本节书摘来自异步社区《ansys fluent 16.0超级学习手册》一书中的第2章,第2.6节,作者: 唐家鹏 更多章节内容可以访问云栖社区“异步社区”公众号查看。
1.网格导入与处理
(1)在“开始”程序菜单中运行fluent主程序,出现fluent launcher对话框,选择3d,其他保持默认设置,单击ok按钮进入fluent界面,读入网格文件jointpipe.msh,如图2-33所示。

https://yqfile.alicdn.com/0bf0c40ee979cce40088debcffd4b62d578e0b64.png" >
(2)检查网格文件。单击general面板中的check按钮,如图2-34所示,对网格进行检查。需保证网格最小单元体积minimum volume不小于0,即没有负体积网格。
(3)设置网格区域尺寸。单击general面板中的scale按钮,在mesh was created in下拉列表中选择mm,在view length unit in下拉列表中选择mm,如图2-35所示。完成后单击close按钮关闭对话框。
(4)执行mesh→reorder→domain命令,如图2-36所示。对计算域内的网格进行重新排序,可加快求解速度。这个操作有可能需要重复操作几次,直至命令窗口中显示bandwidth reduction小于或等于1为止。
(5)将网格转化为polyhedra型网格(多面体型)。执行mesh→polyhedra→convert domain命令,如图2-37所示。转化后的网格如图2-38所示。
2.设置物理模型和材料
(1)设置求解器类型。本例选择基于压力的求解器(pressure based),求解定常流动。设置重力加速度为y轴负方向,大小为9.8 m/s²,如图2-34所示。
(2)选择湍流模型。单击项目树中的models,打开models面板,如图2-39所示。双击viscous laminar项,打开viscous model对话框。
在viscous model对话框中的model列表中选择k-epsilon两方程模型,对话框将自动扩展成图2-40所示的对话框。在k-epsilon model选项区中选择realizable,其他保持默认设置。单击ok按钮关闭对话框。
(3)设置流体材料。本例中的流体是水,为不可压缩流动。单击materials项,打开materials面板,双击面板中的fluid项,打开create/edit materials对话框,如图2-41所示。

https://yqfile.alicdn.com/e4ff2c64ed0b5b552608d2f9dac034b55ee3077f.png" >
面板中的默认流体材料为空气(air),固体材料为铝(aluminum)。单击fluent database按钮,打开fluent database materials对话框,在fluent fluid materials列表中选中water-liquid(h2o),即液态水,单击copy按钮将材料参数复制到当前材料库中,单击close按钮,关闭fluent database materials对话框,如图2-42所示。再单击close按钮,关闭create/edit materials对话框。
单击项目树中的cell zone conditions项,打开cell zone conditions面板。双击zone列表中的branch-fluid项,打开fluid对话框,如图2-43所示。在material name右侧的下拉列表中选择water-liquid,单击ok按钮,关闭fluid对话框。

https://yqfile.alicdn.com/c3ca7c3539dc8f3afe9dfcbd5313325c79dde48d.png" >
3.设置操作环境和边界条件
(1)设置操作环境。执行define→operating conditions命令,打开operating conditions对话框,如图2-44所示。在operating pressure下面的输入框中输入操作压强为0 pa,其他保持默认设置。单击ok按钮,关闭operating conditions对话框。
(2)设置边界条件。单击boundary conditions选项,打开boundary conditions面板。选中zone列表中的inlet。inlet是流量入口边界条件。在inlet被选中的状态下,在type下拉列表中选中mass-flow-inlet边界条件,如图2-45所示。
在inlet被选中的状态下,单击edit按钮,弹出mass-flow-inlet对话框,在该对话框中设置流量入口边界条件的各项参数,如图2-46所示。
在mass flow rate中输入605.6 kg/s,初始表压设置为200 000 pa。
在specification method下拉列表中选择湍流强度的定义方法,此处选择intensity and viscosity ratio,即用湍流强度和黏性比定义湍流,其值采用默认值10。
单击ok按钮,关闭流量入口边界条件设置对话框。
用同样的方法设置出口边界条件。出口边界条件为outflow,即自由出口边界条件,如图2-47所示。保持默认设置,单击ok按钮,关闭outflow对话框。

https://yqfile.alicdn.com/0448670112a64d59771fc0a6e33a71c0e9afa436.png" >
4.设置求解方法和控制参数
(1)设置求解方法。单击项目树中的solution methods选项,打开solution methods面板,对求解方法进行设置,选择求解的方程类型和微分方程离散格式。
在压力-速度耦合方式(pressure-velocity coupling)下拉列表中选择simple,即采用simple算法。
在微分离散格式选项中,梯度(gradient)选择green-gauss cell based,压力(pressure)采用presto!格式,动量方程(momentum)选择二阶迎风格式(second order upwind),湍流脉动能量和湍流耗散率也采用二阶迎风格式。设置完成后的solution methods面板如图2-48所示。
(2)设置求解控制参数。单击项目树中的solution controls选项,打开solution controls面板,对求解过程中的控制参数进行设置,设置各项松弛因子如图2-49所示。

https://yqfile.alicdn.com/288de8e097e54d3d149953b57a4986448a6b2494.png" >
5.设置监视窗口和初始化
(1)设置残差监视窗口。单击项目树中的monitors选项,打开monitors面板,双击monitors面板中的residuals项,打开残差监视器设置对话框。选中print to console和plot选项,在convergence criterion下拉列表中选择none,即不自动进行收敛判断。其他保持默认设置,如图2-50所示。
(2)设置出口速度监视窗口。单击项目树中的monitors项,打开monitors面板,单击monitors面板中surface monitors列表下面的create按钮,打开surface monitor对话框,如图2-51所示。选中print to console、plot和write选项。设置监视器窗口为2号窗口。
报告类型(report type)设置为facet average,即表面平均值。
在场变量(field variable)下拉列表中选择velocity。
在surfaces列表中选中outlet。

https://yqfile.alicdn.com/ec7bca6205bb004e3538b4005b94a4f231d31b7f.png" >
(3)初始化流场。单击项目树中的solution initialization选项,打开solution initialization(流场初始化)面板,如图2-52所示。在compute from下拉列表中选中inlet,表示整个流场中的初始状态与边界inlet上的流场状态相同。单击initialize按钮,完成流场的初始化。

https://yqfile.alicdn.com/61d56375fdbc910d03ba225d3e4fbff112d838fd.png" >
6.求解
(1)开始迭代。单击项目树中的run calculation项,打开run calculation面板,在number of iterations下面的输入框中输入500,即迭代500步,如图2-53所示。单击calculate按钮,开始计算。
(2)迭代开始后,图形窗口中会动态显示残差值和出口平均水流速度随迭代过程的变化曲线。当迭代完500步后,出口速度收敛曲线如图2-54所示。由收敛曲线可以看出,出口平均水流速度基本不变化,可以认为计算已经收敛。

https://yqfile.alicdn.com/d1e91bfe1a92e6d604e12779fe6bfa52259c4075.png" >
7.保存计算结果并退出
(1)执行file→write→case&data命令,保存案例文件和计算结果。
(2)退出fluent。单击fluent界面右上角的图片 56(关闭)按钮,在弹出的warning对话框(见图2-55)中,单击ok按钮,退出fluent。