CFDFoundation.ppt

数值流体仿真基础 No.1 一元对流扩散方程的解法 张 政@Japan 2011年6月25日 流体力学的基础方程 一元对流扩散方程的解法 一元对流扩散方程的解法 一元对流扩散方程的解法 解法流程图 例程的动作确认 什么是扩散 扩散方程的数值稳定性 什么是对流 对流方程的数值稳定性(1) 对流方程的数值稳定性(2) 数值流体仿真基础 No.2 一元Poisson方程的解法 张 政@Japan 2011年6月25日 一元Poisson方程 一元Poisson方程的解法(1) 一元Poisson方程的解法(2) 数值流体仿真基础 No.3 二维非压缩性流体的解法 张 政@Japan 2011年6月25日 计算对象（二维管内层流） 二维非压缩性流体的基础方程 SMAC法的压力Poisson方程的推导(1) SMAC法的压力Poisson方程的推导(2) 边界条件 SMAC法的计算步骤（流程图） 计算对象（二维管内层流） SMAC法的计算步骤（流程图） 计算结果（流入速度一定 1/2） 计算结果（流入速度一定 2/2） 边界条件 计算结果（周期边界条件 1/2） 计算结果（周期边界条件 2/2） 数值流体仿真基础 No.4 k-epsilon紊流模型 张 政@Japan 2011年6月25日 紊流的特征 计算对象（二维高Re流） 计算结果（不使用紊流模型） 雷诺平均方程 紊流模型（雷诺应力模型） 紊流模型（涡粘性模型） 紊流模型（涡粘性模型 0方程模型） 紊流模型（涡粘性模型 2方程模型） 紊流模型（涡粘性模型 2方程模型的壁面处理） 计算对象（二维高Re流 使用紊流模型） k-epsilon紊流模型的离散化 SMAC法的计算步骤（流程图） 边界条件（k 和 epsilon） 计算结果（k-epsilon紊流模型） 紊流计算方法的总结 END t V t V 层流 Re↑ 紊流 雷诺数　　　　　　　的意思？代表长的选择？（質問） 由于涡的存在而发生速度变动 　雷诺数越高、边界层就越不稳定、边界层附近开始产生涡旋（紊乱）、当这个紊乱发展并生长到整个流场区域（如管的全断面，最大涡尺寸：与管径相当）的同时、涡旋互相冲撞、涡旋的尺寸也随之变小（Cascade、最小涡尺寸：Kolmogoroff Scale）、消失。 非常靠壁面的地方由于受到壁面（受分子粘性支配）的影响、难以产生涡旋、并且涡旋也较容易消失、在这个区域速度分布呈线性变化、称为粘性底层。而在边界层外部的区域被称之为外层、平均流速基本均一、紊乱也具有各方同性的特点。外层与粘性底层之间的区域的速度分布遵循对数法则、被称为对数区域。这个区域受到壁面的影响、具有各方非同性的特点。 最小涡尺寸 20 300 粘性底层 对数区域 外层 内层 层流边界层 过渡层 紊流边界层 粘性底层 计算对象 　　如下图所示的22（I方向）×12（j方向）的领域为对象、格子大小为 Δx = Δy = 5 mm 。边界条件为从左侧流速U = 1 m/s 的流入边界、上下壁面边界、右侧自由流出边界。计算区域的流体为常温常压的水（密度: 1000 kg/m^3, 粘性: 0.001Pa?s）。（Re = 50000, 紊流） 12 22 流入 流出 壁面 壁面 均一流 U = 1.0 m/s 結果 　不使用紊流模型的计算所得到的紊流速度分布与层流一样。不能正确反映紊流速度分布、因此需要采取相应的手法来解决。 解决方法： 　１．直接数值仿真（DNS） 　　　NS方程不论层流还是紊流都能够适用， 如果计算格子分割的比最小涡尺寸还要 小的话，就能够再现紊乱所造成的速度 变动 　２．追加紊流模型 即：对紊乱（速度变动）采用适当的 模型来表示 層流 乱流 例程 運動量保存式 质量守恒式 雷诺平均方程式 一般形式 t V 速度变形 平均速度变形 分子粘性应力 紊流应力（雷诺应力） 求解紊流问题需要紊流应力的评价方法（即：紊流模型）。 t V 層流 ↑ 乱流 雷诺平均 运动量守恒式 与层流同样 雷诺应力模型（RSM）使用以下步骤计算紊流应力。 　１．一般形式的运动量守恒式　－　雷诺平均运动量方程式　＝　变动速度的运动量方程式 　２．通过变动速度的运动量方程式导出以下的雷诺应力方程式、计算紊流应力。 乱流応力 （由变动速度构成） 生成项 压力－变形项 散逸项 扩散项 求解紊流问题需要紊流应力的评价方法（即：紊流模型）。 雷诺平均 运动量守恒式 与层流同样 乱流応力 涡粘性模型通过假定涡旋所产生的紊流应力与分子粘性应力具有相似性来计算紊流应力。 分子粘性応力 分子粘性系数 紊流粘性系数 需要紊流粘性系数的计算方法 也就是说、雷诺