Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrocode”

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 19 records

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗

Hydrocode predictions of collisional outcomes: Effects of target size

Traditionally, laboratory impact experiments, designed to simulate asteroid collisions, attempted to establish a predictive capability for collisional outcomes given a particular set of initial conditions. Unfortunately, laboratory experiments are restricted to using targets considerably smaller than the modelled objects. It is therefore necessary to develop some methodology for extrapolating the extensive experimental results to the size regime of interest. Results are reported obtained through the use of two dimensional hydrocode based on 2-D SALE and modified to include strength effects and the fragmentation equations. The hydrocode was tested by comparing its predictions for post-impact fragment size distributions to those observed in laboratory impact experiments.

Ryan, Eileen V.↗

Dynamic fragmentation in impacts - Hydrocode simulation of laboratory impacts

The dynamic fragmentation in impacts into solids is examined using a physical model for the formation and growth of cracks in rocks. The physical model is then inserted into a numerical model (hydrocode) of stress wave propagation and interaction, from which the outcome of a given impact event can be computed. The hydrocode model predicts fragment sizes due to impact in terms of shock waves propagating in a homogeneous elastic medium containing a distribution of crack nucleation centers known as Weibull flaws.

Melosh, H. J.↗

A Novel Porosity Model for Use in Hydrocode Simulations

Introduction: Numerical modeling of impact cratering has reached a high degree of sophistication; however, the treatment of porous materials still poses a large problem in hydrocode calculations. Porosity plays only a minor role in the formation of large craters on most planetary objects, but impacts on comets are believed to be highly affected by the presence of porosity, which may be as much as 80%. The upcoming Deep Impact Mission (launched January 2005) will provide more detailed data about the composition of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of its nucleus. The numerical simulations of such impact events requires an appropriate model for how pore space in the comet is crushed out during the violent initial stage of the impact event. Most hydro-codes compute the pressure explicitly using an "equation of state" (EOS) for each material, which relates changes in density and internal energy to changes in pressure. The added complication introduced by porosity is that changes in a material s density are due to both the closing of pore space (compaction) and compression of the matrix. The amount of resistance to volume change and the amount of irreversible work done during these two processes is very different; it is far easier to compact a porous material sample than to compress a non-porous sample of the same material. As an alternative to existing porosity models, like the Pdot(alpha) model [1], we present a novel approach for dealing with the compaction of porosity in hydrocode calculations.

Wuennemann, K.↗

Verification and Validation of the RAGE Hydrocode in Preparation for Investigation of Impacts into a Volatile-rich Target

Before a hydrocode is used to investigate a question of scientific interest, it should be tested against analogous laboratory experiments and problems with analytical solutions. The Radiation Adaptive Grid Eulerian (RAGE) hydrocode[1], developed by Los Alamos National Laboratory (LANL) and Science Applications International Corporation (SAIC)[2,3] has been subjected to many tests during its development.[4,5] We extend and review this work, emphasizing tests relevant to impact cratering into volatile-rich targets.

Plesko, C. S.↗

Comparison of Damage from Hydrocode Simulations of an Asteroid Airburst or Impact on Land, in Deep, or in Shallow Water

If an asteroid is discovered to be on a collision course with Earth and there is insufficient time for a deflection effort to make it miss Earth completely, should it be redirected to a land or ocean impact? While distance from densely populated areas should obviously be maximized, the differing ability of air blast, seismic waves, and tsunami waves to cause damage at distance does not make the choice between land and ocean impacts an immediately obvious one. More broadly this work is a step towards improving damage models from asteroid impacts. This extended abstract follows the hypothetical scenario of the 2017 IAA Planetary Defense Conference where a 100-250m diameter asteroid is on a potential impact course with Earth. A hydrocode was used to simulate impacts into the most sparsely populated areas along the eastern end of the hypothetical impact corridor- specifically in the Gobi Desert, in the shallow waters of the Sea of Japan, and in the deep waters of the Japan Trench in the Pacific Ocean.

Shallow Water↗

Exploring density and strength variations in asteroid 16 Psyche’s composition with 3D hydrocode modeling of its deepest impact structure

Asteroid 16 Psyche is the largest metallic Main Belt Asteroid and is the subject of a forthcoming NASA mission. The composition of Psyche is still unknown and subject of recent debate. In particular, how much porosity is within Psyche, along with how much of Psyche consists of non-metallic versus metallic materials, are central questions to the issue of Psyche’s composition. If Psyche is indeed predominantly composed of metallic materials, it would need to have considerable porosity (~ 30%–50%) for a composition consistent with its expected bulk density (~ 3.7–4.1 g/cm). In this work, we vary the density and strength of Psyche by including uniform and layered fields of pseudo-microporosity, in addition to investigating the presence of macroscopic voids, i.e., spaces larger than the size of the simulation’s mesh cells, in rubble-pile configurations. Further, all configurations result in bulk densities within the uncertainties of measured values, however the strength of Psyche and the distribution of pseudo-pores are varied. Through 3D computational models of Psyche’s deepest impact structure, we show that Psyche’s composition is unlikely to contain only pseudo-microporosity. Rather, rubble pile structures, which include macroscopic voids, are shown to match the crater’s measured aspect ratio better than simulations of structures that included only pseudo-microporosity.

3D↗

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗

Tethers as Debris: Hydrocode Simulation of Impacts of Tether Fragments on Planar Aerospace Materials

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas, and so tethers are quite susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, and hence have longer working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such encounters is important in assessing the risks to other spacecraft posed by tethers. In this paper I discuss the damage that two types of tethers may produce on planar aerospace materials, as revealed by hyper- velocity impact simulations using the SPHC hydrodynamic code. Tether types considered include a single nylon line and a complex design including metal wires. Target materials considered include the aluminum plates typically used in debris shielding, and solar panels.

Evans, Steven W.↗

Tethers as Debris: Hydrocode Simulation of Impacts of Polymer Tether Fragments on Aluminum Plates

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas. Consequently, tethers are very susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, in order to lengthen their working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris, to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such collisions is important in assessing the risks tethers pose to other spacecraft. This paper discusses the damage that polymer tethers may produce on aluminum plates, as revealed by hypervelocity impact simulations using the SPHC hydrodynamic code.

Evans, Steven W.↗

Hydrocode Simulations of the Chesapeake Bay Impact

The Chesapeake Bay Impact Crater (CBIC) formed about 35 million years ago (late Eocene), in a shallow marine environment (400-600 m water depth). The crater is complex and developed in a multi-layer, rheologically-variable target that comprised 400-1000 meters of soft, water-saturated sediments overlying crystalline basement. Seismic reflection data illustrates that the Chesapeake Bay crater morphology - often described as an "inverted sombrero" - is similar to other marine-target impact craters. It consists of a approx. 1 - 1.5-km deep, highly disturbed central crater, surrounded by a shallower, less deformed basin. The inner crater has a diameter of approx. 40 km; the edge of the outer basin extends to ~85-km diameter. The morphological divide between the inner and outer crater is termed the inner ring or peak ring. Little is known about the nature of the inner ring. Seismic reflection data show that the underlying basement is modestly uplifted; however, it is unclear whether the pristine surface expression of the inner ring was elevated above the floor of the outer crater.

Collins, G. S.↗