Vertical flow map is strange at the rotor disk
I have made a simple model, and computed the flow map at the rotor disk. The result looks fine in the left part, but wrong in the right one. However, everything is fine if I move the plane just 1 mm before or after the rotor.
import numpy as np
import matplotlib.pyplot as plt
from py_wake.deficit_models import BastankhahGaussianDeficit
from py_wake.wind_farm_models import PropagateDownwind
from py_wake.flow_map import YZGrid
from py_wake.site import UniformSite
from py_wake.examples.data.dtu10mw import DTU10MW
windTurbines = DTU10MW()
z_hub = windTurbines.hub_height()
radius = windTurbines.diameter() / 2
site = UniformSite(initial_position=np.array([
[0.0, 0.0],
[8 * radius, 0.0],
]))
x, y = site.initial_position.T
wfm = PropagateDownwind(site, windTurbines, wake_deficitModel=BastankhahGaussianDeficit())
sim_res = wfm(x, y, wd=270, ws=10)
# Plot top view.
plt.figure()
sim_res.flow_map(wd=270).plot_wake_map()
# Plot front view.
x_turbine = x[1]
y_turbine = y[1]
z_grid = z_hub + np.linspace(-1.1 * radius, 1.1 * radius, 30)
y_grid = y_turbine + np.linspace(-1.1 * radius, 1.1 * radius, 30)
fig, ax = plt.subplots(ncols=3, figsize=(14, 8))
for i, dx in zip(range(3), [-0.001, 0.0, +0.001]):
x_plane = x_turbine + dx
ax[i].set_title(f'x = {x_plane:.3f}, y = {y_turbine:.1f}')
sim_res.flow_map(YZGrid(x_plane, y_grid, z_grid)).plot_wake_map(ax=ax[i])
plt.tight_layout()
Edited by Riccardo Riva