Search NASA⌕ Search

SEARCH · Search NASA

Results for “HYDRODYNAMIC EQUATION”

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

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2023

Scissors modes of a Bose-Einstein condensate in a synthetic magnetic field

We study the scissors modes of a harmonically trapped Bose-Einstein condensate under the influence of a synthetic magnetic field, which induces rigid rotational components in the velocity field. Our investigation reveals that the scissors mode, excited in the plane perpendicular to the synthetic magnetic field, becomes coupled to the quadrupole modes of the condensate, giving rise to typical beating effects. Moreover, the two scissors modes excited in the vertical planes are also coupled together by the synthetic magnetic field, resulting in intriguing gyroscope dynamics. Our analytical results, derived from a spinor hydrodynamic theory, are further validated through numerical simulations of the three-dimensional Gross-Pitaevskii equation. These predictions for the condensates subject to a synthetic magnetic field are experimentally accessible with current cold-atom setups and hold promise for potential applications in quantum sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗
DOE OSTI2025

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Lastly, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Chiral hydrodynamics in strong external magnetic fields

We construct the general hydrodynamic description of (3+1)-dimensional chiral charged (quantum) fluids subject to a strong external magnetic field with effective field theory methods. We determine the constitutive equations for the energy-momentum tensor and the axial charge current, in part from a generating functional. Furthermore, we derive the Kubo formulas which relate two-point functions of the energy-momentum tensor and charge current to 27 transport coefficients: 8 independent thermodynamic, 4 independent non-dissipative hydrodynamic, and 10 independent dissipative hydrodynamic transport coefficients. Five Onsager relations render 5 more transport coefficients dependent. We uncover four novel transport effects, which are encoded in what we call the shear-induced conductivity, the two expansion-induced longitudinal conductivities and the shear-induced Hall conductivity. Remarkably, the shear-induced Hall conductivity constitutes a novel non-dissipative transport effect. As a demonstration, we compute all transport coefficients explicitly in a strongly coupled quantum fluid via holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Performance Evaluation of a Two-Dimensional Flood Model on Heterogeneous High-Performance Computing Architectures

This paper describes the implementation of a two-dimensional hydrodynamic flood model with two different numerical schemes on heterogeneous high-performance computing architectures. Both schemes were able to solve the nonlinear hyperbolic shallow water equations using an explicit upwind first-order approach on finite differences and finite volumes, respectively, and were conducted using MPI and CUDA. Four different test cases were simulated on the Summit supercomputer at Oak Ridge National Laboratory. Both numerical schemes scaled up to 128 nodes (768 GPUs) with a maximum 98.2x speedup of over 1 GPU. The lowest run time for the 10 day Hurricane Harvey event simulation at 5 meter resolution (272 million grid cells) was 50 minutes. GPUDirect communication proved to be more convenient than the standard communication strategy. Both strong and weak scaling are shown.

Sharif, Md Bulbul↗
DOE OSTI2022

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatio-temporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Furthermore, our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

An efficient reconstruction algorithm for diffusion on triangular grids using the nodal discontinuous Galerkin method

High-energy-density (HED) hydrodynamics studies such as those relevant to inertial confinement fusion and astrophysics require highly disparate densities, temperatures, viscosities, and other diffusion parameters over relatively short spatial scales. This presents a challenge for high-order accurate methods to effectively resolve the hydrodynamics at these scales, particularly in the presence of highly disparate diffusion. A significant volume of engineering and physics applications use an unstructured discontinuous Galerkin (DG) method developed based on the finite element mesh generation and algorithmic framework. This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection–diffusion equations. However, to the best of the authors’ knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. As a result, the algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.

Computational efficiency↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗