use mesh property n

festim-dev · Oct 17, 2023 · aca8ebb · aca8ebb
1 parent b0901a0
commit aca8ebb
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/festim/hydrogen_transport_problem.py b/festim/hydrogen_transport_problem.py
@@ -267,11 +267,10 @@ def run(self, final_time: float):
 
             cm = self.species[0].solution
             # TODO this should be a property of Mesh
-            n = FacetNormal(self.mesh.mesh)
             D = self.subdomains[0].material.get_diffusion_coefficient(
                 self.mesh.mesh, self.temperature
             )
-            surface_flux = form(D * dot(grad(cm), n) * self.ds(2))
+            surface_flux = form(D * dot(grad(cm), self.mesh.n) * self.ds(2))
             flux = assemble_scalar(surface_flux)
             flux_values.append(flux)
             times.append(t)