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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Dayside Proton Aurora: Comparisons between Global MHD Simulations and Image Observations

The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the iMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model. though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside ionosphere are enhanced. Overall the results of these initial comparisons between global MHD simulation results and IMAGE observations emphasize the interplay between reconnection and dynamic pressure processes at the dayside magnetopause. as well as the intricate connection between the bow shock and the auroral region.

Berchem, J.↗

NASA NTRS1990

Numerical simulations of the rebound shock model for solar spicules

Using time-dependent numerical simulations, the proposed rebound shock mechanism for spicules has been examined. At temperatures above a critical value, T(c), the radiation is characteristic of the conditions in the optically thin corona and near optically thin transition region. When T less than T(c), the atmosphere has a radiative cooling time, tau(rad) characteristic of chromosphere. The spicule is initiated with a quasi-impulsive force in the low chromosphere, which drives a train of upward propagating rebound shocks along the rigid magnetic flux tube. These shocks then move the transition region upward. The material below the displaced transition region has temperatures and densities similar to those of spicules when T(c) = 20,000 K or more and tau(rad) = 500 s or more, but not when T(c) = 10,000 K, and probably not when tau(rad) = 100 s. For all the cases where the cross sectional area diverges rapidly with height, the upward velocity of the transition region is less than that of spicules. Moreover, the maximum height is less than that of average spicules. Taller, higher velocity spicules result when the magnetic field cross sectional area is constant. In all cases, the rebound shock mechanism produces substantial motions and temperature and density variations in chromospheric and transition region material. It is suggested that this may be a partial explanation for the continuous dynamic state of the lower solar atmosphere.

Sterling, Alphonse C.↗

NASA NTRS2001

Vorticity Transport and Wave Emission in the Protoplanetary Nebula

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler/Energy equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the nebula to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived using the simulation data and compared with analytical dispersion relations from the linearized Euler/Energy equations.

Davis, S. S.↗

NASA NTRS2002

Vorticity Transport and Wave Emission In A Protoplanetary Disk

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the disk to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived from the nonlinear simulation data and correlated with analytical dispersion relations from the linearized Euler and energy equations.

Davis, S. S.↗

NASA NTRS1961

Concerning Certain Collisionless Plasma-Shock Wave Models

A shock-wave model proposed by both Kahn and Parker for the collision of two interstellar gas clouds or of the solar wind with the ionized gas around the earth is re-evaluated in the light of recent developments in plasma stability theory. It is shown that the instability which was supposed to arrest the counter-streaming of the ions and to transfer much of their kinetic energy to the electrons does not occur in many important examples. Specifically, it is absent if the initial ion thermal energy exceeds 4.2 x 10 (sup -4) times their translational energy. While other plausible shock structures are known which yield supra-thermal electrons, they depend on magnetic fields, leaving no satisfactory theory in the limit of zero magnetic field.

Noerdlinger, Peter D.↗

NASA NTRS2005

Evidence for a Canonical GRB Afterglow Light Curve in the Swift/XRT Data

We present new observations of the early X-ray afterglows of the first 27 gamma-ray bursts (GRBs) detected with the Swift X-ray Telescope (XRT). The early X-ray afterglows show a canonical behavior, where the light curve broadly consists of three distinct power law segments. These power law segments are separated by two corresponding break times. On top of this canonical behavior of the early X-ray light curve, many events have superimposed X-ray flares, which are most likely caused by internal shocks due to long lasting sporadx activity of the central engine, up to several hours after the GRB. We find that the initial steep decay is consistent with it being the tail of the prompt emission: from photons that are radiated at large angles relative to our line of sight. The first break in the light curve takes place when the forward shock emission becomes dominant, with the intermediate shallow flux decay likely caused by the continuous energy injection into the external shock. When this energy injection stops, a second break is then observed in the light curve. This energy injection increases the energy of the afterglow shock by at least a factor of f greater than or approx. equal to 4, and augments the already severe requirements for the efficiency of the prompt gamma-ray emission.

Nousek, J. A.↗

NASA NTRS1992

Kinetic structure of slow shocks - Effects of the electromagnetic ion/ion cyclotron instability

The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.

Omidi, N.↗

DOE OSTI2021

1D kinetic study of pinch formation in a dense plasma focus: Transition from collisional to collisionless regimes

Here, the pinch-formation stage of a deuterium dense plasma focus, and associated “shock-flash” neutron yield, is studied using 1D kinetic simulations considering a plasma column with initial pressure P, initial radius R, and the compression to be driven by a constant current I. The relative behavior of the compression is shown to be similar for fixed ratios of the characteristic ion mean free path to the radius of the plasma column at stagnation, $λ_{st}/R_{st}$. This dimensionless parameter is shown to scale like $I^4/(P^3R^5)$. The compression ratio, $R/R_{st}$, is found to be a minimum when $λ_{st}/R_{st} ≈$ 1 and is the largest in the collisionless limit where $λ_{st} \gg R_{st}$. This behavior is in contrast to the analogous planar pinch where R/Rst decreases from one constant for λst/Rst$\ll$1 to a smaller constant for $λ_{st}/R_{st} \gg$ 1. The yield in the collisionless regime is shown to fall between the two well-known $I^4$ scaling laws. Furthermore, this regime exhibits qualities that potentially make it appealing for radiography applications, such as increased localization in time and space of the neutron formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2008

CAWSES November 7-8, 2004, Superstorm: Complex Solar and Interplanetary Features in the Post-Solar Maximum Phase

The complex interplanetary structures during 7 to 8 Nov 2004 are analyzed to identify their properties as well as resultant geomagnetic effects and the solar origins. Three fast forward shocks, three directional discontinuities and two reverse waves were detected and analyzed in detail. The three fast forward shocks 'pump' up the interplanetary magnetic field from a value of approx.4 nT to ~44 nT. However, the fields after the shocks were northward, and magnetic storms did not result. The three ram pressure increases were associated with major sudden impulses (SI + s) at Earth. A magnetic cloud followed the third forward shock and the southward Bz associated with the latter was responsible for the superstorm. Two reverse waves were detected, one at the edge and one near the center of the magnetic cloud (MC). It is suspected that these 'waves' were once reverse shocks which were becoming evanescent when they propagated into the low plasma beta MC. The second reverse wave caused a decrease in the southward component of the IMF and initiated the storm recovery phase. It is determined that flares located at large longitudinal distances from the subsolar point were the most likely causes of the first two shocks without associated magnetic clouds. It is thus unlikely that the shocks were 'blast waves' or that magnetic reconnection eroded away the two associated MCs. This interplanetary/solar event is an example of the extremely complex magnetic storms which can occur in the post-solar maximum phase.

solar flares↗

DOE OSTI2020

Effects of High Fuel Loading and CO2 Dilution on Oxy-Methane Ignition Inside a Shock Tube at High Pressure

Abstract Ignition delay times were measured for methane/O2 mixtures in a high dilution environment of either CO2 or N2 using a shock tube facility. Experiments were performed between 1044 K and 1356 K at pressures near 16 ± 2 atm. Test mixtures had an equivalence ratio of 1.0 with 16.67% CH4, 33.33% O2, and 50% diluent. Ignition delay times were measured using OH* emission and pressure time-histories. Data were compared to the predictions of two literature kinetic mechanisms (ARAMCO MECH 2.0 and GRI Mech 3.0). Most experiments showed inhomogeneous (mild) ignition which was deduced from five time-of-arrival pressure transducers placed along the driven section of the shock tube. Further analysis included determination of blast wave velocities and locations away from the end wall of initial detonations. Blast velocities were 60–80% of CJ-Detonation calculations. A narrow high temperature region within the range was identified as showing homogenous (strong) ignition which showed generally good agreement with model predictions. Model comparisons with mild ignition cases should not be used to further refine kinetic mechanisms, though at these conditions, insight was gained into various ignition behavior. To the best of our knowledge, we present first shock tube data during ignition of high fuel loading CH4/O2 mixtures diluted with CO2 and N2.

Energy & Fuels↗

DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2010

Experimental Evaluation of Fatigue Crack Initiation from Corroded Hemispherical Notches in Aerospace Structural Materials

A test program was developed and executed to evaluate the influence of corroded hemispherical notches on the fatigue crack initiation and propagation in aluminum 7075-T7351, 4340 steel, and D6AC steel. Surface enhancements such as shot peening and laser shock peening were also incorporated as part of the test effort with the intent of improving fatigue performance. In addition to the testing, fracture mechanics and endurance limit based analysis methods were evaluated to characterize the results with the objective of challenging typical assumptions used in modeling fatigue cracks from corrosion pits. The results specifically demonstrate that the aluminum and steel alloys behave differently with respect to fatigue crack initiation from hemispherical corrosion pits. The aluminum test results were bounded by the fracture mechanics and endurance limit models while exhibiting a general insensitivity to the residual stress field generated by shot peening. The steel specimens were better characterized by the endurance limit fatigue properties and did exhibit sensitivities to residual stresses from the shot peening and laser shock peening

Garcia, Daniel B.↗

NASA NTRS1994

Igneous inclusions from ordinary chondrites: High temperature cumulates and a shock melt

We report microprobe, instrumental neutron activation analysis, and radiochemical neutron activation analysis data for three large igneous inclusions in the Yamato (Y-)75097, Y-793241, and Y-794046 ordinary chondrites. The inclusions in the first two chondrites are troctolitic cumulates that have undergone appreciable reactions with their hosts either during emplacement and/or cooling. Olivine-spinel Fe-Mg exchange pairs in these two inclusions record equilibration temperatures of about 710 C, and these temperatures are similar to those exhibited by mineral pairs in the Y-75097 and Y-793241 hosts. The inclusion in Y-794046 is texturally unique, consisting of fine-grained, randomly distributed olivines, coarse (approximately 2 mm) fascicular pyroxene laths, and angular pockets of maskelynite/plagioclase feldspar. The phase compositions are readily interpreted as having resulted from extremely rapid, essentially isochemical cooling to temperatures less than 1000 C of a melt with an initial temperature greater than 1670 C. We suggest that this igneous inclusion formed in-situ by shock.

Sack, Richard O.↗

NASA NTRS2018

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for a range of high speed Earth entry experiments conducted in the Electric Arc Shock Tube at NASA Ames. The nonequilibrium radiative transport equations are solved in NEQAIR using flow field variables from the full facility CFD simulations of the EAST shock tube performed by US3D ow solver. These physics-based flow calculations lead to a significantly different post-shock gas state and associated radiation field as compared to that based on a simplified but computationally inexpensive calculation for flow over a blunt-body with appropriate initial conditions. The radiation spectra and radiance profiles are computed for an extensive range of wavelengths, from deep VUV to IR, which are pertinent to the emission characteristics of high enthalpy shock waves in air. The radiation properties of the shocked gas are calculated both in the nonequilibrium region at the shock, and in the equilibrium region behind the shock. Numerical predictions are found to be consistent with the experimental observations.

Simulations↗

NASA NTRS1985

Conical similarity of shock/boundary-layer interactions generated by swept and unswept fins

A parametric experimental investigation has been made of the class of three-dimensional shock wave/turbulent boundary layer interactions generated by swept and unswept leading-edge fins. The fin sweepback angles were 0-65 deg at 5, 9, and 15 deg angles of attack. Two equilibrium two-dimensional turbulent boundary layers with a freestream Mach number of 2.95 and a Reynolds number of 6.3 x 10 to the 7th/m were used as incoming flow conditions. All of the resulting interactions were found to possess conical symmetry of the surface flow patterns and pressures outside of an initial inception zone. Further, these interactions were found to obey a simple conical similarity rule based on inviscid shock wave strength, irrespective of fin sweepback or angle of attack. This is one of the first demonstrations of similarity among three-dimensional interactions produced by geometrically dissimilar shock generators.

Settles, G. S.↗

NASA NTRS2012

Radiation from Accelerated Particles in Shocks and Reconnections

Plasma instabilities are responsible not only for the onset and mediation of collisionless shocks but also for the associated acceleration of particles. We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electrons transverse deflection and, more generally, relativistic acceleration behind the shock. We have calculated, self-consistently, the radiation from electrons accelerated in the turbulent magnetic fields. We found that the synthetic spectra depend on the Lorentz factor of the jet, its thermal temperature and strength of the generated magnetic fields. Our initial results of a jet-ambient interaction with anti-parallelmagnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in a transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K. I.↗

DOE OSTI2020

Variable-density buoyancy-driven turbulence with asymmetric initial density distribution

The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.

42 ENGINEERING↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗