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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Dynamic-overlapped-grid simulation of aerodynamically determined relative motion

Currently, there is a need to develop a means of analyzing and studying unsteady flowfields which involve multiple component configurations with at least one of the components in relative motion with respect to the others. Two of the important phenomena that such analyses can help to understand are the unsteady aerodynamic interference and the boundary-induced component of the flowfield. With this motivation, a computational method is developed which couples the governing equations of the unsteady flowfield and the rigid-body dynamics in six degrees-of-freedom. These equations are solved on composite meshes of overlapped subdomain grids which can move with respect to each other. Initially, several measures that reduce the numerical error are studied and compared with the exact solution of a moving normal shock in a tube. It is concluded that a second-order accurate method, for spatial and temporal discretizations as well as for the moving subdomain interpolations, is needed as a minimum measure. Furthermore, the CFL numbers should be restricted to unity. Then, the method is used to simulate the flowfield history and predict the aerodynamically determined trajectory of a store dropped from its initial position under a wing.

Yen, Guan-Wei↗
NASA NTRS1995

Waves associated with quasi-parallel shocks: Generation, mode conversion and implications

Ions that are energized at quasi-parallel collisionless shocks and move back upstream generate low-frequency waves, largely on the fast/magnetosonic branch. At sufficient Mach number, the waves are convected back into the shock, lead to shock re-formation, and are mode-converted into downstream (magnetosheath) Alfvenic turbulence. Other waves are generated more locally at the interface of the incoming solar wind and the partially thermalized plasma. This paper reviews how recent simulation studies of collisionless shocks in conjunction with linear kinetic theory and proper wave diagnostics have aided in our understanding of the upstream and magnetosheath waves.

Krauss-Varban, D.↗
DOE OSTI2021

Prediction of Detonation-Induced Disturbances Propagating Upstream into Inlets of Rotating Detonation Combustors

Disturbances caused by the detonation wave in a rotating detonation combustor (RDC) propagate upstream through the inlet, and can potentially affect and couple to upstream components, such as turbomachinery or isolators. These disturbances can potentially also affect the operation of the RDC itself. By drawing from the analogy of a detonation wave bounded by an inert gas, the pressure disturbances observed upstream of the inlet are explained as the consequence of the passage of an upstream propagating oblique shock. In this study, the pressure rise in the plenum from the oblique shock is measured in an axial air inlet RDC. The speed of the upstream propagating wave is estimated to be moderately above the acoustic speed of the oxidizer in the plenum. The wave propagates into the plenum despite local regions of choking in the inlet. It is estimated that the time it takes a fluid particle to transit from the plenum to the detonation channel through the inlet is much larger than the rotational time of the detonation wave. This implies that a fluid particle experiences multiple shocks prior to entering the detonation channel. The oblique shock propagating upstream through the inlet area change is modeled by leveraging an analogy with a quasi-1D shock wave moving in a variable area duct with mean (incoming) flow. Due to flow expansion along the area change, fluid particles are found to experience stronger shocks in the inlet than in the plenum, thereby creating different thermodynamic states within the fill region as the oxidizer emerges from the inlet.

Feleo, Alexander↗
NASA NTRS1986

Relaminarization of the boundary layer over a flat plate in shock tube experiments

The relaminarization of the boundary layer over a flat plate in the shock tube was investigated by using the partially reflected shock wave technique. The flow Mach number was approximately 0.14, which corresponds to the inleft flow Mach number for the first row of vanes in a gas turbine. The thin film platinum heat gauges were used to measure the heat transfer rate and the Stanton number was calculated from the oscilloscope voltage traces. The Reynolds number was varied by changing the operation pressure of the shock tube and the values varied from 2.3 x 10 to the 4th to 5.3 x 10 to the 5th. For a Reynolds number range of 7 x 10 to the 4th to 3.5 x 10 to the 5th, the relaminarization of the boundary layer was observed. This phenomenon is due to the decay of the turbulence level in the flow as the reflected shock wave moves upstream from the flat plate. As the Reynolds number increased, the relaminarization was delayed and the delay was related to the turbulence generated by the reflected shock wave.

Hinckel, J. N.↗
NASA NTRS1992

A burst from a thermonuclear runaway on an ONeMg white dwarf

Studies which examine the consequences of accretion, at rates of 10(exp -9) solar mass/yr and 10(exp -10) solar mass/yr, onto an ONeMg white dwarf with a mass of 1.35 solar masses are performed. In these studies, a Lagrangian, hydrodynamic, one-dimensional computer code was used. The code now includes a network with 89 nuclei up to Ca-40, elemental diffusion, new opacities, and new equation of state. The initial abundance distribution corresponded to a mixture that was enriched to either 25, 50, or 75 percent in products of carbon burning. The remaining material in each case is assumed to have a solar composition. The evolution of the thermonuclear runaway in the 1.35 solar mass white dwarf, with M = 10(exp -9) solar mass, produced peak temperatures in the shell source exceeding 300 million degrees. The sequence produced significant amounts of Na-22 from proton captures onto Ne-20 and significant amounts of Al-26 from proton captures on Mg-24. This sequence ejected 5.2 x 10(exp -6) solar mass moving with speeds from approximately 100 km/s to 2300 km/s. When the mass accretion rate was decreased to 10(exp -10) solar mass, the resulting thermonuclear runaway produced a shock that moved through the outer envelope of the white dwarf and raised the surface luminosity to L greater than 10(exp 7) solar luminosity and the effective temperature to values exceeding 10(exp 7) K. The interaction of the material expanding from off of the white dwarf with the accretion disk should produce a burst of gamma-rays.

Starrfield, S.↗
NASA NTRS2005

Shock Location Dominated Transonic Flight Loads on the Active Aeroelastic Wing

During several Active Aeroelastic Wing research flights, the shadow of the over-wing shock could be observed because of natural lighting conditions. As the plane accelerated, the shock location moved aft, and as the shadow passed the aileron and trailing-edge flap hinge lines, their associated hinge moments were substantially affected. The observation of the dominant effect of shock location on aft control surface hinge moments led to this investigation. This report investigates the effect of over-wing shock location on wing loads through flight-measured data and analytical predictions. Wing-root and wing-fold bending moment and torque and leading- and trailing-edge hinge moments have been measured in flight using calibrated strain gages. These same loads have been predicted using a computational fluid dynamics code called the Euler Navier-Stokes Three Dimensional Aeroelastic Code. The computational fluid dynamics study was based on the elastically deformed shape estimated by a twist model, which in turn was derived from in-flight-measured wing deflections provided by a flight deflection measurement system. During level transonic flight, the shock location dominated the wing trailing-edge control surface hinge moments. The computational fluid dynamics analysis based on the shape provided by the flight deflection measurement system produced very similar results and substantially correlated with the measured loads data.

Lokos, William A.↗
NASA NTRS1988

Evolution of recurrent solar wind structures between 14 AU and the termination shock

The solar wind conditions observed from Voyager 2 at approximately 14 AU are extrapolated to the region of the outer heliosphere bounded by the termination shock, using an MHD simulation model. Results from two simulation studies are presented for two sets of nearly recurrent solar wind interaction regions, with initial conditions generated from plasma and magnetic field data observed on March 1984 at 13.8 AU, and on November 1984 at 15.4, respectively. Each simulation describes an idealized recurrent solar wind structure in the supersonic region of the outer heliosphere out to the termination shock far beyond the present reaches of the Pioneer and Voyager spacecraft. It is shown that a collision between the forward shock and the reverse shock occurs approximately every 40 AU. When a forward shock interacts with the termination shock, the latter is weakened and moves outward; the termination shock is strengthened and moves inward when a reverse shock interacts with it.

Whang, Y. C.↗
NASA NTRS1983

A viscous shock-layer flowfield analysis by an explicit-implicit method

An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.

Gupta, R. M.↗
NASA NTRS1983

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗
NASA NTRS1994

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗
NASA NTRS1994

Disintegration of Dust Aggregates in Interstellar Shocks and the Lifetime of Dust Grains in the ISM

Interstellar grains are destroyed by shock waves moving through the ISM. In fact, the destruction of grains may be so effective that it is difficult to explain the observed abundance of dust in the ISM as a steady state between input of grains from stellar sources and destruction of grains in shocks. This is especially a problem for the larger grains. Therefore, the dust grains must be protected in some way. Jones et al. have already considered coatings and the increased post-shock drag effects for low density grains. In molecular clouds and dense clouds, coagulation of grains is an important process, and the largest interstellar grains may indeed be aggregates of smaller grains rather than homogeneous particles. This may provide a means to protect the larger grains, in that, in moderate velocity grain-grain collisions in a shock the aggregates may disintegrate rather than be vaporized. The released small particles are more resilient to shock destruction (except in fast shocks) and may reform larger grains later, recovering the observed size distribution. We have developed a model for the binding forces in grain aggregates and apply this model to the collisions between an aggregate and fast small grains. We discuss the results in the light of statistical collision probabilities and grain life times.

Dominik, C.↗
NASA NTRS1987

Transonic analysis and design of axisymmetric bodies in nonuniform flow

An inviscid nonuniform axisymmetric transonic code was developed for applications in analysis and design. Propfan slipstream effect on pressure distribution for a body with and without sting was investigated. Results show that nonuniformity causes pressure coefficient to be more negative and shock strength to be stronger and more rearward. Sting attached to a body reduced the pressure peak and moves the rear shock forward. Extent and Mach profile shapes of the nonuniformity region appeared to have little effect on the pressure distribution. Increasing nonuniformity magnitude made pressure coefficient more negative and moved the shock rearward. Design study was conducted with the CONMIN optimizer for an ellipsoid and a body with the NACA-0012 counter. For the ellipsoid, the general trend showed that to reduce the pressure drag, the front portion of the body should be thinner and the contour of the rear portion should be flatter than the ellipsoid. For the design of a body with a sharp trailing edge in transonic flow with an initial shape given by the NACA-0012 contour, the pressure drag was reduced by decreasing the nose radius and increasing the thickness in the aft portion. Drag reduction percentages are given.

Chang, Jen-Fu↗
NASA NTRS2011

Voyager Observations of Magnetic Fields and Cosmic Rays in the Heliosheath

The major features of the profile of greater than 70 MeV/nuc cosmic ray intensity (CRI) observed by Voyager 1 (VI) in the heliosheath from 2005.8-2010.24 are described by the empirical "CR-B" relation as the cumulative effect of variations of the magnetic field strength B. The CRI profile observed by Voyager 2 (V2) from 2008.60 to 2010.28 in the heliosheath is also described by the CR-B relation. On a smaller scale, of the order of a hundred days, a sequence on CRI decreases observed by V 1 during 2006 was interpreted as the effect of a propagating interplanetary shock first interacting with the termination shock, then moving past V1, and finally reflecting from the heliopause and propagating back to V1. Our observations show that the second CRI decrease in this sequence began during the passage of a "Global Merged Interaction Region" (GMIR), ~40 days after the arrival of the GMIR and its possible shock. The first and third CRI decreases in the sequence were associated with local enhancements of B. The magnetic field observations associated with the second sequence of 3 cosmic ray intensity decreases observed by V 1 in 2007/2008 are more difficult to reconcile with the scenario of Webber et al. (2009) and the CR-B relation. The discrepancy might indicate the importance of latitudinal effects

Burlaga, L. F.↗
NASA NTRS2011

Voyager Observations of Magnetic Fields and Cosmic Rays in the Heliosheath

The major features of the profile of >70 MeV/nuc cosmic ray intensity (CRI) observed by Voyager 1 (V1) in the heliosheath from 2005.8 - 2010.24 are described by the empirical "CR-B" relation as the cumulative effect of variations of the magnetic field strength B. The CRI profile observed by Voyager 2 (V2) from 2008.60 to 2010.28 in the heliosheath is also described by the CR-B relation. On a smaller scale, of the order of a hundred days, a sequence of 3 CRI decreases observed by V1 during 2006 was interpreted as the effect of a propagating interplanetary shock first interacting with the termination shock, then moving past V1, and finally reflecting from the heliopause and propagating back to V1. Our observations show that the second CRI decrease in this sequence began during the passage of a "Global Merged Interaction Region" (GMIR), approx. 40 days after the arrival of the GMIR and its possible shock. The first and third CRI decreases in the sequence were associated with local enhancements of B. The magnetic field observations associated with the second sequence of 3 cosmic ray intensity decreases observed by V1 in 2007/2008 are more difficult to reconcile with the scenario of Webber et al. and the CR-B relation. The discrepancy might indicate the importance of latitudinal effects.

Burlaga, L. F.↗
NASA NTRS1995

Two-stream modeling of plasmaspheric refilling

Plasmaspheric refilling on an L = 4 flux tube was studied by using a time-dependent, hydrodynamic plasmaspheric flow model in which the ion streams from the two hemispheres are treated as distinct fluids. In the model the continuity, momentum, and energy equations of a two-ion (O(+) and H(+)), quasi-neutral, currentless plasma are solved along a closed geomagnetic field line; diffusive equilibrium is not assumed. collisions between all stream pairs and with neutral species are included. The model includes a corotating, tilted dipole magnetic field and neutral winds. Ionospheric sources and sinks are accounted for in a self-consistent manner. Electrons are assumed to be heated by photoelectrons. The model flux tube extends from a 200-km altitude in one hemisphere to a 200-km altitude in the other hemisphere. Initially, the upwelling streams pass through each other practically unimpeded. When the streams approach the boundary in the conjugate ionosphere, a shock develops there, which moves upward and dissipates slowly; at about the same time a reverse shock develops in the hemisphere of origin, which moves upward. After about 1 hour, large shocks develop in each stream near the equator; these shocks move toward the equator and downward after crossing the equator. However, these shocks are probably artificial, because counterstreaming flows occur in each H(+) fluid, which the model can only handle by creating shocks.

Guiter, S. M.↗
NASA NTRS1987

Shocked relativistic magnetohydrodynamic flows with application to pulsar winds

The time-dependent behavior of a shocked spherically symmetric relativistic fluid with tangential magnetic field is investigated, considering the case where the boundaries of the shocked fluid move at constant velocity so that self-similar solutions exist. The behavior of the fluid in the ultrarelativistic regime is compared to that in the nonrelativistic regime; there is a smooth transition between these limits. If a magnetic field is present, the magnetic pressure becomes increasingly important with distance from the shock wave; the gas pressure vanishes at the contact discontinuity that bounds the flow. Analytic expressions are given which describe the flow. The solutions can be applied to the evolution of shocked relativistic pulsar winds, which are probably observed as Crab-like supernova remnants. A model for the Crab Nebula, based on the steady-state model of Kennel and Coroniti (1984), indicates that sigma = 0.0016, where sigma is twice the ratio of magnetic to particle energy in the wind as measured in the fluid frame. This is about half the value suggested by Kennel and Coroniti and is much smaller than the value that might be expected for a pulsar wind.

Emmering, Robert T.↗
NASA NTRS1976

Table and charts of equilibrium normal shock and shock tube properties for pure argon with velocities to 18 km/sec

Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves in pure argon. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular-weight ratio, isentropic exponent, velocity, and species mole fractions. Incident (moving) shock velocities are varied from 2 to 18 km/sec for a range of initial pressure of 5 N/sq m to 500 kN/sq m. Working charts illustrating shock tube performance with argon test gas and heated helium and hydrogen driver gases are also presented.

Miller, C. G., III↗
NASA NTRS1974

Tables and charts of equilibrium normal shock and shock tube solutions for pure CO2 with velocities to 16 km/second

Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves in pure CO2, representative of Mars and Venus atmospheres. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular weight ratio, isentropic exponent, velocity and species mole fractions. Incident (moving) shock velocities are varied from 1 to 16 km/sec for a range of initial pressure of 5 Newtons per square meter to 500 kilo Newtons per square meter. The present results are applicable to shock tube flows, and to free-flight conditions for a blunt body at high velocities. Working charts illustrating idealized shock-tube performance with CO2 test gas and heated helium and hydrogen driver gases are also presented.

Miller, C. G., III↗