数值模拟计算的整个过程Word格式文档下载.docx
《数值模拟计算的整个过程Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《数值模拟计算的整个过程Word格式文档下载.docx(20页珍藏版)》请在冰豆网上搜索。
1.2几何生成的方法
1.2.1泵的叶轮和导叶部分可以根据各自的木模图使用BLADEGEN较为方便的生成
1.2.2而其他部分则可以通过ProE等三维CAD工具生成,其中离心泵窝壳由窝壳木模图先将各断面绘制成型,再利用扫掠的方法成型。
1.3.几何输出
1.3.1从PRO/E中导出文件时可以选择保存成igs格式,也可以保存成stp格式,在导出时按其默认格式保存,即igs格式的保存成面的形式,stp格式的保存成体和壳的形式。
1.3.2.进出水流道部分(轴流泵),进口端(离心泵)要做适当的延伸。
1.3.3从PRO/E中导出之前可以可以改单位,或者明确几何生成时所用单位,以便导入。
1.3.4各部分的特征位置的坐标要明确,如几何中心,原点,以便各部分导入后的合并。
二.网格划分(软件:
ANSYSICEM)
网格划分主要有以下几部分:
2.1.几何检查及修复通过检查几何命令检查几何并将错误的部分根据实际情况修复(以轴流泵出水流道为例,见图2.1)
图2.1(a)轴流泵出水流道几何检查
图2.1(b)修复后的轴流泵出水流道几何
2.2设置part
图2.2设置part
2.3.建立物质点(生成四面体网格时必须)
图2.3建立物质点
2.4.设置网格大小,生成网格
2.4.1六面体网格的生成。
分块六面体网格生成主要是分块的思想,一般不外乎O型,C型,H型,J型。
就叶轮以及导叶,对包角较大的叶片采用J型,包角小的则采用H型网格,往往为提高质量可以在叶片表面附着一层O型网格;
对于离心泵的窝壳采用C型网格,其割舌部分用C型网格;
对于轴流泵进出口流道的拓朴结构可以根据具体几何形状划分块,可参照附图(图2.4)。
图2.4网格拓扑和生成网格
2.4.2近壁网格处理
近壁面的网格需要加密,以满足y+的要求。
对于
模型,y+值一般要求在100这个量级,Y+的具体要求如下:
模型,y+值一般要求在1这个量级,一般为4~5。
因此需要控制第一层网格的大小及近壁网格的数量。
第一层网格大小的估算公式如下
此外须在近壁区
厚度的范围内添加足够的网格数,要求如下:
其
的估算公式为
设置方法以下图为参考,首先确定通过上面的估算公式计算
及
,进行对块上网格的设置,并生成非结构网格,此时第一层网格的大小为
,然后根据公式
,计算网格增长率
,根据
将此时的第一层网格切成
份即可。
图2.4.1近壁面网格切分
2.5网格检查及修复
对于六面体网格,在未生成非结构化网格之前在blocking选项卡中进行网格质量的检查及光顺,对于生成非结构化网格之后的网格及四面体网格在editmesh选项卡中进行网格质量检查。
四面体网格可以通过网格光顺工具来提高网格质量,但不易调整。
六面体网格需要通过调整块节点的位置来提高网格质量;
块节点的调整需根据几何形状具体调整。
也可以通过editmesh中的选项卡中的工具进行网格修复。
一般而言按照quality方式进行网格质量检查时推荐quality大于0.3。
存在的问题:
1)综合网格质量quality是如何定义的,其值大于0.3的依据尚不明确
2)网格质量判定除了quality之外还有哪些重要的影响因素
图2.5网格质量检查及光顺
2.6.导出网格
首先要将网格转化成非结构网格(仅指六面体),再在output选项卡中选择求解器,再输出网格。
图2.6网格输出
三.CFX的计算
3.1建立一个simulation
打开CFX进入CFX-pre,在file菜单中建立一个newsimulation,选择generalmodel进入CFX-pre的界面。
3.1.1导入网格
进入mesh选项卡,选择importmesh或从file菜单中选择importmesh,进入网格导入界面。
在网格导入界面中,选择definition选项卡,在meshformat中的下拉菜单中选择ICEMCFD,在file中,选取所要引入的网格,并在meshunits中正确选择单位(见图3.1.1)。
一般轴流泵分进水流道,叶轮,导叶以及出水流道四个部分的网格;
离心泵分为叶轮,窝壳以及入口扩充断的三个部分的网格。
图3.1.1网格导入
由于各部分网格在几何造型时,位置安排不一定统一,以及导叶或者叶轮部分网格可通过Turbogrid或者ICEMCFD中生成一个流道的网格,所以网格在导入后需要调整位置,复制叶轮网格等操作。
对需要调整的部分在mesh选项卡中用鼠标右键单击该部分网格选择Transform,此时需要调整的网格部分变为绿色,并在窗口左下方definition选项卡中选择Transformation的类型,对于需要移动的部分选择translation选项,可以通过调整移动坐标来移动该部分;
对于需要旋转的部分选择rotation选项,通过对旋转轴及旋转角度设定可以旋转网格。
对于单一流道的叶片网格在旋转的同时需要重复拷贝该部分网格,需选择Multiplecopy选项,拷贝n-1个,n为叶片或者导叶的数量,并选择GlueMatchingAssumbles,这样可以将拷贝后的各部分网格合并为一个整体,这个必须保证相应边上网格完全一致,否则会产生多余的面。
(见图3.1.2)
图3.1.2网格调整
3.1.2计算定义
在tool菜单中选择Turbomode进入叶轮机计算模式,在basic设置中选择旋转轴及坐标系。
进入componentdefinition,建立一个新的component,在弹出的newcomponent菜单中,根据实际情况选择是静止或旋转的类型,其中叶轮为旋转部分,需定义转速,转速方向定义根据右手法则,其余部分为静止部分。
在meshvolume框体中的下拉菜单里选择相应的网格部分。
如果在网格面定义时命名规范则在turbomodel中regioninformation会自动适配相应的区域,并对该区域的边界条件进行定义。
进入physicaldefinition框体,选择计算类型(定常、非定常,非定常需定义参数)、湍流模式及进出口边界条件,对于非定场计算需设置timestep及totaltime,timestep根据CFL数预先估计。
一般选择massflowinletp-staticoutlet给定入口流量及出口静压值。
图3.1.3计算定义
进入interfacedefinition定义interface,其中转动部分和静止部分之间的类型是frozenrotor,静止部件之间是none。
进入boundarydefinition,定义边界条件,一般系统会自动的根据前面的设置自动定义边界条件,也可根据需要自行定义。
以上部分参见图3.1.3。
进入generalmode,根据需要设置迭代次数及残差。
残差推荐是1e-4。
可以在Solvecontrol中可以选择时间步(timescalcontrol和physicaltimescal,其中physicaltimescal为转速分之一)。
见图3.1.4。
保存设置。
图3.1.4solvercontrol定义
3.2进入计算
在菜单栏中选择writesolverfile按钮,弹出writesolverfile对话框,见图3.1.5,选择ok,进入CFX-solver界面(见图3.1.6),选择startrun进入计算。
图3.1.5writesolverfile对话框
图3.1.6CFX-solver界面
CFX设置中存在问题的部分:
1)在定义interface时,对于动静之间的interface其transformationtype有时需设置为none才可以计算,否则会提示出错。
2)入口湍流强度该如何定义。
3)出口边界为opening和outlet对计算结果的影响。
附件一:
CFX中设置的具体操作如下:
a.创建区域
b.设置边界条件
C.建立交界面
d.设置时间步、收敛精度及迭代次数