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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The BTSbot-nearby Discovery of SN 2024jlf: Rapid, Autonomous Follow-up Probes Interaction in an 18.5 Mpc Type IIP Supernova

We present observations of the Type IIP supernova (SN) SN 2024jlf, including spectroscopy beginning just 0.7 days (∼17 hr) after first light. Rapid follow-up was enabled by the new BTSbot-nearby program, which involves autonomously triggering target-of-opportunity requests for new transients in Zwicky Transient Facility data that are coincident with nearby (D < 60 Mpc) galaxies and identified by the BTSbot machine learning model. Early photometry and nondetections shortly prior to first light show that SN 2024jlf initially brightened by >4 mag day −1 , quicker than ∼90% of Type II SNe. Early spectra reveal weak flash ionization features: narrow, short-lived (1.3 < τ[days] < 1.8) emission lines of Hα, He II , and C IV . Assuming a wind velocity of v w = 50 km s −1 , these properties indicate that the red supergiant progenitor exhibited enhanced mass loss in the last year before explosion. We constrain the mass-loss rate to $1{0}^{-4}\lt \dot{M}\,[{M}_{\odot }\,{\mathrm{yr}}^{-1}]\lt 1{0}^{-3}$ by matching observations to model grids from two independent radiative hydrodynamics codes. BTSbot-nearby automation minimizes spectroscopic follow-up latency, enabling the observation of ephemeral early-time phenomena exhibited by transients.

core-collapse supernovae↗
DOE OSTI2025

The R -process Alliance: Enrichment of r -process Elements in a Simulated Milky Way–like Galaxy

We study the formation of stars with varying amounts of heavy elements synthesized by the rapid neutron-capture process (r-process) based on our detailed cosmological zoom-in simulation of a Milky Way–like galaxy with an N-body/smoothed particle hydrodynamics code, asura. Most stars with no overabundance in r-process elements, as well as the strongly r-process-enhanced (RPE) r-II stars ([Eu/Fe] > +0.7), are formed in dwarf galaxies accreted by the Milky Way within the 6 Gyr after the Big Bang. In contrast, over half of the moderately enhanced r-I stars (+0.3 < [Eu/Fe] ≤ +0.7) are formed in the main in situ disk after 6 Gyr. Our results suggest that the fraction of r-I and r-II stars formed in disrupted dwarf galaxies is larger the higher their [Eu/Fe] is. Accordingly, the most strongly enhanced r-III stars ([Eu/Fe] > +2.0) are formed in accreted components. These results suggest that non-r-process-enhanced stars and r-II stars are mainly formed in low-mass dwarf galaxies that hosted either none or a single neutron star merger, while the r-I stars tend to form in the well-mixed in situ disk. We compare our findings with high-resolution spectroscopic observations of RPE metal-poor stars in the halo and dwarf galaxies, including those collected by the R-Process Alliance. We conclude that observed [Eu/Fe] and [Eu/Mg] ratios can be employed in chemical tagging of the Milky Way’s accretion history.

Chemical enrichment↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2026

The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts

When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.

Astronomy and AstroPhysics↗
DOE OSTI2021

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

(U) Radiation Hydrodynamics in the Ristra Project’s Symphony code: Initial release

One of the recent focuses in the next-generation code (NGC/Ristra)project at LANL is to develop a cell-centered Lagrangian radiation-hydrodynamics code under FleCSI framework. The FleCSI frame-work provides a unified data and execution model in order to insulate production code development from the uncertainty of evolving runtimeand architectures. We have developed the code, Symphony, which in-corporates the hydrodynamics code, FleCSALE, and the radiation transport code, Puno.This memo is intended to describe the coupling of radiation and hydrodynamics within the Ristra Project’s code Symphony.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

New charged-particle transport computational capability: the SIT code An L-4 milestone

We have developed a new high-fidelity code for direct transport of charged particles using the simple integral transport method. The code can be coupled to the outputs of any hydrodynamical simulation code in 1-D, 2-D or 3-D. In this report we summarize the formalism involved in treating complex transport problems. We present physical examples wherein we have used the code to calculate the transport of alpha particles. Future work is planned to study the sensitivity of hydrodynamical mix to charged-particle radiochemistry and reaction-in-flight neutrons for the complex inertial confinement fusion problems encountered at NIF and at the Z-machine.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗
DOE OSTI2020

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement.The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response.With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal.In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

DIII-D↗
DOE OSTI2023

SPH modeling of biomass granular flow: Theoretical implementation and experimental validation

The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.

09 BIOMASS FUELS↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗