Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equations of fluid dynamics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Experimental observation and numerical investigation of filamentary structures in magnetized plasmas

Filamentary structures in low-pressure, low-temperature plasmas are produced when strong magnetic fields are applied parallel to the electric field defined by parallel electrodes. Filamentary structures are regions within the plasma that have distinct properties such as optical brightness and extend along the magnetic field lines. In our experiments, an argon, radio frequency discharge is exposed to a strong background magnetic field in the magnetized dusty plasma experiment at Auburn University. Different forms of filamentary structures emerge in the plasma such as columns, target, or spiral-like structures. To investigate the origin and the characteristics of these patterns, we have developed a three-dimensional fluid model that can reproduce the experimental observations which enables us to investigate the underlying mechanisms of the filamentation process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of ignition and burn propagation in inertial fusion implosions

A detailed analytic model is presented here to investigate the physics of burn propagation in inertially confined plasmas. The onset of ignition and burn propagation occurs when alpha heating of the hot spot causes rapid ablation of shell mass into the hot spot. This allows large energy gains to be achieved since most of the fuel mass is located in the shell. Here, we first present a comprehensive review of previous analytic models that have been used to describe the physics of hot-spot evolution and ignition; we then show that a proper description of a propagating burn wave requires a comprehensive model of hot spot and shell evolution that includes proper mass conservation in the shell, fusion reactivity, and fuel depletion. The analytic theory is in good agreement with detailed radiation-hydrodynamic simulations that predict the onset of burn propagation as occurring when the yield enhancement caused by alpha heating is between 15- and 25-fold, $f_α$ ~ 1.4, where $f_α$ = alpha energy deposited/hot-spot energy at bang time, and the hot-spot burnup fraction is approximately 2%. We show that the definition of ignition is not sensitive to the alpha-particle stopping power nor asymmetries provided that the absorbed fraction of alpha particles $θ_α$ is correctly accounted for. Finally, we use the results of 2-D simulations to show that even when $θ_α$ is small and unknown (as is true in hot spots with mid modes that have significant leakage of alpha particles into the surrounding cold bubbles), one can still relate the experimentally measureable parameter $\chi^{53}_α$ to the yield amplification and the burning-plasma parameter $Q^{\text{hs}}_α$ = alpha energy deposited/total input work delivered to the hot spot

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Generation of excited species in a streamer discharge

At or near atmospheric pressure, most transient discharges, particularly in molecular gases or gas mixture containing molecular gases result in a space-charge dominated transport called a streamer discharge. The excited species generation in such discharges forms the basis for plasma chemistry in most technological applications. In this paper we simulate the propagation of streamers in atmospheric pressure N 2 to understand the energy partitioning in the formation of various excited species and compare the results to a uniform Townsend discharge. The model is fully two-dimensional with azimuthal symmetry. Results show a significantly larger fraction of the energy goes into vibrational excitation of the N 2 ground state in a streamer-type discharge in comparison to a Townsend discharge. For lower applied voltages, the anode-directed (negative) steamer is slightly more efficient in channeling energy for excited species production in comparison to a cathode-directed (positive) streamer. Near 70% overvoltage both types of streamers show very similar energy partitioning but quite different from a Townsend discharge.

42 ENGINEERING↗
DOE OSTI2021

Laser-plasma acceleration beyond wave breaking

Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. Aside from limiting the maximum electric field, this trapping can lead to accelerated electron bunches with large energy spreads. Here we introduce a novel regime of plasma wave excitation and wakefield acceleration that allows for arbitrarily high electric fields while avoiding the deleterious effects of unwanted trapping. Here, the regime, enabled by spatiotemporal shaping of laser pulses, exploits the property that nonlinear plasma waves with superluminal phase velocities cannot trap charged particles and are therefore immune to wave breaking.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scaling of the mean transverse flow and Reynolds shear stress in turbulent plane jet

Here, proper scaling for the mean transverse flow and Reynolds shear stress in a turbulent plane jet is determined using a scaling patch approach. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity equation, a proper scale for the mean transverse flow in a turbulent plane jet is found as V ref = –δdU ctr /dx, where δ is the jet half width and dU ctr /dx is the decay rate of the mean axial velocity at the jet centerline. By seeking an admissible scaling for the mean axial momentum equation, a proper scale for the kinematic Reynolds shear stress is found as R uv,ref =U ctr V ref , which is a mix of the velocity scales in the axial and transverse directions. Approximation functions for the scaled mean transverse flow and Reynolds shear stress are developed and found to agree well with experimental and numerical data. Similarities and differences between the scales of the mean transverse flow and Reynolds shear stress in turbulent plane jets and zero-pressure-gradient turbulent boundary layer flows are clarified.

42 ENGINEERING↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗