Unable to run FAST.Farm simulation

[j-killeen]

Dear @jjonkman

I hope this message finds you well. I am trying to run a wind farm simulation that contains 2 5MW_OC3Spar_DLL_WTurb_WavesIrr offshore fst files. When I run the simulation, I am constantly getting the following error message:

"forrtl: severe (47): write to READONLY file, unit 13, file C:\Users\jamie\Downloads\r-test-548bf2143633a7066f4582167fdf4180951fd7f6 (1)\r-test-548bf2143633a7066f4582167fdf4180951fd7f6\glue-codes\openfast\5MW_OC3Spar_DLL_WTurb_WavesIrr\fort.13
Image PC Routine Line Source
FAST.Farm_x64_OMP 00007FF724AF688F Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF724AB7993 Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF7236A4A3C Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF7236A2AC6 Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF7236A2A5B Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF723698219 Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF72369672C Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF7231052AB Unknown Unknown Unknown
FAST.Farm_x64_OMP 00007FF723263ECA Unknown Unknown Unknown
libiomp5md.dll 00007FFCCEC2A853 Unknown Unknown Unknown
libiomp5md.dll 00007FFCCEB7A617 Unknown Unknown Unknown
libiomp5md.dll 00007FFCCEB7A0DC Unknown Unknown Unknown
libiomp5md.dll 00007FFCCEC025C7 Unknown Unknown Unknown
KERNEL32.DLL 00007FFD275353E0 Unknown Unknown Unknown
ntdll.dll 00007FFD27EE485B Unknown Unknown Unknown"

I don't know if this is a bug issue as I'm using FAST.Farm version 3.5.1. The simulation runs normally without any errors if I remove the second wind turbine, making it only one turbine in the simulation. Any advice or guidance on how to fix this issue would be appreciated.

Here is what I have in my .fstf file in case you need to take a look at that:

"------- FAST.Farm for OpenFAST INPUT FILE -------------------------------------------------
FAST.Farm input file, using two turbines separated by 3D downstream and an offset of 30m, 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"}
90.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]
3.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]
3.0 DT_Low - Time step for low -resolution wind data interpolation; will be used as the global FAST.Farm time step (s) [>0.0]
0.1 DT_High - Time step for high-resolution wind data interpolation (s) [>0.0]
213 NX_Low - Number of low -resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
81 NY_Low - Number of low -resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
30 NZ_Low - Number of low -resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
-100 X0_Low - Origin of low -resolution spatial nodes in X direction for wind data interpolation (m)
-400.0 Y0_Low - Origin of low -resolution spatial nodes in Y direction for wind data interpolation (m)
5.0 Z0_Low - Origin of low -resolution spatial nodes in Z direction for wind data interpolation (m)
10.17 dX_Low - Spacing of low -resolution spatial nodes in X direction for wind data interpolation (m) [>0.0]
10.0 dY_Low - Spacing of low -resolution spatial nodes in Y direction for wind data interpolation (m) [>0.0]
10.0 dZ_Low - Spacing of low -resolution spatial nodes in Z direction for wind data interpolation (m) [>0.0]
16 NX_High - Number of high-resolution spatial nodes in X direction for wind data interpolation (-) [>=2]
16 NY_High - Number of high-resolution spatial nodes in Y direction for wind data interpolation (-) [>=2]
16 NZ_High - Number of high-resolution spatial nodes in Z direction for wind data interpolation (-) [>=2]
"IW.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)
50.0 100.0 0.0 "5MW_OC3Spar_DLL_WTurb_WavesIrr.fst" -30 25.0 5.0 10.17 10.0 10.0
50.0 200.0 0.0 "5MW_OC3Spar_DLL_WTurb_WavesIrr.fst" -30 125.0 5.0 10.17 10.0 10.0
--- WAKE DYNAMICS ---
1 Mod_Wake - Switch between wake formulations {1:Polar, 2:Curl, 3:Cartesian} (-) (switch)
5.0 dr - Radial increment of radial finite-difference grid (m) [>0.0]
40 NumRadii - Number of radii in the radial finite-difference grid (-) [>=2]
136 NumPlanes - Number of wake planes (-) [>=2]
0.17 f_c - Cutoff (corner) frequency of the low-pass time-filter for the wake advection, deflection, and meandering model [recommended=1.28*U0/R] (Hz) [>0.0] or DEFAULT [DEFAULT=12.5/R, R estimated from dr and NumRadii, not recommended]
DEFAULT C_HWkDfl_O - Calibrated parameter in the correction for wake deflection defining the horizontal offset at the rotor (m ) or DEFAULT [DEFAULT= 0.0 ]
DEFAULT C_HWkDfl_OY - Calibrated parameter in the correction for wake deflection defining the horizontal offset at the rotor scaled with yaw error (m/deg) or DEFAULT [DEFAULT= 0.0 if Mod_Wake is 2, 0.3 otherwise]
DEFAULT C_HWkDfl_x - Calibrated parameter in the correction for wake deflection defining the horizontal offset scaled with downstream distance (- ) or DEFAULT [DEFAULT= 0.0 ]
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 ---
False WrDisWind - Write low- and high-resolution disturbed wind data to .Low.Dis.t.vtk etc.? (flag)
0 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]
90.0 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]
63.,441. 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]
3.0 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 ---
True 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]
7 NOutRadii - Number of radial nodes for wake output for an individual rotor (-) [0 to 20]
0, 2, 5, 11, 17, 21, 39 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 ]
63.0 OutDist - List of downstream distances for wake output for an individual rotor (m) [1 to NOutDist ] [unused for NOutDist =0]
0 NWindVel - Number of points for wind output (-) [0 to 9]
95.0 WindVelX - List of coordinates in the X direction for wind output (m) [1 to NWindVel] [unused for NWindVel=0]
200.0 WindVelY - List of coordinates in the Y direction for wind output (m) [1 to NWindVel] [unused for NWindVel=0]
90.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)"

I've tried a few different X and Y co-ordinates for the second turbine but it will either give me the "outside grid bounds" error or the other error I mentioned above.

Thank you for your help.

[jjonkman]

Dear @j-killeen,

I'm not exactly sure what the issue is, but I noticed that you are using the same OpenFAST primary input file for both wind turbines, which goes against the modeling guidance because presumably you need some differences in the input files, e.g., because of the need to have two distinct controller DLLs. I suggest fixing this issue and see if the simulation runs successfully.

FYI: I see that you have only a 100-m spacing between the two turbines, which is less than 1D for the NREL 5-MW baseline wind turbine. I'm not sure why you've placed the two turbines so close together, but the wakes interactions of FAST.Farm are not accurate if one turbine operates in the near wake of the other turbine. FAST.Farm will be much more accurate if the second turbine is in the far wake (> 3-4D) of the other turbine.

Best regards,

[andrew-platt]

Hi @j-killeen,

This might be a file opening issue that we recently fixed. Could you try version 3.5.4 and let us know if the issue persists there? The input files for 3.5.4 are identical to 3.5.1 (all input files in the 3.5.x series are the same).

Regards

[andrew-platt]

Note: I initially suspected this was fixed in 3.5.4 with #2340. However the error above says write not read. Perhaps there is another location that needs a !$OMP critical around some file opening for writing similar to what is done in #2340.

[deslaughter]

@andrew-platt I think you're right that this is another issue with using GetNewUnit in a parallel context. In this case the program gets the same unit number for a file that's been opened for reading and writing in two separate threads. So it fails when it tries to do the write on a file that's only been opened for reading. Wrapping GetNewUnit and the call to open a file for reading in !$OMP critical should fix the issue.

[andrew-platt]

After attempting to fix this issue for our 3.5.5 release with PR #2538, we decided against fixing this until the 4.0.0 release due to the complexity of the code fix. We expect this issue will be fully resolved in 4.0.0