add davis wave speed estimate for upcoming change

trixi-framework · Feb 8, 2024 · 48ab702 · 48ab702
1 parent 3eb1b39
commit 48ab702
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/examples/structured_1d_dgsem/elixir_traffic_flow_lwr_greenlight.jl b/examples/structured_1d_dgsem/elixir_traffic_flow_lwr_greenlight.jl
@@ -7,7 +7,7 @@ using Trixi
 equations = TrafficFlowLWREquations1D()
 
 basis = LobattoLegendreBasis(3)
-surface_flux = flux_hll
+surface_flux = FluxHLL(min_max_speed_davis)
 solver = DGSEM(basis, surface_flux)
 
 coordinates_min = (-1.0,) # minimum coordinate

diff --git a/src/equations/traffic_flow_lwr_1d.jl b/src/equations/traffic_flow_lwr_1d.jl
@@ -88,6 +88,11 @@ end
     return λ_min, λ_max
 end
 
+@inline function min_max_speed_davis(u_ll, u_rr, orientation::Integer,
+                                     equations::TrafficFlowLWREquations1D)
+    min_max_speed_naive(u_ll, u_rr, orientation, equations)
+end
+
 @inline function max_abs_speeds(u, equations::TrafficFlowLWREquations1D)
     return (abs(equations.v_max * (1.0 - 2 * u[1])),)
 end