I have made sure that all the parameters are the same between both AeroDyn (stand alone version) and AeroDyn in OpenFAST. But still I see significant difference between the Aerodynamic Power calculated using standalone Aerodyn vs OpenFAST.
Also, the AeroDyn input file of OpenFAST (V2.40) require some extra inputs like âAA_InputFileâ, âDBEMT_Modâ,âtau1_constâ, which are not present in the standalone AeroDyn input file. This may be because of the version difference. But I made sure that, the options I chose are agnostic to these differences, unless modifications are made to the way calculations are performed even with the same settings between the two (like âAFAeroModâ=1).
Currently I am using stand alone version of âAeroDyn_v15.04.02â. If there is latest version, please let me know where to download it from.
I have attached relevant files. If it is not a big ask, could you please go through the attachment and and let me know if I am making any mistakes? Thanks.
I did a quick comparison between your OpenFAST model and standalone AeroDyn model and donât see any obvious differences that would lead to different power.
Iâm guessing the problem is that you are using two different versions of AeroDyn between the cases. You mention that you are using OpenFAST v2.4 and standalone AeroDyn v15.04. While we donât provide precompiled Windows executable of the standalone modules (like AeroDyn) with each OpenFAST release, the source code needed to compile the standalone AeroDyn driver is included in the OpenFAST v2.4 repository. The AeroDyn driver input file has not changed between AeroDyn v15.04 and OpenFAST v2.4. So, all you should have to do is compile the standalone AeroDyn driver for OpenFAST v2.4, use your AeroDyn driver file together with the AeroDyn input file(s) from your OpenFAST v2.4 model, and rerun the driver-level simulation.
Dear Sina Sir,
I am also facing the same error while running the file Run_OpenLoop.m
Error using Run_OpenLoop (line 30)
Error reported by S-function âFAST_SFuncâ in âOpenLoop/FAST Nonlinear Wind Turbine/S-Functionâ:
FAST_InitializeAll:FAST_Init:FAST_ReadPrimaryFile:OpenFInpFile:The input file, ââŠ..\CertTest\Test01.fstâ, was not found.
Please guide how to solve this.
Thanks and regards,
Shivaji
This error means that the input file âTest01.fstâ does not exist, i.e., the filename is wrong or the path to it is incorrect.
As I told Sina, "In the transition from FAST v8 to OpenFAST, what used to be the CertTest of FAST v8 has been changed to the regression tests (r-test) of OpenFAST: github.com/OpenFAST/r-test. The primary input file names have also been changed.
I wanted to look at the change of the power coefficient Cp of a waked wind turbine. Thus, in FAST.Farm, I considered 2 wind turbines, one 4*D upstream the second one (D is the rotor diameter), and I plot the Tip-Speed Ratio (RtTSR), the rotor speed (RtSpeed), the averaged-disk wind speed (RtVAvgxh), the power coefficient (RtAeroCp), and the output power of the turbines.
I am surprised to have a Cp of the waked wind turbine increasing higher than the Cp_max =0.48, even though I made sure that the controllers are on (pitch & torque), and WakeMod=1. I checked also TSR and RtSpeed, and they seem to be the ones I expected (referring to the plots in the 5MW NREL Definition)
I was expecting a decrease of Cp for the second turbine as the TSR increases (going from Region 2 to Region 1/2). Do you have any idea where this issue may come from?
Are you running a FAST.Farm simulation using OpenFAST models of the NREL 5-MW baseline wind turbine or something else? What inflow conditions are you simulating in FAST.Farm (presumably steady, uniform inflow based on the smooth curves in your plots)? What structural DOFs and aerodynamic settings are enabled in your OpenFAST models?
I would generally expect the NREL 5-MW baseline wind turbine to have a Cpmax of 0.482 at a TSR of 7.55, but that was computed a while ago and different aerodynamic settings in AeroDyn v15 or enabling/disabling structural DOFs may have some effect on that. The peak of the Cp-TSR surface is also fairly flat for small changes in TSR, so, I could see that small variations in TSR will may have a small effect on Cp. However, a Cp of 0.54 sounds a bit higher than I would expect.
Yes, I am using the NREL 5MW baseline. Otherwise, where can I download its files again?
Yes, steady of 10m/s.
For the structural DOFs, only the Generator is enabled as can be seen here:
[code]------- ELASTODYN v1.03.* INPUT FILE -------------------------------------------
NREL 5.0 MW Baseline Wind Turbine for Use in Offshore Analysis. Properties from Dutch Offshore Wind Energy Converter (DOWEC) 6MW Pre-Design (10046_009.pdf) and REpower 5M 5MW (5m_uk.pdf)
---------------------- 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)
---------------------- ENVIRONMENTAL CONDITION ---------------------------------
9.80665 Gravity - Gravitational acceleration (m/s^2)
---------------------- DEGREES OF FREEDOM --------------------------------------
False !JASON:True FlapDOF1 - First flapwise blade mode DOF (flag)
False !JASON:True FlapDOF2 - Second flapwise blade mode DOF (flag)
False !JASON:True EdgeDOF - First edgewise blade mode DOF (flag)
False TeetDOF - Rotor-teeter DOF (flag) [unused for 3 blades]
False !JASON:True DrTrDOF - Drivetrain rotational-flexibility DOF (flag)
True GenDOF - Generator DOF (flag)
False !JASON:True YawDOF - Yaw DOF (flag)
False !JASON:True TwFADOF1 - First fore-aft tower bending-mode DOF (flag)
False !JASON:True TwFADOF2 - Second fore-aft tower bending-mode DOF (flag)
False !JASON:True TwSSDOF1 - First side-to-side tower bending-mode DOF (flag)
False !JASON:True 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 OoPDefl - Initial out-of-plane blade-tip displacement (meters)
0 IPDefl - Initial in-plane blade-tip deflection (meters)
0 BlPitch(1) - Blade 1 initial pitch (degrees)
0 BlPitch(2) - Blade 2 initial pitch (degrees)
0 BlPitch(3) - Blade 3 initial pitch (degrees) [unused for 2 blades]
0 TeetDefl - Initial or fixed teeter angle (degrees) [unused for 3 blades]
0 Azimuth - Initial azimuth angle for blade 1 (degrees)
11.0 !JASON: 12.1 RotSpeed - Initial or fixed rotor speed (rpm)
0 NacYaw - Initial or fixed nacelle-yaw angle (degrees)
0 TTDspFA - Initial fore-aft tower-top displacement (meters)
0 TTDspSS - Initial side-to-side tower-top displacement (meters)
0 PtfmSurge - Initial or fixed horizontal surge translational displacement of platform (meters)
0 PtfmSway - Initial or fixed horizontal sway translational displacement of platform (meters)
0 PtfmHeave - Initial or fixed vertical heave translational displacement of platform (meters)
0 PtfmRoll - Initial or fixed roll tilt rotational displacement of platform (degrees)
0 PtfmPitch - Initial or fixed pitch tilt rotational displacement of platform (degrees)
0 PtfmYaw - Initial or fixed yaw rotational displacement of platform (degrees)
---------------------- TURBINE CONFIGURATION -----------------------------------
3 NumBl - Number of blades (-)
63 TipRad - The distance from the rotor apex to the blade tip (meters)
1.5 HubRad - The distance from the rotor apex to the blade root (meters)
-2.5 PreCone(1) - Blade 1 cone angle (degrees)
-2.5 PreCone(2) - Blade 2 cone angle (degrees)
-2.5 PreCone(3) - Blade 3 cone angle (degrees) [unused for 2 blades]
0 HubCM - Distance from rotor apex to hub mass [positive downwind] (meters)
0 UndSling - Undersling length [distance from teeter pin to the rotor apex] (meters) [unused for 3 blades]
0 Delta3 - Delta-3 angle for teetering rotors (degrees) [unused for 3 blades]
0 AzimB1Up - Azimuth value to use for I/O when blade 1 points up (degrees)
-5.0191 OverHang - Distance from yaw axis to rotor apex [3 blades] or teeter pin [2 blades] (meters)
1.912 ShftGagL - Distance from rotor apex [3 blades] or teeter pin [2 blades] to shaft strain gages [positive for upwind rotors] (meters)
-5 ShftTilt - Rotor shaft tilt angle (degrees)
1.9 NacCMxn - Downwind distance from the tower-top to the nacelle CM (meters)
0 NacCMyn - Lateral distance from the tower-top to the nacelle CM (meters)
1.75 NacCMzn - Vertical distance from the tower-top to the nacelle CM (meters)
-3.09528 NcIMUxn - Downwind distance from the tower-top to the nacelle IMU (meters)
0 NcIMUyn - Lateral distance from the tower-top to the nacelle IMU (meters)
2.23336 NcIMUzn - Vertical distance from the tower-top to the nacelle IMU (meters)
1.96256 Twr2Shft - Vertical distance from the tower-top to the rotor shaft (meters)
87.6 TowerHt - Height of tower above ground level [onshore] or MSL [offshore] (meters)
0 TowerBsHt - Height of tower base above ground level [onshore] or MSL [offshore] (meters)
0 PtfmCMxt - Downwind distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)
0 PtfmCMyt - Lateral distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)
0 PtfmCMzt - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)
0 PtfmRefzt - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform reference point (meters)
---------------------- MASS AND INERTIA ----------------------------------------
0 TipMass(1) - Tip-brake mass, blade 1 (kg)
0 TipMass(2) - Tip-brake mass, blade 2 (kg)
0 TipMass(3) - Tip-brake mass, blade 3 (kg) [unused for 2 blades]
56780 HubMass - Hub mass (kg)
115926 HubIner - Hub inertia about rotor axis [3 blades] or teeter axis [2 blades] (kg m^2)
534.116 GenIner - Generator inertia about HSS (kg m^2)
240000 NacMass - Nacelle mass (kg)
2.60789E+06 NacYIner - Nacelle inertia about yaw axis (kg m^2)
0 YawBrMass - Yaw bearing mass (kg)
0 PtfmMass - Platform mass (kg)
0 PtfmRIner - Platform inertia for roll tilt rotation about the platform CM (kg m^2)
0 PtfmPIner - Platform inertia for pitch tilt rotation about the platform CM (kg m^2)
0 PtfmYIner - Platform inertia for yaw rotation about the platform CM (kg m^2)
---------------------- BLADE ---------------------------------------------------
17 BldNodes - Number of blade nodes (per blade) used for analysis (-)
âNRELOffshrBsline5MW_Blade.datâ BldFile(1) - Name of file containing properties for blade 1 (quoted string)
âNRELOffshrBsline5MW_Blade.datâ BldFile(2) - Name of file containing properties for blade 2 (quoted string)
âNRELOffshrBsline5MW_Blade.datâ BldFile(3) - 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 TeetDmpP - Rotor-teeter damper position (degrees) [used only for 2 blades and when TeetMod=1]
0 TeetDmp - Rotor-teeter damping constant (N-m/(rad/s)) [used only for 2 blades and when TeetMod=1]
0 TeetCDmp - Rotor-teeter rate-independent Coulomb-damping moment (N-m) [used only for 2 blades and when TeetMod=1]
0 TeetSStP - Rotor-teeter soft-stop position (degrees) [used only for 2 blades and when TeetMod=1]
0 TeetHStP - Rotor-teeter hard-stop position (degrees) [used only for 2 blades and when TeetMod=1]
0 TeetSSSp - Rotor-teeter soft-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1]
0 TeetHSSp - Rotor-teeter hard-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1]
---------------------- DRIVETRAIN ----------------------------------------------
100 GBoxEff - Gearbox efficiency (%)
97 GBRatio - Gearbox ratio (-)
8.67637E+08 DTTorSpr - Drivetrain torsional spring (N-m/rad)
6.215E+06 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 ---------------------------------------------------
20 TwrNodes - Number of tower nodes used for analysis (-)
âNRELOffshrBsline5MW_Onshore_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)
0 TStart - Time to begin tabular output (s) (currently unused)
1 DecFact - Decimation factor for tabular output {1: output every time step} (-) (currently unused)
0 NTwGages - Number of tower nodes that have strain gages for output [0 to 9] (-)
10, 19, 28 TwrGagNd - List of tower nodes that have strain gages [1 to TwrNodes] (-) [unused if NTwGages=0]
3 NBlGages - Number of blade nodes that have strain gages for output [0 to 9] (-)
5, 9, 13 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, (-)
âOoPDefl1â - Blade 1 out-of-plane and in-plane deflections and tip twist
âIPDefl1â - Blade 1 out-of-plane and in-plane deflections and tip twist
âTwstDefl1â - Blade 1 out-of-plane and in-plane deflections and tip twist
âBldPitch1â - Blade 1 pitch angle
âHSShftTqâ
âHSShftPwrâ
âYawPznâ
âAzimuthâ - Blade 1 azimuth angle
âRotSpeedâ - Low-speed shaft and high-speed shaft speeds
âGenSpeedâ - Low-speed shaft and high-speed shaft speeds
âTTDspFAâ - Tower fore-aft and side-to-side displacements and top twist
âTTDspSSâ - Tower fore-aft and side-to-side displacements and top twist
âTTDspTwstâ - Tower fore-aft and side-to-side displacements and top twist
âSpn2MLxb1â - Blade 1 local edgewise and flapwise bending moments at span station 2 (approx. 50% span)
âSpn2MLyb1â - Blade 1 local edgewise and flapwise bending moments at span station 2 (approx. 50% span)
âRootFxb1â - Out-of-plane shear, in-plane shear, and axial forces at the root of blade 1
âRootFyb1â - Out-of-plane shear, in-plane shear, and axial forces at the root of blade 1
âRootFzb1â - Out-of-plane shear, in-plane shear, and axial forces at the root of blade 1
âRootMxb1â - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1
âRootMyb1â - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1
âRootMzb1â - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1
âRotTorqâ - Rotor torque and low-speed shaft 0- and 90-bending moments at the main bearing
âLSSGagMyaâ - Rotor torque and low-speed shaft 0- and 90-bending moments at the main bearing
âLSSGagMzaâ - Rotor torque and low-speed shaft 0- and 90-bending moments at the main bearing
âYawBrFxpâ - Fore-aft shear, side-to-side shear, and vertical forces at the top of the tower (not rotating with nacelle yaw)
âYawBrFypâ - Fore-aft shear, side-to-side shear, and vertical forces at the top of the tower (not rotating with nacelle yaw)
âYawBrFzpâ - Fore-aft shear, side-to-side shear, and vertical forces at the top of the tower (not rotating with nacelle yaw)
âYawBrMxpâ - Side-to-side bending, fore-aft bending, and yaw moments at the top of the tower (not rotating with nacelle yaw)
âYawBrMypâ - Side-to-side bending, fore-aft bending, and yaw moments at the top of the tower (not rotating with nacelle yaw)
âYawBrMzpâ - Side-to-side bending, fore-aft bending, and yaw moments at the top of the tower (not rotating with nacelle yaw)
âTwrBsFxtâ - Fore-aft shear, side-to-side shear, and vertical forces at the base of the tower (mudline)
âTwrBsFytâ - Fore-aft shear, side-to-side shear, and vertical forces at the base of the tower (mudline)
âTwrBsFztâ - Fore-aft shear, side-to-side shear, and vertical forces at the base of the tower (mudline)
âTwrBsMxtâ - Side-to-side bending, fore-aft bending, and yaw moments at the base of the tower (mudline)
âTwrBsMytâ - Side-to-side bending, fore-aft bending, and yaw moments at the base of the tower (mudline)
âTwrBsMztâ - Side-to-side bending, fore-aft bending, and yaw moments at the base of the tower (mudline)
END of input file (the word âENDâ must appear in the first 3 columns of this last OutList line)
[/code]
Here is my AeroDyn file as well:
[code]------- AERODYN v15 for OpenFAST INPUT FILE -----------------------------------------------
NREL 5.0 MW offshore baseline aerodynamic input properties, with OC3 Monopile tower
====== General Options ============================================================================
False Echo - Echo the input to â.AD.echâ? (flag)
âdefaultâ DTAero - Time interval for aerodynamic calculations {or âdefaultâ} (s)
1 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? (flag)
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 CompAA - Flag to compute AeroAcoustics calculation [used only when WakeMod = 1 or 2]
âunusedâ AA_InputFile - AeroAcoustics input file [used only when CompAA=true]
====== Environmental Conditions ===================================================================
1.225 AirDens - Air density (kg/m^3)
1.464E-05 KinVisc - Kinematic air viscosity (m^2/s)
335 SpdSound - Speed of sound (m/s)
103500 Patm - Atmospheric pressure ¶ [used only when CavitCheck=True]
1700 Pvap - Vapour pressure of fluid ¶ [used only when CavitCheck=True]
0.5 FluidDepth - Water depth above mid-hub height (m) [used only when CavitCheck=True]
====== Blade-Element/Momentum Theory Options ====================================================== [unused when WakeMod=0 or 3]
2 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/32*pi} (-) [used only when SkewMod=2; unused when WakeMod=0 or 3]
True TipLoss - Use the Prandtl tip-loss model? (flag) [unused when WakeMod=0 or 3]
True HubLoss - Use the Prandtl hub-loss model? (flag) [unused when WakeMod=0 or 3]
True 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]
100 MaxIter - Maximum number of iteration steps (-) [unused when WakeMod=0]
====== Dynamic Blade-Element/Momentum Theory Options ============================================== [used only when WakeMod=2]
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]
âunusedâ OLAFInputFileName - Input file for OLAF [used only when WakeMod=3]
====== Beddoes-Leishman Unsteady Airfoil Aerodynamics Options ===================================== [used only when AFAeroMod=2]
3 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]
====== 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 (-)
8 NumAFfiles - Number of airfoil files used (-)
âAirfoils/Cylinder1.datâ AFNames - Airfoil file names (NumAFfiles lines) (quoted strings)
âAirfoils/Cylinder2.datâ
âAirfoils/DU40_A17.datâ
âAirfoils/DU35_A17.datâ
âAirfoils/DU30_A17.datâ
âAirfoils/DU25_A17.datâ
âAirfoils/DU21_A17.datâ
âAirfoils/NACA64_A17.datâ
====== Rotor/Blade Properties =====================================================================
True UseBlCm - Include aerodynamic pitching moment in calculations? (flag)
âNRELOffshrBsline5MW_AeroDyn_blade.datâ ADBlFile(1) - Name of file containing distributed aerodynamic properties for Blade #1 (-)
âNRELOffshrBsline5MW_AeroDyn_blade.datâ ADBlFile(2) - Name of file containing distributed aerodynamic properties for Blade #2 (-) [unused if NumBl < 2]
âNRELOffshrBsline5MW_AeroDyn_blade.datâ ADBlFile(3) - Name of file containing distributed aerodynamic properties for Blade #3 (-) [unused if NumBl < 3]
====== Tower Influence and Aerodynamics ============================================================= [used only when TwrPotent/=0, TwrShadow=True, or TwrAero=True]
12 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, or TwrAero=True]
TwrElev TwrDiam TwrCd TwrTI (used only with TwrShadow=2)
(m) (m) (-) (-)
0.0000000E+00 6.0000000E+00 1.0000000E+00 1.0000000E-01
8.5261000E+00 5.7870000E+00 1.0000000E+00 1.0000000E-01
1.7053000E+01 5.5740000E+00 1.0000000E+00 1.0000000E-01
2.5579000E+01 5.3610000E+00 1.0000000E+00 1.0000000E-01
3.4105000E+01 5.1480000E+00 1.0000000E+00 1.0000000E-01
4.2633000E+01 4.9350000E+00 1.0000000E+00 1.0000000E-01
5.1158000E+01 4.7220000E+00 1.0000000E+00 1.0000000E-01
5.9685000E+01 4.5090000E+00 1.0000000E+00 1.0000000E-01
6.8211000E+01 4.2960000E+00 1.0000000E+00 1.0000000E-01
7.6738000E+01 4.0830000E+00 1.0000000E+00 1.0000000E-01
8.5268000E+01 3.8700000E+00 1.0000000E+00 1.0000000E-01
8.7600000E+01 3.8700000E+00 1.0000000E+00 1.0000000E-01
====== Outputs ====================================================================================
True SumPrint - Generate a summary file listing input options and interpolated properties to â.AD.sumâ? (flag)
0 NBlOuts - Number of blade node outputs [0 - 9] (-)
2, 9, 16, 4, 5, 7, 11, 13, 14 BlOutNd - Blade nodes whose values will be output (-)
0 NTwOuts - Number of tower node outputs [0 - 9] (-)
1, 2, 3, 4, 5 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, (-)
âRtAeroFxhâ
âRtAeroCtâ
âRtAeroCpâ
âRtSpeedâ
âRtVAvgxhâ
âRtTSRâ
END of input file (the word âENDâ must appear in the first 3 columns of this last OutList line)
[/code]
I agree with your statements and I was expecting the same. Do you think the issue comes from one of the files?
OK, thanks for clarifying. You can find a complete OpenFAST model of the NREL 5-MW baseline wind turbine from the OpenFAST r-test (github.com/OpenFAST/r-test), but I suspect you are using the correct data anyway.
Another thought I have is that the inflow to the second wind turbine is not uniform, but includes a (likely Gaussian-shaped) wake deficit from the upwind turbine. Perhaps the value of the rotor disk-averaged wind speed (RtVAvgxh), which is also used to compute the rotor aerodynamic power coefficient (RtAeroCp), is not computing the disk-averaged wind speed accurately enough for this inflow condition. RtVAvgxh is computed simply by summing the relative velocity at each aerodynamic analysis node along each blade and dividing by the total number of nodes. In reality, the more outboard nodes should have a higher weighting in this summation than the inboard nodes because they sweep a larger area, but this weighting is not included in the spatial averaging. You could look at the inflow to the second turbine from the FAST.Farm outputs and manually compute a more accurate representation of the disk-averaged wind speed and associated rotor aerodynamic power coefficient. I suspect if you do that, youâll get a TSR and Cp that is more consistent with the upstream turbine.
Thank you so much for the reply and guidance. For the maximum power coefficient, I found using the AeroDyn module that it is Cp=0.4862, for a TSR_opt=7.55 and Pitch_opt=0°, as we discussed it in the âNREL 5-MW reference turbine - CP, CQ, CT Coefficientsâ topic [url]Question #1 for turbulence fields (points outside of field)], which is very close to the 0.482 found a while ago. I think the difference comes from the blade nodes considered while computing Cp, and the different aerodynamic settings in AeroDyn v15, and/or enabling/disabling structural DOFs may have some effect, as you highlighted.
Regarding your suggestion for the computation of RtVAvgxh, I compute it now by double integrating the local wind speeds, in the following nodes, on the swept area divided by (R^2)*pi.
Hi Jason,
A basic question as I havenât had much experience with openFAST.
I am working with r-tests. But I am struggling to get around the naming convention of r-tests. ReadMe file doesnât have much information. Also, I didnât get much info from the ânaming scheme from testNNâ. E.g what does â5MW_ITIBarge_DLL_WTurb_WavesIrrâ mean?
I need to run the simulation for a 5MW, three-bladed,variable-speed offshore(fixed monopile) turbine. Which test should I refer to?
Warm regards
The r-test named 5MW_ITIBarge_DLL_WTurb_WavesIrr is a simulation of the NREL 5-MW baseline wind turbine atop the floating offshore ITI Energy barge with the controller DISCON DLL excited by wind turbulence and irregular waves.
