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

Using tabulated NLTE data for Hohlraum simulations

Non-local thermodyamic equilibrium (NLTE) atomic kinetics is necessary in inertial confinement fusion Hohlraum simulations for adequately modeling high-Z walls and dopants but is computationally very expensive. We present here an approach for tabulating NLTE material data in an economical manner. Material properties for an arbitrary radiation field are provided through tabulated data for a limited set of radiation fields plus derivatives with respect to components of those radiation fields. We have implemented this method in a radiation-hydrodynamics code and compare Hohlraum simulations done with inline NLTE calculations with those done with tables constructed from the same atomic data. The results demonstrate that NLTE tables can replace inline calculations in Hohlraum simulations without a significant loss of accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

New Yield Estimates for Nuclear Detonations Over Water

One of the primary methods used by the United States to estimate the yields of atmospheric nuclear tests performed in the 1940’s, 50’s, and 60’s was to film each detonation with a high-speed camera, measure the blast radius as a function of time, and then use Taylor’s spherical blast wave equation to calculate the source energy (i.e., the yield). For surface detonations, the yield in Taylor’s equation was assumed to be twice as large as the actual yield to account for the energy reflected from the surface, assuming the ground was a perfect reflector. Recent simulations of surface detonations over water (referred to as barge shots) have shown that Taylor’s equation must be modified further to account for material entrainment and heat losses. We have computed this lost energy using a high-order radiation-hydrodynamics code to determine the reduction in blast radius vs. time. By analyzing simulation data in the same manner as film data, we have determined that a yield correction factor of 1.27± 0.04 should be applied to the barge shots.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1980

Curvilinear grids for sinuous river channels

In order to effectively analyze the flow in sinuous river channels, a curvilinear grid system was developed for use in the appropriate hydrodynamic code. The CENTERLINE program was designed to generate a two dimensional grid for this purpose. The Cartesian coordinates of a series of points along the boundaries of the sinuous channel represent the primary input to CENTERLINE. The program calculates the location of the river centerline, the distance downstream along the centerline, and both radius of curvature and channel width as a function of such distance downstream. These parameters form the basis for the generation of the curvilinear grid. Based on input values for longitudinal and lateral grid spacing, the corresponding grid system is generated and a file is created containing the appropriate parameters for use in the associated explicit finite difference hydrodynamic programs. Because of the option for a nonuniform grid, grid spacing can be concentrated in areas containing the largest flow gradients.

Tatom, F. B.↗
DOE OSTI2025

Cosmic recombination in the presence of primordial magnetic fields

Primordial magnetic fields (PMFs) may explain observations of magnetic fields on extragalactic scales. They are most cleanly constrained by measurements of cosmic microwave background radiation (CMB) anisotropies. Their effects on cosmic recombination may even be at the heart of the resolution of the Hubble tension. We present the most detailed analysis of the effects of PMFs on cosmic recombination to date. To this end we extend the public magneto-hydrodynamic code ENZO with a new cosmic recombination routine, Monte-Carlo simulations of Lyman-α photon transport, and a Compton drag term in the baryon momentum equation. The resulting code allows us, for the first time, to realistically predict the impact of PMFs on the cosmic ionization history and the clumping of baryons during cosmic recombination. Our results identify the importance of mixing of Lyman-α photons between overdense- and underdense- regions for small PMF strength. This mixing speeds up recombination beyond the speed-up due to clumping. We also investigate the effects of pecuilar flows on the recombination rate and find it to be small for small PMF strengths. For non-helical PMFs with a Batchelor spectrum we find a surprising dependency of results on ultra-violet magnetic modes. We further show that the increase in the ionization fraction at low redshift by hydrodynamic baryon heating due to PMF dissipation is completely compensated by the faster recombination from baryon clumping. In conclusion, the present study shall serve as a theoretical foundation for a future precise comparison of recombination with PMFs to CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

Context. Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims. Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods. We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results. The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

hydrodynamics↗
DOE OSTI2021

Unabsorbed light beamlets for diagnosing coronal density profiles and absorption nonuniformity in direct-drive implosions on OMEGA

Laser light scattered from a target is the most-direct measurement for diagnosing laser absorption in a direct-drive implosion. The 3ωGOI (gated optical imager) beamlets diagnostic images unabsorbed light from all 60 OMEGA beams as distinct “beamlet” spots for each beam. The implosion can be diagnosed using the position and intensity of these beamlets. The position of each beamlet in the image is determined by refraction and can be used to fit the coronal plasma density profile of the implosion. The inferred plasma density profiles are comparable to profiles predicted by the 1-D hydrodynamics code LILAC but suggest the predictions underestimate the density farther out in the corona. The intensity of light in each spot depends on the cumulative effects of absorption and cross-beam energy transfer along the beamlet’s path through the corona. Furthermore, the measured variation in intensity and polarization between similar spots indicates that absorption during OMEGA implosions is less uniform than previously known.

47 OTHER INSTRUMENTATION↗
NASA NTRS1992

Effects of convection electric field on upwelling and escape of ionospheric O(+)

A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.

Cladis, J. B.↗
DOE OSTI2020

FLOP for FLAG Output Plotting

The Los Alamos hydrodynamics code, FLAG, is capable of dumping many types of output for many different variables in an array of formats. While some outputs are best viewed in a multidimensional engineering analysis visualization software, others are best viewed as 1-D “this versus that” curves. During development of a FLAG model, it is frequently required to quickly assess a model’s performance by reviewing such curves, and doing so may involve writing scripts repeatedly, adapting them each time to a specific model’s parameters. This report describes an application developed specifically to improve user efficiency in reviewing 1-D curve dumps from FLAG.

42 ENGINEERING↗
DOE OSTI2022

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Progenitors of Superluminous Type Ia Supernovae

Abstract Recent observations of type Ia supernovae (SNe Ia) have discovered a subclass of “super-Chandrasekhar” SNe Ia (SC SNe Ia) whose high luminosities and low ejecta velocities suggest that they originate from the explosions of white dwarfs (WDs) with masses that exceed the Chandrasekhar mass limit. Different models have been proposed to explain the progenitors of these explosions, including a “magnetized WD” model and a “WD merger” model. To test the robustness of these models, we conduct a 1D numerical parameter survey of WD explosions using these models as initial conditions. We follow the explosions using the hydrodynamics code Castro and then use the radiation transport code SuperNu to create light curves and spectra for the models. We find that while both classes of models fall within the range of SC SNe Ia observations on the light-curve width–luminosity relation, only the WD merger models reproduce the observed low ejecta velocities. The light curves of our merger models are more similar photometrically to observations than our magnetized models. Given this, we discuss possible explanations for the brightest SC SNe Ia observations that cannot be reproduced with our WD merger models. This study provides the basis for future SC SNe Ia observations and higher-dimensional numerical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling Turbulence in Galactic Centers

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way’s Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dust Survival in Galactic Winds

This repository contains three-dimensional volumetric data from an Eulerian hydrodynamical simulation (conducted on a uniform Cartesian grid) generated by the Cholla hydrodynamics code. The datasets contain snapshots (full-grid, projections, and slices) in the HDF5 format of a multi-phase medium in which a hot, diffuse, dust-free background wind accelerates a cool, dense cloud of gas and dust. This scenario is intended to represent a supernova-driven galactic outflow, in which hot supernova winds are thought to accelerate cool interstellar medium material out of the galactic disk into the surrounding circumgalactic medium. There are three separate datasets for simulations corresponding to three cloud evolutionary scenarios: long-term cloud survival (surv), marginal cloud survival (disr), and cloud destruction (dest). Projection and slice images of the simulations are also included in this repository.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
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 OSTI2023

On the importance of three-dimensional modeling for high-energy-density physics experiments

Laser-driven cylindrical implosion experiments enable direct measurements of hydrodynamic instability growth in convergent geometries, providing a wealth of validation data in the high-energy-density regime. These experiments are designed to be nearly axially invariant, allowing for modeling with complementary two-dimensional slices of the cylinder. Two distinct hydrodynamics codes are employed to model a subset of these experiments, and the results are shown to be in very good agreement with each other and the available experimental data. While this 2D modeling approach adequately captures most of the physics of the implosion and ensuing instability growth, there are crucial aspects from the three-dimensional nature of the experiments that are missed in 2D. The first fully 3D simulations of these experiments are presented, and small but significant differences are found to arise from both the axial and azimuthal non-uniformity in the laser drive. Recent experimental results confirming the drive asymmetry are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Staged Z-pinch radiation-hydrodynamic simulations on a 20-MA driver

We report on one-dimensional simulations of a staged Z-pinch (SZP) fusion scheme using the radiation-hydrodynamics code HYDRA [Marinak et al., Phys. Plasmas 8, 2275 (2001)]. A pulsed-power driver, based on the Z-machine at Sandia National Laboratories, with a peak current around 20 MA, is used to implode a high-Z gas liner surrounding a cylindrical DT fuel target. For the best conditions used in this work, excluding some physical processes, such as mix, we find a “clean” fusion yield of 7.05 MJ (2.5 ×1018 neutrons), thus showing the potential for high yields using this scheme. The design achieving this yield was determined through a study of different liner materials, liner masses, and fuel masses. We also address previous simulations of the SZP setup that used ad hoc limits of the total internal energy of the liner that can artificially increase fusion yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗