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 253 records · Page 14

Motion of the heliospheric termination shock. 4: MHD effects

In this paper we generalize earlier gasdynamic analyses of the motion of the heliospheric termination shock in response to upstream disturbances (Barnes, 1993, 1994; Naidu and Barnes, 1994), to include magnetohydrodynamic (MHD) phenomena. We assume that the termination shock is a strong, perpendicular shock and that the initial upstream disturbance is a tangential discontinuity. The resulting configuration after the interaction is very similar to that in the gasdynamic models after an interaction with a contact discontinuity or interplanetary shock, and for an increase (decrease) in dynamic pressure consists of an outward (inward) propagating termination shock and an outward propagating shock (MHD rarefraction wave) that carries the signal of the disturbance into the far downstream plasma. The plasma immediately behind the new termination shock is separated from the downstream signal by a tangential discontinuity. The results of the model show that the speed of the new termination shock depends mainly on the magnitude of the change in dynamic pressure and are typically of order approximately 100 km/s, comparable to the results of the gasdynamic models.

Naidu, Kamcilla↗

Motion of the Heliospheric Termination Shock: MHD Effects - 4

In this paper we generalize earlier gasdynamic analyses of the motion of the heliospheric termination shock in response to upstream disturbances to include magnetohydrodynamic (MHD) phenomena. We assume that the termination shock is a strong, perpendicular shock and that the initial upstream disturbance is a tangential discontinuity. The resulting configuration after the interaction is very similar to that in the gasdynamic models after an interaction with a contact discontinuity or interplanetary shock, and for an increase (decrease) in dynamic pressure consists of an outward (inward) propagating termination shock and an outward propagating shock (MHD rarefaction wave) that carries the signal of the disturbance into the far downstream plasma. The plasma immediately behind the new termination shock is separated from the downstream signal by a tangential discontinuity. The results of the model show that the speed of the new termination shock depends mainly on the magnitude of the change in dynamic pressure and are typically of order -100 km/s, comparable to the results of the gasdynamic models.

Naidu, Kamcilla↗

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pyrotechnic shock at the orbiter/external tank forward attachment

During the initial certification test of the forward structural attachment of the space shuttle orbiter to the external tank, pyrotechnic shock from actuation of the separation device resulted in structural failure of the thermal protection tiles surrounding the attachment. Because of the high shock associated with the separation bolt, the development of alternative low shock separation designs was initiated. Two concepts that incorporate a 5.08 centimeter frangible nut as the release device were developed and tested.

Rogers, W. F.↗

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗

Shock compression of a low-density carbon foam

Porous polymer foams are found in numerous applications where they may be subjected to dynamic compression. The shockwave properties of polymer foams have several features that make shock Hugoniot measurements and interpretation challenging. Here, we present a comparison of dynamic X-ray phase contrast imaging and traditional optical velocimetry approaches to measuring the shocked states for a low-density carbon foam. Overall, dynamic X-ray phase contrast imaging was found to reduce errors associated with the measured Hugoniot states. In this study, the foam exhibits anomalous compression as expected due to its high initial porosity arising from shock heating.

36 MATERIALS SCIENCE↗

Ion reflection and downstream thermalization at the quasi-parallel bow shock

Using the results of ISEE 2 plasma and magnetic field measurements, two features related to the ion thermalization process at high-Mach-number (M above 2) quasi-parallel collisionless shocks are discussed. These are the presence of a coherent secondary beam of ions within the shock layer which is considered to be produced by reflection, and downstream ion distributions which contain both a relatively cold core of directly transmitted ions and a hotter 'shell' of ions, which appear to result from the disruption and scattering of ions initially reflected at the shock. Evidence is presented that coherent ion reflection is an important element of the ion energy dissipation process at high-Mach-number quasi-parallel shocks.

Gosling, J. T.↗

Passage of Bolides Through the Atmosphere

Different fragmentation models are applied to a number of events, including the entry of TC3 2008 asteroid in order to reproduce existing observational data. Keywords meteoroid entry fragmentation modeling 1 Introduction Fragmentation is a very important phenomenon which occurs during the meteoroid entry into the atmosphere and adds more drastic effects than mere deceleration and ablation. Modeling of bolide fragmentation (100 106 kg in mass) may be divided into several approaches. Detail fitting of observational data (deceleration and/or light curves) allows the determination of some meteoroid parameters (ablation and shape-density coefficients, fragmentation points, amount of mass loss) (Ceplecha et al. 1993; Ceplecha and ReVelle 2005). Observational data with high accuracy are needed for the gross-fragmentation model (Ceplecha et al. 1993), which is used for the analysis of European and Desert bolide networks data. Hydrodynamical models, which describe the entry of the meteoroid including evolution of its material, are applied mainly for large bodies (>106 kg) (Boslough et al. 1994; Svetsov et al. 1995; Shuvalov and Artemieva 2002, and others). Numerous papers were devoted to the application of standard equations for large meteoroid entry in the attempts to reproduce dynamics and/or radiation for different bolides and to predict meteorite falls. These modeling efforts are often supplemented by different fragmentation models (Baldwin and Sheaffer, 1971; Borovi.ka et al. 1998; Artemieva and Shuvalov, 2001; Bland and Artemieva, 2006, and others). The fragmentation may occur in different ways. For example, few large fragments are formed. These pieces initially interact through their shock waves and then continue their flight independently. The progressive fragmentation model suggests that meteoroids are disrupted into fragments, which continue their flight as independent bodies and may be disrupted further. Similar models were suggested in numerous papers, beginning with Levin (1956) and initial interaction of fragments started to be taken into account after the paper by Passey and Melosh (1980). The progressive fragmentation model with lateral spreading of formed fragments is widely used (Artemieva and Shuvalov, 1996; Nemtchinov and Popova, 1997; Borovi.ka et al. 1998; Bland and Artemieva, 2006).

Popova, O.↗

Anomalous transport by Kelvin-Helmholtz instabilities

A simulation of magnetohydrodynamic Kelvin-Helmholtz instabilities has been performed for parallel and transverse configurations, modeling high latitude (or downstream flanks) and dayside low latitude magnetospheric boundaries. In the parallel configuration, a super-Alfvenic and transsonic shear flow develops into small eddies, which strongly compresses, twists, and hence amplifies the magnetic field by the dynamo action with an amplification factor of M(A)/2. In the nonlinear stage, however large the initial Alfven mach number M(A) may be, the magnetic field amplified and twisted by the hydromagnetic flow vortices reacts back upon the flow evolution, and the flow vortices cascade into smaller structures. In the transverse configuration, the instability leads to the formation of a fast shock discontinuity from an initially subfast shear flow. Anomalous tangential stress by the instability in the transverse configuration reaches 1 percent of the magnetosheath momentum flux, but for the parallel configuration, the anomalous transport is 2-3 times larger than the anomalous transport in the transverse configuration. The anomalous transport for both configurations satisfies the requirement of the viscous-like interaction at the magnetospheric boundary.

Miura, A.↗

Compensatory reserve detects subclinical shock with more expeditious prediction for need of life‐saving interventions compared to systolic blood pressure and blood lactate

Abstract Introduction We conducted a prospective observational study on 205 trauma patients at a level I trauma facility to test the hypothesis that a compensatory reserve measurement (CRM) would identify higher risk for progression to shock and/or need a life‐saving interventions (LSIs) earlier than systolic blood pressure (SBP) and blood lactate (LAC). Methods A composite outcome metric included blood transfusion, procedural LSI, and mortality. Discrete measures assessed as abnormal (ab) were SBP <90 mmHg, CRM <60%, and LAC >2.0. A graded categorization of shock was defined as: no shock (normal [n] SBP [n‐SBP], n‐CRM, n‐LAC); sub‐clinical shock (ab‐CRM, n‐SBP, n‐LAC); occult shock (n‐SBP, ab‐CRM, ab‐LAC); or overt shock (ab‐SBP, ab‐CRM, ab‐LAC). Results Three patients displayed overt shock, 53 displayed sub‐clinical shock, and 149 displayed no shock. After incorporating lactate into the analysis, 86 patients demonstrated no shock, 25 were classified as sub‐clinical shock, 91 were classified as occult shock, and 3 were characterized as overt shock. Each shock subcategory revealed a graded increase requiring LSI and transfusion. Initial CRM was associated with progression to shock (odds ratio = 0.97; p < .001) at an earlier time than SBP or LAC. Conclusions Initial CRM uncovers a clinically relevant subset of patients who are not detected by SBP and LAC. Our results suggest CRM could be used to more expeditiously identify injured patients likely to deteriorate to shock, with requirements for blood transfusion or procedural LSI.

Convertino, Victor A.↗

Injection and distribution of liquids in the flow fields of blunt shapes at hypersonic speeds

A possible method for reducing the attenuation of radio transmission from hypersonic vehicles is by injection of a liquid into the flow field. The hot plasma layer surrounding the vehicle could presumably be quenched if the latent heat of the liquid and the heat capacity of the vapor could be utilized efficiently. In the case of stagnation-point injection, the momentum of the liquid could be utilized to modify the bow shock and possibly reduce the initial ionization. To attain these objectives of efficient coolant utilization or modification of the bow shock a knowledge of the liquid injection, atomization, spray penetration, and evaporation processes is required. This report presents the results of a study of liquid injection into hypersonic flow fields. A simplified theory for evaporation distances is derived and is shown to be in qualitative agreement with experimental data. Correlation of experimental data for maximum cross-current penetration of the liquid spray for injection near the shoulder pf a hemisphere cylinder, a sphere 9° cone, and a Mercury model are given. These data were obtained in the Langley Mach 8 variable density tunnel with both liquid nitrogen and water injection to provide a large range of liquid parameters. Flow-field disturbances, as indicated by schlieren photographs and surface pressures, increased with the ratio of coolant mass flux to airstream mass flux. For stagnation-point injection, the maximum shock standoff distance was correlated by the ratio of the liquid jet momentum flux to the airstream momentum flux.

Ivan E. Beckwith↗

Study of Unsteady, Sphere-Driven, Shock-Induced Combustion for Application to Hypervelocity Airbreathing Propulsion

A premixed, shock-induced combustion engine has been proposed in the past as a viable option for operating in the Mach 10 to 15 range in a single stage to orbit vehicle. In this approach, a shock is used to initiate combustion in a premixed fuel/air mixture. Apparent advantages over a conventional scramjet engine include a shorter combustor that, in turn, results in reduced weight and heating loads. There are a number of technical challenges that must be understood and resolved for a practical system: premixing of fuel and air upstream of the combustor without premature combustion, understanding and control of instabilities of the shock-induced combustion front, ability to produce sufficient thrust, and the ability to operate over a range of Mach numbers. This study evaluated the stability of the shock-induced combustion front in a model problem of a sphere traveling in a fuel/air mixture at high Mach numbers. A new, rapid analysis method was developed and applied to study such flows. In this method the axisymmetric, body-centric Navier-Stokes equations were expanded about the stagnation streamline of a sphere using the local similarity hypothesis in order to reduce the axisymmetric equations to a quasi-1D set of equations. These reduced sets of equations were solved in the stagnation region for a number of flow conditions in a premixed, hydrogen/air mixture. Predictions from the quasi-1D analysis showed very similar stable or unstable behavior of the shock-induced combustion front as compared to experimental studies and higher-fidelity computational results. This rapid analysis tool could be used in parametric studies to investigate effects of fuel rich/lean mixtures, non-uniformity in mixing, contaminants in the mixture, and different chemistry models.

Axdahl, Erik↗

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The two-dimensional structure of diffuse ions associated with the earth's bow shock

The lateral distribution of upstream, diffuse energetic ions resulting from shock acceleration at an idealized model bow shock is investigated, where the model consists of a uniform planar shock along which steady source ion injection takes place over a finite region of the shock. The solution exhibits growth toward a peak ion concentration whose location beneath the point of initial magnetic connection with the shock is sensitive to the model parameters. It is found that a plateau forms if the source size is large in comparison to the scale length of ion distribution function's growth, and that cross-field diffusion strongly enhances the diffuse ion concentration near the leading edge of the source.

Skadron, G.↗