updated for kestrel with calibrations

NREL · Jun 4, 2024 · 302b5b9 · 302b5b9
1 parent 1d76d1f
commit 302b5b9
Show file tree

Hide file tree

Showing 3 changed files with 8 additions and 3 deletions.
diff --git a/demo/documented/studies/powercurve_validation/single_turbine.yaml b/demo/documented/studies/powercurve_validation/single_turbine.yaml
@@ -25,7 +25,9 @@ turbines:
     # CTprime0: 1.350       # Region 2 Thrust Coeff     | -
     # CTprime0: 1.431       # Region 2 Thrust Coeff     | -
     # CTprime0: 1.639       # Region 2 Thrust Coeff     | -
-    Prated: 3077551.02040 # Specific Rated Power      | W
+    Prated: 2570282.74287 # Specific Rated Power      | W
+    # Prated: 2532257.89912 # Specific Rated Power      | W
+    # Prated: 3077551.02040 # Specific Rated Power      | W
     Trated: 368163.26530  # Specific Thrust           | N
     force: sine           # radial force distribution | -
     HH: 110               # Hub Height                | m

diff --git a/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml b/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml
@@ -25,7 +25,9 @@ turbines:
     # CTprime0: 1.350       # Region 2 Thrust Coeff     | -
     # CTprime0: 1.431       # Region 2 Thrust Coeff     | -
     # CTprime0: 1.639       # Region 2 Thrust Coeff     | -
-    Prated: 3077551.02040 # Specific Rated Power      | W
+    Prated: 2570282.74287 # Specific Rated Power      | W
+    # Prated: 2532257.89912 # Specific Rated Power      | W
+    # Prated: 3077551.02040 # Specific Rated Power      | W
     Trated: 368163.26530  # Specific Thrust           | N
     force: sine           # radial force distribution | -
     HH: 110.0             # Hub Height                | m

diff --git a/windse/turbine_types/ActuatorDiskSimplePowerCurve.py b/windse/turbine_types/ActuatorDiskSimplePowerCurve.py
@@ -195,7 +195,8 @@ def force_gradient(self):
 
     def power(self, u, inflow_angle):
         # adjust for turbine inefficiency
-        return self.calibration_factor_CPprime0*self.mCPprime0/(self.calibration_factor_CTprime0*self.mCTprime0)*dot(-self.tf,u)/1.0e6  # report in megawatts
+        magic_number = 1.1291826687642552  # comes from calibration studies
+        return magic_number*self.calibration_factor_CPprime0*self.mCPprime0/(self.calibration_factor_CTprime0*self.mCTprime0)*dot(-self.tf,u)/1.0e6  # report in megawatts
 
     def thrust(self, u, inflow_angle):
         # adjust for turbine inefficiency