Compressibility effects on fluid entrainment by turbulent mixing layers
Compressibility effects on fluid entrainment by supersonic turbulent mixing layers
SEARCH · Search NASA
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.
Compressibility effects on fluid entrainment by supersonic turbulent mixing layers
Turbulent mixing of compressible air in supersonic and subsonic coaxial jet flows
An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.
The primary objective of this work is to provide accurate numerical solutions for selected flow fields and to compare and evaluate the performance of selected turbulence models with experimental results. Four popular turbulence models have been tested and validated against experimental data often turbulent flows. The models are: (1) the two-equation k-epsilon model of Wilcox, (2) the two-equation k-epsilon model of Launder and Sharma, (3) the two-equation k-omega/k-epsilon SST model of Menter, and (4) the one-equation model of Spalart and Allmaras. The flows investigated are five free shear flows consisting of a mixing layer, a round jet, a plane jet, a plane wake, and a compressible mixing layer; and five boundary layer flows consisting of an incompressible flat plate, a Mach 5 adiabatic flat plate, a separated boundary layer, an axisymmetric shock-wave/boundary layer interaction, and an RAE 2822 transonic airfoil. The experimental data for these flows are well established and have been extensively used in model developments. The results are shown in the following four sections: Part A describes the equations of motion and boundary conditions; Part B describes the model equations, constants, parameters, boundary conditions, and numerical implementation; and Parts C and D describe the experimental data and the performance of the models in the free-shear flows and the boundary layer flows, respectively.
A computational tool named SUPIN is being developed to design and analyze the aerodynamic performance of supersonic inlets. The inlet types available include the axisymmetric pitot, three-dimensional pitot, axisymmetric outward-turning, two-dimensional single-duct, two-dimensional bifurcated-duct, and streamline-traced inlets. The aerodynamic performance is characterized by the flow rates, total pressure recovery, and drag. The inlet flow-field is divided into parts to provide a framework for the geometry and aerodynamic modeling. Each part of the inlet is defined in terms of geometric factors. The low-fidelity aerodynamic analysis and design methods are based on analytic, empirical, and numerical methods which provide for quick design and analysis. SUPIN provides inlet geometry in the form of coordinates, surface angles, and cross-sectional areas. SUPIN can generate inlet surface grids and three-dimensional, structured volume grids for use with higher-fidelity computational fluid dynamics (CFD) analysis. Capabilities highlighted in this paper include the design and analysis of streamline-traced external-compression inlets, modeling of porous bleed, and the design and analysis of mixed-compression inlets. CFD analyses are used to verify the SUPIN results.
A two-equation turbulence model (k-epsilon) has been modified to model compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. Computations are carried out using the Navier-Stokes solver SPARK which includes a finite rate chemistry model for hydrogen-air combustion. Initially, planar compressible mixing layers are simulated. Comparisons of the predictions with available experimental data and the predictions of a compressible Reynolds stress closure indicate that the model is well suited for the study of such flows. The decrease in the growth rate of the mixing layer with increasing convective Mach number is well predicted. The model is then tested on a three-dimensional supersonic flow problem where solid boundaries are present. An axisymmetric nonreacting mixing layer flow involving dissimilar gases is evaluated next using the model and the predictions are compared with available experimental data. Finally, two reacting turbulent flow configurations have been studied and comparison of predictions and experimental data has been done. The comparisons indicate that the model predicts the overall features of the test cases very well.
Two-dimensional nonisoenergetic turbulent jet mixing between two compressible streams at constant pressure
The first experiment is a study of a coaxial jet discharging into stagnant laboratory air, with center jet of a mixture of 5% oxygen and 95% helium by volume and coflow jet of air. The exit flow pressure of both center-jet and coflow nozzles is 1 atmosphere. The presence of oxygen in the center jet is to allow the use of an oxygen flow-tagging technique (RELIEF4) to obtain non-intrusive velocity measurements. Both jets are nominally Mach 1.8, but, because of the greater speed of sound, the center jet velocity is more than twice that of the coflow. The mixing layer which forms between the center jet and the coflow near the nozzle exit is compressible, with a calculated convective Mach number of approximately 0.7. This geometry has several advantages: The streamwise development of the flow is generally dominated by turbulent stresses (rather than pressure forces), and thus calculations are sensitive to turbulence modeling. It includes features present in supersonic combustors, including a compressible mixing layer near the nozzle exit and a light-gas/air plume downstream. Since it is a free jet, it provides easy access for both optical instrumentation and probes. Since it is axisymmetric, it requires fewer experimental measurements to fully characterize, and calculations can be performed with more modest computer resources. However, weak shock waves formed at the nozzle exit strengthen and turn normal as they approach the axis, complicating the flow. Care is thus taken in the design of the facility to provide as near as possible to 1-D flow at the exit of both center and coflow nozzles, and to minimize the strength of waves generated at the nozzle exit. Results from this experiment are compared to CFD solutions obtained by VULCAN, a previously developed code used in engine analysis. The second experiment is a study of a supersonic combustor consisting of a diverging duct with single downstream-angled wall injector. Thus, the geometry is relatively simple and large regions of subsonic recirculating flow are avoided. The nominal entrance Mach number is 2 and the enthalpy of the test gas (hot air "simulant") is nominally that of Mach 7 flight. It was believed, on the basis of calculations performed that this would produce mixing-limited flow, that is to say, one for which chemical reaction to equilibrium proceeds at a much greater rate than mixing. It later proved that this was not the case. The primary experimental technique employed is coherent anti-Stokes Raman spectroscopy, known by its acronym CARS. The species probed is molecular nitrogen and the quantity measured is temperature. Intrusive probes, such as Pitot, total temperature, hot-wire, etc., are not used due to access difficulty and high heat flux in the combustor, and because they may alter the flow. CARS has several advantages over other optical methods. It is a relatively mature and well-understood technique. Signal levels are relatively high and the signal is in the form of a coherent (laser) beam that can be collected through small windows. Incoherent (non-CARS) interferences are rejected by spatial filtering.
The distortion characteristics were investigated at the subsonic diffuser exit of a full-scale, Mach 2.5, axisymmetric, mixed compression inlet. Performance and steady-state distortion characteristics were obtained at zero and maximum angle of attack and during an inlet unstart-restart sequence. For the configuration with no cowl bleed, steady-state distortion P(max)P(min)P(bar) ranged from 0.10 for critical inlet operation at 0 deg angle-of-attack to 0.306 for supercritical inlet operation at 6.84 deg angle-of-attack. Vortex generators provided a 50 percent reduction in steady-state distortion for critical operation. Bleed has a smaller effect on steady-stated distortion.
A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.
The focus of the program was on the use of direct numerical simulations of turbulent flow for study of turbulence physics and modeling. A special interest was placed on turbulent mixing layers. The required data for these investigations were generated from four newly developed codes for simulation of time and spatially developing incompressible and compressible mixing layers. Also of interest were the structure of wall bounded turbulent and transitional flows, evaluation of diagnostic techniques for detection of organized motions, energy transfer in isotropic turbulence, optical propagation through turbulent media, and detailed analysis of the interaction of vortical structures.
A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.
Mixing and combustion of turbulent, compressible free jets
Three computer applications have been developed to solve inviscid compressible fluids problems using interactive computer graphics. The first application is a compressible flow calculator which solves for isentropic flow, normal shocks, and oblique shocks or centered expansions produced by two dimensional ramps. The second application couples the solutions generated by the first application to a more graphical presentation of the results to produce a desk top simulator of three compressible flow problems: 1) flow past a single compression ramp; 2) flow past two ramps in series; and 3) flow past two opposed ramps. The third application extends the results of the second to produce a design tool which solves for the flow through supersonic external or mixed compression inlets. The applications were originally developed to run on SGI or IBM workstations running GL graphics. They are currently being extended to solve additional types of flow problems and modified to operate on any X-based workstation.
Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.
The bifurcated inlet examined in this study (which is one of several being considered in the High Speed Research (HSR) Program) was chosen based upon paper trade studies of axisymmetric, single sided, and bifurcated inlets. For a given compression ratio and mass flow a bifurcated inlet weighs less than a single sided inlet. An axisymmetric inlet has less bleed requirements than 2D inlets but has trouble matching transonic airflow requirements without going to a variable diameter centerbody. The bifurcated inlet was selected as one of the candidates because of its ability to match airflow schedules. The inlet examined in this study, the Two Stage Supersonic Inlet (TSSI), was a candidate mixed compression bifurcated inlet. It has a novel concept to aid in inlet stability. This concept was tested in the 10x10 wind tunnel at NASA Glenn. CFD tools were used to predict and interpret the experimental results.
Dynamic response of supersonic mixed compression inlet coupled to cold pipe or turbojet engine
Wind tunnel tests to increase stable operating range of supersonic mixed compression inlet