Aerodynamics · CFD + FEA + Experiment
Biomimetic winglet retrofit — NREL 5-MW wind turbine blade
BEng dissertation, Liverpool John Moores University · project mark 87% · supervisor: Russell English
- Solver
- ANSYS Fluent 2024 R2
- Model
- k-ω SST + low-Re
- Mesh
- 2.1M cells
- Wall
- y+ ≈ 1
- Scheme
- Coupled · 2nd-order
- Rig
- GUNT ET210
A wind turbine blade sheds a strong vortex off its tip, and the drag that comes with it costs you energy you would otherwise capture. Birds have the same problem. The Andean condor deals with it by spreading its primary feathers so that one strong vortex becomes several weak ones. I wanted to know whether that still works on a 61 metre blade, so I tested it three separate ways: structurally, computationally, and on real hardware.
Design
The blade came first, and rebuilding it properly took longer than the winglets did. The NREL 5-MW definition gives you seventeen radial stations, each with its own chord length, twist angle and aerofoil profile. I placed an aerofoil section at every station, scaled each one to its local chord, rotated it to its prescribed twist, and lofted a surface through the set. The profiles are not all the same: three cylindrical root sections, then DU40 at node 4 running through DU35, DU30, DU25 and DU21 toward mid-span, then NACA 64224 from node 12 out to the tip at 61.63 m. Getting section-to-section continuity right at the profile transitions mattered more than anything else, because a discontinuity there produces numerical artefacts in both the FEA and the CFD rather than an obviously broken model.
Chord falls from root to tip to trade structural stiffness against drag, and twist is set so each section sits near its best angle of attack across the operating range. I implemented both distributions straight from the published data without adjustment, so the baseline keeps its original aerodynamic behaviour and any change in the results belongs to the winglet rather than to me.
On top of that I built three winglets differing only in sweep angle: 10°, 25° and 40°. Height stayed at 1 m throughout, roughly 1.6% of blade length and in line with published full-scale practice. Cant angle, aerofoil profile and junction geometry were fixed across all four configurations, which is what lets the results be attributed to sweep and nothing else.
Structural check first
There is no point improving the aerodynamics if the blade cannot carry the load, so I checked that first. All four configurations were run in SolidWorks Simulation at an extreme aerodynamic load of 36,569 Pa using GFRP properties. Peak von Mises stress came out at the mid-span chord transition rather than at the root or the winglet junction, which makes sense once you consider that this is where a high bending moment meets the widest section.
CFD methodology
The domain is a 150 m hemispherical inlet with a 250 m wake extension and a body of influence zone for refinement. Ten inflation layers bring y+ down to roughly 1, which the low-Reynolds k-ω SST formulation needs, and an 18° curvature setting puts cells where they are needed at the leading and trailing edges without me placing them by hand. My first attempt used 15 inflation layers and produced non-positive volume cells at the winglet junction; dropping to 10 fixed it in every case. For mesh independence I ran three refinement levels and saw under 2% variation in CL and CD between the medium and fine meshes.
I used two convergence criteria rather than one: scaled residuals down to 10⁻⁶ on velocity, and lift and drag monitors holding within 0.1% over the last 50 iterations. Residuals flattening out on their own is not enough, because the forces can still be moving when they do.
Results
Select winglet sweep angle
1.4736
CL
0.6452
CD
2.284
CL / CD
—
Δ efficiency vs baseline
Reported CFD coefficients, steady-state RANS. Values marked ~ are derived from reported CL/CD and CD.
What the winglet actually does
Q-criterion iso-surfaces show what is going on. On the baseline blade there is one tightly wound helical vortex tube. Add the swept winglet and it splits in two: a primary core at the blade-winglet junction and a secondary one at the winglet tip. Two weaker structures instead of one strong one, which is the same thing the condor's feathers are doing. Drag the handle to compare.
Testing it on the rig
CFD results are worth more once something physical agrees with them. I designed three winglet sets for printing, all NACA 0015, 12 mm tall, with sweep offsets of 2.1, 5.6 and 10.1 mm, then FDM printed them and bonded them to a GUNT ET210 rig blade. Print orientation mattered because I wanted the layer lines off the aerodynamic surfaces, and the bonded root joint had to hold up to 3,000 rpm without letting the incidence shift, so this was as much a manufacturing problem as an aerodynamic one. I took 16 steady state points per configuration between 200 and 3,000 rpm at 10 m/s and repeated the runs. Variation between repeats was small compared with the differences between configurations.
Measured improvement vs baseline — 16-node steady-state test, 10 m/s
Mean electrical powerPeak electrical powerBars scaled to the largest gain
| Config | ΣP (W) | Mean P (W) | Peak P (W) | Mean Δ | Peak Δ |
|---|---|---|---|---|---|
| Baseline | 59.565 | 3.723 | 8.426 | — | — |
| 10° winglet | 64.317 | 4.020 | 9.094 | +8.0% | +7.9% |
| 25° winglet | 65.509 | 4.094 | 9.163 | +10.0% | +8.7% |
| 40° winglet | 65.595 | 4.100 | 9.211 | +10.1% | +9.3% |