Rhie-Chow插值有没有对应的连续的微分方程/积分方程?
-
越来越复杂了。
可否认为Rhie-Chow是用额外的3阶项(对压力方程式4阶项)来平滑掉压力波动的,因为是3阶项,所以不会影响FVM的二阶及以下的空间精度。但是这样的话R-C插值算激波绝对是有问题的,激波的FFT展开后频率很宽,微分高阶项比低阶项大,这玩意儿搞插值还是会抖啊。
关于黑体字书上确实这么解释的。同为网格压力棋格子分布很明显就可以看出来,考虑一维度情况并且方程左边不做处理:
\begin{equation}
f(u)_P=\frac{dp}{dx}=p_e-p_w=\frac{p_E+p_P}{2}-\frac{p_P+p_W}{2}=\frac{p_E-p_W}{2}
\end{equation}
可见关于$u_P$的离散后的方程不包含$p_P$,但是却和$p_W,p_E$联系起来了。这就是棋格子分布的产生。不管什么插值,都是要在$u_P$的离散方程添加$p_P$的影响。不太明白,RC插值和那个A,H算子是个什么关系?A, H 算子是simple/piso算法这个层次里用到的东西还是RC差值实现的一个技巧?
UEqn.A()
和UEqn.H()
是和SIMPLE以及SIMPLEC有关的,考虑SIMPLE,UEqn.A()
和UEqn.H()
就是离散后压力方程的对角线系数以及临点的影响。可能说的不是很精准但是大体是这个意思。SIMPLEC算法需要重写UEqn.H()
。这部分太复杂了并且感兴趣的人很少。 -
-
不可压缩流既可以算压力修正方程(pressure correction equation),e.g. SIMPLE - semi-implicit method with pressure-link equation,也可以算压力泊松方程(pressure Poison equation), e.g., PISO - pressure implicit with splitting of operator.
-
可压缩流既可以算压力修正方程(pressure correction equation),e.g. SIMPLE - semi-implicit method with pressure-link equation,也可以算压力泊松方程(pressure Poison equation), e.g., PISO - pressure implicit with splitting of operator.
以上两点都是隐式压力更新,并且压力矩阵的残差在很大程度上控制了质量守恒。算SIMPLE跟PISO只是两个不同的保证质量守恒的方式而已。
那个临点的说法,把很简单的问题复杂化了,已经脱离了Rhie-Chow插值的话题。
-