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 325 records · Page 18

Reconstructing Richtmyer–Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks

We develop an ML-based approach for density reconstruction based on transformer neural networks. This approach is demonstrated in the setting of ICF-like double shell hydrodynamic simulations wherein the parameters related to material properties and initial conditions are varied. The new method can robustly recover the complex topologies given by the Richtmyer-Meshkoff instability (RMI) from a sequence of hydrodynamic features derived from radiographic images corrupted with blur, scatter, and noise. A noise model is developed to characterize errors in extracting features from synthetic radiographs of the simulated density field. The key component of the network is a transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model. This approach is shown to exhibit an excellent ability to accurately recover the RMI growth rates, despite the gas-metal interface being greatly obscured by radiographic noise. Our approach can be applied in a broad array of fields involving shock physics and material science.

47 OTHER INSTRUMENTATION↗

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A renormalization group analysis of two-dimensional magnetohydrodynamic turbulence

The renormalization group (RNG) method is used to study the physics of two-dimensional (2D) magnetohydrodynamic (MHD) turbulence. It is shown that, for a turbulent magnetofluid in two dimensions, no RNG transformation fixed point exists on account of the coexistence of energy transfer to small scales and mean-square magnetic flux transfer to large scales. The absence of a fixed point renders the RNG method incapable of describing the 2D MHD system. A similar conclusion is reached for 2D hydrodynamics, where enstrophy flows to small scales and energy to large scales. These analyses suggest that the applicability of the RNG method to turbulent systems is intrinsically limited, especially in the case of systems with dual-direction transfer.

Liang, Wenli Z.↗

Hydrodynamics of higher-rank gauge theories

We extend recent work on hydrodynamics with global multipolar symmetries — known as “fracton hydrodynamics” — to systems in which the multipolar symmetries are gauged. We refer to the latter as “fracton magnetohydrodynamics”, in analogy to conventional magnetohydrodynamics (MHD), which governs systems with gauged charge conservation. We show that fracton MHD arises naturally from higher-rank Maxwell’s equations and in systems with one-form symmetries obeying certain constraints; while we focus on “minimal” higher-rank generalizations of MHD that realize diffusion, our methods may also be used to identify other, more exotic hydrodynamic theories (e.g., with magnetic subdiffusion). In contrast to semi-microscopic derivations of MHD, our approach elucidates the origin of the hydrodynamic modes by identifying the corresponding higher-form symmetries. Being rooted in symmetries, the hydrodynamic modes may persist even when the semi-microscopic equations no longer provide an accurate description of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MHD modelling of the earth's magnetosphere

Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.

Wu, C. C.↗

Anomalous Behavior of Newtonian Hydrodynamics Coupled with Radiation Transport

This study shows that Newtonian hydrodynamics coupled with radiation transport, using a wide range of methods for treating the material-motion corrections, results in anomalous behavior. In particular, the flow of infinite-medium equilibration will accelerate whenever viewed moving past an observer in the laboratory frame. The acceleration may cause the velocity to increase exponentially in time. An exact fully-relativistic solution is derived to show that there is no acceleration, independent of the reference frame. This solution is analytic and may be used as a validation problem for the simulation of fully-relativistic radiation hydrodynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗

Weighted FFT estimators for 1D and 3D correlations of the Lyman- α forest

Correlations in the Lyman-α (Lyα) forest, both as a function of line of sight separation (1D) and 3D separation, provide a unique window to the distribution of matter at redshifts not accessible by current galaxy surveys. While optimal quadratic estimators have been used to measure 1D correlations, they are computationally expensive and difficult to extend to 3D analyses. On the other hand, estimators based on the Fast Fourier Transform (FFT) are significantly faster, but are affected by missing data in the spectra (masked pixels) and so far have not used pixel weights to reduce the uncertainties in the measurement. In this publication we describe how to compute the window matrix that enables forward-modelling the impact of masked pixels and weights on the FFT-based estimators. Here, we use Gaussian and hydrodynamical simulations with artificially masked pixels to validate the method on the measurement of 1D correlations. Finally, we show that the formalism can be extended to model the impact on 3D correlations, in particular on the cross-spectrum, the correlation of 1D Fourier modes as a function of transverse separation. This work will enable more precise clustering measurements with the Lyα forest dataset recently collected by the Dark Energy Spectroscopic Instrument (DESI).

Lokken, Martine [Univ. Autonoma de Barcelona (Spai↗

Impact of late-time neutrino emission on the diffuse supernova neutrino background

In the absence of high-statistics supernova neutrino measurements, estimates of the diffuse supernova neutrino background (DSNB) hinge on the precision of simulations of core-collapse supernovae. Understanding the cooling phase of protoneutron star (PNS) evolution (≳ 1 s after core bounce) is crucial, since approximately 50% of the energy liberated by neutrinos is emitted during the cooling phase. We model the cooling phase with a hybrid method by combining the neutrino emission predicted by 3D hydrodynamic simulations with several cooling-phase estimates, including a novel two-parameter correlation depending on the final baryonic PNS mass and the time of shock revival. We find that the predicted DSNB event rate at Super-Kamiokande can vary by a factor of ~2–3 depending on the cooling-phase treatment. Here we also find that except for one cooling estimate, the range in predicted DSNB events is largely driven by the uncertainty in the neutrino mean energy. With a good understanding of the late-time neutrino emission, more precise DSNB estimates can be made for the next generation of DSNB searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Diffusion (Mixing) of Material in Numerical Simulations of Hydrodynamics [Slides]

It is often assumed that a material interface between two materials is spread over a few numerical cells in numerical simulations for hydrodynamics. Also, we have the impression that higher order methods introduce less numerical diffusion (mixing) of material. As we know one of the purposes of adaptive mesh refinement (AMR) is to resolve interfaces between materials, but we would like to know how effective AMR is to reduce numerical diffusion of material. We will present our investigation about numerical diffusion (mixing) of material in xRage. The result of the investigation indicates that the assumptions mentioned above are not always valid. In this talk, we will also demonstrate the effectiveness of numerical techniques to reduce numerical diffusion of material, including contact discontinuity steepening, isotropic interface steepening, max interface steepening, material interface reconstruct.

36 MATERIALS SCIENCE↗

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of inducer load and stress, volume 2

A program of analysis, design, fabrication and testing has been conducted to develop computer programs for predicting rocket engine turbopump inducer hydrodynamic loading, stress magnitude and distribution, and vibration characteristics. Methods of predicting blade loading, stress, and vibration characteristics were selected from a literature search and used as a basis for the computer programs. An inducer, representative of typical rocket engine inducers, was designed, fabricated, and tested with special instrumentation selected to provide measurements of blade surface pressures and stresses. Data from the tests were compared with predicted values and the computer programs were revised as required to improve correlation. For Volume 1 see N71-20403. For Volume 2 see N71-20404.

Source record↗

Effects of supply conditions on film thickness in lubricated Hertzian contacts

A generalization of the hydrodynamic expression for Hertzian contacts is described and various methods for calculating the thickness of the oil film winter steady-state, isothermal conditions are given. This is important for engineering applications such as gears and bearings because these results are closer to real operating conditions. Theories of lubrication are discussed, and the mathematics involved are presented using approximately 30 equations and 13 figures. For lubricated, linear, elliptical or point Hertzian contacts it is demonstrated how to calculate the thickness of the oil film at the center of the contact for steady-state isothermal conditions.

Dalmaz, G.↗

Lubrication background

Surface topography, including the various physical methods of measuring surfaces, and the various lubrication regimes (hydrodynamic, elastohydrodynamic, boundary, and mixed) are discussed. The historical development of elastohydrodynamic lubrication is outlined. The major accomplishments in four periods, the pre-1950's, the 1950's, the 1960's, and the 1970's are presented.

Hamrock, B. J.↗

Fragmentation in stellar collapse

A three dimensional Cartesian code was developed to model molecular cloud collapse and protostellar formation. The basic Eulerian scheme employs spatially second-order accurate, finite difference methods to advance the fluid variables in time. The hydrodynamic equations governing the collapse include the effects of self-gravity, rotation, and radiative transfer. Supplementary equations include Poisson's equation which relates the gravitational potential to the density distribution and a mean intensity equation derived with the Eddington approximation for radiative transfer. The gas pressures and internal energies are determined from equations of state derived for a molecular cloud composed of hydrogen, helium, and heavier elements.

Myhill, E. M.↗

qSmoothed particle interpolation

Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of a nonhydrodynamical nature. Under this approach, SPH is viewed as strickly as interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogenous sound speed and (2) Burger's equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite-difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.

Laguna, Pablo↗

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗