Fast.Farm generates errors at low wind speeds

Dear @Jason.Jonkman,
Hello! I hope you’re doing well.

I encountered an issue while running a FAST.Farm simulation to evaluate turbine fatigue due to wake effects at a relatively low wind speed with high turbulence intensity (V: 2 m/s, T.I: 37.1). Initially, I received the following error:

T2:FARM_UpdateStates:WD_UpdateStates:ThomasAlgorithm:Tridiagonal matrix is not diagonally
dominant, i.e., abs(b(i)) <= ( abs(a(i))+abs(c(i)) ). Try reducing the FAST.Farm timestep.

In response, I reduced the value of DT_Low in the .fstf file, which resulted in another error:

T1:FARM_UpdateStates:NearWakeCorrection:Wake model is not valid in the propeller-brake region, i.e., Ct > 2.0.

Upon reviewing the FAST.Farm results, it appears that this error occurs when the turbulence intensity is high, causing the wind speed to drop close to 0.3, resulting in a Ct value exceeding 2. According to the FAST.Farm manual, Ct values are notably high at low wind speeds, but I am struggling to find a suitable solution to this issue.

I would greatly appreciate any advice you can provide.

Thank you.

Dear @Young-Ju.Noh,

I would not expect a wind turbine to operate at such a low wind speed with such a high thrust coefficient. I think it is best to avoid this simulation set-up by not simulating with such low wind speeds or changing the wind turbine controller so that the turbine shuts down when the wind speed gets so low.

Best regards,

Dear @Jason.Jonkman,

Thank you for your prompt response and valuable advice.
To address this issue, I will explore either avoiding simulations with such low wind speeds or modifying the wind turbine controller settings to ensure the turbine shuts down when wind speeds drop to such low levels. Your guidance has been very helpful. Thank you once again for your assistance.

Best regards,

Hi @Jason.Jonkman ,
I ran into the same issue for a few cases like [vhub = 7 m/sec, TI = 14%] and [vhub = 8m/sec and TI = 17%]. Lowering DT_Low to 2 m/sec didn’t help. Is this expected?

Simulation were run on FAST.Farm 3.4.1 with NREL 2.8 kW turbine.

Abhineet

Dear @Abhineet.Gupta,

I’m not sure what you mean when you refer to an “NREL 2.8 kW turbine”. Are you using one of the OpenFAST models from the r-test?

How low is the ambient wind speed dipping in these simulations and are the corresponding rotor speed (tip speed ratio) and blade-pitch angles expected for this situation?

Best regards,

Hi @Jason.Jonkman
Thanks for the reply.

I am using this turbine “openfast-turbine-models/IEA-scaled/legacy_models/NREL-2.8-127 at main · NREL/openfast-turbine-models · GitHub”.

The wind speed dips to around 2.5 m/sec but recovers, However, the TSR drops to around 5.6 just before the error. The rotor speed drops to around 5 RPM.
The following plot shows these values for all 12 turbines in the farm


Here’s a a zoomed in plot focusing on data right before the crash.

Regards

Dear @Abhineet.Gupta,

I’m not really familiar with this OpenFAST model (which is a 2.8-MW wind turbine instead of a 2.8-kW turbine), but I would guess the issue is associated with the wind turbines with high tip-speed ratio (TSR, such as the purple or orange), which would correspond to a high thrust coefficient. It doesn’t look like the TSR dwells at a very high TSR for long, so perhaps adjusting (reducing) the corner frequency (f_c) of the low-pass time filter would help FAST.Farm make it through the periods of high tip-speed ratio. Otherwise, it may be necessary to modify the controller so that the TSR does not get so high.

Best regards,

Dear @Jason.Jonkman ,

When I ran a simulation using FAST.Farm to verify the wake steering effect at low wind speeds, I encountered the same error as others, and the calculation was forcibly terminated midway.

T2:FARM_UpdateStates:NearWakeCorrection:Wake model is not valid in the propeller-brake region,
i.e., Ct > 2.0.

The simulation conditions are as follows:
Average wind speed: 4 m/s
Turbulence intensity: IEC Category C

Turbine: 2x IEA-15-240-RWT on UMaine VolturnUS-S semi-submersible

The .fstf file used in the simulation is as follows:

--- SIMULATION CONTROL ---
False              Echo               - Echo input data to <RootName>.ech? (flag)
FATAL              AbortLevel         - Error level when simulation should abort (string) {"WARNING", "SEVERE", "FATAL"}
4500                TMax               - Total run time (s) [>=0.0]
False              UseSC              - Use a super controller? (flag)
3                  Mod_AmbWind        - Ambient wind model (-) (switch) {1: high-fidelity precursor in VTK format, 2: one InflowWind module, 3: multiple instances of InflowWind module}
2                  Mod_WaveField      - Wave field handling (-) (switch) {1: use individual HydroDyn inputs without adjustment, 2: adjust wave phases based on turbine offsets from farm origin}
0                  Mod_SharedMooring  - Shared mooring system model (switch) {0: None, 3=MoorDyn}}
--- SUPER CONTROLLER --- [used only for UseSC=True]
"unused"           SC_FileName        - Name/location of the dynamic library {.dll [Windows] or .so [Linux]} containing the Super Controller algorithms (quoated string)
--- SHARED MOORING SYSTEM --- [used only for Mod_SharedMoor>0]
""                 SharedMoorFile     - Name of file containing shared mooring system input parameters (quoted string) [used only when Mod_SharedMooring > 0]
0.25	           DT_Mooring         - Time step for farm-level mooring coupling with each turbine (s) [used only when Mod_SharedMooring > 0]
--- AMBIENT WIND: PRECURSOR IN VTK FORMAT --- [used only for Mod_AmbWind=1]
30.0                DT_Low-VTK         - Time step for low -resolution wind data input files; will be used as the global FAST.Farm time step (s) [>0.0]
1.0                DT_High-VTK        - Time step for high-resolution wind data input files (s) [>0.0]
"unused"           WindFilePath       - Path name to VTK wind data files from precursor (string)
False              ChkWndFiles        - Check all the ambient wind files for data consistency? (flag)
--- AMBIENT WIND: INFLOWWIND MODULE --- [used only for Mod_AmbWind=2 or 3]
4.00                    DT_Low      - Time step for low -resolution wind data interpolation; will be used as the global FAST.Farm time step (s) [>0.0]
0.25                    DT_High     - Time step for high-resolution wind data interpolation (s) [>0.0]
69                    NX_Low      - Number of low -resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
61                     NY_Low      - Number of low -resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
25                     NZ_Low      - Number of low -resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
-240.0                 X0_Low      - Origin of low -resolution spatial nodes in X direction for wind data interpolation (m)
-720.0                 Y0_Low      - Origin of low -resolution spatial nodes in Y direction for wind data interpolation (m)
3.0                    Z0_Low      - Origin of low -resolution spatial nodes in Z direction for wind data interpolation (m)
36.0                 dX_Low      - Spacing of low -resolution spatial nodes in X direction for wind data interpolation (m) [>0.0]
24.0                   dY_Low      - Spacing of low -resolution spatial nodes in Y direction for wind data interpolation (m) [>0.0]
24.0                   dZ_Low      - Spacing of low -resolution spatial nodes in Z direction for wind data interpolation (m) [>0.0]
101                     NX_High     - Number of high-resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
49                     NY_High     - Number of high-resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
51                     NZ_High     - Number of high-resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
"InflowWind_01.dat"           InflowFile         - Name of file containing InflowWind module input parameters (quoted string)
--- WIND TURBINES ---
2                  NumTurbines        - Number of wind turbines (-) [>=1]                         [last 6 columns below used only for Mod_AmbWind=2 or 3]
WT_X    WT_Y   WT_Z   WT_FASTInFile     X0_High  Y0_High  Z0_High  dX_High  dY_High  dZ_High
(m)     (m)    (m)    (string)          (m)      (m)      (m)      (m)      (m)      (m)
  0.0  0.0    0.0    "FFTest_WT1_01.fst"   -240  -144   1.0      6.0    6.0     6.0
  1440.0  0.0    0.0    "FFTest_WT2_01.fst"   1200  -144   1.0      6.0    6.0     6.0
--- WAKE DYNAMICS ---
1                      Mod_Wake    - Switch between wake formulations {1:Polar, 2:Curl, 3:Cartesian} (-) (switch)
6.155                  dr          - Radial increment of radial finite-difference grid (m) [>0.0]
60                     NumRadii    - Number of radii in the radial finite-difference grid (-) [>=2]
134                    NumPlanes   - Number of wake planes (-) [>=2]
0.17                   f_c         - Cutoff (corner) frequency of the low-pass time-filter for the wake advection, deflection, and meandering model [recommended=1.28*U0/R] (Hz) [>0.0] or DEFAULT [DEFAULT=12.5/R, R estimated from dr and NumRadii, not recommended]
DEFAULT                C_HWkDfl_O  - Calibrated parameter in the correction for wake deflection defining the horizontal offset at the rotor (m ) or DEFAULT [DEFAULT= 0.0 ]
DEFAULT                    C_HWkDfl_OY - Calibrated parameter in the correction for wake deflection defining the horizontal offset at the rotor scaled with yaw error (m/deg) or DEFAULT [DEFAULT= 0.0 if Mod_Wake is 2, 0.3 otherwise]
DEFAULT                C_HWkDfl_x  - Calibrated parameter in the correction for wake deflection defining the horizontal offset scaled with downstream distance (- ) or DEFAULT [DEFAULT= 0.0 ]
-0.004                  C_HWkDfl_xY - Calibrated parameter in the correction for wake deflection defining the horizontal offset scaled with downstream distance and yaw error (1/deg) or DEFAULT [DEFAULT= 0.0 if Mod_Wake is 2, -0.004 otherwise]
DEFAULT                C_NearWake  - Calibrated parameter for the near-wake correction (-) [>1.0 and <2.5] or DEFAULT [DEFAULT=1.8]
DEFAULT                k_vAmb      - Calibrated parameter for the influence of ambient turbulence in the eddy viscosity (-) [>=0.0] or DEFAULT [DEFAULT=0.05 ]
DEFAULT               k_vShr      - Calibrated parameter for the influence of the shear layer in the eddy viscosity (-) [>=0.0] or DEFAULT [DEFAULT=0.016]
DEFAULT                C_vAmb_DMin - Calibrated parameter in the eddy viscosity filter function for ambient turbulence defining the transitional diameter fraction between the minimum and exponential regions (-) [>=0.0 ] or DEFAULT [DEFAULT= 0.0 ]
DEFAULT                C_vAmb_DMax - Calibrated parameter in the eddy viscosity filter function for ambient turbulence defining the transitional diameter fraction between the exponential and maximum regions (-) [> C_vAmb_DMin ] or DEFAULT [DEFAULT= 1.0 ]
DEFAULT                C_vAmb_FMin - Calibrated parameter in the eddy viscosity filter function for ambient turbulence defining the value in the minimum region (-) [>=0.0 and <=1.0] or DEFAULT [DEFAULT= 1.0 ]
DEFAULT                C_vAmb_Exp  - Calibrated parameter in the eddy viscosity filter function for ambient turbulence defining the exponent in the exponential region (-) [> 0.0 ] or DEFAULT [DEFAULT= 0.01]
DEFAULT                C_vShr_DMin - Calibrated parameter in the eddy viscosity filter function for the shear layer defining the transitional diameter fraction between the minimum and exponential regions (-) [>=0.0 ] or DEFAULT [DEFAULT= 3.0 ]
DEFAULT                C_vShr_DMax - Calibrated parameter in the eddy viscosity filter function for the shear layer defining the transitional diameter fraction between the exponential and maximum regions (-) [> C_vShr_DMin ] or DEFAULT [DEFAULT=25.0 ]
DEFAULT                C_vShr_FMin - Calibrated parameter in the eddy viscosity filter function for the shear layer defining the value in the minimum region (-) [>=0.0 and <=1.0] or DEFAULT [DEFAULT= 0.2 ]
DEFAULT                C_vShr_Exp  - Calibrated parameter in the eddy viscosity filter function for the shear layer defining the exponent in the exponential region (-) [> 0.0 ] or DEFAULT [DEFAULT= 0.1 ]
DEFAULT                Mod_WakeDiam - Wake diameter calculation model (-) (switch) {1: rotor diameter, 2: velocity based, 3: mass-flux based, 4: momentum-flux based} or DEFAULT [DEFAULT=1]
DEFAULT                C_WakeDiam  - Calibrated parameter for wake diameter calculation (-) [>0.0 and <0.99] or DEFAULT [DEFAULT=0.95] [unused for Mod_WakeDiam=1]
DEFAULT                Mod_Meander - Spatial filter model for wake meandering (-) (switch) {1: uniform, 2: truncated jinc, 3: windowed jinc} or DEFAULT [DEFAULT=3]
DEFAULT                C_Meander   - Calibrated parameter for wake meandering (-) [>=1.0] or DEFAULT [DEFAULT=1.9]

I would like to ask you a few questions in order to resolve this situation.

・Is the calculation stopping midway because the TSR of the turbine in question is too high? Or is it because, as indicated by the “NearWakeCorrection” error message, it becomes impossible to calculate the wake near the turbine?

The TSR of T2 (downstream wind turbine) before the simulation was stopped was as shown in the image below.

・To prevent the TSR from increasing due to continued operation at low wind speeds, which parameters should I change to set the system to stop operation when the wind speed falls below 3 m/s?

The following is the DISCON.IN file currently in use.

!------- SIMULATION CONTROL ------------------------------------------------------------
1                   ! LoggingLevel		- {0: write no debug files, 1: write standard output .dbg-file, 2: LoggingLevel 1 + ROSCO LocalVars (.dbg2) 3: LoggingLevel 2 + complete avrSWAP-array (.dbg3)}
0                   ! DT_Out    		  - {Time step to output .dbg* files, or 0 to match sampling period of OpenFAST}
1                   ! Ext_Interface		- (0 - use standard bladed interface, 1 - Use the extened DLL interface introduced in OpenFAST 3.5.0.)
0                   ! Echo		        - (0 - no Echo, 1 - Echo input data to <RootName>.echo)

!------- CONTROLLER FLAGS -------------------------------------------------
2                   ! F_LPFType			  - (1: first-order low-pass filter, 2: second-order low-pass filter), [rad/s] (currently filters generator speed and pitch control signals
0                   ! IPC_ControlMode	- Turn Individual Pitch Control (IPC) for fatigue load reductions (pitch contribution) {0: off, 1: 1P reductions, 2: 1P+2P reductions}
2                   ! VS_ControlMode	- Generator torque control mode in above rated conditions (0- no torque control, 1- k*omega^2 with PI transitions, 2- WSE TSR Tracking, 3- Power-based TSR Tracking)}
0                   ! VS_ConstPower  	- Do constant power torque control, where above rated torque varies, 0 for constant torque}
1                   ! PC_ControlMode  - Blade pitch control mode {0: No pitch, fix to fine pitch, 1: active PI blade pitch control}
0                   ! Y_ControlMode   - Yaw control mode {0: no yaw control, 1: yaw rate control, 2: yaw-by-IPC}
1                   ! SS_Mode         - Setpoint Smoother mode {0: no setpoint smoothing, 1: introduce setpoint smoothing}
0                   ! PRC_Mode        - Power reference tracking mode{0: use standard rotor speed set points, 1: use PRC rotor speed setpoints}
2                   ! WE_Mode         - Wind speed estimator mode {0: One-second low pass filtered hub height wind speed, 1: Immersion and Invariance Estimator, 2: Extended Kalman Filter}
0                   ! PS_Mode         - Pitch saturation mode {0: no pitch saturation, 1: implement pitch saturation}
1                   ! SD_Mode         - Shutdown mode {0: no shutdown procedure, 1: pitch to max pitch at shutdown}
2                   ! Fl_Mode         - Floating specific feedback mode {0: no nacelle velocity feedback, 1: feed back translational velocity, 2: feed back rotational veloicty}
0                   ! TD_Mode         - Tower damper mode (0- no tower damper, 1- feed back translational nacelle accelleration to pitch angle
0                   ! TRA_Mode        - Tower resonance avoidance mode (0- no tower resonsnace avoidance, 1- use torque control setpoints to avoid a specific frequency
0                   ! Flp_Mode        - Flap control mode {0: no flap control, 1: steady state flap angle, 2: Proportional flap control, 2: Cyclic (1P) flap control}
0                   ! OL_Mode         - Open loop control mode {0: no open loop control, 1: open loop control vs. time, 2: rotor position control}
2                   ! PA_Mode         - Pitch actuator mode {0 - not used, 1 - first order filter, 2 - second order filter}
0                   ! PF_Mode         - Pitch fault mode {0 - not used, 1 - constant offset on one or more blades}
0                   ! AWC_Mode        - Active wake control {0 - not used, 1 - complex number method, 2 - Coleman transform method}
0                   ! Ext_Mode        - External control mode {0 - not used, 1 - call external dynamic library}
0                   ! ZMQ_Mode        - Fuse ZeroMQ interface {0: unused, 1: Yaw Control}
0                   ! CC_Mode         - Cable control mode [0- unused, 1- User defined, 2- Open loop control]
0                   ! StC_Mode        - Structural control mode [0- unused, 1- User defined, 2- Open loop control]

!------- FILTERS ----------------------------------------------------------
1.00810             ! F_LPFCornerFreq	  - Corner frequency (-3dB point) in the low-pass filters, [rad/s]
0.70000             ! F_LPFDamping		  - Damping coefficient {used only when F_FilterType = 2} [-]
1                   ! F_NumNotchFilts   - Number of notch filters placed on sensors
3.3550              ! F_NotchFreqs      - Natural frequency of the notch filters. Array with length F_NumNotchFilts
0.0000              ! F_NotchBetaNum    - Damping value of numerator (determines the width of notch). Array with length F_NumNotchFilts, [-]
0.2500              ! F_NotchBetaDen    - Damping value of denominator (determines the depth of notch). Array with length F_NumNotchFilts, [-]
0                   ! F_GenSpdNotch_N   - Number of notch filters on generator speed
0                   ! F_GenSpdNotch_Ind - Indices of notch filters on generator speed
1                   ! F_TwrTopNotch_N   - Number of notch filters on tower top acceleration signal
1                   ! F_TwrTopNotch_Ind - Indices of notch filters on tower top acceleration signal
0.62830             ! F_SSCornerFreq    - Corner frequency (-3dB point) in the first order low pass filter for the setpoint smoother, [rad/s].
0.20944             ! F_WECornerFreq    - Corner frequency (-3dB point) in the first order low pass filter for the wind speed estimate [rad/s].
0.17952             ! F_YawErr          - Low pass filter corner frequency for yaw controller [rad/s].
0.213000  1.000000  ! F_FlCornerFreq    - Natural frequency and damping in the second order low pass filter of the tower-top fore-aft motion for floating feedback control [rad/s, -].
0.01042             ! F_FlHighPassFreq  - Natural frequency of first-order high-pass filter for nacelle fore-aft motion [rad/s].
10.4616 1.0000      ! F_FlpCornerFreq   - Corner frequency and damping in the second order low pass filter of the blade root bending moment for flap control

!------- BLADE PITCH CONTROL ----------------------------------------------
30                  ! PC_GS_n			- Amount of gain-scheduling table entries
0.062094  0.088628  0.109270  0.127148  0.143223  0.158025  0.171764  0.184722  0.197105  0.209011  0.220338  0.231364  0.242063  0.252399  0.262569  0.272384  0.282068  0.291504  0.300785  0.309890  0.318839  0.327654  0.336324  0.344874  0.353301  0.361612  0.369834  0.377928  0.385970  0.393863                ! PC_GS_angles	    - Gain-schedule table: pitch angles [rad].
-1.143968  -0.988219  -0.859873  -0.752282  -0.660788  -0.582031  -0.513524  -0.453388  -0.400177  -0.352760  -0.310240  -0.271895  -0.237139  -0.205490  -0.176551  -0.149987  -0.125517  -0.102903  -0.081942  -0.062459  -0.044302  -0.027341  -0.011462  0.003436  0.017441  0.030631  0.043075  0.054835  0.065966  0.076516                ! PC_GS_KP		- Gain-schedule table: pitch controller kp gains [s].
-0.119726  -0.108178  -0.098661  -0.090684  -0.083900  -0.078060  -0.072981  -0.068522  -0.064576  -0.061061  -0.057908  -0.055065  -0.052488  -0.050141  -0.047995  -0.046026  -0.044211  -0.042535  -0.040980  -0.039536  -0.038189  -0.036932  -0.035754  -0.034650  -0.033611  -0.032633  -0.031711  -0.030839  -0.030013  -0.029231                ! PC_GS_KI		- Gain-schedule table: pitch controller ki gains [-].
0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000                ! PC_GS_KD			- Gain-schedule table: pitch controller kd gains
0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000                ! PC_GS_TF			- Gain-schedule table: pitch controller tf gains (derivative filter)
1.570000000000      ! PC_MaxPit			- Maximum physical pitch limit, [rad].
0.000000000000      ! PC_MinPit			- Minimum physical pitch limit, [rad].
0.034900000000      ! PC_MaxRat			- Maximum pitch rate (in absolute value) in pitch controller, [rad/s].
-0.03490000000      ! PC_MinRat			- Minimum pitch rate (in absolute value) in pitch controller, [rad/s].
0.791680000000      ! PC_RefSpd			- Desired (reference) HSS speed for pitch controller, [rad/s].
0.000000000000      ! PC_FinePit		- Record 5: Below-rated pitch angle set-point, [rad]
0.017450000000      ! PC_Switch			- Angle above lowest minimum pitch angle for switch, [rad]

!------- INDIVIDUAL PITCH CONTROL -----------------------------------------
8.592000  10.740000  ! IPC_Vramp		- Start and end wind speeds for cut-in ramp function. First entry: IPC inactive, second entry: IPC fully active. [m/s]
2                   ! IPC_SatMode		- IPC Saturation method (0 - no saturation (except by PC_MinPit), 1 - saturate by PS_BldPitchMin, 2 - saturate sotfly (full IPC cycle) by PC_MinPit, 3 - saturate softly by PS_BldPitchMin)
0.3                 ! IPC_IntSat		- Integrator saturation (maximum signal amplitude contribution to pitch from IPC), [rad]
0.000e+00 0.000e+00 ! IPC_KP			- Proportional gain for the individual pitch controller: first parameter for 1P reductions, second for 2P reductions, [-]
0.000e+00 0.000e+00 ! IPC_KI			- Integral gain for the individual pitch controller: first parameter for 1P reductions, second for 2P reductions, [-]
0.000000  0.000000  ! IPC_aziOffset		- Phase offset added to the azimuth angle for the individual pitch controller, [rad]. 
0.0                 ! IPC_CornerFreqAct - Corner frequency of the first-order actuators model, to induce a phase lag in the IPC signal {0: Disable}, [rad/s]

!------- VS TORQUE CONTROL ------------------------------------------------
96.55000000000      ! VS_GenEff			- Generator efficiency mechanical power -> electrical power, [should match the efficiency defined in the generator properties!], [%]
19624046.66639      ! VS_ArSatTq		- Above rated generator torque PI control saturation, [Nm]
4500000.000000      ! VS_MaxRat			- Maximum torque rate (in absolute value) in torque controller, [Nm/s].
21586451.33303      ! VS_MaxTq			- Maximum generator torque in Region 3 (HSS side), [Nm].
0.000000000000      ! VS_MinTq			- Minimum generator torque (HSS side), [Nm].
0.523600000000      ! VS_MinOMSpd		- Minimum generator speed [rad/s]
31320452.18335      ! VS_Rgn2K		- Generator torque constant in Region 2 (HSS side). Only used in VS_ControlMode = 1,3
15000000.00000      ! VS_RtPwr			- Wind turbine rated power [W]
19624046.66639      ! VS_RtTq			- Rated torque, [Nm].
0.791680000000      ! VS_RefSpd			- Rated generator speed [rad/s]
1                   ! VS_n				- Number of generator PI torque controller gains
-38016920.26262      ! VS_KP				- Proportional gain for generator PI torque controller [-]. (Only used in the transitional 2.5 region if VS_ControlMode =/ 2)
-4588245.18720      ! VS_KI				- Integral gain for generator PI torque controller [s]. (Only used in the transitional 2.5 region if VS_ControlMode =/ 2)
9.00                ! VS_TSRopt		    - Power-maximizing region 2 tip-speed-ratio. Only used in VS_ControlMode = 2.
0.314000000000      ! VS_PwrFiltF		- Low pass filter on power used to determine generator speed set point.  Only used in VS_ControlMode = 3.

!------- SETPOINT SMOOTHER ---------------------------------------------
1.00000             ! SS_VSGain         - Variable speed torque controller setpoint smoother gain, [-].
0.00100             ! SS_PCGain         - Collective pitch controller setpoint smoother gain, [-].

!------- POWER REFERENCE TRACKING --------------------------------------
2                   ! PRC_n			    -  Number of elements in PRC_WindSpeeds and PRC_GenSpeeds array
0.07854             ! PRC_LPF_Freq   - Frequency of the low pass filter on the wind speed estimate used to set PRC_GenSpeeds [rad/s]
3.0000 25.0000      ! PRC_WindSpeeds   - Array of wind speeds used in rotor speed vs. wind speed lookup table [m/s]
0.7917 0.7917       ! PRC_GenSpeeds   - Array of generator speeds corresponding to PRC_WindSpeeds [rad/s]

!------- WIND SPEED ESTIMATOR ---------------------------------------------
120.970             ! WE_BladeRadius	- Blade length (distance from hub center to blade tip), [m]
1                   ! WE_CP_n			- Amount of parameters in the Cp array
0.0                 ! WE_CP - Parameters that define the parameterized CP(lambda) function
0.0          		  ! WE_Gamma			- Adaption gain of the wind speed estimator algorithm [m/rad]
1.0                 ! WE_GearboxRatio	- Gearbox ratio [>=1],  [-]
312456272.00000     ! WE_Jtot			- Total drivetrain inertia, including blades, hub and casted generator inertia to LSS, [kg m^2]
1.225               ! WE_RhoAir			- Air density, [kg m^-3]
"../IEA-15-240-RWT/Cp_Ct_Cq.IEA15MW.txt"      ! PerfFileName      - File containing rotor performance tables (Cp,Ct,Cq) (absolute path or relative to this file)
36      26          ! PerfTableSize     - Size of rotor performance tables, first number refers to number of blade pitch angles, second number referse to number of tip-speed ratios
60                  ! WE_FOPoles_N      - Number of first-order system poles used in EKF
3.0000 3.2669 3.5338 3.8007 4.0676 4.3345 4.6014 4.8683 5.1352 5.4021 5.6690 5.9359 6.2028 6.4697 6.7366 7.0034 7.2703 7.5372 7.8041 8.0710 8.3379 8.6048 8.8717 9.1386 9.4055 9.6724 9.9393 10.2062 10.4731 10.7400 11.2153 11.6907 12.1660 12.6413 13.1167 13.5920 14.0673 14.5427 15.0180 15.4933 15.9687 16.4440 16.9193 17.3947 17.8700 18.3453 18.8207 19.2960 19.7713 20.2467 20.7220 21.1973 21.6727 22.1480 22.6233 23.0987 23.5740 24.0493 24.5247 25.0000               ! WE_FOPoles_v      - Wind speeds corresponding to first-order system poles [m/s]
-0.02567895 -0.02796349 -0.03024803 -0.03253257 -0.03481711 -0.03710166 -0.03938620 -0.04167074 -0.04395528 -0.04623982 -0.04852436 -0.05080890 -0.05309344 -0.05537798 -0.05766253 -0.05994707 -0.06223161 -0.06451615 -0.06680069 -0.06908523 -0.07136977 -0.07365431 -0.07593885 -0.07822339 -0.08050794 -0.08279248 -0.08507702 -0.08736156 -0.08964610 -0.09059663 -0.06049391 -0.06284948 -0.06999814 -0.07929908 -0.08986026 -0.10164081 -0.11414940 -0.12758646 -0.14144758 -0.15588074 -0.17113587 -0.18650122 -0.20251939 -0.21912926 -0.23583120 -0.25319804 -0.27095717 -0.28905917 -0.30746443 -0.32634101 -0.34558843 -0.36517204 -0.38493017 -0.40510931 -0.42563447 -0.44649211 -0.46766192 -0.48900392 -0.51050507 -0.53219472               ! WE_FOPoles        - First order system poles [1/s]

!------- YAW CONTROL ------------------------------------------------------
0.00000             ! Y_uSwitch		- Wind speed to switch between Y_ErrThresh. If zero, only the second value of Y_ErrThresh is used [m/s]
4.000000  8.000000  ! Y_ErrThresh    - Yaw error threshold/deadbands. Turbine begins to yaw when it passes this. If Y_uSwitch is zero, only the second value is used. [deg].
0.00870             ! Y_Rate			- Yaw rate [rad/s]
0.00000             ! Y_MErrSet		- Integrator saturation (maximum signal amplitude contribution to pitch from yaw-by-IPC), [rad]
0.00000             ! Y_IPC_IntSat		- Integrator saturation (maximum signal amplitude contribution to pitch from yaw-by-IPC), [rad]
0.00000             ! Y_IPC_KP			- Yaw-by-IPC proportional controller gain Kp
0.00000             ! Y_IPC_KI			- Yaw-by-IPC integral controller gain Ki

!------- TOWER CONTROL ------------------------------------------------------
0.00000             ! TRA_ExclSpeed	    - Rotor speed for exclusion [LSS, rad/s]
0.00000             ! TRA_ExclBand	    - Size of the rotor frequency exclusion band [LSS, rad/s]. Torque controller reference will be TRA_ExclSpeed +/- TRA_ExlBand/2
0.00000e+00         ! TRA_RateLimit	    - Rate limit of change in rotor speed reference [LSS, rad/s].  Suggested to be VS_RefSpd/400.
0.00000             ! FA_KI				- Integral gain for the fore-aft tower damper controller,  [rad*s/m]
0.00000             ! FA_HPFCornerFreq	- Corner frequency (-3dB point) in the high-pass filter on the fore-aft acceleration signal [rad/s]
0.00000             ! FA_IntSat			- Integrator saturation (maximum signal amplitude contribution to pitch from FA damper), [rad]

!------- MINIMUM PITCH SATURATION -------------------------------------------
60                  ! PS_BldPitchMin_N  - Number of values in minimum blade pitch lookup table (should equal number of values in PS_WindSpeeds and PS_BldPitchMin)
3.000 3.267 3.534 3.801 4.068 4.334 4.601 4.868 5.135 5.402 5.669 5.936 6.203 6.470 6.737 7.003 7.270 7.537 7.804 8.071 8.338 8.605 8.872 9.139 9.406 9.672 9.939 10.206 10.473 10.740 11.215 11.691 12.166 12.641 13.117 13.592 14.067 14.543 15.018 15.493 15.969 16.444 16.919 17.395 17.870 18.345 18.821 19.296 19.771 20.247 20.722 21.197 21.673 22.148 22.623 23.099 23.574 24.049 24.525 25.000               ! PS_WindSpeeds     - Wind speeds corresponding to minimum blade pitch angles [m/s]
0.060      0.060      0.060      0.060      0.060      0.060      0.056      0.052      0.047      0.041      0.036      0.029      0.022      0.015      0.008      0.001      0.000      0.000      0.000      0.000      0.000      0.000      0.000      0.000      0.000      0.006      0.021      0.033      0.044      0.053      0.064      0.074      0.084      0.093      0.103      0.112      0.121      0.130      0.138      0.147      0.155      0.163      0.172      0.180      0.188      0.196      0.203      0.211      0.219      0.227      0.234      0.242      0.250      0.257      0.265      0.272      0.279      0.287      0.294      0.301                    ! PS_BldPitchMin    - Minimum blade pitch angles [rad]

!------- SHUTDOWN -----------------------------------------------------------
0.698100000000      ! SD_MaxPit         - Maximum blade pitch angle to initiate shutdown, [rad]
0.418880000000      ! SD_CornerFreq     - Cutoff Frequency for first order low-pass filter for blade pitch angle, [rad/s]

!------- Floating -----------------------------------------------------------
1                   ! Fl_n              - Number of Fl_Kp gains in gain scheduling, optional with default of 1
-9.2561             ! Fl_Kp             - Nacelle pitching proportional feedback gain [s]
11.2770             ! Fl_U              - Wind speeds for scheduling Fl_Kp, optional if Fl_Kp is single value [m/s]

!------- FLAP ACTUATION -----------------------------------------------------
0.000000000000      ! Flp_Angle         - Initial or steady state flap angle [rad]
0.00000000e+00      ! Flp_Kp            - Blade root bending moment proportional gain for flap control [s]
0.00000000e+00      ! Flp_Ki            - Flap displacement integral gain for flap control [-]
0.174500000000      ! Flp_MaxPit        - Maximum (and minimum) flap pitch angle [rad]

!------- Open Loop Control -----------------------------------------------------
"unused"            ! OL_Filename       - Input file with open loop timeseries (absolute path or relative to this file)
0                   ! Ind_Breakpoint    - The column in OL_Filename that contains the breakpoint (time if OL_Mode = 1)
  0   0   0         ! Ind_BldPitch      - The columns in OL_Filename that contains the blade pitch (1,2,3) inputs in rad [array]
0                   ! Ind_GenTq         - The column in OL_Filename that contains the generator torque in Nm
0                   ! Ind_YawRate       - The column in OL_Filename that contains the yaw rate in rad/s
0                   ! Ind_Azimuth       - The column in OL_Filename that contains the desired azimuth position in rad (used if OL_Mode = 2)
0.0000 0.0000 0.0000 0.0000        ! RP_Gains - PID gains and Tf of derivative for rotor position control (used if OL_Mode = 2)
0                   ! Ind_CableControl  - The column(s) in OL_Filename that contains the cable control inputs in m [Used with CC_Mode = 2, must be the same size as CC_Group_N]
0                   ! Ind_StructControl - The column(s) in OL_Filename that contains the structural control inputs [Used with StC_Mode = 2, must be the same size as StC_Group_N]

!------- Pitch Actuator Model -----------------------------------------------------
3.140000000000       ! PA_CornerFreq     - Pitch actuator bandwidth/cut-off frequency [rad/s]
0.707000000000       ! PA_Damping        - Pitch actuator damping ratio [-, unused if PA_Mode = 1]

!------- Pitch Actuator Faults -----------------------------------------------------
0.00000000 0.00000000 0.00000000                 ! PF_Offsets     - Constant blade pitch offsets for blades 1-3 [rad]

!------- Active Wake Control -----------------------------------------------------
1                   ! AWC_NumModes       - Number of user-defined AWC forcing modes 
1                   ! AWC_n              - Azimuthal mode number(s) (i.e., the number and direction of the lobes of the wake structure)
1                   ! AWC_harmonic       - Harmonic(s) to apply in the AWC Inverse Coleman Transformation (only used when AWC_Mode = 2)
0.0500              ! AWC_freq           - Frequency(s) of forcing mode(s) [Hz]
1.0000              ! AWC_amp            - Pitch amplitude(s) of individual forcing mode(s) [deg]
0.0000              ! AWC_clockangle     - Initial angle(s) of forcing mode(s) [deg]


・If I want to lower f_c, what value should I set it to? Is it best to use the recommended value of 1.28 × average wind speed ÷ radius?

Thank you for your support.

Best regards,

Dear @Yuki.Ogura,

My understanding is that the very high TSR (above 20) leads to a high thrust that causes FAST.Farm to abort because the near-wake correction model is not valid for high thrust coefficients. My guess is that the high TSR results from winds that fall below the cut-in wind speed due to the wake effect on T2 from T1, but you’d have to plot other things (such as hub-height wind speed, the rotor seed, the pitch angle, etc. of T2) to verify. It may also help if you compare the results of T2 side-by-side with those from T1 to better understand what is happening.

Regarding f_c, how are you calculating this now? For a fixed-bottom wind turbine, we generally recommend setting f_c equal to 1.28*U0/R, but for floating systems, we recommend using a lower value, i.e., settings f_c to a frequency below the rigid-body natural frequencies to avoid having unreasonable streaks in the wake resulting from the floater motion interacting with FAST.Farm’s wake-evolution model.

Regarding ROSCO, are you asking how to set up ROSCO so that the turbine shuts down for wind conditions below cut-in?

Best regards,

Dear @Jason.Jonkman

Thank you for your reply.

I apologize for the delay in my response.

Regarding the value of f_c, I performed the calculation using 1.28 × 4 m/s / 120 m = 0.043, as recommended in the formula, but just like in my previous post, an error message appeared around 500 seconds into the simulation.

I have a question regarding ROSCO and how to configure it to shut down the turbine when wind speeds fall below the cut-in speed.
My understanding is that ROSCO reduces the TSR by feathering the blades to lower the rotational speed when wind speeds fall below the cut-in speed. However, as shown in the figure below, there is almost no blade pitch adjustment even when wind speeds are below the cut-in speed.
Therefore, I suspect there may be an issue with the ROSCO settings.

Thank you for your help.

Best regards,

Hi,

ROSCO does not have built-in logic to shut down at low wind speeds. There are a few ways it could be done; they all require some configuration and validation:

  1. Use power reference control and a create a table that reduces the PRC_R_Pitch at low wind speeds
  2. Use the power reference control mode that will set a very low generator speed set point at low wind speeds.
  3. Add elements to the peak shaving minimum pitch table at lower wind speeds than 3 m/s.
  4. Add a shutdown mode based on the wind speed and/or power threshold.

Option three is likely your lowest effort option that will get you past the possibility that this is causing your issue. Option 4 is the most realistic representation of in-field turbine operation.

I hope this helps.

Best, Dan.

Dear @Daniel.Zalkind and @Jason.Jonkman ,

Thank you for your detailed and thorough response.

@Daniel.Zalkind , thank you for the clarification. Now I understand that, as it stands, ROSCO does not have built-in logic to automatically shut down when wind speeds fall below the cut-in threshold.

I will try Option 3, which you recommended, to see if that resolves the issue.

@Jason.Jonkman, regarding the value of f_c, since it must be lower than the natural frequency of the rigid body, is my understanding correct that it should be set to a value lower than 0.007 Hz—the smallest natural frequency in the surge direction listed in Table 11 of this document?

Definition of the UMaine VolturnUS-S Reference Platform Developed for the IEA Wind 15- Megawatt Offshore Reference Wind Turbine

Thank you for your assistance.

Best regards,

Dear @Yuki.Ogura,

Yes, I agree.

Best regards,