Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock Physics”

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 865 records · Page 48

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗

Physical properties of glasses exposed to Earth-facing and trailing-side environments on LDEF

The exposure of 108 glass samples and 12 glass-ceramic samples to Earth-orbit environments permitted measurements which establish the effects of each environment. Examination of five glass types and one glass ceramic located on both the Earth-facing side and the trailing edge revealed no reduction in strength within experimental limits. Strength measurements subjected less than 5 percent of the sample surface area to stresses above 90 percent of the glass's failure strength. Seven micrometeorite or space debris impacts occurred on trailing edge samples. One of those impacts occurred in a location which was subjected to 50 percent of the applied stress at failure. Micrometeorite or space debris impacts were not observed on Earth-facing samples. The physical shape and structure of the impact sites were carefully examined using stereographic scanning electron microscopy. These impacts induce a stress concentration at the damaged region which influences mechanical strength. The flaw size produced by such damage was examined to determine the magnitude of strength degradation in micrometeorite or space-debris impacted glasses. Scanning electron microscopy revealed topographical details of impact sites which included central melt zones and glass fiber production. The overall crater structure is similar to much larger impacts of large meteorite on the Moon in that the melt crater is surrounded by shocked regions of material which fracture zones and spall areas. Residual stresses arising from shock compression and cooling of the fused zone cannot currently be included in fracture mechanics analyses based on simple flaw size examination.

Wiedlocher, David E.↗

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

Warm Molecular Gas Traced with CO J = 7 --> 6 in the Galaxy's Central 2 Parsecs: Dynamical Heating of the Circumnuclear Disks

We present an 11" resolution map of the central 2 pc of the Galaxy in the CO J = 7 --> 6 rotational transition. The CO emission shows rotation about Sgr A* but also evidence for noncircular turbulent motion and a clumpy morphology. We combine our data set with available CO measurements to model the physical conditions in the disk. We find that the molecular gas in the region is both warm and dense, with T approx. 200-300 K and n(sub H2) approx. (5-7) x 10(exp 4) cm(exp -3). The mass of warm molecular gas we measure in the central 2 pc is at least 2000 M(solar), about 20 times the UV-excited atomic gas mass, ruling out a UV heating scenario for the molecular material. We compare the available spectral tracers with theoretical models and conclude that molecular gas is heated with magnetohydrodynamic shocks with v approx. 10-20 km s(exp -1) and B approx. 0.3- 0.5 mG. Using the conditions derived with the CO analysis, we include the other important coolants, neutral oxygen and molecular hydrogen, to estimate the total cooling budget of the molecular material. We derive a mass-to-luminosity ratio of approx. 2-3 M(solar)(L(solar)exp -1), which is consistent with the total power dissipated via turbulent decay in 0.1 pc cells with v(sub rms) approx. 15 kilometers per second. These size and velocity scales are comparable to the observed clumping scale and the velocity dispersion. At this rate, the material near Sgr A* is dissipating its orbital energy on an orbital timescale and cannot last for more than a few orbits. Our conclusions support a scenario in which the features near Sgr A* such as the circumnuclear disk and northern arm are generated by infalling clouds with low specific angular momentum.

galatic nuclei↗

The interstellar medium on the Gamma Cassiopeiae line of sight

The interstellar medium on the Gamma Cas line of sight is studied through the observations of O I, H2, H I, D I, and Ar I absorption features with the Copernicus satellite. By using the velocity structure of the line of sight previously determined through atomic nitrogen, it is demonstrated that these neutrals are located in the same physical regions and that the O I and Ar I abundances could be solar, on the average, relative to N I in all components. The 1066.660 A Ar I line profile is contaminated by an unidentified line, already detected in the spectrum of Zeta Pup. Furthermore, a peculiar behavior of the Ar I abundance in the different components is reported, possibly associated with high-velocity interstellar gas present on the line of sight. It is concluded that the Gamma Cas line of sight could intercept a dense cloud which has been engulfed and disrupted by a shock.

Ferlet, R.↗

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolution, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolutions, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗

The effect of a type 3 and type 4 shock/shock interaction on heat transfer in the stagnation region

One of the major engineering challenges in designing the National Aerospace Plane, NASP, is to overcome augmented heating on the intake cowl lip from shock/shock interactions. The shock/shock interaction arises when the bow shock from the craft's nose interferes with the bow shock from the cowl lip. Considering only the region immediately around the cowl lip, the problem geometry may be simplified as that of an oblique shock impinging on a bow shock from a circular cylinder. Edney classified six different interference patterns resulting from an oblique-shock/curved bow-shock interaction. Of these six types, type 3 and 4 are most significant in that augmented surface heat transfer may be ten to thirty times greater than the case without the shock/shock interaction. The objective was to begin to develop a mathematical model which is capable of predicting the effect of a type 3 and 4 shock/shock interaction in the stagnation region of an arbitrary 2-D body. This model must be capable of predicting the maximum surface heat flux and the surface stagnation point pressure once the outer (effectively inviscid) flowfield is given. Therefore, it must capture the unsteady physics of the impinging shear layer.

Wilson, Dennis↗

Studies of Interstellar Pickup Ions in the Solar Wind

The work under this grant involves studies of the interaction of interstellar pickup ions with the solar wind, with the goal of a comprehensive model of the particle distributions and wave intensities to be expected throughout the heliosphere, as well as the interactions of those distributions with the solar wind termination shock. In the past year, we have completed a number of projects, including observations and modeling of the effects of a large scattering mean free path on the pickup He(+) seen at AMPTE, an analytical model of anisotropic pickup tons in a steady radial magnetic field, and a derivation of a reduced solar wind Mach number due to increased estimates on the inflowing hydrogen density allowing for a weak termination shock. In the next year, we plan to investigate in more detail the correspondence between our models of anisotropic pickup ions and the data on spectra, variations, and proton-He(+) correlation provided by AMPTE, Ulysses, and our instrument on SOHO. We will model the time-dependent pickup ion density resulting from finite periods of radial magnetic field. We will also incorporate the effects of a large mean free path into our analysis of the He(+) focusing cone, leading to more accurate parameter values for the interstellar helium gas. This progress report also includes a discussion of our Space Physics Educational Outreach activities in the past year and plans for the next year.

Isenberg, Philip A.↗

Geomechanical models of impact cratering: Puchezh-Katunki structure

Impact cratering is a complex natural phenomenon that involves various physical and mechanical processes. Simulating these processes may be improved using the data obtained during the deep drilling at the central mound of the Puchezh-Katunki impact structure. A research deep drillhole (named Vorotilovskaya) has been drilled in the Puchezh-Katunki impact structure (European Russia, 57 deg 06 min N, 43 deg 35 min E). The age of the structure is estimated at about 180 to 200 m.y. The initial rim crater diameter is estimated at about 40 km. The central uplift is composed of large blocks of crystalline basement rocks. Preliminary study of the core shows that crystalline rocks are shock metamorphosed by shock pressure from 45 GPa near the surface to 15-20 GPa at a depth of about 5 km. The drill core allows the possibility of investigating many previously poorly studied cratering processes in the central part of the impact structure. As a first step one can use the estimates of energy for the homogeneous rock target. The diameter of the crater rim may be estimated as 40 km. The models elaborated earlier show that such a crater may be formed after collapse of a transient cavity with a radius of 10 km. The most probable range of impact velocities from 11.2 to 30 km/s may be inferred for the asteroidal impactor. For the density of a projectile of 2 g/cu cm the energy of the impact is estimated as 1E28 to 3E28 erg. In the case of vertical impact, the diameter of an asteroidal projectile is from 1.5 to 3 km for the velocity range from 11 to 30 km/s. For the most probable impact angle of 45 deg, the estimated diameter of an asteroid is slightly larger: from 2 to 4 km. Numerical simulation of the transient crater collapse has been done using several models of rock rheology during collapse. Results show that the column at the final position beneath the central mound is about 5 km in length. This value is close to the shock-pressure decay observed along the drill core. Further improvement of the model needs to take into account the blocky structure of target rocks revealed by drilling.

Ivanov, B. A.↗

Effect of Scattering on Instabilities in Line-driven Stellar Winds

Line driven O-B stellar winds are unstable to perturbations of short spatial wavelength and the growth rates for such instabilities are very rapid. The nature of the nonlinear development of this instability is unknown, but might possibly be one of blobs of gas driven through ambient gas or a quasiregular train of outward moving shocks. In either case the resulting dissipation of mechanical energy might explain the observed anomalous heating n O-B stars as evidenced by their X-ray emission and high ionization state. It also might explain the observed fine structure of the absorption lines and their time variability. The driving due to the absorption of the stellar continuum flux was considered and the effects of the diffuse, scattered radiation field were neglected. It was shown that under certain special conditions the effect of scattering could reduce the instability growth rate to zero. A stability analysis that includes scattering, but that uses the more physically realistic assumption of complete redistribution instead of coherent scattering, and that includes the effects of transverse velocity gradients, which become important as the flow moves away from the stellar surface is presented. It is found that the instability is eliminated right at the base of the wind, but that as the flow moves outward the instability rate rapidly becomes equal to a substantial fraction of the calculated value, the fraction asymptotically reaching 80% at large radii. Since this still implies many e folds in a characteristic outflow time, the primary conclusion that these winds are highly unstable is unchanged.

Owocki, S. P.↗

Transient Three-Dimensional Startup Side Load Analysis of a Regeneratively Cooled Nozzle

The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate.

Wang, Ten-See↗

Aeroacoustics of Turbulent High-Speed Jets

Aeroacoustic noise generation in a supersonic round jet is studied to understand in particular the effect of turbulence structure on the noise without numerically compromising the turbulence itself. This means that direct numerical simulations (DNS's) are needed. In order to use DNS at high enough Reynolds numbers to get sufficient turbulence structure we have decided to solve the temporal jet problem, using periodicity in the direction of the jet axis. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. Therefore in order to answer some questions about the turbulence we will partially compromise the overall structure of the jet. The first section of chapter 1 describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. In the second section we present preliminary work done using a TVD numerical scheme on a CM5. This work is only two-dimensional (plane) but shows very interesting results, including weak shock waves. However this is a nonviscous computation and the method resolves the shocks by adding extra numerical dissipation where the gradients are large. One wonders whether the extra dissipation would influence small turbulent structures like small intense vortices. The second chapter is an extensive discussion of preliminary numerical work using the spectral method to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which are solved in O(N) steps. A very accurate highly resolved DNS of a turbulent jet flow is expected.

Rao, Ram Mohan↗

Numerical Simulations of High-Speed Chemically Reacting Flow

The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗

Laser-Induced Thermal Acoustics Theory and Expected Experimental Errors when Applied to a Scramjet Isolator Model

A scramjet isolator model test apparatus is being assembled in the Isolator Dynamics Research Lab (IDRL) at the NASA Langley Research Center in Hampton, Virginia. The test apparatus is designed to support multiple measurement techniques for investigating the flow field in a scramjet isolator model. The test section is 1-inch high by 2-inch wide by 24-inch long and simulates a scramjet isolator with an aspect ratio of two. Unheated, dry air at a constant stagnation pressure and temperature is delivered to the isolator test section through a Mach 2.5 planar nozzle. The isolator test section is mechanically back-pressured to contain the resulting shock train within the 24-inch isolator length and supports temperature, static pressure, and high frequency pressure measurements at the wall. Additionally, nonintrusive methods including laser-induced thermal acoustics (LITA), spontaneous Raman scattering, particle image velocimetry, and schlieren imaging are being incorporated to measure off-wall fluid dynamic, thermodynamic, and transport properties of the flow field. Interchangeable glass and metallic sidewalls and optical access appendages permit making multiple measurements simultaneously. The measurements will be used to calibrate computational fluid dynamics turbulence models and characterize the back-pressured flow of a scramjet isolator. This paper describes the test apparatus, including the optical access appendages; the physics of the LITA method; and estimates of LITA measurement uncertainty for measurements of the speed of sound and temperature.

Middleton, Troy F.↗