11-计算流体力学理论基础.ppt

2.4流体力学方程组其他形式2.4流体力学方程组其他形式2.4流体力学方程组其他形式2.4流体力学方程组其他形式2.4流体力学方程组其他形式2.4流体力学方程组其他形式2.4流体力学方程组其他形式1.流场（流体运动过程与规律）2.方程组（反映流体介质运动过程中所要遵循的规律）3.控制体与初始条件与边界条件（决定流体运动过程与规律）2.5流体力学方程组数学性质2.5流体力学方程组数学性质2.4.2N-S方程性质对于定常N-S方程，方程为非线性椭圆方程对于非定常N-S方程，方程性质为双曲-抛物性方程2.4.3EULER方程性质对于亚音速流动，定常欧拉方程为椭圆型方程对于超音速流动，定常欧拉方程为双曲型方程对于跨音速流动，定常欧拉方程为椭圆与双曲混合型方程更多形式参考文献[计算流体力学—差分方法的原理与应用；第一章：流体力学的各级近似方程及其数学性质；张涵信，沈孟育著；国防工业出版社]2.6模型方程2.6模型方程2.6模型方程2.6模型方程2.6模型方程2.6模型方程2.6模型方程2.6模型方程2.6模型方程*/573.2.5差分格式的稳定性分析方法4.波动方程差分方程迎风格式把代入格式恒不稳定格式有条件稳定*/573.3差分方程数值特点3.3.1差分方程模型方程3.3.2数值耗散与色散3.4常见的差分格式*/573.3.1差分方程模型方程定义：差分方程所准确逼近的微分方程（即截断误差为零）修正方程是差分方程的微分表达式;修正方程与原方程之差是截断误差;对差分方程稳定性分析有帮助,截断误差中，偶次导数项称为耗散误差项，奇次导数项称为色散误差项*/572、推导过程（2）其有如下差分方程：（3）将和按?Tailer级数展开代入上式有：（4）计算修正方程中对高阶导数项（5）消去修正方程中高阶导数项，得只包含对自变量的各阶导3.3.1差分方程模型方程（1）考虑波动方程：*/573、修正方程性质分析：修正方程是差分方程的微分表达式，对差分方程稳定性分析有帮助修正方程与原方程之差是截断误差截断误差中，偶次导数项称为耗散误差项，奇次导数项称为色散误差项截断误差中，最低阶为偶次，则耗散误差占主导，最低阶为奇次，则色散误差占主导。截断误差中偶次导数项又称为差分粘性项，或者隐式人工粘性项，相对于直接加入到差分方程中的人工粘性项称为显式人工粘性项。3.3.1差分方程模型方程*/573.3.2数值耗散与色散1、修正方程通式：2、误差分析：满足修正方程，则则方程解析解为：设误差分布函数*/57修正方程中，右边偶数项影响振幅，称为数值解的耗散效应修正方程中，右边奇数项影响频率与相位，称为数值解的色散效应则误差增长的增幅随时间增大而减小，格式稳定则误差增长的增幅随时间增大而增大，格式不稳定则误差不减少，频率与相位出现变化，导致数值结果出现锯齿型高频振动，中性稳定，计算也会失败。4、与稳定性的关系5、实用的稳定性判据3.3.2数值耗散与色散3、方程性质分析*/573.4常见的差分格式1、模型方程离散格式稳定性条件（1）显式一阶迎风格式稳定性条件：格式精度：*/573.4常见的差分格式（2）Lax格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（3）Lax-Wendroff格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（4）LeapFrog格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（5）两步Lax-Wendroff格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（6）Mac-Cormack格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（7）Warming-Kutler-Lomax格式稳定性条件：格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（8）时间中心隐式格式稳定性条件：无条件稳定格式精度：1、模型方程离散格式稳定性条件*/573.4常见的差分格式（1）简单显式差分稳定性条件：格式精度：1、