def test_YZGrid_terrain_parallel():
site = ParqueFicticioSite(distance=StraightDistance())
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x, y, wd=0, ws=10)
x = 264000
fm = simulation_result.flow_map(
grid=YZGrid(x, z=110, resolution=20, extend=0))
y = fm.X[0]
x = np.zeros_like(y) + x
z = site.elevation(x, y)
simulation_result.flow_map(XYGrid(extend=0.005)).plot_wake_map()
plt.plot(x, y, '.')
plt.figure()
simulation_result.flow_map(grid=YZGrid(
fm.x.item(), fm.y, z=np.arange(30, 210, 10))).plot_wake_map()
plt.plot(y, z + 110, '.')
plt.plot(y, fm.WS_eff_xylk[:, 0, 0, 0] * 100, label="ws*100")
plt.legend()
if 0:
print(np.round(fm.WS_eff_xylk[:, 0, 0, 0], 2).values.tolist())
plt.show()
plt.close('all')
npt.assert_array_almost_equal(fm.WS_eff_xylk[:, 0, 0, 0], [
4.55, 3.89, 3.15, 2.31, 4.41, 4.3, 7.41, 7.2, 7.03, 7.23, 7.32, 6.77,
6.83, 5.58, 11.01, 11.51, 11.93, 12.17, 11.03, 9.89
], 2)
def test_YZGrid_terrain_perpendicular():
site = ParqueFicticioSite(distance=StraightDistance())
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x, y, wd=270, ws=10)
x = x.max() + 10
fm = simulation_result.flow_map(
grid=YZGrid(x, z=110, resolution=20, extend=0))
y = fm.X[0]
x = np.zeros_like(y) + x
z = site.elevation(x, y)
simulation_result.flow_map(XYGrid(extend=.005)).plot_wake_map()
plt.plot(x, y, '.')
plt.figure()
simulation_result.flow_map(grid=YZGrid(
fm.x.item(), y=fm.y, z=np.arange(30, 210, 10))).plot_wake_map()
plt.plot(y, z + 110, '.')
plt.plot(y, fm.WS_eff_xylk[:, 0, 0, 0] * 100, label="ws*100")
plt.legend()
if 0:
print(np.round(fm.WS_eff_xylk[:, 0, 0, 0], 2).values.tolist())
plt.show()
plt.close('all')
npt.assert_array_almost_equal(fm.WS_eff_xylk[:, 0, 0, 0], [
5.39, 8.48, 8.42, 6.42, 5.55, 11.02, 4.99, 11.47, 5.32, 10.22, 13.39,
8.79, 8.51, 12.4, 5.47, 10.78, 10.12, 6.54, 10.91, 7.18
], 2)
def test_min_ws_eff_line():
site = IEA37Site(16)
x, y = [0, 600, 1200], [0, 0, 0] # site.initial_position[:2].T
windTurbines = IEA37_WindTurbines()
wfm = IEA37SimpleBastankhahGaussian(
site, windTurbines, deflectionModel=JimenezWakeDeflection())
yaw_ilk = np.reshape([-30, 30, 0], (3, 1, 1))
plt.figure(figsize=(14, 3))
fm = wfm(x, y, yaw=yaw_ilk, wd=270, ws=10).flow_map(
XYGrid(x=np.arange(-100, 2000, 10), y=np.arange(-500, 500, 10)))
min_ws_line = fm.min_WS_eff()
if 0:
fm.plot_wake_map()
min_ws_line.plot()
print(np.round(min_ws_line[::10], 2))
plt.show()
npt.assert_array_almost_equal(min_ws_line[::10], [
np.nan, np.nan, 11.6, 21.64, 30.42, 38.17, 45.09, 51.27, -8.65, -18.66,
-27.51, -35.37, -42.38, -48.58, -1.09, -1.34, -1.59, -1.83, -2.07,
-2.31, -2.56
], 2)
def test_aep_map():
site = IEA37Site(16)
x = [0, 0]
y = [0, 200]
windTurbines = IEA37_WindTurbines()
wake_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
aep = AEPCalculator(wake_model)
# print(aep.calculate_AEP([0], [0]).sum())
x_j = np.arange(-150, 150, 20)
y_j = np.arange(-250, 250, 20)
X, Y, Z = aep.aep_map(x_j, y_j, 0, x, y, wd=[0], ws=np.arange(3, 25))
m = 17
if 0:
import matplotlib.pyplot as plt
c = plt.contourf(X, Y, Z, 100) # , np.arange(2, 10, .01))
plt.colorbar(c)
windTurbines.plot(x, y)
plt.plot(X[m], Y[m], '.-r')
plt.show()
# print(np.round(Z[m], 2).tolist()) # ref
ref = [
21.5, 21.4, 21.02, 20.34, 18.95, 16.54, 13.17, 10.17, 10.17, 13.17,
16.54, 18.95, 20.34, 21.02, 21.4
]
npt.assert_array_almost_equal(Z[m].squeeze('h'), ref, 2)
def test_FlowBox():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wf_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
sim_res = wf_model(x, y)
flow_box = sim_res.flow_box(x=np.arange(0, 100, 10), y=np.arange(0, 100, 10), h=np.arange(0, 100, 10))
def test_power_xylk():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
# NOJ wake model
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x, y)
fm = simulation_result.flow_map(grid=HorizontalGrid(resolution=3))
npt.assert_array_almost_equal(fm.power_xylk(with_wake_loss=False)[:, :, 0, 0] * 1e-6, 3.35)
def test_not_implemented_plane():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wf_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
sim_res = wf_model(x, y)
grid = YZGrid(x=-100, y=None, resolution=100, extend=.1)
grid = grid(x, y, windTurbines.hub_height(x * 0), windTurbines.hub_height(x * 0))
with pytest.raises(NotImplementedError):
sim_res.flow_map(grid=grid, wd=270, ws=None).plot_wake_map()
def test_aep_map_no_turbines():
site = IEA37Site(16)
windTurbines = IEA37_WindTurbines()
wake_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
aep = AEPCalculator(wake_model)
x_j = np.arange(-150, 150, 20)
y_j = np.arange(-250, 250, 20)
_, _, Z = aep.aep_map(x_j, y_j, 0, [], [], wd=[0])
expect = (3.35 * 1e6) * 24 * 365 * (1e-9)
npt.assert_array_almost_equal(Z, expect, 2)
def test_YZGrid_plot_wake_map():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wf_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
sim_res = wf_model(x, y)
sim_res.flow_map(grid=YZGrid(x=-100, y=None, resolution=100, extend=.1), wd=270, ws=None).plot_wake_map()
if 0:
plt.show()
plt.close()
def test_Points():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wf_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
sim_res = wf_model(x, y)
flow_map = sim_res.flow_map(Points(x, y,
x * 0 + windTurbines.hub_height()),
wd=0,
ws=None)
if 0:
flow_map.WS_eff.plot()
plt.show()
plt.close('all')
def test_YZGrid():
site = IEA37Site(16)
x, y = site.initial_position.T
m = x < -1000
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x[m], y[m], wd=270)
fm = simulation_result.flow_map(grid=YZGrid(-1000, z=110, resolution=20))
if 0:
simulation_result.flow_map(grid=YZGrid(-1000)).plot_wake_map()
plt.plot(fm.X[0], fm.Y[0], '.')
print(np.round(fm.WS_eff_xylk[0, :, 0, 0], 2).data.tolist())
plt.show()
npt.assert_array_almost_equal(fm.WS_eff_xylk[:, 0, 0, 0],
[9.8, 9.8, 8.42, 5.24, 9.74, 9.8, 9.8, 9.8, 9.76, 7.61, 7.61,
9.76, 9.8, 9.8, 9.8, 9.74, 5.24, 8.42, 9.8, 9.8], 2)
def test_IEA37SimpleBastankhahGaussian_all_ex():
"""check SBG matches IEA37 value in file"""
for n_wt in [9, 16, 36, 64]:
x, y, aep_ref = read_iea37_windfarm(iea37_path + 'iea37-ex%d.yaml' % n_wt)
if 0:
plt.plot(x, y, '2k')
for i, (x_, y_) in enumerate(zip(x, y)):
plt.annotate(i, (x_, y_))
plt.axis('equal')
plt.show()
site = IEA37Site(n_wt)
windTurbines = IEA37_WindTurbines(iea37_path + 'iea37-335mw.yaml')
wake_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
aep_ilk = wake_model(x, y, wd=np.arange(0, 360, 22.5), ws=[9.8]).aep_ilk(normalize_probabilities=True)
aep_MW_l = aep_ilk.sum((0, 2)) * 1000
# test that the result is equal to results provided for IEA task 37
npt.assert_almost_equal(aep_ref[0], aep_MW_l.sum(), 5)
npt.assert_array_almost_equal(aep_ref[1], aep_MW_l, 5)
def test_YZGrid_parallel():
site = IEA37Site(16)
x, y = site.initial_position.T
m = x < -1000
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x[m], y[m], wd=0)
fm = simulation_result.flow_map(grid=YZGrid(-1000, z=110, resolution=20))
if 0:
simulation_result.flow_map(grid=YZGrid(-1000)).plot_wake_map()
plt.plot(fm.X[0], fm.Y[0], '.')
print(np.round(fm.WS_eff_xylk[:, 0, 0, 0], 2).data.tolist())
plt.plot(fm.X[0], fm.WS_eff_xylk[:, 0, 0, 0] * 100, label='ws*100')
plt.legend()
plt.show()
npt.assert_array_almost_equal(fm.WS_eff_xylk[:, 0, 0, 0], [
7.32, 7.02, 6.63, 8.86, 8.79, 8.71, 8.63, 8.53, 8.42, 8.3, 8.16, 7.99,
7.81, 7.59, 7.33, 7.0, 6.52, 9.8, 9.8, 9.8
], 2)
def test_IEA37SimpleBastankhahGaussian_wake_map():
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines(iea37_path + 'iea37-335mw.yaml')
wake_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
x_j = np.linspace(-1500, 1500, 200)
y_j = np.linspace(-1500, 1500, 100)
flow_map = wake_model(x, y, wd=0, ws=9).flow_map(HorizontalGrid(x_j, y_j))
X, Y, Z = flow_map.X, flow_map.Y, flow_map.WS_eff_xylk[:, :, 0, 0]
# test that the result is equal to last run (no evidens that these number are correct)
ref = [
3.32, 4.86, 7.0, 8.1, 7.8, 7.23, 6.86, 6.9, 7.3, 7.82, 8.11, 8.04,
7.87, 7.79, 7.85, 8.04, 8.28
]
if 0:
flow_map.plot_wake_map()
plt.plot(X[49, 100:133:2], Y[49, 100:133:2], '.-')
plt.show()
npt.assert_array_almost_equal(Z[49, 100:133:2], ref, 2)
def test_deflection_model(deflectionModel, dy10d):
site = IEA37Site(16)
x, y = [0], [0]
windTurbines = V80()
D = windTurbines.diameter()
wfm = IEA37SimpleBastankhahGaussian(site, windTurbines, deflectionModel=deflectionModel())
yaw_ilk = np.reshape([-30], (1, 1, 1))
sim_res = wfm(x, y, yaw=yaw_ilk, wd=270, ws=10)
fm = sim_res.flow_map(XYGrid(x=np.arange(-D, 10 * D + 10, 10)))
min_WS_line = fm.min_WS_eff()
if 0:
plt.figure(figsize=(14, 3))
fm.plot_wake_map()
min_WS_line.plot()
plt.plot(10 * D, dy10d * D, '.', label="Ref, 10D")
plt.legend()
plt.show()
npt.assert_almost_equal(min_WS_line.interp(x=10 * D).item() / D, dy10d)
def test_plot_deflection_grid(deflectionModel):
site = IEA37Site(16)
x, y = [0], [0]
windTurbines = V80()
D = windTurbines.diameter()
wfm = IEA37SimpleBastankhahGaussian(site, windTurbines, deflectionModel=deflectionModel())
yaw_ilk = np.reshape([-30], (1, 1, 1))
sim_res = wfm(x, y, yaw=yaw_ilk, wd=270, ws=10)
fm = sim_res.flow_map(XYGrid(x=np.arange(-D, 10 * D + 10, 10)))
plt.figure(figsize=(14, 3))
fm.plot_wake_map()
fm.plot_deflection_grid()
min_WS_line = fm.min_WS_eff()
min_WS_line.plot()
plt.legend()
plt.title(wfm.deflectionModel)
if 0:
plt.show()
plt.close('all')
def main():
if __name__ == '__main__':
from py_wake import NOJ
from py_wake import IEA37SimpleBastankhahGaussian, Fuga
import py_wake
import os
import matplotlib.pyplot as plt
from py_wake.examples.data.hornsrev1 import HornsrevV80, Hornsrev1Site
from py_wake.examples.data.iea37._iea37 import IEA37Site
LUT_path = os.path.dirname(
py_wake.__file__
) + '/tests/test_files/fuga/2MW/Z0=0.03000000Zi=00401Zeta0=0.00E+00/'
wt_x, wt_y = IEA37Site(16).initial_position.T
windTurbines = HornsrevV80()
site = Hornsrev1Site()
wake_models = [
NOJ(site, windTurbines),
IEA37SimpleBastankhahGaussian(site, windTurbines),
Fuga(LUT_path, site, windTurbines)
]
for wake_model in wake_models:
# Calculate AEP
sim_res = wake_model(wt_x, wt_y)
# Plot wake map
plt.figure(wake_model.__class__.__name__)
plt.title('AEP: %.2f GWh' % sim_res.aep().sum())
flow_map = sim_res.flow_map(wd=[0], ws=[9])
flow_map.plot_wake_map()
flow_map.plot_windturbines()
plt.show()
def main():
if __name__ == '__main__':
from py_wake.examples.data.iea37 import IEA37Site, IEA37_WindTurbines
from py_wake import IEA37SimpleBastankhahGaussian
import matplotlib.pyplot as plt
site = IEA37Site(16)
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x, y)
fm = simulation_result.flow_map(wd=30)
fm.plot_wake_map()
plt.figure()
fm.plot(fm.power_xylk().sum(['wd', 'ws']) * 1e-3, "Power [kW]")
fm = simulation_result.flow_map(grid=HorizontalGrid(resolution=50))
plt.figure()
fm.plot(fm.aep_xy(), "AEP [GWh]")
plt.show()
def test_YZGrid_terrain():
site = ParqueFicticioSite(distance=StraightDistance())
x, y = site.initial_position.T
windTurbines = IEA37_WindTurbines()
wind_farm_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
simulation_result = wind_farm_model(x, y, wd=270, ws=10)
x = x.max() + 10
fm = simulation_result.flow_map(grid=YZGrid(x, z=110, resolution=20))
y = fm.X[0]
x = np.zeros_like(y) + x
z = site.elevation(x, y)
simulation_result.flow_map().plot_wake_map()
plt.plot(x, y, '.')
plt.figure()
simulation_result.flow_map(grid=YZGrid(x.max() + 10)).plot_wake_map()
plt.plot(y, z + 110, '.')
if 0:
print(np.round(fm.WS_eff_xylk[0, :, 0, 0], 2).tolist())
plt.show()
plt.close()
npt.assert_array_almost_equal(fm.WS_eff_xylk[:, 0, 0, 0],
[10.0, 10.0, 11.69, 5.8, 8.7, 5.15, 10.63, 4.99, 10.73, 8.16, 13.24, 5.59,
12.4, 7.02, 10.91, 7.21, 8.77, 8.41, 9.99, 10.0], 2)
def test_YZGrid_variables():
site = IEA37Site(16)
x, y = [0], [0]
windTurbines = IEA37_WindTurbines()
wf_model = IEA37SimpleBastankhahGaussian(site, windTurbines)
sim_res = wf_model(x, y)
fm = sim_res.flow_map(grid=YZGrid(x=100, y=None, resolution=100,
extend=.1),
wd=270,
ws=None)
fm.WS_eff.plot()
plt.plot(fm.y[::10], fm.y[::10] * 0 + 110, '.')
if 0:
print(np.round(fm.WS_eff.interp(h=110)[::10].squeeze().values, 4))
plt.show()
plt.close('all')
npt.assert_array_almost_equal(
fm.WS_eff.interp(h=110)[::10].squeeze(), [
9.1461, 8.4157, 7.3239, 6.058, 5.022, 4.6455, 5.1019, 6.182, 7.446,
8.506
], 4)
import multiprocessing
from py_wake.examples.data.hornsrev1 import Hornsrev1Site, wt_x, wt_y
from py_wake import IEA37SimpleBastankhahGaussian
from py_wake.tests.check_speed import timeit
import numpy as np
from py_wake.tests import npt
from py_wake.wind_turbines import WindTurbines
from py_wake.examples.data import wtg_path
wt = WindTurbines.from_WAsP_wtg(wtg_path + "Vestas-V80.wtg")
site = Hornsrev1Site()
wf_model = IEA37SimpleBastankhahGaussian(site, wt)
wd_lst = np.arange(0, 360, 10)
def aep_wd(args):
x, y, wd = args
return wf_model(x, y, wd=wd, ws=None).aep().sum()
def aep_all_multiprocessing(pool, x, y):
return np.sum(pool.map(aep_wd, [(x, y, i) for i in wd_lst]))
def aep_wfm_xy(args):
wfm, x, y = args
return wfm(x, y, wd=wd_lst).aep().sum()
def aep_xy(args):
def get_wfm(grad=True):
wt = WindTurbines.from_WAsP_wtg(wtg_path + "Vestas-V80.wtg", )
site = Hornsrev1Site()
return IEA37SimpleBastankhahGaussian(site, wt)
