Turbine CFD Analysis
Turbine CFD: Theoretical Foundations
Overview
What's the difference between turbine CFD and compressor CFD?
Turbines are on the side that extracts energy from the fluid. Because the flow accelerates, large-scale separation like in compressors is less likely to occur. Instead, blade cooling, secondary flow losses, and transonic shock waves become the main challenges.
Stage Work and Isentropic Efficiency
How is turbine work expressed?
Output and efficiency are defined from the Euler equation.
$h_{02s}$ is the enthalpy after isentropic expansion. For modern design levels, $\eta_{is}=90\sim92\%$ for HP stages and $88\sim90\%$ for LP stages in gas turbines.
Blade Loading Coefficient
How do you evaluate the magnitude of blade loading?
The Zweifel blade loading coefficient is the standard.
$s$: pitch, $c_x$: axial chord. $Z_w \approx 0.8$ has been considered the traditional optimum, but recent high-loading designs are also researching $Z_w > 1.0$.
Software Selection
What software is used for turbine CFD?
Ansys CFX + TurboGrid is the most widely used among aero-engine manufacturers. NUMECA FINE/Turbo is efficient for setting up multi-stage turbines and is used by companies like Rolls-Royce. STAR-CCM+ has strengths in CHT (Conjugate Heat Transfer) analysis for turbine blade cooling.
The Mystery of the Betz Limit 59.3%—Why "100%" is Impossible
The Betz limit (16/27 ≈ 59.3%) derived by Albert Betz in 1919 is the theoretical upper limit of energy that a wind turbine can extract from the wind. Intuitively, one might think "just take all the energy," but doing so would make the downstream wind speed zero, stopping the flow and preventing new wind from entering. The key to maximum output is not completely reducing the wind speed but appropriately "letting it pass." The actual efficiency of modern large wind turbines is 45–50%, and excluding mechanical/electrical losses, they reach about 80–85% of the Betz limit. CFD continues to pursue improvements of a few percent through airfoil optimization.
Computational Methods for Turbine CFD
Importance of Blade Cooling
How is turbine blade cooling handled in CFD?
HP turbine inlet gas temperatures reach 1500–1800°C, far exceeding the heat resistance limit of blade materials (about 1000°C for Ni-based superalloys). Internal cooling passages and film cooling are used to lower the blade surface temperature.
Cooling Model Hierarchy
How do you incorporate cooling into CFD?
There are multiple levels depending on the trade-off between accuracy and cost.
| Level | Model | Computational Cost | Accuracy |
|---|---|---|---|
| L0 | No cooling flow (adiabatic wall) | Lowest | Baseline evaluation without cooling |
| L1 | Source Term (mass/energy injection) | Low | Rough estimate for film cooling |
| L2 | Discrete Hole (individual cooling hole BC) | Medium | Quantitative evaluation of film effectiveness |
| L3 | Resolved Cooling Holes (holes meshed) | High | Highest accuracy but high effort |
| L4 | CHT (fluid + solid conjugate) | Highest | Predicts internal blade temperature distribution |
Are L3 and L4 practical?
L3/L4 for a single blade is practical as a Singleton calculation. STAR-CCM+'s CHT is highly rated for this purpose. L3/L4 for multi-stage is currently at the research level.
Film Cooling Effectiveness
How do you evaluate the effectiveness of film cooling?
It is defined by the adiabatic film cooling effectiveness.
$T_g$: Mainstream gas temperature, $T_{aw}$: Adiabatic wall temperature, $T_c$: Cooling air temperature. $\eta_f = 0$ means no cooling, $\eta_f = 1$ means perfect cooling. In CFD, it is calculated by outputting the adiabatic wall temperature on the blade surface.
100 Years of the Actuator Disk Method—From BEM Theory to LES
The first mathematical model for wind turbine fluid analysis dates back to Betz's momentum theory in the 1920s. The actuator disk method, which treats the rotor as an "infinitely thin disk generating thrust," is still widely used today for wake analysis of entire wind farms (farm CFD). Since the 2000s, the actuator line model (ALM) has emerged, allowing the lift and drag of individual blades to be imposed as volume forces onto an LES flow field, enabling the reproduction of tip vortex generation and breakdown. The evolution story of theory starting from BEM and merging with LES over 100 years is a microcosm of CFD development history.
Turbine CFD in Practice
Turbine Blade Row Mesh
Is the mesh for a turbine blade row the same as for a compressor?
The basic structure is the same, but there are turbine-specific points to note.
- Trailing Edge Thickness: Turbine blades have very thin trailing edges (0.3–0.8mm). Sufficient cells are needed around the trailing edge in the O-grid.
- Cooling Holes: Local refinement around cooling holes is necessary for L2/L3 models.
- Transonic Regions on Blade Surface: Resolving the supersonic patch on the suction side and the trailing edge shock wave.
If the trailing edge is 0.3mm, the mesh must be quite fine, right?
The O-grid at the trailing edge should have at least 10 cells radially, and the wake region immediately behind the trailing edge should also have a fine mesh. TurboGrid's trailing edge cutoff function can control the trailing edge shape.
Transonic Turbine Blade Row
Does turbine flow become supersonic?
In HP turbines, the blade-to-blade Mach number reaches 1.1–1.3. After accelerating to supersonic speed on the suction side, an oblique shock wave is emitted from the trailing edge. Accurate prediction of this Trailing Edge Shock System, where the shock wave impinges on the adjacent blade, is key to CFD accuracy.
How much mesh is needed to resolve shock waves?
It is recommended that the cell size orthogonal to the shock wave direction be less than 0.5% of the chord, with at least 10 cells before and after the shock. Adaptive Mesh Refinement (AMR) to concentrate mesh at the shock location is also effective. AMR functions in Fluent or STAR-CCM+ can be used.
Performance Prediction Accuracy
What is the accuracy of turbine CFD?
| Metric | Accuracy |
|---|---|
| Stage efficiency (multi-stage) | ±0.5–1.5 points |
| Blade surface pressure distribution | Good (qualitatively matches experiment) |
| Blade surface heat transfer coefficient | ±10–20% (depends on turbulence model) |
| Trailing edge shock wave location | ±2% of chord |
Related Topics
Experience the theory firsthand with the interactive simulator for this field
All Simulators