Kenji Miki | Search NASA

NASA NTRS2022

An Investigation Into Some Important Aspects of Droplet Breakup/Vaporization Behavior Caused By a Gas Flow (Both Nonreacting & Detonative Combustion) in a Shock Tube

As a part of the rotating detonation engine (RDE) technology enablement project at NASA Glenn Research Center (GRC), an effort was undertaken to extend our current computational capabilities of OpenNCC in some important ways with the implementation of a modeling approach to account for the droplet breakup caused by a shock-induced gas motion, & a vaporization model valid over a wide range of pressure conditions encountered in multiphase detonation. With the modified code, a study was undertaken to investigate the individual droplet behavior followed by the passage of a shock front. The study is carried out by tracking a sparse group of droplets to gain some understanding of shock induced droplet behavior under various shock strengths & fuel injector conditions. The study also looks into the effect of randomization involved in determining the droplet breakup outcomes. Over a wide range of sub-critical conditions examined, larger droplets are observed to undergo significant changes in droplet behavior following their breakup. However, smaller droplets (10 µm or less ) remain unaffected by any shock induced breakup. In a follow-on work, we investigated the impact of shock and droplet interaction in a detonation study involving both gaseous as well as gas/liquid (droplet clouds) fuel/air stoichiometric mixtures in a simple 3D shock-tube configuration. The droplet clouds are made up of different initial droplet sizes of either 6, 10, or 30 µm. We also investigated the individual droplet behavior followed by the passage of a detonation front. The results represent conditions that lead to both overdriven and C-J (ChapmanJouguet) detonations. Under both test conditions, most of the droplet vaporization is completed within a short distance (duration) behind the detonation front & well within the region of complete combustion observed in a corresponding equivalent gas-phase fuel/air mixture. The impact of the shock-induced droplet breakup is found to be significant in the calculations involving the 30-µm droplets. Subsequent to the breakup, the drop sizes vary from 1 to 10 µm. Another factor that contributed to the observed rapid vaporization is the result of vaporization taking place under supercritical conditions. The overall detonation properties of various droplet clouds (made up of different initial sizes) are similar to those observed in a corresponding gaseous fuel/air mixture. In the calculation involving a gaseous fuel, the calculated C-J detonation velocity is 1822 m/s involving Jet-A/air and φ = 1. In the overdriven detonation, it is 2044 m/s. In the calculations involving droplet clouds, the corresponding detonation velocities are lower. The impact of increased droplet size is primarily seen in a higher reduction in the detonation velocity.

shock/droplet interaction↗

NASA NTRS2022

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of free-stream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Re(ex)=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

LES-LDKM↗

NASA NTRS2021

Large-Eddy Simulation of Variable Speed Power Turbine Cascade with Inflow Turbulence

Numerical results are presented from theNASAGlenn Research Center’sin-house turbomachinery code, Glenn-HT applied to the Variable Speed Power Turbine (VSPT) experiment at the NASA Transonic Turbine Blade Cascade Facility.The main goal of this paper is to implement a digital filtering method to generate turbulence upstream and a sub-grid model (Localized dynamic k-equation model (LDKM)) in the framework of LES in order to investigate the effect of inflow turbulence on the transition seen in the VSPT experimental data at the cruise condition (incidence angle of 40oand Tu = 0.5%, 5%,10%, and 15%). Although the boundary layer on thesuction side and pressure side of the blades is initially laminar due to favorable pressure gradient, the laminar flow can transition to turbulent flow past a separation zoneon the suction sideor by natural or by-pass transition. This process determinesthe total-pressure lossesin the wake. Therefore, it is desirable to develop a reliable prediction tool to accurately capture the transition mechanism in blade rows operated under the conditions of low Reynolds number andat a variety offree stream turbulence conditions. Our numerical studies reveal that the location of separationis rather insensitive to the level of Tu, however the effect ofincreasingTu seems to be in reducing the size and ultimately suppressing the separation bubble. In addition, we performedspectral analysis to identify the peak frequencies in the region where the separation bubble is formed, which provides valuable insights into the transition/separation mechanism.

CFD↗

NASA NTRSDate not supplied

RDE Nozzle Computational Design Methodology Development and Application

The Open National Combustion Code (OpenNCC) is used to simulate the aerospike exhaust nozzle region of a rotating detonation engine (RDE).The main objectives of the current study are twofold. The first goal is to validate the proposed computational methodology using the experimental data. The second goal is to demonstrate how the validated prediction tool can be used to optimize the nozzle geometry. To achieve a significant speed-up of computational time, the computational domain in divided into two parts: a combustion region and a throat-nozzle section. In the combustion region, a validated quasi-two-dimensional in-house code is utilized to generate an unsteady RDE flow field solution just upstream of the combustor throat. Subsequently, the unsteady flow data is fed into the three-dimensional throat-nozzle section as the inflow boundary condition. This facilitates the design optimization process since the unsteady inflow can be reused, and a relatively coarse mesh (i.e., larger time-step) can be used to analyze the flow fields around the nozzle. Five nozzle designs were studied and the predicted performance (i.e., thrust)compared. The optimized nozzle was found to produce 3.2% more overall thrust than a baseline nozzle design. Without the nozzle, there is a large low-pressure region at the throat exit, which significantly reduces the overall performance. This methodology is shown to be a promising approach to explore a wide variety of nozzle geometries in a relatively short amount of time.

CFD↗

NASA NTRSDate not supplied

Computational Study of Modeling Fully-coupled Combustor-Turbine Interactions by the Open National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

NASA NTRSDate not supplied

Computational Study of Modeling Fully-coupled Combustor-Turbine Interactions by the Open National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

NASA NTRSDate not supplied

RDE Nozzle Computational Design Methodology Development and Application

The rotating detonation engine (RDE) is currently under investigation as an approach to achieving pressure gain combustion for propulsion and power systems. The RDE typically consists of an annulus with one end open (or having a throat and/or nozzle) and the other end valved using non-mechanical, fluidic means. In this study, a total of 6 nozzle design variations were studied for the NPS combustor geometry, using the validated quasi-2D to 3-D coupled simulation approach.

CFD↗

NASA NTRSDate not supplied

Numerical Study of Combustor-Turbine Interactions Using Open-National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

NASA NTRSDate not supplied

Large-Eddy Simulation of the Variable Speed Power Turbine Cascade with Inflow Turbulence

The main goal of this is to implement a digital filtering method to generate turbulence upstream and a sub-grid model (Localized dynamic k-equation model (LDKM)) in the framework of LES which together enable reproduction of the effect of inflow turbulence on the transition seen in the VSPT experimental data at the cruise condition (incidence angle of 40o and Tu=0.5%, 5%, 10% and 15%.).

CFD↗

NASA NTRSDate not supplied

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of freestream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Reex=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

CFD↗

NASA NTRSDate not supplied

Improved Prediction of Losses with Large Eddy Simulation in a Low Pressure Turbine

We will present the wake profiles resulting from our previous computations for a range of FSTs. We will address key questions such as how the numerics, mesh density in addition to the FST level affect the wake prediction. We will utilize the spectral analysis and investigate the dynamic behavior of the wake as influenced by the flow variations in locations upstream.

LES↗

NASA NTRSDate not supplied

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

Detonation↗

NASA NTRSDate not supplied

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

detonation↗

NASA NTRSDate not supplied

Embracing the fourth Industrial Revolution - Challenges, Opportunities and Path Forward for Propulsion

The aerospace industry is at a point where components are reaching design maturity and performance improvements are incremental. Aggressive goals to achieve sustainability and the threat of climate change necessitate a new paradigm. An artificial intelligence (AI) approach that enables revolutionary changes in system architecture, mission analysis and performance metrics is needed. The growing interest and development in the field of machine learning presents an opportunity to speed up by 10X or more the discovery, analysis and development of aerospace systems using artificial intelligence. Through IDEAS (Intelligent Design and Engineering of Aerospace Systems) we are embarking on a research and development effort that addresses this opportunity. The objective of IDEAS is to enable design of systems based on requirements. We discuss possible approaches to generating data, training models and applying them to near term applications. Results of recent workshops with industry, academia and other agencies to identify challenges to adopting AI and machine learning are presented.

Machine learning↗

NASA NTRSDate not supplied

Numerical Investigation of the Effect of Trailing Edge Thickness of Simulated CMC Blades on Loss Profiles

To understand such an aerodynamic characteristic of the CMC blades, we numerically simulated three different CMC blades with varying trailing edge thicknesses. We utilized both Large Eddy simulation (LES) with the dynamic subgrid closure as well a recently implemented Reynolds-averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model. We considered two sets of experimental test points with different Reynolds numbers at a given high-free stream turbulence range (Tu=~10-13%). The LES was previously validated by comparing with experimental data. We validated the RANS transition model against the ERCOFTAC experimental zero-pressure gradient smooth flat plate boundary layer at low- and high freestream turbulence offered by the experimental data of the Stripf’s for a turbine blade. We confirmed that the predicted skin friction and the heat flux agree well with the data. Finally, we computed the predicted pressure loading profiles, total loss distributions, and integrated loss of the three different CMC blades and made comparisons against the experimental data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Our paper will present the aforementioned results and discuss the ability of LES and RANS model of our choice to predict salient aerodynamic characteristics of the blades with the various trailing edge thickness effects.

LES↗

NASA NTRSDate not supplied

Embracing the Fourth Industrial Revolution - Challenges, Opportunities and Path Forward for Propulsion

The aerospace industry is at a point where components are reaching design maturity and performance improvements are incremental. Aggressive goals to achieve sustainability and the threat of climate change necessitate a new paradigm. An artificial intelligence (AI) approach that enables revolutionary changes in system architecture, mission analysis and performance metrics is needed. The growing interest and development in the field of machine learning presents an opportunity to speed up by 10X or more the discovery, analysis and development of aerospace systems using artificial intelligence. Through IDEAS (Intelligent Design and Engineering of Aerospace Systems) we are embarking on a research and development effort that addresses this opportunity. The objective of IDEAS is to enable design of systems based on requirements. We discuss possible approaches to generating data, training models and applying them to near term applications. Results of recent workshops with industry, academia and other agencies to identify challenges to adopting AI and machine learning are presented.

Machine learning↗

NASA NTRSDate not supplied

Computational Fluid Dynamic Optimization of an Experimental Rotating Detonation Rocket Engine Nozzle

A parametric optimization study is performed on the nozzle of a laboratory rotating detonation rocket engine (RDRE) using a three-dimensional computational fluid dynamic simulation. The primary optimization objective is maximum nozzle thrust. The basic nozzle configuration is a shrouded, truncated plug. The fluid in the RDRE chamber leading to the nozzle is choked at its exit so that its cyclic behavior is unaffected by any changes to the nozzle design. Optimization is performed for a single operating point. Parameters varied are the overall nozzle area expansion ratio and the fraction of the expansion area that is provided by the shroud. These two parameters indirectly affect the angle of the plug nozzle cone, and the bluff body area associated with its truncation. Nozzle thrust is evaluated as the difference between the thrust of the RDRE chamber-plus-nozzle combination and that of the chamber alone. The nozzle produces approximately 20% of the total engine thrust. The baseline nozzle is found to perform well, yielding 58.1% of the thrust calculated for a notional ideal RDRE nozzle which can instantaneously change shape to allow isentropic expansion of every fluid element. Optimization improves the performance, bringing the nozzle thrust to 70.0% of the notional ideal, and total engine thrust (chamber-plus-nozzle) to 94% of the ideal.

detonation↗

NASA NTRSDate not supplied

Computational Fluid Dynamic Optimization of an Experimental Rotating Detonation Rocket Engine Nozzle

A parametric optimization study is performed on the nozzle of a laboratory rotating detonation rocket engine (RDRE) using a three-dimensional computational fluid dynamic simulation. The primary optimization objective is maximum nozzle thrust. The basic nozzle configuration is a shrouded, truncated plug. The fluid in the RDRE chamber leading to the nozzle is choked at its exit so that its cyclic behavior is unaffected by any changes to the nozzle design. Optimization is performed for a single operating point. Parameters varied are the overall nozzle area expansion ratio and the fraction of the expansion area that is provided by the shroud. These two parameters indirectly affect the angle of the plug nozzle cone, and the bluff body area associated with its truncation. Nozzle thrust is evaluated as the difference between the thrust of the RDRE chamber-plus-nozzle combination and that of the chamber alone. The nozzle produces approximately 20% of the total engine thrust. The baseline nozzle is found to perform well, yielding 58.1% of the thrust calculated for a notional ideal RDRE nozzle which can instantaneously change shape to allow isentropic expansion of every fluid element. Optimization improves the performance, bringing the nozzle thrust to 70.0% of the notional ideal, and total engine thrust (chamber-plus-nozzle) to 94% of the ideal.

detonation↗