Engineering Papers⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Characteristics of coronal shock waves and solar type 2 radio bursts

In the solar corona shock waves generated by flares and/or coronal mass ejections can be observed by radio astronomical methods in terms of solar type 2 radio bursts. In dynamic radio spectra they appear as emission stripes slowly drifting from high to low frequencies. A sample of 25 solar type 2 radio bursts observed in the range of 40 - 170 MHz with a time resolution of 0.1 s by the new radiospectrograph of the Astrophvsikalisches Institut Potsdam in Tremsdorf is statistically investigated concerning their spectral features, i.e, drift rate, instantaneous bandwidth, and fundamental harmonic ratio. In-situ plasma wave measurements at interplanetary shocks provide the assumption that type 2 radio radiation is emitted in the vicinity of the transition region of shock waves. Thus, the instantaneous bandwidth of a solar type 2 radio burst would reflect the density jump across the associated shock wave. Comparing the inspection of the Rankine-Hugoniot relations of shock waves under coronal circumstances with those obtained from the observational study, solar type 2 radio bursts should be regarded to be generated by weak supercritical, quasi-parallel, fast magnetosonic shock waves in the corona.

Mann, G.↗
NASA NTRS1995

Analytic MHD Theory for Earth's Bow Shock at Low Mach Numbers

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven M(A) and sonic M(s) Mach numbers are both much greater than 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first-order perturbation solution is identical to the solution found in the earlier theory. The second-order perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle theta between the normal to the bow shock and magnetic field is not close to 0 deg or 180 deg (the solutions are symmetric about 90 deg). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 deg. In this exceptional case the new solution does not improve on the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta approx. = 90 deg another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) is less than or = 2. MHD and gasdynamic simulations have produced empirical models in which the shock's standoff distance a(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between a(s) and X obtained from MHD simulations, a(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots of X and a(s) predicted for the Earth's bow shock. The plots show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the corresponding phenomenological relations in comparisons with data. Furthermore, significant differences exist between the standoff distances predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.

Grabbe, Crockett L.↗
NASA NTRS1994

Standing shocks in a two-fluid solar wind

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Habbal, Shadia R.↗
NASA NTRS1989

Recurrent solar wind structures in the outer heliosphere

The paper presents recent work on evolution of recurrent solar wind structures in the outer heliosphere. Corotating shocks, corotating interaction regions, and merged interaction regions are studied, and an MHD simulation model in which the jump conditions at all shock crossings satisfy the Rankine-Hugoniot solution is examined. Simulation results which describe the evolution of idealized recurrent solar wind structures between 14 AU and the termination shock are reported which show that merged interaction regions belonging to two to four consecutive solar rotations coalesce with each other between 30 AU and the termination shock to reshape the heliospheric structures. In the outer heliosphere, forward and reverse interplanetary shocks which continuously heat the solar wind plasma are investigated.

Whang, Y. C.↗
NASA NTRS1987

Electron heating in a Monte Carlo model of a high Mach number, supercritical, collisionless shock

Preliminary work in the investigation of electron injection and acceleration at parallel shocks is presented. A simple model of electron heating that is derived from a unified shock model which includes the effects of an electrostatic potential jump is described. The unified shock model provides a kinetic description of the injection and acceleration of ions and a fluid description of electron heating at high Mach number, supercritical, and parallel shocks.

Ellison, Donald C.↗
DOE OSTI2020

Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

In this work, we report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for magneto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six vertices and the center of a hexagon covering approximately one-tenth of the surface area of a 9-ft-diameter spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, and Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging, leading to shock formation rather than potentially more desirable shockless jet merging; and (3) the liner Mach number can remain high (≳10), as required for plasma liners to be an effective compression driver. Secondly, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet mass balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. Ultimately, for both experiments, we present comparisons between experimental and synthetic data from three-dimensional hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Heating of minor ions by the coronal slow shock

The coronal slow shock has been predicted to exist embedded in large coronal holes at 4-10 solar radii. In this paper, a three-fluid model is used to study the jumps in minor ion properties across a slow shock such as the coronal slow shock. The jump conditions are formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. Also obtained is a similarity solution for heavy ions. The results show that, on the downstream side of the slow shock, the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe, in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.

Whang, Y. C.↗
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↗
NASA NTRS1990

Interaction of minor ions with fast and slow shocks

The coronal slow shock was predicted to exist embedded in large coronal holes at 4 to 10 solar radii. A three-fluid model was used to study the jumps in minor ions propertes across the coronal slow shock. The jump conditions were formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. A simularity solution was also obtained for heavy ions. The results show that on the downstream side of the coronal slow shock the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.

Whang, Y. C.↗
NASA NTRS1993

Interaction between perpendicular magnetohydrodynamic shocks

A general analysis is made of the collision and merging of perpendicular shocks as well as the interaction between a shock and a tangential discontinuity. It is found that two head-on shocks diminish both in strength after collisions and a tangential discontinuity forms between them. The property of the discontinuity depends on the relative strength of the two shocks. No discontinuity occurs if the shocks are equal in strength. The emerging of two shocks propagating in the same direction results in a strong shock followed by a tangential discontinuity and a reverse wave. The reverse wave is a rarefaction wave if one or both of the shocks are strong. If the shocks are both weak, a critical adiabatic index (CAI) exists. The reverse wave is a rarefaction wave if the wavelength is less than the CAI and a shock exists if the wavelength is greater than the CAI. As a wake shock enters from a medium of higher wave impedance into that of lower wave impedance, the reflected wave is a rarefaction wave and the total pressure ratio decreases and the velocity jump increases after the shock passes through the border.

Hu, Y. Q.↗
NASA NTRS1983

Multiple spacecraft observations of interplanetary shocks ISEE three-dimensional plasma measurements

ISEE 1 and ISEE 3 three-dimensional solar wind plasma measurements are used together with magnetic field measurements across five previously studied interplanetary shocks to test the accuracy of the mixed-mode shock-normal determination technique and to test whether the shock properties are best approximated with a ratio of specific heats of 5/3 or 2. In the shocks examined, the assumption that the velocity jump was along the normal provided an estimate of the shock normal within 15 deg of the best fit normal 50 percent of the time and within 50 deg, 90 percent of the time. The mixed-mode normals lay within 12 deg of the best fit normal 50 percent of the time and within 36 deg, 90 percent of the time. Part of this deviation may be due to differences in the orientation of the local normal from that of the average normal. Finally, the jump in plasma and field across the shock is better predicted from the Rankine-Hugoniot equations using a ratio of specific heats of 5/3 rather than 2.

Russell, C. T.↗
NASA NTRS1984

The adiabatic energy change of plasma electrons and the frame dependence of the cross-shock potential at collisionless magnetosonic shock waves

The adiabatic energy gain of electrons in the stationary electric and magnetic field structure of collisionless shock waves was examined analytically in reference to conditions of the earth's bow shock. The study was performed to characterize the behavior of electrons interacting with the cross-shock potential. A normal incidence frame (NIF) was adopted in order to calculate the reversible energy change across a time stationary shock, and comparisons were made with predictions made by the de Hoffman-Teller (HT) model (1950). The electron energy gain, about 20-50 eV, is demonstrated to be consistent with a 200-500 eV potential jump in the bow shock quasi-perpendicular geometry. The electrons lose energy working against the solar wind motional electric field. The reversible energy process is close to that modeled by HT, which predicts that the motional electric field vanishes and the electron energy gain from the electric potential is equated to the ion energy loss to the potential.

Goodrich, C. C.↗
NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: The Two Dimensional Time Marching Schemes - 1

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Chang, Sin-Chung↗
NASA NTRS2016

The Efficiency of Magnetic Field Amplification at Shocks by Turbulence

Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks (approx.100 micrograms) and in radio relics in galaxy clusters (approx. micrograms). We present high-resolution magnetohydrodynamic simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be significant. We demonstrate numerically converged field amplification which scales with Alfven Mach number, B/B0 varies as MA, up to MA approx.150.This implies that the post-shock field strength is relatively independent of the seed field. Amplification is dominated by compression at low MA, and stretching (turbulent amplification) at high MA. For high MA, the B-field grows exponentially and saturates at equipartition with turbulence, while the vorticity jumps sharply at the shock and subsequently decays; the resulting field is orientated predominately along the shock normal (an effect only apparent in 3D and not 2D). This agrees with the radial field bias seen in supernova remnants. By contrast, for low MA, field amplification is mostly compressional, relatively modest, and results in a predominantly perpendicular field. The latter is consistent with the polarization seen in radio relics. Our results are relatively robust to the assumed level of gas clumping. Our results imply that the turbulent dynamo may be important for supernovae, but is only consistent with the field strength, and not geometry, for cluster radio relics. For the latter, this implies strong pre-existing B-fields in the ambient cluster outskirts.

lementary particles – magnetic fields – radiat↗
NASA NTRS1982

ISEE-1 and -2 observations of magnetic field strength overshoots in quasi-perpendicular bow shocks

According to the Rankine-Hugoniot jump conditions, the magnetic field strength increases in passing from upstream to downstream of a magnetohydrodynamic shock. At the earth's bow shock the magnetic field often increases significantly more than is required by the Rankine-Hugoniot conditions before decaying in an oscillatory fashion to its downstream average value. Recent OGO-5 and ISEE studies, together with the realization that overshoots also occur in the bow shocks of Venus, Jupiter, and Saturn have focussed attention on these phenomena. It is found that none of the nonlinear dispersive wave theories in the shock literature can account for the overshoot. The present investigation documents the dependence of the overshoot amplitude and thickness on solar wind parameters. Magnetic field overshoots are found to be characteristic of supercritical quasi-perpendicular shocks. The overshoot thickness scales as the ion Larmor radius based upon the solar wind speed and magnetic field.

Livesey, W. A.↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Transonic small disturbance calculations including entropy corrections

Murman's fully conservative mixed type finite-difference operators are first modified. A special sonic point operator with an iterative damping term is introduced which helps the convergence and does not affect the spatial conservative differences. Reliable calculations with second order supersonic schemes are obtained using two sonic operators, the regular sonic point operator followed by a first order supersonic scheme. Also, shock point operator is shown to be equivalent to fitting a locally normal shock terminating the supersonic region. The potential calculations are then modified to account for the non-isentropic jump conditions using a simple shock fitting procedure based on Prandtl relation. The entropy increase across the shock is calculated in terms of the Mach number upstream of the shock and the effect of the generated vorticity is estimated via Crocco relation. Different examples are calculated and extensions to the full potential equation are discussed.

Hafez, M.↗