Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic coupling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-yield implosion modeling using the Frustraum: Assessing and controlling the formation of polar jets and enhancing implosion performance with applied magnetization

Frustraums have a higher laser-to-capsule x-ray radiation coupling efficiency and can accommodate a large capsule, thus potentially generating a higher yield with less laser energy than cylindrical Hohlraums for a given Hohlraum volume [Amendt et al., Phys. Plasmas 26, 082707 (2019]. Frustraums are expected to have less m = 4 azimuthal asymmetries arising from the intrinsic inner-laser-beam geometry on the National Ignition Facility. An experimental campaign at Lawrence Livermore National Laboratory to demonstrate the high-coupling efficiency and radiation symmetry tuning of the Frustraum has been under way since 2021. Simulations benchmarked against experimental data show that implosions using Frustraums can achieve more yield with higher ignition margins than cylindrical Hohlraums using the same laser energy. Hydrodynamic jets in capsules along the Hohlraum axis, driven by radiation-flux asymmetries in a Hohlraum with a gold liner on a depleted uranium (DU) wall, are present around stagnation, and these “polar” jets can cause severe yield degradation. The early-time Legendre mode P4<0 radiation-flux asymmetry is a leading cause of these jets, which can be reduced by using an unlined DU Hohlraum because the shape of the shell is predicted to be more prolate. Magnetization can increase the implosion robustness and reduce the required hotspot ρR for ignition; therefore, magnetizing the Frustraum can maintain the same yield while reducing the required laser energy or increase the yield using the same laser energy—all under the constraint that the ignition margin is preserved. Reducing polar jets is particularly important for magnetized implosions because of the intrinsic toroidal hotspot ion temperature topology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Inverse Bremsstrahlung Absorption

Inverse bremsstrahlung absorption was measured based on transmission through a finite-length plasma that was thoroughly characterized using spatially resolved Thomson scattering. Expected absorption was then calculated using the diagnosed plasma conditions while varying the absorption model components. To match data, it is necessary to account for: (1) the Langdon effect; (2) laser-frequency (rather than plasma-frequency) dependence in the Coulomb logarithm, as is typical of bremsstrahlung theories but not transport theories; and (3) a correction due to ion screening. Radiation-hydrodynamic simulations of inertial confinement fusion implosions have to date used a Coulomb logarithm from the transport literature and no screening correction. Here, we anticipate that updating the model for collisional absorption will substantially revise our understanding of laser-target coupling for such implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2024

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

CHEQUP v0.1

CHEQUP (Castro-based Hofi Expansion with QUasineutral Plasma) is a simulation code for modeling the formation of hydrodynamic optical-field-ionized (HOFI) plasma channels, which are used as waveguides in laser-plasma acceleration experiments. This includes experiments performed at LBNL's BELLA facility as well as other laser facilities across the world. CHEQUP extends the open-source Castro hydrodynamics framework with physics modules tailored for modeling HOFI plasma channels -- including multi-species ionization and three-body recombination for mixtures of hydrogen, nitrogen, helium, and argon ; a two-temperature model tracking electron and heavy-species temperatures separately ; and coupling with other codes of the BLAST ecosystem (https://blast.lbl.gov/) such as WarpX, via the openPMD standard. CHEQUP inherits from Castro the ability to run on modern GPU architectures (NVIDIA CUDA, AMD HIP) and supports adaptive mesh refinement (AMR) for efficient multi-scale resolution. Compared to existing tools, CHEQUP would be, to our knowledge, the first open-source code implementing the full HOFI channel formation physics, and the first implementation capable of running on GPUs. This enables significantly faster, large-scale parameter scans critical for the design of next-generation LPA-based accelerators and light sources.

Lehe, Remi [Lawrence Berkeley National Laboratory ↗
DOE OSTI2024

Assessment of vertical stability for negative triangularity pilot plants

Abstract Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗
DOE OSTI2023

Colliding localized, lumpy holographic shocks with a granular nuclear structure

We apply a recent and simple technique which speeds up the calculation of localized collisions in holography to study more realistic models of the pre-hydrodynamic phase of heavy ion collisions using gauge/gravity duality. Our initial data reflects the lumpy nuclear structure of real heavy ions and our projectiles’ aspect ratio mimics the Lorentz contraction of nuclei during RHIC collisions. At the hydrodynamization time of the central region of the quark gluon plasma developed during the collision, we find that most of the system’s vorticity is located well outside the hydrodynamized part of the plasma. Only the relativistic corrections to the thermal vorticity within the hydrodynamized region are non-negligible. We compare the transverse flow shortly after the collision with previous results which did not use granular initial conditions and determine the proper energy density and fluid velocity in the hydrodynamized subregion of the plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗
DOE OSTI2025

Nonperturbative heavy-flavor transport approach for hot QCD matter

The heavy charm and bottom quarks are unique probes of the transport properties of the quark-gluon plasma (QGP) and its hadronization in high-energy nuclear collisions. A key challenge in this context is to embed the interactions of the heavy quarks in the expanding medium compatible with the strong-coupling nature of the QGP, and thus to unravel the underlying microscopic mechanisms. In the present work we progress toward this goal by combining recent T -matrix interactions for elastic scattering with an effective transport implementation of gluon radiation, and apply these in a Langevin framework in a viscous hydrodynamic evolution. Hadronization of heavy quarks is evaluated using a modern recombination model with 4-momentum conservation, supplemented with fragmentation constrained by data in proton-proton collisions. Deploying this approach to charm-hadron observables in Pb-Pb collisions at the LHC yields fair agreement with experiment while also identifying areas of further systematic improvement of the simulations and its current input.

Heavy-flavor transport↗
DOE OSTI2024

Enhanced sensitivity to target offset when using cross-beam energy transfer mitigation techniques in direct-drive inertial confinement fusion implosions

In direct-drive inertial confinement fusion, target offset from the target chamber center (or center of beam convergence) may lead to significant implosion asymmetry and fusion yield degradation. In addition, cross-beam energy transfer (CBET) has been shown to be a significant source of laser energy scattering and leads to a reduction in implosion velocity and yield. To improve energy coupling and implosion performance, several techniques for CBET mitigation have been proposed. Recent simulations, however, have shown that CBET also substantially mitigates the effect of target offset on implosion asymmetry and yield [Anderson et al., Phys. Plasmas 27, 112713 (2020)]. Furthermore, the inclusion of CBET models in radiation-hydrodynamics codes was shown to greatly improve agreement between simulations and experiments involving substantial target offset distances. This paper explores the intensity dependence of this CBET–offset effect. In addition, it is shown that enhanced sensitivity to target offset can be expected when CBET-mitigation techniques are used in direct-drive implosions. This is shown through simulations of two such CBET-mitigation techniques on the OMEGA laser: (1) decreased beam-to-target radius, and (2) beam-to-beam frequency detuning. For the typical target offset distances (<15 μm) observed in experiments on OMEGA, however, overall yield is still anticipated to be substantially higher when CBET-mitigation techniques are employed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
DOE OSTI2022

The Three-dimensional Collapse of a Rapidly Rotating 16 M ⊙ Star

I report on the three-dimensional (3D) hydrodynamic evolution of a rapidly rotating 16 M ⊙ star to iron core collapse. For the first time, I follow the 3D evolution of the angular momentum (AM) distribution in the iron core and convective shell burning regions for the final 10 minutes up to and including gravitational instability and core collapse. In 3D, convective regions show efficient AM transport that leads to an AM profile that differs in shape and magnitude from MESA within a few shell convective turnover timescales. For different progenitor models, such as those with tightly coupled Si/O convective shells, efficient AM transport in 3D simulations could lead to a significantly different AM distribution in the stellar interior affecting estimates of the natal neutron star or black hole spin. The results suggest that 3D AM transport in convective and rotating shell burning regions are critical components in models of massive stars and could qualitatively alter the explosion outcome and inferred compact remnant properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples—random charge densities and periodic charge densities—with different scaling properties. Additionally, under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

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

Observation of ponderomotively driven bow shock using Thomson scattering

High-energy speckled lasers are known to exert ponderomotive forces on a plasma. This can reduce flow transverse to the propagation of the beams. When coupled with a supersonic flow, this has been shown to lead to the formation of a shock that travels against the flow. Experiments conducted on the OMEGA laser facility have used Thomson scattering to observe density and velocity changes consistent with this ponderomotively driven shock. Comparisons of the data with hydrodynamic simulations with the ponderomotive force, particle-in-cell simulations with a full Maxwell field solver, and hydrodynamic simulations without the ponderomotive force show that this shock feature is only reproduced when accounting for the ponderomotive force.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗