Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hugoniot”

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 325 records · Page 18

Shock experiments on maskelynite-bearing anorthosite

A series of shock recovery experiments over 9.9-60.4 GPa have been carried out on naturally shocked anorthosite from the Mistastin impact structure in Labrador consisting primarily of diaplectic plagioclase glass or maskelynite, An(50), and pyroxene. Petrographic observations of the experimental products indicate that the component minerals and diaplectic glasses generally retained their initial character throughout, the only exception being the increase in fracturing which occurred in the 9.9 GPa shot. Reshocking at pressures higher than the initial shock tends to lower the refractive index of maskelynite. The increase in refractive index of maskelynite reshocked to pressures lower than the initial pressure is interpreted as due to shock densification of the diaplectic glass above the Hugoniot elastic limit and below the mixed phase regime. The data suggest that the low-high-low density transition of maskelynite occurs about 8 GPa below that of the crystal of corresponding composition.

Lambert, P.↗

Structure of the November 12, 1978, quasi-parallel interplanetary shock

The jump in plasma parameters exhibited by the intense interplanetary shock event of Nov. 12, 1978 is analyzed using ISEE 1, 2 and 3 data. Magnetic and electric field measurements indicated that the shock magnetic field profile was similar to the earth bow shock profile. Data on the electron and proton densities, temperatures, bulk velocities and alpha particles showed a steady electron temperature increase across the shock on a 12 earth radii scale. The upstream and downstream flow parameters are found to be within 10 percent of Rankin-Hugoniot jump conditions. The shock moved at 614 km/sec and had three dissipative scales, one a few Larmor radii determined by the magnetic field jump, a second 10 earth radii correlated with the electron equilibrium and the other 30 earth radii connected to the energetic proton foreshock.

Kennel, C. F.↗

Particle velocity experiments in anorthosite and gabbro

Shock wave experiments were conducted in San Gabriel anorthosite and San Marcos gabbra 10 11 GPa using a 40 mm-borne propellant gun. Particle velocities were measured directly at several points in each target by means of electromagnetic gauges. Hugoniot states were calculated by determining shock-transit time from the gauge records. Sound speeds indicate a loss of shear strength upon sock compression for both rocks, with the strength loss persisting upon release to zero stress om the anorthosite. Stress-density release paths in the anorthosite indicate possible transformation of albite to jadeite + (quartz or coesite), with the amount of material transformed increasing as a function of shock stress. Electrical interferene effects in the gabbro precluded the determination of accurate release paths for the rock.

Boslough, M. B.↗

Numerical studies of 2-dimensional flows

A formulation of the lambda scheme for the analysis of two dimensional inviscid, compressible, unsteady transonic flows is presented. The scheme uses generalized Riemann variables to determine the appropriate two point, one sided finite difference approximation for each derivative in the unsteady Euler equations. These finite differences are applied at the predictor and corrector levels with shock updating at each level. The weaker oblique shocks are captured, but strong near normal shocks are fitted into the flow using the Rankine-Hugoniot relations. This code is demonstrated with a numerical example of a duct flow problem with developing normal and oblique shock waves. The technique is implemented in a code which has been made efficient by streamlining to a minimal number of operations and by eliminating branch statements. The scheme is shown to provide an accurate analysis of the flow, including formation, motions, and interactions of shocks; the results obtained on a relatively coarse mesh are comparable to those obtained by other methods on much finer meshes.

Moretti, G.↗

The numerical simulation of steady transonic rotational flow using a dual potential formulation

A finite-difference method is presented that simulates steady transonic rotational flow of an inviscid fluid by representing the velocity field as the sum of scalar and vector potentials. This dual potential velocity decomposition extends the validity of the scalar (full) velocity potential to include vorticity. The inclusion of a vector potential also permits an alternate treatment of lift that does not require a circulation wake cut. This is accomplished by specifying the vector potential as a constant on the airfoil surface in order to satisfy a Kutta condition. The governing equations are solved as iteratively decoupled scalar equations using approximate factorization techniques, and the overall efficiency approaches that of the full potential equation. The governing equations are able to convect entropy and vorticity throughout the flow field and are equivalent to the Euler equations in continuous flow domains, however at shocks the Rankine-Hugoniot entropy jump must be supplied. An entropy correction method is presented and verified with transonic airfoil solutions of the Euler equations.

Chaderjian, N. M.↗

Heating of electrons in turbulence around the Space Shuttle

A model is proposed in which the strong hydrodynamic shock forming around the Space Shuttle, and occasionally detaching, generates strong hydrodynamic turbulence downstream. The level of this turbulence is found to be sufficiently great to account for the heating of electrons that has been observed. A shock formed in front of the Orbiter is shown to satisfy the Rankine-Hugoniot momentum and pressure balance equations. The hydrodynamic turbulence is coupled through ion-neutral collisions to the O2(+) and N2(+) ions of outgas.

Parish, J. L.↗

Application of shock wave data to earth and planetary science

It is pointed out that shock wave data for: (1) low temperature condensable gases H2 and He, (2) H2O, CH4, NH3, CO, CO2, and N2 ices, and (3) silicates, metals, oxides and sulfides have many applications in geophysics and planetary science. The present paper is concerned with such applications. The composition of planetary interiors is discussed, taking into account the division of the major constituent of the planets in three groups on the basis of 'cosmic abundance' arguments, the H-He mixtures in the case of Jupiter and Saturn, shock wave data for hydrogen, and constraints on the internal structure of Uranus and Neptune. Attention is also given to the earth's mantle, shock wave data for mantle materials, the earth's core, impacts on planetary surfaces, elastic wave velocities as a function of pressure along the Hugoniot of iron, and reactions which yield the CO2 bearing atmospheres for Venus, earth, and Mars.

Ahrens, T. J.↗

The effect of heavy ions on the formation and structure of cometary bow shocks

A hybrid simulation model is used to investigate the effects of heavy cometary ions on the formation and structure of a cometary bow shock. The calculations are carried out over various Mach numbers and heavy ion velocity distribution functions. The model is based on previous formulations for phenomena in the solar wind and at the earth's bow shock. The generation of the shock is described in terms of particles injected from one side of the simulation field and reflected from the other end of the field, i.e., a solid wall boundary model. This technique permits the steep buildup of the ion density near the cometary nucleus, followed by coupling of the incident and reflected ion streams to produce a shock. It is shown that at low Mach numbers (up to Mach 2) the shock is transitory and periodically formed by protons, then destroyed by heavy ions (O+). Slightly higher Mach numbers lead to a true stationary shock. An examination of coupling effects between the solar wind and the heavy ions at low Mach numbers by using the Rankine-Hugoniot relations reveals that the ions and the solar wind protons cannot be treated as a single fluid calculating the shock characteristics.

Omidi, N.↗

Magnetohydrodynamic model of the Crab Nebula and its radiation

The interaction of the Crab pulsar with its surrounding nebula is modeled as a steady state, spherically symmetric relativistic wind which is confined by the outer layers of the supernova remnant. The pulsar wind is assumed to contain only positrons and electrons, and is terminated by a strong shock at a distance of 10 arcseconds from the pulsar, at the outer edge of the underluminous region. Down-stream of the shock, the flow decelerates and increases its pressure in order to match the boundary conditions imposed by outer nebula. Consistent flow solutions can be found if the magnetic luminosity of the pulsar wind is about 0.3 percent of its particle luminosity. The immediate postshock plasma is modeled as a power law distribution whose density and pressure moments satisfy the strong shock Rankine-Hugoniot relations. Using the downstream flow solution, the synchrotron continuum is calculated. The parameters of the pulsar wind and assumed power law spectral index are adjusted until the calculated total synchrotron power and the spectral distribution of the nebular continuum from optical to gamma-rays, are in reasonable accord with observations.

Coroniti, F. V.↗

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner has been analyzed using the complete set of ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream have been analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity has then been compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The various measurements have also been compared with values derived from a Rankine-Hugoniot analysis. The results indicate that, inbound, the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained, however, arguments are presented that the bow wave is also likely to be a shock at this location.

Smith, E. J.↗

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner was analyzed using the ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream were analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity was compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The measurements were compared with values derived from a Rankine-Hugoniot analysis. The results indicate that inbound the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained. However, arguments that the bow wave is also likely to be a shock at this location are presented.

Smith, E. J.↗

Unsteady hybrid vortex technique for transonic vortex flows and flutter applications

The well developed surface panel method is extended by adding a volume integral term allowing calculation of the full effect of compressibility. The full effect of compressibility is calculated by using mixed type finite difference scheme consistent with the mixed nature of transonic flow. The solution is obtained through successive iteration cycles for subcritical flows and for critical flows the solution is obtained through satisfying the Rankine-Hugoniot relations across the captured shock in addition to the captured shock. The method is applied to airfoils in subcritical and critical flows and the results are in good agreement with the experimental data and finite-difference solutions of the full-potential and Euler equations.

Kandil, Osama A.↗

Transonic airfoil computation using the integral equation with and without embedded Euler domains

Two transonic computational schemes which are based on the Integral Equation Formulation of the full potential equation were presented. The first scheme is a Shock Capturing-Shock Fitting (SCSF) scheme which uses the full potential equation throughout with the exception of the shock wave where the Rankine-Hugoniot relations are used to cross and fit the shock. The second scheme is an Integral Equation with Embedded Euler (IEEE) scheme which uses the full potential equation with an embedded region where the Euler equations are used. The two schemes are applied to several transonic airfoil flows and the results were compared with numerous computational results and experimental domains with fine grids. The SCSF-scheme is restricted to flows with weak shock, while the IEEE-scheme can handle strong shocks. Currently, the IEEE scheme is applied to other transonic flows with strong shocks as well as to unsteady pitching oscillations.

Kandil, Osama A.↗

Molecular dissociation and shock-induced cooling in fluid nitrogen at high densities and temperatures

Radiative temperatures and electrical conductivities were measured for fluid nitrogen compressed dynamically to pressures of 18-90 GPa, temperatures of 4000-14,000 K, and densities of 2-3 g/cu cm. The data show a continuous phase transition above 30 GPa shock pressure and confirm that (delta-P/delta-T)v is less than 0, as indicated previously by Hugoniot equation-of-state experiments. The first observation of shock-induced cooling is also reported. The data are interpreted in terms of molecular dissociation, and the concentration of dissociated molecules is calculated as a function of density and temperature.

Radousky, H. B.↗

A quarter century of collisionless shock research

This review highlights conceptual issues that have both governed and reflected the direction of collisionless shock research in the past quarter century. These include MHD waves and their steepening, the MHD Rankine-Hugoniot relations, the super-critical shock transition, nonlinear oscillatory wave trains, ion sound anomalous resistivity and the resistive-dispersive transition for subcritical shocks, ion reflection and the structure of supercritical quasi-perpendicular shocks, the earth's foreshock, quasi-parallel shocks, and finally, shock acceleration processes.

Kennel, C. F.↗

Ion reflection, gyration, and dissipation at supercritical shocks

Topics on collisionless plasmas are discussed and include: (1) Ion reflection is the dominant ion dissipation mechanism at nearly perpendicular, supercritical shocks. (2) An increasing fraction of the ions incident on a supercritical shock is reflected as the Mach number increases. The actual fraction reflected can be predicted using the Rankine Hugoniot conservation relations. (3) The effective temperature associated with the dispersion in velocity space associated with ion reflection accounts for a large fraction of the temperature rise observed across supercritical, quasi perpendicular shocks.

Gosling, J. T.↗

A comparative study of plasma heating by ion acoustic and modified two-stream instabilities at subcritical quasi-perpendicular shocks

Plasma heating due to the ion acoustic instability and the modified two-stream instability is examined for quasi-perpendicular subcritical shocks. Electron and ion heating is investigated as a function of upstream electron to ion temperature ratio and plasma beta using second-order heating rates. A simple shock model is employed in which the cross-field electron-ion drift speed is adjusted until the total (adiabatic plus anomalous) heating matches that required by the Rankine-Hugoniot relations. Quantities such as the width of the shock and the maximum electric field fluctuations are also calculated, and the results are compared with the ISEE data set of subcritical bow shock crossings. The observed width of the shock, the amount of plasma heating, and the low-frequency electric field intensity are in reasonably good agreement with the calculations for the modified two-stream instability. On the other hand, the wave intensities at higher frequency are about 4 orders of magnitude smaller than those predicted for the ion acoustic instability at saturation, consistent with the fact that the measured shock widths imply cross-field drift speeds that are below threshold for this instability. It is therefore concluded that the dissipation at these shocks is most likely due to the lowest frequency, modified two-stream instability.

Winske, D.↗