Search NASA⌕ Search

Engineering topics

Drozda, Tomasz G.

Publications and source records attributed to Drozda, Tomasz G..

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to find optimal designs for an interdigitated flushwall injector suitable for scramjet applications at hypervelocity conditions. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. A Latin hypercube sampling design-of-experiments method was used to select design points for RAS. A methodology was developed that automated building geometries and generating grids for each design. The ensuing RAS analysis generated the performance database from which the two objective functions of interest were computed using a one-dimensional performance utility. The data were fitted using four surrogate models: an artificial neural network (ANN) model, a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model obtained an optimal solution set that predicted high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts are required in order to lower the errors and perform more accurate surrogate-based optimization. sed optimization.

Shenoy, Rajiv R.↗

Comparison of Boundary Layer Similarity Transformations for High Mach Number Flows

Reynolds averaged Navier-Stokes (RANS) simulations of the turbulent boundary layer flow on a flat plate are performed across a Mach number range of 0.3 up to 16 using the computational fluid dynamics (CFD) software, VULCAN-CFD. The simulation results are used to evaluate the ability of the Van Driest transformation and the transformation of Trettel and Larsson to collapse boundary layer velocity profiles under various flow conditions across the range of Mach numbers. The flat platewall boundary conditions and leading edge geometry are varied from adiabatic to isothermal and sharp to blunted to reveal the physics of how these effects impact the performance of the chosen flow transformations. Results indicate that the transformation of Trettel and Larsson produces a better collapse of velocity profile data than the Van Driest transformation regardless of Mach number, wall boundary condition, or leading edge geometry. In addition, the velocity gradients obtained from the transformation of Trettel and Larsson match the normalized untransformed velocity gradients much more closely than those obtained using the Van Driest transformation.

DiGregorio, Nicholas J.↗

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.

Drozda, Tomasz G.↗

Comparison of Several Global Mixing Performance Metrics for High-Speed Fuel Injectors

To experimentally assess and compare the mixing performance of high-speed fuel injectors for scramjet engines, quantitative global metrics are needed. The one-dimensional metric most commonly used to assess the degree of mixing completeness at a given downstream station is the mixing efficiency parameter. The experimental determination of the mixing efficiency parameter requires measurement of the spatial distributions of both the fuel mass fraction and the mass flux. Standard in-stream gas sampling techniques can be used to measure the fuel mass fraction distribution, however the mass flux distribution is not easily determined experimentally because it requires the measurement of three independent aerothermodynamic variables in addition to the mixture composition. For this reason, several metrics that can be calculated from the fuel distribution alone are commonly used to assess mixing performance. Because these other metrics do not provide a mass flux-weighted measure of the local degree of mixing completeness, they may not correlate well with the mixing efficiency parameter. Therefore, if the substitute metrics are to be used to compare the mixing performance of candidate fuel injector concepts, it is important to understand their relationships to the mixing efficiency parameter in a representative scramjet combustor flowfield. This work investigates the relationships between the mixing efficiency parameter and several substitute metrics that are able to be measured with the current experimental setup of the Enhanced Injection and Mixing Project at the NASA Langley Research Center for baseline strut and ramp injectors. The results of these comparisons have revealed that it is possible to glean different (i.e., incorrect) conclusions about which injector is the better mixer when the substitute mixing performance metrics are used instead of the mixing efficiency parameter, thereby highlighting the importance of mass flux-weighted mixing performance metrics.

Ground, Cody R.↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Drozda, Tomasz G.↗

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.

Shenoy, Rajiv R.↗

Comparisons Between NO PLIF Imaging and CFD Simulations of Mixing Flowfields for High-Speed Fuel Injectors

The current work compares experimentally and computationally obtained nitric oxide (NO) planar laser induced fluorescence (PLIF) images of the mixing flowfields for three types of high-speed fuel injectors: a strut, a ramp, and a rectangular flushwall. These injection devices, which exhibited promising mixing performance at lower flight Mach numbers, are currently being studied 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, and improve the understanding of underlying physical processes relevant to flight Mach numbers greater than eight. In the experiments, conducted in the NASA Langley Arc-Heated Scramjet Test Facility (AHSTF), the injectors are placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number air flow at the entrance of a scramjet combustor. Helium is used as an inert substitute for hydrogen fuel. Both schlieren and PLIF techniques are applied to obtain mixing flowfield flow visualizations. The experimental PLIF is obtained by using a UV laser sheet to interrogate a plane of the flow by exciting fluorescence from the NO molecules, which are present in the AHSTF air. Consequently, the absence of signal in the resulting PLIF images is an indication of pure helium (fuel). The computational PLIF is obtained by applying a fluorescence model for NO to the results of the Reynolds-averaged simulations (RAS) of the mixing flow field carried out using the VULCAN-CFD solver. This approach is required because the PLIF signal is a nonlinear function of not only NO concentration, but also pressure, temperature, and the flow velocity. This complexity allows additional flow features to be identified and compared with those obtained from the computational fluid dynamics (CFD) simulations, however, such comparisons are only semiquantitative. Three-dimensional image reconstruction, similar to that used in magnetic resonance imaging, is also used to obtain images in the streamwise and spanwise planes from select cross-stream PLIF plane data. Synthetic schlieren is also computed from the RAS data. Good agreement between the experimental and computational results provides increased confidence in the CFD simulations for investigations of injector performance.

Drozda, Tomasz G.↗

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.

Drozda, Tomasz G.↗

Computational and Experimental Characterization of the Mach 6 Facility Nozzle Flow for the Enhanced Injection and Mixing Project at NASA Langley Research Center

Computational fluid dynamics analyses and experimental data are presented for the Mach 6 facility nozzle used in the Arc-Heated Scramjet Test Facility for the Enhanced Injection and Mixing Project (EIMP). This project, conducted at the NASA Langley Research Center, aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics relevant to flight Mach numbers greater than 8. The EIMP experiments use a two-dimensional Mach 6 facility nozzle to provide the high-speed air simulating the combustor entrance flow of a scramjet engine. Of interest are the physical extent and the thermodynamic properties of the core flow at the nozzle exit plane. The detailed characterization of this flow is obtained from three-dimensional, viscous, Reynolds-averaged simulations. Thermodynamic nonequilibrium effects are also investigated. The simulations are compared with the available experimental data, which includes wall static pressures as well as in-stream static pressure, pitot pressure and total temperature obtained via in-stream probes positioned just downstream of the nozzle exit plane.

Drozda, Tomasz G.↗

CFD Analysis of Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The calculations were carried out using the VULCAN-CFD solver and Reynolds-Averaged Simulations (RAS). The high Mach number flow conditions match those proposed for the planned 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 and functional relationships relevant to flight Mach numbers greater than eight. Because of the high Mach number flow considered, the injectors consist of a fuel placement device, a strut; and a fluidic vortical mixer, a ramp. These devices accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow albeit via different strategies. Both of these devices were previously studied at lower flight Mach numbers where they exhibited promising performance in terms of mixing efficiency and total pressure recovery. 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, however, at high flight Mach number conditions, the dynamic pressure needed to induce sufficient fuel penetration may be difficult to achieve along with other requirements such as achieving desired levels of fuel-to-air mixing at the required equivalence ratio. The three injectors represent the baseline configurations planned for the experiments. The current work discusses the mixing flow field behavior and differences among the three fuel injectors, mixing performance as described by the mixing efficiency and the total pressure recovery, and performance considerations based on the thrust potential.

Drozda, Tomasz G.↗

Impact of Flight Enthalpy, Fuel Simulant, and Chemical Reactions on the Mixing Characteristics of Several Injectors at Hypervelocity Flow Conditions

The high total temperatures or total enthalpies required to duplicate the high-speed flight conditions in ground experiments often place stringent requirements on the material selection and cooling needs for the test articles and intrusive flow diagnostic equipment. Furthermore, for internal flows, these conditions often complicate the use of nonintrusive diagnostics that need optical access to the test section and interior portions of the flowpath. Because of the technical challenges and increased costs associated with experimentation at high values of total enthalpy, an attempt is often made to reduce it. This is the case for the Enhanced Injection and Mixing Project (EIMP) currently underway in the Arc-Heated Scramjet Test Facility at the NASA Langley Research Center. The EIMP aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships between mixing performance and losses relevant to flight Mach numbers greater than 8. The experiments will consider a "direct-connect" approach and utilize a Mach 6 nozzle to simulate the combustor entrance flow of a scramjet engine. However, while the value of the Mach number is matched to that expected at the combustor entrance in flight, the maximum value of the total enthalpy for these experiments is limited by the thermal-structural limits of the uncooled experimental hardware. Furthermore, the fuel simulant is helium, not hydrogen. The use of "cold" flows and non-reacting mixtures of fuel simulants for mixing experiments is not new and has been extensively utilized as a screening technique for scramjet fuel injectors. In this study, Reynolds-averaged simulations are utilized (RAS) to systematically verify the implicit assumptions used by the EIMP. This is accomplished by first performing RAS of mixing for two injector configurations at planned nominal experimental conditions. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared to the values obtained from RAS under the true enthalpy conditions and using helium and hydrogen. Finally, the impact of combustion on mixing, often deemed small enough to neglect at hypervelocity conditions, is assessed by comparing the results obtained from the hydrogen-fueled reacting and non-reacting RAS. For reacting flows, in addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also considered. In all of the simulations, the incoming air Mach number and the fuel-to-air ratio are the same, while the total pressure, total enthalpy, and the fuel simulant vary depending on the case considered. It is found that under some conditions the "cold" flow experiments are a good approximation of the flight.

Drozda, Tomasz G.↗

A Priori Analysis of a Compressible Flamelet Model using RANS Data for a Dual-Mode Scramjet Combustor

In an effort to make large eddy simulation of hydrocarbon-fueled scramjet combustors more computationally accessible using realistic chemical reaction mechanisms, a compressible flamelet/progress variable (FPV) model was proposed that extends current FPV model formulations to high-speed, compressible flows. Development of this model relied on observations garnered from an a priori analysis of the Reynolds-Averaged Navier-Stokes (RANS) data obtained for the Hypersonic International Flight Research and Experimentation (HI-FiRE) dual-mode scramjet combustor. The RANS data were obtained using a reduced chemical mechanism for the combustion of a JP-7 surrogate and were validated using avail- able experimental data. These RANS data were then post-processed to obtain, in an a priori fashion, the scalar fields corresponding to an FPV-based modeling approach. In the current work, in addition to the proposed compressible flamelet model, a standard incompressible FPV model was also considered. Several candidate progress variables were investigated for their ability to recover static temperature and major and minor product species. The effects of pressure and temperature on the tabulated progress variable source term were characterized, and model coupling terms embedded in the Reynolds- averaged Navier-Stokes equations were studied. Finally, results for the novel compressible flamelet/progress variable model were presented to demonstrate the improvement attained by modeling the effects of pressure and flamelet boundary conditions on the combustion.

Quinlan, Jesse R.↗

Pre-Test CFD for the Design and Execution of the Enhanced Injection and Mixing Project at NASA Langley Research Center

With the increasing costs of physics experiments and simultaneous increase in availability and maturity of computational tools it is not surprising that computational fluid dynamics (CFD) is playing an increasingly important role, not only in post-test investigations, but also in the early stages of experimental planning. This paper describes a CFD-based effort executed in close collaboration between computational fluid dynamicists and experimentalists to develop a virtual experiment during the early planning stages of the Enhanced Injection and Mixing project at NASA Langley Research Center. This projects aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships relevant to flight Mach numbers greater than 8. The purpose of the virtual experiment was to provide flow field data to aid in the design of the experimental apparatus and the in-stream rake probes, to verify the nonintrusive measurements based on NO-PLIF, and to perform pre-test analysis of quantities obtainable from the experiment and CFD. The approach also allowed for the joint team to develop common data processing and analysis tools, and to test research ideas. The virtual experiment consisted of a series of Reynolds-averaged simulations (RAS). These simulations included the facility nozzle, the experimental apparatus with a baseline strut injector, and the test cabin. Pure helium and helium-air mixtures were used to determine the efficacy of different inert gases to model hydrogen injection. The results of the simulations were analyzed by computing mixing efficiency, total pressure recovery, and stream thrust potential. As the experimental effort progresses, the simulation results will be compared with the experimental data to calibrate the modeling constants present in the CFD and validate simulation fidelity. CFD will also be used to investigate different injector concepts, improve understanding of the flow structure and flow physics, and develop functional relationships. Both RAS and large eddy simulations (LES) are planned for post-test analysis of the experimental data.

Drozda, Tomasz G.↗

A Priori Analysis of Flamelet-Based Modeling for a Dual-Mode Scramjet Combustor

An a priori investigation of the applicability of flamelet-based combustion models to dual-mode scramjet combustion was performed utilizing Reynolds-averaged simulations (RAS). For this purpose, the HIFiRE Direct Connect Rig (HDCR) flowpath, fueled with a JP-7 fuel surrogate and operating in dual- and scram-mode was considered. The chemistry of the JP-7 fuel surrogate was modeled using a 22 species, 18-step chemical reaction mechanism. Simulation results were compared to experimentally-obtained, time-averaged, wall pressure measurements to validate the RAS solutions. The analysis of the dual-mode operation of this flowpath showed regions of predominately non-premixed, high-Damkohler number, combustion. Regions of premixed combustion were also present but associated with only a small fraction of the total heat-release in the flow. This is in contrast to the scram-mode operation, where a comparable amount of heat is released from non-premixed and premixed combustion modes. Representative flamelet boundary conditions were estimated by analyzing probability density functions for temperature and pressure for pure fuel and oxidizer conditions. The results of the present study reveal the potential for a flamelet model to accurately model the combustion processes in the HDCR and likely other high-speed flowpaths of engineering interest.

Quinlan, Jesse R.↗

Characterization of the NASA Langley Arc Heated Scramjet Test Facility Using NO PLIF

The nitric oxide planar laser-induced fluorescence (NO PLIF) imaging was used to characterize the air flow of the NASA Langley Arc Heated Scramjet Test Facility (AHSTF) configured with a Mach 6 nozzle. The arc raises the enthalpy of the test gas in AHSTF, producing nitric oxide. Nitric oxide persists as the temperature drops through the nozzle into the test section. NO PLIF was used to qualitatively visualize the flowfield at different experimental conditions, measure the temperature of the gas flow exiting the facility nozzle, and visualize the wave structure downstream of the nozzle at different operating conditions. Uniformity and repeatability of the nozzle flow were assessed. Expansion and compression waves on the free-jet shear layer as the nozzle flow expands into the test section were visualized. The main purpose of these experiments was to assess the uniformity of the NO in the freestream gas for planned experiments, in which NO PLIF will be used for qualitative fuel-mole-fraction sensitive imaging. The shot-to-shot fluctuations in the PLIF signal, caused by variations in the overall laser intensity as well as NO concentration and temperature variations in the flow was 20-25% of the mean signal, as determined by taking the standard deviation of a set of images obtained at constant conditions and dividing by the mean. The fluctuations within individual images, caused by laser sheet spatial variations as well as NO concentration and temperature variations in the flow, were about 28% of the mean in images, determined by taking standard deviation within individual images, dividing by the mean in the same image and averaged over the set of images. Applying an averaged laser sheet intensity correction reduced the within-image intensity fluctuations to about 10% suggesting that the NO concentration is uniform to within 10%. There was no significant difference in flow uniformity between the low and high enthalpy settings. While not strictly quantitative, the temperature maps show qualitative agreement with the computations of the flow.

Kidd, F. Gray, III↗

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗