Use local limiting in free stream elixir

trixi-framework · May 15, 2024 · 98c2144 · 98c2144
1 parent 2ddd78a
commit 98c2144
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 12 deletions.
diff --git a/examples/structured_2d_dgsem/elixir_euler_free_stream_sc_subcell.jl b/examples/structured_2d_dgsem/elixir_euler_free_stream_sc_subcell.jl
@@ -14,8 +14,9 @@ volume_flux = flux_ranocha
 polydeg = 3
 basis = LobattoLegendreBasis(polydeg)
 limiter_idp = SubcellLimiterIDP(equations, basis;
-                                positivity_variables_cons = ["rho"],
-                                positivity_variables_nonlinear = [pressure])
+                                local_twosided_variables_cons = ["rho"],
+                                local_onesided_variables_nonlinear = [(Trixi.entropy_guermond_etal,
+                                                                       min)])
 
 volume_integral = VolumeIntegralSubcellLimiting(limiter_idp;
                                                 volume_flux_dg = volume_flux,
@@ -61,12 +62,12 @@ analysis_callback = AnalysisCallback(semi, interval = analysis_interval)
 
 alive_callback = AliveCallback(analysis_interval = analysis_interval)
 
-save_solution = SaveSolutionCallback(interval = 10000,
+save_solution = SaveSolutionCallback(interval = 100,
                                      save_initial_solution = true,
                                      save_final_solution = true,
                                      solution_variables = cons2prim)
 
-stepsize_callback = StepsizeCallback(cfl = 0.9)
+stepsize_callback = StepsizeCallback(cfl = 0.7)
 
 callbacks = CallbackSet(summary_callback,
                         analysis_callback, alive_callback,

diff --git a/test/test_structured_2d.jl b/test/test_structured_2d.jl
@@ -556,16 +556,16 @@ end
     @test_trixi_include(joinpath(EXAMPLES_DIR,
                                  "elixir_euler_free_stream_sc_subcell.jl"),
                         l2=[
-                            1.4663777294625118e-15,
-                            2.320054900530864e-14,
-                            3.487555722563465e-14,
-                            2.008802099296406e-14,
+                            2.5709100186298863e-16,
+                            2.184017068634721e-15,
+                            2.1781991538183746e-15,
+                            4.935868571076398e-15,
                         ],
                         linf=[
-                            2.3092638912203256e-14,
-                            2.0623780461193064e-13,
-                            2.6795232699328153e-13,
-                            2.362554596402333e-13,
+                            4.218847493575595e-15,
+                            2.7602919949742954e-14,
+                            3.533284775869561e-14,
+                            5.861977570020827e-14,
                         ],
                         atol=7.0e-13)
     # Ensure that we do not have excessive memory allocations