CFD Analysis of Wind Turbines
Theory and Physics
Overview
Professor, how is CFD analysis used for wind turbine blades?
CFD for wind turbines is used for three main purposes: (1) Aerodynamic design of blades and power output prediction, (2) Wake (wake) analysis for wind farm layout optimization, and (3) Structural load evaluation under extreme wind conditions.
The theoretical upper limit for wind turbine power output is defined by the Betz limit. The key challenge in blade design is how close one can get to this limit.
Betz Limit and Power Coefficient
Wind turbine power coefficient:
Betz limit (theoretical maximum efficiency):
This is the theoretical upper limit where about 59.3% of the wind's energy can be extracted. The actual $C_P$ for large wind turbines (e.g., Vestas V164, Siemens Gamesa SG 14-222) is around 0.45--0.50, reaching about 80--85% of the Betz limit.
85% is quite close to the theoretical limit.
This is achieved through airfoil design, pitch control, and variable-speed operation optimization. Since there is little room for further improvement here, reducing wake losses across the entire wind farm has become the next focus.
Relationship Between BEM Theory and CFD
Blade Element Momentum theory (BEM) is the foundation of wind turbine analysis.
Here, $a$ is the axial induction factor, $a'$ is the tangential induction factor, and $\omega$ is the rotational angular velocity.
If BEM exists, why is CFD necessary?
BEM has limitations.
| Limitations of BEM | Advantages of CFD |
|---|---|
| Cannot account for 3D effects (root/tip vortices) | Solves 3D flow directly |
| Difficult to predict dynamic stall | Reproducible with unsteady CFD |
| Ignores interference between blades | Analyzes all blades simultaneously |
| Wake model is simplified | Directly calculates wake diffusion/merging |
| Ignores nacelle/tower interference | Reproduces tower shadow |
Tip Speed Ratio and Airfoils
Tip Speed Ratio (TSR):
For large wind turbines, the optimal $\lambda$ is approximately 6--9.
Wind turbine airfoils:
| Airfoil Series | Developer | Characteristics |
|---|---|---|
| NACA 63-xxx | NACA | Classical. Proven track record |
| DU (Delft) | Delft University of Technology | Thick airfoil. Used for root sections |
| FFA-W3 | FOI (Sweden) | High $C_{L,max}$. Insensitive to surface roughness |
| DTU-LN1xx | DTU (Denmark) | Latest design. CFD-optimized |
So the airfoil is different at the root and tip of the blade.
The root section uses thick airfoils (relative thickness 30--40%) for structural strength, while the tip section uses thin airfoils (relative thickness 18--24%) for aerodynamic performance. The airfoil shape changes continuously along the span.
The Reason Blades Are Twisted—Aerodynamic Basis of Pitch Angle
If you look closely at a wind turbine blade, you'll notice it twists from the root to the tip. This design is called "pitch distribution (twist angle)," and the reason is that the direction of the relative wind speed changes at each position. At the tip, the circumferential speed is high, resulting in a small angle of attack relative to the incoming wind, so the twist corrects the angle of attack. The basic task in blade design is to meticulously verify this optimal pitch distribution, derived from Betz theory, using CFD. Every "visually strange shape" has a proper fluid dynamic reason behind it.
Physical Meaning of Each Term
- Temporal term $\partial(\rho\phi)/\partial t$: Imagine turning on a faucet. At first, the water comes out unstable and splashing, but after a while, it becomes a steady flow, 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 this significantly reduces computational cost, starting with a steady-state solution 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 corner of a room is also because the "carrier," air, transports heat via convection. Here's the interesting part—this term contains "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 is 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 it naturally mixes. That's molecular diffusion. Next question—honey and 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 "sluggishly." In low Reynolds number flows (slow, viscous), diffusion dominates. Conversely, in high Re number flows, convection overwhelmingly dominates, and diffusion plays a minor role.
- Pressure term $-\nabla p$: When you push the plunger of a syringe, the liquid shoots out forcefully from the needle tip, right? Why? Because the piston 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. "Flow occurs where there is a pressure difference"—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 become strange immediately after switching to compressible analysis, it might be due to mixing up absolute/gauge pressure.
- Source term $S_\phi$: Warmed air rises—why? Because it becomes lighter (lower density) than its surroundings, so it is pushed upward by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat from a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... 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 heated room in winter.
Assumptions and Applicability Limits
- Continuum assumption: Valid for Knudsen number Kn < 0.01 (mean free path of molecules ≪ characteristic length)
- Newtonian fluid assumption: Shear stress and strain rate have a linear relationship (viscosity model required for non-Newtonian fluids)
- Incompressible 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 gas (Kn > 0.1), supersonic/hypersonic flow (shock wave 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 pressure 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 coefficient $\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
Analysis Method Hierarchy
What methods are used for wind turbine CFD?
There is a hierarchy of methods based on the trade-off between computational cost and accuracy.
| Method | Modeling | Cell Count | Application |
|---|---|---|---|
| BEM | 1D cross-section theory | -- | Initial design, annual power generation prediction |
| Actuator Disk (AD) | Represents rotor with volume forces | 1 million--10 million | Wind farm layout |
| ALM (Actuator Line) | Distributes volume forces along lines | 5 million--50 million | Wake analysis (LES) |
| Full Blade RANS | Directly solves 3D blade geometry | 10 million--50 million | Blade aerodynamic design |
| Full Blade LES | 3D blade + LES | 100 million--1 billion | Research purposes |
What is the Actuator Line Model (ALM)?
It's a method that does not physically model the blades but distributes volume forces equivalent to lift and drag along rotating lines. Since there's no need to resolve the blade boundary layer, the mesh count can be drastically reduced. It's a standard method for LES of wind farms.
Full Blade CFD Mesh
Full blade analysis for large wind turbines (rotor diameter ~200m class):
- Rotating domain: Cylinder 1.2 times the rotor diameter. Rotates using Sliding Mesh.
- Stationary domain: Outer boundary at least 10 times the rotor diameter.
- Blade surface: $y^+ < 1$, at least 20 prism layers.
- Tip vortex resolution: Refinement zone near the tip.
- Nacelle/Tower: Included in the same mesh (for tower shadow evaluation).
- Total cell count: 5 million--15 million cells per blade.
So for 3 blades + nacelle + tower, it becomes tens of millions of cells.
It's also possible to use a 1/3 model (1 blade + periodic boundary conditions) by exploiting rotational symmetry. However, evaluating tower shadow requires all three blades.
Turbulence Model
Turbulence model selection for wind turbine CFD:
Model Application Notes SST k-omega Blade steady-state aerodynamics Insufficient for dynamic stall $\gamma$-$Re_\theta$ + SST Transition prediction (blade root section) Transition on thick airfoils is important DDES (SST-based) Dynamic stall, tower shadow Unsteady calculation mandatory LES (ALM) Wind farm wake Atmospheric boundary layer turbulence generation required
Time Step and Rotation Handling
Time step for unsteady analysis of rotating blades:
$$ \Delta t = \frac{\Delta\theta}{\omega} $$
For a typical large turbine (rotational speed 12rpm = 0.2rps) at 1 degree per step:
$$ \Delta t = \frac{1°}{360° \times 0.2} = 0.0139 \text{ s} $$
360 steps per revolution. 10 revolutions would be 3600 steps.
The standard practice is to discard at least the first 5 revolutions to exclude the initial transient state, then take a time average over the subsequent 5--10 revolutions.
Coffee Break Casual Talk
Why "Time Step is Critical" in CFD for Rotating Blades
When setting up a rotating domain for wind turbine CFD, the time step setting is delicate. The angle the turbine rotates per step is typically kept below 1-2°, which for a rated speed of 15rpm results in a time step of about 0.01-0.02 seconds. If it's too coarse, the generation and transport of tip vortices become inaccurate, leading to large errors in power output prediction. "Saving on time steps makes you cry later" is a classic lesson in wind CFD, and in practice, it's a golden rule to perform a time step sensitivity test along with convergence checks at the outset.
Upwind Scheme
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
Second-order accurate, but numerical oscillations occur for Peclet 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 shock waves and steep gradients.
Finite Volume Method vs Finite Element Method
Related Topics
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
Turbulence model selection for wind turbine CFD:
| Model | Application | Notes |
|---|---|---|
| SST k-omega | Blade steady-state aerodynamics | Insufficient for dynamic stall |
| $\gamma$-$Re_\theta$ + SST | Transition prediction (blade root section) | Transition on thick airfoils is important |
| DDES (SST-based) | Dynamic stall, tower shadow | Unsteady calculation mandatory |
| LES (ALM) | Wind farm wake | Atmospheric boundary layer turbulence generation required |
Time step for unsteady analysis of rotating blades:
For a typical large turbine (rotational speed 12rpm = 0.2rps) at 1 degree per step:
360 steps per revolution. 10 revolutions would be 3600 steps.
The standard practice is to discard at least the first 5 revolutions to exclude the initial transient state, then take a time average over the subsequent 5--10 revolutions.
Why "Time Step is Critical" in CFD for Rotating Blades
When setting up a rotating domain for wind turbine CFD, the time step setting is delicate. The angle the turbine rotates per step is typically kept below 1-2°, which for a rated speed of 15rpm results in a time step of about 0.01-0.02 seconds. If it's too coarse, the generation and transport of tip vortices become inaccurate, leading to large errors in power output prediction. "Saving on time steps makes you cry later" is a classic lesson in wind CFD, and in practice, it's a golden rule to perform a time step sensitivity test along with convergence checks at the outset.
Upwind Scheme
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
Second-order accurate, but numerical oscillations occur for Peclet 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 shock waves and steep gradients.
Finite Volume Method vs Finite Element Method
Related Topics
なった
詳しく
報告