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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗
NASA NTRS1993

Numerical modeling of impact erosion of atmospheres: Preliminary results

It is clear from the great diversity of atmospheres among the terrestrial planets that their formation and evolution must have depended on a balance among a number of different processes. One of these processes is atmospheric erosion by impacts, which may have been particularly effective on Mars. The reason is that geomorphic evidence on Mars suggests that this planet had, early in its history, dense enough atmosphere to sustain active precipitation over geologically significant periods of time. Analytic calculations indicate that neither the projectile entering the atmosphere nor the main crater ejecta can cause the lose of significant amounts of atmosphere. The vapor plume that is formed, however, expands rapidly as its internal energy is converted into kinetic energy, and may blow off the overlying atmosphere. A model of this part of the impact/atmosphere interaction predicts Mars could have lost a substantial early atmosphere by impact erosion alone. Although our more detailed calculations, which took into account the anisotropy of the atmosphere with respect to zenith angle, show that the process isn't quite as effective, they still indicate the probability of substantial atmospheric loss from Mars. The first results from 2-D hydrocode runs are discussed. These include two runs which make most of the same simplifying approximations as the analytic models, in order to compare the analytic and numerical results directly, and one run (as yet incomplete) that models the full impact.

Vickery, A. M.↗
NASA NTRS1994

Biospheric effects of a large extraterrestrial impact: Case study of the Cretaceous/Tertiary boundary crater

The Chicxulub Crater in Yucatan, Mexico, is the primary candidate for the impact that caused mass extinctions at the Cretaceous/Tertiary boundary. The target rocks at Chicxulub contain 750 to 1500 m of anhydrite (CaSO4), which was vaporized upon impact, creating a large sulfuric acid aerosol cloud. In this study we apply a hydrocode model of asteroid impact to calculate the amount of sulfuric acid produced. We then apply a radiative transfer model to determine the atmospheric effects. Results include 6 to 9 month period of darkness followed by 12 to 26 years of cooling.

Pope, Kevin O.↗
NASA NTRS1995

Biospheric effects of a large extraterrestrial impact: Case study of the cretaceous/tertiary boundary crater

The Chicxulub impact crater, buried in the Yucatan carbonate platform in Mexico, is the site of the impact purported to have caused mass extinctions at the Cretaceous/Tertiary (K/T) boundary. A recently discovered Chicxulub ejecta deposit in Belize contains evidence of carbonate vaporization and precipitation from the vapor plume. Sulfate clasts are almost absent in the Belize ejecta, but are abundant in the coarse ejecta near the crater rim, hwich may reflect the greater abundance of sulfates deep in the target section. The absence of sulfate precipitates in Belize may indicate that most of the vaporized sulfur was deposited in the upper atmosphere. Hydrocode modeling of the impact indicates that between 0.4 to 7.0 x 10(exp 17) g of sulfur were vaporized by the impact in sulfates. Laser experiments indicate that SO2, SO3, and SO4 are produced, and that complex chemical reactions between plume constituents occur during condensation. The sulfur released as SO3 or SO4 converted rapidly into H2HO4 aerosol. A radiative transfer model coupled with a model of coagulation predicts that the aerosol prolonged the initial blackout period caused by impact dust only if it contained impurities. The sulfur released as SO2 converted to aerosol slowly due to the rate limiting oxidation of SO2. Radiative transfer calculations combined with rates of acid production, coagulation, and diffusion indicate that solar transmission was reduced to 10-20 percent of normal for a period of 8-13 years. This reduction produced a climate forcing (cooling) of -300 Wm(exp -2), which far exceeded the +8 Wm(exp -2) greenhouse warming caused by the CO2 released through the vaporization of carbonates, and therefore produced a decade of freezing and near-freezing temperatures. Several decades of moderate warming followed the decade of severe cooling due to the long residence time of CO2. The prolonged impact winter may have been a major cause of the K/T extinctions.

Pope, Kevin O.↗
NASA NTRS1994

Explosion of comet Shoemaker-Levy 9 on entry into the Jovian atmosphere

We use the astrophysical hydrocode ZEUS to compute high-resolution models of the disruption and deceleration of cometary fragments striking Jupiter. We find that simple analytic and semianalytic models work well for kilometer-size impactors. We show that previous numerical models that placed the explosion much deeper in the atmosphere failed to fully resolve important gasdynamical instabilities. These instabilities tear the comet apart, greatly increase its effective cross section, and bring it to an abrupt halt. A 1 km diameter fragment loses over 90% of its kinetic energy within a single scale height at an atmospheric pressure of order 10 bars. For all practical purposes, it explodes.

Mac Low, Mordecai-Mark↗
NASA NTRS1994

Impact winter and the Cretaceous/Tertiary extinctions: Results of a Chicxulub asteroid impact model

The Chicxulub impact crater in Mexico is the site of the impact purported to have caused mass extinctions at the Cretaceous/Tertiary (K/T) boundary. 2-D hydrocode modeling of the impact, coupled with studies of the impact site geology, indiate that between 0.4 and 7.0 x 10(exp 17) g of sulfur were vaporized by the impact into anhydrite target rocks. A small portion of the sulfur was released as SO3 or SO4, which converted rapidly into H2SO4 aerosol and fell as acid rain. A radiative transfer model, coupled with a model of coagulation indicates that the aerosol prolonged the initial blackout period caused by impact dust only if the aerosol contained impurities. A larger portion of sulfur was released as SO2, which converted to aerosol slowly, due to the rate-limiting oxidation of SO2. Our radiative transfer calculations, combined with rates of acid production, coagulation, and diffusion indicate that solar transmission was reduced to 10-20% of normal for a period of 8-13 yr. This reduction produced a climate forcing (cooling) of -300 W/sq.m, which far exceeded the +8 W/sq.m greenhouse warming, caused by the CO2 released through the vaporization of carbonates, and therefore produced a decade of freezing and near-freezing temperatures. Several decades of moderate warming followed the decade of severe cooling due to the long residence time of CO2. The prolonged impact winter may have been a major cause of the K/T extinctions.

Pope, Kevin O.↗
NASA NTRS1997

Comet Shoemaker-Levy 9 Fragment Size Estimates: How Big was the Parent Body?

The impact of Comet Shoemaker-Levy 9 on Jupiter in July, 1994 was the largest, most energetic impact event on a planet ever witnessed. Because it broke up during a close encounter with Jupiter in 1992, it was bright enough to be discovered more than a year prior to impact, allowing the scientific community an unprecedented opportunity to assess the effects such an event would have. Many excellent observations were made from Earth-based telescopes, the Hubble Space Telescope (HST), and the Galileo spacecraft en route to Jupiter. In this paper, these observations are used in conjunction with computational simulations performed with the CTH shock-physics hydrocode to determine the sizes of the fifteen fragments that made discernible impact features on the planet. To do this, CTH was equipped with a radiative ablation model and a postprocessing radiative ray-trace capability that enabled light-flux predictions (often called the impact flash) for the viewing geometries of Galileo and ground-based observers. The five events recorded by Galileo were calibrated to give fragment size estimates. Compared against ground-based and HST observations, these estimates were extended using a least-squares analysis to assess the impacts of the remaining ten fragments. Some of the largest impacts (L, G, and K) were greater that 1 km in diameter, but the density of the fragments was low, about 0.25 g/cm(exp 3). The volume of the combined fifteen fragments would make a sphere 1.8 km in diameter. Assuming a prebreakup density of 0.5 g/cm(exp 3), the parent body of Shoemaker-Levy 9 had a probable diameter of 1.4 km. The total kinetic energy of all the impacts was equivalent to the explosive yield of 300 Gigatons of TNT.

Crawford, David A.↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
NASA NTRS1997

A comparative modeling of supernova 1993J

The light curve of Supernova 1993J is calculated using two computational radiation transport approaches. The two approaches are represented by the computer codes STELLA and EDDINGTON. The emphasis is on the shock breakout and the photometry in the U, B and V bands during the first 120 days of the supernova. The STELLA model includes implicit hydrodynamics and is able to model early supernova evolution before the expansion is homologous. The STELLA model employs multi-group photonics and is able to follow the radiation as it decouples from the matter. The EDDINGTON code uses an algorithm for integrating the transport equation which assumes homologous expansion and uses a finer frequency resolution. The agreement between the two codes is considered to be satisfactory only in the case where compatible physical assumptions are made concerning the opacity. The assumptions are justified. The continuum spectrum for SN 1993J is predicted near the shock breakout to be superior to that predicted by standard single energy group hydrocodes. The uncertainties involved in current time dependent models of supernova light curves are discussed.

Blinnikov, Sergei↗
NASA NTRS2003

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗
NASA NTRS2002

Impact-Mobilized Dust in the Martian Atmosphere

We consider dust production and entrainment into the atmosphere of Mars by impacts. Numerical simulations based on the multidimensional multimaterial hydrocode were conducted for impactors 1 to 100 m in size and velocities 11 and 20 kilometers per second. The size distribution of particles was based on experimentrr wing TNT explosions. Dust can be mobilized even when the impactor does not reach the ground through the release of energy in the atmosphere, We found that the blast produced winds entrained dust by a mechanism similar to boundary layer winds as determined from the wind-tunnel tests. For a l-m radius stony asteroid releasing its energy in the atmosphere the lifted mass of dust is larger than that in a typical dust devil and could trigger local dust storms, For a 100-m-radius meteoroid the amount of injected dust is comparable with the tota! mass of a global dust storm.

Nemtchinov, I. V.↗
NASA NTRS2003

Numerical Modeling of Shatter Cones Development in Impact Craters

Shatter cones are the characteristic forms of rock fractures in impact structures. They have been used for decades as unequivocal fingerprints of meteoritic impacts on Earth. The abundant data about shapes, apical angles, sizes and distributions of shatter cones for many terrestrial impact structures should provide insights for the determination of impact conditions and characteristics of shock waves produced by high-velocity projectiles in geologic media. However, previously proposed models for the formation of shatter cones do not agree with observations. For example, the widely accepted Johnson-Talbot mechanism requires that the longitudinal stress drops to zero between the arrival of the elastic precursor and the main plastic wave. Unfortunately, observations do not support such a drop. A model has been also proposed to explain the striated features on the surface of shatter cones but can not invoked for their conical shape. The mechanism by which shatter cones form thus remains enigmatic to date. In this paper we present a new model for the formation of shatter cones. Our model has been tested by means of numerical simulations using the hydrocodes SALE 2D enhanced with the Grady-Kipp-Melosh fragmentation model.

Baratoux, D.↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Impacts: Modeling and Observations

This document covers the following topics: Cratering on Titan: Projectiles, Craters and Impact Melt; The Cratering Database: Making Code Jockeys Honest; Popigai Impact Structure Modeling: Morphology and Worldwide Ejecta; Anhydrite EOS and Phase Diagram in Relation to Shock Decomposition; Computational Investigations of the Chesapeake Bay Impact Structure; Hydrocode Simulations of the Chesapeake Bay Impact; Lockne Crater as a Result of Oblique Impact; The Influence of a Deep Shelf Sea on the Excavation and Modification of a Marine-Target Crater, the Lockne Crater, Central Sweden; Pre-Drilling Investigation of the Lake Bosumtwi Impact Crater: Constraints from Geophysics and Numerical Modelling; Central Uplift Formation at the Middlesboro Impact Structure, Kentucky, USA; A SRTM Investigation of Serra da Cangalho Impact Structure, Brazil; Brazilian Impact Craters: A Review; Flynn Creek Impact Structure: New Insights from Breccias, Melt Features, Shatter Cones, and Remote Sensing; The Howell Structure, Lincoln County, Tennessee: A Review of Past and Current Research; After the Chicxulub Impact: Control on Depositional and Diagenetic History of the Cenozoic Carbonate Formations of the Northwestern Yucatan Peninsula, Mexico; Ni Contents by Non-Destructive In-Situ XRF Method of Takamatsu-Kagawa Crater District in Japan; and Akiyoshi Limestone Blocks Transported by the P/T Boundary Event to Japan Islands.

Source record↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Impacts on Mars and Earth

The session "Impacts on Mars and Earth" included the following reports:A Timescale for Major Events in Early Mars Crustal Evolution; Impact Craters in Arabia Terra, Mars; Ancient Giant Basin/Aquifer System in the Arabia Region, Mars; Small Crater Morphology Within Gusev Crater and Isidis Planitia: Evidence for Widespread Secondaries on Mars; Martian Surface Properties: Inferences from Resolved Differences in Crater Geometries; Asymmetric Ramparts of Secondary Craters in Cerberus, Mars: Interpretation and Implications for Local Stratigraphy; Depth to Diameter Relationships of Craters in the High Latitudes (70 -80 ) of Mars: Implications for Geologic History of Those Areas; Impact Cratering on Mars: Search for Target Influence on Morphology; Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling; Impact-induced Hydrothermal System at the Sudbury Crater: Duration, Temperatures, Mechanics, and Biological Implications; Major and Trace Element Variations in Impact Crater Clay from Chicxulub, Lonar, and Mistastin, Implications for the Martian Soil; Stable Isotopes and Hydrothermal Fluid Source in the Yaxcopoil-1 Borehole, Chicxulub Impact Structure, Mexico; and Chicxulub Impact Melts: Geochemical Signatures of Target Lithology Mixing and Post-Impact Hydrothermal Fluid Processes.

Source record↗
NASA NTRS2005

Shock Melting of Permafrost on Mars: Water Ice Multiphase Equation of State for Numerical Modeling and Its Testing

The presence of water/ice/brine in upper layers of Martian crust affects many processes of impact cratering. Modeling of these effects promises better understanding of Martian cratering records. We present here the new ANEOS-based multiphase equation of state for water/ice constructed for usage in hydrocodes and first numerical experiments on permafrost shock melting. Preliminary results show that due to multiple shock compression of ice inclusions in rocks the entropy jump in shocked ice is smaller than in pure ice for the same shock pressure. Hence previous estimates of ice melting during impact cratering on Mars should be re-evaluated. Additional information is included in the original extended abstract.

Ivanov, B. A.↗
NASA NTRS2005

Numerical Modelling of the Deep Impact Mission Experiment

NASA s Deep Impact Mission (launched January 2005) will provide, for the first time ever, insights into the interior of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of a comets nucleus. Although it is usually assumed that comets consist of a very porous mixture of water ice and rock, little is known about the internal structure and in particular the constitutive material properties of a comet. It is therefore difficult to predict the dimensions of the excavated crater. Estimates of the crater size are based on laboratory experiments of impacts into various target compositions of different densities and porosities using appropriate scaling laws; they range between 10 s of meters up to ~250 m in diameter [1]. The size of the crater depends mainly on the physical process(es) that govern formation: Smaller sizes are expected if (1) strength, rather than gravity, limits crater growth; and, perhaps even more crucially, if (2) internal energy losses by pore-space collapse reduce the coupling efficiency (compaction craters). To investigate the effect of pore space collapse and strength of the target we conducted a suite of numerical experiments and implemented a novel approach for modeling porosity and the compaction of pores in hydrocode calculations.

Wuennemann, K.↗