Dear Jason.Jonkman,

Dear All,

I am trying to implement a Lidar system in feed forward CPC. Therefore, I’m trying to extract the wind speed that the turbine actually ‘feels’. Lets call it here rotor effective wind speed (rews). The way I currently do that is based on a simple drive train model (J*OmegaDot=Qa-Qg-Qloss) and leads at the end to the equation: Cp(lambda,beta)/lambda^3 = 2(Q_Lss + J*OmegaDot)/ rho*pi*R^5*Omega^2

To avoid issues with calculating the Cp, I perform a look up table on Cp/lambda^3. (I use a pre calculated Cp surface (attached as .zip), not done by me but seemed reasonable, at least at the beginning )

Since the needed turbine data is known (I extract Q_Lss =Qg+Qloss from strain gauges on the LSS), i go with the known pitch angle (from output channel) to the look-up, compare it with the right side of the equation and so, I find the corresponding TSR. From the TSR the rotor effective wind speed is derived: v_eff = Omega*R/TSR

As I tested this with uniform wind and wind steps (12 to 25 m/s with increment of 1m/s, see attached stepInputAndStepResponse.png) i found that there is a bias between the uniform wind speed and the rews. Now I am trying to find out what causes the bias. I’ve checked everything else twice so one of the possible options left is that I actually have a wrong Cp surface (look up table). Attached a graph with the Cp from FAST, Cp calculated with my lookup table and Cp calculated with my look up table but with the TSR input from FAST (RtTSR). So according to this I have a few questions.

Currently, only the gen DOF is enabled. AFAeroMod =1 (steady aerodynamics) Prandtl tip/hub losses are on. No tower shadow or potential flow.

Would you please tell me the correct configuration of settings to calculate relative precise Cp (TSR,beta) surface. Does the aerodynamics model matter? (read in the forum that Cp above .5 can be computed at steady aerodynamics)

I am also using the AeroDyn output channels (RtAeroCp, RtTSR, RtVAvgxh) as refference but could not find how are these calculated.

How exactly is the averagad x-wind speed calculated?

And at the end Do You think that the method I am using to extract the rews is correct at all or You would do it in a different way?

Thank You massively for Your support!!

Best

Filip

Coefficients5MW.zip (183 KB)