计算流体力学笔记D3

标准的有限差分方法需要均匀网格，非均匀网格要变换为均匀分布的矩形网格。

物理平面与计算平面的变换

链式法则

 

 

 

 关于的导数称为度量，

（1）用封闭解析表达式给出；（2）有限差分，通常是中心差分。

 直接变换，直接度量；逆变换，逆度量

网格拉伸

贴体型网格： 通过求解椭圆型偏微分方程（如Laplace方程），变换到计算网格，称为椭圆型网格生成。这里的椭圆型方程知识用来生成网格的，与流动本身的性质无关。

自适应网格，在流场控制方程的时间推进求解中，网格的调整与按时间步计算流场变量的过程同步。

结构网格：有限差分，均匀与非均匀网格。非结构网格：有限体积法，多用于有限元。

• 显式计算方法

Lax-Wendroff拉克斯-温德罗夫方法，显式有限差分

MacCormack麦考马克方法，Lax-Wendroff的变种，显式有限差分。预估步向前差分，校正步向后差分。

均适用于非定常流NS方程求解，具有抛物型与椭圆型混合的性质，不适合于定常NS，因为是椭圆型偏微分方程。根据经验，对流项的预估步、校正步，空间导数分别用向前差分、向后差分表示，粘性项的预估步和校正步都应使用中心差分。

对于显式求解，沿流向的推进一定要用守恒型方程。

松弛法，适用于求解椭圆型偏微分方程的有限差分方法，可以是显式或隐式，显式松弛法又叫点迭代法。简单迭代法、Gauss-Seidel迭代法、逐次松弛法。

• 数值耗散

差分方程的精确解是原偏微分方程的数值解，解带有离散误差；另一种观点，差分方程的精确解是另一个偏微分方程（修正方程）的精确解，无离散误差。

数值耗散和色散来源于数值过程，没有物理意义，数值耗散项的系数称为人工粘性。数值耗散是修正方程右端项中偶数阶导数的直接结果，数值色散则是奇数阶导数的直接结果。数值耗散的影响会将波抹平，数值色散的影响表现为波前波后出现振荡。

人工粘性降低了解的精度，但通常提高解的稳定性。人工粘性的作用使数值解稳定、光滑。一种新的算法，能够自动的在需要的地方添加合适的人工粘性，如TVD全变差减小。

• 隐式计算方法

ADI（Alternating Direction Implict）交替方向隐式方法，对求解抛物线型偏微分方程尤其合适。

第一步，用 t 计算差分方程的x方向是隐式的，第二步，计算，差分方程的y方向是隐式的。

ADI代表一类格式，除了常见的，还有一种叫做近似因子分解方法。

压力修正法SIMPLE，压力耦合方程的半隐式算法

拉普拉斯算子与散度的关系

假设=常数，T=常数，因为是T的函数，所以=常数，能量方程从方程组中解耦，表明不可压流动的速度场与压力场只需用到连续性方程和动量方程。 

 将中心差分用于不可压NS方程时，一旦出现棋盘式的离散压力分布，差分方程就无法改变，改进方法：迎风差分；交错网格使用中心差分。

交错网格使用中心差分

压力修正法，时间导数向前差分，空间导数中心差分，交错网格。