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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Solar cycle variations in the neutral exosphere inferred from the location of the Venus bow shock

Solar UV and EUV varies significantly during the solar cycle. Pioneer Venus can measure this variation both directly and indirectly. A direct measure of the EUV is obtained from the photoelectron current of the Langmuir probe when the spacecraft is in the solar wind. The indirect measure is by monitoring the location of the Venus bow shock. The UV and EUV both heat the upper atmosphere and ionize it. When solar activity is high, the upper atmosphere should be ionized more rapidly. This effect adds a greater number of planetary ions to the magnetosheath plasma as it flows by Venus. It is this increase in mass flow that causes the Venus bow shock to move away from its solar minimum location. Pioneer Venus has now monitored the location of the bow shock for an entire solar cycle. The bow shock location is well correlated with the variation in EUV flux as measured by the Langmuir probe. The bow shock is farther from Venus than expected from the sunspot number or 10.7 cm solar radio flux, indicating that solar UV radiation may be even stronger at the present time than would be predicted from the relationships determined during the previous solar cycle.

Russell, C. T.↗

NASA NTRS1995

Global aspects of the motion of the heliospheric termination shock: A gasdynamic model

The heliospheric termination shock is expected to move in response to variation in upstream solar wind conditions. Using numerical techniques, we extend an earlier strictly one-dimensional analytic gas dynamic model of shock motion to two dimensions, to investigate the qualitative features of global behavior of shock motion, and the consequences of latitudinal variation in dynamic pressure. The boundary conditions of the calculation are given by the solar wind parameters as a function of latitude and time on an inner spherical boundary, and a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. Density variations, specified at the inner boundary as a function of time, are convected into the termination shock. Immediately after the interaction, the shock moves with speeds given by the earlier analytic model. However, as the termination shock propagates outward (or inward), it begins to slow down. After about 2 to 10 years, depending on details of boundary conditions, the signal from the shock interaction has reached the outer boundary and propagates inward to the position of the termination shock, strongly affecting the behavior of the shock. Assuming no further disturbances in the solar wind, the termination shock will reach its new equilibrium after some tens of years. In reality, large-scale variations in solar wind dynamic pressure occur on time scales short in comparison with the eleven year solar cycle, so that one expects that the termination shock is never in an equilibrium position, but rather oscillates inward and outward; this oscillation will vary with heliographic latitude. The effects of a variety of types of solar wind disturbances are investigated and summarized.

Ratkiewicz, R.↗

NASA NTRS1983

Collisionless electron shocks in electron-beam-plasma systems

One-dimensional Vlasov simulations show that when an electron beam is suddenly injected into a plasma, a fast-moving monotonic shock forms during the early stage of the transient plasma response. In the shock, the ions are nearly immobile. The shock appears to be an electrostatic electron-beam-plasma mode. The shock evolves from an initial positive potential perturbation, which is supported by an ion burst. The steepening of the perturbation into a shock is characterized by an electron-beam-plasma mode.

Singh, N.↗

NASA NTRS1990

Simulated Hypersonic Flows About A Blunt Body

Unsteady and steady flows compared. Report describes computer numerical study of two-dimensional, unsteady, viscous, hypersonic flows of air about blunt body with impinging shock. This kind of flow represents many practical phenomena; for example, interaction of fluctuating bow shock of hypersonic airplane with shocks of leading edge of wing or of lip of cowl at inlet to engine. Such interactions give rise to complicated, moving shock-on-shock patterns.

Kutler, P.↗

NASA NTRS2015

Plume and Shock Interaction Effects on Sonic Boom in the 1-foot by 1-foot Supersonic Wind Tunnel

The desire to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas has led to extensive research at NASA. Restrictions are due to the disturbance of the sonic boom, caused by the coalescence of shock waves formed by the aircraft. A study has been performed focused on reducing the magnitude of the sonic boom N-wave generated by airplane components with a focus on shock waves caused by the exhaust nozzle plume. Testing was completed in the 1-foot by 1-foot supersonic wind tunnel to study the effects of an exhaust nozzle plume and shock wave interaction. The plume and shock interaction study was developed to collect data for computational fluid dynamics (CFD) validation of a nozzle plume passing through the shock generated from the wing or tail of a supersonic vehicle. The wing or tail was simulated with a wedgeshaped shock generator. This test entry was the first of two phases to collect schlieren images and off-body static pressure profiles. Three wedge configurations were tested consisting of strut-mounted wedges of 2.5- degrees and 5-degrees. Three propulsion configurations were tested simulating the propulsion pod and aft deck from a low boom vehicle concept, which also provided a trailing edge shock and plume interaction. Findings include how the interaction of the jet plume caused a thickening of the shock generated by the wedge (or aft deck) and demonstrate how the shock location moved with increasing nozzle pressure ratio.

Shock wave interaction↗

NASA NTRS1985

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities.

Davis, S. F.↗

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.↗

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.↗