Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗
NASA NTRS1993

An investigation of shock wave turbulent boundary layer interaction with bleed through slanted slots

Flow-field characteristics are simulated numerically in an oblique shock wave/turbulent boundary layer interactions with six different bleed slot configurations. The strong conservation-law form of the 2D compressible Navier-Stokes equations and the k-epsilon equations are solved throughout the interaction region and inside the bleed slot. The computed results are presented for a normal and 20-deg slanted bleed slots at three different locations, upstream, across and downstream of the impingement point of an oblique shock of sufficient strength to cause boundary layer separation without bleed. The detailed flow characteristics in the interaction zone and inside the bleed slot are compared for the different bleed slot configurations. The resulting surface pressure and shear stress distributions as well as the boundary layer characteristics downstream of the interaction region are also presented for the mix bleed configurations at different bleed mass flows up to choking.

Hamed, A.↗
NASA NTRS1972

Pioneer 7 observations of the August 29, 1966, interplanetary shock-wave ensemble.

The detailed plasma (except for electrons) and magnetic field data from Pioneer 7 are presented for a proposed double-shock ensemble. The Ames Research Center plasma probe data are compared with that of the MIT probe. These data are compared with MHD theory to show that the discontinuities themselves satisfy the magnetic Hugoniot equations. The data are also compared with the theoretical predictions for the double-shock structure. The theory requires as input only the forward shock velocity, the original ambient solar wind conditions, and a reasonable assumption regarding the density and azimuthal magnetic field in the 'new' solar wind which is detected after passage of the reverse shock. A unique solution is then obtained for the theoretical piston and reverse shock positions as well as for the temporal plasma behavior for both infinite and finite electrical conductivities.

Dryer, M.↗
NASA NTRS1976

Curvature effects on stagnation radiative transfer in shock wave around blunt axisymmetric bodies

During hypersonic entry of a probe into the atmospheres of outer planets, the radiative energy transport in the shock-heated layer is of great significance in computing the heating to the heat shield. The so-called tangential slab approximation is frequently employed to evaluate the radiative flux to the surface of entry probes. A simple assessment of the accuracy of the tangential slab approximation at the stagnation point is given here for blunt axisymmetric bodies of various shock shapes which include spheres, paraboloids, hyperboloids and ellipsoids of different shock parameters.

Poon, P. T. Y.↗
NASA NTRS1984

Simulation of blunt-fin-induced shock wave and turbulent boundary-layer interaction

A supersonic flow over a blunt fin mounted on a flat plate is numerically simulated. The fin shock causes the boundary layer to separate and results in a complicated, three-dimensional shockwave and boundary-layer interaction. The computed result is in good agreement with the measured pressure on the fin and the flat plate. The main features, such as peak pressure on the fin leading edge and a double peak pressure on the plate, are closely predicted. The role of the horseshoe vortex is discussed. The vortex leads to the development of high-speed flow, and, hence, low-pressure regions on the fin and the plate. Different thicknesses of the incoming boundary layer have been studied. Varying the thicknesses by an order of magnitude shows that the size of the horseshoe vortex and therefore the spatial extent of the interaction are inviscid-dominated, and are weakly dependent on the Reynolds number. Colored graphics are used to show details of the interaction flow field.

Hung, C.-M.↗
NASA NTRS1988

Essentially non-oscillatory shock capturing methods applied to turbulence amplification in shock wave calculations

ENO (essentially non-oscillatory) schemes can provide uniformly high order accuracy right up to discontinuities while keeping sharp, essentially non-oscillatory shock transitions. Recently, an efficient implementation of ENO schemes was obtained based on fluxes and TVD Runge-Kutta time discretizations. The resulting code is very simple to program for multi-dimensions. ENO schemes are especially suitable for computing problems with both discontinuities and fine structures in smooth regions, such as shock interaction with turbulence, for which results for 1-D and 2-D Euler equations are presented. Much better resolution is observed by using third order ENO schemes than by using second order TVD schemes for such problems.

Osher, Stanley↗
NASA NTRS1990

Color surface-flow visualization of fin-generated shock wave boundary-layer interactions

Kerosene-lampblack mixtures with addition of a ground colored chalk were used in an experiment on visualizing surface flows of swept shock boundary-layer interactions. The results show that contrasting colors intensify the visualization of different regions of the interaction surface, and help the eye in following the fine streaks to locate the upstream influence. The study confirms observations of the separation occurring at shock strength below accepted values. The superiority of the reported technique over the previous monochrome technique is demonstrated.

Lu, F. K.↗
NASA NTRS1975

Collisionless shock waves in space - A very high beta structure

Measurements from six OGO-5 particle and field experiments are used to examine the structure of the earth's bow shock during a period of extremely high beta (the ratio of plasma thermal to magnetic energy density), as determined from simultaneous measurements of the upstream plasma on board the HEOS satellite. Even though the interplanetary field is nearly perpendicular to the shock normal, the shock is extremely turbulent. Large field increases are observed up to a factor of 20 above the upstream values. Ahead of these large enhancements, smaller magnetic effects accompanied by electrostatic noise, electron heating, and ion deflection are observed for several minutes. These observations suggest that a steady-state shock may not be able to form at very high beta. Further, they show that while the magnetic energy density may be relatively unimportant in the upstream flow, it can become very significant within the shock structure, and hence the magnetic field should not be ignored in theoretical treatments of very high beta shocks.

Formisano, V.↗
NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗
NASA NTRS1973

Interference heating due to shock wave impingement on laminar boundary layers.

Laminar interference heating correlations have been developed based on recent experimental data obtained with wedge/flat plate models for wide ranges of Reynolds number and shock strength. Two correlation techniques were developed using the Eckert reference method. The peak interference Stanton number was first correlated with shock strength, Reynolds number, and Prandtl number based on flow conditions upstream of the interference region. The second approach was made by correlating peak interference Stanton number with only Reynolds number and Prandtl number based on downstream flow conditions. The laminar boundary layer remains laminar when both Reynolds number and shock strength are low but becomes transitional or turbulent when Reynolds number or/and shock strength are increased.

Hung, F. T.↗
NASA NTRS1977

Diffraction of a shock wave by a compression corner. II - Single Mach reflection

The two-dimensional time-dependent Euler equations which govern the flow field resulting from the interaction of a planar shock with a compression corner are solved for initial conditions which result in single Mach reflection of the incident planar shock. The Euler equations are first transformed to include the self-similarity of the flow field. A second transformation is employed to normalize the distances between the ramp and the reflected shock and between the wall and the Mach stem. The resulting equations in strong conservation-law form are solved using a second-order discontinuity-fitting finite-difference approach. The results are compared with experimental interferograms and existing first-order shock-capturing numerical solutions.

Shankar, V.↗
NASA NTRS1990

Comparison of shock wave rise time prediction techniques

Two techniques for predicting sonic-boom rise times are described, and the values of the rise time (defined as the time required for the shock front to go from 10 percent to 90 percent of the maximum overpressure) predicted when using the two approaches are compared. The results are applied to the prediction of the rise time from two hypothetical supersonic transport aircraft. In both cases, the rise time predicted at the ground was much greater than that predicted for a conventional SST.

Raspet, Richard↗
NASA NTRS1990

Structure and evolution of small-amplitude intermediate shock waves

A simplified set of equations is derived that approximates the magnetohydrodynamic (MHD) Navier-Stokes equations for weakly nonlinear disturbances whose speeds are close to the MHD intermediate speed. Its shock structure solutions are then examined. The fast and slow shock solutions are uniquely specified by their Rankine-Hugoniot relations. However, the intermediate shock solutions, which are not unique, are characterized by the integral through the shock of the noncoplanar component of the magnetic field. For situations in which this integral is conserved, the Riemann problem is well defined and predicts the evolution of intermediate shocks. This analysis is substantiated by numerical computations.

Kennel, C. F.↗
NASA NTRS1990

Acceleration of energetic particles by shock waves from large solar flares

Time-intensity profiles for solar proton events are analyzed with respect to the acceleration of energetic particles from the ambient solar wind by an interplanetary shock. The time-intensity profiles are derived from the Helio 1 particle data base for 1974 to 1984. It is found that many of the large proton events display a prolonged plateau. The Fe/O abundances in the large proton events are studied. It is observed that initially the Fe/O ratio is high (about 1.0), but after proton intensities increase to about 100 protons/sq cm sr sec MeV the Fe/O ratio decreases to a value near 0.1.

Reames, Donald V.↗
NASA NTRS1974

Post-shock spikes: A new feature of proton and alpha enhancements associated with an interplanetary shock wave

Abrupt and prolonged enhancements in the intensities of 100 to approximately 2000 keV nucleon protons and alpha particles observed in interplanetary space are interpreted as particle populations confined between an interplanetary shock front and a magnetic field discontinuity. Prominent intensity spikes observed only below approximately 400 keV per charge for both protons and alpha particles several hours behind the shock front suggest that some fraction of the confined particles is accelerated by an energy per charge dependent process.

Gloeckler, G.↗