Fully-coupled CFD engine simulations

Authors: Mateus Teixeira, Luigi Romagnosi, Mohamed Mezine, Yannick Baux, Jan Anker, Kilian Claramunt, Charles Hirsch​

The aerospace industry, like many other industries, is under pressure to drastically reduce its environmental footprint. The Flightpath 2050 goals of the European Union state that by the year 2050 all CO2 emissions per passenger kilometer must be reduced by 75%, NOx by 90% and noise pollution by 65%. One of the ways to achieve these environmental goals is by increasing turbomachinery performance. Improved analysis in the design phase, more specifically development of more reliable predictions, advancement in accuracy, inter-disciplinarity and speed of simulation tools, can add several percentage points to engine efficiency and reduce development cost and time.

Complexity

One of the main challenges in designing an engine is the complexity in terms of geometrical details (combustion chamber features, turbine cooling holes…) and of interaction effects between the components which must be modelled with accuracy and acceptable computation time.

NUMECA’s approach to Full Engine CFD Simulation

A full 3D aerodynamic simulation of a complete gas turbine engine, applied to a micro turbine case has been conducted at NUMECA. The analysis was comprised of a single fully-coupled 3D CFD simulation for the flow of a KJ66 engine redesign. The injection and burning of fuel inside the combustion chamber were modeled with a simplified flamelet model. Using advanced RANS treatment with inputs from Nonlinear Harmonic (NLH) method (available as module for FINE™/Turbo), tangential non-uniformities were captured and the flow physics of the interaction between compressor, combustor and turbine were assessed.

Case description

The selected test case is a redesign version of the KJ66 micro gas turbine (Figure 2). The Figure 2 shows the layout of the redesigned version of the KJ66 micro gas turbine used in the full engine computation. The centrifugal compressor is mounted at the engine entrance and it is composed of an impeller and a bladed diffuser row. The combustion chamber is followed by a high-pressure turbine (HPT), which drives the compressor, and a low pressure turbine (LPT), which would drive a propeller in an independent shaft. An exhaust hood is connected to the last LPT row at the engine exit.

FIGURE 2 : Layout of the KJ66 micro gas turbine

Simulation setup

The computational domain encompassed one blade passage for each turbomachinery blade row, a 60° sector for the combustion chamber (containing one fuel injector) and half exhaust hood. The three-dimensional mesh for the blade rows (compressor and turbine rows) was generated automatically using Autogrid5™, NUMECA’s turbomachinery dedicated full automatic hexahedral block-structured grid generator. The mesh for the combustion chamber and exhaust hood was generated with HEXPRESS/Hybrid™, NUMECA’s unstructured hex-dominant conformal body-fitted mesher for arbitrary complex geometries. The entire mesh has 19.20 million points.

FIGURE 3: Compressor, combustion chamber and turbine 3D mesh view and blade-to-blade layout (shrouds are omitted).

As a first investigation, the steady RANS computation is performed with the Spalart-Allmaras turbulence model. Ambient total quantities are imposed at the engine inlet with specified axial velocity direction and static pressure is fixed at the outlet. The fuel injection is specified with static temperature and axial velocity of -120 m/s. Solid walls are assumed smooth and adiabatic. The convergence history is checked by following the evolution of the mass flow, the pressure ratio, and net torque between the HPT rotor and the impeller. The computation is launched in parallel on 144 processors on a computing cluster. In a second step, the computation is restarted with an improved rotor-stator connection based on the NLH method.

Combustion model

The injection and burning of fuel inside the combustion chamber are modeled with a simplified flamelet model implemented in OpenLabs™. With this model, an additional equation is solved for the mixture fraction f, with the flame temperature as function of the composition.
 

Results

Steady-state computation

The convergence history of mass flow rate error between the engine inlet and outlet drops to less than 0.2% after 25000 iterations in approximately 24 hours of wall clock time. The net torque Mz reaches a final value of +0.21 N.m when the couple impeller-HPT rotor spins at -80000 rpm, meaning that the turbine produces enough torque to drive the compressor and both components are very close to be load balanced.

Figure 4 shows the static pressure, static temperature and absolute Mach number distributions at midspan of the compressor and turbine as well as in the solid walls of the combustion chamber and exhaust hood. The flow fields are continuous across the machine with a gradual raise of Mach number through the impeller and subsequent conversion of the kinetic energy into pressure across the diffuser. At approximately 160 kPa, air reaches the combustion chamber where the simulated combustion process takes place with a relatively small pressure loss. The maximum temperature at the combustion chamber is around 2200K at the combustor inner chamber. The hot gases from the combustion enter the HPT at approximately 983K and are expanded through the downstream blade rows, exiting the machine at 932K.

NLH Computation

 

The results from the RANS computation using a mixing plane treatment at the rotor stator interface can be compared to the results of the NLH analysis. In particular, it is interesting to note the effect of the improved connection approach with respect to the inter-components interactions.

Figures 6 and 7 show similar results on the mass flux for both simulations at the exit of the combustion chamber. A pattern linked to the periodicity of the fuel pipes and to the combustion zones can be observed. The transfer of information from the combustor outlet to the adjacent downstream HPT nozzle is different for the RANS simulation and its NLH counterpart. The mixing plane approach shows no tangential non-uniformities, while the NLH rotor-stator treatment proves to be a low cost first step in capturing tangential non-uniformities.

 

FIGURE 6 : Mass flux at the outlet of the combustion chamber and at the inlet of the HPT for the steady-state (left) and NLH 0 harmonic (right) simulations.

 

FIGURE 7 : Mass flux at the outlet of the combustion chamber and at the inlet of the HPT for the steady-state (left) and NLH 0 harmonic (right) simulations.

 

The Figure 8 shows the results for absolute total enthalpy and mixture fraction at the HPT and LPT inlet connections. The azimuthal averaging used in the mixing plane approach, where tangential non-uniformities are not perceived by downstream components, can be noticed for the RANS simulation. In contrast, the Fourier decomposition together with the local non-reflective boundary treatment used in the NLH method is able to show a slight improvement at the interfaces.

 

FIGURE 8 : Absolute total enthalpy and mixture fraction at the inlet of the HPT and LPT blade rows for the steady-state (top) and NLH 0 harmonic (bottom) simulations.

 

Conclusions

 

A successful simulation of the three-dimensional flow of a complete micro gas turbine engine has been achieved, on the basis of a fully-coupled 3D CFD simulation of a redesigned version of the KJ66 micro gas turbine. Compared to the component-by-component analysis, the fully coupled approach enables the solution of the whole engine in one single simulation. Furthermore it simplifies the simulation workflow as only the engine inlet and outlet pressures, the rotating speed and the fuel mass flow rate need to be prescribed.

 

---

Source: Mateus Teixeira, Luigi Romagnosi, Mohamed Mezine, Yannick Baux, Jan Anker, Kilian Claramunt, Charles Hirsch A Methodology for Fully-Coupled CFD Engine Simulations, Applied to a Micro Gas Turbine Engine. ASME. Turbo Expo: Power for Land, Sea and Air, Volume 2C: Turbomachinery ():V02CT42A047. doi:10.1115/GT2018-76870.