1、数值模拟计算的整个过程数值模拟计算的整个过程主要包括一下几个过程:一. 建立模型(应用软件: CAD 工具如 PRO/E,Bladegen 等)几何生成时应注意的问题主要有以下几个部分:1. 几何生成1.1 几何区域的规划 几何的生成可以是一个整体部分,但是有时为了网格划分时的方便可以把几个分成几个部分生成,例如轴流泵几何的生成可以分为四个部分:进水流道、叶轮、导叶和出水流道(图 1.2) ,离心泵几何分为三个部分:进口端,叶轮,窝壳(图1.2) 。图 1.1 轴流泵几何图 1.2 离心泵几何1.2 几何生成的方法 1.2.1 泵的叶轮和导叶部分可以根据各自的木模图使用 BLADEGEN 较为
2、方便的生成1.2.2 而其他部分则可以通过 Pro E 等三维 CAD 工具生成,其中离心泵窝壳由窝壳木模图先将各断面绘制成型,再利用扫掠的方法成型。1.3.几何输出1.3.1 从 PRO/E 中导出文件时可以选择保存成 igs 格式,也可以保存成 stp 格式,在导出时按其默认格式保存,即 igs 格式的保存成面的形式, stp 格式的保存成体和壳的形式。1.3.2. 进出水流道部分(轴流泵) ,进口端(离心泵)要做适当的延伸。1.3.3 从 PRO/E 中导出之前可以可以改单位,或者明确几何生成时所用单位,以便导入。1.3.4 各部分的特征位置的坐标要明确,如几何中心,原点,以便各部分导入
3、后的合并。二. 网格划分(软件: ANSYS ICEM )网格划分主要有以下几部分:2.1. 几何检查及修复 通过检查几何命令检查几何并将错误的部分根据实际情况修复(以轴流泵出水流道为例,见图 2.1)图 2.1(a ) 轴流泵出水流道几何检查图 2.1(b)修复后的轴流泵出水流道几何2.2 设置 part 图 2.2 设置 part2.3. 建立物质点 (生成四面体网格时必须)图 2.3 建立物质点2.4. 设置网格大小,生成网格2.4.1 六面体网格的生成。分块六面体网格生成主要是分块的思想,一般不外乎 O 型,C 型,H 型,J 型。就叶轮以及导叶,对包角较大的叶片采用 J 型,包角小的
4、则采用 H 型网格,往往为提高质量可以在叶片表面附着一层 O 型网格;对于离心泵的窝壳采用 C 型网格,其割舌部分用 C 型网格;对于轴流泵进出口流道的拓朴结构可以根据具体几何形状划分块,可参照附图(图2.4) 。叶片 J 网格拓扑和生成网格出水流道网格块结构和生成网格图 2.4 网格拓扑和生成网格2.4.2 近壁网格处理近壁面的网格需要加密,以满足 y+的要求。对于 模型,y+值一般要求在 100 这个量级,Y+的具体要求如下:k导叶 H 型网格拓扑和生成网格对于 模型,y+值一般要求在 1 这个量级,一般为 45。k因此需要控制第一层网格的大小及近壁网格的数量。第一层网格大小的估算公式如下
5、此外须在近壁区 厚度的范围内添加足够的网格数,要求如下:其 的估算公式为设置方法以下图为参考,首先确定通过上面的估算公式计算 及 ,进行对块上网格的waly设置,并生成非结构网格,此时第一层网格的大小为 ,然后根据公式 ,normalNwly计算网格增长率 ,根据 将此时的第一层网格切成 份即可。normalN图 2.4.1 近壁面网格切分2.5 网格检查及修复对于六面体网格,在未生成非结构化网格之前在 blocking 选项卡中进行网格质量的检查及光顺,对于生成非结构化网格之后的网格及四面体网格在 edit mesh 选项卡中进行网格质量检查。四面体网格可以通过网格光顺工具来提高网格质量,但
6、不易调整。六面体网格需要通过调整块节点的位置来提高网格质量;块节点的调整需根据几何形状具体调整。也可以通过 edit mesh 中的选项卡中的工具进行网格修复。一般而言按照 quality 方式进行网格质量检查时推荐 quality 大于 0.3。存在的问题:1)综合网格质量 quality 是如何定义的,其值大于 0.3 的依据尚不明确2)网格质量判定除了 quality 之外还有哪些重要的影响因素图 2.5 网格质量检查及光顺2.6. 导出网格首先要将网格转化成非结构网格(仅指六面体) ,再在 output 选项卡中选择求解器,再输出网格。图 2.6 网格输出三. CFX 的计算3.1 建
7、立一个 simulation打开 CFX 进入 CFX-pre,在 file 菜单中建立一个 new simulation,选择 general model进入 CFX-pre 的界面。3.1.1 导入网格 进入 mesh 选项卡,选择 import mesh 或从 file 菜单中选择 import mesh,进入网格导入界面。在网格导入界面中,选择 definition 选项卡,在 mesh format 中的下拉菜单中选择ICEM CFD,在 file 中,选取所要引入的网格,并在 mesh units 中正确选择单位(见图3.1.1) 。一般轴流泵分进水流道,叶轮,导叶以及出水流道四个
8、部分的网格;离心泵分为叶轮,窝壳以及入口扩充断的三个部分的网格。图 3.1.1 网格导入由于各部分网格在几何造型时,位置安排不一定统一,以及导叶或者叶轮部分网格可通过 Turbogrid 或者 ICEMCFD 中生成一个流道的网格,所以网格在导入后需要调整位置,复制叶轮网格等操作。对需要调整的部分在 mesh 选项卡中用鼠标右键单击该部分网格选择 Transform,此时需要调整的网格部分变为绿色,并在窗口左下方 definition 选项卡中选择 Transformation 的类型,对于需要移动的部分选择 translation 选项,可以通过调整移动坐标来移动该部分;对于需要旋转的部分选
9、择 rotation 选项,通过对旋转轴及旋转角度设定可以旋转网格。对于单一流道的叶片网格在旋转的同时需要重复拷贝该部分网格,需选择Multiple copy 选项,拷贝 n-1 个,n 为叶片或者导叶的数量,并选择 Glue Matching Assumbles,这样可以将拷贝后的各部分网格合并为一个整体,这个必须保证相应边上网格完全一致,否则会产生多余的面。 (见图 3.1.2)图 3.1.2 网格调整3.1.2 计算定义在 tool 菜单中选择 Turbo mode 进入叶轮机计算模式,在 basic 设置中选择旋转轴及坐标系。进入 component definition,建立一个新的
10、 component,在弹出的 new component 菜单中,根据实际情况选择是静止或旋转的类型,其中叶轮为旋转部分,需定义转速,转速方向定义根据右手法则,其余部分为静止部分。在 mesh volume 框体中的下拉菜单里选择相应的网格部分。如果在网格面定义时命名规范则在 turbo model 中 region information 会自动适配相应的区域,并对该区域的边界条件进行定义。进入 physical definition 框体,选择计算类型(定常、非定常,非定常需定义参数) 、湍流模式及进出口边界条件,对于非定场计算需设置 time step 及 total time,time step 根据CFL 数预先估计。一般选择 mass flow inlet p-static outlet 给定入口流量及出口静压值。