Merge pull request #481 from firedrakeproject/equation_class_extruded…

…_mesh_check Added check for periodicity
firedrakeproject · Mar 5, 2024 · d4afd7c · d4afd7c
2 parents c49a85b + beb2df9
commit d4afd7c
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/gusto/equations.py b/gusto/equations.py
@@ -375,7 +375,8 @@ def set_no_normal_flow_bcs(self, domain, no_normal_flow_bc_ids):
                 + 'when there is a variable called "u" and none was found')
 
         Vu = domain.spaces("HDiv")
-        if Vu.extruded:
+        # we only apply no normal-flow BCs when extruded mesh is non periodic
+        if Vu.extruded and not Vu.ufl_domain().topology.extruded_periodic:
             self.bcs['u'].append(DirichletBC(Vu, 0.0, "bottom"))
             self.bcs['u'].append(DirichletBC(Vu, 0.0, "top"))
         for id in no_normal_flow_bc_ids: