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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗
DOE OSTI2025

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗
DOE OSTI2023

A numerical study of an atomizing jet in a resonant acoustic field

By applying state-of-the-art, high-fidelity methods for compressible multiphase flows, this work features a parametric study of a turbulent liquid jet atomizing under the influence of standing acoustic waves with five different amplitude values. We perform these simulations in a semi-periodic domain to isolate interactions between acoustic and atomization processes and to emulate a pure liquid fuel spray placed at the pressure node of a resonant acoustic field. Characterization of the flow includes measurements of the mean dimensions, droplet number, and surface area of the liquid distribution, alongside qualitative depictions of the interface evolution, breakup events, and hydrodynamic fields. As acoustic forcing is provided to the flow, the added energy provokes atomization when shear would otherwise be insufficient for breakup. The magnitude of the forcing determines the rate of droplet production, which accelerates as the evolving liquid surface modifies the induced motion of the gas phase. At the nodal plane, the acoustic field produces an organizing effect on the liquid structures.

42 ENGINEERING↗
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗
DOE OSTI2024

Quantum Liouville theorem based on Haar measure

Liouville theorem (L theorem) reveals robust incompressibility of the distribution function in phase space, given arbitrary potentials. However, its quantum generalization, Wigner flow, is compressible, i.e., L theorem is only conditionally true (e.g., for perfect Harmonic potential). Here, we develop quantum L theorem (rigorous incompressibility) for arbitrary potentials (interacting or not) in Hamiltonians. Haar measure, instead of symplectic measure dp$\bigwedge$dq used in Wigner’s scheme, plays a central role. The argument is based on general measure theory, independent of specific spaces or coordinates. Comparison of classical and quantum is made: for instance, here we address why Haar measure and metric preservation do not work in the classical case. Applications of the theorems in statistics, topological phase transition, ergodic theory, etc., are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Additively Manufactured Compliant Hybrid Gas Thrust Bearing for Supercritical Carbon Dioxide Turbomachinery: Experimental Evaluation and Fluid–Structure Model Predictions

This paper presents rotating test results and advances an analytical predictive fluid–structure model for a new type of gas-lubricated thrust bearing fabricated using direct metal laser melting. The bearing concept in this study is a compliant hybrid gas thrust bearing that uses external pressurization to increase load carrying capacity, where the testing campaign in this study was only focused on steady-state static performance. The need for the bearing concept comes from enabling highly efficient supercritical carbon dioxide (sCO2) turbomachinery by replacing oil-lubricated bearings with process gas lubrication. Leveraging the process gas of the turbomachine for bearing lubrication results in lowered bearing power loss, simplified mechanical design, and allows for novel oil-free hermetic drivetrains resulting in an efficient emission-free system. The new concept utilizes hydrostatic pressurization on individual tilting pads flexibly mounted with hermetic squeeze film dampers (HSFDs). This paper focuses on rotating tests of a 173 mm outer diameter gas thrust bearing in air up to 10 krpm and hydrostatic inlet pressures to 365 psi (2.52 MPa). The influence of thrust runner speed and bearing inlet pressure on force deflection characteristics and load carrying capability of the gas film were experimentally evaluated. This work also advances a predictive fluid–structure thrust bearing model using an isothermal ideal-gas-based compressible Reynolds flow equation directly coupled to a lumped stiffness element possessing axial and rotational degrees-of-freedom. The rotating testing demonstrated load capability of 1816 lbs (8.1 kN), which equates to a thrust bearing unit load of 67 psi (0.46 MPa). Gas film force–deflection curves reveal a nonlinear relationship between thrust load and film clearance. Comparison of film thickness values with the predictive model show good agreement under high load and inlet pressure, however deviate as load and pressure decrease. Load capability was shown to increase with increasing hydrostatic inlet pressure, while the increase in thrust runner speed revealed a small decrease in load capacity.

Engineering↗
DOE OSTI2025

Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization

We present an optimized implementation of the recently proposed information geometric regularization (IGR) for unprecedented scale simulation of compressible fluid flows applied to multi-engine spacecraft boosters. We improve upon state-of-the-art computational fluid dynamics (CFD) techniques in terms of computational cost, memory footprint, and energy-to-solution metrics. Unified memory on coupled CPU–GPU or APU platforms increases problem size with negligible overhead. Mixed half/single-precision storage and computation are used on well-conditioned numerics. We simulate flow at 200 trillion grid points and 1 quadrillion degrees of freedom, exceeding the current record by a factor of 20. A factor of 4 wall-time speedup is achieved over optimized baselines. Ideal weak scaling is observed on OLCF Frontier, LLNL El Capitan, and CSCS Alps using the full systems. Strong scaling is near ideal at extreme conditions, including 80% efficiency on CSCS Alps with an 8 node baseline and stretching to the full system.

Wilfong, Benjamin [Georgia Institute of Technology↗
DOE OSTI2022

NUMERICAL INVESTIGATION OF AIR DEHUMIDIFICATION THROUGH WATER DROPLETS DIELECTROPHORESIS

This paper numerically investigated the dehumidification of moist airflow in a converging rectangular duct with electrically enhanced condensation. The aim was to develop a condensation model that predicted water vapor condensation by employing the dielectrophoresis principle. The charged droplets were injected into the computational domain to produce a gradient electric field. The dielectrophoretic interaction between the droplets and vapor molecules of the humid air resulted in local water vapor condensation at the vapor and liquid droplet interface. This phenomenon is described in detail in this paper. A hybrid Eulerian-Lagrangian solver SprayFoam was developed based on OpenFOAM® to simulate condensation in two-phase gas-liquid mixtures. The following developments were made based on the standard compressible multiphase flow solver SprayFoam in OpenFOAM®: (1) Eulerian solver for gas phase (2) multi-component species transport, and (3) Lagrangian solver for gas-droplet two-phase flows and sub-models for liquid droplets. The newly developed numerical model was experimentally validated with data from a series of tests conducted in the authors' laboratory. The results showed that the simulations followed the same trends as the data, and the model predicted the condensation due to the electro-spray injection process in the air stream. The injected water droplets, which are electrically charged to their Rayleigh limits, increased in size while wiping out the humidity from the air. The simulations indicated that the dehumidification was about 1% for 0.5 cubic feet per minute airflow rate. Scaling up to larger flows is a future follow-up work.

Yel Mahi, *Maliha↗
DOE OSTI2020

Using Differential Forms to Find Symmetries in the Noh Problem for an Ideal Gas in a Spherical System

The Euler equations describe inviscid, compressible gas flow and are derived using the laws of conservation of mass, momentum, and energy. From these equations a common test problem, the Noh problem, can be solved for invariant solutions in planar, cylindrical, and spherical geometries. Using Lie group analysis methods, we will solve the Noh Problem for an ideal gas in a spherical system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

SIERRA/Aero User Manual - Version 4.56

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗
DOE OSTI2020

SIERRA/Aero Theory Manual - Version 4.56

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗