Update NREL5MW ALM Rigid case

.website_src/amr-wind/actuator_line/NREL5MW_ALM_RIGID/orphan.rst

@@ -3,15 +3,15 @@
 .. toctree::
    :hidden:
 
-   results/README.md
+   postprocessing/OpenFAST_SectionalLoading.ipynb
    postprocessing/AVGPlanes.ipynb
    postprocessing/InstantaneousAvgPlanes.ipynb
-   setup/SetupALM_OFv402_RIGID.ipynb
+   ../../atmospheric_boundary_layer/convective_abl_nrel5mw/README.md
    performance/README.md
+   postprocessing/README.md
+   setup/SetupALM_OFv402_RIGID.ipynb
    README.md
-   setup/README.md
    postprocessing/OpenFAST_v40_Results.ipynb
-   postprocessing/OpenFAST_SectionalLoading.ipynb
-   postprocessing/README.md
-   ../../atmospheric_boundary_layer/convective_abl_nrel5mw/README.md
+   setup/README.md
+   results/README.md
 

.website_src/amr-wind/actuator_line/NREL5MW_ALM_RIGID/results/README.md

@@ -23,7 +23,7 @@ The time-history and time-averaged turbine performance is extracted and calculat
 
 |BldPitch1    |RotSpeed    |RotThrust  |RotTorq       |Electrical Power|
 |---          | ----       |---        |---           |--- |
-|0.000 deg    |12.100 rpm  |668.802 kN |3240.191 kN-m |3875.757 kW  |
+|0.000 deg    |12.100 rpm  |696.604 kN |3535.028 kN-m |4228.427 kW  |
 
 Note that the electrical power was calculated using this formula
 

.website_src/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/SetupALM_OFv402_RIGID.ipynb

@@ -141,8 +141,9 @@
      "text": [
       "GITHASH = df9ee7d\n",
       "Editing NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/5MW_Land_DLL_WTurb_cpp.fst\n",
-      "Editing /projects/wind_uq/lcheung/exawind-benchmarks.lcheung/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat\n",
-      "Editing /projects/wind_uq/lcheung/exawind-benchmarks.lcheung/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/NRELOffshrBsline5MW_Onshore_AeroDyn.dat\n",
+      "['EDFile', 'AeroFile']\n",
+      "Editing /projects/wind_uq/lcheung/exawind-benchmarks.redoALM/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat\n",
+      "Editing /projects/wind_uq/lcheung/exawind-benchmarks.redoALM/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/NRELOffshrBsline5MW_Onshore_AeroDyn.dat\n",
       "NREL5MW_v402_RIGID.yaml\n",
       "Saved turbine setup to NREL5MW_v402_RIGID.yaml\n",
       "turbinetype_name = 'NREL5MW_v402_ALM_RIGID'\n",
@@ -153,8 +154,8 @@
       "Actuator_hub_height = 90\n",
       "Actuator_num_points_blade = 25\n",
       "Actuator_num_points_tower = 11\n",
-      "Actuator_epsilon = [1.5, 1.5, 1.5]\n",
-      "Actuator_epsilon_tower = [1.5, 1.5, 1.5]\n",
+      "Actuator_epsilon = [2.0, 2.0, 2.0]\n",
+      "Actuator_epsilon_tower = [2.0, 2.0, 2.0]\n",
       "Actuator_openfast_start_time = 0.0\n",
       "Actuator_openfast_stop_time = 2000.0\n",
       "Actuator_nacelle_drag_coeff = 0.0\n",
@@ -189,7 +190,8 @@
       "       0.0 BlPitch(3) - Blade 3 initial pitch (degrees) [unused for 2 blades] [EDITED]\n",
       "      12.1 RotSpeed - Initial or fixed rotor speed (rpm) [EDITED]\n",
       "      30.0 NacYaw - Initial or fixed nacelle-yaw angle (degrees) [EDITED]\n",
-      "         0 Wake_Mod - Wake/induction model (switch) {0=none 1=BEMT 3=OLAF} [Wake_Mod cannot be 2 or 3 when linearizing] [EDITED]\n"
+      "         0 Wake_Mod - Wake/induction model (switch) {0=none 1=BEMT 3=OLAF} [Wake_Mod cannot be 2 or 3 when linearizing] [EDITED]\n",
+      "         0 UA_Mod - Unsteady Aero Model Switch (switch) {0=Quasi-steady (no UA) 2=B-L Gonzalez 3=B-L Minnema/Pierce 4=B-L HGM 4-states 5=B-L HGM+vortex 5 states 6=Oye 7=Boeing-Vertol} [EDITED]\n"
      ]
     }
    ],
@@ -251,7 +253,8 @@
     {
      "data": {
       "text/plain": [
-       "OrderedDict([('io.line_plot_int', '1'),\n",
+       "OrderedDict([('io.KE_int', '-1'),\n",
+       "             ('io.line_plot_int', '1'),\n",
        "             ('CoriolisForcing.turn_off_vertical_force', 'True')])"
       ]
      },
@@ -426,7 +429,7 @@
    "outputs": [],
    "source": [
     "turbbase_xy = [[1800, 1800]]\n",
-    "eps         = 1.5"
+    "eps         = 2.0"
    ]
   },
   {
@@ -442,7 +445,7 @@
       "\n",
       "# CSV file should have columns with\n",
       "# name, x, y, type, yaw, hubheight, options\n",
-      "T0, 1800.000000, 1800.000000, NREL5MW_v402_ALM_RIGID, 240.0, ,FSTparam_DT:0.0008608815426997245 EDparam_YawDOF:False AMRparam_Actuator_epsilon:[1.5;1.5;1.5] AMRparam_Actuator_epsilon_tower:[1.5;1.5;1.5] AMRparam_Actuator_openfast_stop_time:3000.0 \n",
+      "T0, 1800.000000, 1800.000000, NREL5MW_v402_ALM_RIGID, 240.0, ,FSTparam_DT:0.0008608815426997245 EDparam_YawDOF:False AMRparam_Actuator_epsilon:[2.0;2.0;2.0] AMRparam_Actuator_epsilon_tower:[2.0;2.0;2.0] AMRparam_Actuator_openfast_stop_time:3000.0 \n",
       "\n"
      ]
     }
@@ -481,10 +484,11 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "docopy = True from /projects/wind_uq/lcheung/exawind-benchmarks.lcheung/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID to T0_NREL5MW_v402_RIGID\n",
+      "docopy = True from /projects/wind_uq/lcheung/exawind-benchmarks.redoALM/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID to T0_NREL5MW_v402_RIGID\n",
+      "copy /projects/wind_uq/lcheung/exawind-benchmarks.redoALM/amr-wind/actuator_line/NREL5MW_ALM_RIGID/setup/NREL5MW_v402_RIGID failed\n",
       "T0_NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/5MW_Land_DLL_WTurb_cpp.fst\n",
-      "Setting Actuator_epsilon to [1.5, 1.5, 1.5]\n",
-      "Setting Actuator_epsilon_tower to [1.5, 1.5, 1.5]\n",
+      "Setting Actuator_epsilon to [2.0, 2.0, 2.0]\n",
+      "Setting Actuator_epsilon_tower to [2.0, 2.0, 2.0]\n",
       "Setting Actuator_openfast_stop_time to 3000.0\n",
       "OrderedDict([('DT', 0.0008608815426997245)])\n",
       "OrderedDict([('YawDOF', 'False')])\n"
@@ -494,8 +498,8 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "0.0008608815426997245 DT - Recommended module time step (s) [EDITED] [EDITED]\n",
-      "     False YawDOF - Yaw DOF (flag) [EDITED] [EDITED]\n"
+      "0.0008608815426997245 DT - Recommended module time step (s) [EDITED] [EDITED] [EDITED]\n",
+      "     False YawDOF - Yaw DOF (flag) [EDITED] [EDITED] [EDITED]\n"
      ]
     }
    ],
@@ -848,6 +852,7 @@
       "[ PASS] Actuator FST:T0      [T0_NREL5MW_v402_RIGID/openfast-cpp/5MW_Land_DLL_WTurb_cpp/NRELOffshrBsline5MW_Onshore_AeroDyn.dat] exists\n",
       "[ PASS] Actuator FST:T0      WakeMod=0 OK\n",
       "[ PASS] Actuator FST:T0      AirDens=1.000000, matches incflo.density=1.000000\n",
+      "[ SKIP] Actuator FST:T0      Skipping DISCON density check\n",
       "[ PASS] Sampling probes:virtualmast is inside domain\n",
       "[ PASS] Sampling probes:T0_rotorplaneUP \n",
       "[ PASS] Sampling probes:T0_rotorplaneDN \n",
@@ -865,7 +870,7 @@
       "\n",
       "Results: \n",
       " 25 PASS\n",
-      " 0 SKIP\n",
+      " 1 SKIP\n",
       " 0 FAIL\n",
       " 0 WARN\n"
      ]
@@ -1009,7 +1014,6 @@
       "ABL.cutoff_height                        = 50.0                \n",
       "time.plot_interval                       = 10000               \n",
       "io.plot_file                             = plt                 \n",
-      "io.KE_int                                = -1                  \n",
       "\n",
       "#---- tagging defs ----\n",
       "tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone\n",
@@ -1055,8 +1059,8 @@
       "Actuator.T0.hub_height                   = 90.0                \n",
       "Actuator.T0.num_points_blade             = 25                  \n",
       "Actuator.T0.num_points_tower             = 11                  \n",
-      "Actuator.T0.epsilon                      = 1.5 1.5 1.5         \n",
-      "Actuator.T0.epsilon_tower                = 1.5 1.5 1.5         \n",
+      "Actuator.T0.epsilon                      = 2.0 2.0 2.0         \n",
+      "Actuator.T0.epsilon_tower                = 2.0 2.0 2.0         \n",
       "Actuator.T0.openfast_start_time          = 0.0                 \n",
       "Actuator.T0.openfast_stop_time           = 3000.0              \n",
       "Actuator.T0.fllc                         = false               \n",
@@ -1218,6 +1222,7 @@
       "XYdomain_153_.Farm_XYdomain153.offset_vector = 0.0 0.0 0.0         \n",
       "\n",
       "#---- extra params ----\n",
+      "io.KE_int                                = -1                  \n",
       "io.line_plot_int                         = 1                   \n",
       "CoriolisForcing.turn_off_vertical_force  = True                \n",
       "#== END AMR-WIND INPUT ==\n",
@@ -1347,7 +1352,6 @@
       "ABL.cutoff_height                        = 50.0                \n",
       "time.plot_interval                       = 10000               \n",
       "io.plot_file                             = plt                 \n",
-      "io.KE_int                                = -1                  \n",
       "\n",
       "#---- tagging defs ----\n",
       "tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone\n",
@@ -1533,6 +1537,7 @@
       "XYdomain_153_.Farm_XYdomain153.offset_vector = 0.0 0.0 0.0         \n",
       "\n",
       "#---- extra params ----\n",
+      "io.KE_int                                = -1                  \n",
       "io.line_plot_int                         = 1                   \n",
       "CoriolisForcing.turn_off_vertical_force  = True                \n",
       "#== END AMR-WIND INPUT ==\n",
@@ -1662,7 +1667,6 @@
       "ABL.cutoff_height                        = 50.0                \n",
       "time.plot_interval                       = 10000               \n",
       "io.plot_file                             = plt                 \n",
-      "io.KE_int                                = -1                  \n",
       "\n",
       "#---- tagging defs ----\n",
       "tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone\n",

amr-wind/actuator_line/NREL5MW_ALM_RIGID/input_files/NREL5MW_ALM_RIGID_OFv402.inp

@@ -94,7 +94,6 @@ ABL.theta_amplitude                      = 0.8
 ABL.cutoff_height                        = 50.0                
 time.plot_interval                       = 10000               
 io.plot_file                             = plt                 
-io.KE_int                                = -1                  
 
 #---- tagging defs ----
 tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone
@@ -140,8 +139,8 @@ Actuator.T0.rotor_diameter               = 126.0
 Actuator.T0.hub_height                   = 90.0                
 Actuator.T0.num_points_blade             = 25                  
 Actuator.T0.num_points_tower             = 11                  
-Actuator.T0.epsilon                      = 1.5 1.5 1.5         
-Actuator.T0.epsilon_tower                = 1.5 1.5 1.5         
+Actuator.T0.epsilon                      = 2.0 2.0 2.0         
+Actuator.T0.epsilon_tower                = 2.0 2.0 2.0         
 Actuator.T0.openfast_start_time          = 0.0                 
 Actuator.T0.openfast_stop_time           = 3000.0              
 Actuator.T0.fllc                         = false               
@@ -303,6 +302,7 @@ XYdomain_153_.Farm_XYdomain153.axis2     = 0.0 5119.9998 0.0
 XYdomain_153_.Farm_XYdomain153.offset_vector = 0.0 0.0 0.0         
 
 #---- extra params ----
+io.KE_int                                = -1                  
 io.line_plot_int                         = 1                   
 CoriolisForcing.turn_off_vertical_force  = True                
 #== END AMR-WIND INPUT ==

amr-wind/actuator_line/NREL5MW_ALM_RIGID/input_files/NREL5MW_ALM_RIGID_noturb.inp

@@ -94,7 +94,6 @@ ABL.theta_amplitude                      = 0.8
 ABL.cutoff_height                        = 50.0                
 time.plot_interval                       = 10000               
 io.plot_file                             = plt                 
-io.KE_int                                = -1                  
 
 #---- tagging defs ----
 tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone
@@ -280,6 +279,7 @@ XYdomain_153_.Farm_XYdomain153.axis2     = 0.0 5119.9998 0.0
 XYdomain_153_.Farm_XYdomain153.offset_vector = 0.0 0.0 0.0         
 
 #---- extra params ----
+io.KE_int                                = -1                  
 io.line_plot_int                         = 1                   
 CoriolisForcing.turn_off_vertical_force  = True                
 #== END AMR-WIND INPUT ==

amr-wind/actuator_line/NREL5MW_ALM_RIGID/input_files/NREL5MW_BR.inp

@@ -94,7 +94,6 @@ ABL.theta_amplitude                      = 0.8
 ABL.cutoff_height                        = 50.0                
 time.plot_interval                       = 10000               
 io.plot_file                             = plt                 
-io.KE_int                                = -1                  
 
 #---- tagging defs ----
 tagging.labels                           = T0_level_0_zone T0_level_1_zone T0_level_2_zone T0_level_3_zone

amr-wind/actuator_line/NREL5MW_ALM_RIGID/postprocessing/InstantaneousAvgPlanes.py

@@ -3,7 +3,7 @@
 
 # # Instantaneous planes from NREL5MW ALM case
 
-# In[ ]:
+# In[1]:
 
 
 # Add any possible locations of amr-wind-frontend here
@@ -22,7 +22,7 @@
 import postproengine as ppeng
 
 
-# In[ ]:
+# In[2]:
 
 
 # Load ruamel or pyyaml as needed
@@ -54,14 +54,14 @@ def formatplot(fig, ax):
 ppeng.formatplot = formatplot
 
 
-# In[ ]:
+# In[3]:
 
 
-replacedict={'RUNDIR':'/gpfs/lcheung/HFM/exawind-benchmarks/NREL5MW_ALM_RIGID_OFv402/',
+replacedict={'RUNDIR':'/tscratch/lcheung/HFM/exawind-benchmarks/NREL5MW_ALM_RIGID_OFv402.redo/'
             }
 
 
-# In[ ]:
+# In[4]:
 
 
 yamlstring="""
@@ -165,14 +165,14 @@ def formatplot(fig, ax):
 yamldict = Loader(f, **loaderkwargs)
 
 
-# In[ ]:
+# In[5]:
 
 
 # Run the driver
 ppeng.driver(yamldict, verbose=True)
 
 
-# In[ ]:
+# In[6]:
 
 
 # Write out the notebook to a python script

amr-wind/actuator_line/NREL5MW_ALM_RIGID/postprocessing/OpenFAST_SectionalLoading.py

@@ -133,14 +133,14 @@ def makeSecBladeDF(csvfile, rpts, bladekeysdict):
 yamldict = Loader(f, **loaderkwargs)
 
 
-# In[ ]:
+# In[9]:
 
 
 # Run the driver
 ppeng.driver(yamldict, verbose=True)
 
 
-# In[ ]:
+# In[10]:
 
 
 d = makeSecBladeDF(replacedict['RESULTSDIR']+'/NREL5MW_SECLOADS_mean.csv', bladedat[:,0], blistdict)
@@ -149,7 +149,7 @@ def makeSecBladeDF(csvfile, rpts, bladekeysdict):
 
 # ## Plot the blade loading
 
-# In[ ]:
+# In[11]:
 
 
 yamlstring="""
@@ -201,14 +201,14 @@ def makeSecBladeDF(csvfile, rpts, bladekeysdict):
 yamldict = Loader(f, **loaderkwargs)
 
 
-# In[ ]:
+# In[12]:
 
 
 # Run the driver
 ppeng.driver(yamldict, verbose=True)
 
 
-# In[ ]:
+# In[13]:
 
 
 # Write out the notebook to a python script

amr-wind/actuator_line/NREL5MW_ALM_RIGID/postprocessing/OpenFAST_v40_Results.py

@@ -57,7 +57,7 @@ def stringReplaceDict(s, dreplace):
 #rundir='/gpfs/lcheung/HFM/exawind-benchmarks/NREL5MW_ALM_BD_OFv400'
 #RESULTSDIR='../results/OpenFAST_v40_out'
 
-replacedict={'RUNDIR':'/gpfs/lcheung/HFM/exawind-benchmarks/NREL5MW_ALM_RIGID_OFv402/',
+replacedict={'RUNDIR':'/tscratch/lcheung/HFM/exawind-benchmarks/NREL5MW_ALM_RIGID_OFv402.redo/',
              'RESULTSDIR':'../results/OpenFAST_v402_out',
             }