Time depepdent problem in work

Author: niravshah241

demo/time_dependent_diffusion/dolfinx_time_dependent_diffusion.py

@@ -21,7 +21,7 @@
 uD_sol = dolfinx.fem.Function(VD)
 uD_prev = dolfinx.fem.Function(VD)
 # uD_prev.x.array[:] = 0.2 # Initial value
-uD_prev.interpolate(lambda x: 5 * np.sin(x[0] * 2 * np.pi) * np.cos(x[1] * 2 * np.pi))
+uD_prev.interpolate(lambda x: 5. * np.sin(x[0] * 2 * np.pi) * np.cos(x[1] * 2 * np.pi))
 D = ufl.as_tensor([[1., 0., 0.],[0., 1.2, 0.],[0., 0., 3.]])
 f_source = dolfinx.fem.Constant(mesh, PETSc.ScalarType(0.))
 flux = dolfinx.fem.Constant(mesh, PETSc.ScalarType(0.))
@@ -65,7 +65,7 @@
     uD_prev.x.array[:] = uD_sol.x.array
 
     xdmf.write_function(uD_sol, t)
-    print(mesh.comm.allreduce(dolfinx.fem.assemble_scalar(dolfinx.fem.form(ufl.inner(uD_sol.dx(2), uD_sol.dx(2)) * ufl.dx)), op=MPI.MAX))
+    print(mesh.comm.allreduce(dolfinx.fem.assemble_scalar(dolfinx.fem.form(ufl.inner(uD_sol.dx(2), uD_sol.dx(2)) * ufl.dx)), op=MPI.SUM))
 
 xdmf.close()
 

demo/time_dependent_diffusion/mesh_data/mesh.py

@@ -5,13 +5,13 @@
 gmsh.initialize(" ", False)
 gdim = 3
 
-center_x, center_y, center_z, radius = 0., 0., 0., 2.
+center_x, center_y, center_z, radius = 0., 0., -1., 2.
 circle = gmsh.model.occ.addCircle(center_x, center_y, center_z, radius)
 circle_loop = gmsh.model.occ.addCurveLoop([circle])
 surface = gmsh.model.occ.addPlaneSurface([circle_loop])
 gmsh.model.occ.synchronize()
 
-h = 1
+h = 2
 extrusion = gmsh.model.occ.extrude([(2, surface)], 0, 0, h)
 gmsh.model.occ.synchronize()