I would like to start by apologizing for my ignorance, I am sure I am missing something very obvious here.
I am using data for the NREL 5MW reference turbine, and specifically power curve and power coefficient curve. I know that the relationship between the power coefficient (Cp) and the blade pitch angle is not so obvious, but I thought I could use the general approximation formulas:
Where c1, c2, c3, etc… and x are constants, lambda is the tip-speed ratio and beta the blade pitch angle. Since I have all the values for c1, c2, c3, x, etc…, the tip speed ratio and the power coefficient from the NREL report, the only approach I thought about was:
- For each wind speed, get the corresponding “real” tip-speed ratio and Cp from the NREL curve.
- For each wind speed, find beta that minimizes the square difference between the “real” Cp and the one calculated using the equations above
The minimization routine works perfectly, the Cp (from rated wind speed to cut-out) is reproduced to the 8th decimal at all wind speeds after rated (my minimization is labelled as “Mine”):
However, plotting the calculated blade pitch angles against the ones reported by NREL I get this:
You can see that my calculated beta (labelled as “Mine”) diverges almost straight away from the official values, ending up in some ridiculous pitch angle at high wind speed.
Could you please comment on whether the approach appears sensible? Is there maybe a more intelligent (or easier) way to tackle this problem - maybe with other formulas I know nothing about?
Thank you in advance for any suggestions.