Prandtl-meyer expansion for air in local equilibrium and in frozen flows
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Theoretical predictions of weak shocks and Prandtl-Meyer expansions in unsteady flow and three-dimensional steady supersonic flow
Analysis of rarefied gas flow in axisymmetric jet exhausting into vacuum including Prandtl-Meyer expansion
Previous studies of the high-speed viscous inviscid interaction between a turbulent boundary layer and an expansion at a convex corner have noted that surface pressure decreases toward the downstream inviscid value yielded by a Prandtl-Meyer expansion. A downstream influence on the corner is presently identified which is based on the mean surface pressure distribution; a scaling law is proposed for this distance.
The plume flowfield of a helium vent nozzle exhausting into a vacuum is defined by two techniques: the method of characteristics and the direct-simulation Monte Carlo method. The method of characteristics is shown to severely underpredict gas static temperatures due to the assumption of translational equilibrium inherent in the continuum formulation. Mass flow per unit solid angle is shown to be accurately predicted for this flowfield by the method of characteristics until the flow angle approaches within 10 deg of the maximum Prandtl-Meyer expansion angle. Improved treatment of the noncontinuum flow within the nozzle (near the lip) is postulated to have minor effects on the accuracy of this continuum method in comparison with translational nonequilibrium in the external flowfield. Possible treatment of translational nonequilibrium by the method of characteristics is discussed.
We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.
A new linearized-theory pressure-coefficient formulation was studied. The new formulation is intended to provide more accurate estimates of detailed pressure loadings for improved stability analysis and for analysis of critical structural design conditions. The approach is based on the use of oblique-shock and Prandtl-Meyer expansion relationships for accurate representation of the variation of pressures with surface slopes in two-dimensional flow and linearized-theory perturbation velocities for evaluation of local three-dimensional aerodynamic interference effects. The applicability and limitations of the modification to linearized theory are illustrated through comparisons with experimental pressure distributions for delta wings covering a Mach number range from 1.45 to 4.60 and angles of attack from 0 to 25 degrees.
Improved approximations are incorporated into the inviscid fluid method originated by Spreiter et al (1970) in calculations of the altitude of the Venus ionopause. The altitude calculations are then compared with median altitudes measured by the Pioneer Venus retarding potential analyzer. The calculated ionopause shape is found to closely approximate the measured meridian shape in the solar zenith angle (SZA) range of zero-135 deg. The use of improved ionospheric pressure field ionosheath pressure and Prandtl-Meyer expansion approximations lowers the terminator ionopause altitude to approximately half that obtained with the usual Spreiter approximations. It is also determined that the calculated dawn ionopause is about 300 km higher than the dusk ionopause, and that both ionopauses are close to their respective measured meridian values. It is concluded that median ionopause altitude within the SZA angle range may be calculated without inclusion of a viscous interaction in the theory.
Compressible fluid flow problems were solved using an adaptive finite element method. Supersonic flow past a wedge and Prandtl-Meyer expansion were considered. Computational grids are presented.
An integrated computational model of a generic hypersonic vehicle was developed for the purpose of determining the vehicle's performance characteristics, which include the lift, drag, thrust, and moment acting on the vehicle at specified altitude, flight condition, and vehicular configuration. The lift, drag, thrust, and moment are developed for the body fixed coordinate system. These forces and moments arise from both aerodynamic and propulsive sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be determined. The vehicle is assumed to be a lifting body with a single aerodynamic control surface. The body shape and control surface location are arbitrary and must be defined. The aerodynamics are calculated using either 2-dimensional Newtonian or modified Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. Skin-friction drag was also accounted for. The skin-friction drag coefficient is a function of the freestream Mach number. The data for the skin-friction drag coefficient values were taken from NASA Technical Memorandum 102610. The modeling of the vehicle's SCRAMjet engine is based on quasi 1-dimensional gas dynamics for the engine diffuser, nozzle, and the combustor with heat addition. The engine has three variable inputs for control: the engine inlet diffuser area ratio, the total temperature rise through the combustor due to combustion of the fuel, and the engine internal expansion nozzle area ratio. The pressure distribution over the vehicle's lower aft body surface, which acts as an external nozzle, is calculated using a combination of quasi 1-dimensional gas dynamic theory and Newtonian or modified Newtonian theory. The exhaust plume shape is determined by matching the pressure inside the plume, calculated from the gas dynamic equations, with the freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this manner, the pressure distribution along the vehicle after body expansion surface is then determined. The aerodynamic modeling, the engine modeling, and the exhaust plume analysis are described in more detail. A description of the computer code used to perform the above calculations is given and an input/output example is then given. The computer code is available on a Macintosh floppy disk.
An experimental study has been performed with a scramjet combustor and nozzle model to determine the effect of nozzle radii on relaminarization, and the effectiveness of film injection in reducing wall heat flux. The nozzle model was installed downstream of a swept-ramp scramjet combustor model, operated at Mach 2.7 with vitiated air, simulating Mach 6-7 flight enthalpies. Results from tests conducted with three nozzle radii, varying in curvature from a sharp Prandtl-Meyer expansion to a radius of two scramjet combustor exit heights, showed little or no effect on the downstream wall heat flux, which is not in agreement with both subsonic and supersonic relaminarization correlations in the available literature. In contrast, film injection upstream of all three radii produced reductions of 70 percent in the entire nozzle model heating rate, for the nominal film-to-free stream pressure matched condition. Results from film injection, with and without pressure gradient, but without wall curvature, show that the gradients appear to stabilize the film, with the resulting reductions in wall heat transfer rates apparent far downstream of the point of injection.
Experiments were conducted to observe the cross-sectional structure and streamwise growth of round transverse liquid jets injected into a highly accelerated boundary layer in supersonic flow. The accompanying shock structure was also visualized. In one case, a round jet of acetone was injected into a fully turbulent Mach 2.5 boundary layer that was subsequently accelerated and partially laminarized through a sharp Prandtl-Meyer expansion corner. In the second case, a jet was injected into the laminarized Mach 3.2 boundary layer downstream of the expansion corner at the same jet-to-freestream momentum ratio. The jet and shock structure in both cases were visualized using schlieren optics. Wall-flow patterns were visualized using paints. It was found that the lateral spreading of jets injected downstream of the expansion fan was augmented close to the wall and had a cross-sectional structure significantly different from that of the jet injected upstream: the upstream jet spreads rapidly at the expansion corner in both the lateral and vertical directions.
This report, which is a revision and extension of NACA-TN-1428, presents a compilation of equations, tables, and charts useful in the analysis of high-speed flow of a compressible fluid. The equations provide relations for continuous one-dimensional flow, normal and oblique shock waves, and Prandtl-Meyer expansions for both perfect and imperfect gases. The tables present useful dimensionless ratios for continuous one-dimensional flow and for normal shock waves as functions of Mach number for air considered as a perfect gas. One series of charts presents the characteristics of the flow of air (considered a perfect gas) for oblique shock waves and for cones in a supersonic air stream. A second series shows the effects of caloric imperfections on continuous one-dimensional flow and on the flow through normal and oblique shock waves. (author)
The design characteristics of an inlet very much depend on whether the inlet is to be flown at subsonic, supersonic, or hypersonic speed. Whichever the case, the primary function of an inlet is to deliver free-stream air to the engine face at the highest stagnation pressure possible and with the lowest possible variation in both stagnation pressure and temperature. At high speeds, this is achieved by a system of oblique and/or normal shock waves, and possibly some isentropic compression. For both subsonic and supersonic flight, current design practice indicates that the inlet should deliver the air to the engine face at approximately Mach 0.45. As a result, even for flight in the high subsonic regime, the inlet must retard (or diffuse) the air substantially. Second, the design of an inlet is influenced largely by the compromise between high performance and low weight. This compromise involves tradeoffs between the mission requirements, flight trajectory, airframe aerodynamics, engine performance, and weight-all of which, in turn, influence each other. Therefore, to study the effects of some of these influential factors, the Propulsion System Analysis Office of the NASA Lewis Research Center developed the Inlet Performance Analysis Code (IPAC). This code uses oblique shock and Prandtl-Meyer expansion theory to predict inlet performance. It can be used to predict performance for a given inlet geometric design such as pitot, axisymmetric, and two-dimensional. IPAC also can be used to design preliminary inlet systems and to make subsequent performance analyses. It computes the total pressure, the recovery, the airflow, and the drag coefficients. The pressure recovery includes losses associated with normal and oblique shocks, internal and external friction, the sharp lip, and diffuser components. Flow rate includes captured, engine, spillage, bleed, and bypass flows. The aerodynamic drag calculation includes drags associated with spillage, cowl lip suction, wave, bleed, and bypass.
This report provides a user guide for the Compressible Flow Toolbox, a collection of algorithms that solve almost 300 linear and nonlinear classical compressible flow relations. The algorithms, implemented in the popular MATLAB programming language, are useful for analysis of one-dimensional steady flow with constant entropy, friction, heat transfer, or shock discontinuities. The solutions do not include any gas dissociative effects. The toolbox also contains functions for comparing and validating the equation-solving algorithms against solutions previously published in the open literature. The classical equations solved by the Compressible Flow Toolbox are: isentropic-flow equations, Fanno flow equations (pertaining to flow of an ideal gas in a pipe with friction), Rayleigh flow equations (pertaining to frictionless flow of an ideal gas, with heat transfer, in a pipe of constant cross section.), normal-shock equations, oblique-shock equations, and Prandtl-Meyer expansion equations. At the time this report was published, the Compressible Flow Toolbox was available without cost from the NASA Software Repository.
Two nonintrusive laser diagnostics were independently used to demonstrate the measurement of time-averaged and spatially-resolved pressure change across a twodimensional (2-D) shock wave. The first method is Doppler global velocimetry (DGV) which uses water seeding and generates 2-D maps of 3-orthogonal components of velocity. A DGV-measured change in flow direction behind an oblique shock provides an indirect determination of pressure jump across the shock, when used with the known incoming Mach number and ideal shock relations (or Prandtl-Meyer flow equations for an expansion fan). This approach was demonstrated at Mach 2 on 2-D shocks and expansions generated from a flat plate at angles-of-attack approx. equals -2.4deg and +0.6deg, respectively. This technique also works for temperature jump (as well as pressure) and for normal shocks (as well as oblique). The second method, laser-induced thermal acoustics (LITA), is a seedless approach that was used to generate 1-D spatial profiles of streamwise Mach number, sound speed, pressure, and temperature across the same shock waves. Excellent agreement was obtained between the DGV and LITA methods, suggesting that either technique is viable for noninvasive shock-strength measurements.
The objective of this study is to propose a computational methodology that can effectively anchor the base flowfield of a four-engine clustered nozzle configuration. This computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics (CFD) formulation. For efficient CFD calculation, a Prandtl-Meyer solution treatment is applied to the algebraic grid lines for initial plume expansion resolution. As the solution evolves, the computational grid is adapted to the pertinent flow gradients. The CFD model employs an upwind scheme in which second- and fourth-order central differencing schemes with artificial dissipation are used. The computed quantitative base flow properties such as the radial base pressure distributions, model centerline static pressure, Mach number and impact pressure variations, and base pressure characteristic curve agreed reasonably well with those of the measurement.