第二章 控制方程离散化之低阶格式（一）

好久没有更新了，每次敲完公式，都要自己添加\nonumber才不会显示标号，对CSDn也是服了。
最近也不知道忙什么，课也不多，琐事倒是不少，正所谓春水碧于天，划船听雨眠，好了，进入正题
目标

通过本节，了解有限体积法基本原理，通过中心差分、一阶迎风、混合迎风了解这些格式的优缺点

2.1离散化概述

在对CFD计算之前，需要对计算区域离散化，也就是说要对空间连续的计算区域进行划分，把它划分为许多离散子区域，通过确定每个区域中的节点，从而生成网格，说人话就是把偏微分方程转化为代数方程组进行求解，对于瞬态，还要对时间进行离散化，那我们先从空间离散进行讨论吧

2.1.1离散化目的

目的我就不说了，还不是因为偏微分方程那魔性，至今还没有办法求出真解，但是，还好，计算机最不怕的就是迭代，我们可以通过逼近，得到最接近真解的数值解。

2.1.2离散化使用的网格

网格的质量对计算结果影响非常大，至于原因，参考积分怎么求的吧，在二维时，常见的网格有三角形，四边形单元，在三维问题，有四面体、六面体，棱锥， 混合网格等单元

2.1.3 常见离散化方法

1.有限差分法

这个方法在教程中经常提到，通过将求解区域划分为差分网格，用有限个网格节点代替连续的求解域，用差商代替偏导数，这个方法可以用于双曲型和抛物型问题，对于复杂的问题就呵呵了

2.有限元法

有限元法在固体力学占据着主导地位，原因呢，知乎找找就有了，由于其在求解CFD时较慢，在CFD领域用的并不广泛，有限元就是将一个连续的求解区域任意划分为任意形状的许多微小单元，对微小单元构成的方程组进行求解。

3.有限体积法

有限体积法（Finite Volume Method）在流体力学领域用的可就非常广泛了，下面我们开始介绍它

2.2有限体积法及其网格简介

2.2.1 基本思想

有限体积法又叫控制体积法，基础思路为：把计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；将待求解微分方程对每一个控制体积积分，从而得到一组离散的方程。其中的未知数是网格点上的因变量$\phi$。为了求出控制体积的积分，必须假设$\phi$值在网格点之间的变化规律。从积分区域的选取方法看来，有限体积属于加权余量法中的子域法，从未知解的近似方法看来，有限体积法属于局部近似的离散方法。简言之，子域法加离散，就是有限体积的基本方法。

有限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。这就是其优点，对于其他的离散方法，如差分，必须网格及其密时，离散方程才满足积分守恒，而有限体积法即使在粗网格情况下，也显示出准确的积分守恒

2.3一维稳态问题的有限体积法

控制方程

$$\frac{（\text{d}\rho u \phi）}{\text{d}x} =\frac{\text{d}}{\text{d}x}\big(\Gamma \frac{\text{d}\phi}{\text{d}x}\big)+S$$
该方程又称为一维模型方程，方程中包含对流项、扩散项及源项，方程中的 $\phi​$是广义变量，可以为速度。温度或浓度等一些待求的物理量， $\Gamma​$是相应于 $\phi​$的广义扩散系数， $S​$是广义源项，变量 $\phi​$在端点A和B的边界值为已知。

第一步：生成计算网格$\delta$

建立离散方程，对控制体积P进行积分，

$$\int_{\Delta V}\frac{\text{d}\rho u \phi}{\text{d}x}\text{d}V =\int_{\Delta V}\frac{\text{d}}{\text{d}x}\big(\Gamma \frac{\text{d}\phi}{\text{d}x}\big)\text{d}V+\int_{\Delta V}S\text{d}V$$
式中， $\Delta V$是控制体积的体积值，当控制体积很微小时， $\Delta V$可以表示为 $\Delta V+A$，这里 $A$是控制体积界面的体积

这样积分就可以进行计算了，这个简化过程很重要

$$(\rho u \phi A)_e - (\rho u \phi A)_w =\big(\Gamma A \frac{\text{d}\phi}{\text{d}x}\big)_e - \big(\Gamma A \frac{\text{d}\phi}{\text{d}x}\big)_w +S\Delta V \tag{1}$$
目前我们的位置数为界面 $w$和 $e$处的 $\rho$、 $u$、 $\Gamma$、 $\phi$ $\frac{\text{d}\phi}{\text{d}x}$。为了计算界面出的未知物理量，需要知道物理参数在节点间的近似分布。简单的说，如果是线性分布的话，我们采用的方法就叫做 中心差分，对于均匀网格，各个物理量可以用下面的公式进行求解
$$\begin{align} \Gamma_e&=\frac{\Gamma_P+\Gamma_E}{2}\nonumber\\ \Gamma_w&=\frac{\Gamma_W+\Gamma_P}{2}\nonumber\\ \phi_e &= \frac{\phi_P+\phi_E}{2}\nonumber\\ \phi_w &= \frac{\phi_W+\phi_P}{2}\nonumber\\ \Big(\frac{\text{d}\phi}{\text{d}x}\Big)_e&=\frac{\phi_E-\phi_P}{(\delta x)_e}\nonumber\\ \Big(\frac{\text{d}\phi}{\text{d}x}\Big)_w&=\frac{\phi_P-\phi_W}{(\delta x)_e}\nonumber \end{align}$$
$S​$通常为时间和 $\phi​$的函数，为了简化，可以做如下线性处理：