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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
NASA NTRS1975

Anisotropy measurements of nearly 50 KeV solar protons

The Energetic Particles Experiment on IMP-7 measures the angular distribution of 50-200 keV solar protons in 16 sectors. The velocity of 50 keV protons may be less than 5 times that of the solar wind. A generalized nonlinear Compton-Getting point transformation into the co-moving frame that contains no assumptions as to the angular distribution of either the spectrum or intensity is presented. Nearly 50 keV proton data in the spacecraft frame exhibit an anisotropy ratio that is large (not less than 5) and radial throughout the October 29, 1972 event lasting more than 9 days at this energy. This anisotropy argues against impulsive injection and diffusive decay in the inner solar system. Application of the transformation to the data reveals a long lasting residual anisotropy in the co-moving frame with protons streaming from the sun. Differences between the co-moving frame and solar wind frame velocities suggest residual electric fields upstream from the bow shock.

Gold, R. E.↗
NASA NTRS1991

Application of an unsteady Navier-Stokes solver to transonic turbine design

This study presents a numerical evaluation of the performance of the first stage of a new-generation turbine design. The numerical method solves the two-dimensional Navier-Stokes equations using a system of patched grids. Three-dimensional effects of stream-tube contraction are also modeled. The study focuses on the effects of axial gap variation on the unsteady rotor-stator interactions and on stage performance. Results are presented for three different axial gaps. The results indicate that the unsteady interactions can be very large in this design. These interactions affect not only the stage efficiency but also substantially alter the time-averaged features of the flow. In particular, for the case of the smallest axial gap, it was found that there was an unsteady shock on the stator suction surface which spanned the gap region and impinged upon the moving rotor airfoils.

Rangwalla, Akil A.↗
NASA NTRS1994

A model of sonoluminescence

A bubble of air, trapped at the center of a spherical container of water on the surface of which spherical sound waves are maintained by transducers, may emit light, a phenomenon known as sonoluminescence. The surface of the bubble expands and contracts in obedience to the Rayleigh-Lamb equation, which requires knowledge of the gas pressure on the surface of the bubble. In many investigations of bubble pulsations, it is assumed that the air in the bubble moves adiabatically. To understand sonoluminescence, however, it is necessary to allow for the possibility that shocks are generated within the bubble. We couple the Rayleigh-Lamb equation governing the bubble radius to Euler's equations governing the motion of air in the bubble, and solve the two equations simultaneously. The air is modelled by a van der Waals gas. Results are presented for a number of slightly different conditions of excitation, but in which the response of the system is widely different.

Wu, C. C.↗
NASA NTRS2012

Hybrid Simulation of the Shock Wave Trailing the Moon

A standing shock wave behind the Moon was predicted by Michel (1967) but never observed nor simulated. We use 1D hybrid code in order to simulate the collapse of the plasma-free cavity behind the Moon and for the first time to model the formation of this shock. Starting immediately downstream of the obstacle we consider the evolution of plasma expansion into the cavity in the frame of reference moving along with the solar wind. Well-known effects as electric charging of the cavity affecting the plasma flow and counterstreaming ion beams in the wake are reproduced. Near the apex of the inner Mach cone where the plasma flows from the opposite sides of the obstacle meet, a shock wave arises. We expect the shock to be produced at periods of high electron temperature solar wind streams (T(sub i) much less than T(sub e) approximately 100 eV). The shock is produced by the interaction of oppositely directed proton beams in the plane containing solar wind velocity and interplanetary magnetic field vectors. In the direction across the magnetic field and the solar wind velocity, the shock results from the interaction of the plasma flow with the region of the enhanced magnetic field inside the cavity that plays the role of the magnetic barrier. The appearance of the standing shock wave is expected at the distance of approximately 7R(sub M) downstream of the Moon.

Moon↗
NASA NTRS2012

Solar Wind Charge Exchange During Geomagnetic Storms

On March 31st. 2001, a coronal mass ejection pushed the subsolar magnetopause to the vicinity of geosynchronous orbit at 6.6 RE. The NASA/GSFC Community Coordinated Modeling Center (CCMe) employed a global magnetohydrodynamic (MHD) model to simulate the solar wind-magnetosphere interaction during the peak of this geomagnetic storm. Robertson et aL then modeled the expected 50ft X-ray emission due to solar wind charge exchange with geocoronal neutrals in the dayside cusp and magnetosheath. The locations of the bow shock, magnetopause and cusps were clearly evident in their simulations. Another geomagnetic storm took place on July 14, 2000 (Bastille Day). We again modeled X-ray emission due to solar wind charge exchange, but this time as observed from a moving spacecraft. This paper discusses the impact of spacecraft location on observed X-ray emission and the degree to which the locations of the bow shock and magnetopause can be detected in images.

Robertson, Ina P.↗
NASA NTRS1980

Mechanical energy output of the 5 September 1973 flare

The mechanical energy flux of observed macroscopic mass plasma motions in the solar flare of Sept. 5, 1973, is estimated. Consideration is given to the cool eruptive material in the eruptive filament and large surge as revealed by H alpha observations, the moving emission front seen in Ca II as well as H alpha, the piston-driven shock and mass ejection coronal transient observed in radio spectra and flare core motions, and mechanical energy estimates of 5.6 x 10 to the 29th to 8.9 x 10 to the 30th, 9.0 x 10 to the 29th, 2 x 10 to the 30th (thermal) and 10 to the 31st (magnetic), and 9 x 10 to the 24th erg are obtained, respectively, in agreement with previous estimates. It is concluded that the mechanical energy of large-scale mass motions dominates the radiative output of the flare by more than two orders of magnitude, and that a significant portion of the mechanical energy is in the form of magnetic flux delivered to interplanetary space.

Webb, D. F.↗
NASA NTRS1983

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗
NASA NTRS1995

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(exp -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗
NASA NTRS1995

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(sup -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗
NASA NTRS1993

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗
NASA NTRS2013

Interstellar and Ejecta Dust in the Cas A Supernova Remnant

The ejecta of the Cas A supernova remnant has a complex morphology, consisting of dense fast-moving line emitting knots and diffuse X-ray emitting regions that have encountered the reverse shock, as well as more slowly expanding, unshocked regions of the ejecta. Using the Spitzer 5-35 micron IRS data cube, and Herschel 70, 100, and 160 micron PACS data, we decompose the infrared emission from the remnant into distinct spectral components associated with the different regions of the ejecta. Such decomposition allows the association of different dust species with ejecta layers that underwent distinct nuclear burning histories, and determination of the dust heating mechanisms. Our decomposition identified three characteristic dust spectra. The first, most luminous one, exhibits strong emission features at approx. 9 and 21 micron, and a weaker 12 micron feature, and is closely associated with the ejecta knots that have strong [Ar II] 6.99 micron and [Ar III] 8.99 micron emission lines. The dust features can be reproduced by magnesium silicate grains with relatively low MgO-to-SiO2 ratios. A second, very different dust spectrum that has no indication of any silicate features, is best fit by Al2O3 dust and is found in association with ejecta having strong [Ne II] 12.8 micron and [Ne III] 15.6 micron emission lines. A third characteristic dust spectrum shows features that best matched by magnesium silicates with relatively high MgO-to-SiO2 ratio. This dust is primarily associated with the X-ray emitting shocked ejecta and the shocked interstellar/circumstellar material. All three spectral components include an additional featureless cold dust component of unknown composition. Colder dust of indeterminate composition is associated with [Si II] 34.8 micron emission from the interior of the SNR, where the reverse shock has not yet swept up and heated the ejecta. The dust mass giving rise to the warm dust component is about approx. 0.1solar M. However, most of the dust mass is associated with the unidentified cold dust component. Its mass could be anywhere between 0.1 and 1 solar M, and is primarily limited by the mass of refractory elements in the ejecta. Given the large uncertainty in the dust mass, the question of whether supernovae can produce enough dust to account for ISM dust masses in the local and high-z universe remains largely unresolved.

Interstellar↗
NASA NTRS1984

Type II bursts, shock waves, and coronal transients - The event of 1980 June 29, 0233 UT

The metric Type-II solar burst event of June 29, 1980, is characterized on the basis of spatially resolved radioheliograph observations obtained at Culgoora, Australia, and visible-light observations obtained with the coronograph/polarimeter of the SMM satellite. The data are presented in images, diagrams, and graphs and discussed in detail. The Type-II emission is found to arise in the dense moving material behind the transient loops, which have sky-plane width 0.5 solar radius and line-of-sight depth 0.1-0.4 solar radius. A faint arc observed moving ahead of the transient loops at about 900 km/sec and not associated with the Type-II burst is attributed to a shock front, and the compression ratio and Alfven Mach number of the enhanced-density region are estimated as n2/n1 = 1.3-3 and M(A) = 1.2-3. The ambient material at 3 solar radii is determined to have Alfven speed 250-625 km/sec and magnetic-field strength 50-120 mG. The total mass of the event is calculated as 700 Tg; the total magnetic energy of the loops is (1.5-15) x 10 to the 29th ergs.

Gary, D. E.↗
NASA NTRS2014

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗
NASA NTRS1976

Hydrogen molecules in interstellar space

The most important results obtained from Copernicus satellite observations of ultraviolet H2 absorption features are summarized. The properties of the hydrogen molecule are described along with the radiative processes in which it can be involved. Techniques and equipment used in observing H2 absorption features from space are reviewed. The results are discussed in terms of the fraction of hydrogen in molecular form, rotational excitation of H2, physical conditions within interstellar gas clouds, and approaching clouds characterized by high excitation and negative radial velocity. It is shown that: (1) an appreciable fraction of the hydrogen in denser and thicker clouds is molecular, (2) at least some of the interstellar clouds are far from pressure equilibrium, (3) the approaching clouds are apparently dense thick sheets moving toward the sun as viewed in the local standard of rest, (4) these clouds are probably shock waves produced in the interstellar gas by old supernova remnants or expanding H II regions, and (5) the abundances relative to H of such elements as O, Mg, Si, and Cl in these clouds are nearly equal to the solar values, which is in contrast to the depletions observed in many other interstellar clouds.

Spitzer, L.↗
NASA NTRS1988

Multi-body aircraft with an all-movable center fuselage actively controlling fuselage pressure drag

A multi-body aircraft with an all-movable center fuselage which translates relative to two side fuselages is described. At subsonic and transonic flight the center fuselage is in a forward position. At supersonic speeds the center fuselage moves aft so as to ensure optimum aerodynamic interference at particular Mach numbers. This provides an increased shock strength and greater surface areas so the significant reductions in zero-lift wave drag can be achieved. This concept allows for a significant increase in the wing aspect ratio which would improve high-lift performance at all speeds without incurring a significant supersonic zero-lift wave drag penalty. In addition to an improved low-fineness ratio, high-speed performance is achieved at all speeds and for all flight conditions.

Wood, Richard M.↗
NASA NTRS2008

Overview of the LaNCETS Flight Experiment and CFD Analysis. Supplemental Movies

This presentation focuses on nearfield airborne pressure signatures from the Lift and Nozzle Change Effect on Tail Shocks (LaNCETS) flight test experiment. The primary motivation for nearfield probing in the supersonic regime is to measure the shock structure of aircraft in an ongoing effort to overcome the overland sonic boom barrier for commercial supersonic transportation. LaNCETS provides the opportunity to investigate lift distribution and engine plume effects. During Phase 1 flight testing an F-15B was used to probe the F-15 LaNCETS aircraft in order to validate CFD and pre-flight prediction tools. A total of 29 probings were taken at 40,000 ft. altitude at Machs 1.2, 1.4 and 1.6. LaNCETS Phase 1 flight data are presented as a detailed pressure signature superimposed over a picture of the LaNCETS aircraft. The attenuation of the Inlet-Canard shocks with distance as well as its forward propagation and the coalescence of the noseboom shock are illustrated. A detailed CFD study on a simplified LaNCETS aircraft jet nozzle was performed providing the ability to more accurately capture the shocks from the propulsion system and emphasizing how under- and over-expanding the nozzle affects the existence of shock trains inside the jet plume. With Phase 1 being a success preparations are being made to move forward to Phase 2. Phase 2 will fly similar flight conditions, but this time changing the aircraft's lift distribution by biasing the canard positions, and changing the plume shape by under- and over-expanding the nozzle. Nearfield probing will again be completed in the same manner as in Phase 1. An additional presentation focuses on LaNCETS CFD solution methodology. Discussions highlight grid preprocessing, grid shearing and stretching, flow solving and contour plots. Efforts are underway to better capture the flow features via grid modification and flow solution methodology, which will help to achieve better agreement with flight data. An included CD-ROM provides animations of the nearfield probing procedure and of real data from one of the probings integrated with GPS positional and velocity data. An additional in-flight video from the rear seat of the probing aircraft is also provided.

Cliatt, Larry J., II↗
NASA NTRS1973

Motion of the sources for type 2 and type 4 radio bursts and flare-associated interplanetary disturbances

Shock waves are indirectly observed as the source of type 2 radio brusts, whereas magnetic bottles are identified as the source of moving metric type 4 radio bursts. The difference between the expansion speeds of these waves bottles is examined during their generation and propagation near the flare regions. It is shown that, although generated in the explosive phase of flares, the behavior of the bottles is quite different from that of the waves and that the speed of the former is generally much lower. It is shown that the transit times of disturbances between the sun and the earth give information about the deceleration of shock waves to their local speeds observed near the earth's orbit. A brief discussion is given on the relationship among magnetic bottles, shock waves near the sun, and flare-associated disturbances in interplanetary space.

Sakurai, K.↗