关于动网格的求解方法及步骤咨询

[winsway_zero]

想咨询一下动网格的控制方程应该如何处理呀；控制方程里面忽略了X,Y方向的流动(U,V)，主要关注Z方向(W)的流动问题。
流动只考虑z方向，但是粘性项需要考虑x和y方向。由于网格是随着之间变化的，因此控制方程应该需要进行修改，这个部分理论不足，真不知道如何处理？希望能够有这方面的行家交流学习一下，谢谢谢谢~

[李东岳]

没太看明白你的意思，最近在准备OKSS2的求解器，你在icoFoam里面，添加两个信息：

1. mesh.update()函数处理网格变形

2. 通过phi = phi - mesh.phi()将通量变为相对通量代入到动量方程，压力方程不变

就可以了

从算法来讲，静态网格：
$$\begin{array}{}\text{(1)}& \frac{\mathrm{\partial }\rho }{\mathrm{\partial }t}+\mathrm{\nabla }\cdot \rho \mathbf{U}=0,\end{array}$$
$$\begin{array}{}\text{(2)}& \frac{\mathrm{\partial }\mathbf{U}}{\mathrm{\partial }t}+\mathrm{\nabla }\cdot (\mathbf{U}\otimes \mathbf{U})=-\mathrm{\nabla }\frac{p}{\rho }+\mathrm{\nabla }\cdot (\nu \mathrm{\nabla }\mathbf{U}),\end{array}$$
对于动态网格，需要考虑相对速度通量，则变为：
$$\begin{array}{}\text{(3)}& \frac{\mathrm{\partial }\rho }{\mathrm{\partial }t}+\mathrm{\nabla }\cdot \rho (\mathbf{U}-{\mathbf{U}}_b)=0,\end{array}$$
$$\begin{array}{}\text{(4)}& \frac{\mathrm{\partial }\mathbf{U}}{\mathrm{\partial }t}+\mathrm{\nabla }\cdot (\mathbf{U}(\mathbf{U}-{\mathbf{U}}_b))=-\mathrm{\nabla }\frac{p}{\rho }+\mathrm{\nabla }\cdot (\nu \mathrm{\nabla }\mathbf{U}),\end{array}$$
同时附加空间守恒法则：
$$\begin{array}{}\text{(5)}& \frac{\mathrm{\partial }\rho }{\mathrm{\partial }t}-\mathrm{\nabla }\cdot \rho {\mathbf{U}}_b=0,\end{array}$$
把$\text{(5)}$代入到$\text{(3)}$变成了：
$$\begin{array}{}\text{(6)}& \mathrm{\nabla }\cdot \rho \mathbf{U}=0,\end{array}$$
这样，$\text{(6)}$用来构建压力泊松方程，$\text{(4)}$用来构建动量方程，求解即可

[winsway_zero]

@东岳 谢谢东岳老师，我明白如何处理了哈。东岳老师整理的文档和资料都很通俗易懂。