FLUENT6.1全攻略
第二章 FLUENT的计算步骤
本章通过一个稍微复杂一些的算例再次演示FLUENT的求解过程。这个算例的内容是计算一个二维弯管中的湍流流动和热传导过程，在这个算例中可以看到FLUENT计算的标准流程，其中包括：
（1）如何读入网格文件。
（2）如何使用混合的单位制定义几何模型和物质属性。
（3）如何设定边界条件和和物质属性。
（4）如何初始化计算并用残差曲线监视计算进程。
（5）如何用分离求解器计算流场。
（6）如何用FLUENT的图形显示功能检查流场。
（7）如何用二阶精度离散格式获得更高精度的流场。
（8）以温度梯度为基准调整网格以提高对温度场的计算精度。
2.1 问题概述
图2-1 弯管流动图示
如图2-1所示，温度为26℃的冷流体流过弯管，温度为40℃的热流体从转弯处流入，
1
FLUENT6.1全攻略
并与主流中的冷流体混合。管道的尺寸如图2-1所示，单位为英寸，而边界条件和流体材料性质则采用国际单位制。入口处的雷诺数为2.03 x 105，因此必须使用湍流模型。
2.2 处理网格
网格处理包括网格的输入、检查、光顺、比例转换和显示等操作，下面分别进行介绍。
2.2.1读入网格文件
首先启动FLUENT的2D版，然后读入网格文件：
File -> Read -> Case...
这个算例的网格文件可以在FLUENT6.1为用户提供的文档光盘中找到，路径是：
cdrom:\fluent6.1\help\tutfiles\elbow\elbow.msh
2.2.2检查网格
执行下列菜单操作，进行网格检查：
Grid -> Check
此时控制台窗口中会显示与网格有关的信息，包括网格空间范围、体积信息、表面积信息、节点信息等等。网格中存在的任何错误都会出现在这个信息报告中，其中最需要检查的是网格单元的体积不能为负值，否则计算将无法继续下去。
图2-2 Smooth/Swap Grid（光顺/转换网格）面板
2
FLUENT6.1全攻略
2.2.3光顺并转换网格
执行下列菜单操作，打开Smooth/Swap（网格光顺和转换）面板：
Grid -> Smooth/Swap...
光顺网格可以提高网格质量，提高计算精度。在启动光顺/转换网格面板后，点击Smooth （光顺）按钮，然后再点击Swap（转换）按钮，控制台上将报告被转换网格的数量。反复点击Smooth（光顺）按钮和Swap（转换）按钮，直到报告中被转换（Swap）的网格数量降低为零，则网格光顺处理过程结束。
2.2.4按比例调整网格
执行下列菜单操作，打开Scale Grid（按比例调整网格）面板：
Grid -> Scale...
图2-3 Scale Grid（按比例转换网格）面板
在Units Conversion（单位转换）下面的Grid Was Created In列表中选择in（英寸）作为网格创建数据单位。
点击Scale（比例）按钮进行比例计算。注意在计算域范围（Domain Extents）中，网格的单位显示仍然是国际单位制中的m（米）。
点击Change Length Units（改变长度单位）按钮，可以将英寸变为长度单位，同时可以从面板下方计算域范围中看到x、y坐标的最大值改变为64英寸。
3
FLUENT6.1全攻略
图2-4 Grid Display（网格显示）面板
图2-5 网格显示
4
FLUENT6.1全攻略
2.2.5显示网格
执行下列菜单操作，打开Grid Display（网格显示）面板，如图2-4所示：
Display -> Grid...
选择所有表面，然后点击Display（显示）按钮，FLUENT随即打开一个新的显示窗口显示网格，如图2-5所示。
用右键点击图形显示窗口中的一个边界，则与这个边界相关的信息就会显示在控制台窗口中。在网格区域很多时，用这种方法可以比较方便地分辨各个区域。
2.3 计算模型
图2-6 Solver（求解器）面板
5
FLUENT6.1全攻略
2.3.1设置求解器参数
执行下列菜单操作，打开Solver（求解器）面板：
Define -> Models->Solver...
求解器的缺省设置如图2-6所示，本算例中可以保持这些缺省参数设置。
图2-7 粘性模型面板
2.3.2选择湍流模型
执行下列菜单操作，打开Vicous（粘性）面板，如图2-7所示：
6
FLUENT6.1全攻略
7
Define->Models->Viscous...
在Model 列表中选择k-epsilon (2 eqn)，即ε?k 两方程模型。点击OK 按钮接受系统
对ε?k 模型的缺省设置，包括标准模型设置、壁面函数设置和湍流模型常数的设置等等。
2.3.3启用能量方程计算
执行下列菜单操作，打开Energy （能量）面板，并激活能量计算，如图2-8所示：
Define->Models->Energy...
图2-8 能量方程面板
2.4 定义材料性质
创建一种新的物质，命名为water （水），并设置相关参数。
首先执行下列菜单操作，打开Materials （材料）面板，如图2-9所示：
Define -> Materials...
在Name （名称）中键入材料名称water （水），然后设定水的各种物理性质，其中Density （密度）为1000kg/m 3， C p （定压比热）为4216J/kg-K ，thermal conductivity （导热系数）为0.677W/m-K ，viscosity （粘度）为8x10-4kg/m-s 。最后点击Change/Create （修改或创建）按钮完成创建过程。在FLUENT 提示是否用新的材料数据覆盖air 的数据时，点击NO 按钮予以否认。
实际上FLUENT 的材料性质数据库中已经存在水的数据，因此也可以简单地将相关数据直接从Database （数据库）中调入，如果个别数据与前面的数据有出入，可以直接在面板上进行修改。
设置完毕后，关闭Materials （材料）面板。
FLUENT6.1全攻略
图2-9 Materials（材料）面板
2.5 定义边界条件
执行下列菜单命令，打开边界条件面板，如图2-10所示：
Define -> Boundary Conditions...
2.5.1设定流体条件
在Zone（区域）列表中选择fluid-9，Type（类型）列表中相应地反白显示为fluid（流体）。
8
FLUENT6.1全攻略
图2-10 边界条件面板
图2-11 Fluid（流体）面板
9
FLUENT6.1全攻略
点击Set...（设置）按钮打开Fluid（流体）面板，如图2-11所示。在Material Name（材料名称）列表中选择water（水），点击OK按钮结束设置过程。
2.5.2设置速度入口边界条件。
（1）设置velocity-inlet-5上的边界条件。
在Zone（区域）列表中选择velocity-inlet-5，然后点击Set...（设置）按钮，进入速度入口设置面板，如图2-12所示。
图2-12 Velocity Inlet（速度入口）面板
在Velocity Specification Method（速度定义方法）中，选择Components（分量形式）定义速度，然后将X-Velocity（X向速度分量）设为0.2m/s。将Temperature（温度）设为293K。在Turbulence Specification Method（湍流定义方法）中，选择Intensity and Hydraulic Diameter（强度和水力直径）定义湍流，然后将Turbulence Intensity（湍流强度）设为5％，将Hydraulic Diameter（水力直径）设为32in。点击OK按钮结束对velocity-inlet-5的设置。
（2）设置velocity-inlet-6上的边界条件，方法同上，相关参数如下：Y-Velocity（Y向速度分量）为1.0m/s，temperature（温度）为313K，Turbulence Intensity（湍流强度）为5％，Hydraulic Diameter（水力直径）为8in（英寸）。
10
FLUENT6.1全攻略
2.5.3设置压强出口边界条件
pressure-outlet-7上的压强出口边界条件的相关参数为：表压0 pascal，回流总温为300 K，Backflow Direction Specification Method（回流方向定义方法）为Normal to Boundary（垂直于边界），Turbulence Specification Method（湍流定义方法）为Intensity and Hydraulic Diameter（湍流强度与水力直径），Backflow Turbulence Intensity（回流湍流强度）为5％，Backflow Hydraulic Diameter（回流水力直径）为32in。
图2-13 Pressure Outlet（压强出口）面板
2.5.4壁面条件设置
壁面wall-4的边界条件保持缺省设置，其中Heat Flux（热流通量）的值为0。壁面wall-8也保持缺省设置。
2.6 求解过程
参数设置完毕并检查无误后，既可以开始流场计算。流场计算的流程如下：
11
FLUENT6.1全攻略
2.6.1初始化
执行下列菜单操作打开Solution Intialization（求解初始化）面板，如图2-14所示：Solve->Initialize->Initialize...
图2-14 Solution Initialization（初始化解）面板
在Compute From（计算起始于）列表中选择velocity-inlet-5，然后在Y Velocity（Y向速度分量）中加入0.2m/sec，点击Init（初始化）按钮开始初始化，最后点击Close（关闭）按钮关闭初始化面板。
2.6.2设置残差监视器
执行下列菜单操作，打开Residual Monitors（残差监视器）面板：
Solve->Monitors->Residual...
在Options（选项）下选择Plot（绘图），点击OK按钮关闭面板。
2.6.3保存算例文件
执行下列菜单操作保存算例文件：
12
FLUENT6.1全攻略
File->Write->Case...
不用改动Write Binary Files（写二进制文件）的缺省设置，这样算例文件可以用二进制格式记录在硬盘上。
图2-15 Residual Monitors（残差监视器）面板
2.6.4打开迭代面板
执行下列菜单操作，打开Iterate（迭代）面板：
Solve -> Iterate...
将Number of Iterations（迭代次数）设为100，点击Iterate（迭代）按钮开始计算。在开始计算时残差监视器将同时被打开。
2.6.5监视收敛过程
图2-16是计算进行到第60步时残差曲线的走势。因为没有普适的收敛判断标准，所
13
FLUENT6.1全攻略
以在观察残差曲线时，不要仅仅监视残差曲线下降的数量级，最好同时能够监视相关流场变量的变化情况。
图2-16 残差监视器
简单地说，可以用三种方法判断计算是否已经收敛：
（1）观察残差曲线。
可以在残差监视器面板中设置Convergence Criterion（收敛判据），比如设为10-3，则残差下降到小于10-3时，系统既认为计算已经收敛并同时终止计算。
（2）流场变量不再变化。
有时候不论怎样计算，残差都不能降到收敛判据以下。此时可以用具有代表性的流场变量来判断计算是否已经收敛——如果流场变量在经过很多次迭代后不再发生变化，就可以认为计算已经收敛。
（3）总体质量、动量、能量达到平衡。
在Flux Reports（通量报告）面板中检查质量、动量、能量和其他变量的总体平衡情况。通过计算域的净通量应该小于0.1%。Flux Reports（通量报告）面板如图2-17所示，其启动方法为：
Report -> Fluxes
14
FLUENT6.1全攻略
图2-17 Flux Reports（通量报告）面板
2.6.6保存数据文件
执行下列菜单操作保存数据文件：
File -> Write -> Data...
2.7 显示计算结果
2.7.1等值线显示
执行下列菜单操作，启动Contours（等值线）面板，如图2-18所示：
Display -> Contours...
首先显示速度值的填充等值线。在Contours of下的列表中选择Velocity...（速度）和Velocity Magnitude（速度值），然后在Options（选项）下面选择Filled（填充），最后点击Display（显示）按钮显示速度的等值线。在等值线图中用右键点击一点，控制台窗口中将显示该点的速度值。
用同样的方法可以显示温度、压强等变量的等值线图。
15
FLUENT6.1全攻略
图2-18 Contours（等值线）面板
图2-19 速度值的等值线图
16
FLUENT6.1全攻略
2.7.2显示矢量图
除了显示等值线图，还可以用矢量图观察流场中的速度矢量。
首先执行下列菜单命令，启动Vectors（矢量）面板：
Display -> Vectors...
图2-20 Vectors（矢量）面板
在Vectors of下选择Velocity（速度），然后点击Display（显示）按钮，就可以在不改变缺省设置的情况下显示速度矢量图。如果对显示效果不满意，可以在Vectors（矢量）面板中修改显示设置，包括改变Scale（比例）和Style（风格）。
如果计算的是三维问题，则需要先在Surfaces（表面）下选择一个显示面，才能进行显示，显示的结果是速度矢量在这个面上的投影。
17
FLUENT6.1全攻略
图2-21 速度矢量图
图2-22 解变量XY绘图面板
18
FLUENT6.1全攻略
2.7.3 XY曲线显示
出口处的温度变化可以用XY曲线形式予以显示。首先启动Solution XY Plot（解变量XY绘图）面板：
Plot -> XY Plot...
在Y Axis Function（Y轴函数）下面的列表中选择Temperature...（温度）和Static Temperature（静温）作为Y轴变量，在Surfaces（表面）列表中选择pressure-outlet-7，点击Plot（绘制）按钮，则静温在出口处的分布曲线就显示在图形窗口中。用同样方法还可以显示压强等变量在出口、入口等处的分布曲线。
图2-23 温度在出口处的分布曲线
2.7.4定义场函数
定义动压头的场函数。首先启动Custom Field Function Calculator（场函数计算器）面板：
Define -> Custom Field Functions...
将新函数的名称加入New Function Name（新函数名称）栏，然后根据动压头2/2v?ρ
19
FLUENT6.1全攻略
的定义公式进行公式编辑。动压头的输入内容如图2-24所示。
在公式中需要输入流场变量（比如密度、速度）时，就在Field Functions（场函数）列表下进行选择；在使用数字和运算操作符的时候，则直接从面板左部的键盘中点取。输入完毕后，点击Define（定义）按钮完成函数定义。
场函数的显示与其他场变量的显示方法是完全相同的。动压头的显示与速度值的显示方法相同：
（1）启动等值线面板。
（2）然后从Contours of列表中选择Custom Field Function...（修改场函数）和dynam-head（动压头）。
图2-24 场函数计算器面板
（3）设置显示选项，点击Display（显示）按钮。
2.8 启用二阶精度离散格式
前面的计算是用一阶精度格式完成的。一阶精度格式的缺点是耗散性很大，计算稳定性好，但是对流场中梯度比较大区域内的解有比较严重的“抹平”现象，因此为了获得精度更高的结果，可以采用二阶精度格式。因为二阶精度格式的稳定性不如一阶精度，所以在采用二阶精度格式的时候要适当减小亚松弛因子。二阶精度格式的具体设置和求解过程如下：
（1）启动Solution Controls（求解控制）面板：
20
二维模拟： 一、模拟类型： 1、 大区域空间速度场模拟 计算区域大小设置：迎风面是建筑长度的3倍，背风面是建筑长度的12倍，两侧面是建筑宽度的3倍，高度是建筑高度的4倍。 根据相似理论：l C -几何比例尺 速度比例尺：2 10l C C =υ 风量比例尺：2520l l Q C C C C =?=υ 热量比例尺： 250l T Q C C C Cq =?=? 2、 建筑户型温度场、速度场模拟 二、基本操作步骤及注意事项： A gambit 建模 1、 建模： 方法一：直接在GAMBIT 建模； 方法二：CAD 导入gambit ； 1） 在CAD 中用PL 线将户型的基本构造画出来，创建为面域； 2） 输入命令acisoutver ，把‘70’修改为‘30’。 3） “文件”——“输出”——sat 文件 4） 在gambit 中导入Acis 文件 注意：在用PL 线构画户型时，在进口和出口边界（窗户、内户门），要各边界端点连续画线。 2、 划分网格： Interval Size ：50 3、 设置边界条件 内部开口边界（门）设置为internal ，房间相邻墙壁设置为Wall 4、 保存文件，并输出mesh 文件 B 导入fluent 计算： 1、 导入mesh 文件 2、 检查网格 3、 设置单位 gambit 里可以缩小建筑比例建模，在fluent 中设置单位恢复原模型。 4、 选择计算模型 5、 设置材料类型 6、 设置边界条件 7、 设置模拟控制条件 8、 边界初始化
9、设置监视窗口 10、设置迭代次数进行计算 11、结果显示 12、保存文件 三、需解决问题： 1、湍流强度等计算； 2、层流湍流界定问题； 3、壁面湿度设置问题； 四、待提高部分： 1、户型流场模拟时，墙壁考虑采用双钱； 2、南京理工校区原始模型（不简化）模拟； 3、三维模型模拟； 五、
实体入水模拟过程 3.2.1利用GAMBIT建立计算模型 1）启动GAMBIT，打开对话框如图3.2.1选择工作目录为D:\GAMBIT working。 图 3.2.1 2）首先建立等边三角形，单击Geometry Vertex Create Real Vertex，在Create Real Vertex面板的x、y、z坐标输入（0,0,0），单击Apply按钮生成第一个点，按同样的方法建立点（0.4,0,0）。然后单击Geometry Edge Create Straight Edge，在Create Straight Edge面板中选择点1与点2，连接这两点省成线段。如图3.2.2 图3.2.2 3）单击Edge面板中的Move/Copy Edges按钮，打开如图3.2.3的面板，选择线段1，单击copy按钮，并选择Operation为Rotate，在Angle栏输入60，其他保持默认，单击Apply 按钮。即旋转复制生成第二条线段。
图3.2.3 4）剩下的一条线段只需连接右侧两点即可，如图3.2.4所示。 图3.2.4 5）创建三角形面。单击Geometry Face Create Face from Wireframe，在Create Face from Wireframe面板中利用鼠标左键框选等边三角形的三条边，然后单击Apply按钮创建面。 6）由于三角形面域的位置不对，所以还要对其位置进行调整。首先需将其旋转210度。单击Face面板中的Move/Copy Faces按钮，在Move/Copy Faces面板中，选择面1（face.1），单击Move并选择Operation为Rotate，在Angle栏输入210，其他保持默认，单击Apply 按钮。其次，需要将三角形平移，在Move/Copy Edges面板中选择面1（face.1），单击Move 并选择Operation为Translate，在x与y栏分别输入3和8.4，单击Apply按钮完成平移操作，此时的视图窗口如图3.2.5所示。
Fluent数值模拟的主要步骤 使用Gambit划分网格的工作： 首先建立几何模型，再进行网格划分，最后定义边界条件。 Gambit中采用的单位是mm，Fluent默认的长度是m。 Fluent数值模拟的主要步骤： （1）根据具体问题选择2D或3D求解器进行数值模拟； （2）导入网格(File-Read-Case)，然后选择由Gambit导出的msh文件。 （3）检查网格(Grid-Check)，如果网格最小体积为负值，就要重新进行网格划分。（4）选择计算模型(Define-Models-Solver)。（6） （5）确定流体的物理性质(Define-Materials)。 （6）定义操作环境(Define-Operating Conditions)。 （7）指定边界条件(Define-Boundary Conditions )。 （8）求解方法的设置及其控制（Solve-Control-Solution）。 （9）流场初始化（Solve-Initialize）。 （10）打开残插图（Solve-Monitors-Residual）可动态显示残差，然后保存当前的Case和Data文件（File-Writer-Case&Data）。 （11）迭代求解（Solve-Iterate）。 （12）检查结果。 （13）保存结果（File-Writer-Case&Data），后处理等。 在运行Fluent软件包时，会经常遇到以下形式的文件： .jou文件：日志文档，可以编辑运行。 .dbs文件：Gambit工作文件，若想修改网格，可以打开这个文件进行再编辑。 .msh文件：Gambit输出的网格文件。 .cas文件：是.msh文件经过Fluent处理后得到的文件。 .dat文件：Fluent计算数据结果的数据文件。 三维定常速度场的计算实例操作步骤 对于三维管道的速度场的数值模拟，首先利用Gambit画出计算区域，并且对边界条件进行相应的指定，然后导出Mesh文件。接着，将Mesh文件导入到Fluent求解器中，再经过一些设置就得到形影的Case文件，再利用Fluent求解器进行求解。最后，可以将Fluent 求解的结果导入到Tecplot中，并对感兴趣的结果进行进一步的处理。
FLUENT 12 模拟步骤 Problem Setup 读入网格：file read case 选择网格文件（后缀为。Mesh） 1 General 1）Mesh（网格） > Check（点击查看网格的大致情况，如有无负体积等） Maximum volume (m3)（最大体积，不能为负） Minimum volume (m3)（最小体积，不能为负） Total volume (m3)（总体体积，不能为负） > Report Quality（点击报告网格质量） Maximum cell squish（最大单元压扁，如果该值等于1，表示得到了很坏的单元） Maximum cell skewness（最大单元扭曲，该值在0到1之间，0表示最好，1表示最坏） Maximum aspect ratio（最大长宽比，1表示最好） > Scale（点击缩放网格尺寸，FLUENT默认的单位是米） Mesh Was Create In（点选mm →点击Scale按钮且只能点击一次） View Length Unit In（点选mm →直接点击Close按钮不能再点击Scale按钮） > Display（点击显示网格设定）→弹出Mesh Colors窗口 Options（选Edges和Faces） Edge Type（点选All） Surface（点选曲面） →点击Display按钮 点击Colors按钮→弹出Mesh Display窗口 Options（点选Color by ID） →点击Close按钮→再点击Display按钮 2）Solver（求解器） > Pressure-Based（压力基，压力可变，用于低速不可压缩流动） > Density-Based（密度基，密度可变，用于高速可压缩流动） 3）Velocity Formulation（速度格式） > Absolute（绝对速度） > Relative（相对速度） 4）Time（时间） > Steady（稳态） > Transient（瞬态） 5）Units（点击设置变量单位） 点击按钮→弹出Set Units窗口→在Quantities项里点选pressure →在Units项里点选atm →点击
7.19多孔介质边界条件 多孔介质模型适用的范围非常广泛，包括填充床，过滤纸，多孔板，流量分配器，还有管群，管束系统。当使用这个模型的时候，多孔介质将运用于网格区域，流场中的压降将由输入的条件有关，见Section 7.19.2.同样也可以计算热传导，基于介质和流场热量守恒的假设，见Section 7.19.3. 通过一个薄膜后的已知速度/压力降低特性可以简化为一维多孔介质模型，简称为“多孔跳跃”。多孔跳跃模型被运用于一个面区域而不是网格区域，而且也可以代替完全多孔介质模型在任何可能的时候，因为它更加稳定而且能够很好地收敛。见Section 7.22. 7.19.1 多孔介质模型的限制和假设 多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流动阻力。本质上，多孔介质模型仅仅是在动量方程上叠加了一个动量源项。这种情况下，以下模型方面的假设和限制就可以很容易得到： ?因为没有表示多孔介质区域的实际存在的体，所以fluent默认是计算基于连续性方程的虚假速度。做为一个做精确的选项，你可以适用fluent 中的真是速度，见section7.19.7。 ?多孔介质对湍流流场的影响，是近似的。 ?当在移动坐标系中使用多孔介质模型的时候，fluent既有相对坐标系也可以使用绝对坐标系，当激活相对速度阻力方程。这将得到更精确的源项。 相关信息见section7.19.5和7.19.6。 ?当需要定义比热容的时候，必须是常数。 7.19.2 多孔介质模型动量方程 多孔介质模型的动量方程是在标准动量方程的后面加上动量方程源项。源项包含两个部分：粘性损失项（达西公式项，方程7.19-1右边第一项），和惯性损失项（方程7.19-1右边第二项） (7.19-1)
空调间内空气流场的模拟 模拟课题名称：空调间内空气流场的模拟 模拟的步骤： 1：gambit建立空间模型及建立网格 1.1打开gambit，选择文件的存储路径。 为了寻找的方便，我们将文件尊放在d盘fluent文件下 点击run，进入gambit界面； 1.2建立空调间的空间模型 点击体控制面板中的体控制按钮直接生成六面体，选择direction为+x+y+z方向输入房间尺寸，生成房间如图
同样的方法生成空调的六面体，并将六面体通过移动的选项将六面体移动到正确的位置上如图： 因为整个的计算区域为空调间减去空调的空间所以选择选项，将空调的空间去除，经过渲染之后得到下图： 这样计算区域的几何模型就建立起来了。 1.3网格的划分
这里我们直接对体进行网格划分点击，进入网格划分界面，选择该体，采用六面体的网格单元和submap网格划分方法进行划分如图 Interval size输入0.05，建立网格如图 如图，我们可以看到总共建立了359840个网格。接下来需要对网格单元记性检 查，经过检查发现网格符合要求。 1.4定义边界类型
本文定义房间的上下左右先后分别为壁面边界，空调的下边界为速度进口，空调的上边界为速度出口空调的其他面为壁面边界。建立的过程如图 建立完成后得到：
1.5保存文件后安全退出 执行File Export Mash，输出三维的数据文件，不要在2维的地方打勾，为 了以后查找的方便我们将文件的名字改为kongtiaojian，最后执行Exit安全退出。这里注意不要点右上角的，要选择exit键才能安全退出，否则生成的是 碎片文件。 2 fluent模拟 2.1文件的读取和检查 打开fluent三维单精度计算器如图，执行Flie Read Case，读入我们已经划分好的网格文件，当出现 表示网格文件读入完毕。执行Grid Check 检查网格，表示网格检查完毕。执行Grid Scale我们可
一收缩-扩张喷管实例1.1问题描述 本节内容主要依托收缩-扩张喷管内的流动计算展开。喷管外形如图1-1所示，A 为沿轴圆形截面面积，喷管的外形尺寸满足如下条件（单位：m）： 21.0x A +=5 .05.0<<-x 计算求解时可以将模型琪简化为二维轴对称问题，边界条件为：入口压力P m =101325Pa,入口总温T i0=300K，出口静压P 0=3738.9Pa。 图1-1喷管几何示意图 1.2创建几何模型 （1）设定工作目录 File→Change Working Dir，选择文件存储路径。 （2）创建Point，如图1-2所示。 Step 1通过输入坐标的方法创建P_1、P_2。选择Geometry 标签栏中的，单击，选择Create 1point（创建一个点），输入P_1的坐标，单击Apply 按钮确定，如图1-3所示。P_2创建方法与之相似，坐标为（0.5,0,0）。 Step 2创建点集3。因为横截面积为21.0x A +=，因此沿X 轴方向半径的函数为：5.02]/)1.0[()( x x R +=。单击，在Explicit Locations 下拉菜单中选择Create Multiple points，按照如图1-4所示输入数据，单击Apply 按钮确定。单击Apply 按钮确定。
图1-3创建P_1图1-4创建点集3 （3）创建Curve，如图1-5所示。 图1-5创建Curve结果图 Step1选择Geometry标签栏，单击创建Curve。如图1-6所示，单击，再单击，依次选择点集3中的各点连成曲线，创建C_4。 Step2采用Step1的方法创建其余三条Curve。 （4）定义Part。 ICEM中定义Part的名称将会是导出网格后边界的名称，可以简化在求解器