Hello again,
After comparing your rans deflection model to our implementation at TUM, I observed differences in the far wake. The problem lies in equation 7.4 of the paper of Bastankhah and Porté-Agel from 2016.
In rans_deficit_deflection.m line 93 the equation is given:
FW_delta = obj.delta_x0+obj.theta_C0*(obj.rotorRadius/7.35)*sqrt(sqrt(det(obj.C/((diag([obj.ky obj.kz])*obj.ct)^2))))*(2.9+1.3*sqrt(1-obj.ct)-obj.ct)*lnTerm.';
I think it should be written
FW_delta = obj.delta_x0+obj.theta_C0*(obj.rotorRadius/7.35)*sqrt(sqrt(det(obj.C/((diag([obj.ky obj.kz])*sqrt(obj.ct))^2))))*(2.9+1.3*sqrt(1-obj.ct)-obj.ct)*lnTerm.';
that the code correctly represents eq 7.4 with sqrt(obj.ct) instead of just obj.ct.
Please check this equation.
Best,
Andi
Hello again,
After comparing your rans deflection model to our implementation at TUM, I observed differences in the far wake. The problem lies in equation 7.4 of the paper of Bastankhah and Porté-Agel from 2016.
In rans_deficit_deflection.m line 93 the equation is given:
FW_delta = obj.delta_x0+obj.theta_C0*(obj.rotorRadius/7.35)*sqrt(sqrt(det(obj.C/((diag([obj.ky obj.kz])*obj.ct)^2))))*(2.9+1.3*sqrt(1-obj.ct)-obj.ct)*lnTerm.';I think it should be written
FW_delta = obj.delta_x0+obj.theta_C0*(obj.rotorRadius/7.35)*sqrt(sqrt(det(obj.C/((diag([obj.ky obj.kz])*sqrt(obj.ct))^2))))*(2.9+1.3*sqrt(1-obj.ct)-obj.ct)*lnTerm.';that the code correctly represents eq 7.4 with sqrt(obj.ct) instead of just obj.ct.
Please check this equation.
Best,
Andi