在二维 cfd 仿真中，应如何将边界条件应用于压力场？

问：

我正在编写一个 2D 不可压缩非定常纳维-斯托克斯求解器，以计算形状上方和下方的气流。该程序使用MAC方法，在速度的边界条件下使用带有幽灵单元的交错网格，我在边界条件（形状）内将所有值设置为零。问题在于边界条件未正确应用。没有一个迭代在边界内显示为 0 压力（但显示形状外的 x 速度的预期结果），最终所有迭代都显示所有像元的所有值都设置为零。

此外，对于大多数单元，在第 9 次迭代后，不满足连续性。

我尝试将形状内的压力以及所有其他值设置为 0。形状的边界平均了它之前/之后的速度，以强制执行无滑移条件。

除了形状内的细胞外，所有细胞都应该满足连续性，但大多数细胞（仅 30/3300）都没有满足连续性