Hello,
The code seems correct, but removing the source terms from cells near the boundary means the volume over which you apply the source terms depends not on the geometry but on the mesh refinement.
Is this really what you need to do?
Regards,
Yvan
Several uniform heat source
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.