Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK LAYER”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Boundary-layer electron profiles for entry of a blunts slender body at high altitude

New calculations of boundary-layer electron concentration profiles for entry of a blunt-nosed slender body into the earth's atmosphere are compared with previous calculations in which ambipolar diffusion was neglected. The old and new results agree in those flight regimes where ambipolar diffusion is unimportant, but large differences are noted in both peak electron concentration and profile shape at the higher altitudes, where diffusion effects are greatest. The new results are also compared with flight-measured profiles and with calculated profiles for a viscous-shock-layer theory which was recently reported in the literature. The boundary-layer results and the data agree in most respects. Differences which occur between predicted results and the data in the outer parts of the profile are discussed in terms of the effects of aerodynamic heating of the probes.

Evans, J. S.↗
NASA NTRS1986

Computation of three-dimensional flow about aerobrake configurations

Ellipsoid, cone and cylinder aerobrake configurations are analyzed to provide comparison data between experimental and model predictions. An analytical model was devised to account for the shock layer ahead of the body and in the near-wake region in terms of the Navier-Stokes equations expressed in conformal polar and azimuthal-angle coordinates. Using polar coordinates simplified the equations by mapping the body onto a sphere, a procedure which also reduced the magnitude of the discretization errors. The equations are then solved using an alternating direction implicit (ADI) factorization technique. Computations were carried out for Mach 3-10 at various grid resolutions and compared with available wind tunnel data. The model generated pressure distributions, heat transfer coefficients and velocity profile data that agreed relatively well with experimental data at a reduced computational cost. Further work is necessary to identify the location of shocks and to model flows about asymmetric configurations.

Li, C. P.↗
NASA NTRS2013

Mars Science Laboratory Heatshield Aerothermodynamics: Design and Reconstruction

The Mars Science Laboratory heatshield was designed to withstand a fully turbulent heat pulse based on test results and computational analysis on a pre-flight design trajectory. Instrumentation on the flight heatshield measured in-depth temperatures in the thermal protection system. The data indicate that boundary layer transition occurred at 5 of 7 thermocouple locations prior to peak heating. Data oscillations at 3 pressure measurement locations may also indicate transition. This paper presents the heatshield temperature and pressure data, possible explanations for the timing of boundary layer transition, and a qualitative comparison of reconstructed and computational heating on the as-flown trajectory. Boundary layer Reynolds numbers that are typically used to predict transition are compared to observed transition at various heatshield locations. A uniform smooth-wall transition Reynolds number does not explain the timing of boundary layer transition observed during flight. A roughness-based Reynolds number supports the possibility of transition due to discrete or distributed roughness elements on the heatshield. However, the distributed roughness height would have needed to be larger than the pre-flight assumption. The instrumentation confirmed the predicted location of maximum turbulent heat flux near the leeside shoulder. The reconstructed heat flux at that location is bounded by smooth-wall turbulent calculations on the reconstructed trajectory, indicating that augmentation due to surface roughness probably did not occur. Turbulent heating on the downstream side of the heatshield nose exceeded smooth-wall computations, indicating that roughness may have augmented heating. The stagnation region also experienced heating that exceeded computational levels, but shock layer radiation does not fully explain the differences.

Edquist, Karl T.↗
NASA NTRS2001

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗
NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗
NASA NTRS1982

An experimental investigation and a numerical prediction of a transonic normal shock wave/turbulent boundary layer interaction

Detailed pitot, static and wall pressure measurements have been obtained for a transonic normal shock wave/turbulent boundary layer interaction at free-stream Mach numbers of 1.28, 1.37 and 1.48 and at a constant unit Reynolds number of 4.92 x 10 to the 6th per meter in an axisymmetric, internal flow. Measurements have also been obtained at a unit Reynolds number of 9.84 x 10 to the 6th per meter at a free-stream Mach number of 1.29. The interaction depends very strongly on the Mach number. The effect of Reynolds number under study is small. Flow confinement due to the wind-tunnel wall boundary layer produces a weaker interaction and a much larger supersonic tongue than is observed for planar flows. Comparisons are made with solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model (1980). The computations are in agreement with the experimental results.

Om, D.↗
NASA NTRS1989

Experimental and numerical investigation of an oblique shock wave/turbulent boundary layer interaction with continuous suction

An numerical and experimental investigation has been conducted into the interaction of an incident oblique shock wave with a turbulent boundary layer, for the cases of a rough plate and a porous plate with suction, at a nominal Mach number of 2.5 and flow deflection angles of 0, 4, 6, and 8 deg. Attention is given to the pitot pressure profiles, wall static pressures, and porous plate local bleed distributions measured for the two plates. Suction is found to increase the strength of the incident shock required to separate the boundary layer; for all shock strengths tested, separation is completely eliminated.

Benhachmi, Driss↗
NASA NTRS1989

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗
NASA NTRS1992

Research activity at the shock tube facility at NASA Ames

The real gas phenomena dominate the relaxation process occurring in the flow around hypersonic vehicles. The air flow around these vehicles undergoes vibrational excitation, chemical dissociation, and ionization. These chemical and kinetic phenomena absorb energy, change compressibility, cause temperature to fall, and density to rise. In high-altitude, low density environments, the characteristic thicknesses of the shock layers can be smaller than the relaxation distances required for the gas to attain chemical and thermodynamic equilibrium. To determine the effects of chemical nonequilibrium over a realistic hypersonic vehicle, it would be desirable to conduct an experiment in which all aspects of fluid flow are simulated. Such an experiment is extremely difficult to setup. The only practical alternative is to develop a theoretical model of the phenomena and to compute the flow around the vehicle including the chemical nonequilibrium, and compare the results with the experiments conducted in the facilities under conditions where only a portion of the flow phenomena is simulated. Three types of experimental data are needed to assist the aerospace community in this model development process: (1) data which will enhance our phenomenological understanding of the relaxation process, (2) data on rate reactions for the relevant reactions, and (3) data on bulk properties, such as spectral radiation emitted by the gas, for a given set of aerodynamic conditions. NASA Ames is in a process of collecting such data by simulating the required aerothermochemical conditions in an electric arc driven shock tube.

Sharma, Surendra P.↗
NASA NTRS1989

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗
NASA NTRS2000

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗
DOE OSTI2021

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1963

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗
NASA NTRS1965

Investigation to define the propagation characteristics of a finite amplitude acoustic pressure wave, phase 1 final report, 29 jun. 1964 - 29 jul. 1965

The contribution of high entropy production regions to the generation and propagation characteristics of a finite amplitude pressure is considered. Preliminary analysis indicates that, for nozzles where pressure rations are above critical, the predominant contribution may come from the shock layer formation in the exhaust region. Temperature effects indicate high dependence of the forcing function upon the initial temperature of the media.

WAVE PROPAGATION↗