我对 CFD 很陌生，我想从 Chorin 的投影开始。起始问题只是在研究区域中流动的自由射流。我几乎立刻就得到了可怕的振荡，我认为这是由于对中间速度的边界条件处理不当$u^*$. 在这个过程中的某个地方只是一个菜鸟错误：

1. 的计算$u^*$使用：

$$u^* = u^{n}+\Delta t(u^n\cdot\nabla)u^n$$

2. 用源项求解泊松方程$\frac{\Delta t}{\rho}\nabla \cdot u^*$然后更新速度：

$$u^{n+1} = u^* + \frac{\Delta t}{\rho}\nabla p^{n+1}$$

• BC 应该在哪里进行$u^*$介入;涉足？
• 如何使用边界条件$u^*=u^{n+1}$如果我没有$u^{n+1}$然而？

我需要零 Neumann BC 的速度和零 Dirichlet BC 的压力。