Partially connect yaw to foxes

rafmudaf · Mar 21, 2024 · 042df6d · 042df6d
1 parent 89a76e9
commit 042df6d
Show file tree

Hide file tree

Showing 8 changed files with 78 additions and 35 deletions.
diff --git a/...astankhah_deflection/energy_resource.yaml → ...nkhah2016_deflection/energy_resource.yaml b/...astankhah_deflection/energy_resource.yaml → ...nkhah2016_deflection/energy_resource.yaml
diff --git a/...24/bastankhah_deflection/energy_site.yaml → ...astankhah2016_deflection/energy_site.yaml b/...24/bastankhah_deflection/energy_site.yaml → ...astankhah2016_deflection/energy_site.yaml
diff --git a/...ankhah_deflection/wind_energy_system.yaml → ...ah2016_deflection/wind_energy_system.yaml b/...ankhah_deflection/wind_energy_system.yaml → ...ah2016_deflection/wind_energy_system.yaml
@@ -4,20 +4,17 @@ wind_farm: !include wind_farm.yaml
 attributes:
   analyses:
     yaw_angles:
-    - 12.0
+    - 10.0
     wake_model:
       velocity:
         name: bastankhah2016
         parameters:
           alpha: 0.58
           beta: 0.077
-          k: 0.38
+          k: 0.022
       deflection:
         name: bastankhah2016_deflection
         parameters:
-          ad: 0.0
-          bd: 0.0
           alpha: 0.58
           beta: 0.077
-          ka: 0.38
-          kb: 0.004
+          k: 0.022
diff --git a/...2024/bastankhah_deflection/wind_farm.yaml → .../bastankhah2016_deflection/wind_farm.yaml b/...2024/bastankhah_deflection/wind_farm.yaml → .../bastankhah2016_deflection/wind_farm.yaml
diff --git a/examples/torque2024_1turbine.ipynb b/examples/torque2024_1turbine.ipynb
diff --git a/examples/wesc2023_demo.ipynb b/examples/wesc2023_demo.ipynb
diff --git a/wcomp/floris_interface.py b/wcomp/floris_interface.py
@@ -59,12 +59,10 @@
     "bastankhah2016_deflection": {
         "model_ref": "gauss",
         "parameters": {
-            "ad": "ad",
-            "bd": "bd",
             "alpha": "alpha",
             "beta": "beta",
-            "ka": "ka",
-            "kb": "kb",
+            "ka": "k",
+            "kb": "k",
         }
     },
 }

diff --git a/wcomp/foxes_interface.py b/wcomp/foxes_interface.py
@@ -7,6 +7,7 @@
 import matplotlib.pyplot as plt
 import numpy as np
 import pandas as pd
+import foxes
 from foxes.algorithms.downwind import Downwind
 from foxes.core import WindFarm, WakeModel
 from foxes.input.farm_layout import add_from_df
@@ -78,6 +79,7 @@
         "parameters": {
             "alpha": "alpha",
             "beta": "beta",
+            "k": "k",
         }
     },
 }
@@ -314,6 +316,8 @@ def read_analyses(self, analyses, mbook, farm, states, keymap={}, **algo_pars):
             The algorithm
 
         """
+        yaw_angles = np.array([analyses["yaw_angles"]])
+
         wes_analysis = analyses
         _velocity_model_mapping = WAKE_MODEL_MAPPING[wes_analysis["wake_model"]["velocity"]["name"]]
         _velocity_model = _velocity_model_mapping["model_ref"]
@@ -322,6 +326,26 @@ def read_analyses(self, analyses, mbook, farm, states, keymap={}, **algo_pars):
             for k, v in _velocity_model_mapping["parameters"].items()
         }
 
+        if wes_analysis["wake_model"]["deflection"]["name"] is not None:
+            # NOTE: foxes supports only one deflection model, so there's no need to parse the
+            # deflection model settings from windIO.
+            # Check the name, and error if it isn't Bastankhah2016
+            if wes_analysis["wake_model"]["deflection"]["name"] != "bastankhah2016_deflection":
+                raise ValueError("foxes supports only Bastankhah2016 for the deflection model.")
+
+            # _deflection_model_mapping = WAKE_MODEL_MAPPING[wes_analysis["wake_model"]["deflection"]["name"]]
+            # _deflection_model = _deflection_model_mapping["model_ref"]
+            # _deflection_model_parameters = {
+            #     k: wes_analysis["wake_model"]["deflection"]["parameters"][v]
+            #     for k, v in _deflection_model_mapping["parameters"].items()
+            # }
+            # _deflection_model = _deflection_model(**_deflection_model_parameters)
+            mbook.turbine_models["set_yawm"] = foxes.models.turbine_models.SetFarmVars()
+            mbook.turbine_models["set_yawm"].add_var(FV.YAWM, yaw_angles)
+            wake_frame="yawed"
+        else:
+            wake_frame="rotor_wd"
+
         temp_model_name = "this_model"
         mbook.wake_models[temp_model_name] = _velocity_model(
             **_velocity_model_parameters,
@@ -335,6 +359,7 @@ def read_analyses(self, analyses, mbook, farm, states, keymap={}, **algo_pars):
             verbosity=0,
             rotor_model="grid16",
             wake_models=[temp_model_name],
+            wake_frame=wake_frame,
             **algo_pars
         )