FF: Curl/Cartesian wake fixes and AWAE VTK wake-plane output

.gitignore

@@ -64,3 +64,8 @@ varcache
 # Python cache files
 openfast_io/dist/
 openfast_io/openfast_io/_version.py
+
+# docker with VScode config
+*.code-workspace
+docker-compose.yml
+.devcontainer/

glue-codes/fast-farm/src/FAST_Farm_Subs.f90

@@ -211,7 +211,11 @@ SUBROUTINE Farm_Initialize( farm, InputFile, ErrStat, ErrMsg )
    call AllocAry( farm%p%MaxNumPlanes, farm%p%NumTurbines, 'farm%p%MaxNumPlanes', ErrStat2, ErrMsg2);  CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName); if (Failed()) return
    do i=1,farm%p%NumTurbines
       ! Eventually, we will have different settings for different rotors
-      farm%p%MaxNumPlanes(i) = ceiling( 15.0 * ( WD_InitInput%InputFileData%NumDFull + WD_InitInput%InputFileData%NumDBuff ) / AWAE_InitInput%InputFileData%C_Meander )
+      if (WD_InitInput%InputFileData%Mod_Wake == Mod_Wake_Polar) then
+         farm%p%MaxNumPlanes(i) = ceiling( 18.0 * ( WD_InitInput%InputFileData%NumDFull + WD_InitInput%InputFileData%NumDBuff ) / AWAE_InitInput%InputFileData%C_Meander )
+      else
+         farm%p%MaxNumPlanes(i) = ceiling( 54.0 * ( WD_InitInput%InputFileData%NumDFull + WD_InitInput%InputFileData%NumDBuff ) )
+      endif
       farm%p%MaxNumPlanes(i) = max( 2, min( farm%p%MaxNumPlanes(i) , farm%p%n_TMax + 2 ) )
    end do
 
@@ -1818,6 +1822,7 @@ SUBROUTINE Transfer_AWAE_to_WD(farm)
       farm%WD(nt)%u%V_plane(:,0:MaxPln) = farm%AWAE%y%V_plane(:,0:MaxPln,nt)  ! Advection, deflection, and meandering velocity of wake planes, m/s
       farm%WD(nt)%u%Vx_wind_disk        = farm%AWAE%y%Vx_wind_disk(nt)        ! Rotor-disk-averaged ambient wind speed, normal to planes, m/s
       farm%WD(nt)%u%TI_amb              = farm%AWAE%y%TI_amb(nt)              ! Ambient turbulence intensity of wind at rotor disk
+      farm%WD(nt)%u%wakePlaneDomainExit(:,0:MaxPln) = farm%AWAE%y%wakePlaneDomainExit(:,0:MaxPln,nt)  ! Per-dimension domain exit flags
    END DO
 
 END SUBROUTINE Transfer_AWAE_to_WD

modules/awae/CMakeLists.txt

@@ -26,6 +26,7 @@ add_library(awaelib STATIC
   src/AWAE.f90
   src/AWAE_IO.f90
   src/AWAE_Types.f90
+  src/AWAE_vtk.f90
   src/amrex_utils.F90
 )
 target_link_libraries(awaelib awaelib_c ifwlib nwtclibs)

modules/awae/src/AWAE_IO.f90

@@ -80,7 +80,7 @@ subroutine WriteDisWindFiles( n, WrDisSkp1, p, y, m, errStat, errMsg )
    ! TimeStamp
    write(Tstr, '(i' // trim(Num2LStr(p%VTK_tWidth)) //'.'// trim(Num2LStr(p%VTK_tWidth)) // ')') n_out ! TODO use n instead..
 
-   FileName = trim(p%OutFileVTKRoot)//".Low.Dis."//trim(Tstr)//".vtk"
+   FileName = trim(p%OutFileFFvtkRoot)//".Low.Dis."//trim(Tstr)//".vtk"
    call WrVTK_SP_header( FileName, "Low resolution disturbed wind for time = "//trim(num2lstr(t_out))//" seconds.", Un, errStat2, errMsg2 )
       call SetErrStat(errStat2, errMsg2, ErrStat, ErrMsg, RoutineName)
       if (ErrStat >= AbortErrLev) return
@@ -92,7 +92,7 @@ subroutine WriteDisWindFiles( n, WrDisSkp1, p, y, m, errStat, errMsg )
       ! We are only writing out the first of the high res data for a given low res time step
       ! NOTE: y%Vdist_high(nt)%data(:,:,:,:,1) is at T=t_low, and index 0 is at T=t_low-DT_high
 
-      FileName = trim(p%OutFileVTKRoot)//".HighT"//trim(num2lstr(nt))//".Dis."//trim(Tstr)//".vtk"
+      FileName = trim(p%OutFileFFvtkRoot)//".HighT"//trim(num2lstr(nt))//".Dis."//trim(Tstr)//".vtk"
       call WrVTK_SP_header( FileName, "High resolution disturbed wind for time = "//trim(num2lstr(t_out))//" seconds.", Un, errStat2, errMsg2 )
          call SetErrStat(errStat2, errMsg2, ErrStat, ErrMsg, RoutineName)
          if (ErrStat >= AbortErrLev) return

modules/awae/src/AWAE_Registry.txt

@@ -127,13 +127,13 @@ typedef   ^ MiscVarType    IntKi       iPlaneTurbTurb        {:}{:}{:} - -  "Fir
 typedef   ^ MiscVarType    IntKi       iPlaneTurbChunk       {:}{:}{:} - -  "First and Last plane index by source turbine and destination chunk index" -
 typedef   ^ MiscVarType    Logical     LowResChunkHasWake          {:} - -  "Low-res gridFirst and Last plane index by source turbine and destination chunk index" -
 typedef   ^ MiscVarType    ReKi        MaxWakePointSep              -  - -  "Maximum separation between wake points" -
-typedef   ^ MiscVarType    Logical  parallelFlag   {:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     r_s   {:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     r_e   {:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     rhat_s   {:}{:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     rhat_e   {:}{:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     pvec_cs   {:}{:}{:}    - -  "" -
-typedef   ^ MiscVarType    ReKi     pvec_ce   {:}{:}{:}    - -  "" -
+typedef   ^ MiscVarType    Logical  parallelFlag   {:}{:}    - -  "Flag indicating whether adjacent wake planes np and np+1 are (numerically) parallel for each turbine; dims: (plane-pair index np, turbine index nt)" -
+typedef   ^ MiscVarType    ReKi     r_s   {:}{:}    - -  "Perpendicular distance from the start wake-plane center p_plane(:,np,nt) to the intersection line of adjacent (non-parallel) wake planes np and np+1; dims: (plane-pair index np, turbine index nt)" m
+typedef   ^ MiscVarType    ReKi     r_e   {:}{:}    - -  "Perpendicular distance from the end wake-plane center p_plane(:,np+1,nt) to the intersection line of adjacent (non-parallel) wake planes np and np+1; dims: (plane-pair index np, turbine index nt)" m
+typedef   ^ MiscVarType    ReKi     rhat_s   {:}{:}{:}    - -  "Unit vector lying in the start wake plane (np) pointing from the plane-plane intersection line toward the start plane center; dims: (XYZ component, plane-pair index np, turbine index nt)" -
+typedef   ^ MiscVarType    ReKi     rhat_e   {:}{:}{:}    - -  "Unit vector lying in the end wake plane (np+1) pointing from the plane-plane intersection line toward the end plane center; dims: (XYZ component, plane-pair index np, turbine index nt)" -
+typedef   ^ MiscVarType    ReKi     pvec_cs   {:}{:}{:}    - -  "Closest point on the plane-plane intersection line to the start wake-plane center, p_plane(:,np,nt) - r_s*rhat_s; dims: (XYZ component, plane-pair index np, turbine index nt)" m
+typedef   ^ MiscVarType    ReKi     pvec_ce   {:}{:}{:}    - -  "Closest point on the plane-plane intersection line to the end wake-plane center, p_plane(:,np+1,nt) - r_e*rhat_e; dims: (XYZ component, plane-pair index np, turbine index nt)" m
 #
 typedef   ^ MiscVarType    SiKi     outVizXYPlane  {:}{:}{:}{:} - - "An array holding the output data for a 2D visualization slice" -
 typedef   ^ MiscVarType    SiKi     outVizYZPlane  {:}{:}{:}{:} - - "An array holding the output data for a 2D visualization slice" -
@@ -145,7 +145,10 @@ typedef   ^ MiscVarType    InflowWind_OutputType   y_IfW_High {:} - -  "InflowWi
 
 #wake added turbulence
 typedef   ^ MiscVarType    ReKi     V_amb_low_disk       {:}{:} - -  "Rotor averaged ambiend wind speed for each wind turbine (3 x nWT)" m/s
-typedef   ^ MiscVarType    IntKi    planeDomainExit      {:}{:} 0 -  "Value indicates edge number (0: still in domain, +/-1: +/-X, +/-2: +/-Y, +/-3: +/-Z) the plane crossed" -
+typedef   ^ MiscVarType    IntKi    planeDomainExit      {:}{:}{:} 0 -  "Per-dimension flag (0: still in domain, -1: crossed lower bound, +1: crossed upper bound) for each plane [dim,plane,turbine]" -
+
+#visualization
+typedef   ^ MiscVarType    IntKi    WakeVTK_StartN       {:}{:} - - "Time step when wake plane starts - counted by N_dtLow. Indices [wakenum,turbnum]" -
 
 # Low-resolution grid chunk data
 typedef   ^ LRGChunkType     IntKi    iChunk                  {3} -  -  "XYZ index of chunk" -
@@ -225,8 +228,12 @@ typedef  ^  ParameterType  IntKi            NOutDisWindXZ      -  - -   "Number
 typedef  ^  ParameterType  ReKi             OutDisWindY       {:} - -   "Y coordinates of XZ planes for output of disturbed wind data across the low-resolution domain [1 to NOutDisWindXZ]" meters
 typedef  ^  ParameterType  LOGICAL          OutDisWindYvalid  {:}   - - "Valid XZ planes for output of disturbed wind data across the low-resolution domain [1 to NOutDisWindXZ]" -
 typedef  ^  ParameterType  CHARACTER(1024)  OutFileRoot        -  - -   "The root name derived from the primary FAST.Farm input file" -
-typedef  ^  ParameterType  CHARACTER(1024)  OutFileVTKRoot     -  - -   "The root name for VTK outputs" -
+typedef  ^  ParameterType  CHARACTER(1024)  OutFileFFvtkRoot   -  - -   "The root name for VTK outputs" -
+typedef  ^  ParameterType  CHARACTER(1024)  OutFileFFvtkWakeRoot -  - -   "The root name for VTK outputs for wake planes" -
+typedef  ^  ParameterType  CHARACTER(1024)  OutFileFFvtkWakeNullData - - -   "Null data for an unpopulated wake plane" -
 typedef  ^  ParameterType  IntKi            VTK_tWidth         -  - -   "Number of characters for VTK timestamp outputs" -
+typedef  ^  ParameterType  IntKi            VTK_tWidthPlanes   -  0 -   "Number of charactes for the VTK plane numbers" -
+typedef  ^  ParameterType  LOGICAL          WrPlanes           - .false. - "Write plane data out" -
 #typedef  ^  ParameterType  OutParmType      OutParam          {:} - -  "Names and units (and other characteristics) of all requested output parameters"      -
 #wake added turbulence
 typedef  ^  ParameterType  Logical          WAT_Enabled        -  - -   "Switch for turning on and off wake-added turbulence" -
@@ -239,6 +246,7 @@ typedef   ^ OutputType     AWAE_HighWindGrid Vdist_High   {:}     -      -
 typedef   ^ OutputType     ReKi       V_plane      {:}{:}{:}              -      -      "Advection, deflection, and meandering velocity of wake planes for each turbine"      m/s
 typedef   ^ OutputType     ReKi       TI_amb       {:}                    -      -      "Ambient turbulence intensity of wind at rotor disk for each turbine"      -
 typedef   ^ OutputType     ReKi       Vx_wind_disk {:}                    -      -      "Rotor-disk-averaged ambient wind speed, normal to disk, for each turbine"      m/s
+typedef   ^ OutputType     ReKi       wakePlaneDomainExit {:}{:}{:}           0      -      "Per-dimension flag (0: in domain, -1: crossed lower, +1: crossed upper) for each plane [dim,plane,turbine]" -
 
 
 # ..... Inputs ...........................................................................................................