Add HLL numerical flux

doc/misc.rst

@@ -71,3 +71,6 @@ References
     `(DOI) <https://doi.org/10.1007/978-3-642-59721-3_14>`__
 .. [Ihme_2014] Yu Lv and Matthias Ihme (2014) Journal of Computationsl Physics 270 105 \
     `(DOI) <http://dx.doi.org/10.1016/j.jcp.2014.03.029>`__
+.. [Toro_2009] Eleuterio F. Toro (2009), Riemann Solvers and Numerical Methods for Fluid Dynamics, Springer \
+    `(DOI) <http://doi.org/10.1007/978-3-540-49834-6>`__
+

examples/autoignition-mpi.py

@@ -553,6 +553,8 @@ def my_post_step(step, t, dt, state):
             logmgr.tick_after()
         return make_obj_array([cv, fluid_state.temperature]), dt
 
+    from mirgecom.inviscid import inviscid_flux_rusanov
+
     def my_rhs(t, state):
         cv, tseed = state
         from mirgecom.gas_model import make_fluid_state
@@ -561,6 +563,7 @@ def my_rhs(t, state):
         return make_obj_array([
             euler_operator(discr, state=fluid_state, time=t, boundaries=boundaries,
                            gas_model=gas_model,
+                           inviscid_numerical_flux_func=inviscid_flux_rusanov,
                            quadrature_tag=quadrature_tag)
             + eos.get_species_source_terms(cv, fluid_state.temperature),
             0*tseed])

mirgecom/boundary.py

@@ -42,7 +42,10 @@
 from meshmode.mesh import BTAG_ALL, BTAG_NONE  # noqa
 from mirgecom.fluid import make_conserved
 from grudge.trace_pair import TracePair
-from mirgecom.inviscid import inviscid_facial_flux
+from mirgecom.inviscid import (
+    inviscid_facial_flux,
+    inviscid_flux_rusanov
+)
 
 from abc import ABCMeta, abstractmethod
 
@@ -108,21 +111,26 @@ def _boundary_quantity(self, discr, btag, quantity, **kwargs):
         return discr.project(btag, btag.with_dtag("all_faces"), quantity)
 
     def inviscid_divergence_flux(self, discr, btag, gas_model, state_minus,
+                                 numerical_flux_func=inviscid_flux_rusanov,
                                  **kwargs):
         """Get the inviscid boundary flux for the divergence operator."""
         # This one is when the user specified a function that directly
         # prescribes the flux components at the boundary
         if self._inviscid_bnd_flux_func:
-            return self._inviscid_bnd_flux_func(discr, btag, gas_model, state_minus,
-                                                **kwargs)
+            return self._inviscid_bnd_flux_func(
+                discr, btag, gas_model, state_minus,
+                numerical_flux_func=numerical_flux_func, **kwargs)
 
         state_plus = self._bnd_state_func(discr=discr, btag=btag,
                                           gas_model=gas_model,
                                           state_minus=state_minus, **kwargs)
         boundary_state_pair = TracePair(btag, interior=state_minus,
                                         exterior=state_plus)
 
-        return self._inviscid_div_flux_func(discr, state_tpair=boundary_state_pair)
+        return self._inviscid_div_flux_func(discr, state_pair=boundary_state_pair,
+                                            gas_model=gas_model,
+                                            numerical_flux_func=numerical_flux_func,
+                                            **kwargs)
 
 
 class DummyBoundary(PrescribedFluidBoundary):

mirgecom/euler.py

@@ -63,8 +63,10 @@
 
 from mirgecom.inviscid import (
     inviscid_flux,
-    inviscid_facial_flux
+    inviscid_flux_rusanov,
+    inviscid_boundary_flux_for_divergence_operator,
 )
+
 from mirgecom.operators import div_operator
 from mirgecom.gas_model import (
     project_fluid_state,
@@ -81,6 +83,7 @@ class _EulerTseedTag:
 
 
 def euler_operator(discr, state, gas_model, boundaries, time=0.0,
+                   inviscid_numerical_flux_func=inviscid_flux_rusanov,
                    quadrature_tag=None):
     r"""Compute RHS of the Euler flow equations.
 
@@ -176,20 +179,13 @@ def _interp_to_surf_quad(utpair):
 
     # Compute volume contributions
     inviscid_flux_vol = inviscid_flux(vol_state_quad)
+
     # Compute interface contributions
-    inviscid_flux_bnd = (
-
-        # Interior faces
-        sum(inviscid_facial_flux(discr, state_pair)
-            for state_pair in interior_states_quad)
-
-        # Domain boundary faces
-        + sum(
-            boundaries[btag].inviscid_divergence_flux(
-                discr, as_dofdesc(btag).with_discr_tag(quadrature_tag), gas_model,
-                state_minus=boundary_states_quad[btag], time=time)
-            for btag in boundaries)
-    )
+    inviscid_flux_bnd = \
+        inviscid_boundary_flux_for_divergence_operator(
+            discr, gas_model, boundaries, interior_states_quad,
+            boundary_states_quad, quadrature_tag=quadrature_tag,
+            numerical_flux_func=inviscid_numerical_flux_func, time=time)
 
     return -div_operator(discr, dd_quad_vol, dd_quad_faces,
                          inviscid_flux_vol, inviscid_flux_bnd)