Multiphysics (production)

examples/combozzle-mpi.py

@@ -33,7 +33,7 @@
 from meshmode.array_context import PyOpenCLArrayContext
 
 from meshmode.mesh import BTAG_ALL, BTAG_NONE  # noqa
-from grudge.dof_desc import DTAG_BOUNDARY
+from grudge.dof_desc import BoundaryDomainTag
 from grudge.shortcuts import make_visualizer
 from grudge.dof_desc import DISCR_TAG_QUAD
 from mirgecom.discretization import create_discretization_collection
@@ -859,8 +859,8 @@ def get_fluid_state(cv, tseed):
     if not periodic:
         if multiple_boundaries:
             for idir in range(dim):
-                boundaries[DTAG_BOUNDARY(f"+{idir}")] = wall
-                boundaries[DTAG_BOUNDARY(f"-{idir}")] = wall
+                boundaries[BoundaryDomainTag(f"+{idir}")] = wall
+                boundaries[BoundaryDomainTag(f"-{idir}")] = wall
         else:
             boundaries = {BTAG_ALL: wall}
 

examples/doublemach-mpi.py

@@ -31,7 +31,7 @@
 from functools import partial
 
 from meshmode.mesh import BTAG_ALL, BTAG_NONE  # noqa
-from grudge.dof_desc import DTAG_BOUNDARY
+from grudge.dof_desc import BoundaryDomainTag
 from grudge.shortcuts import make_visualizer
 
 
@@ -239,22 +239,22 @@ def main(ctx_factory=cl.create_some_context, use_logmgr=True,
     eos = IdealSingleGas()
     gas_model = GasModel(eos=eos, transport=transport_model)
 
-    def _boundary_state(discr, btag, gas_model, state_minus, **kwargs):
+    def _boundary_state(discr, dd_bdry, gas_model, state_minus, **kwargs):
         actx = state_minus.array_context
-        bnd_discr = discr.discr_from_dd(btag)
+        bnd_discr = discr.discr_from_dd(dd_bdry)
         nodes = actx.thaw(bnd_discr.nodes())
         return make_fluid_state(initializer(x_vec=nodes, eos=gas_model.eos,
                                             **kwargs), gas_model)
 
     boundaries = {
-        DTAG_BOUNDARY("ic1"):
-        PrescribedFluidBoundary(boundary_state_func=_boundary_state),
-        DTAG_BOUNDARY("ic2"):
-        PrescribedFluidBoundary(boundary_state_func=_boundary_state),
-        DTAG_BOUNDARY("ic3"):
-        PrescribedFluidBoundary(boundary_state_func=_boundary_state),
-        DTAG_BOUNDARY("wall"): AdiabaticNoslipWallBoundary(),
-        DTAG_BOUNDARY("out"): AdiabaticNoslipWallBoundary(),
+        BoundaryDomainTag("ic1"):
+            PrescribedFluidBoundary(boundary_state_func=_boundary_state),
+        BoundaryDomainTag("ic2"):
+            PrescribedFluidBoundary(boundary_state_func=_boundary_state),
+        BoundaryDomainTag("ic3"):
+            PrescribedFluidBoundary(boundary_state_func=_boundary_state),
+        BoundaryDomainTag("wall"): AdiabaticNoslipWallBoundary(),
+        BoundaryDomainTag("out"): AdiabaticNoslipWallBoundary(),
     }
 
     if rst_filename:

examples/heat-source-mpi.py

@@ -30,7 +30,7 @@
 from meshmode.mesh import BTAG_ALL, BTAG_NONE  # noqa
 import grudge.op as op
 from grudge.shortcuts import make_visualizer
-from grudge.dof_desc import DTAG_BOUNDARY
+from grudge.dof_desc import BoundaryDomainTag
 from mirgecom.discretization import create_discretization_collection
 from mirgecom.integrators import rk4_step
 from mirgecom.diffusion import (
@@ -126,8 +126,8 @@ def main(actx_class, ctx_factory=cl.create_some_context, use_logmgr=True,
     nodes = actx.thaw(discr.nodes())
 
     boundaries = {
-        DTAG_BOUNDARY("dirichlet"): DirichletDiffusionBoundary(0.),
-        DTAG_BOUNDARY("neumann"): NeumannDiffusionBoundary(0.)
+        BoundaryDomainTag("dirichlet"): DirichletDiffusionBoundary(0.),
+        BoundaryDomainTag("neumann"): NeumannDiffusionBoundary(0.)
     }
 
     u = discr.zeros(actx)
@@ -181,7 +181,7 @@ def rhs(t, u):
             set_dt(logmgr, dt)
             logmgr.tick_after()
     final_answer = actx.to_numpy(op.norm(discr, u, np.inf))
-    resid = abs(final_answer - 0.00020620711665201585)
+    resid = abs(final_answer - 0.0002062062188374177)
     if resid > 1e-15:
         raise ValueError(f"Run did not produce the expected result {resid=}")
 

examples/hotplate-mpi.py

@@ -32,7 +32,7 @@
 from meshmode.mesh import BTAG_ALL, BTAG_NONE  # noqa
 
 from grudge.shortcuts import make_visualizer
-from grudge.dof_desc import DTAG_BOUNDARY
+from grudge.dof_desc import BoundaryDomainTag
 
 from mirgecom.discretization import create_discretization_collection
 from mirgecom.fluid import make_conserved
@@ -223,21 +223,22 @@ def tramp_2d(x_vec, eos, cv=None, **kwargs):
 
     exact = initializer(x_vec=nodes, eos=gas_model.eos)
 
-    def _boundary_state(discr, btag, gas_model, state_minus, **kwargs):
+    def _boundary_state(discr, dd_bdry, gas_model, state_minus, **kwargs):
         actx = state_minus.array_context
-        bnd_discr = discr.discr_from_dd(btag)
+        bnd_discr = discr.discr_from_dd(dd_bdry)
         nodes = actx.thaw(bnd_discr.nodes())
         return make_fluid_state(initializer(x_vec=nodes, eos=gas_model.eos,
                                             **kwargs), gas_model)
 
-    boundaries = {DTAG_BOUNDARY("-1"):
-                  PrescribedFluidBoundary(boundary_state_func=_boundary_state),
-                  DTAG_BOUNDARY("+1"):
-                  PrescribedFluidBoundary(boundary_state_func=_boundary_state),
-                  DTAG_BOUNDARY("-2"): IsothermalWallBoundary(
-                      wall_temperature=bottom_boundary_temperature),
-                  DTAG_BOUNDARY("+2"): IsothermalWallBoundary(
-                      wall_temperature=top_boundary_temperature)}
+    boundaries = {
+        BoundaryDomainTag("-1"): PrescribedFluidBoundary(
+            boundary_state_func=_boundary_state),
+        BoundaryDomainTag("+1"): PrescribedFluidBoundary(
+            boundary_state_func=_boundary_state),
+        BoundaryDomainTag("-2"): IsothermalWallBoundary(
+            wall_temperature=bottom_boundary_temperature),
+        BoundaryDomainTag("+2"): IsothermalWallBoundary(
+            wall_temperature=top_boundary_temperature)}
 
     if rst_filename:
         current_t = restart_data["t"]

examples/lump-mpi.py

@@ -181,9 +181,9 @@ def main(actx_class, ctx_factory=cl.create_some_context, use_logmgr=True,
     from mirgecom.gas_model import GasModel, make_fluid_state
     gas_model = GasModel(eos=eos)
 
-    def boundary_solution(discr, btag, gas_model, state_minus, **kwargs):
+    def boundary_solution(discr, dd_bdry, gas_model, state_minus, **kwargs):
         actx = state_minus.array_context
-        bnd_discr = discr.discr_from_dd(btag)
+        bnd_discr = discr.discr_from_dd(dd_bdry)
         nodes = actx.thaw(bnd_discr.nodes())
         return make_fluid_state(initializer(x_vec=nodes, eos=gas_model.eos,
                                             **kwargs), gas_model)

examples/mixture-mpi.py

@@ -214,9 +214,9 @@ def main(actx_class, ctx_factory=cl.create_some_context, use_logmgr=True,
     initializer = MixtureInitializer(dim=dim, nspecies=nspecies,
                                      massfractions=y0s, velocity=velocity)
 
-    def boundary_solution(discr, btag, gas_model, state_minus, **kwargs):
+    def boundary_solution(discr, dd_bdry, gas_model, state_minus, **kwargs):
         actx = state_minus.array_context
-        bnd_discr = discr.discr_from_dd(btag)
+        bnd_discr = discr.discr_from_dd(dd_bdry)
         nodes = actx.thaw(bnd_discr.nodes())
         return make_fluid_state(initializer(x_vec=nodes, eos=gas_model.eos,
                                             **kwargs), gas_model,