Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moving shock”

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 55 records · Page 3

Spectrophotometry of the Cygnus Loop

New spectrophotometry for 17 filaments in the Cygnus Loop are presented. A wide variety of spectral properties is observed within the remnant. Comparisons with published shock models indicate significant differences between the models and observations with the range of observed line intensities not well described by models with varying shock velocity. Interference-filter photographs of the remnant's northeastern, northwestern, and western nebulae in H-alpha + N II, S II, O III, and O I show radial variations of line emission intensity with a scale much larger than that expected from a shock moving through a uniform medium. It is suggested that the emission variations and apparent ionization structure are the result of the shock front encountering an inhomogeneous interstellar medium.

Fesen, R. A.↗

Computing unsteady shock waves for aeroacoustic applications

The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. The paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third-order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation is quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.

Meadows, Kristine R.↗

Computing unsteady shock waves for aeroacoustic applications

The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. This paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation in quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.

Meadows,, Kristine r.↗

Boundary-layer development and transition due to free-stream exothermic reactions in shock-induced flows

A study of the effect of free-stream thermal-energy release from shock-induced exothermic reactions on boundary-layer development and transition is presented. The flow model is that of a boundary layer developing behind a moving shock wave in two-dimensional unsteady flow over a shock-tube wall. Matched sets of combustible hydrogen-oxygen-nitrogen mixtures and inert hydrogen-nitrogen mixtures were used to obtain transition data over a range of transition Reynolds numbers from 1,100,000 to 21,300,000. The heat-energy is shown to significantly stabilize the boundary layer without changing its development character. A method for application of this data to flat-plate steady flows is included.

Hall, J. L.↗

Overview of Recent EAST Testing, Modeling & Analysis

Shock Layer Radiation at NASA Ames - Background: Complex aerothermal and thermochemical phenomena of planetary entry define convective and radiative heating. A spacecraft's TPS (Thermal Protective System) mitigates heat transfer to substructure. Successful TPS design relies on verifiable characterization of these phenomena in the anticipated flight environment. - Approach: EAST (Electric Arc Shock Tube) [facility at Ames] simulates high-enthalpy, real-gas phenomena encountered by hypersonic vehicles entering planetary atmospheres by spectrally imaging the flow behind a moving shock wave. - Goal: Validate aerothermal models (DPLR (Data Parallel Line Relaxation Code) and NEQAIR (Nonequilibrium Radiative Transport and Spectra Program)), inform model improvements, reduce uncertainty and quantify design uncertainties. - Recent Relevant Projects: MSL (Mars Science Laboratory) & Mars 2020, InSight, OSIRIS-REx (Origins Spectral Interpretation Resource Identification Security Regolith Explorer), Orion EFT-1 (Exploration Flight-1) and EM-1 (Exploration Mission-1) and New Frontiers.

Brandis, Aaron↗

Cometary kilometric radio waves and plasma waves correlated with ion pick-up effect at Comet Halley

Bow-shock movements at Comet Halley are inferred from the discrete spectra of the cometary kilometric radiation (30-195 kHz); the observed emissions can be interpreted as being generated and propagating from the moving shock. The shock motion is possibly associated with the time variation of the solar wind and cometary outgassing. It is concluded that these plasma wave phenomena are manifestations of ion pick-up processes, which occur even in a remote region 7 million to 10 million km from the cometary nucleus.

Oya, H.↗

Grid Convergence of High Order Methods for Multiscale Complex Unsteady Viscous Compressible Flows

Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier-Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analogue of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales.

Sjoegreen, B.↗

Fluidic Fence Flow Control on A 30° Swept Wing in Compressible Freestreams

A 30° swept wing with a NACA 64(3)-618 airfoil cross-section was evaluated in four different freestream conditions with Mach numbers increasing from 0.05 to 0.7 using computational fluid dynamics simulations. Active flow control (AFC) through a streamwise row of vortex generating jets was applied at 70% span, blowing inboard to obstruct the development of spanwise flow and increase lift outboard of the control location. The effectiveness of the AFC was evaluated at α = 5° and 10° as freestream Mach number increased and the AFC mass flow coefficient was held constant. Increases to lift coefficient were reduced at higher Mach numbers, but the ratio of change in lift coefficient and the AFC momentum coefficient remained roughly constant. The fluidic fence caused decreased pressure coefficients on the outboard wing sections, but this region of decreased pressure coefficient moved aft on the chord with increasing Mach number. Mach 0.7 is higher than the critical Mach number of the wing and produced interesting interplay with the AFC, as the suction surface shock moved forward on the inboard wing section but moved aft on the outboard wing section.

Evan J McFadden↗

Fluidic Fence Flow Control on A 30° Swept Wing in Compressible Freestreams

A 30° swept wing with a NACA 64(3)-618 airfoil cross-section was evaluated in four different freestream conditions with Mach numbers increasing from 0.05 to 0.7 using computational fluid dynamics simulations. Active flow control (AFC) through a streamwise row of vortex generating jets was applied at 70% span, blowing inboard to obstruct the development of spanwise flow and increase lift outboard of the control location. The effectiveness of the AFC was evaluated at α = 5° and 10° as freestream Mach number increased and the AFC mass flow coefficient was held constant. Increases to lift coefficient were reduced at higher Mach numbers, but the ratio of change in lift coefficient and the AFC momentum coefficient remained roughly constant. The fluidic fence caused decreased pressure coefficients on the outboard wing sections, but this region of decreased pressure coefficient moved aft on the chord with increasing Mach number. Mach 0.7 is higher than the critical Mach number of the wing and produced interesting interplay with the AFC, as the suction surface shock moved forward on the inboard wing section but moved aft on the outboard wing section.

Evan J McFadden↗

Numerical study of unsteady shockwave reflections using an upwind TVD scheme

An unsteady TVD Navier-Stokes solver was developed and applied to the problem of shock reflection on a circular cylinder. The obtained numerical results were compared with the Schlieren photos from an experimental study. These results show that the present computer code has the ability of capturing moving shocks.

Hsu, Andrew T.↗

Model development of supersonic trough wind with shocks

The time dependent one dimensional hydrodynamic equations describe the evolution of the thermal plasma flow along closed magnetic field lines outside of the plasmasphere. The convection of the supersonic polar wind onto a closed fieldline results in the assumed formation of collisionless plasma shocks. These shocks move earthward as the field line with its frozen-in plasma remains fixed or contracts with time to smaller L coordinates. The high equatorial plasma temperature (of the order of electron volts) produced by the shock process decreases with time if the flow is isothermal but it will increase if the contraction is under adiabatic conditions. Assuming adiabaticity a peak in the temperature forms at the equator in conjunction with a depression in the ion density. After an initial contraction, if the flux tube drifts to higher L coordinates the direction of the shock motion can be reversed so that the supersonic region will expand along the field line towards the state characterizing the supersonic polar wind. A rapid expansion will lower the equatorial density while the temperature decreases with time under adiabatic but not isothermal conditions.

Grebowsky, J. M.↗

Model development of supersonic trough wind with shocks.

The time-dependent one dimensional hydrodynamic equations describe the evolution of the thermal plasma flow along closed magnetic field lines outside of the plasmasphere. The convection of the supersonic polar wind onto a closed field line results in the assumed formation of collisionless plasma shocks. These shocks move earthward as the field line with its 'frozen-in' plasma remains fixed or contracts with time to smaller L coordinates. The high equatorial plasma temperature (of the order of electron volts) produced by the shock process decreases with time if the flow is isothermal, but it will increase if the contraction is under adiabatic conditions. Assuming adiabaticity a peak in the temperature forms at the equator in conjunction with a depression in the ion density. After an initial contraction, if the flux tube drifts to higher L coordinates, the direction of the shock motion can be reversed so that the supersonic region will expand along the field line toward the state characterizing the supersonic polar wind.

Grebowsky, J. M.↗

Weak shock compaction on granular salt

This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exhaust Nozzle Plume Effects on Sonic Boom Test Results for Isolated Nozzles

Reducing or eliminating the operational restrictions of supersonic aircraft over populated areas has led to extensive research at NASA. Restrictions were due to the disturbance of the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Recent work has been performed to reduce the magnitude of the sonic boom N-wave generated by airplane components with focus on shock waves caused by the exhaust nozzle plume. Previous Computational Fluid Dynamics (CFD) analysis showed how the shock wave formed at the nozzle lip interacts with the nozzle boat-tail expansion wave. An experiment was conducted in the 1- by 1-ft Supersonic Wind Tunnel at the NASA Glenn Research Center to validate the computational study. Results demonstrated how the nozzle lip shock moved with increasing nozzle pressure ratio (NPR) and reduced the nozzle boat-tail expansion, causing a favorable change in the observed pressure signature. Experimental results were presented for comparison to the CFD results. The strong nozzle lip shock at high values of NPR intersected the nozzle boat-tail expansion and suppressed the expansion wave. Based on these results, it may be feasible to reduce the boat-tail expansion for a future supersonic aircraft with under-expanded nozzle exhaust flow by modifying nozzle pressure or nozzle divergent section geometry.

Castner, Raymond S.↗

Diffraction of shock wave by moving thin wing.

Flow field pressure distribution due to plane shock wave impinging by thin wing moving in opposite direction, discussing mathematical formulation, analytic solution and applications

Gunzburger, M.↗

A General Relativistic Magnetohydrodynamics Simulation of Jet Formation with a State Transition

We have performed the first fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity sim 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving shock wave. This evolution of jet-disk coupling suggests that the low/hard state of the jet system may switch to the high/soft state with a wind, as the accretion rate diminishes.

Nishikawa, K. I.↗

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗

Effect of Scattering on Instabilities in Line-driven Stellar Winds

Line driven O-B stellar winds are unstable to perturbations of short spatial wavelength and the growth rates for such instabilities are very rapid. The nature of the nonlinear development of this instability is unknown, but might possibly be one of blobs of gas driven through ambient gas or a quasiregular train of outward moving shocks. In either case the resulting dissipation of mechanical energy might explain the observed anomalous heating n O-B stars as evidenced by their X-ray emission and high ionization state. It also might explain the observed fine structure of the absorption lines and their time variability. The driving due to the absorption of the stellar continuum flux was considered and the effects of the diffuse, scattered radiation field were neglected. It was shown that under certain special conditions the effect of scattering could reduce the instability growth rate to zero. A stability analysis that includes scattering, but that uses the more physically realistic assumption of complete redistribution instead of coherent scattering, and that includes the effects of transverse velocity gradients, which become important as the flow moves away from the stellar surface is presented. It is found that the instability is eliminated right at the base of the wind, but that as the flow moves outward the instability rate rapidly becomes equal to a substantial fraction of the calculated value, the fraction asymptotically reaching 80% at large radii. Since this still implies many e folds in a characteristic outflow time, the primary conclusion that these winds are highly unstable is unchanged.

Owocki, S. P.↗