Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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

Shock-induced volatile loss from a carbonaceous chondrite Implications for planetary accretion

Solid recovery impact-induced volatile loss experiments on the Murchison C2M meteorite indicate that for an impact of a given velocity, H2O and total volatiles are driven from the sample in the same proportion as present initially. The primitive surface volatile budget of a planet growing by accretion would have the same bulk elemental composition as the volatiles in the incident planetesimals. Incipient devolatilization of Murchison occurs at an initial shock pressure of about 11 GPa and complete devolatilization occurs at a pressure of about 30 GPa. For the earth, incipient and complete devolatilization of accreting planetesimals would occur when the planet reached approximately 12 percent and 27 percent, respectively, of its present-day radius. Impact-induced devolatilization would profoundly affect the volatile distribution within the accreting planet. Prior to metallic core formation and internal differentiation the growing planet would have a very small core with the same volatile content as the incident material, a volatile depleted mantle, and an extremely volatile rich surface. In the case of the earth, 99.4 wt pct of the total incident volatile material would end up on or near the planetary surface.

Tyburczy, James A.↗

The effect of neutrino transport on the collapse of iron stellar cores

A multigroup flux-limited diffusion approximation to neutrino transport, correct to first order in material velocities, is described. The role of neutrino-electron scattering in determining core deleptonization is studied using a modified Fokker-Planck approximation in which adjacent neutrino energy groups are coupled. The moderate deleptonizaton obtained during infall restricts homologous core masses to less than 0.8 solar mass. With 0.6-0.8 solar mass of outer core to traverse, nuclear dissociation alone stalls the shock waves produced at core bounce. In addition, neutrino energy losses remove at least another 4 x 10 to the 51st ergs. In no case is matter ejected promptly as a direct result of the initial shock wave.

Myra, E. S.↗

Cyclic behavior at quasi-parallel collisionless shocks

Large scale one-dimensional hybrid simulations with resistive electrons have been carried out of a quasi-parallel high-Mach-number collisionless shock. The shock initially appears stable, but then exhibits cyclic behavior. For the magnetic field, the cycle consists of a period when the transition from upstream to downstream is steep and well defined, followed by a period when the shock transition is extended and perturbed. This cyclic shock solution results from upstream perturbations caused by backstreaming gyrating ions convecting into the shock. The cyclic reformation of a sharp shock transition can allow ions, at one time upstream because of reflection or leakage, to contribute to the shock thermalization.

Burgess, D.↗

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗

Enhancements to viscous-shock-layer technique

A solution procedure is presented which considerably improves the computational efficiency of the viscous-shock-layer technique, especially for long slender bodies. The 'predictor-corrector' procedure suggested for obtaining the shock shape beyond the nose region requires only a single global pass. The accuracy of the present method is demonstrated by comparison with globally iterated results over the entire body and with ground- and flight-test data. A good comparison of the results computed with the two methods is shown for different flowfield chemistry models and axisymmetric body shapes. The new procedure results in computer run times 1/3 to 1/2 of the times required for the full-body global iteration procedure. Further, the algebraic expressions used to specify the initial shock shape eliminate the need for a shock shape generated by external means and permit immediate introduction of the full viscous-shock-layer equations. Finally, the present method of solution for the VSL equations provides to the aerothermal designer a very efficient and accurate tool for detailed flowfield as well as future technology studies.

Gupta, Roop N.↗

Solution of a large hydrodynamic problem using the STAR-100 computer

A representative hydrodynamics problem, the shock initiated flow over a flat plate, was used for exploring data organizations and program structures needed to exploit the STAR-100 vector processing computer. A brief description of the problem is followed by a discussion of how each portion of the computational process was vectorized. Finally, timings of different portions of the program are compared with equivalent operations on serial machines. The speed up of the STAR-100 over the CDC 6600 program is shown to increase as the problem size increases. All computations were carried out on a CDC 6600 and a CDC STAR 100, with code written in FORTRAN for the 6600 and in STAR FORTRAN for the STAR 100.

Weilmuenster, K. J.↗

Wind laws for shockless initialization

A system of diagnostic equations for the velocity field, or wind laws, was derived for each of a number of models of large-scale atmospheric flow. The derivation in each case is mathematically exact and does not involve any physical assumptions not already present in the prognostic equations, such as nondivergence or vanishing of derivatives of the divergence. Therefore, initial states computed by solving these diagnostic equations should be compatible with the type of motion described by the prognostic equations of the model and should not generate initialization shocks when inserted into the model. Numerical solutions of the diagnostic system corresponding to a barotropic model are exhibited. Some problems concerning the possibility of implementing such a system in operational numerical weather prediction are discussed.

Ghil, M.↗

Atmospheric model development in support of SEASAT. Volume 5: Test results with fine-mesh (187 by 187 by 10) model

A model (187 x 187 grid with ten layers) was used to produce a 24 hour forecast using initial conditions for 1200Z, 20 May 1976. This forecast was compared to its five layer counterpart on the 187 x 187 grid and to its 63 x 63 coarse mesh counterpart having the same number of layers. Increases in horizontal resolution lead to significant differences in a one day forecast. Many of these differences represent improvements. Increases in vertical resolution tend to produce smaller impacts on the forecast, except in the region near and above the tropopause. The effect on forecast precipitation is in the 10-20% range, as opposed to the 100% for increases in the horizontal resolution. With respect to model energetics (kinetic energy; square vorticity; square divergence), the model tend to group according to horizontal resolution. The time variations of these parameters show that dynamic initialization is needed to: (1) minimize initialization shock; and (2) stabilize the model context to prevent wash-out of small scale information during the adjustment period (first 6-12 forecast hours).

Kesel, P. G.↗

Shock tube and shock wave research; Proceedings of the Eleventh International Symposium, University of Washington, Seattle, Wash., July 11-14, 1977

Papers are presented on the applications of shock-wave technology to the study of hydrodynamics, the use of the pressure-wave machine for charging diesel engines, and measurements of the heat-transfer rate in gas-turbine components. Consideration is given to shock propagation along 90-degree bends, the explosive dissemination of liquids, and rotational and vibrational relaxation behind weak shock waves in water vapor. Shock phenomena associated with expansion flows are described and stratospheric-related research using the shock tube is outlined. Attention is given to shock-wave ignition of magnesium powders, Mach reflection and boundary layers, and transition in the shock-induced unsteady boundary layer on a flat plate. Shock-tube measurements of induction and post-induction rates for low-Btu gas mixtures are presented and shock-initiated ignition in COS-N2O-Ar mixtures is described. Cluster growth rates in supersaturated lead vapor are presented and a study of laser-induced plasma motion in a solenoidal magnetic field is reviewed.

Ahlborn, B.↗

Systematic optimization of a detailed kinetic model using a methane ignition example

An approach to the systematic optimization of a large-scale dynamic model is proposed which consists in parameterization of simulation results as response surfaces. The optimization procedure is carried out using a second-order orthogonal design. The approach proposed here is demonstrated by an example involving the shock-initiated ignition of methane.

Frenklach, M.↗

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗

Prompt arrival of solar energetic particles from far eastern events - The role of large-scale interplanetary magnetic field structure

Intensity-time profiles of solar energetic particle enhancements generally show an asymmetry with respect to the heliolongitude of the associated solar event. Particles arrive promptly from events to the west of an observer because of good magnetic connection, whereas particle enhancements from poorly connected eastern source regions usually show much slower onsets. However, some 15 percent of eastern events do show prompt onsets. Two prompt particle enhancements associated with eastern solar events are studied using data from the Goddard Space Flight Center instruments on the ISEE 3 and IMP 8 spacecraft. In both events the prompt particle onset was observed when the spacecraft were in a postshock plasma region, apparently within a magnetic bottle. It is suggested that the magnetic bottle extended back to the sun and served as a channel for fast particle propagation to the spacecraft. Particles accelerated at an expanding coronal shock initiated by the eastern event could be injected onto field lines in the foot of the bottle.

Richardson, I. G.↗

Melt droplet formation in energetic impacts

Impacts between rocky bodies at velocities exceeding about 15 km/sec are capable of melting or vaporizing both the impacting object and a portion of the target. Geological materials initially shocked to high pressure approach the liquid-vapor phase boundary from the liquid side as they decompress, breaking up into an expanding spray of liquid droplets. A simple theory is presented for estimating the sizes of these droplets as a function of impactor size and velocity. It is shown that these sizes are consistent with observations of microtektites and spherules found in the Cretaceous-Tertiary boundary layer, the Acraman impact structure, Archean beds in South Africa and lunar regolith. The model may also apply to the formation of chondrules.

Vickery, A. M.↗

A loudness calculation procedure applied to shaped sonic booms

Described here is a procedure that can be used to calculate the loudness of sonic booms. The procedure is applied to a wide range of sonic booms, both classical N-waves and a variety of other shapes of booms. The loudness of N-waves is controlled by overpressure and the associated rise time. The loudness of shaped booms is highly dependent on the characteristics of the initial shock. A comparison of the calculated loudness values indicates that shaped booms may have significantly reduced loudness relative to N-waves having the same peak overpressure. This result implies that a supersonic transport designed to yield minimized sonic booms may be substantially more acceptable than an unconstrained design.

Shepherd, Kevin P.↗

LSENS, a general chemical kinetics and sensitivity analysis code for gas-phase reactions: User's guide

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS, are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include static system, steady, one-dimensional, inviscid flow, shock initiated reaction, and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method, which works efficiently for the extremes of very fast and very slow reaction, is used for solving the 'stiff' differential equation systems that arise in chemical kinetics. For static reactions, sensitivity coefficients of all dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters can be computed. This paper presents descriptions of the code and its usage, and includes several illustrative example problems.

Radhakrishnan, Krishnan↗

Global MHD Magnetospheric Simulation of January 10, 1997 Encounter with Magnetic Cloud

We present a global MHD magnetospheric simulation of the encounter of the magnetic cloud with the terrestrial magnetosphere on January 10-11, 1997. The MHD simulation is driven by IMF and solar wind plasma measurements provided by Wind situated about 100 $R_E$ upstream from Earth. Field-aligned currents produced in the model are mapped down to the ionosphere and are directly compared to auroral images from the Polar UVI. Particular attention will be paid to the interval following the the initial shock wave arrival (around 0100 UT on the 1Oth) but preceding the passage of the magnetic cloud proper (commencing at about 0430 UT) for which there was continuous viewing of the entire auroral oval from apogee by Polar. This turbulent period is characterized by numerous dynamic pressure rises and dips and several northward and southward turnings of the IMF, all of which generate dynamic activity in the simulation that is reflected in the mapped field-aligned current patterns. As discussed by Brittnacher et al. in this session, the auroral morphology imaged by UVI during this period includes: shock wave-induced brightening of the oval followed by a pseudo-onset near midnight, several sun-aligned and curled arcs within the polar cap, an expanding polar cap cleared of arcs, and finally substorm onset at 0337 UT. These features will be directly compared to region 1 and 2 current systems as well as cusp currents in the simulation. Spacecraft magnetopause crossings are also predicted by the global simulation and will be compared to observed crossings, including numerous dayside crossings by Geotail on the 1Oth, and geosynchronous crossings early on January 11.

Elsen, R.↗

Realistic Modeling of Interaction of Quiet-Sun Magnetic Fields with the Chromosphere

High­‐resolution observations and 3D MHD simulations reveal intense interaction between the convection zone dynamics and the solar atmosphere on subarcsecond scales. To investigate processes of the dynamical coupling and energy exchange between the subsurface layers and the chromosphere we perform 3D radiative MHD modeling for a computational domain that includes the upper convection zone and the chromosphere, and investigate the structure and dynamics for different intensity of the photospheric magnetic flux. For comparison with observations, the simulation models have been used to calculate synthetic Stokes profiles of various spectral lines. The results show intense energy exchange through small­‐scale magnetized vortex tubes rooted below the photosphere, which provide extra heating of the chromosphere, initiate shock waves, and small­‐scale eruptions.

quiet-sun↗