20章-n-s方程的解举例.ppt

N-S 方程解举例 20.1 引言 ? 粘流解的第三种选择 完整的雷诺平均 N-S 方程数值解 ? 计算流体动力学（ CFD ） 快速发展的技术方法研究 20.2 方法 ? N-S 方程 ? 连续方程 ? 动量方程 能量方程 注意： ? 时间导数在左，空间在右 ? 该形式适合于方程的 时间相关解 20.3 数值解举例 共性 ? 采用 MacCormack 技术 —— 一种时间相关解 法 ? 应用 Baldwin-Lomax 紊流模型，一个 0 方程 湍流模型 —— 湍流计算 ? 应用数千 - 百万网格点，需应用大规模数字 计算机 20.3.1 后向台阶绕流 绕超音速后向台阶的流动 - 速度矢量图 绕超音速后向台阶的流动 等温线 ? 等温线 20.3.2 绕翼型的流动 ? 采用围 绕翼型 非矩形 网格 ? N-S 方 程组变 换到曲 线坐标 系统 计算简介： Wortmann 翼型， 0 攻角，马赫数 M=0.5 ，雷诺数 Re=10000 ? 层流流线图 ? 层流：流动上下翼面均分离 湍流流线图 ? 湍流：无分离 ? 湍流抵抗流动分离的 能力更强 20.3.3 全 机绕流 ? 绕 X-24C 高 超音速试验 机的高超音 速流动 ? 采用 MacCormack 技术——一 种时间相关 解法 20.3.4 激波 / 边界层干扰 激波 / 边界层干扰 ? 激波切入边界层后平板边界层的增长过程： 穿过激波压力升高很大 , 作为很大的逆压梯度作用到边界层； 引起边界层当地分离； 由于逆压梯度通过边界层的低压区向上游传递， 分离发生 在切入激波切入点的前面； 分离边界层诱导一个激波（ 诱导分离激波 ）； 分离边界层向表面偏形成 一系列膨胀激波 ； 最后分离 边界层再附于表面 ，并形成 再附激波 ；再附点边 界层变得很薄，压力变高，为一高的气动加热区； 分离激波和再附激波在离开表面后 合并形成反射激波 （经 典无粘理论结果）； 干扰与流态（层流或湍流）有关：层流干扰更易发生；后 续结果更严重。 数值模拟结果（ Ma=3 ， Re=10 6 ） N-S 方程， B-L 湍流模型 ? 压力在干扰区前增加，穿过分离区达到一个平台， 在再附点前再增加； ? 表面剪应力直降到 0 ，改变方向，经过相当复杂的 变化，在再附点附近恢复到正值。 20.3.5 绕有表面突起的翼型流动 ? 流动条件： NACA0015 翼型， Ma=0.15 ， Re=1.5 × 10 6 ? N-S 方程，时间推进 有限体积代码 (OVERFLOEW,NASA 发 展）， 1 方程湍流模 型； 重叠网格技术 右图为网格结构图 网格全貌 下表面突起的网格 下表面突起的网格（放大） 计算结果速度矢量图