我想通过使用制造的解决方案来验证不可压缩 Navier-Stokes 方程的有限体积求解器（SIMPLE-Algorithm）。我对所有边界处的速度使用狄利克雷边界条件。我使用的速度的制造解决方案被构造为矢量场的旋度，从而满足任意点的连续性方程。

根据我如何选择求解方程的域，我遇到了收敛问题：压力校正方程的残差在初始减小后停滞不前。

我正在使用中点规则近似通过边界面的质量通量，显然这会导致考虑整个问题域的质量净损失/增加。

这是正常行为吗？除了选择在边界处同样消失的制造解决方案之外，还有其他补救措施吗？是否可以选择提供精确的质量通量而不是近似值（不确定这是否允许，因为 Dirichlet 边界条件仅固定速度）？