计算流体力学new.pptVIP

计算流体力学 陈丽萍 第1章 绪论 §1-2有限体积法的基本思想 一.通用变量方程 压力速度耦合方程 第1章 绪论 二.有限体积法的基本思想 1.求解域离散,形成网格.每一网格节点按一定的方式形成一个包围该节点的控制容积V。 2.有限体积法的关键步骤是将控制微分方程式在控制容积内积分,即 第1章 绪论 利用高斯散度定理将对流项和扩散项的体积分转换成关于控制容积V表面A上的面积分。 高斯散度定理 三.有限体积法的特点 （1）积分形式的控制方程，确保守恒特性。 （2）积分方程中每一项都有明确的物理意义。 （3）区域离散的节点网格与进行积分的控制容积分立。 第2章 区域离散和导数的差分格式 § 2-1区域离散的实质与内容 实质：用有限个离散的点来代替原来连续空间。 内容： （1）计算区域划分成许多个互不重迭的子区域 （2）确定每个子区域中的节点位置及该节点所代表的控制容积 （3）确定分界面位置 第2章 区域离散和导数的差分格式 § 2-2两类设置节点的方法 一.外节点法 先确定节点的坐标,再确定界面(两节点间中心位置) 第2章 区域离散和导数的差分格式 二.内节点法 先确定界面位置,再确定节点的坐标(两界面中心位置) 第2章 区域离散和导数的差分格式 三.两种方法比较 第2章 区域离散和导数的差分格式 三.两种方法比较 区域离散：网格参数 一维为例 区域离散：网格参数 第2章 区域离散和导数的差分格式 § 2-3导数的差分格式 一.一阶导数 第2章 区域离散和导数的差分格式 一.二阶导数 第3章 扩散问题的有限体积法 § 3-1一维稳态扩散 § 3-1一维稳态扩散 1.区域离散(内节点法） 3.1 3.1 3.1 § 3-2多维稳态扩散 例3.3 二维受热平板。板厚1cm，导热系数k=1000w/(mk),西侧边界有稳定热流输入，热流强度q=500kw/m2.东侧和南侧边界绝热，北侧边界保持常温，TN=100oC. 求板内温度分布。 第4章 对流扩散问题 例4.1 设某场变量Φ经过对流扩散过程从一维区域的x=0点输运到x=L点，，流体密度ρ=1.0kg/m3,L=1.0m,扩散系数Г=0.1kg/(ms),求分布 1.内节点法，离散成5个节点网格，δx=0.2m 2.方程离散 3.边界节点离散 （1）节点1： w界面上 第5章 差分格式问题 此章讨论对流项离散时界面φ如何用节点φ近似。近似过程中应满足3个特性：守恒性、有界性、输运性。 可证明中心差分满足守恒性 § 5.2一阶差分格式 迎风格式的主要缺点是当流动不是沿着网格线时结果有误差.此时,迎风格式使输运量Φ的分布变得混淆,这样产生的误差有扩散特性,称为伪扩散(false diffusion) 理论上网格加密可解决伪扩散问题,实际中网格加密代价太大,不可行.如此,迎风差分格式不适合精确流动计算. 混合差分格式: 小Pe数(Pe2)用中心差分;大Pe数(Pe=2)用迎风格式.混合差分格式的主要缺点是只有泰勒级数截断误差的一级精度. 指数格式精度高，但费时。 乘方格式精度高。 4.3对流扩散问题的高阶差分格式 二阶迎风QUICK格式 QUICK(The Quadratic Upstream Interpolation for Convective Kinetics )格式 用单元面两侧的两个节点和上游侧的再一个节点的二次插值表示这个单元面的输运量Φ 值. 可证明对均匀网格,在单元面两侧的节点i和i-1,上游节点i-2之间插值,得 对正负流动方向都有效的一般表达式: QUICK格式的不足之处是不能保证主系数(E和W)为正,而且aww,aEE为负,这就有稳定性问题, QUICK格式是有条件稳定的.改进的QUICK格式把带来稳定性问题的系数放进源项,这样主系数就是正的. QUICK格式精度比中心差分和混合格式高,保持迎风加权的特性,总的伪扩散性小,但计算较不稳定.但某些问题如Φ梯度较大时, QUICK格式会出现小的过低和过冲,极端情况会出现非物理行为,如出现负湍流量k,ε. 尽管如此,只要使用时小心判断, QUICK格式能给出非常精确的解. 工况2：u=2.5m/s代入， 解得： 误差： 网格加密误差减少。 ?x (?x)+e (?x)-w (?x)-e (?x)+w (?x)w (?x)e w e W P E 图 1 一维问题空间区域的离散化 一.守恒性 西界面通量 东界面通量 近似过程 守恒性表现为:无源项时, 西界面通量= 东界面通量 二.有界性 扩散问题和对流扩散问题离散后统一形式如下: 扩散问题中 对流扩散问题中 美国学者Scarborough(1958年）提出了由离散方程系数判断上述方程求解时是否收敛于