计算流体力学(中科院力学所)_第7讲-差分方法3.pptVIP

* Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第七讲 差分方法（3） 李新亮 lixl@ ；力学所主楼219； 知识回顾： 守恒型方法 通量分裂—— 逐点分裂与特征投影分裂 知识点： 精确特征方向的守恒型格式——Roe格式 数值解的震荡机制及对策—— TVD格式 * 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理论 ” Copyright by Li Xinliang Copyright by Li Xinliang * 知识回顾 1） 守恒型方程与守恒型格式 守恒型方程： 散度型 守恒型格式： 差量型 习惯写为 仅为记号， 与j+1/2点上的值无关！ 守恒型方程+ 守恒型格式=解守恒 “解守恒”： 数值解的积分误差为0 （如果边界准确） 保证总量（总质量、总动量、总能量）严格守恒 （无误差） 边界点 早期的CFD:极为重视守恒性； 目前CFD: 仍很重视守恒性 难点——复杂非线性系统的守恒性很难保证 中间项全部消去，只剩两端 Copyright by Li Xinliang 2） 通量分裂—— 便于使用迎风格式 方法（1）： 逐点分裂 （Steger-Warming, Van Leer, L-F） 原理： 利用了性质 使得 的Jocabian阵特征值纯正或纯负 优点： 无需矩阵运算，计算量小，使用方便 不足： 仅重新组合，没有做到真正解耦。 原因： 具体方法： Steger-Warming L-F Van Leer or Copyright by Li Xinliang * 方法2） 特征（投影）分裂 在网格基上冻结系数 … j-2 j-1 j j+1 … 在基架点上系数 不变 优点：特征分解，（局部）解耦——耗散小，数值振荡低 缺点：大量矩阵运算，计算量大 每计算一个点的导数，要进行m次矩阵运算 （m为网格基点数） 原理描述 （非守恒，很少采用； 实际上使用下一页的方法） 守恒型方式计算 … j-2 j-1 j j+1 … 在基架点上系数 不变 具体步骤： 假设已知 U， 且针对模型方程（线性单波方程） 已构造出差分格式 （1） 1） 计算出 教材130页的公式(6.1.11-6.1.13), 式中用到各变量在j+1/2的值（例如 ） 可使用j, j+1 点值的算术平均 （如 ） 或Roe平均 ； 由 计算；方法很多，例如前面介绍的 或 * Copyright by Li Xinliang 均可 推荐使用Roe-平均！ 2） 在网格基上计算 … j-2 j-1 j j+1 … 计算fj+1/2用到的点 注意，在该网格基上（例如k=j-1,j,j+1） 保持不变 例如： 3) 利用已构造好的差分格式，计算通量 4） 得到总通量 5） 计算差分 （j点处） 步骤的算法描述 （注意： 实际上是两重循环） do j=1,N do k=j-1,j+1 (网格基，可以是更多或更少点) enddo enddo do j=1,N enddo 需要多次矩阵运算，计算量大 守恒性好，耗散小，数值解质量好 * Copyright by Li Xinliang 通量分裂+迎风差分 ? 引入数值耗散 分裂本身不带来耗散， 但会放大（或减少）差分的耗散 举例： 分裂过程 耗散 如果差分格式无耗散（例如都用中心差分），则通量分裂不带来耗散。 = + 向上平移 向下平移 分裂后的流场越偏离原先流场，则总体耗散越大 如使用低精差分度格式， 则对分裂形式敏感 （推荐使用特征分裂） 如使用高精度格式（低耗散），则对分裂形式不敏感 （可使用逐点分裂） * Copyright by Li Xinliang 概念澄清 耗散放大系数 Copyright by Li Xinliang * § 7.1