Skip to content
Snippets Groups Projects
Commit cda2ab9e authored by Mikkel Friis-Møller's avatar Mikkel Friis-Møller Committed by Mads M. Pedersen
Browse files

Tilt defl

parent 8ac1703d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
py_wake/deflection_models/jimenez.py
+ 2
− 2
View file @ cda2ab9e
......@@ -25,8 +25,8 @@ class JimenezWakeDeflection(DeflectionModel):
nominator_ijxl = (1 + (self.beta / D_src_il)[:, na, na, :] * dw_ijxl)**2
alpha = denominator_ilk[:, na, na] / nominator_ijxl[..., na]
deflection_ijlk = np.trapz(np.sin(alpha), dw_ijxl[..., na], axis=2)
self.hcw_ijlk = hcw_ijl[..., na] + deflection_ijlk * np.cos(theta_deflection_ilk[..., na])
self.dh_ijlk = dh_ijl[..., na] + deflection_ijlk * np.sin(theta_deflection_ilk[..., na])
self.hcw_ijlk = hcw_ijl[..., na] + deflection_ijlk * np.cos(theta_deflection_ilk[:, na])
self.dh_ijlk = dh_ijl[..., na] + deflection_ijlk * np.sin(theta_deflection_ilk[:, na])
return dw_ijl[..., na], self.hcw_ijlk, self.dh_ijlk
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment