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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗
NASA NTRS1992

Heat shield erosion in a dusty Martian atmosphere

The study calculates the effects of dust particle impacts on the erosion of the forebody heat shield for a 26-m diameter aerobraking vehicle entering the Martian atmosphere at 8.6 km/s. An explicit thin-layer Navier-Stokes code is used to compute the dustless flowfield about the vehicle for the actual Martian atmospheric composition. The deceleration and melting of the dust particles within the forebody shock layer was computed for dust spherules having initial diameters from 1 to 19 microns. All the particles began vaporizing shortly after entering the shock layer, but most survived to hit the heat shield surface. The two different heat shield materials considered were Shuttle ceramic tiles with glassy surfaces and a moderate-density ablator. For a vehicle with a ballistic coefficient of 200 kg/sq m, a glassy heat shield surface experienced an average of 7 mm of surface erosion, or 20 times the thickness of the glassy layer on ceramic tiles.

Papadopoulos, Periklis↗
DOE OSTI2021

Improved modeling of the solid-to-plasma transition of polystyrene ablator for laser direct-drive inertial confinement fusion hydrocodes

We report that the target performance of laser direct-drive inertial confinement fusion (ICF) can be limited by the development of hydrodynamic instabilities resulting from the nonhomegeneous laser absorption at the target surface, i.e., the laser imprint on the ablator. To understand and describe the formation of these instabilities, the early ablator evolution during the laser irradiation should be considered. In this work, an improved modeling of the solid-to-plasma transition of a polystyrene ablator for laser direct-drive ICF is proposed. This model is devoted to be implemented in hydrocodes dedicated to ICF which generally assume an initial plasma state. The present approach consists of the two-temperature model coupled to the electron, ion and neutral dynamics including the chemical fragmentation of polystyrene. The solid-to-plasma transition is shown to significantly influence the temporal evolution of both free electron density and temperatures, which can lead to different shock formation and propagation compared with an initial plasma state. The influence of the solid-to-plasma transition on the shock dynamics is evidenced by considering the scaling law of the pressure with respect to the laser intensity. The ablator transition is shown to modify the scaling law exponent compared with an initial plasma state.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Weak shock propagation through a turbulent atmosphere

Consideration is given to the propagation through turbulence of transient pressure waveforms whose initial onset at any given point is an abrupt shock. The work is motivated by the desire to eventually develop a mathematical model for predicting statistical features, such as peak overpressures and spike widths, of sonic booms generated by supersonic aircraft. It is argued that the transient waveform received at points where x greater than 0 will begin with a pressure jump and a formulation is developed for predicting the amount of this jump and the time derivatives of the pressure waveform immediately following the jump.

Pierce, Allan D.↗
NASA NTRS2006

Solar Sentinels: Report of the Science and Technology Definition Team

The goal of NASA s Living With a Star (LWS) program is to develop the scientific understanding necessary to effectively address those aspects of the connected Sun Earth system that directly affect life and society. Along with the other elements of LWS, Solar Sentinels aims to discover, understand, and model the heliospheric initiation, propagation, and solar connection of those energetic phenomena that adversely affect space exploration and life and society here on Earth. The Solar Sentinels mission will address the following questions: (1) How, where, and under what circumstances are solar energetic particles (SEPs) accelerated to high energies and how do they propagate through the heliosphere? And (2) How are solar wind structures associated with these SEPs, like CMEs, shocks, and high-speed streams, initiated, propagate, evolve, and interact in the inner heliosphere? The Sentinels STDT recommends implementing this mission in two portions, one optimized for inner heliospheric in-situ measurements and the other for solar remote observations. Sentinels will greatly enhance the overall LWS science return.

Source record↗
NASA NTRS2004

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗
NASA NTRS1985

Shock drift acceleration in the presence of waves

Attention is given to the initial results of a model designed to study the modification of the scatter-free, shock drift acceleration of energetic test particles by wave activity in the vicinity of a quasi-perpendicular, fast-mode MHD shock. It is emphasized that the concept of magnetic moment conservation is a valid approximation only in the perpendicular and nearly perpendicular regimes, when the angle theta-Bn between the shock normal and the upstream magnetic field vector is in the range from 70 deg to 90 deg. The present investigation is concerned with one step in a program which is being developed to combine the shock drift and diffusive processes at a shock of arbitrary theta-Bn.

Decker, R. B.↗
DOE OSTI2023

Statistical characterization of a shock interacting with an inclined gas column

We present a statistical characterization of the interaction between a planar shock and a finite-diameter, cylindrical column of dense gas based on three-dimensional, large-eddy simulation results. In the simulation, the column of gas is initially inclined at an angle $\alpha _0$ with respect to the shock plane. Effects of the initial column angle on the mixing characteristics are examined at Mach number 2.0 for column incline angles $1^\circ$ , $5^\circ$ , $10^\circ$ and $30^\circ$ . Mean velocity profiles show that the column angle affects the gas velocity components in the vertical plane, but not in the spanwise direction. The gas undergoes higher initial upward acceleration at larger initial column incline angles. With time, the gas motion tends to become one-dimensional in the streamwise direction. Initially, velocity fluctuations are most intense within the interior of the column, but concentrate near the column leading edge over time. At high wavenumbers $\kappa$ , the turbulent kinetic energy spectra follow a power-law scaling of $\kappa ^{-1}$ . The structure functions of the mass fraction do not clearly demonstrate power-law scaling except at early times for $\alpha _0=30^\circ$ , manifesting overall trends very similar to those observed in earlier experiments. Probability distributions of the mass fraction show independence of the mean and the standard deviation of the mixed gas on $\alpha _0$ . The column angle was also found to have little effect on the mixing efficiency characterized by the molecular mixedness. Velocity components in the streamwise and transverse directions tend towards a bimodal distribution for larger $\alpha _{0}$ .

42 ENGINEERING↗
NASA NTRS1992

Ion acceleration at collisionless shock interactions

When two collisionless shocks collide, a population of high energy ions is produced. Using hybrid numerical simulations, it is shown that when the two shocks are quasi-perpendicular, ions with energies up to 15 E sub 0 (where E sub 0 is the initial kinetic energy of upstream ions in the shock frame) are produced by direct electric field acceleration. In the quasi-parallel regime, energies as high as 80 E sub 0 can be obtained. These ions are accelerated by both scattering off the approaching shocks and subsequently in the intense turbulence left behind by the shock collision.

Cargill, Peter J.↗
NASA NTRS1995

Distances to the termination shock and heliopause from a simulation analysis of the 1992-93 heliospheric radio emission event

A new heliospheric radio emission event observed by Voyagers 1 and 2 in mid-1992 is believed to have been produced by the interaction of an interplanetary shock with the heliopause. The shock is thought to have oriented near the Sun during a period of intense solar activity in late-May and early-June, 1991. The observed travel time of the shock to the heliopause is 408 days; the initial speed is estimated to be between 600 and 800 km/s. We use a numerical gasdynamic simulation of an interplanetary shock, propagating through an equilibrium solution of the solar wind/interstellar medium interaction, to compute the distances to the termination shock and the heliopause that are consistent with these observations. For a shock speed of 600 km/s, the termination shock is located at 92 AU, and the heliopause is located at 128 AU. These distances increase to 112 AU and 156 AU when the shock speed is increased to 800 km/s.

Steinolfson, R. S.↗
NASA NTRS1994

Focused interplanetary transport of approximately 1 MeV solar energetic protons through self-generated Alfven waves

We present a model of the focused transport of approximately 1 MeV solar energetic protons through interplanetary Alfven waves that the protons themselves amplify or damp. It is based on the quasi-linear theory but with a phenomenological pitch angle diffusion coefficient in the 'resonance gap.' For initial Alfven wave distributions that give mean free paths greater than approximately 0.5 AU for approximately 1 MeV protons in the inner heliosphere, the model predicts greater than roughly an order of magnitude amplification (damping) in the outward (inward) propagating resonant Alfven waves at less than or approximately equal to o.3 AU heliocentric distance. As the strength of proton source is increased, the peak differential proton intensity at approximately 1 MeV at 1 AU increases to a maximum of approximately 250 particles (/(sq cm)(s)(sr)(MeV)) and then decreases slowly. It may be attenuated by a factor of 5 or more relative to the case without wave evolution, provided that the proton source is sufficiently intense that the resulting peak differential intensity of approximately 1 MeV protons at 1 AU exceeds approximately 200 particles (/(sq cm)(s)(sr)(MeV)). Therefore, in large solar proton events, (1) one may have to take into account self-amplified waves in studying solar particle propagation, (2) the number of accelerated protons escaping from a flare or interplanetary shock may have been underestimated in past studies by a significant factor, and (3) accelerated protons escaping from a traveling interplanetary shock at r less than or approximately equal to 0.3 AU should amplify the ambient hydromagnetic waves siginificantly to make the shock an efficient accelerator, even if initially the mean free path is greater than or approximately equal to 1 AU.

Ng, C. K.↗
NASA NTRS1981

Stall flutter experiment in a transonic oscillating linear cascade

Two dimensional biconvex airfoils were oscillated at reduced frequencies up to 0.5 based on semi-chord and a free stream Mach number of 0.80 to simulate transonic stall flutter in rotors. Steady-state periodicity was confirmed through end-wall pressure measurements, exit flow traverses, and flow visualization. The initial flow visualization results from flutter tests indicated that the oscillating shock on the airfoils lagged the airfoil motion by as much as 80 deg. These initial data exhibited an appreciable amount of scatter; however, a linear fit of the results indicated that the greatest shock phase lag occurred at a positive interblade phase angle. Photographs of the steady-state and unsteady flow fields reveal some of the features of the lambda shock wave on the suction surface of the airfoils.

Boldman, D. R.↗
DOE OSTI2021

Particle Curtain Experiments at Sandia

Efforts at Sandia National Laboratories have focused on fundamental experiments to understand the dispersal of dense particle distributions in high-speed compressible flow. The experiments are conducted in shock tube facilities where the flow conditions and the initial conditions of the particle distributions are well controlled and well characterized. An additional advantage of the shock tube is that it is more readily able to accommodate advanced measurement diagnostics in comparison to explosive field tests.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Pressure distributions and shock shapes for a bent-nose biconic at incidence

Initial results of a NASA study of the lift-drag characteristics of a 12.84/7 deg biconic model intended for airbraking during atmospheric entry of probes to Mars, Venus, Saturn, and Titan are reported. Pressure distributions and shock shapes were measured in the Langley 20 in. Mach 6 tunnel with the spherically blunted bent-nose model set at angles from 0-25 deg. Pressure distributions and shock shapes where computed using the STEIN flowfield code, which features a MacCormack scheme to integrate the three-dimensional Euler equations, the Rankine-Hugoniot jump conditions to model shock waves as discontinuities, and requires a supersonic condition at every step. A comparison was made between measured and predicted values. The leeward shock angle was found to be predictable to within 3% for all angles of attack, while parabolized Navier-Stokes equations are regarded as offering more accurate results than the STEIN code for surface pressure distributions.

Miller, C. G., III↗
DOE OSTI2021

Post-explosion Evolution of Core-collapse Supernovae

Here we investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 explosion models during several seconds. We present a broad study based on three progenitors (11.2, 15, and 27 M ⊙ ), different neutrino heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind (NDW) that survives for several seconds. This indicates that NDWs are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 M ⊙ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.

79 ASTRONOMY AND ASTROPHYSICS↗