算流体力学1.pptVIP

高正阳 计算流体力学 CFD 差分解 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 假扩散—人工粘性 一、概述 1、离散误差的一种 2、研究目的 3、假扩散系数 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 二、原因 1、一阶偏导差分采用一阶截差格式 2、流动与网格交叉 3、非常数源项 高正阳 计算流体力学 CFD 三、解决方法 1、提高格式精度 2、考虑流动与网格交叉 3、考虑非常数源项 高正阳 计算流体力学 CFD 三、Quick 格式的几种处理方式 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 格式 对流项向上游传递扰动 对流项向下游传递扰动 稳定性条件 中心差分 一阶迎风差分 0 绝对稳定 二阶迎风差分 0 绝对稳定 三阶迎风差分 Leonard迎风差分 0 绝对稳定 不同格式的稳定性 高正阳 计算流体力学 CFD 格式 对流项向上游传递扰动 对流项向下游传递扰动 稳定性条件 最佳迎风差分 绝对稳定 超级迎风差分 绝对稳定 Quick格式 指数格式 绝对稳定 不同格式的稳定性 高正阳 计算流体力学 CFD §2.3数学描述 一、微分方程的种类 1.原始变量法（非守恒形式） 高正阳 计算流体力学 CFD 非守恒形式 高正阳 计算流体力学 CFD 2.涡量-流函数法 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 二、湍流 1、层流向湍流的转捩 （1）流体动力学稳定性 层流，动力学稳定 扰动放大条件 速度分布的影响 （2）转捩 射流： Re10； 平板：Re91000 管流： Re2000 高正阳 计算流体力学 CFD 二、湍流结构 三维、非稳态、涡 大涡：边界决定 从主流获得能量 尺寸与流场尺寸相当 低频脉动 各向异性 小涡：从大涡获得能量 趋向各向同性 高正阳 计算流体力学 CFD 三、计算方法 1、直接模拟 节点量大，1cm3内需要105个节点，计算量大 数学描述简单 高正阳 计算流体力学 CFD 2、大涡模拟 大涡对时均值影响大，能量大 计算量较大 3、雷诺时均方程法 应力方程法 粘性系数法 高正阳 计算流体力学 CFD 四、雷诺时均方程 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 高正阳 计算流体力学 CFD 雷诺应力： 高正阳 计算流体力学 CFD 多为常数 高正阳 计算流体力学 CFD 五、湍流模型 1.零方程模型 高正阳 计算流体力学 CFD 2.单方程模型 高正阳 计算流体力学 CFD 3.双方程模型 高正阳 计算流体力学 CFD 评价促 C1，C2对计算影响比较大 各系数通用性有限制 适用于高雷诺数流动 发展的K-e模型 高正阳 计算流体力学 CFD K方程的推导 (1) (1)-(2) (2) (3) 高正阳 计算流体力学 CFD (4) 对上式作时均 高正阳 计算流体力学 CFD 4.雷诺应力方程模型 高正阳 计算流体力学 CFD 节点方程的建立 计算区域的离散 节点方程建立方法 差分格式的特性 高正阳 计算流体力学 CFD x y i j i, j i+1, j i-1, j i, j-1 i, j+1 计算区域的离散 高正阳 计算流体力学 CFD x y i j P E W S N e n w s 高正阳 计算流体力学 CFD 节点方程建立方法 直接差分 泰勒级数法 多项式拟合法 控制容积积分法 控制容积平衡法 高正阳 计算流体力学 CFD 差分格式的特性 1、误差 截断误差：差分方程与微分间的误差 舍入误差：计算机字长限制导致误差 离散误差：精确解与差分方程组精确解 间的误差 高正阳 计算流体力学 CFD 2、相容性 时间、空间步长趋近于零时，截差趋于零 3、收敛性 时间、空间步长趋近于零时，离散误差趋于零 4、减小误差的方法 高正阳 计算流体力学 CFD 　 1 2 3 4 5 精确解 0 0.0473 0.135 0.3679 1 二阶格式 0 0.0582 0.1552 0.3944 1 四阶格式 0 0.0505 0.1348 0.3918 1 截差对误差的影响 高正阳 计算流体力学 CFD 网格的影响 加密前 0.0582 0.1552 0.3944 加密后 0.0502 0.1404 0.3752 精确解 0.0473 0.135 0.3679