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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

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.↗
NASA NTRS1991

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.↗
NASA NTRS1992

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.↗
NASA NTRS1993

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.↗
NASA NTRS1993

Collisional and dynamical history of Gaspra

Interpretation of the impact record on Gaspra requires understanding of the effects of collisions on a target body of Gaspra's size and shape, recognition of impact features that may have different morphologies from craters on larger planets, and models of the geological processes that erase and modify impact features. Crater counts on the 140 sq km of Gaspra imaged at highest resolution by the Galileo spacecraft show a steep size-frequency distribution (cumulative power-law index near -3.5) from the smallest resolvable size (150 m diameter) up through the large feature (1.5 km diameter crater) of familiar crater-like morphology. In addition, there appear to be as many as eight roughly circular concavities with diameters greater than 3 km visible on the asteroid. If we restrict our crater counts to features with traditionally recognized crater morphologies, these concavities would not be included. However, if we define craters to include any concave structures that may represent local or regional damage at an impact size, then the larger features on Gaspra are candidates for consideration. Acceptance of the multi-km features as craters has been cautious for several reasons. First, scaling laws (the physically plausible algorithms for extrapolating from experimental data) indicate that Gaspra could not have sustained such large-crater-forming impacts without being disrupted; second, aside from concavity, the larger structures have no other features (e.g. rims) that can be identified with known impact craters; and third, extrapolation of the power-law size distribution for smaller craters predicts no craters larger than 3 km over the entire surface. On the other hand, recent hydrocode modeling of impacts shows that for given impact (albeit into a sphere), the crater size is much larger than given by scaling laws. Gaspra-size bodies can sustain formation of up to 8-km craters without disruption. Besides allowing larger impact craters, this result doubles the lifetime since the last catastrophic fragmentation event up to one billion years. Events that create multi-km craters also globally damage the material structure, such that regolith is produced, whether or not Gaspra 'initially' had a regolith, contrary to other models in which initial regolith is required in order to allow current regolith. Because the globally destructive shock wave precedes basin formation, crater size is closer to the large size extrapolated from gravity-scaling rather than the strength-scaling that had earlier been assumed for such small bodies. This mechanism may also help explain the existence of Stickney on Phobos. Moreover, rejection of the large concavities as craters based on unfamiliar morphology would be premature, because (aside from Stickney) we have no other data on such large impact structures on such a small, irregular body. The eight candidate concavities cover an area greater than that counted for smaller craters, because they are most apparent where small craters cannot be seen: on low resolution images and at the limb on high resolution images. We estimate that there are at least two with diameter greater than 4 km per 140 sq km, which would have to be accounted for in any model that claims these are impact craters.

Greenberg, R.↗
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.↗