Penetration of an initially radial shock wave through an aluminum-glass interface
Propagation of spherical shock wave through aluminum-glass interface
INTERFACE↗
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.
Initial operation of the UTA shock tunnel
Initial experience in operating a small UTA (University of Texas at Arlington) shock tunnel at high Reynolds numbers for perfect gas simulation of hypersonic flow and boundary layer studies is described. Particular attention is given to tradeoffs between various constraints for achieving the test requirements, and methods for obtaining precise control of test conditions and protecting low-range pressure transducers from overpressure damage.
Lu, Frank K.↗
Magma ocean formation due to giant impacts
The thermal effects of giant impacts are studied by estimating the melt volume generated by the initial shock wave and corresponding magma ocean depths. Additionally, the effects of the planet's initial temperature on the generated melt volume are examined. The shock pressure required to completely melt the material is determined using the Hugoniot curve plotted in pressure-entropy space. Once the melting pressure is known, an impact melting model is used to estimate the radial distance melting occurred from the impact site. The melt region's geometry then determines the associated melt volume. The model is also used to estimate the partial melt volume. Magma ocean depths resulting from both excavated and retained melt are calculated, and the melt fraction not excavated during the formation of the crater is estimated. The fraction of a planet melted by the initial shock wave is also estimated using the model.
Tonks, W. B.↗
A finite-difference solution for unsteady wave interactions with an application to the expansion tube
An analytical study of the unsteady wave interaction of a shock wave overtaking a previously initiated shock tube flow has been conducted for a real gas by solving a finite-difference representation of the one-dimensional Navier-Stokes equations. This analysis has been used to determine the effect on the available test time of opening the secondary diaphragm in the expansion-tube operating cycle prior to the arrival of the incident shock wave. Special emphasis is placed on determining expansion-tube test-time limitations resulting from a nonuniform entropy distribution occurring during the early test-gas flow due to the varying Mach number of the incident shock wave as it traverses the expansion fan caused by the preopening process. Time distributions of velocity, temperature, pressure, and entropy at various distances from the diaphragm are presented as a function of the time required to open the diaphragm.
Weilmuenster, K. J.↗
Role of Sea‐Surface Salinity in Simulating Historical Decadal Variations of Atlantic Meridional Overturning Circulation in a Coupled Climate Model
Abstract Extending climate prediction from seasonal to decadal timescales requires realistic initialization of not only upper ocean heat content but also the Atlantic meridional overturning circulation (AMOC). However, it remains a major challenge to realistically initialize AMOC in a coupled system while also maintaining a balanced atmosphere‐ocean initial state. This study demonstrates the feasibility of generating fully coupled historical states with realistic AMOC variability. Employing a forced ocean—sea‐ice (FOSI) model simulation as the “truth,” we show reproducibility of key features of historical AMOC decadal variability in a fully coupled model by restoring sea‐surface salinity, in addition to sea‐surface temperature restoring widely used in seasonal prediction. The atmospheric state of the restored coupled model solution is much closer to that of the free coupled simulation than to the observations used in FOSI, pointing to potential advantages of using this approach for initializing decadal predictions with reduced ocean initialization shock.
54 ENVIRONMENTAL SCIENCES↗
THEORY OF INITIATION OF DETONATION IN SOLID AND LIQUID EXPLOSIVES
Differences in growth of shock initiated detonation in liquid and solid explosives - energy release rate model for solid granular explosives
ENERGY CONVERSION↗
Numerical investigation of separated nozzle flows
A numerical study of axisymmetric overexpanded nozzle is presented. The flow structure of the startup and throttle-down processes are examined. During the impulsive startup process, observed flow features include the Mach disk, separation shock, Mach stem, vortex core, contact surface, slip stream, initial shock front, and shocklet. Also the movement of the Mach disk is not monotonical in the downstream direction. For a range of pressure ratios, hysteresis phenomenon occurs; different solutions were obtained depending on different processes. Three types of flow structures were observed. The location of separation point and the lower end turning point of hysteresis are closely predicted. A high peak of pressure is associated with the nozzle flow reattachment. The reversed vortical structure and affects engine performance.
Chen, C. L.↗
Shock, Post-Shock Annealing, and Post-Annealing Shock in Ureilites
The thermal and shock histories of ureilites can be divided into four periods: 1) formation, 2) initial shock, 3) post-shock annealing, and 4) post-annealing shock. Period 1 occurred approx.4.55 Ga ago when ureilites formed by melting chondritic material. Impact events during period 2 caused silicate darkening, undulose to mosaic extinction in olivines, and the formation of diamond, lonsdaleite, and chaoite from indigenous carbonaceous material. Alkali-rich fine-grained silicates may have been introduced by impact injection into ureilites during this period. About 57% of the ureilites were unchanged after period 2. During period 3 events, impact-induced annealing caused previously mosaicized olivine grains to become aggregates of small unstrained crystals. Some ureilites experienced reduction as FeO at the edges of olivine grains reacted with C from the matrix. Annealing may also be responsible for coarsening of graphite in a few ureilites, forming euhedral-appearing, idioblastic crystals. Orthopyroxene in Meteorite Hills (MET) 78008 may have formed from pigeonite by annealing during this period. The Rb-Sr internal isochron age of approx.4.0 Ga for MET 78008 probably dates the annealing event. At this late date, impacts are the only viable heat source. About 36% of ureilites experienced period 3 events, but remained unchanged afterwards. During period 4, approx.7% of the ureilites were shocked again, as is evident in the polymict breccia, Elephant Moraine (EET) 83309. This rock contains annealed mosaicized olivine aggregates composed of small individual olivine crystals that exhibit undulose extinction. Ureilites may have formed by impact-melting chondritic material on a primitive body with heterogeneous O isotopes. Plagioclase was preferentially lost from the system due to its low impedance to shock compression. Brief melting and rapid burial minimized the escape of planetary-type noble gases from the ureilitic melts. Incomplete separation of metal from silicates during impact melting left ureilites with relatively high concentrations of trace siderophile elements.
Rubin, Alan E.↗
Shock experiments on maskelynite-bearing anorthosite
A series of shock recovery experiments over 9.9-60.4 GPa have been carried out on naturally shocked anorthosite from the Mistastin impact structure in Labrador consisting primarily of diaplectic plagioclase glass or maskelynite, An(50), and pyroxene. Petrographic observations of the experimental products indicate that the component minerals and diaplectic glasses generally retained their initial character throughout, the only exception being the increase in fracturing which occurred in the 9.9 GPa shot. Reshocking at pressures higher than the initial shock tends to lower the refractive index of maskelynite. The increase in refractive index of maskelynite reshocked to pressures lower than the initial pressure is interpreted as due to shock densification of the diaplectic glass above the Hugoniot elastic limit and below the mixed phase regime. The data suggest that the low-high-low density transition of maskelynite occurs about 8 GPa below that of the crystal of corresponding composition.
Lambert, P.↗
Prediction and measurement of heat transfer rates for the shock-induced unsteady laminar boundary layer on a flat plate
The unsteady laminar boundary layer induced by the flow-initiating shock wave passing over a flat plate mounted in a shock tube was theoretically and experimentally studied in terms of heat transfer rates to the plate for shock speeds ranging from 1.695 to 7.34 km/sec. The theory presented by Cook and Chapman for the shock-induced unsteady boundary layer on a plate is reviewed with emphasis on unsteady heat transfer. A method of measuring time-dependent heat-transfer rates using thin-film heat-flux gages and an associated data reduction technique are outlined in detail. Particular consideration is given to heat-flux measurement in short-duration ionized shocktube flows. Experimental unsteady plate heat transfer rates obtained in both air and nitrogen using thin-film heat-flux gages generally agree well with theoretical predictions. The experimental results indicate that the theory continues to predict the unsteady boundary layer behavior after the shock wave leaves the trailing edge of the plate even though the theory is strictly applicable only for the time interval in which the shock remains on the plate.
Cook, W. J.↗
The size distributions of fragments ejected at a given velocity from impact craters
The mass distribution of fragments that are ejected at a given velocity for impact craters is modeled to allow extrapolation of laboratory, field, and numerical results to large scale planetary events. The model is semi-empirical in nature and is derived from: (1) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) observed ejecta blanket particle size distributions, (3) an empirical relationship between maximum ejecta fragment size and crater diameter, (4) measurements and theory of maximum ejecta size versus ejecta velocity, and (5) an assumption on the functional form for the distribution of fragments ejected at a given velocity. This model implies that or planetary impacts into competent rock, the distribution of fragments ejected at a given velocity is broad, e.g., 68% of the mass of the ejecta at a given velocity contains fragments having a mass less than 0.1 times a mass of the largest fragment moving at that velocity. The broad distribution suggests that in impact processes, additional comminution of ejecta occurs after the upward initial shock has passed in the process of the ejecta velocity vector rotating from an initially downward orientation. This additional comminution produces the broader size distribution in impact ejecta as compared to that obtained in simple brittle failure experiments.
Okeefe, J. D.↗
The size distributions of fragments ejected at a given velocity from impact craters
The mass distribution of fragments that are ejected at a given velocity for impact craters is modeled to allow extrapolation of laboratory, field, and numerical results to large scale planetary events. The model is semi-empirical in nature and is derived from: (1) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) observed ejecta blanket particle size distributions, (3) an empirical relationship between maximum ejecta fragment size and crater diameter, (4) measurements and theory of maximum ejecta size versus ejecta velocity, and (5) an assumption on the functional form for the distribution of fragments ejected at a given velocity. This model implies that for planetary impacts into competent rock, the distribution of fragments ejected at a given velocity is broad, e.g., 68 percent of the mass of the ejecta at a given velocity contains fragments having a mass less than 0.1 times a mass of the largest fragment moving at that velocity. The broad distribution suggests that in impact processes, additional comminution of ejecta occurs after the upward initial shock has passed in the process of the ejecta velocity vector rotating from an initially downward orientation. This additional comminution produces the broader size distribution in impact ejecta as compared to that obtained in simple brittle failure experiments.
O'Keefe, John D.↗
Molecule formation and infrared emission in fast interstellar shocks. I Physical processes
The paper analyzes the structure of fast shocks incident upon interstellar gas of ambient density from 10 to the 7th per cu cm, while focusing on the problems of formation and destruction of molecules and infrared emission in the cooling, neutral post shock gas. It is noted that such fast shocks initially dissociate almost all preexisting molecules. Discussion covers the physical processes which determine the post shock structure between 10 to the 4 and 10 to the 2 K. It is shown that the chemistry of important molecular coolants H2, CO, OH, and H2O, as well as HD and CH, is reduced to a relatively small set of gas phase and grain surface reactions. Also, the chemistry follows the slow conversion of atomic hydrogen into H2, which primarily occurs on grain surfaces. The dependence of this H2 formation rate on grain and gas temperatures is examined and the survival of grains behind fast shocks is discussed. Post shock heating and cooling rates are calculated and an appropriate, analytic, universal cooling function is developed for molecules other than hydrogen which includes opacities from both the dust and the lines.
Hollenbach, D.↗
Cascading economic losses from port disruptions under capacity constrained multimodal freight networks
This study quantifies how throughput disruptions at major seaports cascade through capacity-constrained multimodal freight networks and interregional production systems. We couple an agent-based model (ABM) multimodal freight simulation that resolves rerouting, terminal queueing, and inventory drawdown under binding modal and facility capacities with a multiregional output loss input-output (MRIIM) model that propagates realized delivery shortfalls across regions and sectors. The framework is demonstrated for the Port of Los Angeles using Freight Analysis Framework flows and Bureau of Economic Analysis input-output accounts and is evaluated over a 52-week horizon under deterministic sector targeted shocks and stochastic disruption realizations with uncertain severity and duration. Results indicate nonlinear amplification: realized national losses concentrate in manufacturing and transportation/warehousing even when exogenous port shocks are dispersed, suggesting that congestion spillback and limited short-run substitution can dominate the initial shock allocation. We further evaluate a tabular reinforcement-learning (Q-learning) intervention layer that selects among a small set of implementable system level levers (truck-to-rail and truck-to-barge shift settings) without overriding shipper routing, finding that such interventions reduce total losses for moderate disruptions but yield diminishing returns once substitute modes approach capacity. By linking operational freight behavior to system wide impacts under uncertainty, the proposed ABM-MRIIM pipeline provides a reusable workflow for port disruption stress testing, identification of structurally critical sectors/corridors, and evaluation of resilience interventions under realistic capacity limits.
42 ENGINEERING↗
Configurations of corotating shocks in the outer heliosphere
Configurations of corotating shocks in the outer heliosphere are computed on the assumption that the shocks move at constant speeds between collisions. The basic physical process is the interaction of a forward shock with a reverse shock. The points at which this interaction takes place are determined algebraically in terms of the initial shock positions and speeds. Between the points at which the shocks interact, the shocks have the form of Archimedian spirals. Shock configurations are determined for the case of two corotating shock pairs originating at 2.5 AU and for a single shock pair at 2.5 AU. In both cases the heliosphere is divided into several distinct regions distinguished by the number of times the particles have passed through a shock. Since the shock strengths decrease when a forward shock interacts with a reverse shock, corotating shocks are expected to be weaker at larger distances from the sun.
Burlaga, L. F.↗
An approximate viscous shock layer technique for calculating chemically reacting hypersonic flows about blunt-nosed bodies
An approximate axisymmetric method was developed which can reliably calculate fully viscous hypersonic flows over blunt nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. The approach is applicable to perfect gas, equilibrium, and nonequilibrium flowfields. Since the method is fully viscous, the problems associated with a boundary layer solution with an inviscid layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed before extending the method to nonequilibrium calculations. Perfect gas (laminar and turbulent), equilibrium, and nonequilibrium solutions were generated for airflows over several analytic body shapes. Surface heat transfer, skin friction, and pressure predictions are comparable to VSL results. In addition, computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher order procedures which require starting solutions.
Cheatwood, F. Mcneil↗
An approximate viscous shock layer approach to calculating hypersonic flows about blunt-nosed bodies
An approximate axisymmetric method has been developed which can reliably calculate fully viscous hypersonic flows over blunt-nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. Since the method is fully viscous, the problems associated with coupling a boundary-layer solution with an inviscid-layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed. Surface heat transfer and pressure predictions are comparable to both VSL results and experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher-order procedures which require starting solutions.
Cheatwood, F. MCN.↗
An approximate viscous shock layer technique for calculating nonequilibrium hypersonic flows about blunt-nosed bodies
An approximate axisymmetric method has been developed which can reliably calculate nonequilibrium fully viscous hypersonic flows over blunt-nosed bodies. By substituting Maslen's second-order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. This procedure is significantly faster than the parabolized Navier-Stokes and VSL solvers and would be useful in a preliminary design environment. Solutions have been generated for air flows over several analytic body shapes. Surface heat transfer and pressure predictions are comparable to VSL results. Computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher-order procedures which require starting solutions.
Cheatwood, F. M.↗