Ali Ameri | Search NASA

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 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

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

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

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

CMC↗

NASA NTRSDate not supplied

Combustor-Turbine Interactions By Using the Open National Combustion Code (OpenNCC) and the Glenn-HT Code

We develop the methodology of coupling two different computational fluid dynamic (CFD) codes, OpenNCC and Glenn-HT, in order to investigate the unsteady flow fields inside the combustor and around the first-stage stator of a high-pressure turbine (HPT) from the energy efficient engine (E 3 ) program. In our coupling strategy, OpenNCC, a time accurate unstructured mesh multi-phase combustion CFD code, is applied to the combustor region and Glenn-HT, a structured mesh CFD code designed for turbine applications, is utilized inside the HPT region. As a proof of concept, we first simulate the three-dimensional airflow over the backward facing step and compared the predicted pressure coefficient against experimental data. We then considered a Jet-A/Air combustor with the NASA single learn direct injector (LDI) as a reacting flow test case. Here, Glenn-HT is used in the downstream of the combustor and assumes that the working fluid is a single gas (i.e., air) where the properties of the gas are taken from time- and area-averaged values at the combustor exit. A fully coupled combustor-turbine simulation using both OpenNCC and Glenn-HT was conducted on the E 3 combustor and a detailed analysis done pertaining to how a realistic combustor outflow affects migration of hot-streaks and its impact on the aerodynamic behavior inside the HPT. It is found that the combustion dynamics (i.e., the switching main and pilot flame strength) significantly alters the aerodynamic behavior inside the HPT including hot-streak migration. Variations of the temperature and the velocity magnitude in the passage could be up to ±75 and ±15. These large variations are an important part of the combustor-turbine interactions, which are overlooked by a single component simulation of the HPT while imposing the steady inflow boundary condition at the HPT inlet. Finally, it is observed the strong density gradient associated with the hot-streaks and the pressure gradient at the passage significantly enhances the baroclinic torque, affecting the vorticity field.

CFD hot-streaks↗

NASA NTRSDate not supplied

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

LES↗

NASA NTRSDate not supplied

LES Analysis of Loss Anomaly in a Turbine Blade with Large Trailing Edge Radius

There is great interest in the usage of Ceramic Matrix Composites (CMC) as a turbine blade material. CMCs can tolerate higher temperatures than alloyed metallic blades, resulting in better thermal efficiency of gas turbine engines. However, depending on the manufacturing process of the CMC, blades may have a larger trailing edge thickness. The design space therefore needs to be updated due to the resulting flow physics. Recently, experimental results acquired at NASA Glenn Transonic Turbine Blade Cascade Rig showed that a loss measure generally increased with increasing trailing edge thickness. For one case however, the losses initially decreased with increasing Reynolds number but rose and peaked at Re ~1.24x10 6 and then subsequently dropped to the original pre-peak trend. A possible cause of this intriguing phenomenon was speculated to be transonic vortex shedding, which is the mechanism of vortex shedding promoted by reflected shed pressure waves at the trailing edge at relatively high Reynolds numbers and higher Mach numbers. A Reynolds-averaged Navier-Stokes (RANS) analysis or a low-resolution large eddy simulation (LES) does not reproduce this apparent anomaly and thus it is worth performing a wall-resolved LES (the total mesh count of ~290 million cells) to investigate the aerodynamics of the CMC blade in the context of transonic vortex shedding. Our numerical results at Re=1,246,350 show that the pressure waves generated by the vortex shedding in the wake travel upstream and significantly influence the transition and separation on the suction side thus enhancing the said vortex shedding in the wake. This feedback seems not to hold under a low-Re condition (Re=621,900). The Reynolds number dependence was also examined by numerical perturbation of the pressure waves in the wake and by examining how such perturbation attenuates or endures. It is confirmed that the perturbation of the pressure waves is quickly damped.

LES↗

NASA NTRSDate not supplied

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.

Aerodynamics↗