The Effect of Turbulence Modeling on the Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented on the effects of turbulence modeling choices on the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The analyses were carried out with the VULCAN-CFD solver using Reynolds-Averaged Simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The three injectors are the baseline configurations used in the experiments and represent three categories of injectors typically considered individually or in combination for fueling high-speed propulsive devices. The current work discusses the impact of the turbulence model and the turbulent Schmidt number on the mixing flow field behavior and the mixing performance as described by the one-dimensional values of the Mach number, total pressure recovery, and the mixing efficiency. Because planar laser induced fluorescence (PLIF) images are available from the EIMP experiments, the sensitivity of the synthetic LIF signal to turbulence modeling choices is also examined to determine whether PLIF can be extended beyond its intended qualitative visualization purpose and used to guide CFD turbulence model and parameter selections. It is found that the mixing performance, as quantified using mixing efficiency, exhibits a strong sensitivity to both turbulence model choice and turbulent Schmidt number value. However, the synthetic LIF signal only demonstrates a modest level of sensitivity, which suggests that PLIF is of limited use for guiding CFD turbulence model and parameter selections

Hybrid Reynolds-Averaged/Large-Eddy Simulations of a Co-Axial Supersonic Free-Jet Experiment

Reynolds-averaged and hybrid Reynolds-averaged/large-eddy simulations have been applied to a supersonic coaxial jet flow experiment. The experiment utilized either helium or argon as the inner jet nozzle fluid, and the outer jet nozzle fluid consisted of laboratory air. The inner and outer nozzles were designed and operated to produce nearly pressure-matched Mach 1.8 flow conditions at the jet exit. The purpose of the computational effort was to assess the state-of-the-art for each modeling approach, and to use the hybrid Reynolds-averaged/large-eddy simulations to gather insight into the deficiencies of the Reynolds-averaged closure models. The Reynolds-averaged simulations displayed a strong sensitivity to choice of turbulent Schmidt number. The baseline value chosen for this parameter resulted in an over-prediction of the mixing layer spreading rate for the helium case, but the opposite trend was noted when argon was used as the injectant. A larger turbulent Schmidt number greatly improved the comparison of the results with measurements for the helium simulations, but variations in the Schmidt number did not improve the argon comparisons. The hybrid simulation results showed the same trends as the baseline Reynolds-averaged predictions. The primary reason conjectured for the discrepancy between the hybrid simulation results and the measurements centered around issues related to the transition from a Reynolds-averaged state to one with resolved turbulent content. Improvements to the inflow conditions are suggested as a remedy to this dilemma. Comparisons between resolved second-order turbulence statistics and their modeled Reynolds-averaged counterparts were also performed

Comparison of Mixing Characteristics for Several Fuel Injectors at Mach 8, 12, and 15 Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at flight Mach numbers of 8, 12, and 15. The Reynolds-averaged simulations (RAS) were carried out using the VULCAN-CFD solver. The high Mach number flow conditions match those of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies relevant to flight Mach numbers greater than 8. The injectors include a fuel placement device, a strut, and a fluidic vortical mixer, a ramp. These fuel injectors accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow, albeit via different strategies. For comparison, a flush-wall injector is also included. This type of injector generally represents the simplest method of introducing fuel into a scramjet combustor. The three injectors represent the baseline configurations of the EIMP experiments. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are computed from the RAS and compared for the three flight conditions and injector configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations. Plotting the total pressure recovery and thrust potential as a function of mixing efficiency provides added insight into critical aspects of combustor performance as the flight condition and injector type are varied

Laser-Induced Thermal Acoustic Measurements in a Highly Back-Pressured Scramjet Isolator Model: A Research Plan

Under the Propulsion Discipline of NASA s Fundamental Aeronautics Program s Hypersonics Project, a test apparatus, for testing a scramjet isolator model, is being constructed at NASA's Langley Research Center. The test apparatus will incorporate a 1-inch by 2-inch by 15-inch-long scramjet isolator model supplied with 2.1 lbm/sec of unheated dry air through a Mach 2.5 converging-diverging nozzle. The planned research will incorporate progressively more challenging measurement techniques to characterize the flow field within the isolator, concluding with the application of the Laser-Induced Thermal Acoustic (LITA) measurement technique. The primary goal of this research is to use the data acquired to validate Computational Fluid Dynamics (CFD) models employed to characterize the complex flow field of a scramjet isolator. This paper describes the test apparatus being constructed, pre-test CFD simulations, and the LITA measurement technique

Comparison of Mixing Characteristics for Several Fuel Injectors on an Open Plate and in a Ducted Flowpath Configuration at Hypervelocity Flow Conditions

In order to reduce the cost and complexity associated with fuel injection and mixing experiments for high-speed flows, and to further enable optical access to the test section for nonintrusive diagnostics, the Enhanced Injection and Mixing Project (EIMP) utilizes an open flat plate configuration to characterize inert mixing properties of various fuel injectors for hypervelocity applications. The experiments also utilize reduced total temperature conditions to alleviate the need for hardware cooling. The use of "cold" flows and non-reacting mixtures for mixing experiments is not new, and has been extensively utilized as a screening technique for scramjet fuel injectors. The impact of reduced facility-air total temperature, and the use of inert fuel simulants, such as helium, on the mixing character of the flow has been assessed in previous numerical studies by the authors. Mixing performance was characterized for three different injectors: a strut, a ramp, and a flushwall. The present study focuses on the impact of using an open plate to approximate mixing in the duct. Toward this end, Reynolds-averaged simulations (RAS) were performed for the three fuel injectors in an open plate configuration and in a duct. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared for the two configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations

The NASA Langley Isolator Dynamics Research Lab

The Isolator Dynamics Research Lab (IDRL) is under construction at the NASA Langley Research Center in Hampton, Virginia. A unique test apparatus is being fabricated to support both wall and in-stream measurements for investigating the internal flow of a dual-mode scramjet isolator model. The test section is 24 inches long with a 1-inch by 2-inch cross sectional area and is supplied with unheated, dry air through a Mach 2.5 converging-diverging nozzle. The test section is being fabricated with two sets (glass and metallic) of interchangeable sidewalls to support flow visualization and laser-based measurement techniques as well as static pressure, wall temperature, and high frequency pressure measurements. During 2010, a CFD code validation experiment will be conducted in the lab in support of NASA s Fundamental Aerodynamics Program. This paper describes the mechanical design of the Isolator Dynamics Research Lab test apparatus and presents a summary of the measurement techniques planned for investigating the internal flow field of a scramjet isolator model

Shuttle Gaseous Hydrogen Venting Risk from Flow Control Valve Failure

This paper describes a series of studies to assess the potential risk associated with the failure of one of three gaseous hydrogen flow control valves in the orbiter's main propulsion system during the launch of Shuttle Endeavour (STS-126) in November 2008. The studies focused on critical issues associated with the possibility of combustion resulting from release of gaseous hydrogen from the external tank into the atmosphere during assent. The Shuttle Program currently assumes hydrogen venting from the external tank will result in a critical failure. The current effort was conducted to increase understanding of the risk associated with venting hydrogen given the flow control valve failure scenarios being considered in the Integrated In-Flight Anomaly Investigation being conducted by NASA

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements

Numerical modeling and simulation of supersonic flows in propulsion systems by open-source solvers

Two open-source solvers, Eilmer and hyFoam, are here considered for their performance in simulating high-speed flows in different flow conditions and geometric configurations typical of propulsive systems at supersonic speeds. The goal is to identify the open-source platform providing the best compromise between accuracy, flexibility and computational cost to eventually simulate the flow fields inside ramjet and scramjet engines. The differences in terms of discretization and solution methods of the selected solvers are discussed in terms of their impact on solution accuracy and computational efficiency and in view of the aerothermodynamic analysis and design of future trans-atmospheric propulsive systems. In this work steady state problems are considered. Numerical results of two scramjet type engines demonstrated a similar predictive capability of both codes in non-reacting conditions. These results highlight their potential to be considered for further characterization of overall engine performance