离心压气机转子数值模拟培训教程

1、离心压气机转子数值模拟培训教程,NUMECA FINE/TurboTM 7.4_1,目 录,基础知识 离心压气机简介 CFD技术简介 FINE/Turbo软件简介 增压器离心压气机转子数值模拟 网格生成 边界设定及求解 计算结果显示,离心压气机简介,增压器中压气机的任务是吸收涡轮的机械功,以提高空气的压力,然后经发动机进气管将空气导入发动机,以达到增压目的. 离心压气机由进气道、工作轮、无叶（叶片）扩压器、蜗壳组成。 离心压气机的主要参数有：转速、压比、效率和流量等。 离心压气机内部流动的一般特征： 三维，有粘，非定常 可压缩，跨音速 多相流动，化学反应 传热传质，流热耦合，流固耦合，气动噪声

2、,CFD技术简介,计算流体动力学（CFD）是二十世纪六十年代开始、随着电子计算机的出现和应用而得以蓬勃发展的一门新兴学科。 CFD是在经典流体力学、数值计算方法和计算机技术的基础上建立起来的。 CFD技术弥补了理论分析和实验研究的不足，并促进了理论和实验的发展。 应用CFD技术，在分析、掌握离心压气机内部流动的基础上，对其通流部件气动性能进行优化，是提高产品竞争力的一个有效方法。,CFD分析步骤,采用数值方法对具体的流动现象进行模拟的过程主要分为以下几个步骤： 根据具体流动情况和所关心的问题，建立或选择对应的数学模型； 如欧拉方程，雷诺平均的NS方程，湍流模型等 根据所选择的模型，对计算区域进

3、行网格剖分； 选择对应的流体属性； 如理想气体、实际气体、水 采用合理的边界条件； 根据具体条件或实验，设定进、出口条件、固壁边界 流场求解； 选择数值方法、收敛准则等 计算结果处理和显示。 制做各种定性、定量物理量图线,CFD分析步骤,1. 网格生成,工业问题,2. 给定流体特性 3. 给定边界条件,4. 流场求解,5. 数据处理及演示,Backflow,CFD求解器数学模型,DNS：在非定常条件下求解NS方程及所有尺度的湍流涡旋网格非常密（约为RANS的105倍） LES：在非定常条件下求解“过滤”后的NS方程及最大的涡旋而模拟小的涡旋网格较密（约为RANS的100倍） RANS：求解“平

4、均的”NS方程而模拟所有的湍流涡旋，网格的数量级为几十万到几百万 （可在较好的PC机上实现）,势流、（两类相对流面理论）、无粘/粘性流、定常/非定常、雷诺平均NS方程+湍流模型（RANS），大涡模拟（LES），NS方程直接解（DNS）,Fine/Turbo软件组成,IGG 交互式分区结构化网格生成器 可高质量地生成任意复杂区域的网格 AutoGrid 柱形回转面结构(如搅拌器及透平机械)自动网格生成器 可以方便、快速地生成包括离心、轴流压缩机、泵、风机、透平、燃气轮机、航空发动机、增压器、内燃机等多种对象的旋转机械叶片网格 EURANUS 全二阶精度结构化块网格求解器 基于密度求解方法保证了快

5、速的流场收敛速度 采用多重网格技术，显著提高计算速度 提供了针对透平机械专用的边界及初场给定方法 CFView 功能强大的流场显示后处理器 提供定性、定量参数分析，并提供针对透平机械专用的参数显示功能（如子午平均、叶片叶片）,培训内容,带分流叶片离心压气机叶轮数值模拟，主要包括： 网格生成：AutoGridTM 计算设定和求解：FINETM 结果处理：CFViewTM,准备工作,用三维绘图软件PRO/E、UG、SOLIDWORKS.等绘制要计算的叶轮三维图或子午面图，并将这些图存为.igs或.iges格式的文件。 没有三维绘图功能的CAD无法绘制或打开.igs文件，要通过proe、UG等转换成

6、DWG格式才能打开。Fine只能打开.igs或.iges格式的文件。Iges（initial graphics exchange specification）是业界通用三维数据（点，线，面等）的一种交换标准。但转换时有时容易破面。,准备工作,2. 打开FINE界面(开始程序NUMECA software fine74_1 FINE) 3. 选择创建一个新的工程文件（Create a New Project），命名为compressor.iec。 4. 在弹出的网格选择对话框中，选择Create Grid File来打开IGG窗口。,第一部分 网格生成,带分流叶片离心压气机网格生成 操作及示例,

7、网格生成步骤,Step1：准备几何数据,几何数据的准备有两种方法，一种是通过AutoGrid专用的几何文件格式准备为GeomTurbo文件，另外一种方法是通过IGG从三维模型中提取所需的几何型线及面。一般工程中都是从图纸开始准备，因此这里着重介绍第二种方法。,基本思路：,在CAD作出的图形中，叶片型面可能是由很多个小面组成，组成面的弧线通常也会由若干段弧线拼成。而用在AutoGrid中的几何，希望叶片不同截面(根部、顶部)的型线控制点数尽可能相同或相似，这样可以减小或消除几何上的误差。因此，通常需要对已有的型线进行重新描述。常用的方法是贴网格面，通过控制网格点数及分布来控制几何线的控制点数及分

8、布。,由于叶片通道是周期性的，每个流道内的流动情况都相同(只计算叶轮通道，不考虑下游涡壳不对称性对叶轮内流动的影响)，因此计算的时候只需计及一个通道即可。故可选任一长叶片及其相邻的短叶片作为取型及计算对象。,Step1：准备几何数据,将准备好的叶轮文件（.igs格式）输入IGG， (File Import IGES DATA)，然后关闭输入窗口(CLOSE)。 IGG中每条线都对应有一个名称。如果原始的名称相同，在IGG中调入后，会弹出一个对话框，此时只要选择“Auto Rename”即可,IGG会对其重新命名。再选择“Apply to all”，,Step1：准备几何数据,点击右下脚显示按钮

9、，输入的图形就会显示出来，如图。右下角有一排按钮，可以对显示的图形进行旋转、放大等操作。,注意：AutoGrid中，叶轮的旋转轴必须是Z轴，如果在CAD作图时选的轴不是Z轴，需要进行旋转操作，以符合AutoGrid的要求。,Step1：准备几何数据,重新打开一个只画有一片叶片的文件，如图。点击屏幕左侧上方快捷面板中的“Insert New Face”。做出一个网格面。按图示顺序确定该网格面的两个顶点。刚画完时按ESC键可以取消。,附：网格的有关名词,每个网格块由六个面组成。 一个面（Face）可以由一个或多个片（patch）组成。 一个面由四个边（edge）组成。 一个边可以由一个或多个线段（

10、Segment）,一个网格块内可以插入若干个内部网格面 一个网格面内可以插入若干个内部网格线,Step1：准备几何数据,将网格面的顶点拖至长叶片压力面的四个定点上。点一下网格面，线会变粗，顶点会显示为方块。点一下顶点，松开左键，顶点就会随着鼠标移动，拖至合适的地方，点一下即可。此时不满意时怎样取消该面？随时保存 如果叶片沿某一方向是由一条线确定的，那么网格面的边会自动贴附到几何线上。但是在本算例中，叶片绘图时是由几段线段组成，因此，需要在每一个线段连接点处插入点来使网格边和几何边贴合。选择屏幕左上方快捷面板中的”Insert vertex”(插入顶点)，会有一个顶点出现，将其放至网格面的边上，

11、点一下，然后松开左键，这个点就会随着鼠标移动，当其移动到线段连接点时，点一下松开即可。为了方便，可以把图形放大再操作。加的位置？顶点的数目？数目对划分网格和计算结果有否影响？对边加的数目是否需要一致并对称？ 贴好的网格面如图所示。,Step1：准备几何数据,Step1：准备几何数据,调整网格数目（目的：控制几何线的描述）： 左键点击选中需要调整的网格边（该边变为黄色），右键选择SegmentSet Number of Points,在弹出对话框中将网格数目设置为61，点击Apply后关闭对话框。 用同样的方法调整另一条边数目为21。 这样，该长叶片压力面就由一个6121的网格面表述出来。,St

12、ep1：准备几何数据,查看生成的网格面：快捷面板 View/Grid/Toggle Face Grid 投影网格面至原始几何面（目的：使网格面和原始几何面完全贴合） 9.1 用选面模式选中原始叶轮压力面（由若干个面组成，选中后面的四条边显示为红色） 9.2 点击快捷面板中的“Project Face”（Generate下面） 9.3 点击Apply确定,9.3,附：如何激活选面模式,在IGG中，面的选择模式可以通过以下几种方法激活： 下拉式菜单GeometrySelectSurfaces 快捷键“Ctrl s”（IGG中区分大小写，在按下s之前，请确认CAPSLOCK健未被激活） 在IGG窗口

13、的空白处，点击鼠标右键，弹出式菜单中点击“Select Surfaces” 当选面模式被激活时，IGG窗口下侧会出现黄色提示文字： 在选面模式下，用户可以通过鼠标点选相应的面。当鼠标移动至该面的位置时，面会呈现出蓝色的网状面。鼠标点击后，面的轮廓呈红色则被选中，面呈青色则未被选中。 在选面模式下，用户可以连续选择/取消选择多个面。选面完毕后，点击鼠标右键释放选面功能，则IGG底部的提示文字消失。 在选面模式下，按下小写字母“a”可全部选择或全部取消选择所有显示的面。每按一次“a”，则全部面的被选状态全部发生改变。,Step1：准备几何数据,用贴好的网格面生成几何面：在左侧的快捷面板中，点击Fr

14、om Face Grid，则会生成一个和网格面重合的几何面。 取消网格面及边的显示。,Step1：准备几何数据,在选面模式（ctrls）下，按下小写字母a，释放所有选中的面，然后鼠标点击选中刚才生成的面（可在该面位置用鼠标中键切换不同的面）。 用GeometryModify Surface Representation（画像，表示法），将选中的面沿叶高方向描述为5条曲线。在下方的输入栏中输入（5 2），回车确定。如果这时分不清哪个是叶高方向，可试着输入，如果得到的不是沿叶高方向分布的5条线，可继续执行该操作，将输入栏中的数字顺序颠倒即可）。 取消面的选择，通过右键激活选线模式（Select C

15、urves），从叶根至叶顶，依次选中这5条线，使用Geometry Create Surface lofted命令，生成loft形式的面。,Step1：准备几何数据,通过Geometry View Surfaces，显示刚才生成的面（lofted_2） 通过选面模式选中该面，用FileExportGeometry Selection将选中的面输出为PS1.dat,作为长叶片的压力面。 用同样的方法将长叶片的吸力面及分流叶片的压力面、吸力面分别输出为SS1.dat，Splitter_PS.dat, Splitter_SS.dat。（每次选择面前，需要按释放先前所选的面。）,叶片的取型到此完成。下

16、面需要取出子午面的型线。,通过通过Geometry View Surfaces，显示所有几何面 用同样的方法将长叶片的吸力面及分流叶片的压力面、吸力面分别输出为SS1.dat，Splitter_PS.dat, Splitter_SS.dat。（每次选择面前，需要按释放先前所选的面。）,Step1：准备几何数据,重新打开一个IGG工程（ File New ），输入子午面文件（.IGES格式）。,由于CAD中的端壁线不完整，需要对HUB（中心）和Shroud（覆盖物）进行修补。具体的方法是沿原始端壁的走向用C样条曲线做出所需的端壁线。,Step1：准备几何数据,点击左侧快捷面板Geometry C

17、urve CSpline，按照图示的起点开始画起，通过若干个中间控制点，画出一条和原始Hub线重合的线条（黄色所示）。 通过 FileExportGeometry Control Points将该线输出为hub.dat,作为轮毂线（默认处于激活状态）,HUB,Step1：准备几何数据,通过选线模式取消HUB线的选择。 用同样的方法画出Shroud线。 通过 FileExportGeometry Control Points将该线输出为shroud.dat,作为轮盖线,SHROUD,附：如何激活选线模式,在IGG中，线的选择模式可以通过以下几种方法激活： 下拉式菜单GeometrySelectC

18、urves 快捷键“s”（IGG中区分大小写，在按下s之前，请确认CAPSLOCK健未被激活） 在IGG窗口的空白处，点击鼠标右键，弹出式菜单中点击“Select Curves” 当选线模式被激活时，IGG窗口下侧会出现黄色提示文字： 在选线模式下，用户可以通过鼠标点选相应的线。此时，在线上点击鼠标后，如果线的颜色为黄色，则代表该线被选中。如果线为青色，则代表该线未被选择。 在选线模式下，用户可以连续选择/取消选择多条线。选线完毕后，点击鼠标右键释放选线功能，则IGG底部的提示文字消失。 在选线模式下，按下小写字母“a”可全部选择或全部取消选择所有显示的线。每按一次“a”，则全部线的被选状态全

19、部发生改变。,Step2：导入几何参数,从IGG切换至AutoGrid模块,启动AutoGrid后，会出现一个对话框，记录了以前保存过的文件。用户可以直接选择打开。如果列表中没有所需的文件，可以通过按钮2选择打开。如果是第一次启动AutoGrid，那么这个区域则是空白的。,26.1 关闭该对话框。,2,Step2：导入几何参数,AutoGrid 5 主界面：,Step2：导入几何参数,快捷面板：,Step2：导入几何参数,在左侧快捷面板Geometry Definition中打开Import CAD窗口。用户可以在这个窗口中定义Hub，Shroud，Blade（Pressure Side、Su

20、ction Side）以及子午效应、3D效应等。,打开先前生成的hub.dat, shroud.dat, PS.dat, SS.dat, Splitter_PS.dat, Splitter_SS.dat。 （File Open ）,在叶片排参数定义过程中，必须保证该叶片排处于激活状态，即显示为图示红色！,Step2：导入几何参数,定义HUB，SHROUD。,29.1 IGG中的线段是有方向的，要使AutoGrid能正确调用，每一条线的方向必须是从进口指向出口，从叶根指向叶顶。选中所需要的线（Geometry Select Curves,黄色为选中，青色为未选中），查看线段的方向（+o），如果发

21、现有方向不对的线，用Geometry Modify Curve Reverse工具将该线反向。本例中在画HUB和SHROUD线时都是从进口指向出口，因此不存在反向问题。,Step2：导入几何参数,29.2 左键点击任意空白处，释放两条线的选择。 29.3 点击选中Shroud线（变为黄色），点击鼠标右键，在弹出菜单中选择Link to Shroud 。在AutoGrid主界面子午视图中会出现Shroud 的型线。,注意：Shroud（Hub）也可用多条线段定义，但在选择指定的时候，选择的顺序必须是从进口到出口！,Step2：导入几何参数,29.4 点击任意空白处释放Shroud线的选择 29.

22、5 点击选中Hub线，在右键弹出菜单中，选择Link to Hub。,Hub和Shroud定义完后，AutoGrid会自动将两条线的端点连接，并将所围成的区域设定为计算域，以蓝色表示。,Step2：导入几何参数,定义长叶片。 30.1 将鼠标移近长叶片压力面，当出现蓝色网格面的时候，左键点击。这时该面的四条变为黄色，表示被选中。点击鼠标右键，在弹出菜单中选择Link to Pressure Side。这一操作将在AutoGrid主界面子午视图中定义出叶片的压力面。定义完成后，左键点击任意空白处，面的四条边恢复成青色，表示该面的选择被释放。,30.2 用同样的方法定义长叶片吸力面(Link to

23、 Suction Side)。,注意：只有在叶片前缘和尾缘都是钝头的情况下需要分别link pressure side和suction side，只要有一个是圆头，则只需同时选中吸力面和压力面，link to blade！,Step2：导入几何参数,定义前缘、尾缘位置。,31.1 选择长叶片进口处任意一条线（ Geometry Select Curves ），右键选择Link to Leading Edge,31.2 用同样的方法选择长叶片出口处任意一条线，右键选择Link to Trailing Edge,到此为止，大叶片的定义全部完成。由于该叶轮是半开式，叶片顶部和轮盖间存在间隙，故叶片和

24、轮盖并不相交。这时，AutoGrid会弹出一个右图所示的提示框。关闭即可。,Step2：导入几何参数,定义分流叶片。,32.1 在左侧快捷面板中，左键激活row 1，右键打开弹出菜单。选择Add Blade。 32.2 激活splitter 1，按照长叶片的方法，在Import CAD窗口中定义分流叶片的吸力面、压力面、叶片前缘、尾缘。,同大叶片一样，在小叶片定义完后，会弹出一个提示框，提示叶片和端壁没有相交。,定义几何单位。,33.1 在左侧快捷面板几何定义菜单下，将几何单位设置为Millimeters,Step2：导入几何参数,叶片延伸,34.1 激活Main Blade，在右键弹出菜单中

25、选择Expand Geometry。 34.2 在Hub栏下选择expand，设置Cut offset=0.5，Extent offset=0.5。表示将Hub线先向内移动0.5mm，保证Hub和叶片有交线；再向外移动0.5mm，将叶片延伸至该处。 34.3 在Shroud栏下，分别设置Cut offset1，Extent offset1。 34.4 点击Apply确定。 34.5 激活splitter1，用同样的参数将小叶片延伸。,AutoGrid在生成网格时，要求叶片和端壁必须有交线，否则网格不能正常生成。因此，必须对叶片进行延伸。这一操作只是为了满足网格生成的需要，延伸出的部分会被端壁及

26、间隙所切割，不会改变真实几何。,Step2：导入几何参数,所有的几何定义完全后，在计算域的前后会出现两条线，表示计算域的进口位置和出口位置。如果定义不完全，那么这两条线将不会出现，因此可通过这种方法检验几何定义是否完全以及正确。,Step2：导入几何参数,关闭Import CAD窗口。 按照图示调整进出口位置。,默认的进出口位置是hub和shroud两端点的连线，用户可根据需要或实际的物理情况对其位置及形状进行调整。,控制点控制方法： 鼠标拖动: 将鼠标移至进(出)口上,当进(出)口线变为黄色时,左键确定.这时会出现若干个控制点,鼠标移至控制点时,该点会变成黄色,左键点击该点后,即可自由拖动。

27、 精确定位：当控制点变成黄色后，点击鼠标右键，则会弹出文本框，用户可输入其（R,Z）坐标进行定位。 当进(出)口线变为黄色时，点击鼠标右键，在弹出菜单中选择Properties，在Shape栏中定义所需位置及形状。,Step2：定义叶片参数,左键激活row1（变为红色），在右键单击弹出的菜单中选择Properties。 在弹出的对话框中，输入 Periodicity（叶片数）6，Rotational Speed（转速）60000 rpm， 选择Impeller，Centrifugal。 关闭该对话框。,NUMECA中，转速有正负之分。可用右手法则来判断。四指的旋向和叶轮转动方向相同，大拇指的指

28、向和Z轴正方向相同，则转速为正，反之则为负。,定义叶顶间隙。,40.1 激活row1，在右键弹出的菜单中选择Define Shroud Gap。 40.2 在弹出的对话框中，按下图数值定义叶顶间隙尺寸。每输入一个数值，必须按确定。 40.3 关闭该对话框。,间隙设定： 间隙既可以是根部间隙，也可以是顶部间隙，分别对应Hub Gap和Shroud Gap. 间隙的尺寸单位和叶型数据单位相同。用户在进行其他算例的网格生成时，需额外注意叶型尺寸的单位，这样方可确定应当给定的间隙尺寸的单位是米、厘米还是毫米。 间隙中展向网格的数目包含在展向总网格数中。 AutoGrid中支持三种方式定义间隙： a.间

29、隙尺寸不变。此时只需指定前缘和尾缘的间隙尺寸为相同值即可。 b.间隙尺寸线性变化。此时需分别指定前缘和尾缘的间隙尺寸，AutoGrid会根据这两个尺寸值进行插值确定间隙尺寸的变化规律。 c.间隙尺寸是随意的变化形式。此时可通过Defined Shape来引入定义间隙根部的型线即可。,Step2：定义叶片参数,Step3：生成默认拓扑结构,在左侧快捷面板Mesh Control下，选择Coarse。（一般来说，为了保证数值模拟的精确性，网格越密越好，通常都选择Medium作为基准网格。本例只做演示使用，为了节省时间，故选择了最粗糙的网格Coarse） 计算壁面第一层网格尺度,ywall为壁面第一

30、层网格大小，单位m； Vref为参考速度，单位m/s； Lref为参考长度，单位m； 为流体的运动粘性，单位m2/s; Y+为无量纲量，对应不同的湍流模型，有不同 的取值范围。,这里，Vref取叶轮出口速度177.5 m/s；Lref取叶轮半径0.02825 m；取1.4310-5m2/s; Y+取5，计算得到ywall110-5 m。由于本例中的几何单位是mm，故这里需要将单位转换为mm。因此，最后得到的第一层网格尺度为0.01mm。,42.1 将Cell Width设置为0.01，确定,Step3：生成默认拓扑结构,点击（Re）set Default Topology,AutoGrid会根

31、据选择的网格层数及第一层网格大小，自动生成一个网格拓扑结构。一般来说，这个默认的网格拓扑结构都可以获得一个比较满意的网格。当然，用户也可在此基础上根据需要对网格进行调整。调整的内容包括两个方面，一个是子午面，一个是叶片叶片（B2B）。,附：网格拓扑结构介绍（一）,H网格,H-I网格,H-O-H网格,附：网格拓扑结构介绍（二）,O 网格,嵌套网格,附：网格拓扑结构介绍（三）,H网格,蝶型网格,叶顶间隙,附：网格拓扑结构介绍（四）,Skin网格,子午面控制包括调整叶高方向网格数目，分布，以及间隙内的径向网格数目及分布。这里可不进行任何调整，使用默认拓扑结构设定的参数。,Step4：子午面控制,St

32、ep5：叶片叶片控制,叶片叶片控制包括调整网格拓扑形式、周向、流向网格数目以及间隙内O网格数目。这里也不做任何修改，使用默认拓扑结构设定的参数。,Step6：生成三维网格,点击工具栏中的Generate 3D，在弹出的窗口中点击yes确认三维网格生成。,最终生成的网格，不能存在负网格，否则计算无法进行。 原则上，最小正交性角度越接近90越好；最大网格长宽比越接近1越好；最大网格延展比越接近1越好。但实际中，很难得到三者兼得得网格质量，所以一般推荐： 最小网格正交性角度5 最大网格长宽比5000 最大网格延展比10 应当尽可能地提高网格质量，以避免不必要的数值误差。 对于网格正交性而言，在某些算

33、例中，可能甚至会出现角度小于1的情况。但一般而言，只要不存在负网格，计算就可以进行。,Step7：检查网格质量,三维网格生成后，会自动弹出一个网格质量检测报告。通过这个报告，用户可以粗略地了解所生成网格的最小正交性，最大长宽、最大延展比及其分布。,附：网格质量参数定义,正交性（Orthogonality）：网格相邻两条边之间相互垂直的程度。二维标准，范围0-90。正交性表征网格面两条边之间的最小角度。如果两条边的夹角大于90，则正交性按（180-实际角度）确定。 长宽比（Aspect Ratio)：如图示。 二维标准，范围1-10000。如果实际值超出这 个范围，重置为10000。 延展比（E

34、xpansion Ratio)：如图示。 三维标准，范围 1-100。表征相邻两个网格 的大小变化。和方向有关。如果实际值超出 这个范围，重置为100。 如果在某个方向只 有一个网格，那么这个参数没有任何意义。,Step7：检查网格质量,除了质量报告外，AutoGrid还提供了详细的网格检测工具，包括各个block的正交性，长宽比，延展比，以及边界条件，负网格检测等。用户可根据需要选取相应的菜单命令。,查看各个block中的网格质量。,在Type下拉菜单中可以选择想要查看的量，并在Block下选择相应的块，回车确定后，可用该页面下的Show chart和More info查看详细信息。,Ste

35、p7：检查网格质量,查看边界条件（Boundary Conditions）,三维网格生成结束后,AutoGrid会自动设定边界类型,其中包括:进出口边界(INL,OUT),周期性边界(PER/PERM)，连接边界(CON/NMB)，固壁边界(SOL)等.,在IGG/AutoGrid中,支持用户人工设定的边界条件类型包含: UND:未定义条件(Undefined); SNG:奇线边界条件(Singular); EXT:远场边界条件(External),用于外流; SOL:固壁边界条件(Solid); INL:进口边界条件(Inlet); OUT:出口边界条件(Outlet); ROT:转静子交接

36、面边界(Rotor/Stator); MIR:镜像边界条件(Mirror); FNMB: 完全非匹配边界条件(Full Non-matching Boundary);,在IGG/AutoGrid中,可自动搜索出的边界类型包含: CON:匹配连接边界(Connection); NMB:非匹配连接边界(Non Matching Boundary); PER:匹配周期边界(Periodicity); PERNM:非匹配周期边界(Periodicity Non-Matching ).,附：边界条件连接边界,直接连接（CON）,非匹配连接（NMB）,全非匹配连接（FNMB）,附：边界条件周期边界,附：边

37、界条件对称边界,对称边界（MIRror 边界） 几何对称和流场对称 例如直列叶栅的中叶展面,对称面,附：边界条件远场边界,远场边界（EXT边界） 主要用于外流分析 绕弹体网格外边界均为远场边界,附：边界条件奇面边界,奇面边界（SNG 边界） 一个边界面退化成一条线,奇面边界,附：转/静子交界面,转/静子交界面（ROT边界）,Step8：保存网格及模板,File Save project,这一操作可将网格及模板同时保存。在用户目录中，生成12个文件，包含了所有的网格信息。主要的几个有： *.geomTurbo：几何模板，记录了所有的几何信息。 *.trb：网格模板文件，记录了所有的网格控制参数，

38、用户在下次需要生成网格的时候可以直接打开该文件。 *.igg：网格文件，为FINE计算所需。,如果不再对网格质量进行一步的改进，即可跳过下面的网格调整部分，返回FINE界面，进行计算设定。,网格调整：子午面控制,调整展向网格数目及分布,48.1 将Flow Paths Number调整为49。这49个节点是展向的总网格点，包括间隙内网格点。 48.2 按照右图数目调整展向网格分布。（Flow Path Control） 48.3 点击Generate后，关闭该对话框。,调整间隙内网格数目及分布,49.1 在Main Blade下激活Tip gap，在右键弹出菜单中选择Properties。在弹

39、出对话框中设定长叶片间隙内网格数目及分布。,49.2 激活splitter1下的间隙，按照上述方法设定短叶片间隙内的网格数目及分布（参数同长叶片）。,49.3 点击Generate Flow Paths后关闭该对话框。,网格调整：叶片叶片控制,在row1下激活Main blade。以下的操作都是针对长叶片。打开B2B Mesh Topology Control窗口。,50.1 打开Topology项。,这一页面可以设置匹配/非匹配周期边界及进出口形式。对于本算例，AutoGrid自动生成了匹配的周期边界（matching）。为了改善网格质量，将其改为非匹配（non-matching）。进出口形

40、式都为Normal。用户可根据需要调整。,50.2 打开Grid Points项。,按照右图相应的数目调整网格数量。为了保证计算时可以采用至少三层的多重网格，网格数目必须为2n+1（n2）。,50.3 激活Mesh项。,将间隙内的网格数（Number of Point in Clearance O-mesh）设置为5。,50.4 关闭该窗口。,在row1下激活Splitter1。打开B2B Mesh Topology Control窗口。用同样的步骤按右图参数调整短叶片的网格。,在工具栏中点击Generate B2B按钮，更新B2B网格。,网格调整：叶片叶片控制,NUMECA为了加速收敛，采用

41、了多重网格技术（详细介绍参见FINE计算设定中数值模型，即Numerical Model部分）。而多重网格的层数是通过网格数目来确定的。因此，为了满足计算时能够采用多重网格，在设置网格数目时就要满足一定的要求。,如：17241，min（n）4，即满足5重多重网格。 61252423221，min（n）2，即满足3重多重网格。,附：网格数目调整原则,（n 2）,多重网格的层数为：min（n）1,在IGG/AutoGrid中，用户可以方便地通过网格输入框右侧的箭头选择网格数目，以保证其符合多重网格的要求。,网格调整：叶片叶片控制,调整网格光顺步数。,光顺步数越大，网格过渡越圆滑，但是网格生成所需的

42、时间也越长。 用户可尝试改变光顺步数（0，50），更新B2B网格后观察网格的变化。 本例使用默认值100。,重新生成三维网格。,可以看到经过调整，网格质量有所提高。,重新保存网格。,至此，网格生成部分已全部完成，下面进行计算设定。,第二部分：FINETM,带分流叶片离心压气机计算设定 操作及示例,FINE/Turbo主界面,从AutoGrid中切换至FINE界面：,You are now back to FINE graphical user interface！,计算设定步骤,Step1：链接网格,点击选择网格文件按钮，指定前面生成的网格文件（compressor.igg）。 设定网格属性（

43、Millimeter，Cylindrical，3 Dimensions）。,如果在设定完成后，需要变更网格属性，可通过菜单Mesh Properties修改。,如果需要变更网格文件，可通过文件夹按钮来指定。,Step2：创建或更改计算名称,2.1 为计算更名： 点击计算名（变为蓝色高亮度框） 单击鼠标右键，选择Rename或者在顶部菜单中单击Rename 输入新的计算名并回车确定,2.2 创建或复制当前计算： 点击计算名（变为蓝色高亮度框） 单击鼠标右键，选择New-Duplicate或者在顶部菜单中单击New 输入新的计算名并回车确定,2.3 删除当前计算： 点击计算名（变为蓝色高亮度框）

44、单击鼠标右键，选择Remove或者在顶部菜单中单击Remove OK确认。,该操作可得到一个与原有计算完全相同的计算，对于进行变工况非常适用，用户只需要改变几个参数即可，其他完全不用改变。,Step3：选择工作介质,Parameters Configuration Fluid Model (双击),3.1 选择AIR（Perfect）作为工作介质（单击，变为蓝色高亮度框为选中）,用户可以在已有的流体列表中选择相应的工作介质。如果列表中没有需要的流体，可以通过Add New Fluid来定义新的流体。,Step4：设定流动模型,Parameters Configuration Flow Mode

45、l (双击),4.1 选择Steady选项（定常计算） 4.2 选择Turbulent Navier-Stokes方程组及Spalart-Allmaras湍流模型 4.3 定义参考长度：Length=0.02825m, Velocity=177.5m/s, Density=1.2kg/m3, Temperature=293K, Pressure=101300Pa.,参考长度是用来计算雷诺数（Reynolds）的，对收敛性和计算结果等没有任何影响。为了简便起见，这里的取值和网格生成时计算y+所用的值保持一致。,对层流流动： 未知量为：密度r、三个速度分量和温度。 压力可由状态方程求出。理想气体的

46、状态方程为,附：湍流模型,控制方程,六个未知数、六个方程，方程封闭,对湍流流动： 假设任何流动参数其时均值其脉动值 对NS方程进行时均处理后 RANS方程新的未知量:雷诺应力张量,方程不再封闭,采用湍流模型封闭方程,零方程模型 BaldwinLomax（简称BL模型） 特点：简单易用，可用于附着流动和小的分离流动。往往过度模拟分离流动。,附：湍流模型,一方程模型 Spalart-Allmaras模型（简称SA模型） 特点：加入一个控制方程，与BL模型比较，求解速度慢、空间需求大，计算精度高。可模拟流动分离和边界层转捩。,两方程模型 k-模型 特点：加入两个控制方程，与SA模型比较，求解速度慢、

47、空间需求大；计算精度近似，可模拟小尺度流动分离。可得到湍动能和耗散率的分布。,Step5：设定转动部件,Parameters Configuration Rotating Machinery (双击),在这一页面中，FINE根据AutoGrid中定义的叶片排，自动将block组合，并按照AutoGrid/Properties中设定的转速，自动设置Rotational Speed。因此，这里不需要用户做任何修改。,Step6：设定边界条件（BC）,Parameters Boundary Conditions (双击),6.1 打开网格 Mesh/View On/Off 6.2 拖动网格区域左上角

48、的按钮，同时显示边界设定区域和网格。,附：进口边界给定方式,亚音速 （绝对马赫数的子午分量1） 要给定所有的未知量的值,附：出口边界给定方式,亚音速 （绝对马赫数的子午分量1） 不用给定任何条件 。,Step6：设定边界条件（BC）,6.3 进口边界条件（INLET）,选择INLET，选中所有patch，将其组合为Inlet。 给定进口边界： 轴向进气 总压101325 Pa 总温293 K 进口湍流粘性3e-5 m2/s,Step6：设定边界条件（BC）,6.4 出口边界条件（OUTLET）,选择OUTLET，选中所有patch，将其组合为Outlet 给定出口边界条件 平均静压105000

49、 Pa,Step6：设定边界条件（BC）,6.3 周期性边界条件（PERIODIC） 无需做任何更改,Step6：设定边界条件（BC）,6.4 固壁边界条件（SOLID）,Fine会根据AutoGrid中设定的叶片参数，自动将转动部分（Hub&Blade）组合，并设置好转速，因此用户如果没有特殊要求，该页面也不需要设定。,选择SOLID，查看其设定。 将其余的patch组合为Shroud，使用默认的设置。,取消网格显示 Mesh/View On/Off,Step7：设定数值参数,Parameters Numerical Model (双击),这一页面中的参数可使用默认值，不需要改变。,附：CF

50、L数,CFL (Courant-Friedrich-Levy)数或库朗特数 是一个控制时间步长的参数。CFL越大，时间步长就大，计算的时间就短。 但由于稳定性的限制，CFL不能很大。一般在0.510的范围之内。常用值为24。 对于复杂问题，CFL要小；接近喘振/失速工况，CFL要小。 细网格计算的CFL要比粗网格时小。,近似的：,附：多重网格法,一种加速收敛的方法,细网格 0,粗网格 1,粗网格 2,对于三维问题： 网格000：在I， J， K方向的网格点数为所用网格的最大点数（如33，65，49） 网格111：在I， J， K方向的网格点数分别比网格000时约少一半（如17，33，25） 网

51、格222：在I， J， K方向的网格点数分别比网格111时约少一半（如9，17，13） 对于网格012，其网格点数：33，33，13,附：多重网格法,V形多重网格循环,在从细网格0粗网格1粗网格2 粗网格2 粗网格1细网格的工程中，在每个网格层上都要进行若干从RungeKutta计算或称为sweep（扫掠）。 缺省的sweep次数为：网格0：1次；网格1：2次；网格3：4次； 。 为了提高计算速度和计算的稳定性，也可采用：网格0：1次；网格1：4次；网格2：8次； 。,附：全多重网格（Full Multi Grid）,第3层,第2层,V-Cycle 全多重网格示意图,第1层,全多重网格循环,第

52、0层,N次单一网格计算,N次三重网格计算,在粗网格上计算，速度快。 残差降到-3左右即可,三层多重网格,Step7：设定数值参数,在窗口右上角的下拉框中选择专家模式（Expert Mode）,为了加速收敛，可增大在粗网格上的循环次数，以得到较好的初场。,7.1 在Scheme definition中，选择User Defined。 7.2 将Number of sweeps on successive grid levels中设为“1 4 8”。,增大粗网格上的循环次数可增强收敛稳定性，并显著提高收敛速度。但是循环次数越大，每一步迭代所需的时间也越长。,Step8：设定求解初场,Paramet

53、ers Initial Solution (双击),For turbomachinery是特别针对叶轮机械的一种初场给定方式。这里只需要给定进口的压力（如果有转静子面的话，还需给出转静子面处的压力）。,8.1 设定Estimated static pressure90000Pa,初场对最终的收敛结果没有影响。但一个较好的初场可以大大提高收敛速度，而一个不合理的初场很可能导致计算在一开始就发散。,Step9：输出变量设定,Parameters Outputs (双击),输出变量中有几个量是必选的：静温、静压、密度、相对速度。其他的热力学参数都可以在CFView中通过这几个量计算得到，无需在FIN

54、E中输出，这样可以节省存储空间及时间。为了简化设定，FINE中已将这几个量自动选定。 其他几个子页面中，常用的几个量也已作为默认设定勾选，因此如没有特别要求，这里不需要进行设定。,Step10：控制参数设定,Parameters Computation Steering Control Variables (双击),10.3 Save Solution Every (100) Iterations:中间结果保存设定 10.4 Minimum output：最小输出选项（不激活） 10.5 Set the required memory:内存需求设定(无需更改) 10.6 专家参数（无需更改）,

55、到此为止，所有的计算设定已全部完成！,10.7 保存文件（FileSave）,10.1 Maximum Number of Iterations：在细网格上的迭代次数(500) 10.2 Convergence Criteria:收敛量级（- 6）,Step11：开始计算,11.1 在菜单中选择SolverStart,或在图标栏中点击开始图标，开始计算。,11.2 双击Parameters Computation Steering Convergence History，观察收敛历史。,用户也可在Monitor中查看详细的收敛史，包括每一个block的残差收敛史，进出口流量收敛史，压比收敛史等

56、等。关于Monitor详细的说明请参见FINE用户手册。 12.1 在Windows下：开始程序NUMECA software fine74_1 Monitor 12.2 通过Add，加载当前计算的res文件。,Step12：查看收敛史,附：计算残差定义,计算残差 收敛后为计算残差值：,理想值,实际可达到值,全局残差：建议全局残差下降三个量级以上。 各块中残差：由于各块中网格质量以及流动特性的不同，每一块中的残差下降幅度也会不同。但仍然推荐计算中,每一块中的残差下降三个量级以上。 进出口流量：收敛准则中最重要的一个参数.建议进出口流量相对误差小于0.5%，且流量步再发生变化。对于有大分离涡的流

57、动(尤其在进出口处)，流量收敛曲线会发生振荡，此时由于迭代中分离涡的位置和强度都会发生不同程度的变化，呈现非定常特性，因此流量也会随之发生变化(但这种变化近似为周期性)。在这种情况下，也可认为计算收敛。 流场当地值：计算迭代收敛时，流场每一点处的参数值不应当再发生变化，或者对于有分离涡的情况，涡内某一点的参数应当为周期性变化。用户可以在FINE介面中跟踪某一特性点的参数，并观察其变化参数。 总体参数：对于定常计算,所有的总体性能(效率、扭矩、推力等）都应当变为恒定值，不再随迭代步数而发生变化。对于有大分离的情况，这些参数则会呈现周期性变化。这两种情况下都可认为计算收敛。对于非定常计算，所有的参

58、数都应当呈现近似周期性变化。,附：收敛标准,13.1 计算挂起：Solver/Suspend，求解器会在当前计算步完成后，存储结果并停止。,Step13：挂起计算,14.1 新建若干个计算，依次修改其边界条件及初场，以计算工况线。对于压气机计算，可通过改变出口背压来计算工况线。 14.2 保存.run文件,Step14：多工况点计算,将已经收敛的结果做为初场附给新的计算，可大大缩短收敛时间。,14.3 打开Task Manager 14.4 添加子任务，依次指定相 应的.run文件,Step14：多工况点计算,14.5 选择Run After Previous Subtask。 14.6 点击

59、Start开启计算。Fine会在第一个subtask完成后自动开启第二个。,Step14：多工况点计算,该功能在进行工况线计算时非常方便，用户只需提前设定好边界条件，在任务栏中添加任务，Fine会自动按顺序依次执行，而不需要在一个计算完成后人为开启下一个。,附：主要计算结果文件,.cgns 求解结果，存储流场内所有离散点的参数值，用于续算或者后处理（二进制存储格式） .mf 存储进出口平均值及总体性能（文本格式） .res 存储计算过程中的迭代残差（文本格式） .log 记录当前计算的所有信息，包括错误提示、警告信息等（文本格式） .steering 记录任务管理器中显示的跟踪参数信息（文本格式） .std 记录当前计算迭代过程的所有信息，包括挂起、删除、迭