Engineering Papers⌕ Search

Engineering topics

Menon, Suresh

Publications and source records attributed to Menon, Suresh.

Validation of High-Fidelity CFD Simulations for Rocket Injector Design

Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by evaluating the sensitivity of performance and injector-driven thermal environments to the details of the injector geometry and key operational parameters. Methodical verification and validation efforts on a range of coaxial injector elements have shown the current production CFD capability must be improved in order to quantitatively impact the injector design process. This paper documents the status of a focused effort to compare and understand the predictive capabilities and computational requirements of a range of CFD methodologies on a set of single element injector model problems. The steady Reynolds-Average Navier-Stokes (RANS), unsteady Reynolds-Average Navier-Stokes (URANS) and three different approaches using the Large Eddy Simulation (LES) technique were used to simulate the initial model problem, a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants. While one high-fidelity LES result matches the experimental combustion chamber wall heat flux very well, there is no monotonic convergence to the data with increasing computational tool fidelity. Systematic evaluation of key flow field regions such as the flame zone, the head end recirculation zone and the downstream near wall zone has shed significant, though as of yet incomplete, light on the complex, underlying causes for the performance level of each technique. 1 Aerospace Engineer and Combustion CFD Team Leader, MS ER42, NASA MSFC, AL 35812, Senior Member, AIAA. 2 Professor and Director, Computational Combustion Laboratory, School of Aerospace Engineering, 270 Ferst Dr., Atlanta, GA 30332, Associate Fellow, AIAA. 3 Reilly Professor of Engineering, School of Mechanical Engineering, 585 Purdue Mall, West Lafayette, IN 47907, Fellow, AIAA. 4 Principal Member of Technical Staff, Combustion Research Facility, 7011 East Avenue, MS9051, Livermore, CA 94550, Associate Fellow, AIAA. 5 J. L. and G. H. McCain Endowed Chair, Mechanical Engineering, 104 Research Building East, University Park, PA 16802, Fellow, AIAA. American Institute of Aeronautics and Astronautics 1

Tucker, P. Kevin↗

An Approach to Improved Credibility of CFD Simulations for Rocket Injector Design

Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by simulating the sensitivity of performance and injector-driven thermal environments to. the details of the injector geometry and key operational parameters. Methodical verification and validation efforts on a range of coaxial injector elements have shown the current production CFD capability must be improved in order to quantitatively impact the injector design process.. This paper documents the status of an effort to understand and compare the predictive capabilities and resource requirements of a range of CFD methodologies on a set of model problem injectors. Preliminary results from a steady Reynolds-Average Navier-Stokes (RANS), an unsteady Reynolds-Average Navier Stokes (URANS) and three different Large Eddy Simulation (LES) techniques used to model a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants are presented. Initial observations are made comparing instantaneous results, corresponding time-averaged and steady-state solutions in the near -injector flow field. Significant differences in the flow fields exist, as expected, and are discussed. An important preliminary result is the identification of a fundamental mixing mechanism, accounted for by URANS and LES, but missing in the steady BANS methodology. Since propellant mixing is the core injector function, this mixing process may prove to have a profound effect on the ability to more correctly simulate injector performance and resulting thermal environments. Issues important to unifying the basis for future comparison such as solution initialization, required run time and grid resolution are addressed.

Tucker, Paul K.↗

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗

Large-Eddy/Lattice Boltzmann Simulations of Micro-blowing Strategies for Subsonic and Supersonic Drag Control

This report summarizes the progress made in the first 8 to 9 months of this research. The Lattice Boltzmann Equation (LBE) methodology for Large-eddy Simulations (LES) of microblowing has been validated using a jet-in-crossflow test configuration. In this study, the flow intake is also simulated to allow the interaction to occur naturally. The Lattice Boltzmann Equation Large-eddy Simulations (LBELES) approach is capable of capturing not only the flow features associated with the flow, such as hairpin vortices and recirculation behind the jet, but also is able to show better agreement with experiments when compared to previous RANS predictions. The LBELES is shown to be computationally very efficient and therefore, a viable method for simulating the injection process. Two strategies have been developed to simulate multi-hole injection process as in the experiment. In order to allow natural interaction between the injected fluid and the primary stream, the flow intakes for all the holes have to be simulated. The LBE method is computationally efficient but is still 3D in nature and therefore, there may be some computational penalty. In order to study a large number or holes, a new 1D subgrid model has been developed that will simulate a reduced form of the Navier-Stokes equation in these holes.

Menon, Suresh↗

Premixed Turbulent Flame Propagation in Microgravity

A combined numerical-experimental study has been carried out to investigate the structure and propagation characteristics of turbulent premixed flames with and without the influence of buoyancy. Experimentally, the premixed flame characteristics are studied in the wrinkled regime using a Couette flow facility and an isotropic flow facility in order to resolve the scale of flame wrinkling. Both facilities were chosen for their ability to achieve sustained turbulence at low Reynolds number. This implies that conventional diagnostics can be employed to resolve the smallest scales of wrinkling. The Couette facility was also built keeping in mind the constraints imposed by the drop tower requirements. Results showed that the flow in this Couette flow facility achieves full-developed turbulence at low Re and all turbulence statistics are in good agreement with past measurements on large-scale facilities. Premixed flame propagation studies were then carried out both using the isotropic box and the Couette facility. Flame imaging showed that fine scales of wrinkling occurs during flame propagation. Both cases in Ig showed significant buoyancy effect. To demonstrate that micro-g can remove this buoyancy effect, a small drop tower was built and drop experiments were conducted using the isotropic box. Results using the Couette facility confirmed the ability to carry out these unique reacting flow experiments at least in 1g. Drop experiments at NASA GRC were planned but were not completed due to termination of this project.

Menon, Suresh↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

Turbulent Premixed Flames in Microgravity

The experimental cold-flow facility is now full operational and is currently being used to obtain baseline turbulence data in a Couette flow. The baseline turbulence data is necessary to confirm the capability of the chosen device to generate and maintain the required turbulence intensity. Subsequent reacting flow studies will assume that a similar turbulent flow field exists ahead of the premixed flame. Some modifications and refinements had to be made to enable accurate measurements. It consists of two rollers, one (driven by a motor) which drives a continuous belt and four smaller rollers used to set the belt spacing and tension to minimize belt flutter. The entire assemble is enclosed in a structure that has the dimensions to enable future drop tower experiments of the hot facility. All critical dimensions are the same as the original plans except for the pulley ratio which has been changed to enable a wider operating regime in terms of the Reynolds number. With the current setup, Reynolds numbers as low as 100 and as high as 14,000 can be achieved. This is because the in-between belt spacing can be varied from 1 cm to 7.6 cm, and the belt speed can be accurately varied from .15 m/sec to 3.1 m/sec.

Menon, Suresh↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

Simulation of Unsteady Combustion in a Ramjet Engine Using a Highly Parallel Computer

Combustion instability in ramjets is a complex phenomenon that involve nonlinear interaction between acoustic waves, vortex motion and unsteady heat release in the combustor. To numerically simulate this 3-D, transient phenomenon, enormous computer resources (time, memory and disk storage) are required. Although current generation vector supercomputers are capable of providing adequate resources for simulations of this nature, their high cost and limited availability, makes such machines less than satisfactory for routine use. The primary focus of this study is to assess the feasibility of using highly parallel computer systems as a cost-effective alternative for conducting such unsteady flow simulations. Towards this end, a large-eddy simulation model for combustion instability was implemented on the Intel iPSC/860 and a careful study was conducted to determine the benefits and the problems associated with the use of such machines for transient simulations. Details of this study along with the results obtained from the unsteady combustion simulations carried out on the iPSC/860 are discussed in this paper.

Menon, Suresh↗

A new unsteady mixing model to predict NO(x) production during rapid mixing in a dual-stage combustor

An advanced gas turbine engine to power supersonic transport aircraft is currently under study. In addition to high combustion efficiency requirements, environmental concerns have placed stringent restrictions on the pollutant emissions from these engines. A combustor design with the potential for minimizing pollutants such as NO(x) emissions is undergoing experimental evaluation. A major technical issue in the design of this combustor is how to rapidly mix the hot, fuel-rich primary zone product with the secondary diluent air to obtain a fuel-lean mixture for combustion in the second stage. Numerical predictions using steady-state methods cannot account for the unsteady phenomena in the mixing region. Therefore, to evaluate the effect of unsteady mixing and combustion processes, a novel unsteady mixing model is demonstrated here. This model has been used to study multispecies mixing as well as propane-air and hydrogen-air jet nonpremixed flames, and has been used to predict NO(x) production in the mixing region. Comparison with available experimental data show good agreement, thereby providing validation of the mixing model. With this demonstration, this mixing model is ready to be implemented in conjunction with steady-state prediction methods and provide an improved engineering design analysis tool.

Menon, Suresh↗

A new unsteady mixing model to predict NO(x) production during rapid mixing in a dual-stage combustor

An advanced gas turbine engine to power supersonic transport aircraft is currently under study. In addition to high combustion efficiency requirements, environmental concerns have placed stringent restrictions on the pollutant emissions from these engines. A dual-stage combustor with the potential for minimizing pollutants such as NO(x) emissions is undergoing experimental evaluation. A major technical issue in the design of this combustor is how to rapidly mix the hot, fuel-rich primary stage product with the secondary diluent air to obtain a fuel-lean mixture for combustion in the secondary stage. Numerical design studies using steady-state methods cannot account for the unsteady phenomena in the mixing region. Therefore, to evaluate the effect of unsteady mixing and combustion processes, a novel unsteady mixing model is demonstrated here. This model has been used in a stand-alone mode to study mixing and combustion in hydrogen-air nonpremixed jet flames. NO(x) production in these jet flames was also predicted. Comparison of the computed results with experimental data show good agreement thereby providing validation of the mixing model.

Menon, Suresh↗

Numerical simulation of pressure oscillations in a ramjet combustor

Large-eddy simulations of compressible cold flow in a ramjet combustor configuration were performed. The objectives were to investigate the mechanisms for vortex-acoustic interaction in such a device and to develop a simulation method that can be extended to cases with combustion heat release to study combustion instability. From the simulation, it was found that the separated shear layer rolls up into concentrated vortices that merge to form large coherent structures. These vortices interact with the choked nozzle downstream to produce an axial acoustic dipole. The spectrum of the pressure fluctuation at the base of the backward-facing step shows that there are two types of oscillations: an acoustic resonant mode and a vortex-acoustic coupled mode. Based on the flow physics observed in the simulations, a simple one-dimensional model for the vortex-acoustic coupled mode was proposed. The eigenvalue problem based on this model was solved to obtain the frequency of the coupled mode.

Jou, Wen-Huei↗

Numerical investigation of V/STOL jet induced interactions

Direct numerical simulation using the full three dimensional, time dependent Navier-Stokes equations is used to investigate V/STOL jet induced interactions. The objective of this numerical simulation is to compute accurately the details of the flow field and to achieve a better understanding of the physics of the flow, including the role of initial turbulence in the jet, the influence of forward motion on hover aerodynamics, the collision zone and fountain characteristics. Preliminary results are presented.

Rizk, Magdi H.↗

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗