examples.diffusion.steadyState.mesh50x50.inputΒΆ

This input file again solves a 1D diffusion problem as in examples.diffusion.mesh1D. The difference being that the mesh is two dimensional.

The result is again tested in the same way:

>>> DiffusionTerm().solve(var)
>>> Lx = nx * dx
>>> x = mesh.cellCenters[0]
>>> analyticalArray = valueLeft + (valueRight - valueLeft) * x / Lx
>>> print(var.allclose(analyticalArray, rtol = 1e-9))
1
Last updated on Jun 26, 2024. Created using Sphinx 7.1.2.