Skip to content

FAST.Farm wake origin moving discontinuously and offset from turbine vtp #2661

Closed
@ReillanS

Description

@ReillanS

Hello, I've been modeling a small turbine in openfast (0.6m rotor diameter) but when I simulate it in FAST.Farm there is an issue with the wake. The wake does not seem to originate from the correct y position and jumps around discontinuously over time. This issue is shown visually in the following images, which are taken using ParaView to visualize an XY plane of disturbed wind data along with the turbine's vtp files:

Image
Image

The vtp files seem to show the turbine in the correct spot from the descriptions of the coordinate system in the documentation.
Any help in understanding this issue would be greatly appreciated. Below is the fstf file used and inflow wind file, I am new to FAST.Farm, so let me know if there is any other information needed.

FAST.Farm version: 3.5.4
Windows 10

fstf file:
------- FAST.Farm for OpenFAST INPUT FILE -------------------------------------------------
FAST.Farm input file, using two turbines, with a turbulent inflow given by one turbsim box
--- SIMULATION CONTROL ---
True Echo - Echo input data to .ech? (flag)
FATAL AbortLevel - Error level when simulation should abort (string) {"WARNING", "SEVERE", "FATAL"}
10.0 TMax - Total run time (s) [>=0.0]
False UseSC - Use a super controller? (flag)
2 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.04 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]
1.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]
0.1 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]
0.3 DT_Low - Time step for low -resolution wind data interpolation; will be used as the global FAST.Farm time step (s) [>0.0]
0.05 DT_High - Time step for high-resolution wind data interpolation (s) [>0.0]
200 NX_Low - Number of low -resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
101 NY_Low - Number of low -resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
51 NZ_Low - Number of low -resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
-10 X0_Low - Origin of low -resolution spatial nodes in X direction for wind data interpolation (m)
-4.95 Y0_Low - Origin of low -resolution spatial nodes in Y direction for wind data interpolation (m)
0.11 Z0_Low - Origin of low -resolution spatial nodes in Z direction for wind data interpolation (m)
0.12 dX_Low - Spacing of low -resolution spatial nodes in X direction for wind data interpolation (m) [>0.0]
0.097 dY_Low - Spacing of low -resolution spatial nodes in Y direction for wind data interpolation (m) [>0.0]
0.098 dZ_Low - Spacing of low -resolution spatial nodes in Z direction for wind data interpolation (m) [>0.0]
26 NX_High - Number of high-resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
41 NY_High - Number of high-resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
40 NZ_High - Number of high-resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
"Inflow_Wind\InflowWind.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 1.5 0.0 "G06.fst" -1 -1 0.11 0.08 0.075 0.12
3.0 1.5 0.0 "G06.fst" 2 -1 0.11 0.08 0.075 0.12
--- WAKE DYNAMICS ---
3 Mod_Wake - Switch between wake formulations {1:Polar, 2:Curl, 3:Cartesian} (-) (switch)
0.05 dr - Radial increment of radial finite-difference grid (m) [>0.0]
55 NumRadii - Number of radii in the radial finite-difference grid (-) [>=2]
3 NumPlanes - Number of wake planes (-) [>=2]
20.8 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 ]
DEFAULT 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]
--- CURLED-WAKE PARAMETERS [only used if Mod_Wake=2 or 3] ---
DEFAULT Swirl - Switch to include swirl velocities in wake (-) (switch) [DEFAULT=TRUE]
DEFAULT k_VortexDecay - Vortex decay constant for curl (-) [DEFAULT=0]
DEFAULT NumVortices - The number of vortices in the curled wake model (-) [DEFAULT=100]
DEFAULT sigma_D - The width of the vortices in the curled wake model non-dimesionalized by rotor diameter (-) [DEFAULT=0.2]
DEFAULT FilterInit - Switch to filter the initial wake plane deficit and select the number of grid points for the filter {0: no filter, 1: filter of size 1} or DEFAULT [DEFAULT=1] (switch)
DEFAULT k_vCurl - Calibrated parameter for scaling the eddy viscosity in the curled-wake model (-) [>=0] or DEFAULT [DEFAULT=2.0 ]
DEFAULT Mod_Projection - Switch to select how the wake plane velocity is projected in AWAE {1: keep all components, 2: project against plane normal} or DEFAULT [DEFAULT=1: if Mod_Wake is 1 or 3, or DEFAULT=2: if Mod_Wake is 2] (switch)
--- VISUALIZATION ---
True WrDisWind - Write low- and high-resolution disturbed wind data to .Low.Dis.t.vtk etc.? (flag)
4 NOutDisWindXY - Number of XY planes for output of disturbed wind data across the low-resolution domain to .Low.DisXY<n_out>.t.vtk (-) [0 to 99]
1.5, 1.9, 2.3, 2.5 OutDisWindZ - Z coordinates of XY planes for output of disturbed wind data across the low-resolution domain (m) [1 to NOutDisWindXY] [unused for NOutDisWindXY=0]
0 NOutDisWindYZ - Number of YZ planes for output of disturbed wind data across the low-resolution domain to /Low.DisYZ<n_out>.t.vtk (-) [0 to 99]
0 OutDisWindX - X coordinates of YZ planes for output of disturbed wind data across the low-resolution domain (m) [1 to NOutDisWindYZ] [unused for NOutDisWindYZ=0]
0 NOutDisWindXZ - Number of XZ planes for output of disturbed wind data across the low-resolution domain to /Low.DisXZ<n_out>.t.vtk (-) [0 to 99]
0.0 OutDisWindY - Y coordinates of XZ planes for output of disturbed wind data across the low-resolution domain (m) [1 to NOutDisWindXZ] [unused for NOutDisWindXZ=0]
0.3 WrDisDT - Time step for disturbed wind visualization output (s) [>0.0] or DEFAULT [DEFAULT=DT_Low or DT_Low-VTK] [unused for WrDisWind=False and NOutDisWindXY=NOutDisWindYZ=NOutDisWindXZ=0]
--- OUTPUT ---
False SumPrint - Print summary data to .sum? (flag)
99999.9 ChkptTime - Amount of time between creating checkpoint files for potential restart (s) [>0.0]
0.0 TStart - Time to begin tabular output (s) [>=0.0]
1 OutFileFmt - Format for tabular (time-marching) output file (switch) {1: text file [.out], 2: binary file [.outb], 3: both}
True TabDelim - Use tab delimiters in text tabular output file? (flag) {uses spaces if False}
"ES10.3E2" OutFmt - Format used for text tabular output, excluding the time channel. Resulting field should be 10 characters. (quoted string)
DEFAULT OutAllPlanes - Output all wake planes at all time steps. [DEFAULT=False]
3 NOutRadii - Number of radial nodes for wake output for an individual rotor (-) [0 to 20]
0, 2, 5 OutRadii - List of radial nodes for wake output for an individual rotor (-) [1 to NOutRadii] [unused for NOutRadii=0]
1 NOutDist - Number of downstream distances for wake output for an individual rotor (-) [0 to 9 ]
5.0 OutDist - List of downstream distances for wake output for an individual rotor (m) [1 to NOutDist ] [unused for NOutDist =0]
1 NWindVel - Number of points for wind output (-) [0 to 9]
0.0 WindVelX - List of coordinates in the X direction for wind output (m) [1 to NWindVel] [unused for NWindVel=0]
1.5 WindVelY - List of coordinates in the Y direction for wind output (m) [1 to NWindVel] [unused for NWindVel=0]
2.0 WindVelZ - List of coordinates in the Z direction for wind output (m) [1 to NWindVel] [unused for NWindVel=0]
OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels (quoted string)
"RtVAmbT1"
"RtVAmbT2"
"RtVRelT2"
"YawErrT1"
"YawErrT2"
"CtT1N02 , CtT1N03 , CtT1N04 , CtT1N05 , CtT1N06"
WkAxsXT2D1"
WkAxsYT2D1"
WkAxsZT2D1"
WkPosXT2D1"
WkPosYT2D1"
WkPosZT2D1"
"WkDfVxT1N01D1, WkDfVxT1N02D1, WkDfVxT1N03D1, WkDfVxT1N04D1, WkDfVxT1N05D1, WkDfVxT1N06D1, WkDfVxT1N07D1"
"WkDfVrT1N01D1, WkDfVrT1N02D1, WkDfVrT1N03D1, WkDfVrT1N04D1, WkDfVrT1N05D1, WkDfVrT1N06D1, WkDfVrT1N07D1"
"WkDfVxT2N01D1, WkDfVxT2N02D1, WkDfVxT2N03D1, WkDfVxT2N04D1, WkDfVxT2N05D1, WkDfVxT2N06D1, WkDfVxT2N07D1"
"WkDfVrT2N01D1, WkDfVrT2N02D1, WkDfVrT2N03D1, WkDfVrT2N04D1, WkDfVrT2N05D1, WkDfVrT2N06D1, WkDfVrT2N07D1"
END of input file (the word "END" must appear in the first 3 columns of this last OutList line)

inflow wind file:
------- InflowWind v3.01.* INPUT FILE -------------------------------------------------------------------------
Generated with AeroElasticSE FAST driver

False Echo - Echo input data to .ech (flag)
3 WindType - switch for wind file type (1=steady; 2=uniform; 3=binary TurbSim FF; 4=binary Bladed-style FF; 5=HAWC format; 6=User defined)
0.0 PropagationDir - Direction of wind propagation (meteoroligical rotation from aligned with X (positive rotates towards -Y) -- degrees)
0.0 VFlowAng - Upflow angle (degrees) (not used for native Bladed format WindType=7)
False VelInterpCubic - Use cubic interpolation for velocity in time (false=linear, true=cubic) [Used with WindType=2,3,4,5,7]
1 NWindVel - Number of points to output the wind velocity (0 to 9)
0.0 0 0 0 0 0 0 0 0 WindVxiList - List of coordinates in the inertial X direction (m)
0 -90 -60 -30 0 30 60 90 130 WindVyiList - List of coordinates in the inertial Y direction (m)
1 110 110 110 110 110 110 110 110 WindVziList - List of coordinates in the inertial Z direction
================== Parameters for Steady Wind Conditions [used only for WindType = 1] =========================
8.5 HWindSpeed - Horizontal windspeed (m/s)
0.64 RefHt - Reference height for horizontal wind speed (m)
0.2 PLexp - Power law exponent (-)
================== Parameters for Uniform wind file [used only for WindType = 2] ============================
"Wind/IEAonshore_10mps.inp" Filename_Uni - Filename of time series data for uniform wind field. (-)
110.0 RefHt_Uni - Reference height for horizontal wind speed (m)
125.88 RefLength - Reference length for linear horizontal and vertical sheer (-)
================== Parameters for Binary TurbSim Full-Field files [used only for WindType = 3] ==============
"TurbSim_U7_T5_5x10.bts" Filename_BTS - Name of the Full field wind file to use (.bts)
================== Parameters for Binary Bladed-style Full-Field files [used only for WindType = 4] =========
"TurbSim" FilenameRoot - Rootname of the full-field wind file to use (.wnd, .sum)
False TowerFile - Have tower file (.twr) (flag)
================== Parameters for HAWC-format binary files [Only used with WindType = 5] =====================
"./Wind/basic_5u.bin" FileName_u - name of the file containing the u-component fluctuating wind (.bin)
"./Wind/basic_5v.bin" FileName_v - name of the file containing the v-component fluctuating wind (.bin)
"./Wind/basic_5w.bin" FileName_w - name of the file containing the w-component fluctuating wind (.bin)
64 nx - number of grids in the x direction (in the 3 files above) (-)
32 ny - number of grids in the y direction (in the 3 files above) (-)
32 nz - number of grids in the z direction (in the 3 files above) (-)
16.0 dx - distance (in meters) between points in the x direction (m)
3.0 dy - distance (in meters) between points in the y direction (m)
3.0 dz - distance (in meters) between points in the z direction (m)
110.0 RefHt_Hawc - reference height; the height (in meters) of the vertical center of the grid (m)
------------- Scaling parameters for turbulence ---------------------------------------------------------
2 ScaleMethod - Turbulence scaling method [0 = none, 1 = direct scaling, 2 = calculate scaling factor based on a desired standard deviation]
1.0 SFx - Turbulence scaling factor for the x direction (-) [ScaleMethod=1]
1.0 SFy - Turbulence scaling factor for the y direction (-) [ScaleMethod=1]
1.0 SFz - Turbulence scaling factor for the z direction (-) [ScaleMethod=1]
1.2 SigmaFx - Turbulence standard deviation to calculate scaling from in x direction (m/s) [ScaleMethod=2]
0.8 SigmaFy - Turbulence standard deviation to calculate scaling from in y direction (m/s) [ScaleMethod=2]
0.2 SigmaFz - Turbulence standard deviation to calculate scaling from in z direction (m/s) [ScaleMethod=2]
------------- Mean wind profile parameters (added to HAWC-format files) ---------------------------------
8.0 URef - Mean u-component wind speed at the reference height (m/s)
2 WindProfile - Wind profile type (0=constant;1=logarithmic,2=power law)
0.2 PLExp_Hawc - Power law exponent (-) (used for PL wind profile type only)
0.03 Z0 - Surface roughness length (m) (used for LG wind profile type only)
0 XOffset - Initial offset in +x direction (shift of wind box) (-)
================== LIDAR Parameters ===========================================================================
0 SensorType - Switch for lidar configuration (0 = None, 1 = Single Point Beam(s), 2 = Continuous, 3 = Pulsed)
0 NumPulseGate - Number of lidar measurement gates (used when SensorType = 3)
30 PulseSpacing - Distance between range gates (m) (used when SensorType = 3)
0 NumBeam - Number of lidar measurement beams (0-5)(used when SensorType = 1)
-200 FocalDistanceX - Focal distance co-ordinates of the lidar beam in the x direction (relative to hub height) (only first coordinate used for SensorType 2 and 3) (m)
0 FocalDistanceY - Focal distance co-ordinates of the lidar beam in the y direction (relative to hub height) (only first coordinate used for SensorType 2 and 3) (m)
0 FocalDistanceZ - Focal distance co-ordinates of the lidar beam in the z direction (relative to hub height) (only first coordinate used for SensorType 2 and 3) (m)
0.0 0.0 0.0 RotorApexOffsetPos - Offset of the lidar from hub height (m)
17 URefLid - Reference average wind speed for the lidar[m/s]
0.25 MeasurementInterval - Time between each measurement [s]
False LidRadialVel - TRUE => return radial component, FALSE => return 'x' direction estimate
1 ConsiderHubMotion - Flag whether to consider the hub motion's impact on Lidar measurements
====================== OUTPUT ==================================================
True SumPrint - Print summary data to .IfW.sum (flag)
OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-)
"WindHubVelXY" - horizontal component of wind at (moving) hub point (m/s)
END of input file (the word "END" must appear in the first 3 columns of this last OutList line)

Metadata

Metadata

Assignees

No one assigned

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions