Engineering Papers⌕ Search

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.

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

At least 379 records · Page 21

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

Sm-Nd, Rb-Sr, and Mn-Cr Ages of Yamato 74013

Yamato 74013 is one of 29 paired diogenites having granoblastic textures. The Ar-39 - Ar-40 age of Y-74097 is approximately 1100 Ma. Rb-Sr and Sm-Nd analyses of Y-74013, -74037, -74097, and -74136 suggested that multiple young metamorphic events disturbed their isotopic systems. Masuda et al. reported that REE abundances were heterogeneous even within the same sample (Y-74010) for sample sizes less than approximately 2 g. Both they and Nyquist et al. reported data for some samples showing significant LREE enrichment. In addition to its granoblastic texture, Y-74013 is characterized by large, isolated clots of chromite up to 5 mm in diameter. Takeda et al. suggested that these diogenites originally represented a single or very small number of coarse orthopyroxene crystals that were recrystallized by shock processes. They further suggested that initial crystallization may have occurred very early within the deep crust of the HED parent body. Here we report the chronology of Y-74013 as recorded in chronometers based on long-lived Rb-87 and Sm-147, intermediate- lived Sm-146, and short-lived Mn-53.

Nyquist, L. E.↗

Shaping the Dust Mass-Star-Formation Rate Relation

There is a remarkably tight relation between the observationally inferred dust masses and star-formation rates (SFRs) of Sloan Digital Sky Survey galaxies, M dust ∝ SFR 1.11(exp). Here we extend the M dust-SFR relation to the high end and show that it bends over at very large SFRs (i.e., dust masses are lower than predicted for a given SFR). We identify several distinct evolutionary processes in the diagram: (1) a star-bursting phase in which dust builds up rapidly at early times. The maximum attainable dust mass in this process is the cause of the bend-over of the relation. A high dust-formation efficiency, a bottom-light initial mass function, and negligible supernova shock dust destruction are required to produce sufficiently high dust masses. (2) A quiescent star-forming phase in which the subsequent parallel decline in dust mass and SFR gives rise to the M dust-SFR relation, through astration and dust destruction. The dust-to-gas ratio is approximately constant along the relation. We show that the power-law slope of the M dust-SFR relation is inversely proportional to the global Schmidt-Kennicutt law exponent (i.e., ~0.9) in simple chemical evolution models. (3) A quenching phase which causes star formation to drop while the dust mass stays roughly constant or drops proportionally. Combined with merging, these processes, as well as the range in total baryonic mass, give rise to a complex population of the diagram which adds significant scatter to the original M dust-SFR relation. (4) At very high redshifts, a population of galaxies located significantly below the local relation is predicted.

Shaping↗

Development of Langmuir and Spatially Resolved OES Diagnostics for Arc Jet Shock and Boundary Layer Measurements Relevant to Atmospheric Re-Entry Environments

The second-generation miniature Arc Jet Research Chamber (mARC) ground testing facility at NASA Ames requires physical characterization of its supersonic plasma flow to verify numerical models and material sample testing. Two types of plasma diagnostics will be utilized: Langmuir probes and optical emission spectroscopy (OES). Several single and triple Langmuir probes have been constructed and will be inserted parallel to the flow via a sweep arm to measure post-shock electron temperature and number density along radial profiles at different axial distances. Due to the high heat flux (200-3000 W/cm2) sweep durations are necessarily short (< 0.1 s). Measuring these properties near the probe surface provide characterization of the post-shock environment. This is an important metric for validating coupled CFD-material codes, and will provide insight into longstanding thermocouple anomalies attributed to charging of heatshield surfaces. An OES system has also been assembled to obtain spatially resolved density and temperature for various flow species. The system consists of a convex lens that forms an image of the plasma onto a 16-channel linear fiber array, which will simultaneously capture the spectra from sixteen different locations. The light emission from the plasma is then focused by a spherical mirror into a spectrometer and recorded with a camera. These measurements non-intrusively provide a new spatially resolved flow characterization capability. This will initially be used to resolve issues with shock radiation interference in pyrometer measurements of material surface temperature. The design of each diagnostic and any initial results will be discussed.

Ethan Main-Loam Leong↗

Experimental observation of transition behavior on a flat plate

In studying transition behavior a shock tube and tunnel were used to produce high temperatures, and thin-film platinum heat gauges were used to measure local heat flux as well as to detect the transition of the laminar boundary layer over a flat plate and a cone. Initial investigations were conducted in the hypersonic shock tunnel to obtain high-temperature information for the development of an ICBM nose cone. Shock Mach numbers as large as 50 with a temperature of 15,000 K after the incident wave were produced in the driven tube. Shock tubes are used to investigate the heat transfer over various surfaces to 2500 K for the development of future gas turbines.

Nagamatsu, H. T.↗

The Formation of Relativistic Jets from Kerr Black Holes

We have performed the first fully three-dimensional general relativistic magnetohydrodynamics (GRMHD) simulation for Schwarzschild and Kerr black holes with a free falling corona and thin accretion disk. The initial simulation results with a Schwarzschild metric show that a jet is created as in the previous axisymmetric simulations with mirror symmetry at the equator. However, the time to form the jet is slightly longer than in the 2-D axisymmetric simulation. We expect that the dynamics of jet formation are modified due to the additional freedom in the azimuth dimension without axisymmetry with respect to the Z axis and reflection symmetry respect to the equatorial plane. The jet which is initially formed due to the twisted magnetic fields and shocks becomes a wind at the later time. The wind flows out with a much wider angle than the initial jet. The twisted magnetic fields at the earlier time were untwisted and less pinched. The accretion disk became thicker than the initial condition. Further simulations with initial perturbations will provide insights for accretion dynamics with instabilities such as magneto-rotational instability (MRI) and accretion-eject instability (AEI). These instabilities may contribute to variabilities observed in microquasars and AGN jets.

Nishikawa, K.-I.↗

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

Hybrid Simulations of Interstellar Pickup Ion Acceleration at the Solar Wind Termination Shock

Hybrid (kinetic ion/fluid electron) simulations have been used to study self-consistently the injection and initial stages of acceleration of interstellar pickup ions at the solar wind termination shock. Results are presented from one-dimensional simulations of high Mach number oblique (40 deg-50 deg) shocks with a 10% population of interstellar pickup hydrogen. In these simulations, the pickup ions, the solar wind ions, the shock fields and the waves are all treated self-consistently. Pickup ions reflected by the shock excite large amplitude (DELTA.B / B approximately 0.3) upstream magnetosonic waves. These waves, in turn, scatter the pickup hydrogen, as expected in the diffusive shock acceleration process. The spectrum of excited waves broadens in time. We find that, for the parameters studied, the termination shock efficiently injects and accelerates the interstellar pickup hydrogen, and thus this work represents a microscopic verification of this anomalous cosmic ray injection mechanism. Comparison of injected fluxes from the simulations with Voyager fluxes extrapolated to the termination shock suggests that injection at the termination shock is probably greater than injection by interplanetary shocks. A study of accelerated pickup ion orbits shows that the energy gain comes predominantly from shock drift acceleration in the shock front, with the upstream waves aiding the acceleration by allowing multiple encounters with the shock.

Liewer, P. C.↗

High-temperature shock formation of N2 and organics on primordial Titan

Theoretical models suggest that the initial form of nitrogen in Titan's atmosphere may have been NH3. The possible importance of strong shocks produced during high-velocity impacts accompanying the late stages of accretion are investigated as a method for converting NH3 to N2, of which Titan's atmosphere is now primarily composed. The focused beam of a high-power laser is used to simulate the effects of an impact in Titan's atmosphere. For mixtures of 10, 50, and 90 percent NH3, yields of 0.25, 1, and 6 x 20 to the 17th molecules of N2 per joule, respectively, were obtained. It is also found that the yield of HCN is comparable to that for N2. Several other hydrocarbons are produced, many with yields in excess of theoretical high-temperature-equilibrium models. The above yields, when combined with models of the satellite's accretion, result in a total N2 production comparable to that present in TItan's atmosphere and putative ocean.

Mckay, Christopher P.↗

Spall and subsequent recompaction of copper under shock loading

Research was undertaken to study spall and subsequent recompaction of oxygen-free high-thermal-conductivity copper using a single-stage large-bore light gas gun capable of planar impacts. Gun experiments were conducted that produced an initial spall in the target with a subsequent recompaction of the spall damage/layer by the use of layered impactors. Symmetric spall experiments at similar conditions were also conducted as a control to the recompaction experiments. Photonic Doppler velocimetry was used to obtain the velocity history of the target’s back surface; these velocity histories were analyzed, and the results were compared to numerical simulations. After recovery of selected samples, the microstructure was analyzed using optical imaging microscopy and electron backscatter diffraction. Recompaction waves were clearly observed in the time-resolved data obtained at the back surface of the targets. Two different series of experiments were performed to understand the effect of damage morphology on the stress required for spall recompaction. In the initial series, the peak pressure of the first shock was systematically increased to alter the amount of damage in the material. While in the second series, the peak pressure of the second shock was increased, keeping the magnitude of the first shock relatively constant to understand the effect of peak stress on the recompaction process. Findings clearly show that recompaction results in a perturbed band of the microstructure at the location of the expected spall plane in all cases except where the peak stress of the second shock is lower than the critical value.

, recompaction↗

Fuel convergence sensitivity in indirect drive implosions

In inertial confinement fusion experiments at the National Ignition Facility, a spherical shell of deuterium–tritium fuel is imploded in an attempt to reach the conditions needed for fusion, self-heating, and eventual ignition. Since theory and simulations indicate that ignition efficacy in 1D improves with increasing imploded fuel convergence ratio, it is useful to understand the sensitivity of the scale-invariant fuel convergence on all measurable or inferable 1D parameters. In this paper, we develop a simple isobaric and isentropic compression scaling model incorporating sensitivity to the in-flight adiabat inferred from shock strengths, to measured implosion velocity, and to known initial ablator and fuel aspect ratio and mass ratio. The model is first benchmarked to 1D implosion simulations spanning a variety of relevant implosion designs. We then use the model to compare compressibility trends across all existing indirect-drive layered implosion data from the facility spanning three ablators [CH, carbon (C), and Be], for which in-flight fuel adiabats varied from 1.6 to 5 by varying the number of drive shocks from 2 to 4, peak implosion velocities varied by 1.4×, capsule radii by 50%, and initial fuel aspect ratios by 1.4×. We find that the strength of the first shock is the dominant contributor setting the maximum fuel convergence. We also observe additional sensitivities to successive shock strengths and fuel aspect ratios that improve the agreement between the expected and measured compression for carbon and Be designs with adiabats above 3. A principal finding is that the adiabat 2.5 C-shell designs exhibit less convergence than CH-shell designs of similar inferred in-flight adiabat.

Landen, O. L. (ORCID:0000000214998217)↗

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

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

Miller, C. G., III↗

Unusual locations of Earth's bow shock on September 24 - 25, 1987: Mach number effects

International Sun Earth Explorer 1 (ISEE 1) and Interplanetary Monitoring Platform 8 (IMP 8) data are used to identify 19 crossings of Earth's bow shock during a 30-hour period following 0000 UT on September 24, 1987. Apparent standoff distances for the shock are calculated for each crossing using two methods and the spacecraft location; one method assumes the average shock shape, while the other assumes a ram pressure-dependent shock shape. The shock's apparent standoff distance, normally approximately 14 R(sub E), is shown to increase from near 10 R(sub E) initially to near 19 R(sub E) during an 8-hour period, followed by an excursion to near 35 R(sub E) (where two IMP 8 shock crossings occur) and an eventual return to values smaller than 19 R(sub E). The Alfven M(sub A) and fast magnetosonic M(sub ms). Mach numbers remain above 2 and the number density above 4/cu cm for almost the entire period. Ram pressure effects produce the initial near-Earth shock location, whereas expansions and contractions of the bow shock due to low Mach number effects account, qualitatively and semiquantitatively, for the timing and existence of almost all the remaining ISEE crossings and both IMP 8 crossings. Significant quantitative differences exist between the apparent standoff distances for the shock crossings and those predicted using the observed plasma parameters and the standard model based on Spreiter et al.'s (1966) gasdynamic equation. These differences can be explained in terms of either a different dependence of the standoff distance on Mach number at low M(sub A) and M(sub ms), or variations in shock shape with M(sub A) and M(sub ms) (becoming increasingly "puffed up" with decreasing M(sub A) and M(sub ms), as expected theoretically), or by a combination of both effects.

Cairns, Iver H.↗

Heterogeneous shock-induced thermal radiation in minerals

A 500-channel optical imaging intensifying and spectral digital recording system is used for recording the shock-induced radiation emitted from 406 to 821 nm from transparent minerals during the time interval that a shock wave propagates through the sample. The initial results obtained for single crystals of gypsum, calcite and halite in the 30 to 40 GPa (300 to 400 kbar) pressure range reveal grey-body emission spectra corresponding to temperatures in the 3000 to 4000 K range and emissivities ranging from 0.003 to 0.02. With gypsum and calcite, distinctive line spectra are superimposed on the thermal radiation. The observed color temperatures are greater than the Hugoniot temperature by a factor of 2 to 10; this is calculable on the basis of continuum thermodynamics and equation of state models for the shock states achieved in the three minerals. These observed high temperatures are thought to be real. It is concluded that a large number of closed spaced high temperature shear-band regions are being detected immediately behind the shock front.

Kondo, K.-I.↗

Mass-loss rates, ionization fractions, shock velocities, and magnetic fields of stellar jets

In this paper we calculate emission-line ratios from a series of planar radiative shock models that cover a wide range of shock velocities, preshock densities, and magnetic fields. The models cover the initial conditions relevant to stellar jets, and we show how to estimate the ionization fractions and shock velocities in jets directly from observations of the strong emission lines in these flows. The ionization fractions in the HH 34, HH 47, and HH 111 jets are approximately 2%, considerably smaller than previous estimates, and the shock velocities are approximately 30 km/s. For each jet the ionization fractions were found from five different line ratios, and the estimates agree to within a factor of approximately 2. The scatter in the estimates of the shock velocities is also small (+/- 4 km/s). The low ionization fractions of stellar jets imply that the observed electron densities are much lower than the total densities, so the mass-loss rates in these flows are correspondingly higher (approximately greater than 2 x 10(exp -7) solar mass/yr). The mass-loss rates in jets are a significant fraction (1%-10%) of the disk accretion rates onto young stellar objects that drive the outflows. The momentum and energy supplied by the visible portion of a typical stellar jet are sufficient to drive a weak molecular outflow. Magnetic fields in stellar jets are difficult to measure because the line ratios from a radiative shock with a magnetic field resemble those of a lower velocity shock without a field. The observed line fluxes can in principle indicate the strength of the field if the geometry of the shocks in the jet is well known.

Hartigan, Patrick↗

Initial dynamic load estimates during configuration design

This analysis includes the structural response to shock and vibration and evaluates the maximum deflections and material stresses and the potential for the occurrence of elastic instability, fatigue and fracture. The required computations are often performed by means of finite element analysis (FEA) computer programs in which the structure is simulated by a finite element model which may contain thousands of elements. The formulation of a finite element model can be time consuming, and substantial additional modeling effort may be necessary if the structure requires significant changes after initial analysis. Rapid methods for obtaining rough estimates of the structural response to shock and vibration are presented for the purpose of providing guidance during the initial mechanical design configuration stage.

Schiff, Daniel↗

The Study on the Shock Wave Propagation Rule of a Gas Explosion in a Gas Compartment

Combined with the k-ε turbulence model of general application, a refined finite element model of a utility tunnel’s gas compartment filled with the methane/air mixture is developed. A series of analyses are made by using the powerful industry-leading computational fluid dynamics (CFD) software flame acceleration simulator (FLACS) to study the shock wave propagation rule in the gas compartment. The longitudinal and transversal distribution laws of the explosion shock wave are gained taking into consideration the spatial characteristics of the gas compartment. The influences of a few parameters, such as initial conditions and section size of the gas compartment, on the shock wave propagation rule are further discussed. The basic procedure for predicting the peak pressure of the blast wave is provided by considering the initial conditions and the gas compartment, and the corresponding injury effect of the explosion wave on the living beings is assessed. The investigation demonstrates that the peak pressure by the coupled effect between the initial conditions is significantly influenced, especially at the upper and lower gas explosion limits. The peak pressure increases gradually as the width or height increases, and both basically meet the linear relation. The proposed method can forecast the peak pressure of the explosion shock wave in the gas compartment accurately. According to the peak pressure longitudinal and transversal distributions of the blast wave, the peak pressure is far greater than the killing pressure threshold in the underground and closed space; consequently, it is not safe for the living beings in the gas compartment.

Lin, Shu-Chao↗