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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2025

Comparison of tungsten versus molybdenum for double shell capsules using machine learning design optimization

Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Simulations of high-velocity impacts on metal in preparation for the Psyche mission

In 2026, the NASA Discovery mission Psyche will orbit the asteroid (16) Psyche, the largest known metal-rich asteroid in the main belt. To estimate relative ages of the surface, identify re-surfacing events and better constrain Psyche's history, impact craters will be counted and characterized. No spacecraft has ever visited a metal-rich small body; therefore, laboratory-scale impact experiments and numerical simulations will play an important role in the interpretation of the mission's data. However, the planetary applications of high-velocity impacts have so far mostly been studied for silicate targets. Limited attention has been drawn to planetary objects predominantly made of metal, and more laboratory experiments and numerical calibrations are needed. As part of this effort, we present a suite of numerical simulations using an adaptative smoothed particles hydrodynamics numerical code (ASPH) reproducing a high-velocity impact experiment conducted on a steel target. This work primarily focuses on the influence of the chosen equation of state and initial distribution of flaws in the material on the estimated crater dimensions, damage and temperature. Additionally, we find that changing the EOS and initial flaw distribution affects the crater dimensions, though for a vast majority of cases the dimensions remain within 20% of the experimental values. The target is in most cases only locally weakened but not fully damaged, independently from the EOS chosen. Finally, temperatures at the impact point and around the forming crater can reach values above the melting point of iron at <100 GPa, which is in agreement with experimental observations. These results allow us to speculate on the differences expected between the surfaces of visited silicate-rich asteroids and that of the metal-rich target of the Psyche mission.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Probing ion species separation and ion thermal decoupling in shock-driven implosions using multiple nuclear reaction histories

Simultaneously measured DD, DT, and D 3 He reaction histories are used to probe the impacts of multi-ion physics during the shock phase of inertial confinement fusion implosions. In these relatively hydrodynamiclike (burn-averaged Knudsen number &lang;N K &rang; ~0.3) shock-driven implosions, average-ion hydrodynamic DUED simulations are able to reasonably match burnwidths, nuclear yields, and ion temperatures. However, kinetic-ion FPION simulations are able to better simulate the timing differences and time-resolved reaction rate ratios between DD, DT, and D 3 He reactions. FPION simulations suggest that the D 3 He/DT reaction rate ratio is most directly impacted by ion species separation between the 3 He and T ions, whereas the D 3 He/DD reaction rate ratio is affected by both ion species separation and ion temperature decoupling effects.

Hydrodynamic codes↗

DOE OSTI2020

Effect of cross-beam energy transfer on target-offset asymmetry in direct-drive inertial confinement fusion implosions

The unintentional mispositioning of inertial confinement fusion (ICF) capsules from the center of laser beam convergence has long been shown in simulations to generate large ℓ = 1 asymmetry and significantly degrade implosion symmetry and fusion yields. Experimental yields on the OMEGA Laser System, however, have shown much less sensitivity to this initial target offset. This paper presents simulations of offset ICF implosions improved by including a physics model of cross-beam energy transfer (CBET), a mechanism of laser energy scattering from one beam to another. Room-temperature OMEGA implosion experiments with prescribed target offsets are simulated with and without CBET, illustrating that CBET mitigates the ℓ = 1 implosion asymmetry from target offset. Here, comparison of simulations to multiple complementary experimental observables indicates the addition of CBET physics in offset simulations is necessary to match experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Laser transport and backscatter in low-density SiO 2 and Ta2O 5 foams

Experiments using a single 527 nm wavelength beam interacting with sub- and supercritical density SiO 2 and Ta 2 O 5 foams examined laser propagation and backscatter from laser–plasma instabilities such as Stimulated Brillouin Scattering (SBS). Two densities of each material were examined, and multiple diagnostics were used to characterize the propagation and backscatter. For 5 mg/cc SiO 2 (n e /n c = 0.375), the laser propagation distance was well approximated by treating the foam as a gas. However, for the 2 mg/cc SiO 2 foam (n e /n c = 0.15), the same model over-predicts the propagation distance by ~40%. Furthermore, existing analytical theories on propagation through subcritical foams were able to account for this difference. The laser heat wave propagated ~1/2 as far in Ta 2 O 5 than SiO 2 foams with similar electron density. We showed that this difference is due to the increased radiation losses in the higher Z foam. The fraction of backscattered light scales linearly with incident laser intensity for the range of intensities examined. Ta 2 O 5 foams had significantly lower levels of backscatter (1–3%) than the SiO 2 (4–8%), which is consistent with estimates of large Landau damping due to the presence of the oxygen atoms. The measured fraction of SBS backscattered laser energy for a 2 mg/cc SiO 2 foam shot was ~4 times lower than predicted by simulations assuming a gas-like foam. We found that we needed to assume increased ion heating such that T i /T e ~1.2–1.5 in the plasma to agree with the measured SBS reflectivity. Analytical models of laser-heated foams predict preferential heating of the ions as has been observed in previous experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Understanding and controlling capsule symmetry in near vacuum hohlraums at the National Ignition Facility

The near vacuum hohlraum platform is an inertial confinement fusion design at the National Ignition Facility (NIF) that uses the lowest practical density of helium gas of 30 μg/cc to fill the hohlraum, which is ten times lower than now used routinely. This has several advantages, such as high laser coupling; however, the inability to understand and simulate the symmetry of the imploded capsule has limited the use of this platform. This work presents the first simulations that are able to accurately capture the highly prolate implosion seen experimentally without unphysical, ad hoc model changes. While previous investigations attributed this asymmetry to multi-species interpenetration in the hohlraum, we find that this alone has little effect on symmetry. Instead, it is the presence of crossed-beam energy transfer (CBET), occurring with no applied wavelength shift between the laser beams, that increases the laser power to the inner cones and causes a more prolate implosion. The effect of CBET is increased in the simulation model when the hohlraum laser entrance hole hardware is included. Using this understanding, CBET is exploited by shifting the inner-beam wavelength by –0.75 Å (at 1ω) with respect to the outer-beams. This transfers laser power to the outer-beams in contrast to positive wavelength shifts as done routinely on NIF and produces a round capsule implosion in our simulations. Furthermore, this work shows the possibility of the near vacuum hohlraum as a viable experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Hybrid direct drive with a two-sided ultraviolet laser

This paper presents a “hybrid” approach to direct drive inertial confinement fusion that can exploit a high-energy gas laser with two opposed beams. The target and driver are asymmetric, much like experiments performed on the National Ignition Facility, but have been designed to benefit from scale and their particular compatibility with a fusion power plant. The imploded masses (and areal densities) are increased by a factor of 12 (3) and provide a path to high-gain implosions that robustly ignite. The design also mitigates common concerns such as laser imprint and cross-beam energy transfer. We discuss the rationales for a hybrid target, the methods used to control implosion symmetry, and the implication(s) for inertial fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Dynamics of implosion, bounce, and re-implosion of a shell directly driven by a time-shaped laser pulse

Gas-filled plastic spherical shells have been imploded using the 60-beam OMEGA laser. The same laser pulse was used as in a previous experiment, which demonstrated the dynamic formation of a shell by irradiating a homogeneous plastic foam sphere. In the present experiment, the shell initially implodes, bounces, and expands; a new shell then forms and undergoes a second implosion. Analysis of x-ray self-emission-streaked images and sequences of 2D x-ray self-emission framing camera images show shell evolution that is in reasonable agreement with 1D simulations. Shell integrity and symmetry are well maintained until the formation of the new shell. The final implosion, instead, is affected by substantial asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗