Fix, update some parameters in the tutorials.

loganoz · Apr 4, 2024 · e5cfcd6 · e5cfcd6
1 parent 6341555
commit e5cfcd6
Show file tree

Hide file tree

Showing 8 changed files with 18 additions and 11 deletions.
diff --git a/tutorials/Cylinder/Cylinder_Re100.control b/tutorials/Cylinder/Cylinder_Re100.control
@@ -1,23 +1,28 @@
 Flow equations        = "NS"
 mesh file name        = "MESH/cyl_circ.msh"
 Polynomial order      = 3
-Number of time steps  = 10000
-Output Interval       = 100
-autosave mode         = "iteration"
-autosave interval     = 2500
+Number of time steps  = 1000000
+Output Interval       = 5000
 Convergence tolerance = 1.d-10
 cfl                   = 0.6
 dcfl                  = 0.6
 mach number           = 0.2
 Reynolds number       = 100.0
 AOA theta             = 0.0
 AOA phi               = 0.0
+
 solution file name    = "RESULTS/Cylinder_Re100.hsol"
 save gradients with solution = .false.
 restart               = .true.
 restart file name     = "RESULTS/Cylinder_Re40.hsol"
 riemann solver        = roe
 simulation type       = time-accurate
+final time            = 200
+
+autosave mode         = "time"
+autosave interval     = 10.0
+
+
 time integration      = explicit
 explicit method       = RK3
 
@@ -57,4 +62,4 @@ explicit method       = RK3
    Name = wake_u
    Position = [2.0,0.0,0.5]
    Variable = u
-#end
+#end
diff --git a/tutorials/Cylinder/Tutorial_Cylinder.docx b/tutorials/Cylinder/Tutorial_Cylinder.docx
diff --git a/tutorials/Cylinder/Tutorial_Cylinder.pdf b/tutorials/Cylinder/Tutorial_Cylinder.pdf
diff --git a/tutorials/Cylinder/mesh_gmsh/Custom_mesh.geo b/tutorials/Cylinder/mesh_gmsh/Custom_mesh.geo
@@ -102,7 +102,10 @@ Physical Surface("inlet", 33) = {8};
 Physical Surface("outlet", 34) = {14};
 Physical Surface("up", 35) = {18};
 Physical Surface("down", 36) = {10};
-Physical Surface("walls", 37) = {13, 2, 17, 3, 4, 1, 9, 20};
-Physical Surface("cylinder", 38) = {19, 6, 16, 12};
+Physical Surface("wall1", 37) = {13, 17, 9, 20};
+Physical Surface("wall2", 38) = {2, 3, 4, 1};
+Physical Surface("cylinder", 39) = {19, 6, 16, 12};
 //The whole volume is limited by Surfaces 1, 2, 3 and 4.
-Physical Volume("fluid", 39) = {4, 3, 2, 1};
+Physical Volume("fluid", 40) = {4, 3, 2, 1};
+// uncoment next line to assing a order 3
+// Mesh.ElementOrder = 3; 
diff --git a/tutorials/Wind_Turbine/ActuatorDef/Act_ActuatorDef.dat b/tutorials/Wind_Turbine/ActuatorDef/Act_ActuatorDef.dat
@@ -278,8 +278,8 @@
 ############################
 # NUMERICAL PARAMETERS
 ############################
+2.5
 # Gaussian parameter epsilon
-0.0753
 # tip correction cosntants c1 c2
 0.125 21.0
 
diff --git a/tutorials/Wind_Turbine/Tutorial_Wind_Turbine.docx b/tutorials/Wind_Turbine/Tutorial_Wind_Turbine.docx
diff --git a/tutorials/Wind_Turbine/Tutorial_Wind_Turbine.pdf b/tutorials/Wind_Turbine/Tutorial_Wind_Turbine.pdf
diff --git a/tutorials/Wind_Turbine/wind_turbine.control b/tutorials/Wind_Turbine/wind_turbine.control
@@ -20,10 +20,9 @@ restart file name     = "RESULTS/actuatorP1.hsol"
 
 ! ACTUATOR LINE
 use actuatorline = .true.
-actuator epsilon type = 0
+actuator epsilon type = 2
 actuator calculate with projection = .false.
 actuator save average = .true.
-actuator initial azimuthal = .true.
 
 riemann solver         = "lax-friedrichs"