DTU10MW_RWT turbine

py_wake/examples/data/DTU10MW_RWT.py

+import numpy as np
+from py_wake.wind_turbines import OneTypeWindTurbines
+power_curve = np.array([[4, 280.2],
+                        [5, 799.1],
+                        [6, 1532.7],
+                        [7, 2506.1],
+                        [8, 3730.7],
+                        [9, 5311.8],
+                        [10, 7286.5],
+                        [11, 9698.3],
+                        [12, 10639.1],
+                        [13, 10648.5],
+                        [14, 10639.3],
+                        [15, 10683.7],
+                        [16, 10642],
+                        [17, 10640],
+                        [18, 10639.9],
+                        [19, 10652.8],
+                        [20, 10646.2],
+                        [21, 10644],
+                        [22, 10641.2],
+                        [23, 10639.5],
+                        [24, 10643.6],
+                        [25, 10635.7],
+                        ])
+ct_curve = np.array([[4, 0.923],
+                     [5, 0.919],
+                     [6, 0.904],
+                     [7, 0.858],
+                     [8, 0.814],
+                     [9, 0.814],
+                     [10, 0.814],
+                     [11, 0.814],
+                     [12, 0.577],
+                     [13, 0.419],
+                     [14, 0.323],
+                     [15, 0.259],
+                     [16, 0.211],
+                     [17, 0.175],
+                     [18, 0.148],
+                     [19, 0.126],
+                     [20, 0.109],
+                     [21, 0.095],
+                     [22, 0.084],
+                     [23, 0.074],
+                     [24, 0.066],
+                     [25, 0.059],
+                     ])
+
+
+class DTU10MW(OneTypeWindTurbines):
+    '''
+    Data from:
+    Christian Bak, Frederik Zahle, Robert Bitsche, Taeseong Kim, Anders Yde, Lars Christian Henriksen, Anand Natarajan,
+    Morten Hartvig Hansen.“Description of the DTU 10 MW Reference Wind Turbine” DTU Wind Energy Report-I-0092, July 2013. Table 3.5
+
+    '''
+
+    def __init__(self):
+        OneTypeWindTurbines.__init__(
+            self,
+            'DTU10MW',
+            diameter=178.3,
+            hub_height=119,
+            ct_func=self._ct,
+            power_func=self._power,
+            power_unit='kW')
+
+    def _ct(self, u):
+        return np.interp(u, ct_curve[:, 0], ct_curve[:, 1])
+
+    def _power(self, u):
+        return np.interp(u, power_curve[:, 0], power_curve[:, 1])
+
+
+DTU10WM_RWT = DTU10MW
+
+
+def main():
+    wt = DTU10MW()
+    print('Diameter', wt.diameter())
+    print('Hub height', wt.hub_height())
+    ws = np.arange(3, 25)
+    import matplotlib.pyplot as plt
+    plt.plot(ws, wt.power(ws), '.-')
+    plt.show()
+
+
+if __name__ == '__main__':
+    main()