NREL / floris

A controls-oriented engineering wake model.
http://nrel.github.io/floris
BSD 3-Clause "New" or "Revised" License
198 stars 152 forks source link

[BUG] Wake added turbulence not working in TurbOPark model #909

Open misi9170 opened 1 month ago

misi9170 commented 1 month ago

This issue follows discussion #906 , where @prith-gs pointed out that the TurbOPark model is not correctly adding turbulence to the wake (even when running with the Crespo-Hernandez turbulence model). I believe the issue stems from this TODO item. Essentially, the area_overlap computed in the turbopark_solver is all zeros, so the subsequent lines that compute the wake added turbulence also produce zeros.

Steps to reproduce

import numpy as np
from floris import FlorisModel

fmodel = FlorisModel("inputs/turbopark.yaml")

fmodel.set(
    layout_x=[0, 1000],
    layout_y=[0, 0],
    wind_speeds=[9],
    wind_directions=[270],
)

fmodel.run()

print(fmodel.get_turbine_TIs())

produces

[[0.06 0.06]]

which is simply the ambient turbulence intensity at both turbines.

Expected output

According to the Crespo-Hernandez turbulence model, we'd expect that the turbulence intensity at the downstream turbine would be higher than the ambient turbulence intensity, i.e.

[[0.06 a]]

with a > 0.06.

Related

Bartdoekemeijer commented 1 month ago

I believe #907 may already solve this.