任立峰双圆柱绕流解读.docVIP

任立峰双圆柱绕流解读 任立峰双圆柱绕流解读 PAGE / NUMPAGES 任立峰双圆柱绕流解读 二维并排双圆柱绕流数值模拟 纲要：为研究不行压缩流动中二维的并列双圆柱在不一样间距和流速下的流动状况和影响要素。 选用间距分别为 3m，1m 两个距离成立模型在选用速度为 0.2m/s；0.6 m/s；1 m/s 进行模拟，对不一样距离状况下的速度云图速度轮廓图得出了不一样距离状况下的流动状况， 同时研究了同一模型下不一样速度对流动状况的影响。并经过查阅文件考证了用 fluent 模拟的正确性。重点字： 二维双圆柱；数值模拟；不行压缩流动；数值传热学 1.前言 数值传热学， 又称计算传热学， 是指对描绘流动与传热问题的控制方程采纳数值方法， 经过 计算机求解的一门传热学与数值方法相联合的交错学科。 数值传热学的基本思想是把本来在 空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等） ，用一系列有限个离 散点上的值的会合来取代，经过必定的原则成立起这些失散点变量值之间关系的代数方程 （称为失散方程） 。求解所成立起来的代数方程已获取求解变量的近似值。 当前，比较流行的数值模拟剖析应用软件有 FLUENT 、CFX 、STAR-CD 、和 PHOENICS 等， 而 FLUENT 是国内外比较流行的商用 CFD 软件包，该软件以其市场据有率高、计算正确、界面友善、使用简单、应用领域广、物理模型多而获取较高的市场据有率和用户的必定。 成立模型 2.成立模型 ： 以下均采纳二维非耦共计算模式， 采纳非定常计算模式； 动量方程失散模式为二阶迎风 模式；残差控制为。 关于不行压缩粘性流体，在直角坐标系下，其运动规律可用 N-S 方程来描绘，连续性 方程和动量方程分别为 ? uj = 0 （ 1-1） ?x j ? u i + ? 1 ? P ? ?u j ) （ 1-2） u u = - +(v ? t ρ? xi ?x j ?x j ?x j 界限条件 u = U, v = 0 （ 1-3） 采纳双圆柱并排摆列的方式选择两园间距 1 米和 3 米即 1 倍直径和 3 倍直径成立模型此模型 长为圆的直径的 32 倍宽为圆的 16 倍直径圆心距进口距离为 6 米其模型以下列图。 图 1 间距为 1 倍直径的模型 图 2 间距为 3 倍直径的模型 3 用 gambit 成立模型并进行网格区分和设定界限条件 3.1 网格的区分 对圆柱的四周用线进行切割线距离圆柱的圆心的距离为 2m 即 2D，在圆柱四周围城的 矩形框内对网格进行独自的区分，在圆柱四周采纳 pava 进行区分，对圆柱的矩形框外的使 用 map 进行区分。对边采纳等距离进行区分每个节点间距 0.05m 对两圆间距 1 倍直径的进 行网格区分共获取 12868 个网格如图 3。同理对间距为 3 倍直径的模型进行网格区分共获取 12846 个网格如图 4。 图 3 间距为 1m 的模型 图 4 间距为 1m 的模型 3.2 界限条件的定义 对计算地区进行界限条件定义，选中左侧线为进口边线，在 name 文本框中输入 inlet ， 将 type 栏选为 Velocity-inlet ，在 Entity 栏宣威 Edges。选用右边线为出口边线，在 name 栏 目中输入 Outlet ，将 Type 栏选为 Outlet ，在 Entity 中选用 Edges。而后输出网格文件。 4 使用 fluent 进行解算 4.1 对间距为 3m 的模型进行解算 采纳的水作为流体为流经计算地区的流体，从进口的速度分别设定为 0.2m/s， 0.6m/s 和 1m/s 分三个速度对这个模型。 第一导入并检查网格，而后单击 display— grid 显示网格。第二步选择计算模型，在 define— models— viscous 选择 k-epsilon 模型，设置求解器是翻开 solver 对话框，在 Time 下 选择 Unsteady （非定常），在 Unsteady Formulation 中选择 2nd-Order-Implicit 。在 Porous Formulation 中选择 Superficial Velocity 。保持其余默认设置不变。 将流体资料设置为水， Fluent 中有这一资料，直接复制即可。设置界限条件时先设置 fluid 流体地区的界限条件，再设置 进口速度界限条件，将进口速度 velocity magnitude 栏内输入 0.2，右边栏内选择 constant。 并保存出口和墙壁的默认设置。 翻开残差设置对话框， 选择 option 下边的 plot ，convergence 下边对应的数值均为 0.001