编写2D不可压缩NS求解器时压力场边界条件及连续性问题求助

1. 固体区域边界条件的错误实现 #

• 直接将外形内部所有值设为0的逻辑不成立：MAC网格下固体区域需要分离速度约束与压力求解屏蔽：
  • 固体内部的交错速度单元（u、v位置）应直接固定为0（无滑移），禁止迭代更新覆盖；
  • 压力泊松方程求解时，需将固体内部压力单元从求解域中剔除——比如将对应位置的系数矩阵设为单位值，右端项等于当前压力值，避免修正值传入固体区域。你观察到的固体内部压力非零，就是因为压力求解没有屏蔽这些单元，迭代修正依然会渗透进去。

2. 无滑移边界的Ghost Cell处理误区 #

• 边界前后速度取平均的方式不符合MAC网格的无滑移要求，正确处理需区分速度分量位置：
  • 垂直壁面附近的u速度单元：ghost cell的u值应设为 -u_邻域单元值（镜像法，保证壁面处合成速度为0）；
  • 水平壁面附近的v速度单元：ghost cell的v值设为 -v_邻域单元值；
  • 斜壁面需通过几何关系推导ghost cell速度，确保壁面法向/切向速度均为0，而非简单平均。错误的边界处理会引入非物理通量，直接破坏连续性条件。

3. 连续性条件不满足的核心原因 #

• 迭代后期大部分单元不满足连续性，本质是压力修正未有效消除速度散度：
  • 检查压力泊松方程的边界条件：固体壁面处应施加法向压力梯度为0的约束（或根据速度边界推导），而非直接设0；
  • 确认压力求解器的收敛性：若压力方程求解不充分，散度误差会持续累积，最终导致数值崩溃；
  • 核对速度更新公式：MAC方法中速度更新依赖压力梯度，需确保网格步长、1/(ρΔt)等系数计算准确，避免数值误差放大。

4. 全场数值归零的排查方向 #

• 全场归零是典型的数值发散表现，除边界条件外，需检查：
  • 时间步长：不可压缩流需满足Courant数<1，步长过大会引发数值震荡；
  • 压力松弛因子：松弛因子过大可能导致压力修正过度，引发迭代崩溃；
  • 初始条件：外形外区域需设置合理的初始速度/压力，避免初始时刻大散度触发发散。

内容的提问来源于stack exchange，提问作者Akin83