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 505 records · Page 28

Flux-corrected transport techniques for transient calculations of strongly shocked flows

New flux-corrected transport algorithms are described for solving generalized continuity equations. These techniques were developed by requiring that the finite difference formulas used ensure positivity for an initially positive convected quantity. Thus FCT is particularly valuable for fluid-like problems with strong gradients or shocks. Repeated application of the same subroutine to mass, momentum, and energy conservation equations gives a simple solution of the coupled time-dependent equations of ideal compressible fluid dynamics without introducing an artificial viscosity. FCT algorithms span Eulerian, sliding-rezone, and Lagrangian finite difference grids in several coordinate systems. The latest FCT techniques are fully vectorized for parallel/pipeline processing.

Boris, J. P.↗

Finite-amplitude pressure waves in the radial mode of a cylinder

A numerical study of finite-strength, isentropic pressure waves transverse to the axis of a circular cylinder was made for the radial resonant mode. The waves occur in a gas otherwise at rest, filling the cylinder. A method of characteristics was used for the numerical solution. For small but finite amplitudes, calculations indicate the existence of waves of permanent potential form. For larger amplitudes, a shock is indicated to occur. The critical value of the initial amplitude parameter in the power series is found to be 0.06 to 0.08, under various types of initial conditions.

Kubo, I.↗

A supersonic, three-dimensional code for flow over blunt bodies: User's manual

A computer code is described which may be used to calculate the steady, supersonic, three-dimensional, inviscid flow over blunt bodies. The theoretical and numerical formulation of the problem is given (shock-capturing, downstream marching), including exposition of the boundary and initial conditions. The overall flow logic of the program, its usage, accuracy, and limitations are discussed.

Chaussee, D. S.↗

Type 2 solar radio events observed in the interplanetary medium. Part 1: General characteristics

Twelve type 2 solar radio events were observed in the 2 MHz to 30 kHz frequency range by the radio astronomy experiment on the ISEE-3 satellite over the period from September 1978 to December 1979. These data provide the most comprehensive sample of type 2 radio bursts observed at kilometer wavelengths. Dynamic spectra of a number of events are presented. Where possible, the 12 events were associated with an initiating flare, ground based radio data, the passage of a shock at the spacecraft, and the sudden commencement of a geomagnetic storm. The general characteristics of kilometric type 2 bursts are discussed.

Cane, H. V.↗

Computerized three-dimensional aerodynamic design of a lifting rotor blade

A three-dimensional, inviscid, full-potential lifting rotor code was used to demonstrate that pressure distributions on both advancing and retreating blades could be significantly improved by perturbing local airfoil sections. The perturbations were described by simple geometric shape functions. To illustrate the procedure, an example calculation was made at a forward flight speed of 85 m/sec (165 knots) and an advance ratio of 0.385. It was found that a minimum of three shape functions was required to improve the pressures without producing undesirable secondary effects in high-speed forward flight on a hypothetical modern rotor blade initially having an NLR-1 supercritical airfoil. Reductions in the shock strength on the advancing blade could be achieved, while simultaneously lessening leading-edge pressure gradients on the retreating blade. The major blade section modifications required were blunting of the upper surface leading edge and some reshaping of the blade's upper surface resulting in moderately thicker airfoils.

Tauber, M. E.↗

Influence of numerical dissipation in computing supersonic vortex-dominated flows

Steady supersonic vortex-dominated flows are solved using the unsteady Euler equations for conical and three-dimensional flows around sharp- and round-edged delta wings. The computational method is a finite-volume scheme which uses a four-stage Runge-Kutta time stepping with explicit second- and fourth-order dissipation terms. The grid is generated by a modified Joukowski transformation. The steady flow solution is obtained through time-stepping with initial conditions corresponding to the freestream conditions, and the bow shock is captured as a part of the solution. The scheme is applied to flat-plate and elliptic-section wings with a leading edge sweep of 70 deg at an angle of attack of 10 deg and a freestream Mach number of 2.0. Three grid sizes of 29 x 39, 65 x 65 and 100 x 100 have been used. The results for sharp-edged wings show that they are consistent with all grid sizes and variation of the artificial viscosity coefficients. The results for round-edged wings show that separated and attached flow solutions can be obtained by varying the artificial viscosity coefficients. They also show that the solutions are independent of the way time stepping is done. Local time-stepping and global minimum time-steeping produce same solutions.

Kandil, O. A.↗

A solar-wind 'trigger' for the outer heliosphere radio emissions and the distance to the terminal shock

The solar wind data from the plasma science experiment on the Voyager spacecraft are examined to search for the source of the 2 to 3 kHz radio noise detected in the outer heliosphere. It is found that two anomalous high speed streams passed Voyager 2 before the noise was initially observed. It is suggested that the interaction of these streams with the terminal shock is responsible for the more intense emission. Using a time-of-flight argument, the distance to the shock is estimated at about 70 AU-140 AU. The larger value is consistent with an estimate using Voyager data and a standard pressure-balance argument.

Mcnutt, Ralph L., Jr.↗

Mineralogy and Ar-39 - Ar-40 of an old pristine basalt: Thermal history of the HED parent body

Previous investigations of mineral chemistry and Rb-Sr and Sm-Nd ages indicated that clast,84 from eucrite Yamato 75011 had preserved the pristine nature of its initial crystallization during an early stage of the HED parent body. Microscale mineralogy and Ar-39-Ar-40 ages of this clast, however, revealed local disturbance of microtextures and partially reset ages. This evidence suggests that, in addition to initial crystallization and rapid cooling, the Y75011,84 clast experienced shock deformation, reheating of short duration at higher temperature, and brecciation. These characteristics suggest two or more impact events. Fe-rich olivine filling fractures in pyroxene may have been introduced during the accompanying shock fracturing. The inferred Ar-39-Ar-40 degassing ages for Y75011 matrix and clast, 84 are 3.94 +/- 0.04 Ga and 3.98 +/- 0.03 Ga, respectively. The suggested degassing age for a clast from Y790020, believed to be paired with Y75011, is approximately 4.03 Ga, but could be younger. We consider it likely that all three samples experienced a common degassing event 3.95 +/- 0.05 Ga ago, but we cannot rule out two or more events spaced over a approximately 0.1 Ga interval. Higher temperature extractions of the two clast samples show significantly older apparent ages up to approximately 4.5 Ga and suggest that the time/temperature regime of this event was not sufficient to degas Ar totally. Most likely, the K-Ar ages were reset by thermal metamorphism associated with one or more impact events associated with shock fracturing, formation of Fe-rich olivine veins, and/or meteorite brecciation. The pyroxene annealing that commonly occurs in many eucrites is likely to be a much earlier process than the impact-produced textural changes and reset K-Ar ages observed in these meteorites. The existence of mineralogical and chronological evidence for metamorphism in an otherwise pristine eucrite suggests that the HED parent body experienced an extensive degree of early cratering.

Takeda, Hiroshi↗

Ion Chemistry in Atmospheric and Astrophysical Plasmas

There are many differences and also remarkable similarities between the ion chemistry and physics of planetary ionospheres and the ion chemistry and physics of astronomical environments beyond the solar system. In the early Universe, an expanded cooling gas of hydrogen and helium was embedded in the cosmic background radiation field and ionized by it. As the Universe cooled by adiabatic expansion, recombination occurred and molecular formation was driven by catalytic reactions involving the relict electrons and protons. Similar chemical processes are effective in the ionized zones of gaseous and planetary nebulae and in stellar winds where the ionization is due to radiation from the central stars, in the envelopes of supernovae where the ionization is initiated by the deposition of gamma-rays, in dissociative shocks where the ionization arises from electron impacts in a hot gas and in quasar broad-line region clouds where the quasar is responsible for the ionization. At high altitudes in the atmospheres of the Jovian planets, the main constituents are hydrogen and helium and the ion chemistry and physics is determined by the same processes, the source of the ionization being solar ultraviolet radiation and cosmic rays. After the collapse of the first distinct astronomical entities to emerge from the uniform flow, heavy elements were created by nuclear burning in the cores of the collapsed objects and distributed throughout the Universe by winds and explosions. The chemistry and physics became more complicated. Over 90 distinct molecular species have been identified in interstellar clouds where they are ionized globally by cosmic ray impacts and locally by radiation and shocks associated with star formation and evolution. Complex molecules have also been found in circumstellar shells of evolved stars. At intermediate and low altitudes in the Jovian atmospheres, the ion chemistry is complicated by the increasing abundance of heavy elements such as carbon, and an extensive array of complex molecules has been predicted. Reactions involving heavy elements dominate the structure of the ionspheres of the terrestrial planets and the satellites Titan and Triton.

Dalgarno, A.↗

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Equation of State of CO 2 Shock Compressed to 1 TPa

Equation-of-state (pressure, density, temperature, internal energy) and reflectivity measurements on shock-compressed CO 2 at and above the insulating-to-conducting transition reveal new insight into the chemistry of simple molecular systems in the warm-dense-matter regime. CO 2 samples were precompressed in diamond-anvil cells to tune the initial densities from 1.35 g/cm 3 (liquid) to 1.74 g/cm 3 (solid) at room temperature and were then shock compressed up to 1 TPa and 93 000 K. Lastly, variation in initial density was leveraged to infer thermodynamic derivatives including specific heat and Gruneisen coefficient, exposing a complex bonded and moderately ionized state at the most extreme conditions studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Adiabatic invariants in stellar dynamics, 3: Application to globular cluster evolution

The previous two companion papers demonstrate that slowly varying perturbations may not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with R(sub g) less than or approximately 8 kpc: a concentration c = 1.5 cluster with R(sub g) kpc loses up to 95% of its initial mass in 15 Gyr. Only those with concentration c greater than or approximately 1.3 survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will decrease this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, x, which may be considerably larger than their initial values.

Weinberg, Martin D.↗

Propagation of magnetically guided acoustic shocks in the solar chromosphere

The propagation of a train of acoustic shocks guided by diverging magnetic fields through a static model of the solar chromospheric network and transition region is investigated. For initial flux densities of 1,000,000 ergs/sq cm/sec in the lower chromosphere, the local efficiency of acoustic transmission into the corona can be much higher than that calculated for a plane parallel atmosphere. Acoustic energy will tend to be deposited at higher chromospheric levels in diverging magnetic fields, and magnetic guiding may influence the temperature profile of the network and plages. The total flux that can be transmitted into the corona along such diverging fields is severely limited since the magnetic elements occupy a small fractional area of the photosphere and the transmission efficiency is a rapidly decreasing function of initial flux density. Diverging magnetic fields and a varying ratio of specific heats are not likely to allow high frequency shocks to dissipate high enough in a static atmosphere to contribute significantly to the coronal energy balance.

Foukal, P.↗

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗

Investigation of warm dense matter by combining static and dynamic compression

Shock-compressing material samples that are precompressed with diamond-anvil cells (DACs) enables benchmarking equation of state measurements over a wide range of high-pressure conditions, in particular, warm dense matter relevant to planetary science, astrophysics, and inertial confinement fusion. These targets provide optical access to the shocked sample, so optical velocity interferometry and pyrometry are typical diagnostics. We present a scaling for the coupling of the laser drive to the shock pressure in DACs based on over 150 prior experiments that can be used to design future experiments. For initial static pressures ≤ 5 GPa, using a sapphire window mitigates the effect of photoionization compared to diamond windows and allows the measurement of shock pressures up to 1 TPa with optical diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A shock surface geometry: The 15-16 February 1967 event

The flare-associated interplanetary shock observed by Explorer 33 and Pioneer 7 is analyzed to yield an estimation of the ecliptic plane geometry of the shock surface near 1 AU. These spacecraft were separated by 23 deg in the heliocentric longitude and Pioneer 7 was at a distance of 1.12 AU from the sun. Although a data gap occurred at the apparent time of passage of the disturbance at Pioneer 6, which was 85 deg in heliocentric longitude from Pioneer 7 and at 0.83 AU, the recovered data did suggest such a passage. A consistent picture of the shock propagation is given to explain the difference in arrival times at Pioneers 6, 7, and Explorer 33 and the difference of the shock normals observed by Pioneer 7 and Explorer 33. The average shock speed from the sun to each spacecraft and the local speed at Explorer 33 and their relations to the position of the initiating solar flare are obtained and discussed. In the region of space between the earth and Pioneer 7 the shock surface radius of curvature in the ecliptic plane was 0.4 AU or less.

Lepping, R. P.↗