【CFD理论】对流项-01

基本对流格式

练习：
Divergence scheme example

Divergence schemes

$\nabla \cdot Q$

• 如果Q是矢量，等式可以写成：
  $\nabla \cdot Q=\frac{\partial Q_x}{\partial x}+\frac{\partial Q_y}{\partial y}+\frac{\partial Q_z}{\partial z}$

divSchemes
{
    default         none;
    div(Q)          Gauss <interpolation scheme>;
}

• 对于对流项的典型使用，传输一个属性Q，在速度场的影响 $\varphi$，特别的定义：

divSchemes
{
    default         none;
    div(phi,Q)      Gauss <interpolation scheme>;
}

$\varphi$ 这个关键词典型的用于代表通量（流）穿过cell面，例如

• 体积通量：$\varphi ={\mathbf{u}}_f\cdot {\mathbf{S}}_f$
• 质量通量：$\varphi ={\rho }_f({\mathbf{u}}_f\cdot {\mathbf{S}}_f)$

1. Upwind divergence scheme

• 一阶
• 有界
• 面心值是根据上游值来确定
  ${\varphi }_f={\varphi }_c$

divSchemes
{
    default         none;
    div(phi,U)      Gauss upwind;
}

2. Linear divergence scheme

• 二阶
• 无界
• 对于LES计算是很好的选择，因为低耗散
  ${\varphi }_f=0.5({\varphi }_c+{\varphi }_d)$

divSchemes
{
    default         none;
    div(phi,U)      Gauss linear;
}

3. Linear-upwind divergence scheme

• 迎风插值权重，用一个显式的修正基于local cell gradient
• 二阶精度
• 无界

divSchemes
{
    default         none;
    div(phi,U)      Gauss linearUpwind grad(U);
}

4. MUSCL divergence scheme

• 守恒定律的单调上游中心方案
• 二阶
• 无界

divSchemes
{
    default         none;
    div(phi,U)      Gauss MUSCL;
}

5. Limited linear divergence scheme

• 一/二阶
• 无界

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinear <coeff>;
}

其他参考网址：

• 6.2 Numerical schemes
• 3.16 Gradient limiting
• book-2-An Introduction to Computational Fluid Dynamics, The Finite Volume Method (chapter 5)

一维定常对流扩散方程

解析解

方程：
$\frac{d(\rho u\varphi )}{dx}=\frac{d}{x}(\Gamma \frac{d\varphi }{dx})$

解析解：
$\rho u\varphi -(\Gamma \frac{d\varphi }{dx})=c$
$\frac{d\varphi }{dx}=\frac{\rho u}{\Gamma }\varphi -\frac{c}{\Gamma }$
$where\varphi =\frac{\rho u}{\Gamma }\varphi -\frac{c}{\Gamma }$
$\frac{d\varphi }{\varphi }=\frac{\rho u}{\Gamma }dx\Rightarrow \frac{\rho u}{\Gamma }x+c_3\Rightarrow \varphi =\frac{C\lambda e^{\frac{\rho u}{\lambda }x}+c}{\rho u}$

$\frac{\varphi -{\varphi }_W}{{\varphi }_E-{\varphi }_W}=\frac{e^{P{e}_L\frac{x-x_W}{L}}-1}{e^{P{e}_L}-1}$

佩克莱数 Peclet number,简称Pe数:
$P{e}_L=\frac{\rho uL}{\Gamma }(对流/扩散),L=x_E-x_W$

• 当$P{e}_L=0$时，是纯扩散型
• 一般$P{e}_L\le 2$

数值解

三维 对流-扩散
${\int }_{V_P}[\nabla \cdot (\rho \mathbf{U}\varphi )-\nabla \cdot (\Gamma \nabla \varphi )]dV=0$
${\int }_{V_P}[\nabla \cdot (\rho \mathbf{U}\varphi )-\nabla \cdot (\Gamma \nabla \varphi )]dV={\int }_{\partial V_P}(\rho \mathbf{U}\varphi -\Gamma \nabla \varphi )\cdot d\mathbf{S}={\int }_{\partial V_P}[\rho \mathbf{U}\varphi \mathbf{i}-\Gamma \frac{d\varphi }{dx}\mathbf{i}]\cdot d\mathbf{S}=0$

• CD
• upwind
• second order upwind

Central Difference (CD) Scheme

中心差分格式

• 公式：
  $\varphi (x)=k_0+k_1(x-x_C)$
• e处的插值格式：
  ${\varphi }_e={\varphi }_C+\frac{({\varphi }_E-{\varphi }_C)}{(x_E-x_C)}(x_e-x_C)$
• 均匀网格：
  ${\varphi }_e=\frac{{\varphi }_C+{\varphi }_E}{2}$

数值解
对流-扩散e处插值方程
$(\rho u\Delta y{)}_e-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_e=(\rho u\Delta y{)}_e\frac{{\varphi }_C+{\varphi }_E}{2}-(\Gamma \frac{\Delta y}{\delta x}{)}_e({\varphi }_E-{\varphi }_C)$
对流-扩散w处插值方程
$-[(\rho u\Delta y{)}_w-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_w]=-[(\rho u\Delta y{)}_w\frac{{\varphi }_C+{\varphi }_E}{2}-(\Gamma \frac{\Delta y}{\delta x}{)}_w({\varphi }_C-{\varphi }_W)]$
$\because [(\rho u\Delta y{)}_e-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_e]-[(\rho u\Delta y{)}_w-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_w]=0$
$$\begin{gathered} \therefore [(\rho u\Delta y{)}_e-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_e]=[(\rho u\Delta y{)}_w-(\Gamma \frac{d\varphi }{dx}\Delta y{)}_w] \\ (\frac{(\rho u\Delta y{)}_e}{2}+{\Gamma }_e\frac{\Delta y_e}{\delta x_e})+[(\frac{(\rho u\Delta y{)}_w}{2}+{\Gamma }_w\frac{\Delta y_w}{\delta x_w})]{\varphi }_C+[\frac{(\rho u\Delta y{)}_e}{2}-{\Gamma }_e\frac{\Delta y_e}{\delta x_e}]{\varphi }_E+[\frac{(\rho u\Delta y{)}_w}{2}-{\Gamma }_w\frac{\Delta y_w}{\delta x_w}]{\varphi }_W=0 \\ \frac{(\rho u\Delta y{)}_e}{2}-{\Gamma }_e\frac{\Delta y_e}{\delta x_e}>0 \\ \Rightarrow Pe=\frac{(\rho u{)}_e\delta x_e}{{\Gamma }_e}>2 \end{gathered}$$

• 中心差分格式 $\to$容易振荡（发散）
• 为了减少邻值误差放大，通过减小网格尺寸，一般设置为$Pe\le 2$。

example 5.2

• table 5.2是绝对对角占优
• 纯扩散方程求解是对称矩阵，为什么对流-扩散方程不是对称矩阵？
• $a_W=D+\frac{F}{2},a_E=D-\frac{F}{2}$