第五章 工程应用实例.docVIP

第五章 工程应用实例 5.1 圆柱绕流数值模拟 一、问题描述 把一个圆柱体放在静止的流体中，然后流体以很低的速度（雷诺数很小）绕流圆柱体。此时流动与理想流体绕流圆柱一样，流体在前驻点速度为零，而后沿圆柱体两侧流动，流动在圆柱体的前半部分是降压，速度逐渐增大到最大值，而后半部分是升压，速度逐渐下降，到后驻点重新等于零。增大来流速度（雷诺数增大），圆柱体后半部分的压强梯度增加，以致引起边界层的分离。随着来流速度的不断增加，圆柱体后半部分边界层中的流体微团受到更大的阻滞，分离点一直向前移动。当雷诺数（Re）增大到40左右时，在圆柱体的后面便会产生一对旋转方向相反的对称旋涡。当雷诺数超过40后，对称旋涡不断增长且出现摆动，到雷诺数大约等于60时，这对不稳定的旋涡分裂，最后形成有规则的、旋转方向相反的交替旋涡，这就是卡门涡街。 本节将模拟一直径为2cm的圆柱在雷诺数Re为200时卡门涡阶脱落情况。本节涉及以下内容： 1、二维结构化网格的划分； 2、非稳态计算； 3、动画设置； 4、后处理：压力云图、速度矢量图； 5、数据保存 二、GAMBIT建模 第1步：创建点 操作：点击——，如图1。在X，Y，Z上分别输入流体域的边界点的值。在这里取x最小值为-6，最大值为20，y的最小值为-4，最大值为4。 对于二维圆柱的创建只需要输入圆心坐标和不在同一条直线上的二个点坐标值即可。创建完成后如图2所示。 为了帮助画结构化网格，我们需要在圆柱外面加上一个长为4的正方形，取四个点（-2，2），（-2，-2），（2，-2），（2，2）以及点（-2，4），（-2，-4），（2，4），（2，-4）。点创建完成后如图2所示。 图1 GAMBIT工具栏 图2 点创建完成 第2步 由点创线 操作：点击—右键单击会出现下拉菜单，在下拉菜单中选择，打开“Creat Full Circle”对话框如图3所示。 图3 创建圆对话框 按照提示在Center 后选择圆心点，在End-Point 后面选择圆周上的两个点，点击Apply。 图4 创建圆 右键点击，在出现的下拉菜单中选择，创建流体域的外围边界线。依次点击外围各点连接成直线。如图5所示。 图5 创建线 第3步 由线建面 3.1 四等分圆 为了下一步对圆划分结构化网格需要将圆分成四段。此时需要连接圆心与正方形的四个顶点。连接完成后如图6所示。 图6 创建等分线 打断线操作：——，在图7的菜单栏Split with 后面选择Edge。第一个Edge选择需要打断的线，第二个Edge选择去打断的线。在这里，我们选择圆为需要打断的线，依次选择四条斜线将其打断为四段。打断后的图形见图8。依次连接各打断点和对应的四个角点。 图7 打断对话框图 图8 圆被打断 3.2 由线创建面 由GAMBIT进行线创建面的操作中，每个面都需要用几条封闭的曲线构成。因此在创建时需要注意选择的线段封闭。 操作：—— 在Edge后面选择需要构成面的几条线段。注意：圆由四个面构成。 图9 创建面对话框 图10 由线构建面 第4步 划分网格 4.1 线网格划分 1、划分圆的网格。 操作：—— 出现图11所示对话框。 图11 线网格对话框 在Edge 后的黄色框里选择一段圆弧，点击Spacing后面的Interval size 将其改成Interval Count，然后输入30；然后依次将各段圆弧划分网格，内部网格节点均输入30。 图12 一段圆弧划分网格 每段圆弧对应了一段正方形的边，将每段正方形的边网格节点数也设置为30，以划分结构化网格。 图13 正方形网格划分 3、斜边网格 对于四条斜边的网格划分需要局部进行局部加密，即在圆柱表面网格要密，相对远一些的地方网格可以稀疏一些。参数设置如图14所示，选择一段斜线，在Ratio后输入1.05，即表示网格节点以1.05的比例增长。结果如图15所示。 图14 比例操作 图15 比例网格 四条斜边均按1.05的比例进行网格划分，网格数为25。 4、其余边网格 别的边的网格划分要相对简单一些，只要依据网格节点之间的距离合适，对边的网格节点数相等就可以。线网格划分完成以后如图16所示。 图16 线网格划分完成 4.2 面网格划分 操作：—— 选择一个面，ELEMENTS选择QUAD（四边形网格），type选择MAP(结构化网格)，点击APPLY。此时划分出来的网格为结构化四边形网格。依次划分各个面，均划分为结构化网格。如图18所示。 图17 面网格划分对话框