Implemented GLM divergence cleaning for multi-ion MHD

amrueda · Dec 19, 2023 · 31d27ad · 31d27ad
1 parent a908fd7
commit 31d27ad
Show file tree

Hide file tree

Showing 3 changed files with 94 additions and 29 deletions.
diff --git a/examples/tree_2d_dgsem/elixir_mhdmultiion_ec.jl b/examples/tree_2d_dgsem/elixir_mhdmultiion_ec.jl
@@ -40,12 +40,17 @@ save_solution = SaveSolutionCallback(dt = 0.1, # interval=100,
                                      save_final_solution = true,
                                      solution_variables = cons2prim)
 
-stepsize_callback = StepsizeCallback(cfl = 0.5)
+cfl = 0.5
+
+stepsize_callback = StepsizeCallback(cfl = cfl)
+
+glm_speed_callback = GlmSpeedCallback(glm_scale = 0.5, cfl = cfl)
 
 callbacks = CallbackSet(summary_callback,
                         analysis_callback, alive_callback,
                         save_solution,
-                        stepsize_callback)
+                        stepsize_callback,
+                        glm_speed_callback)
 
 ###############################################################################
 # run the simulation

diff --git a/src/callbacks_step/glm_speed_dg.jl b/src/callbacks_step/glm_speed_dg.jl
@@ -7,7 +7,8 @@
 
 function calc_dt_for_cleaning_speed(cfl::Real, mesh,
                                     equations::Union{AbstractIdealGlmMhdEquations,
-                                                     AbstractIdealGlmMhdMulticomponentEquations},
+                                                     AbstractIdealGlmMhdMulticomponentEquations,
+                                                     IdealMhdMultiIonEquations2D},
                                     dg::DG, cache)
     # compute time step for GLM linear advection equation with c_h=1 for the DG discretization on
     # Cartesian meshes
@@ -20,7 +21,8 @@ end
 
 function calc_dt_for_cleaning_speed(cfl::Real, mesh,
                                     equations::Union{AbstractIdealGlmMhdEquations,
-                                                     AbstractIdealGlmMhdMulticomponentEquations},
+                                                     AbstractIdealGlmMhdMulticomponentEquations,
+                                                     IdealMhdMultiIonEquations2D},
                                     dg::DGMulti, cache)
     rd = dg.basis
     md = mesh.md