Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗
DOE OSTI2020

High-performance computing in water resources hydrodynamics

In this work, we present a vision of future water resources hydrodynamics codes that can fully utilize the strengths of modern high-performance computing. The advances to computing power, formerly driven by the improvement of central processing unit processors, now focus on parallel computing and, in particular, the use of graphics processing units (GPUs). However, this shift to a parallel framework requires refactoring the code to make efficient use of the data as well as changing even the nature of the algorithm that solves the system of equations. These concepts along with other features such as the precision for the computations, dry regions management, and input/output data are analyzed in this paper. A 2D multi-GPU flood code applied to a large-scale test case is used to corroborate our statements and ascertain the new challenges for the next-generation parallel water resources codes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A boundary condition for Guderley’s converging shock problem

The Guderley model of a self-similar imploding shock based on the group invariance of the flow equations is a powerful tool in understanding the behavior of converging shock waves. Two modifications described here improve the predictions of observable quantities in spherical-shock wave experiments. First, a non-infinite boundary condition is established by the isentropic release of the outer pressure. Second, a two-temperature system of ions and electrons allows description of higher temperatures while conserving energy and without perturbing the overall hydrodynamics of the solution. Furthermore, these modifications of the Guderley model improve the prediction of the observables in laser driven spherical shock experiments in reference to a one dimensional (1-D) hydrodynamics code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

In this work, we report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for magneto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six vertices and the center of a hexagon covering approximately one-tenth of the surface area of a 9-ft-diameter spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, and Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging, leading to shock formation rather than potentially more desirable shockless jet merging; and (3) the liner Mach number can remain high (≳10), as required for plasma liners to be an effective compression driver. Secondly, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet mass balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. Ultimately, for both experiments, we present comparisons between experimental and synthetic data from three-dimensional hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Laser material interactions in tamped materials on picosecond time scales in aluminum

Here, a 100 ps laser is used to probe the pressure generation, depth of the non-solid ablator, and the non-linear optical effects through tamper materials. Samples consisted of an aluminum ablator with tampers of sapphire and coverslip glass. In general, the sapphire tamped sample achieves higher pressures at lower laser intensities as compared to the coverslip glass tamped sample. Attempts to model the details of this set of experimental data with standard available radiation coupled hydrodynamic codes make clear that more physics is needed in these simulations to accurately predict the impact of the tamper material on the pressure generation and the depth of non-solid aluminum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development of uncertainty-aware equation-of-state models: Application to copper

Sophisticated hydrodynamic codes are commonly used to understand and predict events relevant to natural and applied sciences. The degree to which these simulations reflect reality, however, is dependent on how well we understand the materials and underlying physics involved. These research communities need material models that communicate the uncertainty in the physical properties, which at their basest form comes from the uncertainty in the underlying experimental measurements. Additionally, we have constructed a new framework for using both experimental measurements and the associated experimental uncertainties to build equation-of-state models that reflect not only current best measurements but also the accuracy of those measurements. This method had been used to construct an ensemble of equation-of-state models for copper that communicates the experimental uncertainties in the data through the equation-of-state model, which is available for application in any simulation metric of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Stationary Bragg reflection of laser light in inhomogeneous absorbing plasmas inside inertial confinement fusion Hohlraums

Laser-produced plasma in inertial confinement fusion (ICF) Hohlraums are marked with density non-uniformity whose length scale can go down to micrometers. This scale is of the order of the laser wavelength. The WKB approximation, which is classically used in radiation-hydrodynamic codes to compute the laser trajectory, cannot correctly take into account such small-scale inhomogeneity of the plasma. Going beyond this approximation, we predict a novel mechanism for the laser reflection. Here we show that an electromagnetic plane wave with wave number k resonates with the k B =2 k Fourier component of a multimode perturbation of the background density and generates a reflected wave. It is the first time that this reflection is considered for stationary inhomogeneous ICF plasmas, and the energy absorption is taken into account. This mechanism, which is a form of Bragg reflection, can occur away from the critical surface and generate a drift of the location of the laser absorption. Furthermore, this absorption will be periodically modulated with a k B wave number. The stationary Bragg reflection can explain ongoing discrepancies between experimental and numerical data about laser trajectory and absorption in ICF Hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗
DOE OSTI2024

A review on charged-particle transport modeling for laser direct-drive fusion

Inertial confinement fusion (ICF) with the laser-indirect-drive scheme has recently made a tremendous breakthrough recently after decades of intensive research effort. Taking this success to the next step, the ICF community is coming to a general consensus that laser direct-drive (LDD) fusion might be the viable way for enabling inertial fusion energy (IFE) and high-gain targets for other applications. Designing and understanding LDD fusion targets heavily rely on radiation-hydrodynamic code simulations, in which charged-particle transport plays an essential role in modeling laser-target energy coupling and bootstrap heating of fusion-produced α-particles. To better simulate charged-particle transport in LDD targets, over the past four decades the plasma physics community has advanced transport calculations from simple plasma physics models to sophisticated computations based on first-principles methods. In this review, we give an overview of the current status of charged-particle transport modeling for LDD fusion, including what challenges we still face and the possible paths moving forward to advance transport modeling for ICF simulations. We hope this review will provide a summary of exciting challenges to stimulate young minds to enter the field, facilitate further progress in understanding warm-dense matter physics, and ultimately bridge toward the success of reliable LDD fusion designs for IFE and other high-gain ICF applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hot-electron generation at direct-drive ignition-relevant plasma conditions at the National Ignition Facility

Laser–plasma interaction instabilities can be detrimental for direct-drive inertial confinement fusion by generating high-energy electrons that preheat the target. Moreover, an experimental platform has been developed and fielded on the National Ignition Facility to investigate hot-electron production from laser–plasma instabilities at direct-drive ignition-relevant conditions. The radiation-hydrodynamic code DRACO has been used to design planar-target experiments that generate plasma and interaction conditions comparable to direct-drive ignition designs: I L ~10 15 W/cm 2 , T e > 3 keV, and density-gradient scale lengths of L n ~ 600 μ m in the quarter-critical density region. The hot-electron properties were inferred by comparing the experimentally observed hard x-ray spectra to Monte Carlo simulations of hard x-ray emission from hot electrons depositing energy in the target. Hot-electron temperatures of ~40 keV to 60 keV and the fraction of laser energy converted to hot electrons of ~0.5% to 5% were inferred in plastic targets for laser intensities at the quarter-critical density surface of (~4 to 14) × 10 14 W/cm 2 . The use of silicon ablators was found to mitigate the hot-electron preheat by increasing the threshold laser intensity for hot-electron generation from ~3.5 × 10 14 W/cm 2 in plastic to ~6 × 10 14 W/cm 2 in silicon. The overall hot-electron production is further reduced in silicon ablators when the intensity threshold is exceeded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reaction-in-flight neutrons as a diagnostic for hydrodynamical mixing in double shell inertial confinement fusion capsules

We examine reaction-in-flight (RIF) neutrons as diagnostics for hydrodynamical mixing of high-Z shell material into the hotspot of double shell capsules that are designed for the National Ignition Facility. In particular, we consider the effects of different levels of mixing of tungsten shell material into the DT gas on RIF spectra. Using a set of 1D simulations from the radiation hydrodynamic code xRAGE to determine the temperature and density profiles of the mixed W-DT gas, we find that increasing the mass of mix systematically reduces the ratio of RIF neutrons to primary 14 MeV neutrons. The shape of the RIF spectrum also changes with mix, with the predicted spectrum softening in energy with increasing mix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The impact of varying inhomogeneous reionization histories on metrics of Ly α opacity

The epoch of hydrogen reionization is complete by z = 5⁠, but its progression at higher redshifts is uncertain. Measurements of Ly α forest opacity show large scatter at z < 6⁠, suggestive of spatial fluctuations in neutral fraction, temperature, or ionizing background, either individually or in combination. However, there are degeneracies in the impact of such fluctuations, necessitating careful modelling. We develop a framework for modelling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. We generate several reionization histories using seminumerical code AMBER, and implement them in the Nyx cosmological hydrodynamics code to examine the impact on the evolution of gas within the simulation and the associated metrics of the Ly α forest opacity. We find that the pressure smoothing scale within the intergalactic medium is strongly correlated with the adiabatic index of the temperature–density relation. We find that while models with 20 000 K photoheating at reionization are better able to reproduce the shape of the observed z = 5 1D flux power spectrum than colder ones, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at z < 6.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Radiative transfer in the dynamic atmospheres of long period variable stars

An iterative procedure is presented for determining the thermal structure and dynamics of Mira-type stellar atmospheres, where the non-LTE radiative transfer code PANDORA is used in conjunction with the Bowen hydrodynamics code of Iowa State University. Preliminary results are reported for an atmospheric model of a pulsating AGB star of 1 solar mass, 240 solar radii, Teff = 3000 K, and a period of 320 days. At the present time, H, H(-), Mg I, and Mg II radiative transfer calculations have been completed and synthetic spectra are shown for H-alpha. The radiative transfer calculations demonstrate that cooling in the innermost shock of the original Bowen model is underestimated due to the omission of various hydrogen transitions. These initial results suggest that the main shock of the Bowen models are too hot and/or too deep.

Luttermoser, Donald G.↗