Open and run in Colab (interactive) Edit on Gitlab

Experiment: Validation

In this notebook, you can compare and validate your own WindFarmModel (i.e. combination of engineering submodels) with data from RANS, LES and measurements.

Install PyWake if needed

[1]:
# Install PyWake if needed
try:
    import py_wake
except ModuleNotFoundError:
    !pip install git+https://gitlab.windenergy.dtu.dk/TOPFARM/PyWake.git

Now we install some basic Python elements

[2]:
import numpy as np
import matplotlib as plt
plt.rcParams.update({'figure.max_open_warning': 0})

Now we also import all available models

[3]:
from py_wake.deficit_models import *
from py_wake.wind_farm_models import *
from py_wake.rotor_avg_models import *
from py_wake.superposition_models import *
from py_wake.deflection_models import *
from py_wake.turbulence_models import *
from py_wake.ground_models import *
[4]:
# prepare for the model combination tool
from py_wake.utils.model_utils import get_models, get_signature
from ipywidgets import interact
from IPython.display import HTML, display, Javascript

# Fix ipywidget label width
display(HTML('''<style>.widget-label { min-width: 20ex !important; }</style>'''))

def print_signature(windFarmModel, **kwargs):
    s = """# insert windFarmModel code below
wfm = %s

validation.add_windFarmModel('MyModelName', wfm)"""% get_signature(windFarmModel, kwargs, 1)

    # Write windFarmModel code to cell starting "# insert windFarmModel code below"
    display(Javascript("""
for (var cell of IPython.notebook.get_cells()) {
    if (cell.get_text().startsWith("# insert windFarmModel code below")){
        cell.set_text(`%s`)
    }
}"""%s))

# setup list of models
models = {n:[(getattr(m,'__name__',m), m) for m in get_models(cls)]
          for n,cls in [('windFarmModel', WindFarmModel),
                        ('wake_deficitModel', WakeDeficitModel),
                        ('rotorAvgModel', RotorAvgModel),
                        ('superpositionModel', SuperpositionModel),
                        ('blockage_deficitModel', BlockageDeficitModel),
                        ('deflectionModel',DeflectionModel),
                        ('turbulenceModel', TurbulenceModel),
                        ('groundModel', GroundModel)
                        ]}
[5]:
from py_wake.validation.validation import Validation, ValidationSite, ValidationWindTurbines
site, windTurbines = ValidationSite(), ValidationWindTurbines()

Setup Validation

Instantiate new validation. This cell removes previously added WindFarmModel

[6]:
validation = Validation()

Add WindFarmModels

Add as many windFarmModels as you wish with the function

validation.add_windFarmModel(name, windFarmModel, line_style='-')

for example:

wfm = PropagateDownwind(site, windTurbines, wake_deficitModel=NOJDeficit(k=0.04))
validation.add_windFarmModel("NOJ(k=0.04)", wfm, ':')

You can use the dropdown boxes here to update the code cell below. Note, that some models needs manual specification of some non-optional arguments

[7]:
_ = interact(print_signature, **models)

In the cell below:

  • Replace MyModelName with a name for the windFarmModel

  • Set unspecified arguments, if any

  • Run the cell below to add the windFarmModel to the validation

[8]:
# insert windFarmModel code below
wfm = PropagateDownwind(
    site,
    windTurbines,
    wake_deficitModel=NOJDeficit(
        k=0.1,
        use_effective_ws=False,
        rotorAvgModel=AreaOverlapAvgModel(),
        groundModel=None),
    superpositionModel=LinearSum(),
    deflectionModel=None,
    turbulenceModel=None,
    rotorAvgModel=None)

validation.add_windFarmModel('MyModelName', wfm)
[9]:
# print name of added wind farm models
print(validation.windFarmModel_dict.keys())
dict_keys(['MyModelName'])

Single wake deficit validation

[10]:
validation.plot_deficit_profile()
../../_images/notebooks_exercises_Validation_18_0.png
../../_images/notebooks_exercises_Validation_18_1.png
../../_images/notebooks_exercises_Validation_18_2.png
../../_images/notebooks_exercises_Validation_18_3.png
../../_images/notebooks_exercises_Validation_18_4.png
../../_images/notebooks_exercises_Validation_18_5.png

Single wake integrated momentum deficit validation

[11]:
validation.plot_integrated_deficit()
../../_images/notebooks_exercises_Validation_20_0.png
../../_images/notebooks_exercises_Validation_20_1.png
../../_images/notebooks_exercises_Validation_20_2.png
../../_images/notebooks_exercises_Validation_20_3.png
../../_images/notebooks_exercises_Validation_20_4.png
../../_images/notebooks_exercises_Validation_20_5.png

Wind farm power validation

[12]:
validation.plot_multiwake_power()
../../_images/notebooks_exercises_Validation_22_0.png
../../_images/notebooks_exercises_Validation_22_1.png
../../_images/notebooks_exercises_Validation_22_2.png
../../_images/notebooks_exercises_Validation_22_3.png
../../_images/notebooks_exercises_Validation_22_4.png
../../_images/notebooks_exercises_Validation_22_5.png
../../_images/notebooks_exercises_Validation_22_6.png
../../_images/notebooks_exercises_Validation_22_7.png
../../_images/notebooks_exercises_Validation_22_8.png
../../_images/notebooks_exercises_Validation_22_9.png
../../_images/notebooks_exercises_Validation_22_10.png
../../_images/notebooks_exercises_Validation_22_11.png