Dear Jason,
I tried to implement a BEM code in MATLAB, using the Distributed Blade Aerodynamic Properties of the 5MW NREL wind turbine (Table 3-1 in the 5MW NREL reference pdf: nrel.gov/docs/fy09osti/38060.pdf) and the airfoils data. I find the maximum power coefficient of Cp_max=0.486 at TSR=7.82 and Pitch=0.0°, which are slightly off what is in the pdf reference (Cp_max=0.482 and TSR=7.55). Here is a plot of Cp with TSR and Pitch angle.
I have two questions:
_ Is this tiny difference comes from the number of blade nodes taken?
_ Are those results that I found make sense to you?