Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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 487 records · Page 27

Imaging the initial condition of heavy-ion collisions and nuclear structure across the nuclide chart

High-energy nuclear collisions encompass three key stages: the structure of the colliding nuclei, informed by low-energy nuclear physics, the initial condition , leading to the formation of quark–gluon plasma (QGP), and the hydrodynamic expansion and hadronization of the QGP, leading to final-state hadron distributions that are observed experimentally. Recent advances in both experimental and theoretical methods have ushered in a precision era of heavy-ion collisions, enabling an increasingly accurate understanding of these stages. However, most approaches involve simultaneously determining both QGP properties and initial conditions from a single collision system, creating complexity due to the coupled contributions of these stages to the final-state observables. To avoid this, we propose leveraging established knowledge of low-energy nuclear structures and hydrodynamic observables to independently constrain the QGP’s initial condition. By conducting comparative studies of collisions involving isobar-like nuclei—species with similar mass numbers but different ground-state geometries—we can disentangle the initial condition’s impacts from the QGP properties. This approach not only refines our understanding of the initial stages of the collisions but also turns high-energy nuclear experiments into a precision tool for imaging nuclear structures, offering insights that complement traditional low-energy approaches. Opportunities for carrying out such comparative experiments at the Large Hadron Collider and other facilities could significantly advance both high-energy and low-energy nuclear physics. Additionally, this approach has implications for the future electron-ion collider. While the possibilities are extensive, we focus on selected proposals that could benefit both the high-energy and low-energy nuclear physics communities. Originally prepared as input for the long-range plan of U.S. nuclear physics, this white paper reflects the status as of September 2022, with a brief update on developments since then.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the convergence of an IEQ-based first-order semi-discrete scheme for the Beris-Edwards system

We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.

Weber, Franziska↗

Initial data for first-order causal viscous conformal fluids in general relativity

We solve the Einstein constraint equations for a first-order causal viscous relativistic hydrodynamic theory in the case of a conformal fluid. For such a theory, a direct application of the conformal method does not lead to a decoupling of the equations, even for constant-mean curvature initial data. We combine the conformal method applied to a background perfect fluid theory with a perturbative argument in order to obtain the result.

Disconzi, Marcelo (ORCID:0000000234497778)↗

Finite amplitude instability of second-order fluids in plane Poiseuille flow.

The hydrodynamic stability of plane Poiseuille flow of second-order fluids to finite amplitude disturbances is examined using the method of Stuart and Watson as extended by Reynolds and Potter. For slightly non-Newtonian fluids subcritical instabilities are predicted. No supercritical equilibrium states are expected if the entire spectrum of disturbance wavelengths is present. Possible implications with respect to the Toms phenomenon are discussed.

Mcintire, L. V.↗

Theoretical mean colors for RR Lyrae variables

A hydrodynamically pulsating 0.6 solar mass model of a typical RR Lyrae variable was studied with a radiation transport-hydrodynamic computer program to predict theoretical T sub 3 and colors at many phases and to find the proper methods for getting mean colors and the consequent mean effective temperatures. The variable Eddington radiation approximation method was used with gray and with multifrequency absorption coefficients to represent the radiation flow in the outer optically thin layers. Comparison between observed and computed B-V colors indicate that these low Z population 2 models are reasonably accurate using King 1A composition opacities. The well known Oke, Giver, and Searle relation between B-V and T sub e reproduced. Mean colors were found by four different averaging methods. The method that gives a mean color and the mean T sub e closest to the nonpulsating model was the separate intensity means of B and V.

Davis, C. G.↗

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗

The Formation of Preplanetary Disks from the Collapse of Rotating Molecular Cloud Cores

Solutions that describe the collapse of a molecular cloud core that is initially in unstable equilibrium, embedded within an envelope of uniform density, and rotating at the same rate as the envelope are given. Hydrodynamic equations, including self gravity, are deduced to a set of ordinary differential equations, which are solved by the method of matched asymptotic expansions. Results of these calculations are: (1) the range of stellar masses derived seems to correspond to realistic ranges of observed stellar masses and interstellar cloud parameters, (2) the proper measure of dissipation rate is the ratio of accretion time to viscous diffusion time, and (3) the pressure distribution on the surface of an accreting protostar is nonuniform in a way that favors the channeling of a stellar wind into a bipolar flow directed along the rotation axis.

Cassen, P.↗

Problems of space biology. Volume 50: Nystagmometry for evaulation of the status of the vestibular function

Various aspects of nystagmometry are studied, primarily those in which the study of hystagmus serves as a means to learn about the vestibular apparatus. Along with exhaustive published material, the monograph presents data from many years of research on the physioloigical mechanisms of nystagmus, the features of nystagmus when vestibular stimulation is combined with optokinetic, the pole of vertibular afferentation asymmetry in the asymmetry of reactions to optokinetic stimulus, a nystagmometric approach to studying the hydrodynamic interaction among semicircular canals, as well as several other questions. A great deal of attention is given to methods of recording nystagmus, calibrating nystagmograms, quantitative evaluation of nystagmographic material, new nystagmometric characteristics and diagnostic techniques. A diagnostic model is proposed which makes it possible to obtain important information on the condition of the vestibular system from results of vestibular testing.

Levashov, M. M.↗

A procedure for computing surface wave trajectories on an inhomogeneous surface

Equations are derived for computing surface waves on smooth surfaces, including surfaces with a nonuniform wave speed. The prior literature dealt primarily with the theoretical development with little consideration given to computational methods, and examples were limited to waves on surfaces of simple analytic description, such as cones, spheres, and cylinders. The computational procedure presented is a relatively general method. Sample calculations illustrate the procedure for a class of practical shapes of the type that include aerodynamic and hydrodynamic surfaces. Equations are also included for computing the spreading of rays into a surrounding medium that will support waves.

Barger, Raymond L.↗

Potential flow in engine valves

The extensive applicability of the hydrodynamic theory to the problems of engine construction is clearly shown in the following attempt to determine by exact methods the nature of the flow in valves under variously restricted conditions. Observation shows that two principal kinds of flow occur in simple flat-seated valves. For small valve lifts, the flow is along the horizontal wall and is therefore deflected 90 degrees, but for greater valve lifts the flow separates and forms a free stream, whose angle of deflection naturally increases with increasing lift. Both these kinds of flow can, in fact, be theoretically explained

THEORIES - FLOW - KIRCHOFF AND HELMHOLTZ↗

Load Capacity Estimation of Foil Air Journal Bearings for Oil-Free Turbomachinery Applications

This paper introduces a simple "Rule of Thumb" (ROT) method to estimate the load capacity of foil air journal bearings, which are self-acting compliant-surface hydrodynamic bearings being considered for Oil-Free turbo-machinery applications such as gas turbine engines. The ROT is based on first principles and data available in the literature and it relates bearing load capacity to the bearing size and speed through an empirically based load capacity coefficient, D. It is shown that load capacity is a linear function of bearing surface velocity and bearing projected area. Furthermore, it was found that the load capacity coefficient, D, is related to the design features of the bearing compliant members and operating conditions (speed and ambient temperature). Early bearing designs with basic or "first generation" compliant support elements have relatively low load capacity. More advanced bearings, in which the compliance of the support structure is tailored, have load capacities up to five times those of simpler designs. The ROT enables simplified load capacity estimation for foil air journal bearings and can guide development of new Oil-Free turbomachinery systems.

DellaCorte, Christopher↗

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Theory of the Unsteady Motion of an Airfoil

The paper presents a systematical analysis of the problem of the determination of the unsteady motion about an airfoil moving in an infinite fluid that contains a system of vortices and the determination of the hydrodynamical forces acting on the airfoil. The hydrodynamical problem is reduced to the determination of the function f (xi) which transforms conformally the external region of the airfoil into the interior of a circle. The proposed methods of determining the irrotational motion of a fluid that is produced by any motion of the airfoil are especially simple and effective if the function f (xi) is rational. As an example the flow is determined for the case of an arbitrary motion of an airfoil of the Joukowsky type. The formulas obtained for the determination of the hydrodynamical forces by means of contour integration are similar to those given by S. Chaplygin. These formulas are used to determine the force acting on the airfoil in the cases where the unsteady motion is potential throughout and the circulation about the airfoil is constant and also when the fluid contains a system of vortices. A full discussion is given of the concept of virtual masses together with practical formulas for computing the virtual mass coefficients.

Sedov, L. I.↗

Shock structure in jets of arbitrary exit geometry

The shock structure and hydrodynamic stability of single and twin jets with arbitrary exit geometry are studied. Both vortex sheets and realistic mean profiles are considered for the jet shear layer. The boundary element method is used to predict the shock spacing and screech tones in a vortex sheet model of a single jet. Results are given for circular, elliptic and rectangular jets and the results compared with analysis and experiment. A finite difference technique is used to examine the effect of finite mixing layer thickness for a single jet. The hydrodynamic stability of twin circular jets is also considered. The variation of the stability characteristics are presented for the vortex sheet case for several separation distances. The formulation of the problem to account for finite thickness effects is also given.

Morris, P. J.↗

Reconstructing Ly$α$ Fields from Low-resolution Hydrodynamical Simulations with Deep Learning

Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to resolve density fluctuation in the IGM puts a stringent requirement on the resolution of such simulations, which in turn limits the volumes that can be modeled, even on the most powerful supercomputers. In this work, we present a novel modeling method that combines physics-driven simulations with data-driven generative neural networks to produce outputs that are qualitatively and statistically close to the outputs of hydrodynamical simulations employing eight times higher resolution. We show that the Ly$α$ flux field, as well as the underlying hydrodynamic fields, have greatly improved statistical fidelity over a low-resolution simulation. Importantly, the design of our neural network allows for sampling multiple realizations from a given input, enabling us to quantify the model uncertainty. Using test data, we demonstrate that this model uncertainty correlates well with the true error of the Ly$α$ flux prediction. Ultimately, our approach allows for training on small simulation volumes and applying it to much larger ones, opening the door to producing accurate Ly$α$ mock skies in volumes of Hubble size, as will be probed with DESI and future spectroscopic sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗

Alpha heating of indirect-drive layered implosions on the National Ignition Facility

We report in order to understand how close current layered implosions in indirect-drive inertial confinement fusion are to ignition, it is necessary to measure the level of alpha heating present. To this end, pairs of experiments were performed that consisted of a low-yield tritium–hydrogen–deuterium (THD) layered implosion and a high-yield deuterium–tritium (DT) layered implosion to validate experimentally current simulation-based methods of determining yield amplification. The THD capsules were designed to reduce simultaneously DT neutron yield (alpha heating) and maintain hydrodynamic similarity with the higher yield DT capsules. The ratio of the yields measured in these experiments then allowed the alpha heating level of the DT layered implosions to be determined. The level of alpha heating inferred is consistent with fits to simulations expressed in terms of experimentally measurable quantities and enables us to infer the level of alpha heating in recent high-performing implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗