MRF的发展

最新推荐文章于 2021-02-26 22:36:16 发布

寒冰彻骨

最新推荐文章于 2021-02-26 22:36:16 发布

阅读量964

点赞数 1

分类专栏：控制方程文章标签： CFD OpenFOAM

原文链接：http://openfoamwiki.net/index.php/See_the_MRF_development

版权

控制方程专栏收录该内容

13 篇文章 11 订阅

订阅专栏

本节描述了旋转框架中不可压缩Navier-Stokes公式的发展。

文章目录

1.加速度

首先，我们来看看旋转框架的加速度项 $\vec \Omega:$
$\begin{pmatrix} 0 & -\omega_z & -\omega_y\\ -\omega_z & 0 & -\omega_x \\ -\omega_y & -\omega_x & 0 \\ \end{pmatrix}$

符号：I：惯性，R：旋转

对于一般向量：
$\left [ \frac{d \vec A}{dt} \right ]_I = \left [ \frac{d \vec A}{dt} \right ]_R + \vec \Omega \times \vec A$
对于位置向量:
$\left [ \frac{d \vec r}{dt} \right ]_I = \left [ \frac{d \vec r}{dt} \right ]_R + \vec \Omega \times \vec r$
$\vec u_I = \vec u_R + \vec \Omega \times \vec r$
加速度可以表示为：
$\left [ \frac{d \vec u_I}{dt} \right ]_I = \left [ \frac{d \vec u_I}{dt} \right ]_R + \vec \Omega \times \vec u_I$

$\left [ \frac{d \vec u_I}{dt} \right ]_I = \left [ \frac{d \left [ \vec u_R + \vec\Omega \times \vec r \right ] }{dt} \right ]_R + \vec \Omega \times \left [ \vec u_R + \vec \Omega \times \vec r \right ]$

$\left [ \frac{d \vec u_I}{dt} \right ]_I = \left [ \frac{d \vec u_R}{dt} \right ]_R + \frac{d \vec \Omega}{dt} \times \vec r + \vec \Omega \times \underbrace{ \left [ \frac{d \vec r}{dt} \right ]_R }_{\vec u_R} + \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r$

$\left [ \frac{d \vec u_I}{dt} \right ]_I = \left [ \frac{d \vec u_R}{dt} \right ]_R + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r \tag1$

2. 惯性坐标系系下Navier-Stokes方程的绝对速度

惯性系中分子粘度恒定的不可压缩Navier-Stokes方程为：

$\begin{cases} \frac {D \vec u_I}{D t} = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ \nabla \cdot \vec u_I = 0 \end{cases} \tag2$

$\begin{cases} \frac {\partial \vec u_I}{\partial t} + \vec u_I \cdot \nabla \vec u_I = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ \nabla \cdot \vec u_I = 0 \end{cases}$

$\begin{cases} \frac {\partial \vec u_I}{\partial t} + \nabla \cdot (\vec u_I \otimes \vec u_I) - \underbrace{( \nabla \cdot \vec u_I )}_{0} \vec u_I= - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ \nabla \cdot \vec u_I = 0 \end{cases}$

$\begin{cases} \frac {\partial \vec u_I}{\partial t} + \nabla \cdot (\vec u_I \otimes \vec u_I) = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ \nabla \cdot \vec u_I = 0 \end{cases} \tag3$

3.相对坐标系下Navier-Stokes方程的相对速度

让我们看看方(2)的动量方程的左边，通过考虑方程(1)的加速度项：
$\frac {D \vec u_I}{D t} = \frac{D \vec u_R}{Dt} + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r$

$\frac {D \vec u_I}{D t} = \frac {\partial \vec u_R}{\partial t} + \vec u_R \cdot \nabla \vec u_R + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r$

$\frac {D \vec u_I}{D t} = \frac {\partial \vec u_R}{\partial t} + \nabla \cdot (\vec u_R \otimes \vec u_R) + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r \tag4$
因为 $\nabla \cdot \vec u_R = \nabla \cdot \vec u_I = 0$ ,
$\begin{align}\nabla \cdot \vec u_I & = \nabla \cdot \left [ \vec u_R + \vec \Omega \times \vec r \right ] = 0 \& = \nabla \cdot \vec u_R + \underbrace{\nabla \cdot \left [ \vec \Omega \times \vec r \right ]}_{0} = 0 \& = \nabla \cdot \vec u_R = 0\end{align}$$
另外，可以注意到：
$\begin{align}\nabla \cdot \nabla (\vec u_I) & = \nabla \cdot \nabla \left [ \vec u_R + \vec \Omega \times \vec r \right ] \& = \nabla \cdot \nabla (\vec u_R) + \nabla \cdot \underbrace{\nabla (\Omega \times \vec r)}_{0} \& = \nabla \cdot \nabla (\vec u_R )\end{align}$

方程(3)可以写成:

$\begin{cases} \frac {\partial \vec u_R}{\partial t} + \frac{d \vec \Omega}{dt} \times \vec r + \nabla \cdot (\vec u_R \otimes \vec u_R) + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_R) \\ \nabla \cdot \vec u_R = 0 \end{cases} \tag5$

4.相对坐标系下Navier-Stokes方程的绝对速度

方程(5)用旋转速度表示旋转框架中不可压缩的Navier-Stokes方程。

其左侧第三项可以描述为:
$\nabla \cdot (\vec u_R \otimes \vec u_R) & = \nabla \cdot ( \vec u_R \otimes \left [ \vec u_I - \vec \Omega \times \vec r \right ] ) \& = \nabla \cdot (\vec u_R \otimes \vec u_I) - \underbrace{\nabla \cdot \vec u_R}_{0} (\vec \Omega \times \vec r) - \underbrace{\vec u_R \cdot \nabla(\vec \Omega \times \vec r)}_{\vec \Omega \times \vec u_R} \& = \nabla \cdot (\vec u_R \otimes \vec u_I) - \vec \Omega \times \vec u_R\end{alignat}$
因此，式(5)左边的稳定项可以写成:
在这里插入图片描述
方程(5)可以用绝对速度表示：

$\begin{cases} \frac {\partial \vec u_R}{\partial t} + \frac{d \vec \Omega}{dt} \times \vec r + \nabla \cdot (\vec u_R \otimes \vec u_I) + \vec \Omega \times \vec u_I = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ \nabla \cdot \vec u_I = 0 \end{cases} \tag6$