Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Solar wind dynamic pressure variations and possible ground signatures of flux transfer events

The effects of a series of solar wind dynamic pressure pulses are followed through the solar wind, across the bow shock and magnetosheath, into the magnetosphere, and down to the dayside ionosphere. Each pressure pulse had an intrinsic solar wind origin, i.e., was not generated by processes occurring at the earth's bow shock. Both upstream and downstream of the earth's bow shock, the pulses were identified by depressed magnetic field strengths and enhanced plasma densities. The pulses compressed the magnetosheath and magnetosphere, causing satellites in the magnetosphere to briefly observe enhanced magnetic field strengths, and/or enter the low-latitude boundary layer or magnetosheath. Satellites already in the magnetosheath observed a brief burst of enhanced magnetosheath flow caused by the relative inward motion of magnetosheath flow patterns. Alfven and fast mode compressional waves propagated rapidly down to the polar ionosphere, where they produced transient ionospheric flows at latitudes equatorward of the convection reversal boundary. The flow patterns moved westward around the auroral oval.

Sibeck, D. G.↗
NASA NTRS1993

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gas dynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel. Initial conditions are prescribed by assuming homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. The results suggest possible criteria for assessing the consequences of passage rotation on the wave processes, and they may therefore be applicable to pressure-exchange wave rotors. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface as well as the hub-to-tip radius ratio. Rarefaction fronts, shocks, and contact interfaces are observed to propagate faster with increasing wheel Mach number.

Larosiliere, L. M.↗
NASA NTRS1993

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gas dynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel. Initial conditions are prescribed by assuming homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. The results suggest possible criteria for assessing the consequences of passage rotation on the wave processes, and they may therefore be applicable to pressure-exchange wave rotors. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface as well as the hub-to-tip radius ratio. Rarefaction fronts, shocks, and contact interfaces are observed to propagate faster with increasing wheel Mach number.

Larosiliere, Louis M.↗
NASA NTRS1995

A procedure for automating CFD simulations of an inlet-bleed problem

A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Chyu, Wei J.↗
NASA NTRS1995

Discontinuities formed by steepening the nonlinear Alfven waves

Montgomery (1959) first studied how Alfven waves steepen to form collisionless shocks by use of a two-fluid model. He solved the coupled Lorentz equations of motion for electrons and ions and Maxwell equations by a perturbation technique and showed that linearized polarized Alfven waves in the course of time steepen to form shocks. We have extended this work and have obtained exact simple wave solutions and show that in the limit of slow fluid time scales, we recover the earlier result of Montgomery. For an arbitrary fluid time scale, the solutions include elliptically polarized Alfven waves which also steepen and evolve into shocks. These nonlinear Alfven waves are important in the study of the dynamic solar wind discontinuities.

Lee, N. C.↗
NASA NTRS2009

Shock Structure and Magnetic Fields Generation Associated with Relativistic Jets Unmagnetized Pair Plasma

Using 3D and 2D particle-in-cell simulations we investigate a shock structure, magnetic field generation, and particle acceleration associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized pair plasma. The simulations use long computational grids which allow to study the formation and dynamics of the system in a spatial and temporal way. We find for the first time a relativistic shock system comparable to a predicted magnetohydrodynamic shock structure consisting of leading and trailing shocks separated by a contact discontinuity. Strong electromagnetic fields resulting from the Weibel two-stream instability are generated in the trailing shock where jet matter is thermalized and decelerated. We analyze the formation and nonlinear development through saturation and dissipation of those fields and associated particle acceleration. In the AGN context the trailing shock corresponds to the jet shock at the head of a relativistic astrophysical jet. In the GRB context this trailing shock can be identified with the bow shock driven by relativistic ejecta. The strong electromagnetic field region in the trailing shock provides the emission site for the hot spot at the leading edge of AGN jets and for afterglow emission from GRBs.

Niemiec, J.↗
NASA NTRS2000

Applications of Low Density Flow Techniques and Catalytic Recombination at the Johnson Space Center

The talk presents a brief background on defInitions of catalysis and effects associated with chemically nonequilibrium and low-density flows of aerospace interest. Applications of catalytic recombination on surfaces in dissociated flow are given, including aero heating on reentry spacecraft thermal protection surfaces and reflection of plume flow on pressure distributions associated with the space station. Examples include aero heating predictions for the X-38 test vehicle, the inlet of a proposed gas-sampling probe used in high enthalpy test facilities, and a parabolic body at angle of attack. The effect of accommodation coefficients on thruster induced pressure distributions is also included. Examples of tools used include simple aero heating formulas based on boundary layer solutions, an engineering approximation that uses axisymmetric viscous shock layer flow to simulate full three dimensional flow, full computational fluid dynamics, and direct simulation Monte-Carlo calculations. Methods of determining catalytic recombination rates in arc jet flow are discus ed. An area of catalysis not fully understood is the formation of single-wall carbon nanotubes (SWNT) with gas phase or nano-size metal particles. The Johnson Space Center is making SWNTs using both a laser ablation technique and an electric arc vaporization technique.

Scott, Carl D.↗
NASA NTRS2017

The Mars Crustal Magnetic Field Control of Plasma Boundary Locations and Atmospheric Loss: MHD Prediction and Comparison with MAVEN

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

Global↗
NASA NTRSDate not supplied

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

Aerothermodynamic prediction is important for the design and analysis of many aerospace vehicles. Computational Fluid Dynamics (CFD) tools used for these predictions commonly rely on structured shock-aligned grids due to the strong shocks exhibited at high speeds. This work details the incorporation of an improved HLLE++ scheme for both perfect gas and thermochemical nonequilibrium flows which robustly captures shock waves on both non-shock-aligned and shock-aligned grids. The HLLE++ scheme reduces the impact of eigenvalue limiting on flow solutions which is required for Roe's scheme for high Mach number flows. An adapted grid approach is also demonstrated using mixed-element unstructured grids to improve heat prediction. Thin prismatic boundary layers are used with tetrahedra in the farfield and for capturing shock waves. Various results are presented for hypersonic flows over blunt bodies including cylinders, spheres, and capsules.

CFD↗
NASA NTRSDate not supplied

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

Aerothermodynamic prediction is important for the design and analysis of many aerospace vehicles. Computational Fluid Dynamics (CFD) tools used for these predictions commonly rely on structured shock-aligned grids due to the strong shocks exhibited at high speeds. This work details the incorporation of an improved HLLE++ scheme for both perfect gas and thermochemical nonequilibrium flows which robustly captures shock waves on both non-shock-aligned and shock-aligned grids. The HLLE++ scheme reduces the impact of eigenvalue limiting on flow solutions which is required for Roe's scheme for high Mach number flows. An adapted grid approach is also demonstrated using mixed-element unstructured grids to improve heat prediction. Thin prismatic boundary layers are used with tetrahedra in the farfield and for capturing shock waves. Various results are presented for hypersonic flows over blunt bodies including cylinders, spheres, and capsules.

CFD↗
DOE OSTI2024

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗
NASA NTRS1985

The central engine of quasars and AGNs: A relativistic proton radiative shock

Active galactic nuclei (AGNs) and quasars (QSOs) appear to emit roughly equal energy per decade from radio to gamma-ray energies (e.g. Ramaty and Ligenfelter 1982). This argues strongly for a nonthermal radiation mechanism (see Rees 1984). In addition, statistical studies have indicated that the spectra of these objects in the IR-UV and 2 to 50 keV X-ray band, can be fitted very well with power laws of specific indices. These spectral indices do not seem to depend on the luminosity or morphology of the objects (Rothschild et al. 1983; Malkan 1984), and any theory should account for them in a basic and model independent way. If shocks accelerate relativistic protons via the first-order Fermi mechanism (e.g. Axfor 1981), the radiating electrons can be produced as secondaries throughout the source by proton-proton (p-p) collisions and pion decay, thus eliminating Compton losses (Protheroe and Kazanas 1983). As shown by Kazanas (1984), if relativistic electrons are injected at high energies, e+-e- pair production results in a steady state electron distribution that is very similar to that observed in AGNs, independent of the details of injection and the dynamics of the source. The conditions required by this mechanism are met in the shock model of Eichler (1984) and Ellison and Eichler (1984) which allows the self-consistent calculation of the shock acceleration efficiency.

Kazanas, D.↗
NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗
DOE OSTI2024

Quantifying motion blur by imaging shock front propagation with broadband and narrowband X-ray sources

Time-integrated radiography using MeV Bremsstrahlung X-ray sources is the norm for imaging during system-level testing of components and structures under dynamic condition. One source of error in the analysis of the time-integrated radiography data sets stems from motion blur which smears out sharp interfaces to a greater degree with longer exposure times, which become necessary to provide sufficient signal-to-noise with low X-ray penetration of objects of interest. To quantify motion blur, a 1D shock wave through PMMA was investigated experimentally at The Dynamic Compression Sector at The Advanced Photon Source (DCS@APS) with tapered broadband and 25.46 ± 1.06 keV narrowband X-rays. Four cameras with different exposure times were used for each experiment to compare the effect that exposure time has on motion blur. In addition, our methodology to accurately simulate motion blur in terms of transmission and shape is presented and compared to our experimental results and quantified. There is a high level of agreement between the experimental and simulation results across the range of data sets investigated in this study with a percent difference range of 0.29–1.31% for the four shots. The methodology of this work serves as a steppingstone towards a physically validated model that could be used in conjunction with experimental results to deconvolve physical parameters, densities, and interfaces of interest in a way that would not be possible with experimental results alone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1995

Electron acceleration to high energies at quasi-parallel shock waves in the solar corona

In the solar corona shock waves are generated by flares and/or coronal mass ejections. They manifest themselves in solar type 2 radio bursts appearing as emission stripes with a slow drift from high to low frequencies in dynamic radio spectra. Their nonthermal radio emission indicates that electrons are accelerated to suprathermal and/or relativistic velocities at these shocks. As well known by extraterrestrial in-situ measurements supercritical, quasi-parallel, collisionless shocks are accompanied by so-called SLAMS (short large amplitude magnetic field structures). These SLAMS can act as strong magnetic mirrors, at which charged particles can be reflected and accelerated. Thus, thermal electrons gain energy due to multiple reflections between two SLAMS and reach suprathermal and relativistic velocities. This mechanism of accelerating electrons is discussed for circumstances in the solar corona and may be responsible for the so-called 'herringbones' observed in solar type 2 radio bursts.

Mann, G.↗
NASA NTRS1980

Effects of atmospheric breakup on crater field formation

This paper investigates the physics of meteoroid breakup in the atmosphere and its implications for the observed features of strewn fields. There are several effects which cause dispersion of the meteoroid fragments: gravity, differential lift of the fragments, bow shock interaction just after breakup, centripetal separation by a rotating meteroid, and possibly a dynamical transverse separation resulting from the crushing deceleration in the atmosphere. Of these, it is shown that gravity alone can produce the common pattern in which the largest crater occurs at the downrange end of the scatter ellipse. The average lift-to-drag ratio of the tumbling fragments must be less than about 0.001, otherwise small fragments would produce small craters downrange of the main crater, and this is not generally observed. The cross-range dispersion is probably due to the combined effects of bow shock interaction, crushing deceleration, and possibly spinning of the meteoroid. A number of terrestrial strewn fields are discussed in the light of these ideas, which are formulated quantitatively for a range of meteoroid velocities, entry angles, and crushing strengths. It is found that when the crater size exceeds about 1 km, the separation between the fragments upon landing is a fraction of their own diameter, so that the crater formed by such a fragmented meteoroid is almost indistinguishable from that formed by a solid body of the same total mass and velocity.

Passey, Q. R.↗
NASA NTRS1983

Cloud-particle galactic gas dynamics and star formation

Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.

Roberts, W. W., Jr.↗
DOE OSTI2020

Experimental determination of shock speed versus exciter speed in a two-dimensional dusty plasma

A shock that is continuously driven by a moving exciter will propagate at a speed that depends on the exciter speed. We obtained this dependence experimentally, in a strongly coupled dusty plasma that was prepared as a single two-dimensional layer of charged microparticles. Attaining this result required an experimental advance, developing a method of driving a shock continuously, which we did using an exciter moving at a constant supersonic speed, analogous to a piston in a cylinder. The resulting compressional pulse was a shock that propagated steadily without weakening, ahead of the moving exciter. In this work, we compare our experimental results to an empirical form $M_{\text{shock}} = 1 + sM_{\text{exciter}}$, and to the prediction of a recent simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗