Use wind_speed to set URef

WISDEM · Sep 19, 2024 · e838be8 · e838be8
1 parent aaf6d2a
commit e838be8
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 7 deletions.
diff --git a/examples/02_run_openfast_cases/modeling_options.yaml b/examples/02_run_openfast_cases/modeling_options.yaml
@@ -49,12 +49,12 @@ DLC_driver:
         - DLC: "6.1"
           analysis_time: 10.
           transient_time: 0.1
+          wind_speed: [46.789]
           turbulent_wind:
               HubHt: 142
               RefHt: 142
               GridHeight: 275
               GridWidth: 275
-              URef: 46.789
               PLExp: 0.14
               IECturbc: 0.12
 

diff --git a/weis/aeroelasticse/openmdao_openfast.py b/weis/aeroelasticse/openmdao_openfast.py
@@ -1814,7 +1814,7 @@ def run_FAST(self, inputs, discrete_inputs, fst_vt):
                 # TODO AG: dlc driver (dlc_generator) does not seem to have IECturbc parameter, Check the logic below
                 if dlc_generator.cases[i_case].IECturbc > 0:    # use custom TI for DLC case
                     dlc_generator.cases[i_case].IECturbc = str(dlc_generator.cases[i_case].IECturbc)
-                    dlc_generator.cases[i_case].IEC_WindType = 'NTM'
+                    dlc_generator.cases[i_case].IEC_WindType = 'NTM'        # must use NTM for custom TI
                 else:
                     dlc_generator.cases[i_case].IECturbc = wt_class
                 # Reference height for wind speed

diff --git a/weis/dlc_driver/dlc_generator.py b/weis/dlc_driver/dlc_generator.py
@@ -906,7 +906,8 @@ def generate_6p1(self, dlc_options):
         if 'yaw_misalign' not in dlc_options:
             dlc_options['yaw_misalign'] = [-8,8]
 
-        dlc_options['wind_speed'] = [self.V_e50]
+        if not dlc_options['wind_speed']:
+            dlc_options['wind_speed'] = [self.V_e50]
 
         # parked options
         dlc_options['turbine_status'] = 'parked-idling'
@@ -946,7 +947,8 @@ def generate_6p2(self, dlc_options):
         if 'yaw_misalign' not in dlc_options:
             dlc_options['yaw_misalign'] = np.arange(-180+15,180+15,15).tolist()     # -180 is not valid in OF
 
-        dlc_options['wind_speed'] = [self.V_e50]   # placeholder, could be anything as long as the length is 1, since the EWM50 is just a single speed that turbsim will determine 
+        if not dlc_options['wind_speed']:
+            dlc_options['wind_speed'] = [self.V_e50]
 
         # parked options
         dlc_options['turbine_status'] = 'parked-idling'
@@ -986,8 +988,9 @@ def generate_6p3(self, dlc_options):
         else: # default
             dlc_options['yaw_misalign'] = [-20.,20.]
 
-        dlc_options['wind_speed'] = [self.V_e1]   # placeholder, could be anything as long as the length is 1, since the EWM50 is just a single speed that turbsim will determine 
-
+        if not dlc_options['wind_speed']:
+            dlc_options['wind_speed'] = [self.V_e1]
+
         # parked options
         dlc_options['turbine_status'] = 'parked-idling'
         dlc_options['wake_mod'] = 0

diff --git a/weis/inputs/modeling_schema.yaml b/weis/inputs/modeling_schema.yaml
@@ -3361,7 +3361,7 @@ properties:
                                     type: number
                                     unit: m/s
                                     default: -1
-                                    description: Mean (total) velocity at the reference height [m/s] (or 'default' for JET velocity profile) [must be 1-hr mean for API model; otherwise is the mean over AnalysisTime seconds]
+                                    description: URef is an input to TurbSim, but it's not honored in WEIS.  Please use the wind_speed input to each DLC instead.
                                 IECturbc:
                                     type: number
                                     unit: '(-)'