Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK WAVE INTERACTION”

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 523 records · Page 29

Parabolized Navier-Stokes analysis of ducted turbulent mixing problems with finite-rate chemistry

A hybrid implicit/explicit approach for analyzing 2D ducted turbulent mixing problems with finite-rate chemistry is presented. The approach combines a fully-implicit parabolic mixing algorithm, an explicit viscous-characteristic based shock-capturing algorithm, and linearized implicit chemical kinetic algorithm. The resultant model provides spatial marching solutions of the parabolized Navier-Stokes (PNS) equations for supersonic combustion and viscous nozzle flowfields. Specialized procedures are incorporated to deal with the near wall sublayer and with small central pockets of subsonic flow. A parabolic option is available which can be utilized in the direct or inverse mode for design applications. The ability of the model to treat shock waves and wave/mixing layer interactions is assessed via comparisons of predictions with those of well established explicit shock-capturing Euler (SCIPPY) and PNS (SCIPVIS) models. Applications to a variety of supersonic combustion flowfield problems involving tangential or moderately inclined fuel injection, and, to a viscous nozzle flow problem, are presented. This paper serves to exhibit overall capabilities of the model developed and no comparisons with supersonic combustion data are presented. Such comparisons serve mainly to validate the modeling of the turbulence and turbulence/chemical interactions and will be the subject of a future paper.

Sinha, N.↗

Time dependent heat transfer rates in high Reynolds number hypersonic flowfields

Time dependent heat transfer rates have been calculated from time dependent temperature measurements in the vicinity of shock-wave boundary-layer interactions due to conical compression ramps on an axisymmetric body. The basic model is a cylindrical body with a 10 degree conical nose. Four conical ramps, 20, 25, 30, and 35 degrees serve as shock wave generators. Flowfield surveys have been made in the vicinity of the conical ramp vertex, the separation point, and the reattachment point. A significant effort was made to characterize the natural frequencies and relative powers of the resulting fluctuations in heat transfer rates. This research effort, sponsored jointly by NASA and the Air Force, was conducted in the Air Force Flight Dynamics Directorate High Reynolds Facility. The nominal freestream Mach number was 6, and the freestream Reynolds numbers ranged from 2.2 million/ft to 30.0 million/ft. Experimental results quantify temperature response and the resulting heat transfer rates as a function of ramp angle and Reynolds number. The temperature response within the flowfield appears to be steady-state for all compression ramp angles and all Reynolds numbers, and hence, the heat transfer rates appear to be steady-state.

Flanagan, Michael J.↗

An experimental study of fluctuating pressure loads beneath swept shock/boundary-layer interactions

A database is established on the fluctuating pressure loads produced on aerodynamic surfaces beneath 3-D shock wave/boundary layer interactions. Such loads constitute a fundamental problem of critical concern to future supersonic and hypersonic flight vehicles. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins. Fin angles from 5 to 25 deg at freestream Mach numbers between 2.5 and 4 produce a variety of interaction strengths from weak to very strong. Miniature Kulite pressure transducers mounted in the flat plate were used to measure interaction-induced wall pressure fluctuations. These data will be correlated with proposed new optical data on the fluctuations of the interaction structure, especially that of the lambda-shock system and its associated high-speed jet impingement.

Settles, Gary S.↗

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

Turbulent wake/shock interaction at Mach 6

The paper presents measurements of mean pitot, static pressure and total temperature profiles for a two-dimensional turbulent wake in the vicinity and downstream of intersection with a shock/expansion wave system. The wake centerline Mach number, velocity, and total temperature distributions indicated that the shock increased turbulent mixing, while the expansion field that followed the shock decreased mixing. The overall effect of the shock/expansion wave interaction was to increase mixing for shock angles 15 deg and less for the shock generator geometry used in the present study. Wake growth rates determined from the total temperature profiles confirmed the increased turbulent mixing effect of the shock. These data could be used to validate nonquilibrium turbulence modeling and numerical solution of the time-averaged Navier-Stokes equations.

Walsh, M. J.↗

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗

Shock interacting with a random array of stationary particles underwater

Accurately predicting the force experienced by particles within a random bed during shock traversal is a challenging problem. In this paper, particle-resolved inviscid simulations of shock propagation over randomly distributed beds are performed, with the goal of quantifying the force on individual particles. Water was considered the continuum material, and then results were compared with previous research completed in air and the differences resulting from these different mediums were examined. Simulations were conducted for four different combinations of incident Mach number and bed volume fraction with the particles remaining stationary. Additionally, time-resolved streamwise and transverse force coefficients were calculated for each particle in the bed. It was observed that while the average force coefficient was similar to that of an isolated particle, the force on individual particles was substantially different. An important observation was that the inviscid drag force on some particles remained consistently positive, but others consistently negative even long after the passage of the shock. This persistent force on the particle was attributed to the potential flow fields that result from neighbors' specific locations upstream and downstream. Similar behavior was observed in the transverse force, with particles experiencing sustained force in the same lateral direction. These persistent forces are inviscid quasisteady contributions which with the viscous counterpart will play an essential role in the particle's long-term dispersion if allowed to move in response to the force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interactive phenomena in supersonic jet mixing problems. II - Numerical studies

Applications of the numerical techniques discussed in Part I of this article are presented. Fundamental wave/shear-layer interaction calculations are performed exhibiting the waves produced by high-speed turbulent mixing processes and the interactions of an expansion fan and shock wave with a developing shear layer. Calculations for an underexpanded supersonic jet into a supersonic external stream at different velocity ratios are presented which delineate betweeen weakly interactive and strongly interactive processes in the jet near field. Calculations for the multiple-cell shock structure of an underexpanded jet into still air exhibit the damping of wave intensities by the turbulence and compare quite favorably with available data.

Dash, S. M.↗

A validation study of four Navier-Stokes codes for high-speed flows

A code validation study has been conducted for four different codes for solving the compressible Navier-Stokes equations. Computations for a series of nominally two-dimensional high-speed laminar separated flows were compared with detailed experimental shock-tunnel results. The shock wave-boundary layer interactions considered were induced by a compression ramp in one case and by an externally-generated incident shock in the second case. In general, good agreement was reached between the grid-refined calculations and experiment for the incipient- and small-separation conditions. For the most highly separated flow, three-dimensional calculations which included the finite-span effects of the experiment were required in order to obtain agreement with the data. The finite-span effects were important in determining the extent of separation as well as the time required to establish the steady-flow interaction. The results presented provide a resolution of discrepancies with the experimental data encountered in several recent computational studies.

Rudy, David H.↗

Wall pressure fluctuations beneath swept shock wave/boundary layer interactions

An experimental research program providing basic knowledge and establishing a database on the fluctuating pressure loads produced on aerodynamic surfaces beneath 3D shock wave/boundary layer interactions is presented. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins at angle of attack. Measurements are made for the first time in the aft areas of these interactions, showing fluctuating pressure levels as high as 160 dB.

Garg, S.↗

Three-Dimensional Navier-Stokes Simulations with Two-Equation Turbulence Models of Intersecting Shock-Waves/Turbulent Boundary Layer at Mach 8.3

An investigation of the numerical simulation with two-equation turbulence models of a three-dimensional hypersonic intersecting (SWTBL) shock-wave/turbulent boundary layer interaction flow is presented. The flows are solved with an efficient implicit upwind flux-difference split Reynolds-averaged Navier-Stokes code. Numerical results are compared with experimental data for a flow at Mach 8.28 and Reynolds number 5.3x10(exp 6) with crossing shock-waves and expansion fans generated by two lateral 15 fins located on top of a cold-wall plate. This experiment belongs to the hypersonic database for modeling validation. Simulations show the development of two primary counter-rotating cross-flow vortices and secondary turbulent structures under the main vortices and in each corner singularity inside the turbulent boundary layer. A significant loss of total pressure is produced by the complex interaction between the main vortices and the uplifted jet stream of the boundary layer. The overall agreement between computational and experimental data is generally good. The turbulence modeling corrections show improvements in the predictions of surface heat transfer distribution and an increase in the strength of the cross-flow vortices. Accurate predictions of the outflow flowfield is found to require accurate modeling of the laminar/turbulent boundary layers on the fin walls.

Bardina, J. E.↗

Relatively stable, large-amplitude Alfvenic waves seen at 2.5 and 5.0 AU

Pioneer 11 and 10 observations of the wave structure seen in a corotating interaction region at 2.5 AU on day 284 of 1973 and 8 days later at 5 AU reveal large-amplitude Alfvenic structures with many detailed correlations seen between their features at the two radial distances. Hodogram analysis suggests the dominance of near plane polarized, transverse Alfvenic mode fluctuations with periods between 2 min and one hour or more. Some wave evolution close to the Corotating Interaction Region (CIR) shock is noticed, but waves towards the center of the compression seem to propagate with little damping between the spacecraft observation positions.

Mavromichalaki, H.↗

An experimental study of the sources of fluctuating pressure loads beneath swept shock/boundary-layer interactions

An experimental research program providing basic knowledge and establishing a database on the fluctuating pressure loads produced on aerodynamic surfaces beneath three dimensional shock wave/boundary layer interactions is described. Such loads constitute a fundamental problem of critical concern to future supersonic and hypersonic flight vehicles. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins at angle of attack. Fin angles from 10 to 20 deg at freestream Mach numbers of 3 and 4 produce a variety of interaction strengths from weak to very strong. Miniature Kulite pressure transducers flush-mounted in the flat plate are used to measure interaction-induced wall pressure fluctuations. The distributions of properties of the pressure fluctuations, such as their ring levels, amplitude distributions, and power spectra, are also determined. Measurements were made for the first time in the aft regions of these interactions, revealing fluctuating pressure levels as high as 160 dB. These fluctuations are dominated by low frequency (0-5 kHz) signals. The maximum ring levels in the interactions show an increasing trend with increasing interaction strength. On the other hand, the maximum ring levels in the forward portion of the interactions decrease linearly with increasing interaction sweep back. These ring pressure distributions and spectra are correlated with the features of the interaction flowfield. The unsteadiness of the off-surface flowfield is studied using a new, non-intrusive technique based on the shadow graph method. The results indicate that the entire lambda-shock structure generated by the interaction undergoes relatively low-frequency oscillations. Some regions where particularly strong fluctuations are generated were identified. Fluctuating pressure measurements are also made along the line of symmetry of an axisymmetric jet impinging upon a flat plate at an angle. This flow was chosen as a simple analog to the impinging jet region found in the rear portion of the shock wave/boundary layer interactions under study. It is found that a sharp peak in ring pressure level exists at or near the mean stagnation point. It is suggested that the phenomena responsible for this peak may be active in the swept interactions as well, and may cause the extremely high fluctuating pressures observed in the impinging jet region in the present experimental program.

Settles, G. S.↗

Computational considerations for the simulation of shock-induced sound

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order-accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order-accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks in a quasi-one-dimensional nozzle flow. The design order of accuracy is achieved in the smooth regions of a steady-state test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order-accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustics simulations of flows with shocks that are raised by these results are discussed.

Casper, Jay↗

Understanding the plasticity contributions during laser-shock loading and spall failure of Cu microstructures at the atomic scales

A hybrid atomic-scale and continuum modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline (nc) Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. Here, the simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.

36 MATERIALS SCIENCE↗

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗

High speed interference heating loads and pressure distributions resulting from elevon deflections

Effects of elevon-induced three-dimensional shock-wave turbulent boundary-layer interactions on hypersonic aircraft surfaces are analyzed. Detailed surface pressure and heating rate distributions obtained on wing-elevon-fuselage models representative of aft portions of hypersonic aircraft are compared with analytical and experimental results from other sources. Examples are presented that may be used to evaluate the adequacy of current theoretical methods for estimating the effects of three-dimensional shock-wave turbulent boundary-layer interactions on hypersonic aircraft surfaces.

Johnson, C. B.↗

Estimation and comparison of quasilinear electron heating in the shock foot at Jupiter and earth

A simple quasilinear model is developed to estimate the electron heating that occurs in the foot of a supercritical, quasiperpendicular shock through interactions with electrostatic waves generated by reflected ions. At earth the increase in electron thermal energy calculated using the measured wave amplitudes is negligible, while at Jupiter it is comparable with the observed temperature gain across the shock. The anisotropic quasilinear heating should destabilize whistler mode waves in the foot. These have been detected by the plasma wave instrument on Voyager with amplitudes sufficient to isotropize the electron distribution.

Moses, S. L.↗