Query about TurbSim error

Hi @Jason.Jonkman,
I am getting an error in a TurbSim model that I am trying to simulate.
I am using a user-defined wnd shear profile. I would like to create a wind-field for the 5MW DeepCwind turbine. The hub-height is 90 m and so I have set GridHeight to be 175 m and GirdWidth to be 300 m. I have set NumGrid_Z and NumGirid_Y to be 15 x 15. I use the general coherence model with CohExp of 0.85. I have shown the TurbSim file and the Shear Profile below.
I am getting the following error:
Reading the input file “90m_wind_shearJCS1-12.inp”.

  • Calculating the spectral and transfer function matrices:*
    forrtl: severe (157): Program Exception - access violation
    Image PC Routine Line Source
    ntdll.dll 00007FFD15895BB6 Unknown Unknown Unknown
    ntdll.dll 00007FFD158947B1 Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FFE580C0 Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FFE580FA Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FF2D2679 Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FF2D1BBE Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FF2D0F32 Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FF92D2E6 Unknown Unknown Unknown
    TurbSim_x64.exe 00007FF6FFE57E9C Unknown Unknown Unknown
    KERNEL32.DLL 00007FFD15257614 Unknown Unknown Unknown
    ntdll.dll 00007FFD158C26A1 Unknown Unknown Unknown

Can you help me to fix this please ?
Thank you.

Regards,
AOAW

TurbSim File:
---------TurbSim v2.00.* Input File------------------------

---------Runtime Options-----------------------------------
False Echo - Echo input data to .ech (flag)
1 RandSeed1 - First random seed (-2147483648 to 2147483647)
“RANLUX” RandSeed2 - Second random seed (-2147483648 to 2147483647) for intrinsic pRNG, or an alternative pRNG: “RanLux” or “RNSNLW”
False WrBHHTP - Output hub-height turbulence parameters in binary form? (Generates RootName.bin)
True WrFHHTP - Output hub-height turbulence parameters in formatted form? (Generates RootName.dat)
False WrADHH - Output hub-height time-series data in AeroDyn form? (Generates RootName.hh)
True WrADFF - Output full-field time-series data in TurbSim/AeroDyn form? (Generates RootName.bts)
False WrBLFF - Output full-field time-series data in BLADED/AeroDyn form? (Generates RootName.wnd)
False WrADTWR - Output tower time-series data? (Generates RootName.twr)
False WrFMTFF - Output full-field time-series data in formatted (readable) form? (Generates RootName.u, RootName.v, RootName.w)
False WrACT - Output coherent turbulence time steps in AeroDyn form? (Generates RootName.cts)
True Clockwise - Clockwise rotation looking downwind? (used only for full-field binary files - not necessary for AeroDyn)
0 ScaleIEC - Scale IEC turbulence models to exact target standard deviation? [0=no additional scaling; 1=use hub scale uniformly; 2=use individual scales]

--------Turbine/Model Specifications-----------------------
15 NumGrid_Z - Vertical grid-point matrix dimension
15 NumGrid_Y - Horizontal grid-point matrix dimension
0.05 TimeStep - Time step [seconds]
3600 AnalysisTime - Length of analysis time series [seconds] (program will add time if necessary: AnalysisTime = MAX(AnalysisTime, UsableTime+GridWidth/MeanHHWS) )
3600 UsableTime - Usable length of output time series [seconds] (program will add GridWidth/MeanHHWS seconds unless UsableTime is “ALL”)
90.0 HubHt - Hub height [m] (should be > 0.5GridHeight)
175 GridHeight - Grid height [m]
300 GridWidth - Grid width [m] (should be >= 2
(RotorRadius+ShaftLength))
0 VFlowAng - Vertical mean flow (uptilt) angle [degrees]
0 HFlowAng - Horizontal mean flow (skew) angle [degrees]

--------Meteorological Boundary Conditions-------------------
“IECKAI” TurbModel - Turbulence model (“IECKAI”,“IECVKM”,“GP_LLJ”,“NWTCUP”,“SMOOTH”,“WF_UPW”,“WF_07D”,“WF_14D”,“TIDAL”,“API”,“USRINP”,“TIMESR”, or “NONE”)
“unused” UserFile - Name of the file that contains inputs for user-defined spectra or time series inputs (used only for “USRINP” and “TIMESR” models)
“1-ED3” IECstandard - Number of IEC 61400-x standard (x=1,2, or 3 with optional 61400-1 edition number (i.e. “1-Ed2”) )
7.38 IECturbc - IEC turbulence characteristic (“A”, “B”, “C” or the turbulence intensity in percent) (“KHTEST” option with NWTCUP model, not used for other models)
“NTM” IEC_WindType - IEC turbulence type (“NTM”=normal, “xETM”=extreme turbulence, “xEWM1”=extreme 1-year wind, “xEWM50”=extreme 50-year wind, where x=wind turbine class 1, 2, or 3)
“default” ETMc - IEC Extreme Turbulence Model “c” parameter [m/s]
“USR” WindProfileType - Velocity profile type (“LOG”;“PL”=power law;“JET”;“H2L”=Log law for TIDAL model;“API”;“USR”;“TS”;“IEC”=PL on rotor disk, LOG elsewhere; or “default”)
“windshear12km.dat” ProfileFile - Name of the file that contains input profiles for WindProfileType=“USR” and/or TurbModel=“USRVKM” [-]
90.0 RefHt - Height of the reference velocity (URef) [m]
“default” URef - Mean (total) velocity at the reference height [m/s] (or “default” for JET velocity profile) [must be 1-hr mean for API model; otherwise is the mean over AnalysisTime seconds]
“default” ZJetMax - Jet height [m] (used only for JET velocity profile, valid 70-490 m)
“default” PLExp - Power law exponent [-] (or “default”)
“default” Z0 - Surface roughness length [m] (or “default”)

--------Non-IEC Meteorological Boundary Conditions------------
“default” Latitude - Site latitude [degrees] (or “default”)
0.05 RICH_NO - Gradient Richardson number [-]
“default” UStar - Friction or shear velocity [m/s] (or “default”)
“default” ZI - Mixing layer depth [m] (or “default”)
“default” PC_UW - Hub mean u’w’ Reynolds stress [m^2/s^2] (or “default” or “none”)
“default” PC_UV - Hub mean u’v’ Reynolds stress [m^2/s^2] (or “default” or “none”)
“default” PC_VW - Hub mean v’w’ Reynolds stress [m^2/s^2] (or “default” or “none”)

--------Spatial Coherence Parameters----------------------------
“General” SCMod1 - u-component coherence model (“GENERAL”, “IEC”, “API”, “NONE”, or “default”)
“General” SCMod2 - v-component coherence model (“GENERAL”, “IEC”, “NONE”, or “default”)
“General” SCMod3 - w-component coherence model (“GENERAL”, “IEC”, “NONE”, or “default”)
“default” InCDec1 - u-component coherence parameters for general or IEC models [-, m^-1] (e.g. “10.0 0.3e-3” in quotes) (or “default”)
“default” InCDec2 - v-component coherence parameters for general or IEC models [-, m^-1] (e.g. “10.0 0.3e-3” in quotes) (or “default”)
“default” InCDec3 - w-component coherence parameters for general or IEC models [-, m^-1] (e.g. “10.0 0.3e-3” in quotes) (or “default”)
0.85 CohExp - Coherence exponent for general model [-] (or “default”)

--------Coherent Turbulence Scaling Parameters-------------------
“/home/lkilcher/work/nwtc/turbsim/trunk/Test/EventData/” CTEventPath - Name of the path where event data files are located
“Random” CTEventFile - Type of event files (“LES”, “DNS”, or “RANDOM”)
true Randomize - Randomize the disturbance scale and locations? (true/false)
1 DistScl - Disturbance scale [-] (ratio of event dataset height to rotor disk). (Ignored when Randomize = true.)
0.5 CTLy - Fractional location of tower centerline from right [-] (looking downwind) to left side of the dataset. (Ignored when Randomize = true.)
0.5 CTLz - Fractional location of hub height from the bottom of the dataset. [-] (Ignored when Randomize = true.)
30 CTStartTime - Minimum start time for coherent structures in RootName.cts [seconds]

====================================================
! NOTE: Do not add or remove any lines in this file!

Wind shear profile:

---------TurbSim v2.00.* Profile Input File------------------------
Example file using completely made up profiles
-------- User-Defined Profiles (Used only with USR wind profile or USRVKM spectral model) -------------
38
1.092 StdScale1 - u-component scaling factor for the input standard deviat ion
1.0 StdScale2 - v-component scaling factor for the input standard deviat ion
0.534 StdScale3 - w-component scaling factor for the input standard deviat ion
------------------------------------------------------------------------- --------- -
Height Wind Speed Wind Direction Standard Deviation Length Sc ale
(m) (m/s) (deg, cntr-clockwise ) (m/s) (m)
------------------------------------------------------------------------- --------- -
9.477484385 68.96623387 0
14.62264151 72.40591398 0
17.02928068 73.67001496 0
20.53387105 75.51082109 0
28.1133688 78.50395199 0
31.76020612 79.94408192 0
38.34483475 81.52138174 0
42.86451582 82.60403841 0
45.40941838 83.00115007 0
50.55211484 83.80362664 0
55.9819332 84.4805294 0
62.26655144 85.26399491 0
65.40949825 85.53955518 0
71.78453118 86.09849107 0
74.37421127 86.30573511 0
82.0754717 86.92204301 0
86.90775706 87.15757248 0
98.485148 87.72186382 0
100 87.79569892 0
110.2817846 88.17935743 0
119.8113208 88.53494624 0
128.9663903 88.71554873 0
144.6028718 89.02401033 0
150.4716981 89.13978495 0
153.9116787 89.16031547 0
161.3660325 89.20480463 0
168.8203863 89.24929378 0
174.8861839 89.28549573 0
178.6620821 89.30803109 0
182.7056028 89.33216367 0
184.7522404 89.34074077 0
188.2845322 89.35554399 0
192.7181674 89.37412459 0
196.5671475 89.39025499 0
201.1225859 89.40934604 0
206.128696 89.43032578 0
209.9654958 89.44640514 0
216.9811321 89.47580645 0

Could you try this with a newer version of TurbSim? This case works for me when using the current executable from the main branch of OpenFAST.

Just note that the options in the first section are different:

---------Runtime Options-----------------------------------
False  Echo - Echo input data to .ech (flag)
1 RandSeed1 - First random seed (-2147483648 to 2147483647)
"RANLUX" RandSeed2 - Second random seed (-2147483648 to 2147483647) for intrinsic pRNG, or an alternative pRNG: "RanLux" or "RNSNLW"
False WrBHHTP - Output hub-height turbulence parameters in binary form? (Generates RootName.bin)
True WrFHHTP - Output hub-height turbulence parameters in formatted form? (Generates RootName.dat)
False WrADHH - Output hub-height time-series data in AeroDyn form? (Generates RootName.hh)
True WrADFF - Output full-field time-series data in TurbSim/AeroDyn form? (Generates RootName.bts)
False WrBLFF - Output full-field time-series data in BLADED/AeroDyn form? (Generates RootName.wnd)
False WrADTWR - Output tower time-series data? (Generates RootName.twr)
False         WrHAWCFF        - [Envision addition] Output full-field time-series data in HAWC form?  (Generates RootName-u.bin, RootName-v.bin, RootName-w.bin, RootName.hawc)
False WrFMTFF - Output full-field time-series data in formatted (readable) form? (Generates RootName.u, RootName.v, RootName.w)
False WrACT - Output coherent turbulence time steps in AeroDyn form? (Generates RootName.cts)
0 ScaleIEC - Scale IEC turbulence models to exact target standard deviation? [0=no additional scaling; 1=use hub scale uniformly; 2=use individual scales]

I did have to modify the coherence exponent because it was giving an error about the coherence not being positive definite. If you need to use the exponent of 0.85, you should check the other coherence parameters to make sure they make sense.

default 0.85 CohExp - Coherence exponent for general model [-] (or "default")

Thanks @Bonnie.Jonkman.
Do you have the link for latest version of TurbSim?
Trying to find but no successful yet.

Regards,
AOAW

Hi @Bonnie.Jonkman,
I see TurbSim v2.0 on NREL’s NWTC website but nothing on GitHub. That version from NREL is the one I had used, so if there is a version that links to openfast v3.4.1 can you please direct me to the input files.

Thank you.

Regards,
AOAW

Here is an executable file: https://github.com/OpenFAST/openfast/releases/download/v3.4.1/TurbSim_x64.exe

I don’t know if NREL has been keeping up with the TurbSim input file documentation, but here are example input files: openfast/docs/source/user/turbsim/examples at main · OpenFAST/openfast · GitHub

@Bonnie.Jonkman thank you.

Regards,
AOAW