SEARCH · Engineering Papers

Results for “shock initiation”

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 1,279 records · Page 71

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Experimental and Computational Studies of the Flow Over a Sting Mounted Planetary Probe Configuration

This paper summarizes the results of a series of experimental studies in the LENS shock tunnel and computations with DSMC and Navier Stokes codes which have been made to examine the aerothermal and flowfield characteristics of the flow over a sting-supported planetary probe configuration in hypervelocity air and nitrogen flows. The experimental program was conducted in the LENS hypervelocity shock tunnel at total enthalpies of 5and 10 MJkg for a range of reservoir pressure conditions from 70 to 500 bars. Heat transfer and pressure measurements were made on the front and rear face of the probe and along the supporting sting. High-speed and single shot schlieren photography were also employed to examine the flow over the model and the time to establish the flow in the base recirculation region. Predictions of the flowfield characteristics and the distributions of heat transfer and pressure were made with DSMC codes for rarefied flow conditions and with the Navier-Stokes solvers for the higher pressure conditions where the flows were assumed to be laminar. Analysis of the time history records from the heat transfer and pressure instrumentation on the face of the probe and in the base region indicated that the base flow was fully established in under 4 milliseconds from flow initiation or between 35 and 50 flow lengths based on base height. The measurements made in three different tunnel entries with two models of identical geometries but with different instrumentation packages, one prepared by NASA Langley and the second prepared by CUBRC, demonstrated good agreement between heat transfer measurements made with two different types of thin film and coaxial gage instrumentation. The measurements of heat transfer and pressure to the front face of the probe were in good agreement with theoretical predictions from both the DSMC and Navier Stokes codes. For the measurements made in low density flows, computations with the DSMC code were found to compare well with the pressure and heat transfer measurements on the sting, although the computed heat transfer rates in the recirculation region did not exhibit the same characteristics as the measurements. For the 10MJkg and 500 bar reservoir match point condition, the measurements and heat transfer along the sting from the first group of studies were in agreement with the Navier Stokes solutions for laminar conditions. A similar set of measurements made in later tests where the model was moved to a slightly different position in the test section indicated that the boundary layer in the reattachment compression region was close to transition or transitional where small changes in the test environment can result in larger than laminar heating rates. The maximum heating coefficients on the sting observed in the present studies was a small fraction of similar measurements obtained at nominally the same conditions in the HEG shock tunnel, where it is possible for transition to occur in the base flow, and in the low enthalpy studies conducted in the NASA Langley high Reynolds number Mach 10 tunnel where the base flow was shown to be turbulent. While the hybrid Navier- StokedDMSC calculations by Gochberg et al. (Reference 1) suggested that employing the Navier- Stokes calculations for the entire flowfield could be seriously in error in the base region for the 10 MJkg, 500 bar test case, similar calculations performed by Cornell, presented here, do not.

Holden, Michael S.↗

NASA NTRS1985

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.↗

NASA NTRS1985

Transonic rotor tip design using numerical optimization

The aerodynamic design procedure for a new blade tip suitable for operation at transonic speeds is illustrated. For the first time, 3 dimensional numerical optimization was applied to rotor tip design, using the recent derivative of the ROT22 code, program R22OPT. Program R22OPT utilized an efficient quasi-Newton optimization algorithm. Multiple design objectives were specified. The delocalization of the shock wave was to be eliminated in forward flight for an advance ratio of 0.41 and a tip Mach number of 0.92 at psi = 90 deg. Simultaneously, it was sought to reduce torque requirements while maintaining effective restoring pitching moments. Only the outer 10 percent of the blade span was modified; the blade area was not to be reduced by more than 3 percent. The goal was to combine the advantages of both sweptback and sweptforward blade tips. A planform that featured inboard sweepback was combined with a sweptforward tip and a taper ratio of 0.5. Initially, the ROT22 code was used to find by trial and error a planform geometry which met the design goals. This configuration had an inboard section with a leading edge sweep of 20 deg and a tip section swept forward at 25 deg; in addition, the airfoils were modified.

Tauber, Michael E.↗

NASA NTRS1990

CFD for hypersonic propulsion

An overview is given of research activity on the application of computational fluid dynamics (CDF) for hypersonic propulsion systems. After the initial consideration of the highly integrated nature of air-breathing hypersonic engines and airframe, attention is directed toward computations carried out for the components of the engine. A generic inlet configuration is considered in order to demonstrate the highly three dimensional viscous flow behavior occurring within rectangular inlets. Reacting flow computations for simple jet injection as well as for more complex combustion chambers are then discussed in order to show the capability of viscous finite rate chemical reaction computer simulations. Finally, the nozzle flow fields are demonstrated, showing the existence of complex shear layers and shock structure in the exhaust plume. The general issues associated with code validation as well as the specific issue associated with the use of CFD for design are discussed. A prognosis for the success of CFD in the design of future propulsion systems is offered.

Povinelli, Louis A.↗

NASA NTRS1991

CFD for hypersonic propulsion

Povinelli, Louis A.↗

NASA NTRS1995

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

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

Chyu, Wei J.↗

NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 22

The Lunar and Planetary Science XXXVI, Part 22 is presented. The topics include: 1) Pressure Histories from Thin and Thick Shock-induced Melt Veins in Meteorites; 2) Nano-structured Minerals as Signature of Microbial Activity; 3) The Insoluble Carbonaceous Material of CM Chondrites as Possible Source of Discrete Organics During the Asteroidal Aqueous Phase; 4) Discovery of Abundant Presolar Silicates in Subgroups of Antarctic Micrometeorites; 5) Characteristics of a Seismometer for the LUNAR-A Penetrator; 6) Heating Experiments of the HaH 262 Eucrite and Implication for the Metamorphic History of Highly Metamorphosed Eucrites; 7) Measurements of Ejecta Velocity Distribution by a High-Speed Video Camera; 8) Petrological Comparison of Mongolian Jalanash Ureilite and Twelve Antarctic Ureilites; 9) Metallographic Cooling Rate of IVA Irons Revisited; 10) Inhomogeneous Temperature Distribution in Chondrules in Shock-Wave Heating Model; 11) Subsurface Weathering of Rocks and Soils at Gusev Crater; 12) Extinct Radioactivities in the Early Solar System and the Mean Age of the Galaxy; 13) Correlation of Rock Spectra with Quantitative Morphologic Indices: Evidence for a Single Rock Type at the Mars Pathfinder Landing Site; 14) Silicon Isotopic Ratios of Presolar Grains from Supernovae; 15) Current Status and Readiness on In-Situ Exploration of Asteroid Surface by MINERVA Rover in Hayabusa Mission; 16) Long Formation Period of Single CAI: Combination of O and Mg Isotope Distribution; 17) Supra-Canonical Initial 26Al/27Al Indicate a 105 Year Residence Time for CAIs in the Solar Proto-Planetary Disk; 18) Evolution of Mercury's Obliquity; 19) First Results from the Huygens Surface Science Package; 20) Polyhedral Serpentine Grains in CM Chondrites; 21) Mountainous Units in the Martian Gusev Highland Region: Volcanic, Tectonic, or Impact Related? 22) Petrography of Lunar Meteorite MET 01210, A New Basaltic Regolith Breccia; 23) Earth-Moon Impacts at ~300 Ma and ~500 Ma Ago; 24) Chemical and Isotopic Analyses of Apollo 16 Glasses: An Integrated Approach; 25) Measuring the Surfaces on Mars; 26) Melting in Martian Snowbanks; 27) The Effect of Ni on Fe-FeS Phase Relations: Implications for the Chemistry of the Martian Core; 28) Shorelines in the Western United States as Analogs for Hypothesized Shoreline Features on Mars; 29) Al and Ti Isotopic Ratios of Presolar SiC Grains of Type Z; 30) Physicochemical Conditions of Clinopyroxene Crystallization in the Parental Bodies of Ordinary Chondrites; 31) Metamorphosed CM and CI Carbonaceous Chondrites Could be from the Breakup of the Same Earth-crossing Asteroid; 32) Aqueous Alteration and Hydrogen Generation on Parent Bodies of Unequilibrated Ordinary Chondrites: Thermodynamic Modeling for the Semarkona Composition; and 33) Preliminary Analysis of Relative Abundances of Hydrothermal Alteration Products in the C1-N10, Y6-N19, and Yax-1_863.51 Impact Melt Samples, Chicxulub Structure, Mexico

Source record↗

NASA NTRS1986

Methane oxidation behind reflected shock waves: Ignition delay times measured by pressure and flame band emission

Ignition delay data were recorded for three methane-oxygen-argon mixtures (phi = 0.5, 1.0, 2.0) for the temperature range 1500 to 1920 K. Quiet pressure trances enabled us to obtain delay times for the start of the experimental pressure rise. These times were in good agreement with those obtained from the flame band emission at 3700 A. The data correlated well with the oxygen and methane dependence of Lifshitz, but showed a much stronger temperature dependence (phi = 0.5 delta E = 51.9, phi = 1.0 delta = 58.8, phi = 2.0 delta E = 58.7 Kcal). The effect of probe location on the delay time measurement was studied. It appears that the probe located 83 mm from the reflecting surface measured delay times which may not be related to the initial temperature and pressure. It was estimated that for a probe located 7 mm from the reflecting surface, the measured delay time would be about 10 microseconds too short, and it was suggested that delay times less than 100 microsecond should not be used. The ignition period was defined as the time interval between start of the experimental pressure rise and 50 percent of the ignition pressure. This time interval was measured for three gas mixtures and found to be similar (40 to 60 micro sec) for phi = 1.0 and 0.5 but much longer (100 to 120) microsecond for phi = 2.0. It was suggested that the ignition period would be very useful to the kinetic modeler in judging the agreement between experimental and calculated delay times.

Brabbs, T. A.↗

DOE OSTI2021

Combined microwave and laser Rayleigh scattering diagnostics for pin-to-pin nanosecond discharges

In this work, the temporal decay of electrons produced by an atmospheric pin-to-pin nanosecond discharge operating in the spark regime was measured via a combination of microwave Rayleigh scattering (MRS) and laser Rayleigh scattering (LRS). Due to the initial energy deposition of the nanosecond pulse, a variation in the local gas density occurs on the timescale of electron decay. Thus, the assumption of a constant collisional frequency is no longer applicable when electron number data are extracted from MRS measurements. Furthermore, to recalibrate MRS measurements throughout the electron decay period, temporally resolved LRS measurements of the local gas density were performed over the event duration. The local gas density was calculated to be 30% of the ambient level during the later stages of electron decay, and it recovers at about 1 ms after discharge. A shock front traveling approximately 500 m/s was additionally observed. Coupled with plasma volume calibration via temporally resolved intensified charge-coupled device imaging, the corrected decay curves of the electron number and electron number density are presented with a measured peak electron number density of 4.5 × 1015 cm -3 and a decay rate of ~(0.1–0.35) × 10 7 s -1 . A hybrid MRS and LRS diagnostic technique can be applied for a broad spectrum of atmospheric-pressure microplasmas where a variation in the gas number density is expected due to energy deposition in the discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1982

A first-order estimate of shock heating and vaporization in oceanic impacts

The vaporization of water in oceanic impacts of asteroids or comets of multikilometer dimensions is estimated by a semianalytical modeling of impact heating and shock isobar geometry that is based on computer code calculations. The mass of water vaporized in an infinitely deep ocean by the impact of a 10 km diameter asteroid at 25 km/sec (these values have been proposed for the Cretaceous/Tertiary extinction bolide) is approximately equal to the total mass of water vapor present in the earth's atmosphere, and 3-4 orders of magnitude larger than the mass of water vapor in the stratosphere. For projectiles of this size, however, the finite depth of the ocean becomes significant and may considerably reduce the amount of water vapor initially generated by the impact. Climatological models and extinction scenarios invoking the effects of impact-generated water vapor may critically depend on the a priori ambiguous details of the hypothesized impact.

Croft, S. K.↗

NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

NASA NTRS1973

The initial stage of development of type IV radio bursts and the relation to expanding magnetic bottles

Using the observed data for wideband type IV solar radio bursts, the onset time differences between the microwave and metric frequencies and the peak flux intensities of the metric component are analyzed as a function of the longitudinal position of the associated flares on the solar disk. It is shown that this time difference is dependent on the position of the associated flare and that the peak flux intensity reaches maximum when a flare occurs in the region from 10 to 40 deg west of the central meridian of the solar disk. These results are explained by taking into account the eastward expansion of magnetic bottles which trap mildly relativistic electrons responsible for type IV bursts. The relation between these magnetic bottles and shock waves which excite type II radio bursts is discussed.

Sakurai, K.↗

DOE OSTI2023

The effect of microstructure on the dynamic shock response of 1045 steel

Herein, the effect of microstructure on the shock response of 1045 steel is investigated via plate impact experiments and postmortem characterization. Three unique microstructures are explored: ferrite-pearlite, martensite, and ferrite with spheroidal cementite (i.e. spheroidized). Two spall recovery experiments, at approximate peak pressures of 3.2 and 3.5 GPa, are conducted to assess the Hugoniot elastic limit (HEL), spall strength, and damage morphology of the various microstructures. The ferrite-pearlite and martensite microstructures exhibit continuous yielding at both quasi-static and dynamic rates, while the spheroidized condition displays discontinuous yielding. Discontinuous yielding of the spheroidized microstructure is attributed to a combined low initial dislocation density coupled with a low dislocation nucleation rate. The spall strength of ferrite-pearlite is consistently lower than the spheroidized microstructure, attributed to elongated cementite that is more susceptible to cracking than more spherical cementite precipitates. Despite a high density of boundaries, martensite exhibits the highest spall strength. A large percentage of the boundaries within the martensite microstructure are found to be low energy (i.e. Σ3 or low angle), and are thus less susceptible to spall damage. Overall, the high spall strength of martensite is likely linked to traditional strengthening mechanisms that limit dislocation motion.

36 MATERIALS SCIENCE↗

NASA NTRS1982

Cometary and meteorite swarm impact on planetary surfaces

The impact-induced deformation from hypothetical cometary objects having initial densities in the 0.01 to 1 g/cu cm range and heats of vaporization in the approximately 2 kJ/g (corresponding to water) to approximately 10 to the 7th J/g range is examined for impacts in the 5 to 45 km/s range. Even though the direct effect of an atmosphere is neglected, the atmosphere may in fact cause a cometary object to break up into a shower or equivalent very porous impactor. Besides examining the partitioning of impact energy into internal energy of the impacted planet and impacting cometary material, calculations are made of the relative efficiency of shock-induced melting and vaporization by comets on planetary surface materials and the mass loss from a given planet for various escape velocities.

Okeefe, J. D.↗

NASA NTRS1985

Complex structure of the supernova remnant HB 3

HB 3 is an old, large (84 pc diameter) supernova remnant associated with the W3 H II region/molecular cloud complex. Observations of the imaging proportional counter (IPC) onboard the Einstein X-ray astronomy satellite have been reprocessed to yield a contour map of X-ray brightness and spectra of various regions of this remnant. The measured IPC flux is 2.4 x 10 to the -11th ergs per sq cm per s, giving a 0.2-4 keV luminosity of 1.6 x 10 to the 35th ergs/s for a column densityof 6 x 10 to the 21st per sq cm. The measured X-ray temperatures reveal a decrease from center to limb of the remnant of 1-0.3 keV. HB 3 is in the late adiabatic blast-wave phase of evolution, 30,000 to 50,000 yr old and with an initial blast energy of 3 x 10 to the 50th ergs. The X-ray map is compared with available radio and optical images. In X-rays, HB 3 has two components - a diffuse emission inside the 84 pc radio remnant and a ring of emission at the center of 30 pc in diameter. The diffuse emission is similar to that from other supernova remnants which are moderately obscured (column density, nH approximately 10 to the 22nd per sq cm). Three possibilities for the origin of the ring are explored: (1) a second supernova remnant, (2) a shocked shell in the interstellar medium surrounding HB 3, and (3) reverse-shock heated ejecta. There is no hot neutron star within the remnant.

Leahy, D. A.↗

NASA NTRS2011

Coronal Mass Ejections and Their Heliospheric Consequences

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. The Solar and Heliospheric Observatory (SOHO) mission has accumulated a wealth of information on CMEs with unprecedented extent, clarity, and uniformity. The CME data base has been useful in confirming results from past missions and obtaining a number of new results on the initiation, propagation, and heliospheric consequences of CMEs. This paper highlights some of the recent results on the morphological, physical, and kinematic properties of CMEs observed during solar cycle 23. White-light signatures of CMEs and their utility in inferring coronal magnetic fields near the Sun will be discussed. Presence of coronal holes near eruption regions have been shown to significantly modify the trajectories of CMEs. New results on the interaction of CMEs with large-scale structures such as coronal holes and other CMEs will be presented. Finally, the CME consequences in the heliosphere such as interplanetary shocks, type II radio bursts, energetic particles, geomagnetic storms, and cosmic ray modulation will be discussed

Gopalswamy, Nat↗

NASA NTRS1987

Computation of asymmetric supersonic flows around cones at large incidence

The Schiff-Steger parabolized Navier-Stokes (PNS) code has been modified to allow computation of conical flowfields around cones at high incidence. The improved algorithm of Degani and Schiff has been incorporated with the PNS code. This algorithm adds the cross derivative and circumferential viscous terms to the original PNS code and modifies the algebraic eddy viscosity turbulence model to take into account regions of so called cross-flow separation. Assuming the flowfield is conical (but not necessarily symmetric) a marching stepback procedure is used: the solution is marched one step downstream using improved PNS code and the flow variables are then scaled to place the solution back to the original station. The process is repeated until no change in the flow variables is observed with further marching. The flow variables are then constant along rays of the flowfield. The experiments obtained by Bannik and Nebbeling were chosen as a test case. In these experiments a cone of 7.5 deg. half angle at Mach number 2.94 and Reynolds number 1.372 x 10(7) was tested up 34 deg. angle of attack. At high angle of attack nonconical asymmetric leeward side vortex patterns were observed. In the first set of computations, using an earlier obtained solution of the above cone for angle of attack of 22.6 deg. and at station x=0.5 as a starting solution, the angle of attack was gradually increased up to 34 deg. During this procedure the grid was carfully adjusted to capture the bow shock. A stable, converged symmetric solution was obtained. Since the numerical code converged to a symmetric solution which is not the physical one, the stability was tested by a random perturbation at each point. The possible effect of surface roughness or non perfect body shape was also investigated. It was concluded that although the assumption of conical viscous flows can be very useful for certain cases, it can not be used for the present case. Thus the second part of the investigation attempted to obtain a marching (in space) solution with the PNS method using the conical solution as initial data. Finally, the solution of the full Navier-Stokes equations was carried out.

Degani, David↗

NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗