Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6

A Short Survey of Current Reactive Burn Model Capabilities as of Mid-Year 2025

The authors were tasked with writing a concise memo (now a short report) that adequately addresses the current ``state-of-the-art'' in the field of reactive flow modeling and burn models for mid-year 2025, along with identifying some of the modeling gaps. It is assumed that the reader has experience with burn models and running hydrocode simulations, and thus every effort is made to ensure brevity, clarity, and utility, so that this document will serve as a quick and useful reference. Note that “mid-year” was used in the title, as it is possible that advancements in artificial intelligence (AI) may significantly impact the field before the year is over.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predicting Asteroid Material Properties from a DART-like Kinetic Impact

Abstract NASA’s Double Asteroid Redirection Test (DART) mission is the first full-scale test of the kinetic impactor method for asteroid deflection, in which a spacecraft intentionally impacts an asteroid to change its trajectory. DART represents an important first step for planetary defense technology demonstration, providing a realistic assessment of the effectiveness of the kinetic impact approach on a near-Earth asteroid. The momentum imparted to the asteroid is transferred from the impacting spacecraft and enhanced by the momentum of material ejected from the impact site. However, the magnitude of the ejecta contribution is dependent on the material properties of the target. These properties, such as strength and shear modulus, are unknown for the DART target asteroid, Dimorphos, as well as most asteroids since such properties are difficult to characterize remotely. This study examines how hydrocode simulations can be used to estimate material properties from information available post-impact, specifically the asteroid size and shape, the velocity and properties of the impacting spacecraft, and the final velocity change imparted to the asteroid. Across >300 three-dimensional simulations varying seven material parameters describing the asteroid, we found many combinations of properties could reproduce a particular asteroid velocity. Additional observations, such as asteroid mass or crater size, are required to further constrain properties like asteroid strength or outcomes like the momentum enhancement provided by impact ejecta. Our results demonstrate the vital importance of having as much knowledge as possible prior to an impact mission, with key material parameters being the asteroid’s mass, porosity, strength, and elastic properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Dimorphos’s Material Properties and Estimates of Crater Size from the DART Impact

On 2022 September 26, the Double Asteroid Redirection Test (DART) spacecraft intentionally collided with Dimorphos, the moon of the binary asteroid system 65803 Didymos. This collision provided the first full-scale test of a kinetic impactor for planetary defense. Images from DART’s DRACO camera revealed Dimorphos to be an oblate spheroid covered in boulders of varying sizes and shapes. Very little was known about Dimorphos prior to DART’s impact, including its shape, structure, and material properties. Approach observations and those following the DART impact have provided crucial knowledge that narrows the parameter space relevant to modeling the impact into Dimorphos. Here we present the results of a suite of hydrocode simulations of the DART impact on Dimorphos. Despite remaining uncertainties, initial models of DART’s kinetic impact provide important information about the results of DART (e.g., potential crater size and morphology, ejecta mass) and the properties of Dimorphos. Simulations here suggest that Dimorphos has near-surface strength ranging from a few Pascals to tens of kPa, which corresponds to crater sizes of ~40–60 m. Simulated crater sizes provide a crucial comparison metric for the European Space Agency Hera mission when it arrives at the Didymos system. Hera’s measurement of crater size in combination with measurement of Dimorphos’s mass will allow us to assess our simulations and provide the information needed to make the DART impact experiment both the first test of a planetary defense mitigation mission and the first full-scale planetary defense simulation validation exercise.

36 MATERIALS SCIENCE↗

Nonlinear calculation of a radial mode of oscillation of a hot white dwarf

A linearly unstable first harmonic radial mode of oscillation (period 0.86 s) of a hot white dwarf with a hydrogen envelope is analyzed in the nonlinear regime by application of the Stellingwerf technique to a nonlinear hydrocode. The mass and equilibrium luminosity of the hot white dwarf are 1.3 solar masses, and 1.6 x 10 to the 37th erg/s, respectively. A stable nonlinear mode is found to exist, in which appreciable shock damping occurs at and near the surface - the shock damping being the dominant stabilizing effect. This enhances the possibility that such an object can produce a pulsating X-ray source.

Degregoria, A. J.↗

Constraints on extragalactic-jet characteristics from lobe observations

Very Large Array spectral and polarimetric observations of the lobes of luminous, double radio galaxies are discussed. These observations can provide information on the jets responsible for the sources. Spectral steepening is usually observed in the lobes, with the radio spectral index increasing with distance from the hot spot. These data can be used to infer a 'speed of separation' of the hot spot and lobe material. These speeds are in agreement with hydrocode models of jets. Polarimetric observations indicate that in at least some sources, there is an undectably small amount of internal Faraday rotation, indicating upper limits to the thermal-plasma density of a few times 10 to the -5th/cu cm or less. These measurements are also in agreement with the numerical beam models, provided that the beam density is substantially less than that of the background medium. It is concluded that the lobe observations indicate that these sources are powered by light, high-Mach-number beams.

Spangler, S. R.↗

Collisional stripping of Mercury's mantle

A three-dimensional smoothed-particle hydrocode is the basis of the present numerical simulations of conditions under which a giant collision between a proto-Mercury and a planet one-sixth its size would lead to the loss of most of the silicate mantle of Mercury and thereby account for its anomalously high density. A head-on collision at 20 km/sec, and an off-axis impact parameter of half the radius of the proto-Mercury at 35 km/sec, are approximately equal in damage yielded; both will yield a remnant whose characteristics are those of the present Mercury.

Benz, Willy↗

Application and implementation of transient algorithms in computer programs

A brief introduction is given to the nonlinear finite element programs developed at Lawrence Livermore National Laboratory. The four programs are DYNA3D and DYNA2D, which are explicit hydrocodes, and NIKE3D and NIKE2D, which are implicit programs. The main emphasis is on DYNA3D with asides about the other programs. During the past year several new features were added to DYNA3D, and major improvements were made in the computational efficiency of the shell and beam elements. Most of these new features and improvements will eventually make their way into the other programs. The emphasis in the computational mechanics effort was always, and continues to be, efficiency. To get the most out of the supercomputers, all Crays, the programs were vectorized where possible. Several of the more interesting capabilities of DYNA3D will be described and the impact on efficiency will be discussed. Some of the recent work on NIKE3D and NIKE2D will also be presented. In the belief that a single example is worth a thousand equations, the theory is skipped entirely and the examples presented.

Benson, David J.↗

Shielding requirements for the Space Station habitability modules

The design, analysis, development, and tests of the total meteoroid/debris protection system for the Space Station Freedom habitability modules, such as the habitation module, the laboratory module, and the node structures, are described. Design requirements are discussed along with development efforts, including a combination of hypervelocity testing and analyses. Computer hydrocode analysis of hypervelocity impact phenomena associated with Space Station habitability structures is covered and the use of optimization techniques, engineering models, and parametric analyses is assessed. Explosive rail gun development efforts and protective capability and damage tolerance of multilayer insulation due to meteoroid/debris impact are considered. It is concluded that anticipated changes in the debris environment definition and requirements will require rescoping the tests and analysis required to develop a protection system.

Avans, Sherman L.↗

Impact processes in the Solar System: New understandings through numerical modeling

A collision of two rocky objects circling the sun in space, each roughly the size and mass of a large mountain range, was modeled. A fragmentation hydrocode was developed to perform dynamical computations of collisional outcomes. Explosive framentation and fluid dynamics were used and drawn together into a single application. To model a solid, certain material parameters, such as density, elasticity, rigidity, and energies of melting and vaporization were input. These parameters are well-known for a variety of important materials, such as ice, iron, granite, and basalt. Another important parameter used is the distribution of initial flaws within the material.

Asphaug, E.↗

Freedom Station wall design using hydrodynamic modelling

The paper outlines selected outputs from an ongoing Marshall Space Flight Center/US Army Corps of Engineers parametric study of the effects of meteoroid/space debris impacts on Space Station Freedom module wall protection structures. The advantages and limitations of HULL hydrocode computer simulation are discussed, and examples of meteoroid/ debris impacts on various wall structural designs are presented. Trends in the terminal effects of particle sizes, velocities, shapes, and impact angles upon structural design parameters (wall and bumper thickness and spacing) are summarized and depicted through selected simulation runs. Throughout the paper, proposed structural modifications to the baseline module protection structure will be examined which are capable of mitigating damage from particle impacts.

Williamsen, Joel E.↗

Energy coupling in catastrophic collisions

The prediction of events leading to the catastrophic collisions and disruption of solar system bodies is fraught with the same difficulties as are other theories of impact events; since one simply cannot perform experiments in the regime of interest. In the catastrophic collisions of asteroids that regime involves bodies of a few tons to hundred of kilometers in diameter, and velocities of several kilometers pre second. For hundred kilometer bodies, gravitational stresses dominate material fracture strengths, but those gravitational stresses are essentially absent for laboratory experiments. Only numerical simulations using hydrocodes can in principle analyze the true problems, but they have their own major uncertainties about the correctness of the physical models and properties. The question of the measure of the impactor and its energy coupling is investigated using numerical code calculations. The material model was that of a generic silicate rock, including high pressure melt and vapor phases, and includes material nonlinearity and dissipation via a Mie-Gruniesen model. A series of calculations with various size ratios and impact velocities are reported.

Holsapple, K. A.↗

Continuum modeling of catastrophic collisions

A two dimensional hydrocode based on 2-D SALE was modified to include strength effects and fragmentation equations for fracture resulting from tensile stress in one dimension. Output from this code includes a complete fragmentation summary for each cell of the modeled object: fragment size (mass) distribution, vector velocities of particles, peak values of pressure and tensile stress, and peak strain rates associated with fragmentation. Contour plots showing pressure and temperature at given times within the object are also produced. By invoking axial symmetry, three dimensional events can be modeled such as zero impact parameter collisions between asteroids. The code was tested against the one dimensional model and the analytical solution for a linearly increasing tensile stress under constant strain rate.

Ryan, Eileen V.↗

Analysis of energy dissipation and deposition in elastic bodies impacting at hypervelocities

A series of impact problems were analyzed using the Eulerian hydrocode CTH. The objective was to quantify the amount of energy dissipated locally by a projectile-infinite plate impact. A series of six impact problems were formulated such that the mass and speed of each projectile were varied in order to allow for increasing speed with constant kinetic energy. The properties and dimensions of the plate were the same for each projectile impact. The resulting response of the plate was analyzed for global Kinetic Energy, global momentum, and local maximum shear stress. The percentage of energy dissipated by the various hypervelocity impact phenomena appears as a relative change of shear stress at a point away from the impact in the plate.

Medina, David F.↗

Preliminary design of a meteoroid/orbital debris shield system for a Mars mission spacecraft

A preliminary design of a spacecraft shield system to defeat meteoroids and orbital debris during a Mars mission is presented. The mission scenario is first defined in terms of stage times which include LEO, transit, low-Mars orbit, and Mars surface periods. The environment definitions for earth-orbital space debris, planetary meteoroids, and interplanetary meteoroids are next introduced. Shield system design approaches incorporate stochastic simulation, hydrocode analyses, hypervelocity impact testing, and optimization techniques. Structural design trades presented include spacecraft configuration, mission schedule, penetration risk, and total standoff distance between the bumper and wall.

Robinson, J. H.↗

Momentum transfer in the one-dimensional impact of spaced plates

The two-dimensional hydrocode CSQ III was used to calculate the fraction of momentum transferred from a flyer plate to a target of two spaced plates. The effect of the vaporization phase transition, as calculated with the ANEOS analytical complete three-phase equation of state, was estimated. Application of these results to the protection of spacecraft from meteoroids and orbital debris is discussed.

Abbott, R. D.↗

Hypervelocity impact simulations of Whipple shields

The problem associated with protecting space vehicles from space debris impact is described. Numerical simulation is espoused as a useful complement to experimentation: as a means to help understand and describe the hypervelocity impact phenomena. The capabilities of a PC-based hydrocode, ZeuS, are described, for application to the problem of hypervelocity impact. Finally, results of ZeuS simulations, as applied to the problem of bumper shield impact, are presented and compared with experimental results.

Segletes, Steven B.↗

Ejecting basaltic achondrites from Vesta: Hydrodynamical impact models

Vesta is a large asteroid whose crust is mostly basaltic. Spectral heterogeneity suggests a sizable olivine feature which may be explained as impact excavation. The spectral data probably show a localized feature approximately 200 km in diameter or a diffuse feature approximately 400 km in diameter. Lightcurve irregularities suggest heterogeneity on a similar scale. This heterogeneity may represent the crater bowl, the extent of its ejecta deposit, or indeed something unrelated to cratering. In any case, drawing direct inferences about the state of Vesta's surface on the basis of these observations involves substantial speculation. We are still far from understanding the boundary between cratering and catastrophic disruption, particularly on targets for which strength and self-gravity both matter. But we are now able to model the underlying physical process -- dynamic fragmentation -- accurately with fragmentation hydrocodes such as SALE 2D and SPH3D. We present several impact scenarios for Vesta; our study is similar to a previous impact model for the formation of Stickney crater on Phobos. We illustrate the effects of impactors of different sizes and velocities, and the effect of gravity and structural inhomogeneity.

Asphaug, E.↗