Smooth Particle Hydrodynamic Code Predictions for Meteoroid Damage to Thermal Protection Systems Shielded By Composite Structures
Explore the source record for details and available documents.
Brooke Corbett↗
SEARCH · Search NASA
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.
Laser-Plasma Interaction Modeling in LLNL Radiation-Hydrodynamic Codes
Explore the source record for details and available documents.
Physics - Plasma physics↗
High-Fidelity Simulations of Objects in Accretion Disks by Coupling REBOUND to a Hydrodynamics Code
Explore the source record for details and available documents.
Dempsey, Adam Mahlon↗
Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives
A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.
97 MATHEMATICS AND COMPUTING↗
Development of an Enhanced Radiation Physics Toolset for Modeling Spectral and Imaging Signatures in the Warm Dense Matter Experiments
Radiative and atomic processes in plasmas play a critical role in a wide variety of high energy density laboratory plasma (HEDLP) experiments. The emission, absorption, and transport of radiation can strongly affect the overall energetics and evolution of such plasmas. In addition, radiation-based diagnostics – including imaging, spectroscopy, and absolute flux measurements – are widely used to determine key features of HEDLPs. To advance our understanding of HEDLP science, it is vital to have high-fidelity computational physics tools that have well-tested radiation physics modeling, and that are readily accessible to researchers in the HEDLP community. Simulations play an extremely important role for planning and designing the experiments, as well as for post-experiment data analysis. Prism Computational Sciences develops software that is used by National Laboratories and universities (including five members of LaserNetUS network). Prominent examples of such research efforts include z-pinch and short-pulse laser experiments designed to study the basic physics of photoionized plasmas and photoionization fronts, as well as their application to astrophysical plasmas. The main effort was dedicated to the development of non-equilibrium equation-of-state (EOS) models within the HELIOS-CR code, a hydrodynamics code with inline collisional-radiative atomic kinetics. Gas cell experiments on Z and Omega demonstrated the importance of non-equilibrium effects on atomic kinetics in photoionized plasmas. Recent proof-of-principle experiments on Omega EP confirmed the advantages of using a short-pulse laser to create an intense radiation drive, leading to additional experiments being proposed. Photoionization front experiments at LLE also emphasizes the importance of radiation and atomic physics. In both studies, HELIOS-CR simulations played a crucial role in computing non-equilibrium opacities and ionization distributions. A newly developed non-LTE EOS model will help addressing possible non-equilibrium effects, on for example specific heat, and their influence in plasma evolution. Prism also implemented support for open-source atomic data generated by the Flexible Atomic Code. This allows researchers to generate custom atomic tables and use them within the complex framework of simulation tools developed by Prism. The ability to use open-source atomic data would be extremely valuable for hydrodynamics and spectroscopic simulations that include high-Z materials, e.g., picosecond x-ray pulse generation experiments. Support for new atomic structures was fully implemented, and the data can be used by all simulation tools developed at Prism: radiation-hydrodynamics, imaging and spectroscopy, EOS and opacity. The development resulted in a significant fidelity enhancement to the simulations tools developed by Prism that are currently used in other cutting-edge experiments including: opacity measurement experiments performed to both understand the basic radiative and atomic properties of plasmas as well as provide data for more accurately modeling the internal structure of the Sun and other stars, high-intensity short-pulse laser experiments performed to develop short-wavelength light sources for use as backlighters and to investigate fast ignition concepts for inertial fusion energy; capsule implosion experiments designed to develop inertial fusion as an energy source, etc.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Comparison of SPHC Hydrocode Results with Penetration Equations and Results of Other Codes
The SPHC hydrodynamic code was used to simulate impacts of spherical aluminum projectiles on a single-wall aluminum plate and on a generic Whipple shield. Simulations were carried out in two and three dimensions. Projectile speeds ranged from 2 kilometers per second to 10 kilometers per second for the single-wall runs, and from 3 kilometers per second to 40 kilometers per second for the Whipple shield runs. Spallation limit results of the single-wall simulations are compared with predictions from five standard penetration equations, and are shown to fall comfortably within the envelope of these analytical relations. Ballistic limit results of the Whipple shield simulations are compared with results from the AUTODYN-2D and PAM-SHOCK-3D codes presented in a paper at the Hypervelocity Impact Symposium 2000 and the Christiansen formulation of 2003.
Evans, Steven W.↗
The special relativistic shock tube
The shock-tube problem has served as a popular test for numerical hydrodynamics codes. The development of relativistic hydrodynamics codes has created a need for a similar test problem in relativistic hydrodynamics. The analytical solution to the special relativistic shock-tube problem is presented here. The relativistic shock-jump conditions and rarefaction solution which make up the shock tube are derived. The Newtonian limit of the calculations is given throughout.
Thompson, Kevin W.↗
Coupling laser physics to radiation-hydrodynamics
In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.
47 OTHER INSTRUMENTATION↗
A hydrodynamic approach to cosmology - Methodology
The present study describes an accurate and efficient hydrodynamic code for evolving self-gravitating cosmological systems. The hydrodynamic code is a flux-based mesh code originally designed for engineering hydrodynamical applications. A variety of checks were performed which indicate that the resolution of the code is a few cells, providing accuracy for integral energy quantities in the present simulations of 1-3 percent over the whole runs. Six species (H I, H II, He I, He II, He III) are tracked separately, and relevant ionization and recombination processes, as well as line and continuum heating and cooling, are computed. The background radiation field is simultaneously determined in the range 1 eV to 100 keV, allowing for absorption, emission, and cosmological effects. It is shown how the inevitable numerical inaccuracies can be estimated and to some extent overcome.
Cen, Renyue↗
Hydrodynamic models of a cepheid atmosphere
A method for including the solution of the transfer equation in a standard Henyey type hydrodynamic code was developed. This modified Henyey method was used in an implicit hydrodynamic code to compute deep envelope models of a classical Cepheid with a period of 12(d) including radiative transfer effects in the optically thin zones. It was found that the velocity gradients in the atmosphere are not responsible for the large microturbulent velocities observed in Cepheids but may be responsible for the occurrence of supersonic microturbulence. It was found that the splitting of the cores of the strong lines is due to shock induced temperature inversions in the line forming region. The adopted light, color, and velocity curves were used to study three methods frequently used to determine the mean radii of Cepheids. It is concluded that an accuracy of 10% is possible only if high quality observations are used.
Karp, A. H.↗
TIA: A forward model and analyzer for Talbot interferometry experiments of dense plasmas
Interferometry is one of the most sensitive and successful diagnostic methods for plasmas. However, owing to the design of most common interferometric systems, the wavelengths of operation and, therefore, the range of densities and temperatures that can be probed are severely limited. Talbot–Lau interferometry offers the possibility of extending interferometry measurements to x-ray wavelengths by means of the Talbot effect. While there have been several proof-of-concept experiments showing the efficacy of this method, it is only recently that experiments to probe High Energy Density (HED) plasmas using Talbot–Lau interferometry are starting to take place. To improve these experimental designs, we present here the Talbot-Interferometry Analyzer (TIA) tool, a forward model for generating and postprocessing synthetic x-ray interferometry images from a Talbot–Lau interferometer. Although TIA can work with any two-dimensional hydrodynamic code to study plasma conditions as close to reality as possible, this software has been designed to work by default with output files from the hydrodynamic code FLASH, making the tool user-friendly and accessible to the general plasma physics community. Here, the model has been built into a standalone app, which can be installed by anyone with access to the MATLAB runtime installer and is available upon request to the authors.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Semi-analytic solutions to the Noh problem with a black box EoS
The objective of this paper is to derive a method of constructing semi-analytic solutions to the Noh problem when the equation of state is a black box. Such solutions can be used for verification tests of hydrodynamics codes. We present the underlying theory, the method for finding solutions, and several examples of derived semi-analytic solutions. We end by performing a classic verification convergence test comparing numerical results from a hydrodynamics code against a non-trivial semi-analytic solution.
97 MATHEMATICS AND COMPUTING↗
Validation of ray-based cross-beam energy transfer models
Ray-based cross-beam energy transfer (CBET) models have become a common feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion experiments. Here, these models are necessary for achieving better agreement with experimental measurements, but their detailed implementation can vary widely between the codes and often rely on artificial multipliers. To address this, a series of 2D and 3D test cases has been developed with validated solutions from wave-based calculations. Comparisons of various ray-based CBET models to the wave-based calculations highlight the essential physics that is required for accurate ray-based CBET modeling. Quantitative comparison metrics and/or field data from the wave-based calculations have been made available for use in the validation of other ray-based CBET codes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
An accelerated approach to inline non-LTE modeling
A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications
Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Thermonuclear 28 P(p, γ ) 29 S reaction rate and astrophysical implication in ONe nova explosion
An accurate 28 P(p, γ) 29 S reaction rate is crucial to defining the nucleosynthesis products of explosive hydrogen burning in ONe novae. Using the recently released nuclear mass of 29 S, together with a shell model and a direct capture calculation, we reanalyzed the 28 P(p, γ) 29 S thermonuclear reaction rate and its astrophysical implication. We focus on improving the astrophysical rate for 28 P(p, γ) 29 S based on the newest nuclear mass data. Our goal is to explore the impact of the new rate and associated uncertainties on the nova nucleosynthesis. We evaluated this reaction rate via the sum of the isolated resonance contribution instead of the previously used Hauser-Feshbach statistical model. The corresponding rate uncertainty at different energies was derived using a Monte Carlo method. Nova nucleosynthesis is computed with the 1D hydrodynamic code SHIVA. The contribution from the capture on the first excited state at 105.64 keV in 28 P is taken into account for the first time. We find that the capture rate on the first excited state in 28 P is up to more than 12 times larger than the ground-state capture rate in the temperature region of 2.5 × 10 7 K to 4 × 10 8 K, resulting in the total 28 P(p, γ) 29 S reaction rate being enhanced by a factor of up to 1.4 at ~1 × 10 9 K. In addition, the rate uncertainty has been quantified for the first time. It is found that the new rate is smaller than the previous statistical model rates, but it still agrees with them within uncertainties for nova temperatures. The statistical model appears to be roughly valid for the rate estimation of this reaction in the nova nucleosynthesis scenario. Using the 1D hydrodynamic code SHIVA, we performed the nucleosynthesis calculations in a nova explosion to investigate the impact of the new rates of 28 P(p, γ) 29 S. Our calculations show that the nova abundance pattern is only marginally affected if we use our new rates with respect to the same simulations but statistical model rates. Finally, the isotopes whose abundance is most influenced by the present 28 P(p, γ) 29 S uncertainty are 28 Si, 33,34 S, 35,37 Cl, and 36 Ar, with relative abundance changes at the level of only 3% to 4%.
Astronomy & Astrophysics↗
Physics of the Serrano Sn Transport Package
Serrano is a software package for simulating photon transport in moving material. It presently supports 1-D, unstructured 2-D for a wide variety of element types, and 3-D for tetrahedra and hexahedra. The package is designed to be a component of larger codes, such as radiation hydrodynamics codes. Only the radiation transport equations are solved by Serrano itself.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A novel, magnetically driven convergent Richtmyer–Meshkov platform
Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗