Merge branch 'sdg_linear_transport_circle' of https://github.com/bcub…

…e-project/BcubeTutorials.jl into sdg_linear_transport_circle
bcube-project · May 30, 2024 · 59184fa · 59184fa
2 parents 008f8aa + 24469ea
commit 59184fa
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/src/example/transport_hypersurface/transport_hypersurface.jl b/src/example/transport_hypersurface/transport_hypersurface.jl
@@ -1072,8 +1072,8 @@ function scalar_torus(;
 
     b = Bcube.allocate_dofs(U)
     du = similar(b)
-    progressBar && (progress = Progress(nitemax))
-    for ite in 1:nitemax
+    progressBar && (progress = Progress(nite))
+    for ite in 1:nite
         b .= 0.0
 
         # Apply limitation
@@ -1097,7 +1097,7 @@ function scalar_torus(;
         progressBar && next!(progress)
 
         # Output results
-        if ite % (nitemax ÷ _nout) == 0
+        if ite % (nite ÷ _nout) == 0
             append_vtk(vtk, u, lim_u, t)
         end
     end
@@ -1139,15 +1139,15 @@ end
 # )
 @time scalar_torus(;
     degree = 1,
-    CFL = 0.5, # d=0, CFL=0.4 OK
+    CFL = 0.2, # d=0, CFL=0.4 OK
     rint = 1.0,
     rext = 1.5,
     lc = 0.08,
     ϕ = 0.5 * π / 2,
     C = 1.0,
-    tmax = 10.0,
+    tmax = 11.0, # 11s  = 1 turn for rint=1, rext=1.5
     nout = 100,
-    nitemax = 10000, # d = 4, n = 400 OK,
+    nitemax = 100000, # d = 4, n = 400 OK,
     isLimiterActive = false,
     progressBar = true,
     meshOrder = 2,