SEARCH · Engineering Papers

Results for “shock jumps”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Effects of highly intermittent magnetic fields and plasma flows on solar wind dynamics

Filamentary magnetic fields and intermittent mass flows with highly variant physical parameters as observed in coronal holes provide, from a theoretical point of view, natural conditions for strongly nonlinear dynamics. The presence of sheared mass flows along fine scale magnetic structures results in strong nonlinear instability, most important of which is the explosive instability. We specify the physical conditions for several different manifestations of the onset of explosive instability and its further evolution: (1) fully developed explosive instability - explosive release of the energy; (2) shock formation - stabilization of instability by small scale spatial inhomogeneities leads to formation of subsequent shocks having a number of peculiarities that is determined by the interplay of thermal and viscous losses (for example, in predominance of thermal losses the isothermal jump occurs); and (3) solitary waves - stabilization of explosive instability by nonlinear dispersion effect leads to formation of a 'gas' of solitons which are later either damped away with characteristic time and energy input or evolve to solitons with explosively growing amplitudes. Each scenario is completely determined by the physical parameters of the medium, thus producing a quite uneven distribution of energy in a coronal hole and, respectively, an uneven outward propagation of the energy flux.

Ryutova, M.↗

NASA NTRS1978

Solar wind stream interfaces

Results are presented for a superposed epoch analysis of discontinuous solar wind interfaces. The average time-space profiles of stream interfaces are discussed with reference to fluid properties (flow speed, pressure ridge, density, electron and proton temperatures) and kinetic properties (electron core and halo, flow speed fluctuations, electron heat flux, alpha particles). Other aspects of stream interfaces are described, such as the persistence of individual interfaces, shock associations, the sector boundaries of the interplanetary magnetic field, and sudden impulses in the geomagnetic field. Interface position is considered in terms of the observed temperature jump. A conceptual model of high-speed stream evolution is proposed.

Gosling, J. T.↗

NASA NTRS2002

Numerical Investigation of 'Transonic Resonance' with a Convergent-Divergent Nozzle

At pressure ratios lower than the design value, convergent-divergent (C-D) nozzles often undergo a flow resonance accompanied by the emission of acoustic tones. The phenomenon, driven by the unsteady shock within the divergent section of the nozzle, has been studied experimentally by Zaman et al. In this paper, the space-time conservation element solution element (CE/SE) method is employed to numerically investigate the phenomenon. The computations are performed for a given nozzle geometry for several different pressure ratios. Sustained 'limit cycle' oscillations are encountered in all cases. The oscillation frequencies, their variation with pressure ratio including a 'stage jump', agree well with the experimental results. The unsteady flow data confirm that stage 1 of the resonance (fundamental) involves a one-quarter standing wave while stage 2 (third harmonic) involves a three-quarter standing wave within the divergent section of the nozzle. Details of the shock motion, and the flow and near acoustic field, are documented for one case each of stages 1 and 2.

Loh, Ching Y.↗

NASA NTRS2021

Magnetic Field and Plasma Density Observations of a Pressure Front by Voyager 1during 2020 in the Very Local Interstellar Medium

Voyager 1 has been moving through the very local interstellar medium (VLISM) from the time that it crossed the heliopause on 2012/DOY 238 to 2020/DOY 292. Three notable objects in the magnetic field of the VLISM have been observed: two shocks and one pressure front. This paper reports the observation of a fourth object observed near 2020/DOY 147. There were no upstream electron plasma oscillations of the type often observed ahead of shocks, abrupt increases in energetic particles, or fluctuations in the 48 s increments of the magnetic field associated with this feature, suggesting that it was probably not a shock. This feature was associated with a relatively large increase in the magnetic field strength(B2/B1=1.35)and in the electron density determined by the Plasma Wave Science experiment (N2/N1=1.36) using a new method described in this paper. This feature appears to be a pressure front associated with a compressive wave in the VLISM. The two shocks and the two pressure fronts were associated with the four largest maxima observed in B (t) between 2012/DOY 238 and 2020/DOY 292. Each feature was associated with a jump–ramp structure. The jump–ramp structures were separated by long relatively undisturbed quiet intervals.

Interstellar magnetic fields↗

NASA NTRS1983

Quasi-steady solar wind dynamics

Progress in understanding the large scale dynamics of quasisteady, corotating solar wind structure was reviewed. The nature of the solar wind at large heliocentric distances preliminary calculations from a 2-D MHD model are used to demonstrate theoretical expectations of corotating structure out to 30 AU. It is found that the forward and reverse shocks from adjacent CIR's begin to interact at about 10 AU, producing new shock pairs flanking secondary CIR's. These sawtooth secondary CIR's interact again at about 20 AU and survive as visible entities to 30 AU. The model predicts the velocity jumps at the leading edge of the secondary CIR's at 30 AU should be very small but there should still be sizable variations in the thermodynamic and magnetic parameters. The driving dynamic mechanism in the distant solar wind is the relaxation of pressure gradients. The second topic is the influence of weak, nonimpulsive time dependence in quasisteady dynamics. It is suggested that modest large scale variations in the coronal flow speed on periods of several hours to a day may be responsible for many of the remaining discrepancies between theory and observation. Effects offer a ready explanation for the apparent rounding of stream fronts between 0.3 and 1.0 AU discovered by Helios.

Pizzo, V. J.↗

NASA NTRS1989

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗

DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗

NASA NTRS1966

Nylon shock absorber prevents injury to parachute jumpers

Nylon shock absorbers reduce the canopy-opening shock of a parachute to a level that protects the wearer from injury. A shock absorber is mounted on each of the four risers between the shroud lines and the harness. Because of their size and location, they pose no problem in repacking the chute and harness after a jump.

Mandel, J. A.↗

DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1986

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗

NASA NTRS2017

Feature Detection and Curve Fitting Using Fast Walsh Transforms for Shock Tracking: Applications

Walsh functions form an orthonormal basis set consisting of square waves. Square waves make the system well suited for detecting and representing functions with discontinuities. Given a uniform distribution of 2p cells on a one-dimensional element, it has been proven that the inner product of the Walsh Root function for group p with every polynomial of degree < or = (p - 1) across the element is identically zero. It has also been proven that the magnitude and location of a discontinuous jump, as represented by a Heaviside function, are explicitly identified by its Fast Walsh Transform (FWT) coefficients. These two proofs enable an algorithm that quickly provides a Weighted Least Squares fit to distributions across the element that include a discontinuity. The detection of a discontinuity enables analytic relations to locally describe its evolution and provide increased accuracy. Time accurate examples are provided for advection, Burgers equation, and Riemann problems (diaphragm burst) in closed tubes and de Laval nozzles. New algorithms to detect up to two C0 and/or C1 discontinuities within a single element are developed for application to the Riemann problem, in which a contact discontinuity and shock wave form after the diaphragm bursts.

Gnoffo, Peter A.↗

NASA NTRS1990

Model atmospheres for the early stages of novae in outburst

For the photospheres of classical novae during the early stages of their outbursts, continuum and line blanketed models are presented. The expanding envelopes are characterized by small density gradients that lead to very large geometrical extensions and large temperature differences between inner and outer parts. The spectra show large excesses in the IR and small Balmer jumps which may either be in absorption or emission. For the parameters considered most lines are in absorption if powerlaw density distributions and radiative equilibrium are assumed. The effects of modifications in the temperature structure (e.g., heating from shock fronts) and in the density stratifications as well as changes in the abundances are discussed.

Wehrse, R.↗

NASA NTRS1974

Observation of a driver gas-tangential discontinuity

A complete analysis of an interplanetary disturbance of Nov. 19, 1970 using the Apollo 12-SIDE (Suprathermal Ion Detector Experiment) is presented. The SIDE detectors were pointing at 26.3 degrees from the normal solar-wind direction during the observations. The data were least-squares fitted (using a parabolic hypersurface approximation) to a convected Maxwell-Boltzmann distribution function. The results of the fit combined with two other experiments showed a drastic change in the wind speed (from an steady 352 km/sec down to 219 km/sec), direction, and temperature. Except for a delta-function increase at the onset, the density remained constant. There was a considerable enhancement in the abundance of He and probably of heavier elements. The interplanetary magnetic field exhibited a jump of 21 gamma with a change in latitude from -56 to -76 degrees in solar ecliptic coordinates. It is concluded that the disturbance was due to the driver gas-tangential discontinuity of a solar flare-induced shock wave. The characteristic of the tangential discontinuity fit well with theoretical prediction.

Medrano, R. A.↗

NASA NTRS1988

Ion equation of state in quasi-parallel shocks - A simulation result

Ion equation of state in the quasi-parallel collisionless shock is deduced from simulation results. The simulations were performed for theta(bn) = 10 deg, beta = 0.5 and M sub A in the range from 1.2 to 8, where M sub A is the Alfven Mach number, beta is the upstream ratio of plasma pressure to magnetic pressure, and theta(bn) is the angle between the shock normal and the upstream magnetic field. The equation of state can be approximated by a power law with different exponents in the upstream and downstream sides of the shock transition region. The exponent in the upstream side of the transition region is much greater than the adiabatic value of 5/3 and increases with M sub A. The exponent in the downstream side of the transition region is slightly less than 5/3. The results show that ion heating in the quasi-parallel shock is highly nonadiabatic with a large increase in entropy and in temperature ratio in the upstream side of the transition region, while the heating is highly isentropic with a large increase in temperature difference across the principal density jump in the downstream side of the transition region.

Mandt, M. E.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

Vijayaraghavan, D.↗

DOE OSTI2021

Understanding the evolution of liquid and solid microjets from grooved Sn and Cu samples using radiography

Experiments were performed on grooved Sn and Cu samples to study the temporal evolution of microjets. Jets were generated by the impact of gun-launched flyer plates against the back of grooved targets made from either Cu or Sn (groove depth of ~250 μm). The Hugoniot states in the various Sn targets encompassed conditions where solid phases are maintained throughout (7 and 16 GPa) and also conditions where melting occurs upon the release of compression (25 and 34 GPa); the transition occurs near a Hugoniot pressure of 23 GPa. Cu targets at 27 and 56 GPa provide comparisons in which the jets move at similar speeds but remain solid. In all cases, the spatial distribution of mass within the microjets was measured using high-speed synchrotron radiography. The result is a time history of the jet thickness profile from which quantities like total jet mass and jet velocity can be derived. In both the solid and liquid states, we generally observe that an increase in the shock strength leads to an increase in jet mass. However, this trend breaks down for Hugoniot states near the transition from continuously solid to melted-on-release. This is evidenced by the observation that there was no difference in the rate of mass flow in Sn jets at 16 and 25 GPa, while similar pressure jumps on either side of this range caused substantial changes in the jet mass. This contrasts with the behavior of smaller polishing defects that were present on the same samples (~1 μm deep). From these, no ejecta mass was detected below the melt boundary, but obvious microjets were generated once melting occurred. This indicates that crossing the bulk melt-on-release threshold can alternately promote or inhibit the flow of mass into microjets based on the amplitude of the initial perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗