Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock Physics”

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 901 records · Page 50

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗

Uncertainty Assessment of Hypersonic Aerothermodynamics Prediction Capability

The present paper provides the background of a focused effort to assess uncertainties in predictions of heat flux and pressure in hypersonic flight (airbreathing or atmospheric entry) using state-of-the-art aerothermodynamics codes. The assessment is performed for four mission relevant problems: (1) shock turbulent boundary layer interaction on a compression corner, (2) shock turbulent boundary layer interaction due a impinging shock, (3) high-mass Mars entry and aerocapture, and (4) high speed return to Earth. A validation based uncertainty assessment approach with reliance on subject matter expertise is used. A code verification exercise with code-to-code comparisons and comparisons against well established correlations is also included in this effort. A thorough review of the literature in search of validation experiments is performed, which identified a scarcity of ground based validation experiments at hypersonic conditions. In particular, a shortage of useable experimental data at flight like enthalpies and Reynolds numbers is found. The uncertainty was quantified using metrics that measured discrepancy between model predictions and experimental data. The discrepancy data is statistically analyzed and investigated for physics based trends in order to define a meaningful quantified uncertainty. The detailed uncertainty assessment of each mission relevant problem is found in the four companion papers.

Bose, Deepak↗

Perspectives on hypersonic viscous and nonequilibrium flow research

An attempt is made to reflect on current focuses in certain areas of hypersonic flow research by examining recent works and their issues. Aspects of viscous interaction, flow instability, and nonequilibrium aerothermodynamics pertaining to theoretical interest are focused upon. The field is a diverse one, and many exciting works may have either escaped the writer's notice or been abandoned for the sake of space. Students of hypersonic viscous flow must face the transition problems towards the two opposite ends of the Reynolds or Knudsen number range, which represents two regimes where unresolved fluid/gas dynamic problems abound. Central to the hypersonic flow studies is high-temperature physical gas dynamics; here, a number of issues on modelling the intermolecular potentials and inelastic collisions remain the obstacles to quantitative predictions. Research in combustion and scramjet propulsion will certainly be benefitted by advances in turbulent mixing and new computational fluid dynamics (CFD) strategies on multi-scaled complex reactions. Even for the sake of theoretical development, the lack of pertinent experimental data in the right energy and density ranges is believed to be among the major obstacles to progress in aerothermodynamic research for hypersonic flight. To enable laboratory simulation of nonequilibrium effects anticipated for transatmospheric flight, facilities capable of generating high enthalpy flow at density levels higher than in existing laboratories are needed (Hornung 1988). A new free-piston shock tunnel capable of realizing a test-section stagnation temperature of 10(exp 5) at Reynolds number 50 x 10(exp 6)/cm is being completed and preliminary tests has begun (H. Hornung et al. 1992). Another laboratory study worthy of note as well as theoretical support is the nonequilibrium flow experiment of iodine vapor which has low activation energies for vibrational excitation and dissociation, and can be studied in a laboratory with modest resources (Pham-Van-Diep et al. 1992).

Cheng, H. K.↗

Fifteen Years of Laboratory Astrophysics at Ames

Tremendous strides have been made in our understanding of interstellar material over the past fifteen years thanks to significant, parallel developments in two closely related areas: observational astronomy and laboratory astrophysics. Fifteen years ago the composition of interstellar dust was largely guessed at, the concept of ices in dense molecular clouds ignored, and the notion of large, abundant, gas phase, carbon-rich molecules widespread throughout the interstellar medium (ISM) considered impossible. Today the composition of dust in the diffuse ISM is reasonably well constrained to cold refractory materials comprised of amorphous and crystalline silicates mixed with an amorphous carbonaceous material containing aromatic structural units and short, branched aliphatic chains. In the dense ISM, these cold dust particles are coated with mixed-molecular ices whose compositions are very well known. Lastly, the signature of carbon-rich polycyclic aromatic hydrocarbons (PAHs), shockingly large molecules by early interstellar chemistry standards, is widespread throughout the ISM. This great progress has only been made possible by the close collaboration of laboratory experimentalists with observers and theoreticians, all with the goal of applying their skills to astrophysical problems of direct interest to NASA programs. Such highly interdisciplinary collaborations ensure fundamental, in depth coverage of the wide-ranging challenges posed by astrophysics. These challenges include designing astrophysically focused experiments and data analysis, tightly coupled with astrophysical searches spanning 2 orders of magnitude in wavelength, and detailed theoretical modeling. The impact of our laboratory has been particularly effective as there is constant cross-talk and feedback between quantum theorists; theoretical astrophysicists and chemists; experimental physicists; organic, physical and petroleum chemists; and infrared and UV/Vis astronomers. In this paper, two examples of the Ames Program will be given. We have been involved in identifying 9 out of the 14 interstellar pre-cometary ice species known, determined their abundances and the physical nature of the ice structure. Details on our ice work are given in the paper by Sandford et al. Our group is among the pioneers of the PAH model. We built the theoretical framework, participated in the observations and developed the experimental techniques needed to test the model. We demonstrated that the ubiquitous infrared emission spectrum associated with many interstellar objects can be matched by laboratory spectra of neutral and positively charged PAHs and that PAHs were excellent candidates for the diffuse interstellar band (DIB) carriers. See Salama et al. and Hudgins et al.

Allamandola, L. J.↗

Coronal Mass Ejections: a Summary of Recent Results

Coronal mass ejections (CMEs) have been recognized as the most energetic phenomenon in the heliosphere, deriving their energy from the stressed magnetic fields on the Sun. This paper highlights some of the recent results on CMEs obtained from the Solar and Heliospheric Observatory (SOHO) and the Solar Terrestrial Relations Observatory (STEREO) missions. The summary of the talk follows. SOHO observations revealed that the CME rate is almost a factor of two larger than previously thought and varied with the solar activity cycle in a complex way (e.g., high-latitude CMEs occurred in great abundance during the solar maximum years). CMEs were found to interact with other CMEs as well as with other large-scale structures (coronal holes), resulting in deflections and additional particle acceleration. STEREO observations have confirmed the three-dimensional nature of CMEs and the shocks surrounding them. The EUV signatures (flare arcades, corona) dimming, filament eruption, and EUV waves) associated with CMEs have become vital in the identification of solar sources from which CMEs erupt. CMEs with speeds exceeding the characteristic speeds of the corona and the interplanetary medium drive shocks, which produce type II radio bursts. The wavelength range of type II bursts depends on the CME kinetic energy: type II bursts with emission components at all wavelengths (metric to kilometric) are due to CMEs of the highest kinetic energy. Some CMEs, as fast as 1600 km/s do not produce type II bursts, while slow CMEs (400 km/s) occasionally produce type II bursts. These observations can be explained as the variation in the ambient flow speed (solar wind) and the Alfven speed. Not all CME-driven shocks produce type II bursts because either they are subcritical or do not have the appropriate geometry. The same shocks that produce type II bursts also produce solar energetic particles (SEPs), whose release near the Sun seems to be delayed with respect to the onset of type II bursts. This may indicate a subtle difference in the acceleration of the ions and 10 keV electrons needed to produce type II bursts. Surprisingly, some shocks lacking type II bursts are associated with energetic storm particle events (ESPs) pointing to the importance of electron escape from the shock for producing the radio emission. CMEs slow down or accelerate in the interplanetary medium because of the drag force, which modifies the transit time of CMEs and shocks. Halo CMEs that appear to surround the occulting disk were known before the SOHO era as occasional events. During the SOHO era, they became very prominent because of their ability to impact Earth and producing geomagnetic storms. Halo CMEs are generally more energetic than ordinary CMEs, which means they can produce severe impact on Earth's magnetosphere. Their origin close to the disk center of the Sun ensures direct impact on the magnetosphere, although their internal magnetic structure is crucial in causing storms. The solar sources of CMEs that produce SEP events at Earth, on the other hand, are generally in the western hemisphere because of the magnetic connectivity. Thus, CMEs are very interesting from the point of view of plasma physics as well as practical implications because of their space weather impact.

Gopalswamy, Nat↗

High-resolution SOFIA/EXES Spectroscopy of SO2 Gas in the Massive Young Stellar Object MonR2 IRS3: Implications for the Sulfur Budget

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the ν3 rovibrational band of SO2 obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO2 gas in the massive YSO MonR2 IRS3. This yields a SO2/H abundance lower limit of 5.6 ± 0.5 × 10(exp -7), or >4% of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of b < 3.20 km s(exp -1). The small line widths and high temperature (Tex = 234 ± 15 K) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov–Gerasimenko), which is followed by SO2 formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO2 from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

Infrared ISM↗

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗

Chandra Imaging of the X-Ray Nebula Powered by Pulsar B1509-58

We present observations with the Chandra X-Ray Observatory of the pulsar wind nebula (PWN) powered by the energetic young pulsar B1509-58. These data confirm the complicated morphology of the system indicated by previous observations, and in addition reveal several new components to the nebula. The overall PWN shows a clear symmetry axis oriented at a position angle 150" +- 5" (north through east), which we argue corresponds to the pulsar spin axis. We show that a previously identified radio feature matches well with the overall extent of the X-ray PWN, and propose the former as the long-sought radio nebula powered by the pulsar. We further identify a bright collimated feature, at least 4' long, lying along the nebula's main symmetry axis; we interpret this feature as a physical outflow from the pulsar, and infer a velocity for this jet of greater than 0.2c. The lack of any observed counterjet implies that the pulsar spin axis is inclined at -30" to the line of sight, contrary to previous estimates made from lower resolution data. We also identify a variety of compact features close to the pulsar. A pair of semicircular X-ray arcs lie 17" and 30" to the north of the pulsar; the latter arc shows a highly polarized radio counterpart. We show that these features can be interpreted as ion compression wisps in a particle-dominated equatorial flow, and use their properties to infer a ratio of electromagnetic to particle energy in pairs at the wind shock 0.005, similar to that seen in the Crab Nebula. We further identify several compact knots seen very close to the pulsar; we use these to infer cr < 0.003 at a separation from the pulsar of 0.1 pc.

Gaensler, B. M.↗

The physics of the oscillating lightning

The lightning is proposed to start in a limited region of electric breakdown field. A growing charge pool in a cumulonimbus cloud is likely to produce such a field in its outer layers. When the field is not aligned with the potential surfaces, the electric force pushes positive charge in one direction and negative in the opposite. Self inductance develops and an embryo open electric circuit is born. Oscillations swiftly grow in strength when potential energy of the cloud is released. The channel extends in steps during those cycle stages where charge is crammed into the channel terminals. Electrostatic consequences of a narrow conductor among the cloud charges are computed. The conductor is charged by influence to such a degree that the charge is bound to leak out to the ambient air and appear as a charge sheath some distance off the channel. Ample energy is released to sustain the oscillations. The period increases with channel length and agrees well with observed time between steps. A ground strike leads to a shock pulse travelling up the channel.

Storebo, Per B.↗

The highly ionized regions of NGC 6302 and NGC 6537

We present (Ne V) 24.32 micrometer and 3426 A emission-line fluxes for NGC 6302 and NGC 6537. These fluxes along with those of (Ne V) 14.3 micrometer allow us to determine the physical conditions within the highly ionized regions of these nebulae. The electron densities found via the (Ne V) line ratios are similar to those found in studies of the low-ionization regions. However, the electron temperatures determined from the (Ne V) line ratios are much greater than those in the low-ionization regions. However, the electron temperatures determined from the (Ne V) line ratios are much greater than those in the low-ionization regions. We use photoionization models to show that these high derived electron temperatures are inconsistent with the assumption that the nebular ionization is produced solely by radiation from the central star. Much of the intense (Ne V) 3426 A emission from these objects must arise from highly ionized shocked material; a conclusion which is in accord with other evidence. The fact that the (Ne V) infrared lines are less affected by the shock means that we can estimate the relative volumes of the shocked and radiatively excited portions of these nebulae. In the most likely scenario, the volume of the highly ionized shocked region is less than 10% that of the highly ionized region formed by photoionization. We note that under such conditions, infrared emission lines are more useful than optical or UV emission lines for abundance studies because they are less affected by shocks. With appropiate corrections to the optical and UV emission-line fluxes, these nebulae can still be described by photoionization models.

Rowlands, N.↗

Plausible mechanisms for rapid acceleration of protons during solar flares

The distinctive physical features of 20 gamma-ray line flares observed during Solar Maximum Mission are discussed. The average energy deposition by nonthermal electrons in the impulsive phase of the flare event of April 27, 1980 is calculated, on the basis of spectral data from the Hinotori satellite. Analysis of the hard X-ray spectral evolution of the flare events of June 21, 1980 and April 27, 1981, showed a very hot (about 10 to the 8th K) thermal component with an emission measure of 3 x 10 to the 47th per cu cm. Some possible mechanisms for the particle acceleration observed during the flare events are discussed, including: the delay of high-energy hard X-rays; first-order Fermi acceleration by shocks propagating in a closed flare loop; and the compression of accelerated particles by upstream and downstream scattering centers around a shock front. In addition to the above mechanisms, sotchastic acceleration by turbulence is identified as a possible mechanism for rapid proton acceleration during solar flare events.

Bai, T.↗

Multi-dimensional MHD simple waves

In this paper we consider a formalism for multi-dimensional simple MHD waves using ideas developed by Boillat. For simple wave solutions one assumes that all the physical variables (the density rho, gas pressure p, fluid velocity V, gas entropy S, and magnetic induction B in the MHD case) depend on a single phase function phi(r,t). The simple wave solution ansatz and the MHD equations then require that the phase function has the form phi = r x n(phi) - lambda(phi)t, where = n(phi) = Delta phi / (absolute value of Delta phi) is the wave normal and lambda(phi) = omega/k = -phi t / (absolute value of Delta phi) is the normal speed of the wave front. The formalism allows for more general simple waves than that usually dealt with in which n(phi) is a constant unit vector that does not vary along the wave front. The formalism has implications for shock formation for multi-dimensional waves.

Webb, G. M.↗

Application of Nonstandard Analysis to the Study of the Shock Structure in a Viscous Heat Conducting Fluid

The study of the shock structure in a viscous heat conducting fluid is an old problem. We study this problem from a novel mathematical point of view. A new class of generalized functions is defined where multiplication of any two functions is allowed with the usual properties. A Heaviside function in this class has the unit jump at occurring on an infinitesimal interval of the nonstandard analysis (NSA) in the halo of . This jump has a smooth microstructure over the infinitesimal interval . From this point of view, we have a new class of Heaviside functions, and their derivatives the Dirac delta functions, which are equivalent when viewed as continuous linear functionals over the test function space of Schwartz. However, they differ in their microstructures which in applications are determined from physics of the problem as shown in our presentation.

F. Farassat↗

The ISM in the M82 starburst

We have observed (O 1) (63 microns) and (Si 2) (35 microns) in the central 700 pc of the starburst galaxy M82. The luminosities in these transitions are 7.1 x 10(exp 7) solar luminosity and 6.2 x 10(exp 7) solar luminosity, respectively, which are each approx. 0.15% of the bolometric luminosity from this region. The ratios of (O 1) line luminosity to (O 3), (Si 2) (35 microns), and to bolometric luminosities in M82 are similar to those in M42, M17, and Sgr A. These similarities, and the association of the bulk of the (O 1) and (Si 2) emission with the ionized emission, suggest that the dominant emission mechanism for (O 1) and (Si 2) in M82 is the same as in these Galactic regions, namely warm gas photodissociated by UV flux from the OB stars responsible for the nearby H 2 regions. We argue that shock or x ray heated gas or H 2 plasma is a minor contributor to the intensities of these fine structure lines. Both the (O 1) (53 microns) and the (Si 2) (35 microns) spectrum show an asymmetric line profile indistinguishable in shape from those of the (O 3) (52 and 88 microns) and (N 3) (57 microns) lines and similar to that of the more extended (C 2) 158 micron line measured previously in M82. We detect two distinct velocity components, which we attribute to emission from two regions at either end of the central bar, where the bar connects to an orbiting torus of neutral gas seen in H 1 and CO J = 1-0. We model separately the two velocity components and derive the physical conditions in these two regions. The clouds in these regions are small, R approx. 1-2 pc, have warm neutral gas surfaces, T approx. 200 K, and are concentrated with volume filling factors of approx. 0.02 and area filling factors of 1-5. The entire central region (R approx. 700 pc) is characterized by a large number, approx. 5 x 10(exp 4), of 2 x 10(exp 3) solar mass clouds with surface densities of approx. 3 x 10(exp 4) cm(exp -3), illuminated by FUV fluxes 10(exp 4) times the average local interstellar value for the Milky Way. These clouds reside in the harsh conditions of a starburst nucleus, with photoevaporation times of 10(exp 6) yr, and collision timescales only about an order of magnitude longer.

Lord, Steven↗

Transient Three-Dimensional Side Load Analysis of a Film Cooled Nozzle

Transient three-dimensional numerical investigations on the side load physics for an engine encompassing a film cooled nozzle extension and a regeneratively cooled thrust chamber, were performed. The objectives of this study are to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet history based on an engine system simulation. Ultimately, the computational results will be provided to the nozzle designers for estimating of effect of the peak side load on the nozzle structure. Computations simulating engine startup at ambient pressures corresponding to sea level and three high altitudes were performed. In addition, computations for both engine startup and shutdown transients were also performed for a stub nozzle, operating at sea level. For engine with the full nozzle extension, computational result shows starting up at sea level, the peak side load occurs when the lambda shock steps into the turbine exhaust flow, while the side load caused by the transition from free-shock separation to restricted-shock separation comes at second; and the side loads decreasing rapidly and progressively as the ambient pressure decreases. For the stub nozzle operating at sea level, the computed side loads during both startup and shutdown becomes very small due to the much reduced flow area.

Wang, Ten-See↗

Metamorphism of cosmic dust: Processing from circumstellar outflows to the cometary regolith

Metamorphism of refractory particles continues in the interstellar medium (ISM) where the driving forces are sputtering by cosmic ray particles, annealing by high energy photons, and grain destruction in supernova generated shocks. Studies of the depletion of the elements from the gas phase of the interstellar medium tell us that if grain destruction occurs with high efficiency in the ISM, then there must be some mechanism by which grains can be formed in the ISM. Most grains in a cloud which collapses to form a star will be destroyed; many of the surviving grains will be severely processed. Grains in the outermost regions of the nebula may survive relatively unchanged by thermal processing or hydration. It is these grains which one hopes to find in comets. However, only those grains encased in ice at low temperature can be considered pristine since a considerable degree of hydrous alteration might occur in a cometary regolith if the comet enters the inner solar system. The physical, chemical and isotopic properties of a refractory grain at each stage of its life cycle will be discussed.

Nuth, Joseph A., III↗

Advanced ion thruster and electrochemical launcher research

The theoretical model of orificed hollow cathode operation predicted experimentally observed cathode performance with reasonable accuracy. The deflection and divergence characteristics of ion beamlets emanating from a two grid optics system as a function of the relative offset of screen and accel grids hole axes were described. Ion currents associated with discharge chamber operation were controlled to improve ion thruster performance markedly. Limitations imposed by basic physical laws on reductions in screen grid hole size and grid spacing for ion optics systems were described. The influence of stray magnetic fields in the vicinity of a neutralizer on the performance of that neutralizer was demonstrated. The ion current density extracted from a thruster was enhanced by injecting electrons into the region between its ion accelerating grids. Theoretical analysis of the electrothermal ramjet concept of launching space bound payloads at high acceleration levels is described. The operation of this system is broken down into two phases. In the light gas gun phase the payload is accelerated to the velocity at which the ramjet phase can commence. Preliminary models of operation are examined and shown to yield overall energy efficiences for a typical Earth escape launch of 60 to 70%. When shock losses are incorporated these efficiencies are still observed to remain at the relatively high values of 40 to 50%.

Wilbur, P. J.↗

A more accurate nonequilibrium air radiation code - NEQAIR second generation

Two experiments, one an equilibrium flow in a plasma torch at Stanford, the other a nonequilibrium flow in a SDIO/IST Bow-Shock-Ultra-Violet missile flight, have provided the basis for modifying, enhancing, and testing the well-known radiation code, NEQAIR. The original code, herein termed NEQAIR1, lacked computational efficiency, accurate data for some species and the flexibility to handle a variety of species. The modified code, herein termed NEQAIR2, incorporates recent findings in the spectroscopic and radiation models. It can handle any number of species and radiative bands in a gas whose thermodynamic state can be described by up to four temperatures. It provides a new capability of computing very fine spectra in a reasonable CPU time, while including transport phenomena along the line of sight and the characteristics of instruments that were used in the measurements. Such a new tool should allow more accurate testing and diagnosis of the different physical models used in numerical simulations of radiating, low density, high energy flows.

Moreau, Stephane↗