Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct Numerical Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Effectiveness of multi-stage cooling processes in improving the CH 4 -hydrate saturation uniformity in sandy laboratory samples

Laboratory-created samples of methane hydrate (MH)-bearing media are a necessity because of the rarity and difficulty of obtaining naturally-occurring samples. The hypothesis that the inevitable heterogeneity in the phase saturations of the laboratory samples may lead to unreliable and non-repeatable results provided the impetus for us to undertake this study, which aimed to determine the conditions under which maximum uniformity can be achieved. To that end, we designed four experiments involving different multi-stage cooling regimes (in terms of their duration and number of stages) to induce MH formation under excess-water conditions. In the absence of direct visualization capabilities, we analysed the experimental results by means of numerical simulation, which provided high-resolution predictions of the spatial distributions of the phase saturations in the cores and enabled the estimation of the parameters controlling the kinetic MH-formation behaviour through history-matching. Analysis of the numerical results indicated that, under the conditions of the experiments and with the design of the reactor, significant heterogeneities in phase saturation distributions were observed in all cases, leading to the conclusion that it is not possible to obtain cores with uniform phase saturation. Additionally, contrary to expectations, heterogeneities increased with the number of cooling stages and the duration of cooling, and this was attributed to imperfect insulation of the upper part of the reactor. A set of simulations involving perfect insulation of the reactor top confirmed the validity of this assumption: (a) predicting the formation of high-uniformity MH-bearing cores that became more homogeneous as the number of cooling stages and the length of the cooling period increased; and (b) providing important information for the improvement of the standard design of the experimental apparatus for the laboratory creation of MH-bearing cores using the excess water method.

42 ENGINEERING↗
DOE OSTI2023

Sensitivity of Precipitation Displacement of a Simulated MCS to Changes in Land Surface Conditions

Abstract This study investigates the role of the land surface on the precipitation produced by an elevated mesoscale convective system (MCS) in Iowa between 24–25 June 2015 during the Plains Elevated Convection at Night (PECAN) field campaign. Previous studies have shown a strong effect of low‐level atmospheric moisture on the location of this MCS. A series of semi‐idealized and realistic simulations with irrigation are conducted to understand the effect of moisture perturbations on the MCS precipitation displacement. In general, numerical simulations place the MCS east of the observed location. Adding moisture directly in the low‐level atmosphere in the semi‐idealized experiments reduces this displacement error. However, experiments with perturbed soil moisture result in drying over Iowa induced by moisture flux divergence from cooler low‐level temperatures and higher surface pressure, causing the MCS to move further to the east. The irrigation impact on low‐level moisture is highly dependent on the length of simulation period. Shorter simulations on the order of days generate similar drying over Iowa but the opposite is found for month‐long simulations. Despite the lack of low‐level moistening in the perturbed soil moisture and short‐term irrigation experiments, the sensitivity to low‐level moisture is similar in all runs. More low‐level moisture generates a more convectively unstable environment with less inhibition and a lower level of free convection that leads to more rapid MCS development and a change in MCS location.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulating lossy Gaussian boson sampling with matrix-product operators

Gaussian boson sampling, a computational model that is widely believed to admit quantum supremacy, has already been experimentally demonstrated and is claimed to surpass the classical simulation capabilities of even the most powerful supercomputers today. However, whether the current approach limited by photon loss and noise in such experiments prescribes a scalable path to quantum advantage is an open question. Here, to understand the effect of photon loss on the scalability of Gaussian boson sampling, we analytically derive the asymptotic operator entanglement entropy scaling, which relates to the simulation complexity. As a result, we observe that efficient tensor network simulations are likely possible under the N out ∝ √N scaling of the number of surviving photons N out in the number of input photons N. We numerically verify this result using a tensor network algorithm with U⁡(1) symmetry, and we overcome previous challenges due to the large local Hilbert-space dimensions in Gaussian boson sampling with hardware acceleration. Additionally, we observe that increasing the photon number through larger squeezing does not increase the entanglement entropy significantly. Finally, we numerically find the bond dimension necessary for fixed accuracy simulations, providing more direct evidence for the complexity of tensor networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Probing Dark Matter Using Precision Measurements of Stellar Accelerations

Dark matter comprises the bulk of the matter in the Universe, but its particle nature and cosmological origin remain mysterious. Knowledge of the dark matter density distribution in the Milky Way Galaxy is crucial both to our understanding of the standard cosmological model and for grounding direct and indirect searches for the particles comprising dark matter. Current measurements of Galactic dark matter content rely on model assumptions to infer the forces acting upon stars from the distribution of observed velocities. Here, we propose to apply the precision radial velocity method, optimized in recent years for exoplanet astronomy, to measure the change in the velocity of stars over time, thereby providing a direct probe of the local gravitational potential in the Galaxy. Using numerical simulations, we develop a realistic strategy to observe the differential accelerations of stars in our Galactic neighborhood with next-generation telescopes, at the level of 10 –8 cm/s 2 . Our simulations show that detecting accelerations at this level with an ensemble of 10 3 stars requires the effect of stellar noise on radial velocity measurements to be reduced to < 10 cm/s. Furthermore, the measured stellar accelerations may then be used to extract the local dark matter density and morphological parameters of the density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Evaluation of Machine Learning Approaches to Estimate Aerosol Mixing State Metrics in Atmospheric Models

Aerosol mixing state describes how aerosol compositions are distributed among atmospheric aerosol particles in a population. Oversimplified assumptions of aerosol mixing state in atmospheric modeling can introduce errors in estimations of weather and climate-relevant aerosol microphysical properties. A more comprehensive representation of the aerosol mixing state can be achieved in principle with a Particle-resolved Monte Carlo (PartMC) model but at added computational cost that may be prohibitive for direct invocation in operational numerical weather prediction or multi-year climate simulations. The aim of our research is to explore the machine learning (ML) methodologies for estimating aerosol mixing state metrics, which we define here in three different ways: with respect to (a) hygroscopicity; (b) optical properties; and (c) chemical species abundance. We adopted a data-driven approach, leveraging deep learning and statistical learning techniques, to take advantage of massive PartMC model simulations. First, we performed particle-resolved simulations by PartMC to create a series of scenarios considering a range of global environmental conditions. Each scenario consists of aerosol populations and corresponding mixing state metrics. The gas concentration, aerosol mass concentration, environmental variables and mixing state metrics of each population constitute the datasets for machine learning implementations. We have adopted and evaluated various configurations of machine learning methodologies in this investigation, embracing deep learning, Extreme Gradient Boosting (XGBoost) algorithm, and the ensemble approaches. After a rigorous model selection process, we identified an appropriate model to derive estimates of aerosol mixing state metrics. We used the computational and data resources of the Oak Ridge Leadership Computing Facility (OLCF) and the ORNL Compute and Data Environment for Science (CADES). The NVIDIA DGX-1 hardware was used for the prototyping of the ML models. Our approach has allowed us to gain a new understanding of how machine learning methodologies can be applied to improve the representation of aerosol mixing state in atmospheric models and benefit the atmospheric research community. Next, we plan to extend our research and methodology to quantify some of the aerosol-related uncertainties in the E3SM Atmospheric Model (EAM) Version 1.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Study of incompressible turbulent channel flow by large eddy simulation

Incompressible turbulent channel flow is investigated by large eddy simulation using improved numerical methods and boundary conditions. In downstream and spanwise directions, cyclic boundary conditions are imposed for velocity and pressure, and two types of boundary conditions near the wall are used and compared. One type is based on the logarithmic law of the mean velocity near the wall and has a slip boundary condition where the molecular-viscous term is neglected. The other type is based on a no-slip boundary condition, where fine mesh spacing near the wall is used to take account of the molecular viscosity. Although the present study employs a coarse mesh (16 x 16 x 21), its results are in good agreement with those of Moin and Kim (1981).

Horiuti, K.↗
NASA NTRS1988

Collisionless shock formation and the prompt acceleration of solar flare ions

The formation mechanisms of collisionless shocks in solar flare plasmas are investigated. The priamry flare energy release is assumed to arise in the coronal portion of a flare loop as many small regions or 'hot spots' where the plasma beta locally exceeds unity. One dimensional hybrid numerical simulations show that the expansion of these 'hot spots' in a direction either perpendicular or oblique to the ambient magnetic field gives rise to collisionless shocks in a few Omega(i), where Omega(i) is the local ion cyclotron frequency. For solar parameters, this is less than 1 second. The local shocks are then subsequently able to accelerate particles to 10 MeV in less than 1 second by a combined drift-diffusive process. The formation mechanism may also give rise to energetic ions of 100 keV in the shock vicinity. The presence of these energetic ions is due either to ion heating or ion beam instabilities and they may act as a seed population for further acceleration. The prompt acceleration of ions inferred from the Gamma Ray Spectrometer on the Solar Maximum Mission can thus be explained by this mechanism.

Cargill, P. J.↗
NASA NTRS1990

Magnetic reconnection driven by current repulsion

The evolution of an equilibrium consisting of two magnetic islands with oppositely directed currents embedded in a strong magnetic field is investigated, using numerical simulation methods. The rapid development of an ideal magnetohydrodynamic instability is observed, which first rotates and then expels the islands. The growth rate is on the order of the inverse of the Alfven transit time and is much higher than that for magnetic island coalescence. In the nonlinear stage, resistivity becomes important as the reconnection process ensues and dissipates the magnetic energy. The growth rate of the instability is a weak function of the plasma beta and other plasma parameters such as S, the magnetic Reynolds number. An energy principle analysis, based on eigenfunctions obtained from the simulation, confirms the existence of the instability.

Richard, R. L.↗
NASA NTRS1991

Quasi-static shape estimation and control of adaptive truss structures

Methods for estimating the deformation of adaptive truss structures are proposed which employ internal displacement sensors to measure changes in the length of selected truss members. Based on the measured data from the instrumented truss member, the total truss deformation pattern can be estimated through direct interpolation. To verify the validity of the methods presented here, numerical simulations are carried out for simple plane trusses, a beam truss, and a tetrahedral truss.

Kuwao, Fumihiro↗
DOE OSTI2020

Joining Techniques for Novel Metal Polymer Hybrid Heat Exchangers

In the United States, over 50% of the unrecovered energy from industrial processes is in the form of low-grade heat (< 220°C). Materials and maintenance costs of common heat exchangers are typically too high to justify their usage. Polymers, though more affordable, are usually unsuitable for HX applications due to their low thermal conductivity (~0.2 W/mK). Here, we show that metal-polymer hybrids may be attractive from both performance and cost perspectives. The use of polymers further increases the resistance to corrosion by sulfuric and carbonic acids often present in flue gases. An ongoing work explores different configurations of layered polyimide-copper macroscale hybrids for heat exchanger applications using numerical simulations. This paper explores a manufacturing pathway for producing such layered hybrid tubes that involves directly rolling and bonding tapes made of polymer and copper foil into tubes. A critical problem in the fabrication process is the bonding of metal and polymers. We explore approaches involving adhesives (epoxy, acrylic and silicone) for metal/polymer interfaces and direct welding (ultrasonic) for metal/metal interfaces that can be integrated into the manufacturing process. We report characterizations of the thermomechanical properties of these joining processes. This work paves the way for realizing cost-effective manufacturing of heat exchangers for low grade waste heat recovery.

Kuntumalla, Gowtham↗
DOE OSTI2020

Ultra-high aspect ratio pores milled in diamond via laser, ion and electron beam mediated processes

Microfabrication techniques are critical for the rapid prototyping and development of applications for cutting edge materials. Recently diamond has gained considerable interest for quantum photonic, biosensing, inertial confinement fusion and magnetometer applications. In this article, ultra-high aspect ratio milling of diamond micropores by photon, ion and electron based methods is reported. A multiphoton absorption laser ablation approach is revealed to rapidly produce sub-10 μm diameter micropores in diamond with an aspect ratio of 14:1 and a tapered profile at the surface interface. Chemically active, oxygen focused ion beam milling produces high-aspect ratio pores in diamond with an aspect ratio of 65:1 and minimal tapering over the length of the pore, overcoming the physical interaction volume limitations imposed in conventional gallium based focused ion beam milling and laser ablation methods. Oxygen-mediated electron beam induced etching is revealed to negate the limitations imposed by physical sputtering mechanisms utilized in focused ion beam milling via the direct initiation of chemical reactions at the receding surface, producing aspect ratios on the order of 200:1. Numerical simulations reveal the physical basis for the superior aspect-ratio pore milling of the oxygen focused ion beam milling and electron beam induced etching methods. Our results demonstrate direct-write methods for the fabrication of ultra-high aspect micropores in diamond and provide insight into the underlying mechanisms of these physical processes. The three methods demonstrated here can be interchanged for applications based on the desired characteristic aspect ratio and process throughput.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Observational and Numerical Diagnostics of Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r(sub 200) and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to several recent results, we observe a steepening of the density profiles beyond approximately 0.3r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or gas clumping are in better agreement with the observed gas distribution. We note a systematic difference between cool-core and non-cool core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. The general trend of steepening density around the virial radius indicates that the shallow density profiles found in several recent works were probably obtained along particular directions (e.g., filaments) and are not representative of the typical behavior of clusters. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Therefore, a detailed treatment of gas cooling, star formation, clumping, and AGN feedback is required to construct realistic models of cluster outer regions.

Eckert, D.↗
NASA NTRS1986

Numerical simulation of transitional flow

The applicability of active control of transition by periodic suction-blowing is investigated via direct simulations of the Navier-Stokes equations. The time-evolution of finite-amplitude disturbances in plane channel flow is compared in detail with and without control. The analysis indicates that, for relatively small three-dimensional amplitudes, a two-dimensional control effectively reduces disturbance growth rates even for linearly unstable Reynolds numbers. After the flow goes through secondary instability, three-dimensional control seems necessary to stabilize the flow. An investigation of the temperature field suggests that passive temperature contamination is operative to reflect the flow dynamics during transition.

Biringen, Sedat↗
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
NASA NTRS2000

DSMC Simulations of Shock Interactions About Sharp Double Cones

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 flow about sharp double cones. Computations are made by using the direct simulation Monte Carlo (DSMC) method of Bird. The sensitivity and characteristics of the interactions are examined by varying flow conditions, model size, and configuration. The range of conditions investigated includes those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel.

Moss, James N.↗
NASA NTRS1997

Parallelization of Rocket Engine Simulator Software (PRESS)

Parallelization of Rocket Engine System Software (PRESS) project is part of a collaborative effort with Southern University at Baton Rouge (SUBR), University of West Florida (UWF), and Jackson State University (JSU). The second-year funding, which supports two graduate students enrolled in our new Master's program in Computer Science at Hampton University and the principal investigator, have been obtained for the period from October 19, 1996 through October 18, 1997. The key part of the interim report was new directions for the second year funding. This came about from discussions during Rocket Engine Numeric Simulator (RENS) project meeting in Pensacola on January 17-18, 1997. At that time, a software agreement between Hampton University and NASA Lewis Research Center had already been concluded. That agreement concerns off-NASA-site experimentation with PUMPDES/TURBDES software. Before this agreement, during the first year of the project, another large-scale FORTRAN-based software, Two-Dimensional Kinetics (TDK), was being used for translation to an object-oriented language and parallelization experiments. However, that package proved to be too complex and lacking sufficient documentation for effective translation effort to the object-oriented C + + source code. The focus, this time with better documented and more manageable PUMPDES/TURBDES package, was still on translation to C + + with design improvements. At the RENS Meeting, however, the new impetus for the RENS projects in general, and PRESS in particular, has shifted in two important ways. One was closer alignment with the work on Numerical Propulsion System Simulator (NPSS) through cooperation and collaboration with LERC ACLU organization. The other was to see whether and how NASA's various rocket design software can be run over local and intra nets without any radical efforts for redesign and translation into object-oriented source code. There were also suggestions that the Fortran based code be encapsulated in C + + code thereby facilitating reuse without undue development effort. The details are covered in the aforementioned section of the interim report filed on April 28, 1997.

Cezzar, Ruknet↗
DOE OSTI2018

Calculating electronic stopping power in materials from first principles

From the early models of electronic stopping power to the current first principles simulations, the techniques evolved to increase the range of validity and to reduce empiricism. Thanks to a combination of theoretical advances provided by Time Dependent Density Functional Theory and the development of numerical codes, it is currently possible to predict electronic stopping power for realistic materials by performing direct simulations of the electron excitation processes beyond linear response, and including electronic band structure effects. Electronic stopping power is an important quantity used to predict and understand the effects of particle radiation in matter. First principles calculations of electronic stopping power can be applied to any atomistic system, solids, liquids and alloys. This review aims here to help graduate level students and researchers immerse themselves into state-of-the-art techniques to computationally model and calculate electronic stopping power.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Digital Twin Approach to Study Additive Manufacturing Processing Using Embedded Optical Fiber Sensors and Numerical Modeling

One of the major challenges for metal-powder-based additive manufacturing is measuring and mitigating residual strain induced during the manufacturing processes. This article reports distributed fiber optic sensors embedded in Inconel alloy components as experimental means to validate numerical models of additive manufacturing process. Electroplating was used to deposit a metal protective jacket onto standard telecom single-mode fibers for strain measurements, Fiber sensors were embedded in an Inconel alloy substrate using the laser engineered net shaping (LENS) process. Here by using a Rayleigh-scattering optical frequency domain reflectometer (OFDR), temperature changes, and residual strain in the metal substrate were monitored with 5 mm spatial resolution during the LENS process. Using finite element analysis, temperature and strain profiles induced by the LENS deposition processes were also numerically studied. Discrepancies between the simulated temperature and strain profiles and those measured directly were less than 10%. Results presented in this article demonstrates a digital twin approach to fuse modeling results with distributed fiber sensor measurement data to study additive manufacturing process toward design and fabrication process optimization.

36 MATERIALS SCIENCE↗