我正在求解 2D 中的稳态不可压缩斯托克斯方程：

$$\begin{equation} \frac{\partial u_x}{\partial x} + \frac{\partial u_y}{\partial y} = 0, \end{equation}$$

$$\begin{equation} \mu\left[\frac{\partial^2 u_x}{\partial x^2} + \frac{\partial^2 u_x}{\partial y^2}\right] - \frac{\partial p}{\partial x} = 0, \end{equation}$$ 和 $$\begin{equation} \mu\left[\frac{\partial^2 u_y}{\partial x^2} + \frac{\partial^2 u_y}{\partial y^2}\right] - \frac{\partial p}{\partial y} = 0. \end{equation}$$

我正在使用具有交错网格方法的有限差分（在面定义的速度和在单元中心的压力）。

我正在努力寻找如何在左面和底面实现压力边界条件。我的压力模板是“L”模板：我使用$p_{i,j+1}$,$p_{ij}$， 和$p_{i+1,j}$. 因此，在左侧和左下角的体积中，我不使用“幽灵单元”中的压力值。有关如何在此设置中实施压力 BC 的任何线索？例如，我想通过平行板模拟（左右面的规定压力，顶面和底面的无滑移）来验证我的代码。