计算科学 - 仅在外力作用下的空腔流动。为什么会流传？ - 吾爱随笔录

使用来自http://math.mit.edu/~gs/cse/codes/mit18086_navierstokes.pdf的代码，我模拟了具有外力的空腔流动，所有 4 侧的边界条件为和。的步骤中向 u 方向添加外力到。所以代码只改变了两行，从 $u=v=0$ $\frac{\partial p}{\partial n}=0$

(1) \frac{U^{*} - U^{n}}{Δ t} = - ((U^{n})^{2})_{x} - (U^{n} V^{n})_{y}

$\text{(1)} \quad \frac{U^*-U^n}{\Delta t} = -((U^n)^2)_x - (U^nV^n)_y$

(1) \frac{U^{*} - U^{n}}{Δ t} = - ((U^{n})^{2})_{x} - (U^{n} V^{n})_{y} + f

$\text{(1)} \quad \frac{U^*-U^n}{\Delta t} = -((U^n)^2)_x - (U^nV^n)_y +f$

uN = x*0+1
U = U-dt*(UVy(2:end-1,:)+U2x);

到

uN = x*0
U = U-dt*(UVy(2:end-1,:)+U2x+10);

，我得到了这个结果。

我对这个结果感到困惑。这种外力就像重力，所以我想，当只施加重力时，腔内的流体不会循环。它就像玻璃杯里的水。为什么会流传？这是物理上正确的还是因为数值模拟方法？