NREL / ROSCO

A Reference Open Source Controller for Wind Turbines
https://rosco.readthedocs.io/en/latest/
Apache License 2.0
112 stars 92 forks source link

Error in running 02_ccblade.py #378

Closed GBA217 closed 1 month ago

GBA217 commented 1 month ago

I'm trying to tune the controller for a 80m diameter 2MW turbine and when I run 02_ccblade, I get this error:

(rosco-envir) c:\Users\xxxx\ROSCO\Examples>python 02_ccblade.py Using ofTools in rosco.toolbox...

Loading wind turbine data for NREL's ROSCO tuning and simulation processeses

Loading FAST model: Turb10ms.fst Loading rotor performace data from text file: c:\NREL-2p8-127_Cp_Ct_Cq.txt Loading rotor performance data from CC-Blade. Running CCBlade aerodynamic analysis, this may take a minute... UserWarning: C:\Users\xxxx\AppData\Local\anaconda3\envs\rosco-envir\lib\site-packages\wisdem\ccblade\ccblade.py:970 error. check input values.NaNs at 0/41: 1.5707963267948966 0.0 3.141591653589793 NaNs at 0/41: 1.5707963267948966 0.0 3.141591653589793 NaNs at 0/41: 1.5707963267948966 0.0 3.141591653589793 CCBlade aerodynamic analysis run successfully. Writing rotor performance text file: c:\Users\xxxx\ROSCO\Examples\examples_out\02_Cp_Ct_Cq.Ex03.txt

I have a few questions regarding this.

  1. If I have not generated a cpctcq curve before for the model, what do I give as input for "rotor_performance_filename" in the .yaml file? I provided the path to NREL-2p8-127_Cp_Ct_Cq.txt. Will that work? Or can I leave that input blank?
  2. Is aerofoil corrdinates file mandatory when I'm tuning th econtroller for the first time? It is mentioned that Numcoords can be set to zero if coordinates not included. Does that mean that it is not necessary for the simulation?
  3. What could cause the above mentioned error?

.yaml file that I used: path_params: {FAST_InputFile: Turb10ms.fst, FAST_directory: c:\01_OpenFAST_MM82, rotor_performance_filename: c:\01_OpenFAST_MM82\02_Cp_Ct_Cq.Ex03.txt} controller_params: LoggingLevel: 2 F_LPFType: 1 F_NotchType: 0 IPC_ControlMode: 1 VS_ControlMode: 3 PC_ControlMode: 1 Y_ControlMode: 0 SS_Mode: 1 WE_Mode: 2 PS_Mode: 1 SD_Mode: 0 Fl_Mode: 0 Flp_Mode: 0 zeta_pc: [2.0062413749856227] omega_pc: [0.13653659182284006] zeta_vs: 0.465 omega_vs: 0.11 twr_freq: 0.458 ptfm_freq: 0.2 ps_percent: 0.8 sd_maxpit: 0.4363 DISCON: {PC_MinPit: -0.1745329, PC_FinePit: 0.0, DLL_FileName: libdiscon, DLL_InFile: discon, DLL_ProcName: DISCON} TD_Mode: 0 PwC_Mode: 0 ZMQ_Mode: 0 ZMQ_UpdatePeriod: 2 PA_Mode: 0 Ext_Mode: 0 U_pc: [12] interp_type: sigma max_pitch: 1.57 min_pitch: 0.0 vs_minspd: 0.0 ss_vsgain: 1.0 ss_pcgain: 0.001 flp_maxpit: 0.1745 WS_GS_n: 60 PC_GS_n: 30 tune_Fl: true max_torque_factor: 1.2 IPC_Kp1p: 0.0 IPC_Kp2p: 0.0 IPC_Ki1p: 5.2e-9 IPC_Ki2p: 0.0 IPC_Vramp: [0.0, 0.0] filter_params: {f_we_cornerfreq: 0.20944, f_fl_highpassfreq: 0.01042, f_ss_cornerfreq: 0.6283, f_yawerr: 0.17952, f_sd_cornerfreq: 0.41888} open_loop: {flag: false, filename: unused, OL_Ind_Breakpoint: 1, OL_Ind_BldPitch: 0, OL_Ind_GenTq: 0, OL_Ind_YawRate: 0} PA_CornerFreq: 3.14 PA_Damping: 0.707 linmodel_tuning: type: none linfile_path: none lintune_outpath: lintune_outfiles load_parallel: false stability_margin: 0.1 omega_pc: {} flag: true tuning_yaml: /home/dzalkind/Tools/ROSCO/Tune_Cases/NREL2p8.yaml n_pitch: 30 n_tsr: 30 n_U: 1 omega_flp: 0.0 zeta_flp: 0.0 turbine_params: {rotor_inertia: 5375828.000, rated_rotor_speed: 1.7854, v_min: 3.0, v_rated: 8.4, v_max: 25.0, max_pitch_rate: 0.1745, max_torque_rate: 11000.0, rated_power: 2000000.0, bld_edgewise_freq: 1.663, bld_flapwise_freq: 0.905}

Thanks in advance!

dzalkind commented 1 month ago
  1. Any input for the performance_filename should work. This example will generate a new one.
  2. I don't think the airfoil coordinates are necessary.
  3. I'm not sure what causes this error. I will ask the WISDEM team.
ptrbortolotti commented 1 month ago

hello, the error is most likely caused by some issue with the input files. can you maybe share your setup?

GBA217 commented 1 month ago

Hello @ptrbortolotti, Yes, I figured that the issue is coming from the brent function calculation as the residual BEM solver residual is failing to reach zero in the iterative solver (please correct me if I'm wrong). Below is the setup:

Elastodyn: ------- ELASTODYN v1.03.* INPUT FILE ------------------------------------------- xxxx; GBA; 2024-02-29 ---------------------- SIMULATION CONTROL -------------------------------------- False Echo - Echo input data to ".ech" (flag) 3 Method - Integration method: {1: RK4, 2: AB4, or 3: ABM4} (-) "default" DT Integration time step (s) ---------------------- DEGREES OF FREEDOM -------------------------------------- True FlapDOF1 - First flapwise blade mode DOF (flag) True FlapDOF2 - Second flapwise blade mode DOF (flag) True EdgeDOF - First edgewise blade mode DOF (flag) False TeetDOF - Rotor-teeter DOF (flag) [unused for 3 blades] False DrTrDOF - Drivetrain rotational-flexibility DOF (flag) True GenDOF - Generator DOF (flag) False YawDOF - Yaw DOF (flag) False TwFADOF1 - First fore-aft tower bending-mode DOF (flag) False TwFADOF2 - Second fore-aft tower bending-mode DOF (flag) False TwSSDOF1 - First side-to-side tower bending-mode DOF (flag) False TwSSDOF2 - Second side-to-side tower bending-mode DOF (flag) False PtfmSgDOF - Platform horizontal surge translation DOF (flag) False PtfmSwDOF - Platform horizontal sway translation DOF (flag) False PtfmHvDOF - Platform vertical heave translation DOF (flag) False PtfmRDOF - Platform roll tilt rotation DOF (flag) False PtfmPDOF - Platform pitch tilt rotation DOF (flag) False PtfmYDOF - Platform yaw rotation DOF (flag) ---------------------- INITIAL CONDITIONS -------------------------------------- 0.0 OoPDefl - Initial out-of-plane blade-tip displacement (meters) 0.0 IPDefl - Initial in-plane blade-tip deflection (meters) 0.0 BlPitch(1) - Blade 1 initial pitch (degrees) 0.0 BlPitch(2) - Blade 2 initial pitch (degrees) 0.0 BlPitch(3) - Blade 3 initial pitch (degrees) [unused for 2 blades] 0.0 TeetDefl - Initial or fixed teeter angle (degrees) [unused for 3 blades] 0.0 Azimuth - Initial azimuth angle for blade 1 (degrees) 7.223281620282807 RotSpeed - Initial or fixed rotor speed (rpm) 0.0 NacYaw - Initial or fixed nacelle-yaw angle (degrees) 0.0 TTDspFA - Initial fore-aft tower-top displacement (meters) 0.0 TTDspSS - Initial side-to-side tower-top displacement (meters) 0.0 PtfmSurge - Initial or fixed horizontal surge translational displacement of platform (meters) 0.0 PtfmSway - Initial or fixed horizontal sway translational displacement of platform (meters) 0.0 PtfmHeave - Initial or fixed vertical heave translational displacement of platform (meters) 0.0 PtfmRoll - Initial or fixed roll tilt rotational displacement of platform (degrees) 0.0 PtfmPitch - Initial or fixed pitch tilt rotational displacement of platform (degrees) 0.0 PtfmYaw - Initial or fixed yaw rotational displacement of platform (degrees) ---------------------- TURBINE CONFIGURATION ----------------------------------- 3 NumBl - Number of blades (-) 41.1 TipRad - The distance from the rotor apex to the blade tip (meters) 1.1 HubRad - The distance from the rotor apex to the blade root (meters) -3.5 PreCone(1) - Blade 1 cone angle (degrees) -3.5 PreCone(2) - Blade 2 cone angle (degrees) -3.5 PreCone(3) - Blade 3 cone angle (degrees) [unused for 2 blades] 0.1 HubCM - Distance from rotor apex to hub mass [positive downwind] (meters) 0.0 UndSling - Undersling length [distance from teeter pin to the rotor apex] (meters) [unused for 3 blades] 0.0 Delta3 - Delta-3 angle for teetering rotors (degrees) [unused for 3 blades] 0.0 AzimB1Up - Azimuth value to use for I/O when blade 1 points up (degrees) -3.2 OverHang - Distance from yaw axis to rotor apex [3 blades] or teeter pin [2 blades] (meters) 1.1 ShftGagL - Distance from rotor apex [3 blades] or teeter pin [2 blades] to shaft strain gages [positive for upwind rotors] (meters) -5.0 ShftTilt - Rotor shaft tilt angle (degrees) 1.25 NacCMxn - Downwind distance from the tower-top to the nacelle CM (meters) 0.0 NacCMyn - Lateral distance from the tower-top to the nacelle CM (meters) 2.45 NacCMzn - Vertical distance from the tower-top to the nacelle CM (meters) 0.0 NcIMUxn - Downwind distance from the tower-top to the nacelle IMU (meters) 0.0 NcIMUyn - Lateral distance from the tower-top to the nacelle IMU (meters) 0.0 NcIMUzn - Vertical distance from the tower-top to the nacelle IMU (meters) 1.6 Twr2Shft - Vertical distance from the tower-top to the rotor shaft (meters) 56.65 TowerHt - Height of tower above ground level [onshore] or MSL [offshore] (meters) 0.0 TowerBsHt - Height of tower base above ground level [onshore] or MSL [offshore] (meters) 0.0 PtfmCMxt - Downwind distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters) 0.0 PtfmCMyt - Lateral distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters) 0.0 PtfmCMzt - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters) 0.0 PtfmRefzt - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform reference point (meters) ---------------------- MASS AND INERTIA ---------------------------------------- 0.0 TipMass(1) - Tip-brake mass, blade 1 (kg) 0.0 TipMass(2) - Tip-brake mass, blade 2 (kg) 0.0 TipMass(3) - Tip-brake mass, blade 3 (kg) [unused for 2 blades] 17000 HubMass - Hub mass (kg) 10000 HubIner - Hub inertia about rotor axis [3 blades] or teeter axis [2 blades] (kg m^2) 96 GenIner - Generator inertia about HSS (kg m^2) 66000 NacMass - Nacelle mass (kg) 6234818 NacYIner - Nacelle inertia about yaw axis (kg m^2) 0.0 YawBrMass - Yaw bearing mass (kg) 0.0 PtfmMass - Platform mass (kg) 0.0 PtfmRIner - Platform inertia for roll tilt rotation about the platform CM (kg m^2) 0.0 PtfmPIner - Platform inertia for pitch tilt rotation about the platform CM (kg m^2) 0.0 PtfmYIner - Platform inertia for yaw rotation about the platform CM (kg m^2) ---------------------- BLADE --------------------------------------------------- 41 BldNodes - Number of blade nodes (per blade) used for analysis (-) "ElastoDyn_blade.dat" BldFile1 - Name of file containing properties for blade 1 (quoted string) "ElastoDyn_blade.dat" BldFile2 - Name of file containing properties for blade 2 (quoted string) "ElastoDyn_blade.dat" BldFile3 - Name of file containing properties for blade 3 (quoted string) [unused for 2 blades] ---------------------- ROTOR-TEETER -------------------------------------------- 0 TeetMod - Rotor-teeter spring/damper model {0: none, 1: standard, 2: user-defined from routine UserTeet} (switch) [unused for 3 blades] 0.0 TeetDmpP - Rotor-teeter damper position (degrees) [used only for 2 blades and when TeetMod=1] 0.0 TeetDmp - Rotor-teeter damping constant (N-m/(rad/s)) [used only for 2 blades and when TeetMod=1] 0.0 TeetCDmp - Rotor-teeter rate-independent Coulomb-damping moment (N-m) [used only for 2 blades and when TeetMod=1] 0.0 TeetSStP - Rotor-teeter soft-stop position (degrees) [used only for 2 blades and when TeetMod=1] 0.0 TeetHStP - Rotor-teeter hard-stop position (degrees) [used only for 2 blades and when TeetMod=1] 0.0 TeetSSSp - Rotor-teeter soft-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1] 0.0 TeetHSSp - Rotor-teeter hard-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1] ---------------------- DRIVETRAIN ---------------------------------------------- 94 GBoxEff - Gearbox efficiency (%) 105.4 GBRatio - Gearbox ratio (-) 1.58E+08 DTTorSpr - Drivetrain torsional spring (N-m/rad) 1061776 DTTorDmp - Drivetrain torsional damper (N-m/(rad/s)) ---------------------- FURLING ------------------------------------------------- False Furling - Read in additional model properties for furling turbine (flag) [must currently be FALSE) "unused" FurlFile - Name of file containing furling properties (quoted string) [unused when Furling=False] ---------------------- TOWER --------------------------------------------------- 29 TwrNodes - Number of tower nodes used for analysis (-) "ElastoDyn_tower.dat" TwrFile - Name of file containing tower properties (quoted string) ---------------------- OUTPUT -------------------------------------------------- True SumPrint - Print summary data to ".sum" (flag) 1 OutFile - Switch to determine where output will be placed: {1: in module output file only; 2: in glue code output file only; 3: both} (currently unused) True TabDelim - Use tab delimiters in text tabular output file? (flag) (currently unused) "ES10.3E2" OutFmt - Format used for text tabular output (except time). Resulting field should be 10 characters. (quoted string) (currently unused) 50 TStart - Time to begin tabular output (s) (currently unused) 10 DecFact - Decimation factor for tabular output {1: output every time step} (-) (currently unused) 3 NTwGages - Number of tower nodes that have strain gages for output [0 to 9] (-) 1, 5, 29 TwrGagNd - List of tower nodes that have strain gages [1 to TwrNodes] (-) [unused if NTwGages=0] TwrGagNd - List of tower nodes that have strain gages [1 to TwrNodes] (-) [unused if NTwGages=0] TwrGagNd - List of tower nodes that have strain gages [1 to TwrNodes] (-) [unused if NTwGages=0] 4 NBlGages - Number of blade nodes that have strain gages for output [0 to 9] (-) 7, 15, 23, 31 BldGagNd - List of blade nodes that have strain gages [1 to BldNodes] (-) [unused if NBlGages=0] BldGagNd - List of blade nodes that have strain gages [1 to BldNodes] (-) [unused if NBlGages=0] BldGagNd - List of blade nodes that have strain gages [1 to BldNodes] (-) [unused if NBlGages=0] OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) "Azimuth" "BldPitch1" "BldPitch2" "BldPitch3" "GenSpeed" "RotPwr" "TipClrnc1 TipClrnc2 TipClrnc3" "RootMEdg1" "RootMEdg2" "RootMEdg3" "RootMFlp1" "RootMFlp2" "RootMFlp3" END of input file (the word "END" must appear in the first 3 columns of this last OutList line) ---------------------- NODE OUTPUTS -------------------------------------------- 1 BldNd_BladesOut - Blades to output 99 BldNd_BlOutNd - Blade nodes on each blade (currently unused) OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx, ElastoDyn_Nodes tab for a listing of available output channels, (-) END of input file (the word "END" must appear in the first 3 columns of this last OutList line)

AeroDyn15: ------- AERODYN v15.03. INPUT FILE ------------------------------------------------ Turbine, Onshore; GBA; 2024-03-05 ====== General Options ============================================================================ True Echo - Echo the input to ".AD.ech"? (flag) default DTAero - Time interval for aerodynamic calculations {or "default"} (s) 0 WakeMod - Type of wake/induction model (switch) {0=none, 1=BEMT, 2=DBEMT, 3=OLAF} [WakeMod cannot be 2 or 3 when linearizing] 1 AFAeroMod - Type of blade airfoil aerodynamics model (switch) {1=steady model, 2=Beddoes-Leishman unsteady model} [AFAeroMod must be 1 when linearizing] 0 TwrPotent - Type tower influence on wind based on potential flow around the tower (switch) {0=none, 1=baseline potential flow, 2=potential flow with Bak correction} 0 TwrShadow - Calculate tower influence on wind based on downstream tower shadow (switch) {0=none, 1=Powles model, 2=Eames model} False TwrAero - Calculate tower aerodynamic loads? (flag) False FrozenWake - Assume frozen wake during linearization? (flag) [used only when WakeMod=1 and when linearizing] False CavitCheck - Perform cavitation check? (flag) [AFAeroMod must be 1 when CavitCheck=true] False Buoyancy - Include buoyancy effects? (flag) False CompAA - Flag to compute AeroAcoustics calculation [only used when WakeMod=1 or 2] AeroAcousticsInput.dat AA_InputFile - AeroAcoustics input file [used only when CompAA=true] ====== Environmental Conditions =================================================================== "default" AirDens - Air density (kg/m^3) "default" KinVisc - Kinematic air viscosity (m^2/s) "default" SpdSound - Speed of sound (m/s) "default" Patm - Atmospheric pressure (Pa) [used only when CavitCheck=True] "default" Pvap - Vapour pressure of fluid (Pa) [used only when CavitCheck=True] ====== Blade-Element/Momentum Theory Options ====================================================== [used only when WakeMod=1] 1 SkewMod - Type of skewed-wake correction model (switch) {1=uncoupled, 2=Pitt/Peters, 3=coupled} [unused when WakeMod=0 or 3] default SkewModFactor - Constant used in Pitt/Peters skewed wake model {or "default" is 15/32pi} (-) [used only when SkewMod=2; unused when WakeMod=0 or 3] False TipLoss - Use the Prandtl tip-loss model? (flag) [unused when WakeMod=0 or 3] False HubLoss - Use the Prandtl hub-loss model? (flag) [unused when WakeMod=0 or 3] False TanInd - Include tangential induction in BEMT calculations? (flag) [unused when WakeMod=0 or 3] False AIDrag - Include the drag term in the axial-induction calculation? (flag) [unused when WakeMod=0 or 3] False TIDrag - Include the drag term in the tangential-induction calculation? (flag) [unused when WakeMod=0,3 or TanInd=FALSE] default IndToler - Convergence tolerance for BEMT nonlinear solve residual equation {or "default"} (-) [unused when WakeMod=0 or 3] 500 MaxIter - Maximum number of iteration steps (-) [unused when WakeMod=0] ====== Dynamic Blade-Element/Momentum Theory Options ====================================================== [used only when WakeMod=1] 2 DBEMT_Mod - Type of dynamic BEMT (DBEMT) model {1=constant tau1, 2=time-dependent tau1} (-) [used only when WakeMod=2] 4 tau1_const - Time constant for DBEMT (s) [used only when WakeMod=2 and DBEMT_Mod=1] ====== OLAF -- cOnvecting LAgrangian Filaments (Free Vortex Wake) Theory Options ================== [used only when WakeMod=3] ../IEA-15-240-RWT-OLAF/IEA-15-240-RWT_OLAF.dat OLAFInputFileName - Input file for OLAF [used only when WakeMod=3] ====== Beddoes-Leishman Unsteady Airfoil Aerodynamics Options ===================================== [used only when AFAeroMod=2] 2 UAMod - Unsteady Aero Model Switch (switch) {1=Baseline model (Original), 2=Gonzalez's variant (changes in Cn,Cc,Cm), 3=Minnema/Pierce variant (changes in Cc and Cm)} [used only when AFAeroMod=2] True FLookup - Flag to indicate whether a lookup for f' will be calculated (TRUE) or whether best-fit exponential equations will be used (FALSE); if FALSE S1-S4 must be provided in airfoil input files (flag) [used only when AFAeroMod=2] 0.1 UAStartRad - Starting radius for dynamic stall (fraction of rotor radius) [used only when AFAeroMod=2] 1.0 UAEndRad - Ending radius for dynamic stall (fraction of rotor radius) [used only when AFAeroMod=2] ====== Airfoil Information ========================================================================= 1 AFTabMod - Interpolation method for multiple airfoil tables {1=1D interpolation on AoA (first table only); 2=2D interpolation on AoA and Re; 3=2D interpolation on AoA and UserProp} (-) 1 InCol_Alfa - The column in the airfoil tables that contains the angle of attack (-) 2 InCol_Cl - The column in the airfoil tables that contains the lift coefficient (-) 3 InCol_Cd - The column in the airfoil tables that contains the drag coefficient (-) 4 InCol_Cm - The column in the airfoil tables that contains the pitching-moment coefficient; use zero if there is no Cm column (-) 0 InCol_Cpmin - The column in the airfoil tables that contains the Cpmin coefficient; use zero if there is no Cpmin column (-) 17 NumAFfiles - Number of airfoil files used (-)" "01_15p4.dat" AFNames - Airfoil file names (NumAFfiles lines) (quoted strings) "02_15p5.dat" "03_15p8.dat" "04_17p3.dat" "05_20p0.dat" "06_24p0.dat" "07_27p6.dat" "08_31p0.dat" "09_35p0.dat" "10_40p5.dat" "11_45p0.dat" "12_50p0.dat" "13_58p0.dat" "14_70p9.dat" "15_84p0.dat" "16_94p6.dat" "17_100p0.dat" ====== Rotor/Blade Properties ===================================================================== True UseBlCm - Include aerodynamic pitching moment in calculations? (flag) "AeroDyn15_blade.dat" ADBlFile(1) - Name of file containing distributed aerodynamic properties for Blade #1 (-) "AeroDyn15_blade.dat" ADBlFile(2) - Name of file containing distributed aerodynamic properties for Blade #2 (-) [unused if NumBl < 2] "AeroDyn15_blade.dat" ADBlFile(3) - Name of file containing distributed aerodynamic properties for Blade #3 (-) [unused if NumBl < 3] ====== Hub Properties ============================================================================== [used only when Buoyancy=True] 0.0 VolHub - Hub volume (m^3) 0.0 HubCenBx - Hub center of buoyancy x direction offset (m) ====== Nacelle Properties ========================================================================== [used only when Buoyancy=True] 0.0 VolNac - Nacelle volume (m^3) 0,0,0 NacCenB - Position of nacelle center of buoyancy from yaw bearing in nacelle coordinates (m) ====== Tail fin Aerodynamics ======================================================================== False TFinAero - Calculate tail fin aerodynamics model (flag) "unused" TFinFile - Input file for tail fin aerodynamics [used only when TFinAero=True] ====== Tower Influence and Aerodynamics ============================================================= [used only when TwrPotent/=0, TwrShadow/=0, or TwrAero=True] 29 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, or TwrAero=True] TwrElev TwrDiam TwrCd TwrTI TwrCb !TwrTI used only with TwrShadow=2, TwrCb used only with Buoyancy=True (m) (m) (-) (-) (-) 0.00 3.800 0.64 0.1 0.0 2.00 3.800 0.64 0.1 0.0 4.95 3.800 0.64 0.1 0.0 4.96 3.800 0.64 0.1 0.0 7.90 3.800 0.64 0.1 0.0 10.85 3.800 0.64 0.1 0.0 10.86 3.800 0.64 0.1 0.0 13.55 3.660 0.64 0.1 0.0 16.25 3.520 0.64 0.1 0.0 18.95 3.380 0.64 0.1 0.0 21.65 3.240 0.64 0.1 0.0 24.40 3.098 0.64 0.1 0.0 27.15 2.955 0.64 0.1 0.0 27.16 2.955 0.64 0.1 0.0 30.10 2.955 0.64 0.1 0.0 33.05 2.955 0.64 0.1 0.0 33.06 2.955 0.64 0.1 0.0 36.00 2.955 0.64 0.1 0.0 36.01 2.955 0.64 0.1 0.0 38.95 2.955 0.64 0.1 0.0 41.90 2.955 0.64 0.1 0.0 41.91 2.955 0.64 0.1 0.0 44.85 2.955 0.64 0.1 0.0 44.86 2.955 0.64 0.1 0.0 47.80 2.955 0.64 0.1 0.0 50.75 2.955 0.64 0.1 0.0 53.70 2.955 0.64 0.1 0.0 53.71 2.955 0.64 0.1 0.0 56.65 2.955 0.64 0.1 0.0 ====== Outputs ==================================================================================== False SumPrint - Generate a summary file listing input options and interpolated properties to ".AD.sum"? (flag) 9 NBlOuts - Number of blade node outputs [0 - 9] (-) 6, 11, 16, 21, 25, 30, 35, 39, 41 BlOutNd - Blade nodes whose values will be output (-) 0 NTwOuts - Number of tower node outputs [0 - 9] (-) 0 TwOutNd - Tower nodes whose values will be output (-) OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) "RtFldFxh" "RtFldFyh" "RtFldFzh" "RtFldMxh" "RtFldMyh" "RtFldMzh" "RtVAvgxh" "RtFldCp" "RtFldCt" "RtArea" "RtSpeed" "RtTSR" "B1N1Cl" "B1N2Cl" "B1N3Cl" "B1N4Cl" "B1N5Cl" "B1N6Cl" "B1N7Cl" "B1N8Cl" "B1N9Cl" "B2N1Cl" "B2N2Cl" "B2N3Cl" "B2N4Cl" "B2N5Cl" "B2N6Cl" "B2N7Cl" "B2N8Cl" "B2N9Cl" "B3N1Cl" "B3N2Cl" "B3N3Cl" "B3N4Cl" "B3N5Cl" "B3N6Cl" "B3N7Cl" "B3N8Cl" "B3N9Cl" END of input file (the word "END" must appear in the first 3 columns of this last OutList line) ====== Outputs for all blade stations (same ending as above for B1N1.... =========================== [optional section] 1 BldNd_BladesOut - Number of blades to output all node information at. Up to number of blades on turbine. (-) "All" BldNd_BlOutNd - Future feature will allow selecting a portion of the nodes to output. Not implemented yet. (-) OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) Vrel Alpha Vindx Vindy Cn Ct Fn Ft Fx Fy TnInd AxInd END (the word "END" must appear in the first 3 columns of this last OutList line in the optional nodal output section)

ServoDyn: ------- SERVODYN v1.05.* INPUT FILE -------------------------------------------- Generated with AeroElasticSE FAST driver ---------------------- SIMULATION CONTROL -------------------------------------- False Echo - Echo input data to .ech (flag) default DT - Communication interval for controllers (s) (or "default") ---------------------- PITCH CONTROL ------------------------------------------- 5 PCMode - Pitch control mode {0: none, 3: user-defined from routine PitchCntrl, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) 0.0 TPCOn - Time to enable active pitch control (s) [unused when PCMode=0] 99999.0 TPitManS(1) - Time to start override pitch maneuver for blade 1 and end standard pitch control (s) 99999.0 TPitManS(2) - Time to start override pitch maneuver for blade 2 and end standard pitch control (s) 99999.0 TPitManS(3) - Time to start override pitch maneuver for blade 3 and end standard pitch control (s) [unused for 2 blades] 7.0 PitManRat(1) - Pitch rate at which override pitch maneuver heads toward final pitch angle for blade 1 (deg/s) 7.0 PitManRat(2) - Pitch rate at which override pitch maneuver heads toward final pitch angle for blade 2 (deg/s) 7.0 PitManRat(3) - Pitch rate at which override pitch maneuver heads toward final pitch angle for blade 3 (deg/s) [unused for 2 blades] 90.0 BlPitchF(1) - Blade 1 final pitch for pitch maneuvers (degrees) 90.0 BlPitchF(2) - Blade 2 final pitch for pitch maneuvers (degrees) 90.0 BlPitchF(3) - Blade 3 final pitch for pitch maneuvers (degrees) [unused for 2 blades] ---------------------- GENERATOR AND TORQUE CONTROL ---------------------------- 5 VSContrl - Variable-speed control mode {0: none, 1: simple VS, 3: user-defined from routine UserVSCont, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) 1 GenModel - Generator model {1: simple, 2: Thevenin, 3: user-defined from routine UserGen} (switch) [used only when VSContrl=0] 93.94772808556249 GenEff - Generator efficiency [ignored by the Thevenin and user-defined generator models] (%) True GenTiStr - Method to start the generator {T: timed using TimGenOn, F: generator speed using SpdGenOn} (flag) True GenTiStp - Method to stop the generator {T: timed using TimGenOf, F: when generator power = 0} (flag) 99999.0 SpdGenOn - Generator speed to turn on the generator for a startup (HSS speed) (rpm) [used only when GenTiStr=False] 0.0 TimGenOn - Time to turn on the generator for a startup (s) [used only when GenTiStr=True] 99999.0 TimGenOf - Time to turn off the generator (s) [used only when GenTiStp=True] ---------------------- SIMPLE VARIABLE-SPEED TORQUE CONTROL -------------------- 99999.0 VS_RtGnSp - Rated generator speed for simple variable-speed generator control (HSS side) (rpm) [used only when VSContrl=1] 99999.0 VS_RtTq - Rated generator torque/constant generator torque in Region 3 for simple variable-speed generator control (HSS side) (N-m) [used only when VSContrl=1] 99999.0 VS_Rgn2K - Generator torque constant in Region 2 for simple variable-speed generator control (HSS side) (N-m/rpm^2) [used only when VSContrl=1] 99999.0 VS_SlPc - Rated generator slip percentage in Region 2 1/2 for simple variable-speed generator control (%) [used only when VSContrl=1] ---------------------- SIMPLE INDUCTION GENERATOR ------------------------------ 99999.0 SIG_SlPc - Rated generator slip percentage (%) [used only when VSContrl=0 and GenModel=1] 99999.0 SIG_SySp - Synchronous (zero-torque) generator speed (rpm) [used only when VSContrl=0 and GenModel=1] 99999.0 SIG_RtTq - Rated torque (N-m) [used only when VSContrl=0 and GenModel=1] 99999.0 SIG_PORt - Pull-out ratio (Tpullout/Trated) (-) [used only when VSContrl=0 and GenModel=1] ---------------------- THEVENIN-EQUIVALENT INDUCTION GENERATOR ----------------- 99999.0 TEC_Freq - Line frequency [50 or 60] (Hz) [used only when VSContrl=0 and GenModel=2] 0 TEC_NPol - Number of poles [even integer > 0] (-) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_SRes - Stator resistance (ohms) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_RRes - Rotor resistance (ohms) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_VLL - Line-to-line RMS voltage (volts) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_SLR - Stator leakage reactance (ohms) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_RLR - Rotor leakage reactance (ohms) [used only when VSContrl=0 and GenModel=2] 99999.0 TEC_MR - Magnetizing reactance (ohms) [used only when VSContrl=0 and GenModel=2] ---------------------- HIGH-SPEED SHAFT BRAKE ---------------------------------- 0 HSSBrMode - HSS brake model {0: none, 1: simple, 3: user-defined from routine UserHSSBr, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) 99999.0 THSSBrDp - Time to initiate deployment of the HSS brake (s) 99999.0 HSSBrDT - Time for HSS-brake to reach full deployment once initiated (sec) [used only when HSSBrMode=1] 99999.0 HSSBrTqF - Fully deployed HSS-brake torque (N-m) ---------------------- NACELLE-YAW CONTROL ------------------------------------- 0 YCMode - Yaw control mode {0: none, 3: user-defined from routine UserYawCont, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) 99999.0 TYCOn - Time to enable active yaw control (s) [unused when YCMode=0] 0.0 YawNeut - Neutral yaw position--yaw spring force is zero at this yaw (degrees) 496739375.65069634 YawSpr - Nacelle-yaw spring constant (N-m/rad) 187442.5509148155 YawDamp - Nacelle-yaw damping constant (N-m/(rad/s)) 99999.0 TYawManS - Time to start override yaw maneuver and end standard yaw control (s) 0.25 YawManRat - Yaw maneuver rate (in absolute value) (deg/s) 0.0 NacYawF - Final yaw angle for override yaw maneuvers (degrees) ---------------------- AERODYNAMIC FLOW CONTROL -------------------------------- 0 AfCmode - Airfoil control mode {0: none, 1: cosine wave cycle, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) 0 AfC_Mean - Mean level for cosine cycling or steady value (-) [used only with AfCmode==1] 0 AfC_Amp - Amplitude for for cosine cycling of flap signal (-) [used only with AfCmode==1] 0 AfC_Phase - Phase relative to the blade azimuth (0 is vertical) for for cosine cycling of flap signal (deg) [used only with AfCmode==1] ---------------------- STRUCTURAL CONTROL -------------------------------------- 0 NumBStC - Number of blade structural controllers (integer) "unused" BStCfiles - Name of the files for blade structural controllers (quoted strings) [unused when NumBStC==0] 0 NumNStC - Number of nacelle structural controllers (integer) "unused" NStCfiles - Name of the files for nacelle structural controllers (quoted strings) [unused when NumNStC==0] 0 NumTStC - Number of tower structural controllers (integer) "unused" TStCfiles - Name of the files for tower structural controllers (quoted strings) [unused when NumTStC==0] 0 NumSStC - Number of substructure structural controllers (integer) "unused" SStCfiles - Name of the files for substructure structural controllers (quoted strings) [unused when NumSStC==0] ---------------------- CABLE CONTROL ------------------------------------------- 0 CCmode - Cable control mode {0: none, 4: user-defined from Simulink/Labview, 5: user-defined from Bladed-style DLL} (switch) ---------------------- BLADED INTERFACE ---------------------------------------- [used only with Bladed Interface] "libdiscon.dll" DLL_FileName - Name/location of the dynamic library {.dll [Windows] or .so [Linux]} in the Bladed-DLL format (-) [used only with Bladed Interface] "DISCON.IN" DLL_InFile - Name of input file sent to the DLL (-) [used only with Bladed Interface] "DISCON" DLL_ProcName - Name of procedure in DLL to be called (-) [case sensitive; used only with DLL Interface] default DLL_DT - Communication interval for dynamic library (s) (or "default") [used only with Bladed Interface] False DLL_Ramp - Whether a linear ramp should be used between DLL_DT time steps [introduces time shift when true] (flag) [used only with Bladed Interface] 99999.0 BPCutoff - Cuttoff frequency for low-pass filter on blade pitch from DLL (Hz) [used only with Bladed Interface] 0.0 NacYaw_North - Reference yaw angle of the nacelle when the upwind end points due North (deg) [used only with Bladed Interface] 1 Ptch_Cntrl - Record 28: Use individual pitch control {0: collective pitch; 1: individual pitch control} (switch) [used only with Bladed Interface] 0.0 Ptch_SetPnt - Record 5: Below-rated pitch angle set-point (deg) [used only with Bladed Interface] 0.0 Ptch_Min - Record 6: Minimum pitch angle (deg) [used only with Bladed Interface] 0.0 Ptch_Max - Record 7: Maximum pitch angle (deg) [used only with Bladed Interface] 0.0 PtchRate_Min - Record 8: Minimum pitch rate (most negative value allowed) (deg/s) [used only with Bladed Interface] 0.0 PtchRate_Max - Record 9: Maximum pitch rate (deg/s) [used only with Bladed Interface] 0.0 Gain_OM - Record 16: Optimal mode gain (Nm/(rad/s)^2) [used only with Bladed Interface] 0.0 GenSpd_MinOM - Record 17: Minimum generator speed (rpm) [used only with Bladed Interface] 0.0 GenSpd_MaxOM - Record 18: Optimal mode maximum speed (rpm) [used only with Bladed Interface] 0.0 GenSpd_Dem - Record 19: Demanded generator speed above rated (rpm) [used only with Bladed Interface] 0.0 GenTrq_Dem - Record 22: Demanded generator torque above rated (Nm) [used only with Bladed Interface] 0.0 GenPwr_Dem - Record 13: Demanded power (W) [used only with Bladed Interface] ---------------------- BLADED INTERFACE TORQUE-SPEED LOOK-UP TABLE ------------- 0 DLL_NumTrq - Record 26: No. of points in torque-speed look-up table {0 = none and use the optimal mode parameters; nonzero = ignore the optimal mode PARAMETERs by setting Record 16 to 0.0} (-) [used only with Bladed Interface] GenSpd_TLU GenTrq_TLU
(rpm) (Nm)
---------------------- OUTPUT -------------------------------------------------- False SumPrint - Print summary data to .sum (flag) (currently unused) 1 OutFile - Switch to determine where output will be placed: {1: in module output file only; 2: in glue code output file only; 3: both} (currently unused) True TabDelim - Use tab delimiters in text tabular output file? (flag) (currently unused) "ES10.3E2" OutFmt - Format used for text tabular output (except time). Resulting field should be 10 characters. (quoted string) (currently unused) 0.0 TStart - Time to begin tabular output (s) (currently unused) OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) "GenPwr" "GenTq" END of input file (the word "END" must appear in the first 3 columns of this last OutList line)

DISCON.IN:

! Controller parameter input file for the wind turbine ! - File written using ROSCO version 2.8.0 controller tuning logic on 01/05/24

!------- SIMULATION CONTROL ------------------------------------------------------------ 2 ! 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 .echo)

!------- CONTROLLER FLAGS ------------------------------------------------- 1 ! F_LPFType - (1: first-order low-pass filter, 2: second-order low-pass filter), [rad/s] (currently filters generator speed and pitch control signals 1 ! IPC_ControlMode - Turn Individual Pitch Control (IPC) for fatigue load reductions (pitch contribution) {0: off, 1: 1P reductions, 2: 1P+2P reductions} 3 ! 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} 1 ! PS_Mode - Pitch saturation mode {0: no pitch saturation, 1: implement pitch saturation} 0 ! SD_Mode - Shutdown mode {0: no shutdown procedure, 1: pitch to max pitch at shutdown} 0 ! 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} 0 ! 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 ---------------------------------------------------------- 2.07080 ! F_LPFCornerFreq - Corner frequency (-3dB point) in the low-pass filters, [rad/s] 0.00000 ! F_LPFDamping - Damping coefficient {used only when F_FilterType = 2} [-] 0 ! F_NumNotchFilts - Number of notch filters placed on sensors 0.0000 ! 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.0000 ! 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 0 ! F_TwrTopNotch_N - Number of notch filters on tower top acceleration signal 0 ! 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.000000 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]. 0.0000 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.075990 0.104503 0.127446 0.147544 0.164940 0.181678 0.197073 0.211624 0.225753 0.238993 0.252021 0.264458 0.276515 0.288229 0.299596 0.310932 0.321646 0.332112 0.342584 0.353105 0.363012 0.373072 0.382733 0.392195 0.401619 0.411107 0.420504 0.429096 0.437752 0.446357 ! PC_GS_angles - Gain-schedule table: pitch angles [rad]. -0.023307 -0.020910 -0.018908 -0.017210 -0.015752 -0.014487 -0.013378 -0.012398 -0.011527 -0.010746 -0.010043 -0.009406 -0.008827 -0.008298 -0.007813 -0.007366 -0.006953 -0.006571 -0.006216 -0.005885 -0.005577 -0.005288 -0.005017 -0.004762 -0.004522 -0.004296 -0.004083 -0.003881 -0.003689 -0.003508 ! PC_GS_KP - Gain-schedule table: pitch controller kp gains [s]. -0.000804 -0.000732 -0.000672 -0.000621 -0.000577 -0.000539 -0.000505 -0.000476 -0.000450 -0.000426 -0.000405 -0.000386 -0.000369 -0.000353 -0.000338 -0.000325 -0.000312 -0.000301 -0.000290 -0.000280 -0.000271 -0.000262 -0.000254 -0.000246 -0.000239 -0.000232 -0.000226 -0.000220 -0.000214 -0.000209 ! 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.17453000000 ! PC_MinPit - Minimum physical pitch limit, [rad]. 0.174500000000 ! PC_MaxRat - Maximum pitch rate (in absolute value) in pitch controller, [rad/s]. -0.17450000000 ! PC_MinRat - Minimum pitch rate (in absolute value) in pitch controller, [rad/s]. 188.5096700000 ! 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 ----------------------------------------- 9.120000 11.400000 ! 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, [-] 5.200e-09 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 ------------------------------------------------ 93.94773000000 ! VS_GenEff - Generator efficiency mechanical power -> electrical power, [should match the efficiency defined in the generator properties!], [%] 24248.54567000 ! VS_ArSatTq - Above rated generator torque PI control saturation, [Nm] 1500000.000000 ! VS_MaxRat - Maximum torque rate (in absolute value) in torque controller, [Nm/s]. 26673.40024000 ! VS_MaxTq - Maximum generator torque in Region 3 (HSS side), [Nm]. 0.000000000000 ! VS_MinTq - Minimum generator torque (HSS side), [Nm]. 80.81562000000 ! VS_MinOMSpd - Minimum generator speed [rad/s] 1.654680000000 ! VS_Rgn2K - Generator torque constant in Region 2 (HSS side). Only used in VS_ControlMode = 1,3 2000000.000000 ! VS_RtPwr - Wind turbine rated power [W] 24248.54567000 ! VS_RtTq - Rated torque, [Nm]. 188.5096700000 ! VS_RefSpd - Rated generator speed [rad/s] 1 ! VS_n - Number of generator PI torque controller gains -600.450450000 ! VS_KP - Proportional gain for generator PI torque controller [-]. (Only used in the transitional 2.5 region if VS_ControlMode =/ 2) -85.3230300000 ! VS_KI - Integral gain for generator PI torque controller [s]. (Only used in the transitional 2.5 region if VS_ControlMode =/ 2) 8.25 ! 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 --------------------------------------------- 63.457 ! 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] 97.0 ! WE_GearboxRatio - Gearbox ratio [>=1], [-] 66347470.49793 ! 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] "NREL-2p8-127_Cp_Ct_Cq.txt" ! PerfFileName - File containing rotor performance tables (Cp,Ct,Cq) (absolute path or relative to this file) 30 30 ! 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.2897 3.5793 3.8690 4.1586 4.4483 4.7379 5.0276 5.3172 5.6069 5.8966 6.1862 6.4759 6.7655 7.0552 7.3448 7.6345 7.9241 8.2138 8.5034 8.7931 9.0828 9.3724 9.6621 9.9517 10.2414 10.5310 10.8207 11.1103 11.4000 11.8533 12.3067 12.7600 13.2133 13.6667 14.1200 14.5733 15.0267 15.4800 15.9333 16.3867 16.8400 17.2933 17.7467 18.2000 18.6533 19.1067 19.5600 20.0133 20.4667 20.9200 21.3733 21.8267 22.2800 22.7333 23.1867 23.6400 24.0933 24.5467 25.0000 ! WE_FOPoles_v - Wind speeds corresponding to first-order system poles [m/s] -0.01246409 -0.01366752 -0.01487095 -0.01607438 -0.01727781 -0.01848123 -0.01968466 -0.02088809 -0.02209152 -0.02329495 -0.02449838 -0.02570181 -0.02690524 -0.02810867 -0.02931210 -0.03051553 -0.03171896 -0.03292239 -0.03412582 -0.03532924 -0.03653267 -0.03773610 -0.03893953 -0.04014296 -0.03876534 -0.03527029 -0.03014051 -0.02417865 -0.01714772 -0.00815401 -0.02209033 -0.02677213 -0.03094759 -0.03700125 -0.04188805 -0.04820630 -0.05390840 -0.06043728 -0.06624359 -0.07295380 -0.07986166 -0.08646363 -0.09412176 -0.10174417 -0.10913837 -0.11728247 -0.12497467 -0.13271621 -0.14189107 -0.15083268 -0.15981675 -0.16948290 -0.17870908 -0.18789884 -0.19750668 -0.20786563 -0.21788705 -0.22660607 -0.23612135 -0.24537696 ! 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.290 3.579 3.869 4.159 4.448 4.738 5.028 5.317 5.607 5.897 6.186 6.476 6.766 7.055 7.345 7.634 7.924 8.214 8.503 8.793 9.083 9.372 9.662 9.952 10.241 10.531 10.821 11.110 11.400 11.853 12.307 12.760 13.213 13.667 14.120 14.573 15.027 15.480 15.933 16.387 16.840 17.293 17.747 18.200 18.653 19.107 19.560 20.013 20.467 20.920 21.373 21.827 22.280 22.733 23.187 23.640 24.093 24.547 25.000 ! PS_WindSpeeds - Wind speeds corresponding to minimum blade pitch angles [m/s] 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.019 0.031 0.041 0.050 0.058 0.066 0.074 0.086 0.096 0.107 0.117 0.127 0.137 0.147 0.156 0.165 0.174 0.183 0.192 0.201 0.210 0.218 0.227 0.236 0.244 0.253 0.262 0.270 0.279 0.287 0.295 0.303 0.312 0.320 0.328 0.337 0.345 ! PS_BldPitchMin - Minimum blade pitch angles [rad]

!------- SHUTDOWN ----------------------------------------------------------- 0.436300000000 ! 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 0.0000 ! Fl_Kp - Nacelle velocity proportional feedback gain [s] 0.0000 ! 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]

!------- External Controller Interface ----------------------------------------------------- "unused" ! DLL_FileName - Name/location of the dynamic library in the Bladed-DLL format "unused" ! DLL_InFile - Name of input file sent to the DLL (-) "DISCON" ! DLL_ProcName - Name of procedure in DLL to be called (-)

!------- ZeroMQ Interface --------------------------------------------------------- "tcp://localhost:5555" ! ZMQ_CommAddress - Communication address for ZMQ server, (e.g. "tcp://localhost:5555") 1.000000 ! ZMQ_UpdatePeriod - Update period at zmq interface to send measurements and wait for setpoint [sec.] 0 ! ZMQ_ID - Integer identifier of turbine

!------- Cable Control --------------------------------------------------------- 1 ! CC_Group_N - Number of cable control groups 0 ! CC_GroupIndex - First index for cable control group, should correspond to deltaL 20.000000 ! CC_ActTau - Time constant for line actuator [s]

!------- Structural Controllers --------------------------------------------------------- 1 ! StC_Group_N - Number of cable control groups 0 ! StC_GroupIndex - First index for structural control group, options specified in ServoDyn summary output

It is mostly derived from NREL2p8MW turbine setup as I'm trying to modify it to a smaller rotor with 2MW rated power.

ptrbortolotti commented 1 month ago

Hello, I'm sorry but I don't really understand what's going on from this text. I talked to @dzalkind and our recommendation is to go back to a working example and start changing things one by one. It looks like you've changed the airfoil polars. Maybe start there, as things could easily go wrong in those tables. Please let us know if and when you find the issue

GBA217 commented 1 month ago

Hello @ptrbortolotti, the issue was indeed with the polars list. The order was messed up in the aerodyn input file. Thanks