Film Cooling
Theory and Physics
Film Cooling Principle
Professor, how exactly does film cooling work to cool the blade?
Cooling air is ejected from inside the turbine blade onto the blade surface, forming a low-temperature "film" between it and the high-temperature main gas flow. This film acts as an insulating layer, reducing the heat flux to the blade surface. It's an essential technology in modern gas turbines where inlet temperatures exceed 1500°C.
Is there an index to represent cooling efficiency?
The adiabatic film cooling effectiveness $\eta$ is the fundamental index.
Here $T_{\infty}$ is the mainstream temperature, $T_{aw}$ is the adiabatic wall temperature, and $T_c$ is the cooling air temperature. $\eta = 1$ means the wall is at the cooling air temperature, and $\eta = 0$ means no cooling effect.
Blowing Ratio and Momentum Ratio
Is it better to blow more cooling air?
It's not that simple. The blowing ratio $M$ is an important parameter.
If $M$ is too small, the cooling air gets swallowed by the mainstream and the effect is weak. If $M$ is too large, the cooling jet lifts off from the wall, actually reducing cooling efficiency. For circular holes, the optimal range is said to be around $M \approx 0.5$ to $1.0$.
That depends on the hole shape too, right?
Of course. Shaped holes (fan-shaped, laidback holes) have an expanded exit, which reduces the momentum of the cooling air, making lift-off less likely. Data shows that for the same $M$, $\eta$ can improve by 30-50% compared to circular holes. Engine manufacturers like GE and Rolls-Royce have even patented their own hole shapes.
The Birth of Film Cooling—The Dawn of Jet Engines in the 1940s
The concept of film cooling dates back to the late 1940s, when Rolls-Royce in the UK attempted to eject cooling air from a perforated plate onto the combustor liner of the Derwent engine. Initially called "Air Film Cooling," the mechanism relied on empirical rules. The theoretical foundation was established with the formulation of the cooling effectiveness equation (η = (T∞ - Taw)/(T∞ - Tc)) by Goldstein (1971), and over the next 50 years, it evolved to the elucidation of vortex structures around cooling holes using LES (Large Eddy Simulation). Modern jet engine turbine inlet temperatures exceed 1700°C, surpassing the melting point of nickel superalloys (about 1350°C), making film cooling indispensable for even a second of operation.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out in an unstable, spluttering manner, but after a while, the flow becomes steady, right? This "period of change" is described by the temporal term. The pulsation of blood flow from a heartbeat, or the flow fluctuation each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning this term is set to zero. Since computational cost drops significantly, solving first in steady-state is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far end of a room is also because the "carrier," air, transports heat via convection. Here's the interesting part—this term includes "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They are completely different! Convection is carried by flow, conduction is transmitted by molecules. There's an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while, they naturally mix. That's molecular diffusion. Now, a question—honey or water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move in a "thick" manner. In low Reynolds number flows (slow, viscous), diffusion dominates. Conversely, in high Re number flows, convection overwhelmingly dominates, and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push the plunger of a syringe, liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are densely packed? That's right, strong winds blow. "Where there is a pressure difference, flow is generated"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results go wrong immediately after switching to compressible analysis, it might be due to mixing up absolute/gauge pressure.
- Source Term $S_\phi$: Heated air rises—why? Because it becomes lighter (lower density) than its surroundings, so it's pushed up by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat generated by a gas stove flame, Lorentz force applied to molten metal by an electromagnetic pump in a factory... These are all actions that "inject energy or force into the fluid from the outside," expressed by the source term. What happens if you forget the source term? In natural convection analysis, forgetting buoyancy means the fluid doesn't move at all—a physically impossible result where warm air doesn't rise in a room with the heater on in winter.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (molecular mean free path ≪ characteristic length)
- Newtonian Fluid Assumption: Shear stress and strain rate have a linear relationship (non-Newtonian fluids require viscosity models)
- Incompressibility Assumption (for Ma < 0.3): Density is treated as constant. For Mach numbers above 0.3, compressibility effects must be considered.
- Boussinesq Approximation (Natural Convection): Density variation is considered only in the buoyancy term; constant density is used in other terms.
- Non-applicable Cases: Rarefied gases (Kn > 0.1), supersonic/hypersonic flow (shock capturing required), free surface flow (VOF/Level Set, etc., required)
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units. |
| Pressure $p$ | Pa | Distinguish between gauge and absolute pressure. Use absolute pressure for compressible analysis. |
| Density $\rho$ | kg/m³ | Air: approx. 1.225 kg/m³ @20°C, Water: approx. 998 kg/m³ @20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity $\nu = \mu/\rho$ [m²/s] |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Indicator for laminar/turbulent transition. |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time-step stability. |
Numerical Methods and Implementation
RANS vs LES Selection
Is RANS sufficient for film cooling CFD?
RANS is the mainstream approach but has limitations. The standard k-ε model tends to overestimate lateral diffusion of the cooling jet, often overpredicting $\eta$. Realizable k-ε or SST k-ω models are recommended, but even they struggle to accurately capture jet reattachment after lift-off.
Does LES work well then?
LES can physically and correctly reproduce the mixing of the cooling jet and mainstream, so the $\eta$ distribution is closer to experiments than RANS. However, computational cost becomes 100 to 1000 times that of RANS. In practice, a realistic approach is to use RANS in the design stage and LES for final verification. Recently, Hybrid RANS/LES (DES or SBES models) are used as a compromise.
Mesh Strategy
How should I mesh around the cooling holes?
Use the cooling hole diameter $D$ as a reference. Aim for at least 20 divisions around the hole circumference and about 20 divisions along the hole's internal length. The first wall layer should target $y^+ < 1$. Ideally, maintain a mesh density of $\Delta x / D \approx 0.1$ downstream of the hole exit for at least $x/D = 30$.
If there are hundreds of holes across the entire blade surface, meshing them all is unrealistic, right?
Exactly. In practice, there are three approaches. (1) Mesh only representative holes in detail and simulate others with mass flow boundary conditions. (2) Use a source term model to introduce cooling effects without individually meshing holes. (3) Use dedicated features like the Film Cooling Model in Ansys Fluent. STAR-CCM+ also has a Cooling Film function.
What about OpenFOAM?
OpenFOAM doesn't have a standard dedicated film cooling function, but you can use codedFixedValue to prescribe velocity/temperature profiles at hole exits. Another method is to connect a separate region with fine mesh only around the holes using AMI (Arbitrary Mesh Interface).
Coffee Break Trivia
Film Cooling Hole Mesh—"Insufficient Resolution" is the Biggest Enemy
The most computationally expensive aspect of film cooling CFD is the difficulty of "sufficiently resolving the interface between the cooling hole interior and the mainstream." Hole diameters are typically very small, around 0.5–2 mm, but the turbulent boundary layer inside the hole and the mixing region of the jet and mainstream at the hole exit require a resolution of at least 10–20 divisions relative to the hole diameter. Actual blades have hundreds to thousands of cooling holes, so modeling all holes with high resolution can easily exceed 100 million mesh cells. Therefore, in practice, a two-scale analysis approach is adopted: "detailed analysis of a few representative holes to understand local characteristics, and solving the entire blade with modeled boundary conditions." LES offers significantly higher accuracy for mixing near the hole exit, but computational cost is tens of times higher than RANS.
Upwind Scheme (Upwind)
First-order upwind: Large numerical diffusion but stable. Second-order upwind: Improved accuracy but risk of oscillations. Essential for high Reynolds number flows.
Central Differencing (Central Differencing)
Second-order accurate, but numerical oscillations occur for Pe number > 2. Suitable for low Reynolds number diffusion-dominated flows.
TVD Schemes (MUSCL, QUICK, etc.)
Maintain high accuracy while suppressing numerical oscillations via limiter functions. Effective for capturing shocks and steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD. FEM: Advantageous for complex shapes and multiphysics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit methods: CFL ≤ 1 is the stability condition. Implicit methods: Stable even for CFL > 1, but affects accuracy and iteration count. LES: CFL ≈ 1 is recommended. Physical meaning: Information should not travel more than one cell per time step.
Residual Monitoring
Convergence is typically judged when residuals for the continuity equation, momentum, and energy each drop by 3–4 orders of magnitude. The mass conservation residual is particularly important.
Relaxation Factors
Typical initial values: Pressure: 0.2–0.3, Velocity: 0.5–0.7. Reduce factors if diverging. Increase after convergence to accelerate.
Internal Iterations for Unsteady Calculations
Iterate within each time step until a steady solution converges. Internal iteration count: 5–20 iterations is a guideline. If residuals fluctuate between time steps, reconsider the time step size.
Analogy for the SIMPLE Method
The SIMPLE method is an "alternating adjustment" technique. First, velocity is tentatively determined (predictor step), then pressure is corrected so that mass conservation is satisfied with that velocity (corrector step), and velocity is revised using the corrected pressure—this back-and-forth is repeated to approach the correct solution. It resembles two people leveling a shelf: one adjusts the height, the other balances it, and they repeat this alternately.
Analogy for the Upwind Scheme
The upwind scheme is a method that "stands in the river flow and prioritizes upstream information." A person in the river looking downstream cannot tell where the water comes from—it reflects the physics that upstream information determines downstream conditions. Although it's first-order accurate, it is highly stable because it correctly captures flow direction.
Is RANS sufficient for film cooling CFD?
RANS is the mainstream approach but has limitations. The standard k-ε model tends to overestimate lateral diffusion of the cooling jet, often overpredicting $\eta$. Realizable k-ε or SST k-ω models are recommended, but even they struggle to accurately capture jet reattachment after lift-off.
Does LES work well then?
LES can physically and correctly reproduce the mixing of the cooling jet and mainstream, so the $\eta$ distribution is closer to experiments than RANS. However, computational cost becomes 100 to 1000 times that of RANS. In practice, a realistic approach is to use RANS in the design stage and LES for final verification. Recently, Hybrid RANS/LES (DES or SBES models) are used as a compromise.
How should I mesh around the cooling holes?
Use the cooling hole diameter $D$ as a reference. Aim for at least 20 divisions around the hole circumference and about 20 divisions along the hole's internal length. The first wall layer should target $y^+ < 1$. Ideally, maintain a mesh density of $\Delta x / D \approx 0.1$ downstream of the hole exit for at least $x/D = 30$.
If there are hundreds of holes across the entire blade surface, meshing them all is unrealistic, right?
Exactly. In practice, there are three approaches. (1) Mesh only representative holes in detail and simulate others with mass flow boundary conditions. (2) Use a source term model to introduce cooling effects without individually meshing holes. (3) Use dedicated features like the Film Cooling Model in Ansys Fluent. STAR-CCM+ also has a Cooling Film function.
What about OpenFOAM?
OpenFOAM doesn't have a standard dedicated film cooling function, but you can use codedFixedValue to prescribe velocity/temperature profiles at hole exits. Another method is to connect a separate region with fine mesh only around the holes using AMI (Arbitrary Mesh Interface).
Film Cooling Hole Mesh—"Insufficient Resolution" is the Biggest Enemy
The most computationally expensive aspect of film cooling CFD is the difficulty of "sufficiently resolving the interface between the cooling hole interior and the mainstream." Hole diameters are typically very small, around 0.5–2 mm, but the turbulent boundary layer inside the hole and the mixing region of the jet and mainstream at the hole exit require a resolution of at least 10–20 divisions relative to the hole diameter. Actual blades have hundreds to thousands of cooling holes, so modeling all holes with high resolution can easily exceed 100 million mesh cells. Therefore, in practice, a two-scale analysis approach is adopted: "detailed analysis of a few representative holes to understand local characteristics, and solving the entire blade with modeled boundary conditions." LES offers significantly higher accuracy for mixing near the hole exit, but computational cost is tens of times higher than RANS.
Upwind Scheme (Upwind)
First-order upwind: Large numerical diffusion but stable. Second-order upwind: Improved accuracy but risk of oscillations. Essential for high Reynolds number flows.
Central Differencing (Central Differencing)
Second-order accurate, but numerical oscillations occur for Pe number > 2. Suitable for low Reynolds number diffusion-dominated flows.
TVD Schemes (MUSCL, QUICK, etc.)
Maintain high accuracy while suppressing numerical oscillations via limiter functions. Effective for capturing shocks and steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD. FEM: Advantageous for complex shapes and multiphysics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit methods: CFL ≤ 1 is the stability condition. Implicit methods: Stable even for CFL > 1, but affects accuracy and iteration count. LES: CFL ≈ 1 is recommended. Physical meaning: Information should not travel more than one cell per time step.
Residual Monitoring
Convergence is typically judged when residuals for the continuity equation, momentum, and energy each drop by 3–4 orders of magnitude. The mass conservation residual is particularly important.
Relaxation Factors
Typical initial values: Pressure: 0.2–0.3, Velocity: 0.5–0.7. Reduce factors if diverging. Increase after convergence to accelerate.
Internal Iterations for Unsteady Calculations
Iterate within each time step until a steady solution converges. Internal iteration count: 5–20 iterations is a guideline. If residuals fluctuate between time steps, reconsider the time step size.
Analogy for the SIMPLE Method
The SIMPLE method is an "alternating adjustment" technique. First, velocity is tentatively determined (predictor step), then pressure is corrected so that mass conservation is satisfied with that velocity (corrector step), and velocity is revised using the corrected pressure—this back-and-forth is repeated to approach the correct solution. It resembles two people leveling a shelf: one adjusts the height, the other balances it, and they repeat this alternately.
Analogy for the Upwind Scheme
The upwind scheme is a method that "stands in the river flow and prioritizes upstream information." A person in the river looking downstream cannot tell where the water comes from—it reflects the physics that upstream information determines downstream conditions. Although it's first-order accurate, it is highly stable because it correctly captures flow direction.
Practical Guide
Comparison with Experimental Data
What experimental data should I use to validate film cooling CFD?
Classics include Goldstein et al.'s single circular hole film cooling experiment (1968) and, more recently, shaped hole experimental data from the University of Michigan and Ohio State University. The standard validation procedure is to compare with 2D $\eta$ distributions obtained using temperature-sensitive paint (PSP/TSP) or IR thermography.
How much discrepancy between CFD and experiment is acceptable?
For RANS, $\eta$ on the centerline within 20% of experimental values, and laterally averaged $\eta$ within 15%, is often considered "reasonable." For LES, aim for within 10%. However, discrepancies tend to be larger in lift-off regions and at jet edges.
Design Parameters and Sensitivity
What are the important parameters in design?
Let's organize the main parameters and their effects.
| Parameter | Typical Range | Effect on $\eta$ |
|---|---|---|
| Blowing Ratio $M$ | 0.3–2.0 | Maximum around $M \approx 0.5$–$1.0$ |
| Hole Spacing Ratio $P/D$ | 3–6 | Smaller improves lateral coverage |
| Injection Angle $\alpha$ | 20–45 degrees | Smaller improves wall attachment |
| Hole Shape | Circular, Fan-shaped, Layback | Significant improvement with Shaped holes |
| Density Ratio $DR$ | 1.0–2.0 | In real engines $DR \approx 2$, simulated with $CO_2$ in experiments |
An injection angle of 20 degrees is quite shallow. Can it be manufactured?
Related Topics
なった
詳しく
報告