Skip to content

Fixing bug occurring with combination of time series with turbulence intensity...

Mikkel Friis-Møller requested to merge no_wake_aep into master

Fixing bug occurring with combination of time series with turbulence intensity and no_wake calculation, as reproduced by:

from py_wake import NOJ
import numpy as np
from py_wake.examples.data.iea37 import IEA37Site, IEA37_WindTurbines
from py_wake.examples.data import example_data_path
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
d = np.load(example_data_path + "/time_series.npz")
n_days=30
wd, ws, ws_std = [d[k][:6*24*n_days] for k in ['wd', 'ws', 'ws_std']]
ti = np.minimum(ws_std/ws,.5)
time_stamp = np.arange(len(wd))/6/24
wf_model = NOJ(site, windTurbines, k = 0.04)
sim_res_time = wf_model(x, y, # wind turbine positions
                        wd=wd, # Wind direction time series
                        ws=ws, # Wind speed time series
                        time=time_stamp, # time stamps
                        TI=ti, # turbulence intensity time series
                  )
total_AEP_nowake_time = sim_res_time.aep(with_wake_loss=False).sum().data```
Edited by Mikkel Friis-Møller

Merge request reports