Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW, COMPRESSIBLE - TWO-DIMENSIONAL”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Observations on the deformation of metal microspheres in shock-driven polymer flows

We report that solid particles can be fragmented by a fast-moving fluid if their velocity difference is great enough, such as during the atmospheric entry of meteoroids or the shock compression of engineered particulate composites. The extent of particle deformation and breakup in such systems is poorly understood because the necessary extreme conditions make observation difficult and data scarce. To meet this need, experiments combining ultrafast synchrotron-based radiography with plate impact loading were performed at the dynamic compression sector at the advanced photon source. Metal microspheres of several densities and strengths (Au, Ta, and W) were placed inside a polymer matrix. A planar shock wave was then produced in the polymer by the impact of a gun-launched flyer plate. X-ray images of the resulting flow were collected at ~150ns intervals. These images document the progression of particle deformation across a range of flow conditions and particle materials. They show that the extent of deformation is sensitive to the ratio of drag stress to particle strength. The deforming particle's shape is determined by the initial shock–particle interaction, fluid stagnation pressure, and vorticity, each acting on its own timescale. A set of scaling relationships is presented to capture these observations and enable comparison with prior hydrodynamic data. The result is a framework for predicting the conditions under which strong particles are severely deformed by a shock-driven flow.

36 MATERIALS SCIENCE↗
DOE OSTI2022

X-ray diffraction from shock driven Sn microjets

Here, in situ x-ray diffraction was performed on shock-generated microjets composed of Sn and Sn–4Ag. Under low pressure drives (~21 GPa), a significant fraction of the jet volume was found to be in the β-Sn phase, and these crystallites were much smaller than the initial grain size of the material. Significant quantities of amorphous (molten) material were observed for higher drive pressures (~25–35 GPa). The extent of melting at these pressures was greater than would be predicted for uniaxial shock loading. Diffraction patterns from the Sn–4Ag alloy showed a peak that is consistent with the expected Ag 3 Sn intermetallic phase. This peak remained evident under drive conditions where the sample was otherwise fully amorphous. This indicates a slushy or a mixed phase of liquid Sn and solid Ag 3 Sn. Given the eutectic character of this alloy, this observation is attributed to a kinetic limitation on the dissolution of Ag 3 Sn. This implies that a much broader range of drive conditions will lead to mixed phase jets and ejecta than would be predicted from the equilibrium melt boundary of such alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
DOE OSTI2023

Latest developments in the MOOSE fluid properties module

The fluid properties module in MOOSE serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, SAM~\cite{sam}, the MOOSE thermal hydraulics module, RELAP-7~\cite{relap7} and subchannel~\cite{subchannel}. It is used for coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled simulations~\cite{anl_sam_pgh}. The module offers a consistent set of interfaces to implement to create a new fluid property. There are numerous fluid properties of interest in the entirety of all fields of fluid flow simulations, and this is exacerbated by the use of different variable sets depending on the compressibility of the fluid. For single-phase fluids, the following variable sets may be used to compute fluid properties: (pressure, temperature) and (specific volume, specific internal energy). Some properties may also be computed using the (pressure, density) or the (specific volume, specific enthalpy) variable sets. In order to reduce the challenge of adding a new fluid property, properties may be implemented partially.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
DOE OSTI2023

Discrete element modeling of granular hopper flow of irregular-shaped deformable particles

Many natural and engineered granular materials have relatively deformable particles. Besides particle size and shape, particle deformability is another salient factor that significantly impacts the material’s flow behavior. Here, in this work, the flow of irregular-shaped deformable particles in a wedge-shaped hopper is investigated using discrete element simulations. A bonded-sphere model is developed to simultaneously capture irregular particle shapes and particle-wise deformations (e.g., compression, deflection, and distortion). Quantitative analysis of the effects of irregular shapes and particle deformations shows that the increase in particle stiffness tends to increase initial packing porosity and decrease the flow rate in the hopper. Rigid particles tend to have clogging issues, whereas deformable particles have less chance to, indicating particle deformation reduces the critical bridging width in the hopper flow. Detailed analysis of stress fields is also conducted to provide insights into the mechanism of particle flow and clogging. Stresses and discharge rates calculated from numerical simulations are compared and show good agreement with Walker’s theory and the extended Beverloo formula. Simulations with various particle shape combinations are also performed and show that the initial packing porosity decreases with an increasing percentage of fibers while the discharge rate has a complex dependency on particle shapes.

09 BIOMASS FUELS↗
DOE OSTI2021

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗
DOE OSTI2020

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Nonidealities in Rotating Detonation Engines

A rotating detonation engine (RDE) is a realization of pressure-gain combustion, wherein a traveling detonation wave confined in a chamber provides shock-based compression along with chemical heat release. Due to the high wave speeds, such devices can process high mass flow rates in small volumes, leading to compact and unconventional designs. RDEs involve unsteady and multiscale physics, and their operational characteristics are determined by an equilibrium between large- and small-scale processes. While RDEs can provide a significant theoretical gain in efficiency, achieving this improvement requires an understanding of the multiscale coupling. Specifically, unavoidable nonidealities, such as unsteady mixing, secondary combustion, and multiple competing waves associated with practical designs, need to be understood and managed. The secondary combustion processes arise from fuel/air injection and unsteady and incomplete mixing, and can create spurious losses. In addition, a combination of multiple detonation and secondary waves compete and define the dynamical behavior of mixing, heat release distribution, and the overall mode of operation of the device. This review discusses the current understanding of such nonidealities and describes the tools and techniques used to gain insight into the extreme unsteady environment in such combustors.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Thermal field reconstruction and compressive sensing using proper orthogonal decomposition

Model order reduction allows critical information about sensor placement and experiment design to be distilled from raw fluid mechanics simulation data. In many cases, sensed information in conjunction with reduced order models can also be used to regenerate full field variables. In this paper, a proper orthogonal decomposition (POD) inferencing method is extended to the modeling and compressive sensing of temperature, a scalar field variable. The method is applied to a simulated, critically stable, incompressible flow over a heated cylinder (Re = 1000) with Prandtl number varying between 0.001 and 50. The model is trained on pressure and temperature data from simulations. Field reconstructions are then generated using data from selected sensors and the POD model. Finally, the reconstruction error is evaluated across all Prandtl numbers for different numbers of retained modes and sensors. The predicted trend of increasing reconstruction accuracy with decreasing Prandtl number is confirmed and a Prandtl number/sensor count error matrix is presented.

Matulis, John↗
DOE OSTI2020

Discrete element modeling of switchgrass particles under compression and rotational shear

Switchgrass is a perennial herbaceous plant regarded as a biomass energy crop in the United States for its highadaptability and yield potential. Processing and handling of switchgrass particles are challenging due to the erratic mechanical and flow behavior originating from their intrinsic particulate properties. Here, we present a bonded-sphere discrete element model designed specifically for switchgrass particles. The model simultaneously captures three key particulate features, i.e., fibrous particle shapes, a wide range of particle sizes, and particle deformability. Realistic yet computationally efficient particle shape templates are created based on the image analysis data of switchgrass specimens. A fitting procedure is proposed to ensure both the particle width and length distributions are captured, a unique requirement for fibrous particles. Two full-scale numerical models, i.e., a uniaxial compression model and a Schulze ring shear model, are developed using information fromphysical experiments. The model is calibrated using experimental data of chopped-small switchgrass specimens, and then, is validated using data of chopped-large specimens in both compression and ring-shear tests. Numerical results show that the numerical models capture bulk densities accurately (with an error of 3%) while slightly underestimate the bulk friction angle. Furthermore, an extensive sensitivity analysis reveals that (1) switchgrass particles with rougher edges (due to different processing techniques) exhibit a higher shear strength and a lower flowability; (2) stiffer particles yield a lower bulk density (up to 21% lower) compared to more deformable particles, indicating particle deformability should be incorporated when modeling biomass flow in a preprocessing system.

09 BIOMASS FUELS↗
DOE OSTI2020

An Evaluation of Energy Storage Cost and Performance Characteristics

The energy storage industry continues to expand globally as costs continue to fall and opportunities in consumer, transportation and grid applications are defined. As the rapid evolution of the industry continues, it has become increasingly important to understand how varying technologies compare in terms of cost and performance. This paper defines and evaluates cost and performance parameters of six battery energy storage technologies (BESS) - lithium-ion batteries, lead-acid batteries, redox flow batteries, sodium-sulfur batteries, sodium metal halide batteries, and zinc-hybrid cathode batteries - four non-BESS storage technologies - pumped storage hydropower, flywheels, compressed air energy storage, and ultracapacitors - and combustion turbines. Cost and performance information was compiled based on an extensive literature review, conversations with vendors and stakeholders, and costs of systems procured at sites across the United States. Detailed cost and performance estimates are presented for 2018 and projected out to 2025. Annualized costs were also calculated for each technology.

Mongird, Kendall↗
DOE OSTI2025

Ductility evaluation of nitrided Zircaloy-4 cladding by ring compression testing

Here, this study evaluated the embrittlement of Zircaloy-4 cladding by dry air and dry nitrogen environments as a function of process temperature and duration of exposure. A Netzsch thermogravimetric analyzer was used to measure the mass gain of cladding sections over time with specified temperature and gas flow programs. Oxidation of Zircaloy in dry air produced a brittle ZrO 2 layer with ductile unreacted metal. In contrast, nitridation of Zircaloy had a slower reaction rate but produced a brittle cladding that was more easily breached during ring compression tests. The ring compression tests demonstrated that ductility of Zircaloy cladding decreases as nitrogen pickup increases; a ductile-to-brittle transition occurs near a mass gain of 2.1 wt %. Nitridation showed minimal temperature dependency in ductility degradation, whereas oxidation showed significant temperature dependency in ductility degradation, with severe embrittlement occurring at 1200 °C. Optical micrographs confirmed the formation of uniform ZrN layers (approximately 12 µm thick) on the cladding surface, indicating that cladding embrittlement is primarily the result of nitrogen diffusion rather than metal layer thinning.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ultralow-temperature cryogenic transmission electron microscopy using a new helium flow cryostat stage

Advances in cryogenic electron microscopy have opened new avenues for probing quantum phenomena in correlated materials. This study reports the installation and performance of a new side-entry condenZero cryogenic cooling system for JEOL (Scanning) Transmission Electron Microscopes (S/TEM), utilizing compressed liquid helium (LHe) and designed for imaging and spectroscopy at ultra-low temperatures. The system includes an external dewar mounted on a vibration-damping stage and a pressurized, low-noise helium transfer line with a remotely controllable needle valve, ensuring stable and efficient LHe flow with minimal thermal and mechanical noise. Performance evaluation demonstrates a stable base temperature of 4.37 K measured using a Cernox bare chip sensor on the holder with temperature fluctuations within ±0.004 K. Complementary in-situ electron energy-loss spectroscopy (EELS) via aluminum bulk plasmon analysis was used to measure the local specimen temperature and validate cryogenic operation during experiments. The integration of cryogenic cooling with other microscopy techniques, including electron diffraction and Lorentz TEM, was demonstrated by resolving charge density wave (CDW) transitions in NbSe2 using electron diffraction, and imaging nanometric magnetic skyrmions in MnSi via Lorentz TEM. In conclusion, this platform provides reliable cryogenic operation below 7 K, establishing a low-drift route for direct visualization of electronic and magnetic phase transformations in quantum materials.

Charge density wave↗