Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK WAVE PROFILE”

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 217 records · Page 12

Shock Radiation Tests for Saturn and Uranus Entry Probes

This paper describes a test series in the Electric Arc Shock Tube at NASA Ames Research Center with the objective of quantifying shock-layer radiative heating magnitudes for future probe entries into Saturn and Uranus atmospheres. Normal shock waves are measured in Hydrogen-Helium mixtures (89:11 by volume) at freestream pressures between 13-66 Pa (0.1-0.5 Torr) and velocities from 20-30 kms. No shock layer radiation is detected within measurement limits below 25 kms, a finding consistent with predictions for Uranus entries. Between 25-30 kms, radiance is quantified from the Vacuum Ultraviolet through Near Infrared, with focus on the Lyman-a and Balmer series lines of Hydrogen. Shock profiles are analyzed for electron number density and electronic state distribution. The shocks do not equilibrate over several cm, and in many cases the state distributions are non-Boltzmann. Radiation data are compared to simulations of Decadal Survey entries for Saturn and shown to be as much as 8x lower than predicted with the Boltzmann radiation model. Radiance is observed in front of the shock layer, the characteristics of which match the expected diffusion length.

Hydrogen↗

Unsteady, one-dimensional gas dynamics computations using a TVD type sequential solver

The efficacy of high resolution convection schemes to resolve sharp gradient in unsteady, 1D flows is examined using the TVD concept based on a sequential solution algorithm. Two unsteady flow problems are considered which include the problem involving the interaction of the various waves in a shock tube with closed reflecting ends and the problem involving the unsteady gas dynamics in a tube with closed ends subject to an initial pressure perturbation. It is concluded that high accuracy convection schemes in a sequential solution framework are capable of resolving discontinuities in unsteady flows involving complex gas dynamics. However, a sufficient amount of dissipation is required to suppress oscillations near discontinuities in the sequential approach, which leads to smearing of the solution profiles.

Thakur, Siddharth↗

Contribution of Spicules to Solar Coronal Emission

Recent high-resolution imaging and spectroscopic observations have generated renewed interest in spicules' role in explaining the hot corona. Some studies suggest that some spicules, often classified as type II, may provide significant mass and energy to the corona. Here we use numerical simulations to investigate whether such spicules can produce the observed coronal emission without any additional coronal heating agent. Model spicules consisting of a cold body and hot tip are injected into the base of a warm (0.5 MK) equilibrium loop with different tip temperatures and injection velocities. Both piston- and pressure-driven shocks are produced. We find that the hot tip cools rapidly and disappears from coronal emission lines such as Fe xii 195 and Fe xiv 274. Prolonged hot emission is produced by preexisting loop material heated by the shock and by thermal conduction from the shock. However, the shapes and Doppler shifts of synthetic line profiles show significant discrepancies with observations. Furthermore, spatially and temporally averaged intensities are extremely low, suggesting that if the observed intensities from the quiet Sun and active regions were solely due to type II spicules, one to several orders of magnitude more spicules would be required than have been reported in the literature. This conclusion applies strictly to the ejected spicular material. We make no claims about emissions connected with waves or coronal currents that may be generated during the ejection process and heat the surrounding area.

Shanwlee Sow Mondal↗

Time-dependent coronal shock acceleration of energetic solar flare particles

A global time-dependent model is presented for the coronal and interplanetary shock acceleration and propagation of energetic solar flare particles. The calculations are carried out to help prove that coronal shock acceleration of solar flare particles is responsible for energetic solar flare event data gathered in interplanetary space. The model is based on the theory of diffusive shock acceleration, and requires particle speeds to be much greater than bulk velocities. Also, sufficient scattering must occur upstream and downstream of the shock for the particle scattering mean free path to be smaller than the characteristic scale lengths, which causes the same particles to encounter the shock repeatedly. A spherically symmetric shock wave is assumed, which leads to the same emission configuration for impulsively and monoenergetically emitted particles. Consideration is given to acceleration by compression at the shock front, adiabatic deceleration in the divergent downstream flow, the temporal evolution of the shock and the three-dimensional geometry of the corona. The model is used to generate normalized proton omnidirectional distributions at 1 AU and at the shock front. The spectral exhibit trends similar to those in observational data, especially proton acceleration times and the proton distribution profiles at 1 AU.

Lee, M. A.↗

Shock-wave structure in a fully ionized gas.

The structure of a steady plane shock in a fully ionized gas has been investigated using the Boltzmann equation with the Gross-Krook type model as the governing equation and the discrete-ordinate method as a tool. The present results agree well with the results obtained by the continuum approach for weak shocks. For strong shocks the present approach gives results that are considerably different from those yielded by the continuum approach, particularly in the high pressure region. Also, the present method gives smooth and continuous shock profiles for high Mach numbers. On the other hand, the results from the continuum approach contain discontinuities in slope. The ion and electron distribution functions across the shock are found to be singly peaked for all cases considered and the effects of the induced electric field is small.

Lu, C. S.↗

Measurements of the boundary layer and near wake of a supercritical airfoil at cruise conditions

The flow behavior within the upper-surface boundary layer and near wake of a supercritical airfoil operating at cruise conditions is discussed. Experimental results obtained from wind-tunnel tests are presented which provide a more detailed description of the flow in these regions than was previously available. Mean streamwise velocity profiles measured by pitot-pressure-probe and laser-velocimeter techniques were found to be in excellent agreement. Other mean-flow properties obtained by the laser-velocimeter technique were the local flow angles in the viscous layers and the static pressures at the edges of the boundary layer and wake. The data set also includes measurements of the turbulence intensity and turbulent Reynolds-stress distributions as obtained by the laser-velocimeter technique. To assess the effects of the shock wave, a less extensive set of measurements was realized at a subcritical test condition. The two test conditions (Mach number at free-stream conditions = 0.72, airfoil section lift coefficient = 0.76 and Mach number of free-stream conditions = 0.5, airfoil section lift coefficient = 0.75) provide a good test for state-of-the-art prediction methods because the upper-surface-boundary layer is separated just upstream of the trailing edge in both cases.

Johnson, D. A.↗

Investigation of a hypersonic crossing shock wave/turbulent boundary layer interaction

A combined theoretical and experimental study is presented for the interaction between crossing shock waves generated by (10 deg, 10 deg) sharp fins and a flat plate turbulent boundary layer at Mach 8.3. The theoretical model is the full 3D mean compressible Reynolds-averaged Navier-Stokes (RANS) equations incorporating the algebraic turbulent eddy viscosity model of Baldwin and Lomax (1978). A grid refinement study indicated that adequate resolution of the flow field has been achieved. Computed results agree well with experiment for surface pressure and surface flow patterns and for pitot pressure and yaw angle profiles in the flow field. The computations, however, significantly overpredict surface heat transfer. Analysis of the computed flow field results indicates the formation of complex streamline and wave structures within the interaction region.

Narayanswami, N↗

Investigation of CAST-10-2/DOA 2 airfoil in NAE high Reynolds number two-dimensional test facility

A common airfoil model with the CAST 10-2/DOA-2 profile and 228 mm (9 inches) chord length was tested. The tests performed in NAE covered the Mach numbers from 0.3 to 0.8 and Reynolds numbers from 10 to 30 million. The model was tested with transition free and with transition fixed at 5 percent chord for both the upper and the lower surfaces. The data obtained were analyzed for the effects of Reynolds number, transition fixing and Mach number. The role of the boundary layer on the displacement effect, the interaction with the shock wave and the trailing edge separation are examined. The results are summarized as follows: (1) the airfoil performance depends strongly on Reynolds number and transition fixing; (2) with transition fixed, the aerodynamic quantities such as lift, pitching moment and drag show a monotonic variation with Reynolds number; (3) with transition free, the aerodynamic quantities vary less regularly with Reynolds number and a slight parametric dependency is shown. The weak dependency is due to the compensatory effect of the forward shift of the transition position and the thinning of the turbulent boundary layer as Reynolds number increases; (4) the shock Mach number and the shock position are weakly dependent on Reynolds number; and (5) the long extent of the laminar boundary layer at transonic speeds reduces the drag appreciably at low Reynolds numbers. The drag bucket around the design Mach number can be observed below Reynolds number 15 million.

Chan, Y. Y.↗

Universal limiter for transient interpolation modeling of the advective transport equations: The ULTIMATE conservative difference scheme

A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.

Leonard, B. P.↗

Optimizing Facility Configurations and Operating Conditions for Improved Performance in the NASA Ames 24 Inch Shock Tube

The Ames Electric Arc Shock Tube (EAST) is a shock tube wherein the driver gas can be heated by an electric arc discharge. The electrical energy is stored in a 1.2 MJ capacitor bank. Four inch and 24 inch diameter driven tubes are available. The facility is described and the need for testing in the 24 inch tube to better simulate low density NASA mission profiles is discussed. Three test entries, 53, 53B and 59, are discussed. Tests are done with air or Mars gas (95.7% CO2/2.7% N2/1.6% Ar) at pressures of 0.01 to 0.14 Torr. Velocities spanned 6.3-9.2 km/s, with a nominal center of 7 km/s. Many facility configurations are studied in an effort to improve data quality. Various driver and driven tube configurations and the use of a buffer section between the driver and the driven tube are studied. Diagnostics include test times, time histories of the shock light pulses and tilts of the shock wave off the plane normal to the tube axis. The report will detail the results of the various trials, give the best configuration/operating conditions found to date and provide recommendations for further improvements. Finally, diaphragm performance is discussed.

Shock Tube↗

Two-equation turbulence modeling for 3-D hypersonic flows

An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.

Bardina, J. E.↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

CW laser strategies for simultaneous, multi-parameter measurements in high-speed gas flows

Strategies utilizing continuous wave (CW) lasers are considered which are capable of simultaneously measuring the flow parameters of velocity, temperature, and pressure at sampling rates exceeding 3 kHz. Velocity is determined from the Doppler shift of the spectral profile, temperature is extracted from intensity ratios of multiple lines, and pressure is measured from either the collision of broadening or the magnitude of absorption. Distinctions between strategies concern the specifics of probe spacies (NO, OH, O2, and H2O) in terms of nominal probe wavelength, equipment, and detection scheme. CW lasers were applied to path-integrated absorption measurements of transient shock-tube flows and spatially resolved laser-induced fluorescence measurements of underexpanded jets.

Di Rosa, M. D.↗

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Kinetics of the reaction of nitric oxide with hydrogen

The reaction of nitric oxide with hydrogen has been studied in the temperature range 2400-4500 K using a shock-tube technique. Mixtures of NO and H2 diluted in argon or krypton were heated by incident shock waves, and the infrared emission from the fundamental vibration-rotation band of NO at 5.3 microns was used to monitor the time-varying NO concentration. The decomposition of nitric oxide behind the shock was found to be modeled well by a fifteen-reaction system. A principal result of the study was the determination of the rate constant k1 for the reaction H + NO yields N + OH, which may be the rate-limiting step for NO removal in some combustion systems. Experimental values of k1 were obtained for each test through comparisons of measured and numerically predicted NO profiles. The data are fit closely by the expression k1 = 1.34 times 10 to the fourteenth power exp(-49 200/RT) cu cm/mole-sec. These data appear to be the first available for this rate constant.

Flower, W. L.↗

Continuous wave laser absorption techniques for gasdynamic measurements in supersonic flows

Line-of-sight measurements of velocity, temperature, pressure, density, and mass flux were performed in a transient shock tube flow using three laser absorption schemes. All methods employed an intracavity-doubled ring dye laser tuned to an OH transition at 306 nm. In the first scheme, the gas was labeled by 193.3-nm excimer photolysis of H2O, and the passage of the generated OH was detected downstream. In the second method, the laser was tuned at a rate of 3 kHz over the R1(7) and R1(11) line pair, and absorption was simultaneously monitored at 90 and 60 deg with respect to the flow. Velocity was determined from the Doppler shift of the profiles and the temperature from the intensity ratio of the lines. Pressure was determined from both the magnitude of absorption and the collisional broadening. In the third method, the laser wavelength was fixed at a single frequency, and a continuous measurement of velocity and pressure was obtained using the signals from the two beam paths. All methods gave results which compare favorably to calculated values.

Davidson, David F.↗

The hypersonic shock wave-turbulent boundary layer interaction generated by a sharp fin at Mach 8.2

A combined experimental and numerical study has been conducted on the hypersonic shock-wave turbulent-boundary layer interaction at Mach 8.2 generated by a single fin of angles alpha = 10 and 15 deg. Three models are considered: (1) the 3D compressible Reynolds-averaged Navier-Stokes (RANS) equations using the k-epsilon turbulence model, (2) the 3D RANS using the Rodi turbulence model, and (3) the conical RANS using the Baldwin-Lomax algebraic turbulence model. The computations are compared with various experimental data. The computations using models (1) and (2) show quantitatively very similar results and very good agreement with experimental data for surface pressure and skin friction. Comparison with boundary layer profiles of pitot pressure and yaw angle are also generally good, but the peak surface heat transfer is overestimated by up to 48 percent. The effect of the laminar boundary layer on the fin is restricted to the immediate vicinity of the fin surface. Conical calculations using model (3) show substantially poorer agreement with experiment.

Knight, Doyle D.↗