Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 469 records · Page 26

Accounting for the Influence of Decorrelation in Microphone Phased Array Deconvolution Methods

Microphone phased arrays are a common tool for use in aeroacoustic wind tunnel testing. The analysis of acquired array data is known to suffer from decorrelation effects, where the coherence of an acoustic wave measured by a pair of microphones is degraded as the wave passes through a turbulent free shear layer or boundary layer. This paper describes, in detail, how to mitigate the influence of decorrelation effects when processing array data with deconvolution methods. This is done using the DAMAS algorithm as an example, applied to recent airframe noise test data acquired in the NASA Langley 14- by 22-Foot Subsonic Tunnel. Two ways of handling the turbulent propagation modeling, both assuming plane wave propagation, are described. Results show that while turbulence model fit parameters may differ, both methods output extremely similar deconvolution results. Further improvements likely require more accurate mean shear layer data prior to developing more involved turbulence models.

microphone array↗

Afterbody Heating Characteristics of a Proposed Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 deg sideslip, and normal-shock density, ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75% of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three-dimensional separated wake flows.

Horvath, Thomas J.↗

DIAMOND - an Architecture for Persistent Space Platforms

We define DIAMOND, an architecture for generalized, persistent platforms in space, providing support infrastructure for multiple payloads with diverse missions, and able to accept the integration of payloads delivered to the platform after it is established in space. Such platforms would create a new space ecosystem with opportunities for much larger platforms than currently exist, for adaptation to changing technological and market conditions over long periods of time, and for tenants of the platform to launch only their unique payloads without having to carry all the equipment needed to support them. The DIAMOND architecture uses a novel “Architectural Shearing Layer” approach combining long-lived but replaceable structural and support infrastructure components with shorter-lived, higher-value service and tenant facilities. The structural aspect of the proposed architecture is a regular lattice constructed of tetrahedral cells, constructed of only a few types of simple structural components, patterned according to the cubic lattice of natural diamond crystals. The cubic lattice is arranged on larger scales to be fractal sponge, providing a disproportionately large number of surface attachment points while maintaining overall low mass of the entire structure. Generalized architectural shearing layers are also provided that permit easy separation between structural, power, thermal, propulsion/attitude control, pointing/vibration control, computing/data storage, and communication, and various operational aspects of the architecture. The combined result ensures longevity by enabling space platforms to undergo many cycles of change with little peril and at low cost, changing any part of the platform at different times and rates while leaving the remaining parts of the platform undisturbed.

Breidenthal, Julian↗

Separated and Recovering Turbulent Boundary Layer Flow Behind a Backward Facing Step For Different Reynolds Numbers

Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Jovic, Srba↗

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗

An experimental study of natural and forced modes in an axisymmetric jet

The experiment consisted of two parts. The first part was an effort to study naturally occurring instability modes in the axisymmetric shear layer of a high Reynolds number turbulent jet. Untripped (transitional) and tripped (turbulent) nozzle exit conditions, both with 0.1 percent core turbulence, were studied. For the turbulent nozzle exit boundary layer case, the core turbulence, was varied systematically from 0.1 to 5 percent of the jet exit velocity. The region up to the end of the potential core was dominated by the axisymmetric mode. The azimuthal modes grew rapidly but dominated only downstream of the potential core region. For the jet excited by natural disturbances, the energy content of the higher order modes (m is less than 1) was significantly lower than that of the axisymmetric and m = +/-1 modes. The initial boundary layer had a profound effect on the natural jet evolution and its excitability. The shorter potential core allowed the jet in the transitional case to support helical disturbances closer the nozzle exit than the turbulent case. The natural jet evolution was found to remain unaffected for varying initial core turbulence over the range of 0.15 to 5 percent of the jet velocity. Target modes for efficient excitation of the jet were determined from the results of the naturally occurring jet instability mode experiments. The second part of this work describes an effort to control the axisymmetric shear layer by artificially exciting target modes. Under optimum conditions, two-frequency excitation is indeed more effective than single frequency plane wave excitation as far as jet mixing is concerned. At high amplitudes of fundamental and subharmonic forcing, the subharmonic augmentation and the axial location of the peak were independent of the initial phase difference. Two-frequency excitation also has its limitations, since axisymmetric waves are damped beyond the potential core. Higher spreading rates are obtained when multi-modal forcing is applied.

Raman, Ganesh↗

Axisymmetric computational fluid dynamics analysis of Saturn V/S1-C/F1 nozzle and plume

An axisymmetric single engine Computational Fluid Dynamics calculation of the Saturn V/S 1-C vehicle base region and F1 engine plume is described. There were two objectives of this work, the first was to calculate an axisymmetric approximation of the nozzle, plume and base region flow fields of S1-C/F1, relate/scale this to flight data and apply this scaling factor to a NLS/STME axisymmetric calculations from a parallel effort. The second was to assess the differences in F1 and STME plume shear layer development and concentration of combustible gases. This second piece of information was to be input/supporting data for assumptions made in NLS2 base temperature scaling methodology from which the vehicle base thermal environments were being generated. The F1 calculations started at the main combustion chamber faceplate and incorporated the turbine exhaust dump/nozzle film coolant. The plume and base region calculations were made for ten thousand feet and 57 thousand feet altitude at vehicle flight velocity and in stagnant freestream. FDNS was implemented with a 14 species, 28 reaction finite rate chemistry model plus a soot burning model for the RP-1/LOX chemistry. Nozzle and plume flow fields are shown, the plume shear layer constituents are compared to a STME plume. Conclusions are made about the validity and status of the analysis and NLS2 vehicle base thermal environment definition methodology.

Ruf, Joseph H.↗

Enhanced mixing of a rectangular supersonic jet by natural and induced screech

The influence of shear layer excitation on the mixing of supersonic rectangular jets was studied experimentally. Two methods of excitation were used to control the jet mixing. The first used the natural screech of an underexpanded supersonic jet from a converging nozzle. The level of the screech excitation was controlled by the use of a pair of baffles located to block the acoustic feedback path between the downstream shock structure and the nozzle lip. A screech level variation of over 30 decibels was achieved and the mixing was completely determined by the level of screech attained at the nozzle lip. The second form of self-excitation used the induced screech caused by obstacles or paddles located in the shear layers on either long side of the rectangular jet. With sufficient immersion of the paddles intense jet mixing occurred and large flapping wave motion was observed using a strobed focused Schlieren system. Each paddle was instrumented with a total pressure tap and strain gages to determine the pressure and drag force on the square cross-section paddle. Considerable drag was observed in this initial exploratory study. Future studies using alternate paddle geometries will be conducted to maximize jet mixing with minimum drag.

Rice, Edward J.↗

Transition and Turbulence Modeling for Blunt-Body Wake Flows

Aerobraking has been proposed as an efficient means of decelerating spacecraft for planetary missions. Most current aerobrake designs feature a blunt forebody shielding the payload from the intense heat generated during atmospheric entry. Although this forebody will absorb the largest portion of the heat pulse, accurate prediction of heating in the near wake is of great importance, since large local heating values can occur at points of shear-layer impingement. In order to address the various issues associated with these blunt-body wake flowfields, the Advisory Group for Aerospace Research and Development (AGARD) formed Working Group 18 in 1992. One of the objectives of this activity was to examine real-gas effects in high-speed flow fields around a 70 deg. blunted cone. To date, many researchers have conducted experiments using this geometry in various facilities, such as the Large Energy National Shock (LENS) tunnel at Cubric/Calspan and the HEG shock tunnel at DLR-Goettingen. Several computational studies have also been conducted in concert with these tests. Many of the experimental results have indicated the possible presence of a transitional shear layer through a large increase in heat transfer downstream of the reattachment point. The presence of transition could in fact lead to much higher peak heating than if the separated flow is entirely laminar or turbulent. In the shock-tunnel tests, however, it is difficult to separate such viscous-flow phenomena from real-gas effects. In order to help make this distinction, Horvath et al. recently conducted a set of experiments in the NASA Langley 20-Inch Mach 6 Tunnel, and compared the results to laminar Navier-Stokes calculations. They found heat-transfer distributions similar to those obtained in the high-enthalpy facilities, with the measured peak heating along the sting support markedly greater than that predicted by the laminar computations. These trends point to the need to find transitional and turbulent computational solutions for these flowfields.

Nance, Robert P.↗

Stable Lobed Mixer With Combustion Demonstrated and Measured

The NASA Lewis Research Center collaborated with the Massachusetts Institute of Technology (MIT) on an experiment to study the use of lobed mixers to improve the fuel-air mixing process and increase combustion intensity in combustors with minimal pressure loss. This experiment is the first known stable combusting flow studied for this device, and the data show a much faster and much more uniform combustion process than for flat-plate mixers. Several potential benefits may be realized from this study in future combustors, including a reduction in NO_x emissions because of the more uniform temperature distribution. The experiment was done in Lewis' Planar Reacting Shear Layer facility, which was adapted to accept a lobed mixer in addition to the original planar tip. A graduate student at MIT provided the mixer design concept, and Lewis provided the engineering, operations, and research expertise. The experiment used hydrogen-nitrogen mixtures to react with vitiated hot air at 920 K. A flow speed of about 120 m/sec and a speed ratio of 0.5 were used. Flow diagnostics consisted of traversing fine-wire thermocouples and pitot probes for flow mapping. Supplementary fluorescence images were taken with a charged coupled device (CCD) camera to show the location and temporal behavior of the reaction zone. The data showed that the lobed mixer consumed the reactants between 3 to 10 times faster than a corresponding planar shear layer. The figure shows the dramatic difference in the measured temperature distribution with and without the lobed mixer. The increased mixing rate was due to a larger interfacial area as well as to the secondary flow from the streamwise vortices off the tips of the lobes. In addition, the fluorescence images showed that the lobes acted as flame stabilizers.

Source record↗

Green's Function Applicable to Turbofan Exhaust Noise in Jets with an External Center-Body

The problem of propagation of sound across the shear layer in a turbofan jet exhaust with an external center-body is discussed. The wave equation of interest is compressible Rayleigh equation. Two forms of the equation are considered, and the Green's function solutions subject to appropriate surface conditions on the center-body and flight condition in the ambient are presented. Directivity studied in a heated exhaust at temperature ratio of 2.0 and Mach number 0.90 indicate that a rigid center-body tends to increase the sound propagation at forward angles relative to an exhaust without a center-body, while application of suitable surface liner may significantly reduce this enhancement. A general form of the far-field solution to the propagation equation in a parallel flow under supersonic conditions and in the neighborhood of the shear layer singularity is also discussed.

Noise↗

Predicting the stability of a compressible periodic parallel jet flow

It is known that mixing enhancement in compressible free shear layer flows with high convective Mach numbers is difficult. One design strategy to get around this is to use multiple nozzles. Extrapolating this design concept in a one dimensional manner, one arrives at an array of parallel rectangular nozzles where the smaller dimension is omega and the longer dimension, b, is taken to be infinite. In this paper, the feasibility of predicting the stability of this type of compressible periodic parallel jet flow is discussed. The problem is treated using Floquet-Bloch theory. Numerical solutions to this eigenvalue problem are presented. For the case presented, the interjet spacing, s, was selected so that s/omega =2.23. Typical plots of the eigenvalue and stability curves are presented. Results obtained for a range of convective Mach numbers from 3 to 5 show growth rates omega(sub i)=kc(sub i)/2 range from 0.25 to 0.29. These results indicate that coherent two-dimensional structures can occur without difficulty in multiple parallel periodic jet nozzles and that shear layer mixing should occur with this type of nozzle design.

Miles, Jeffrey H.↗

Studies of shock/shock interaction on smooth and transpiration-cooled hemispherical nosetips in hypersonic flow

A program of experimental research and analysis was conducted to examine the heat transfer and pressure distributions in regions of shock/shock interaction over smooth and transpiration-cooled hemispherical noseshapes. The objective of this investigation was to determine whether the large heat transfer generated in regions of shock/shock interaction can be reduced by transpiration cooling. The experimental program was conducted at Mach numbers of 12 to 16 in the Calspan 48-Inch Shock Tunnel. Type 3 and type 4 interaction regions were generated for a range of freestream unit Reynolds numbers to provide shear layer Reynolds numbers from 10 exp 4 to 10 exp 6 to enable laminar and turbulent interaction regions to be studied. Shock/shock interactions were investigated on a smooth hemispherical nosetip and a similar transpiration-cooled nosetip, with the latter configuration being examined for a range of surface blowing rates up to one-third of the freestream mass flux. While the heat transfer measurements on the smooth hemisphere without shock/shock interaction were in good agreement with Fay-Riddell predictions, those on the transpiration-cooled nosetip indicated that its intrinsic roughness caused heating-enhancement factors of over 1.5. In the shock/shock interaction studies on the smooth nosetip, detailed heat transfer and pressure measurements were obtained to map the variation of the distributions with shock-impingement position for a range of type 3 and type 4 interactions. Such sets of measurements were obtained for a range of unit Reynolds numbers and Mach numbers to obtain both laminar and turbulent interactions. The measurements indicated that shear layer transition has a significant influence on the heating rates for the type 4 interaction as well as the anticipated large effects on type 3 interaction heating. In the absence of blowing, the peak heating in the type 3 and type 4 interaction regions, over the transpiration-cooled model, did not appear to be influenced by the model's rough surface characteristics. The studies of the effects of the transpiration cooling on type 3 and type 4 shock/shock interaction regions demonstrated that large surface blowing rates had significant effect on the structure of the flowfield, enlarging the shock layer and moving the region of peak-heating interaction around the body.

Holden, Michael S.↗

Turbulance boundary conditions for shear flow analysis, using the DTNS flow solver

The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Mizukami, M.↗

Axisymmetric Flow Properties for Magnetic Elements of Differing Strength

Aspects of the structure and dynamics of the flows in the Sun's surface shear layer remain uncertain and yet are critically important for understanding the observed magnetic behavior. In our previous studies of the axisymmetric transport of magnetic elements we found systematic changes in both the differential rotation and the meridional flow over the course of Solar Cycle 23. Here we examine how those flows depend upon the strength (and presumably anchoring depth) of the magnetic elements. Line of sight magnetograms obtained by the HMI instrument aboard SDO over the course of Carrington Rotation 2097 were mapped to heliographic coordinates and averaged over 12 minutes to remove the 5-min oscillations. Data masks were constructed based on the field strength of each mapped pixel to isolate magnetic elements of differing field strength. We used Local Correlation Tracking of the unmasked data (separated in time by 1- to 8-hours) to determine the longitudinal and latitudinal motions of the magnetic elements. We then calculated average flow velocities as functions of latitude and longitude from the central meridian for approx 600 image pairs over the 27-day rotation. Variations with longitude indicate and characterize systematic errors in the flow measurements associated with changes in the signal from disk center to limb. Removing these systematic errors reveals changes in the axisymmetric flow properties that reflect changes in flow properties with depth in the surface shear layer.

Rightmire-Upton, Lisa↗

Evaluation of Aircraft Platforms for SOFIA by Computational Fluid Dynamics

The selection of an airborne platform for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is based not only on economic cost, but technical criteria, as well. Technical issues include aircraft fatigue, resonant characteristics of the cavity-port shear layer, aircraft stability, the drag penalty of the open telescope bay, and telescope performance. Recently, two versions of the Boeing 747 aircraft, viz., the -SP and -200 configurations, were evaluated by computational fluid dynamics (CFD) for their suitability as SOFIA platforms. In each configuration the telescope was mounted behind the wings in an open bay with nearly circular aperture. The geometry of the cavity, cavity aperture, and telescope was identical in both platforms. The aperture was located on the port side of the aircraft and the elevation angle of the telescope, measured with respect to the vertical axis, was 500. The unsteady, viscous, three-dimensional, aerodynamic and acoustic flow fields in the vicinity of SOFIA were simulated by an implicit, finite-difference Navier-Stokes flow solver (OVERFLOW) on a Chimera, overset grid system. The computational domain was discretized by structured grids. Computations were performed at wind-tunnel and flight Reynolds numbers corresponding to one free-stream flow condition (M = 0.85, angle of attack alpha = 2.50, and sideslip angle beta = 0 degrees). The computational domains consisted of twenty-nine(29) overset grids in the wind-tunnel simulations and forty-five(45) grids in the simulations run at cruise flight conditions. The maximum number of grid points in the simulations was approximately 4 x 10(exp 6). Issues considered in the evaluation study included analysis of the unsteady flow field in the cavity, the influence of the cavity on the flow across empennage surfaces, the drag penalty caused by the open telescope bay, and the noise radiating from cavity surfaces and the cavity-port shear layer. Wind-tunnel data were also available to compare to the CFD results; the data permitted an assessment of CFD as a design tool for the SOFIA program.

Klotz, S. P.↗

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗