Fix_ElixirDGMulti_ShuOsher (#2618)

Co-authored-by: Hendrik Ranocha <[email protected]>

Author: DanielDoehring

examples/dgmulti_1d/elixir_euler_shu_osher_gauss_shock_capturing.jl

@@ -32,14 +32,7 @@ end
 
 initial_condition = initial_condition_shu_osher
 
-# Up to version 0.13.0, `max_abs_speed_naive` was used as the default wave speed estimate of
-# `const flux_lax_friedrichs = FluxLaxFriedrichs(), i.e., `FluxLaxFriedrichs(max_abs_speed = max_abs_speed_naive)`.
-# In the `StepsizeCallback`, though, the less diffusive `max_abs_speeds` is employed which is consistent with `max_abs_speed`.
-# Thus, we exchanged in PR#2458 the default wave speed used in the LLF flux to `max_abs_speed`.
-# To ensure that every example still runs we specify explicitly `FluxLaxFriedrichs(max_abs_speed_naive)`.
-# We remark, however, that the now default `max_abs_speed` is in general recommended due to compliance with the 
-# `StepsizeCallback` (CFL-Condition) and less diffusion.
-surface_flux = FluxLaxFriedrichs(max_abs_speed_naive)
+surface_flux = flux_lax_friedrichs
 volume_flux = flux_ranocha
 
 polydeg = 3
@@ -92,7 +85,7 @@ save_solution = SaveSolutionCallback(interval = 100,
                                      solution_variables = cons2prim)
 
 # handles the re-calculation of the maximum Δt after each time step
-stepsize_callback = StepsizeCallback(cfl = 0.2)
+stepsize_callback = StepsizeCallback(cfl = 0.13)
 
 # collect all callbacks such that they can be passed to the ODE solver
 callbacks = CallbackSet(summary_callback, analysis_callback, save_solution,

test/test_dgmulti_1d.jl

@@ -57,16 +57,15 @@ end
 @trixi_testset "elixir_euler_shu_osher_gauss_shock_capturing.jl " begin
     @test_trixi_include(joinpath(EXAMPLES_DIR,
                                  "elixir_euler_shu_osher_gauss_shock_capturing.jl"),
-                        cells_per_dimension=(64,), tspan=(0.0, 1.0),
                         l2=[
-                            1.6967151731067875,
-                            6.018445633981826,
-                            21.77425594743242
+                            1.696712726264938,
+                            6.018435800027037,
+                            21.774221602026298
                         ],
                         linf=[
-                            3.2229876650556477,
-                            10.702690533393842,
-                            38.37424900889908
+                            3.2229930998898952,
+                            10.702431261492814,
+                            38.37420027341893
                         ])
     # Ensure that we do not have excessive memory allocations
     # (e.g., from type instabilities)