数值模拟计算整个过程
- 格式:doc
- 大小:1.96 MB
- 文档页数:17
SPCC板材冲压成形过程的数值模拟SPCC板材是一种优质的低碳钢板材,常用于制作汽车钣金、电器外壳等产品。
其具有良好的加工性能,但冲压成形过程中容易出现变形和裂纹等问题。
为了提高产品质量和生产效率,数值模拟成为一种重要的工具,可以帮助设计者优化冲压工艺和工装结构。
数值模拟的基本原理是根据材料力学特性和变形规律,通过计算机模拟来预测材料在冲压成形过程中的变形、应力和应变等参数,从而评估产品的质量和可行性。
下面我们以SPCC板材的冲压成形为例,介绍其数值模拟的过程和方法。
一、建立模型数值模拟的第一步是建立三维几何模型,包括产品形状、工艺参数、材料性质等要素。
在SPCC板材的冲压成形中,常见的产品包括箱体、盖板、底板等,需要根据不同产品的尺寸、形状和工艺要求来建立几何模型。
同时,需要设置材料参数,如杨氏模量、泊松比、屈服强度等。
二、网格划分建立几何模型后,需要将其转化为数值网格模型,也称为离散化。
这一过程主要是将复杂的几何体划分成数个小体积或小面积,以便于数值计算和求解。
网格划分的密度和精度决定了数值模拟的精度和速度。
通常,高精度的模拟需要更密的网格划分,但也会增加模拟时间和计算量。
三、建立物理模型在网格划分后,需要建立材料的物理模型,包括材料的本构模型和材料的本构参数。
本构模型是描述材料变形和应力关系的数学模型,常见的本构模型包括各向同性模型、非各向同性模型和弹塑性模型等。
用来描述SPCC板材的本构模型可以使用各向同性的屈服准则或者分段屈服的本构模型。
四、施加载荷和边界条件载荷是指在冲压成形过程中施加在材料表面的力和压力,这些载荷包括空气压力、上下模压力和侧向引伸力等。
边界条件是在数值模拟中必须考虑的约束条件,它们可以限制材料在变形过程中的位移、速度和加速度等。
通常,边界条件包括定位销、卡位、滚动辅助和止口等。
五、数值计算和分析完成模型建立和载荷设置后,就可以进行数值计算和分析。
该过程一般采用有限元分析(FEA)软件来完成,例如ABAQUS、ANSYS等。
基于CFX的离心泵内部流场数值模拟基于CFX的离心泵内部流场数值模拟随着计算流体力学和计算机技术的快速发展,泵内部的流动特征成为热点研究方向,目前应用CFX 软件的科研人员还较少,所以将CFX使用的基本过程加以整理供初学者参考。
如有不对之处敬请指教。
一、CFX数值计算的完整流程二、基于ICEM CFD的离心泵网格划分导入几何模型修整模型创建实体创建PRAT设置全局参数划分网格检查网格质量并光顺网格导出网格-选择求解器导出网格三、CFX-Pre 设置过程基本步骤新建文件导入网格定义模拟类型创建计算域指定边界条件建立交界面定义求解控制定义输出控制写求解器输入文件定义运行计算过程四、CFX-Post后处理计算泵的扬程和效率云图矢量图流线图导入几何模型在ICEM CFD软件界面内,单击File→Imort Geometry→STEP/IGES(一般将离心泵装配文件保存成STEP格式),将离心泵造型导入ICEM,如图3所示。
图3 导入几何模型界面修整模型单击Geometry→Repair Geometry→Build Topology,设置Tolerence,然后单击Apply,如图4所示。
拓扑分析后生成的曲线颜色指示邻近表面的关系:green = 自由边,yellow = 单边,red = 双边,blue =多边,线条颜色显示的开/关Model tree →Geometry → Curves → Color by count,Red curves 表示面之间的间隙在容差之内, 这是需要的物理模型,Yellow edges 通常是一些需要修补的几何。
图4 修整模型界面2-3 创建实体单击Geometry→Creade Body,详细过程如图5所示。
图5 创建实体界面创建PRAT创建PART,是为了设置边界时使用,在模型树中,右键点击Part,在出现菜单中选择Create Part。
以此创建各个部件的part,如图6所示。
数值模拟过程及历史拟合方法数值模拟过程及历史拟合方法是科学研究中常用的方法之一,它可以通过计算机模拟来探索不同系统的行为和变化规律。
在很多领域,数值模拟已经成为理论和实验研究的重要补充,它可以模拟各种物理场景、复杂的自然现象以及社会经济系统等。
数值模拟的过程一般包括以下几个步骤:1.定义问题和建立模型:首先需要明确研究问题的具体内容和边界条件,然后建立数学模型来描述问题。
模型的建立通常基于已知的理论和现象,可以是常微分方程、偏微分方程、代数方程等形式。
2.离散化:将连续的物理空间或时间离散化为有限的网格或时间步长。
这个过程通常需要将物理量转化为离散的数值,可以使用有限差分法、有限元法、谱方法等。
3.运用数值方法:利用数值方法求解离散后的问题。
常用的数值方法包括常微分方程数值解法、偏微分方程数值解法、随机数生成方法等。
4.模拟过程:根据所建立的数值模型和数值方法,通过计算机进行模拟运算。
在计算过程中,可以进行参数敏感性分析、收敛性分析等来确保结果的准确性和可靠性。
5.分析和解释:根据模拟结果进行分析和解释,得出结论并与实际情况进行对比。
通过与实验数据、观测数据等进行比较,可以验证模拟结果的合理性。
数值模拟的历史拟合方法是指利用已知的历史数据来拟合数学模型中的参数,以使模拟结果与观测结果尽可能吻合。
其中一个常用的历史拟合方法是最小二乘法。
最小二乘法是通过最小化实际观测值与模型预测值之间的残差平方和来确定模型的参数。
通过求解最小二乘问题的正规方程或使用迭代求解方法,可以得到最优的参数估计值。
另外一个常用的历史拟合方法是最大似然估计。
最大似然估计是假设观测数据来自于一些概率分布,在给定观测数据的条件下,寻找使得观测数据的概率最大化的参数估计值。
通过最大化似然函数或对数似然函数,可以得到最优的参数估计值。
历史拟合方法还包括遗传算法、粒子群优化算法等启发式算法。
这些方法通过模拟生物进化和群体行为的过程,来最优的参数组合。
数值模拟摘要:数值模拟是一种通过计算机模拟方法来研究和分析现实世界中的物理现象、工程问题和自然现象的方法。
本文将探讨数值模拟的原理、步骤和应用场景,并讨论其优点和限制。
1. 引言数值模拟是一种基于计算机技术的仿真方法,可用于模拟和研究各种自然和工程现象。
它通过利用数值计算方法解决传统试验无法解决或者很难解决的问题。
2. 数值模拟的原理和步骤数值模拟的基本原理是将问题转化为数学模型,并通过计算方法求解该模型。
它通常包括以下步骤:2.1 问题建模在数值模拟中,首先需要对待解问题进行建模。
建模的目的是将实际问题转化为数学模型,包括确定问题的边界条件、初值条件和物理方程等。
2.2 离散化离散化是将连续的问题转化为离散的数值问题。
例如,在求解连续介质力学问题时,可以通过将物理空间离散为网格点,并对网格点上的物理量进行离散化处理。
2.3 数值求解数值求解是数值模拟的核心步骤,涉及到使用数值方法和算法对离散化后的问题进行求解。
常用的数值方法包括有限差分法、有限元法、边界元法等。
2.4 结果分析数值模拟的最终结果需要进行分析和验证。
分析结果可以通过与理论分析、实验结果或其他已有数据进行比对来验证其准确性和可靠性。
3. 数值模拟的应用场景数值模拟广泛应用于各个领域,包括物理学、化学、生物学、工程学和计算机科学等。
3.1 天气预报数值模拟在天气预报中有着重要的应用。
通过对大气物理方程进行离散化和数值求解,可以对天气系统进行模拟预测,并提供准确的天气预报。
3.2 污染扩散模拟污染扩散模拟是评估污染物排放对环境影响的重要手段。
通过模拟和计算污染物在大气、水体或土壤中的传输和扩散过程,可以评估污染物的浓度分布和危害程度。
3.3 车辆碰撞模拟车辆碰撞模拟可以通过数值模拟来研究交通事故的发生机理和影响因素。
通过建立车辆和人体的力学模型,并对碰撞过程进行数值求解,可以评估碰撞对车辆和人体的影响。
4. 数值模拟的优点和限制数值模拟作为一种研究方法具有以下优点:4.1 成本低廉相对于传统试验方法,数值模拟不需要大量的实验设备和人力资源,能够在计算机上进行模拟和求解,降低了成本。
一、名词解释差分近似的收敛性:在模式的计算域内,在固定时间n∆t 后,当∆t,∆x→0时,在整个计算域内max{|u i n−u(i∆x,n∆t)|}→0。
混淆误差:指由于有限网格不能正确地分解短波而造成的误差。
计算的稳定性:给定初始值,计算的解随着时间步长(趋于∞时),是有界的。
差分近似的一致性(相容性):当∆t,∆x→0时,差分系统应等同于微分方程。
模式中的动力框架:包括建立相适应的大气动力方程组、确定垂直坐标系和选择合理的离散化的方法(差分还是谱分析)。
物理过程:主要包括3个,即云物理过程参数化、边界层物理过程参数化和太阳辐射过程参数化。
资料同化:即把不同时刻、不同地区、不同类型的观测资料与常规资料融合,通过一定的预报模式,使之在动力和热力上协调起来,求得质量场和流场基本平衡的理想的初始场。
客观分析:利用模式构造一个背景场,再把测站资料插值到该背景场的格点上,然后做质量检验并对其进行初始化。
是整合检错,质量控制,合理插值的综合过程。
Lax等价性原理:对于一个适定的初值问题和它的一个具有相容性的差分格式,其计算稳定性是收敛性的充分必要条件。
CFL条件:显式格式稳定性的必要要件,c∆t/∆x≤1,其中C为所考虑的最快波的波速。
不同离散方法中模式的精度:对于差分方法,格距越小,分辨率越高;对于谱方法,截断系数越大,分辨率越高。
二、数值模式由哪几个部分组成?1、前处理:包括客观分析等,给主模式一个合理的适应模式的初始条件。
2、动力框架、物理过程3、后处理:模式输出的资料的处理,以供应用。
MM5模式有哪些优点?(1)多重嵌套和移动嵌套功能;(2) 非静力动力模式(3) 四维同化能力;(4) 齐全的物理选项;(5)可以多平台、多系统运行(6) 易操作性。
MM5模式有几个模块组成,各自的主要功能是什么?MM5共有8个模块。
(1)TERRAIN模块:对经纬度网格的地形资料和路面特征参数资料进行插值。
(2)REGRID模块:对分析场资料进行插值,形成标准等压面上的初值场。
四、非线笥问题迭代式解法的收敛性每一层次上满足迭代法求解的收敛条件+相邻次间代数方程的系数变化不太大(亦即未知量的变化不太大J多数情形下非线性问题迭代式解法是可以收敛的)。
使相邻两层次间未知量变化不太大的措施:1欠松弛迭代常用逐次欠弛线迭法(SLUR):一组临时系数下逐线迭代求解+对所得的解施以欠松弛,再用欠松弛后的解去计算新的系数,常数,以进入下一层次的迭代。
实施:常把欠松弛处理纳入迭代过程,而不是在一个层次迭代完成后再行欠松弛。
.(n 1)川).'a n bt n bt p =t p (t p )ap(先)t p n1) = 7an b t n b b (1一•)屯t p n)co oa'p t p n 9 、a n bt n b b'a'p -a^ ■, b' = b (^ )(a p )t p n),用交替方向线迭代法求解这一方程,就实现了SLUR的迭代求解。
为一般化起见,上式中t n b上没有标以迭代层次的符号(J, GS时不相同)。
2、采用拟非稳态法前面已指出,稳态问题的迭代解法与非稳态问题的步进法十分相似。
对于非线性稳态问题,从代数方程的一组临时系数进入到另一组临时系数亦好象非稳态问题前进了一个时间层,非稳态问题的物理特性:系数热惯性越大(a; = PM v/也I ),温度变化越慢,仿此,对稳态非线性问题,可在离散方程中加入拟非稳态项,以减小未知量托两个层次间的变化,即由(=a n b -S p:V)t p n。
= Ua n bt n b b=(3a n b - S p:V a;)t p n。
=二a n bt n b b a p tfZa n bt n b - b - a;t p n)(n 1)t p oEa n b -S p心V +a p一直进行到t p,t n b收敛,虚拟时间步的大小通过计算实践确定。
3、采用Jacobi点迭代法中止迭代的判据(该层次迭代)除前述变化率判据外,还可以规定迭代的轮数,例如规定进行4-6次ADI线迭代就结束该层次上的计算。
数值模拟计算的整个过程数值模拟计算的整个过程主要包括一下几个过程:一.建立模型(应用软件: CAD 工具如 PRO/E Bladegen 等)几何生成时应注意的问题主要有以下几个部分:1.几何生成1.1几何区域的规划几何的生成可以是一个整体部分,但是有时为了网格划分时的方 便可以把几个分成几个部分生成,例如轴流泵几何的生成可以分为四个部分:进水流道、叶轮、导叶和出水流道(图 1.2),离心泵几何分为三个部分:进口端,叶轮,窝壳(图 1.2 )。
图1.1轴流泵几何图1.2离心泵几何fl 9U 1,痂(m)--------------- ------------------------------- 1 1 ________________ 1 1出水流道叶轮1.2几何生成的方法1.2.1泵的叶轮和导叶部分可以根据各自的木模图使用BLADEGE较为方便的生成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各部分的特征位置的坐标要明确,如几何中心,原点,以便各部分导入后的合并。
二.网格划分(软件:ANSYS ICEM)网格划分主要有以下几部分:2.1.几何检查及修复通过检查几何命令检查几何并将错误的部分根据实际情况修复(以轴流泵出水流道为例,见图 2.1)图2.1 (a)轴流泵出水流道几何检查图2.1 (b )修复后的轴流泵出水流道几何2.2 设置 part2.4.设置网格大小,生成网格 2.4.1六面体网格的生成。
Lagrange、Euler、ALE三种方法的简单介绍ALE、Lagrange、Euler是数值模拟中处理连续体的广泛应用的三种方法。
Lagrange方法多用于固体结构的应力应变分析,这种方法以物质坐标为基础,其所描述的网格单元将以类似“雕刻”的方式划分在用于分析的结构上,即是说采用Lagrange方法描述的网格和分析的结构是一体的,有限元节点即为物质点。
采用这种方法时,分析结构的形状的变化和有限单元网格的变化完全是一致的(因为有限元节点就为物质点),物质不会在单元与单元之间发生流动。
这种方法主要的优点是能够非常精确的描述结构边界的运动,但当处理大变形问题时,由于算法本身特点的限制,将会出现严重的网格畸变现象,因此不利于计算的进行。
Euler方法以空间坐标为基础,使用这种方法划分的网格和所分析的物质结构是相互独立的,网格在整个分析过程中始终保持最初的空间位置不动,有限元节点即为空间点,其所在空间的位置在整个分析过程始终是不变的。
很显然由于算法自身的特点,网格的大小形状和空间位置不变,因此在整个数值模拟过程中,各个迭代过程中计算数值的精度是不变的。
但这种方法在物质边界的捕捉上是困难的。
多用于流体的分析中。
使用这种方法时网格与网格之间物质是可以流动的。
ALE方法最初出现于数值模拟流体动力学问题的有限差分方法中。
这种方法兼具Lagrange方法和Euler方法二者的特长,即首先在结构边界运动的处理上它引进了Larange方法的特点,因此能够有效的跟踪物质结构边界的运动;其次在内部网格的划分上,它吸收了Euler 的长处,即是使内部网格单元独立于物质实体而存在,但它又不完全和Euler网格相同,网格可以根据定义的参数在求解过程中适当调整位置,使得网格不致出现严重的畸变。
这种方法在分析大变形问题时是非常有利的。
使用这种方法时网格与网格之间物质也是可以流动的。
固体结构分析中一般都选用lagrange坐标,实际上lagrange euler法在有限元中体现的节点意义正如楼主所述,但是本质牵扯的是参考什么样的坐标来描述应力应变关系。
数值模拟计算的整个过程数值模拟计算的整个过程主要包括一下几个过程:一.建立模型(应用软件:CAD工具如PRO/E,Bladegen等)几何生成时应注意的问题主要有以下几个部分:1. 几何生成1.1 几何区域的规划几何的生成可以是一个整体部分,但是有时为了网格划分时的方便可以把几个分成几个部分生成,例如轴流泵几何的生成可以分为四个部分:进水流道、叶轮、导叶和出水流道(图1.2),离心泵几何分为三个部分:进口端,叶轮,窝壳(图1.2)。
图1.1 轴流泵几何图1.2 离心泵几何1.2几何生成的方法1.2.1泵的叶轮和导叶部分可以根据各自的木模图使用BLADEGEN较为方便的生成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各部分的特征位置的坐标要明确,如几何中心,原点,以便各部分导入后的合并。
二.网格划分(软件: 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型,包角小的则采用H型网格,往往为提高质量可以在叶片表面附着一层O型网格;对于离心泵的窝壳采用C型网格,其割舌部分用C型网格;对于轴流泵进出口流道的拓朴结构可以根据具体几何形状划分块,可参照附图(图2.4)。
出水流道网格块结构和生成网格叶片J 网格拓扑和生成网格导叶H 型网格拓扑和生成网格图2.4 网格拓扑和生成网格2.4.2 近壁网格处理近壁面的网格需要加密,以满足y+的要求。
对于k ε-模型,y+值一般要求在100这个量级,Y+的具体要求如下:对于k ω-模型,y+值一般要求在1这个量级,一般为4~5。
因此需要控制第一层网格的大小及近壁网格的数量。
第一层网格大小的估算公式如下此外须在近壁区δ厚度的范围内添加足够的网格数,要求如下:其δ的估算公式为设置方法以下图为参考,首先确定通过上面的估算公式计算δ及wall y ,进行对块上网格的设置,并生成非结构网格,此时第一层网格的大小为δ,然后根据公式normall N wall y δα=,计算网格增长率α,根据α将此时的第一层网格切成normall N 份即可。
图2.4.1近壁面网格切分2.5 网格检查及修复对于六面体网格,在未生成非结构化网格之前在blocking选项卡中进行网格质量的检查及光顺,对于生成非结构化网格之后的网格及四面体网格在edit mesh选项卡中进行网格质量检查。
四面体网格可以通过网格光顺工具来提高网格质量,但不易调整。
六面体网格需要通过调整块节点的位置来提高网格质量;块节点的调整需根据几何形状具体调整。
也可以通过edit mesh中的选项卡中的工具进行网格修复。
一般而言按照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菜单中建立一个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)。
一般轴流泵分进水流道,叶轮,导叶以及出水流道四个部分的网格;离心泵分为叶轮,窝壳以及入口扩充断的三个部分的网格。
图3.1.1网格导入由于各部分网格在几何造型时,位置安排不一定统一,以及导叶或者叶轮部分网格可通过Turbogrid或者ICEMCFD中生成一个流道的网格,所以网格在导入后需要调整位置,复制叶轮网格等操作。
对需要调整的部分在mesh选项卡中用鼠标右键单击该部分网格选择Transform,此时需要调整的网格部分变为绿色,并在窗口左下方definition选项卡中选择Transformation的类型,对于需要移动的部分选择translation选项,可以通过调整移动坐标来移动该部分;对于需要旋转的部分选择rotation选项,通过对旋转轴及旋转角度设定可以旋转网格。
对于单一流道的叶片网格在旋转的同时需要重复拷贝该部分网格,需选择Multiple copy选项,拷贝n-1个,n为叶片或者导叶的数量,并选择Glue Matching Assumbles,这样可以将拷贝后的各部分网格合并为一个整体,这个必须保证相应边上网格完全一致,否则会产生多余的面。
(见图3.1.2)图3.1.2网格调整3.1.2计算定义在tool菜单中选择Turbo mode进入叶轮机计算模式,在basic设置中选择旋转轴及坐标系。
进入component definition,建立一个新的component,在弹出的new component菜单中,根据实际情况选择是静止或旋转的类型,其中叶轮为旋转部分,需定义转速,转速方向定义根据右手法则,其余部分为静止部分。
在mesh volume框体中的下拉菜单里选择相应的网格部分。
如果在网格面定义时命名规范则在turbo model中region information 会自动适配相应的区域,并对该区域的边界条件进行定义。
进入physical definition框体,选择计算类型(定常、非定常,非定常需定义参数)、湍流模式及进出口边界条件,对于非定场计算需设置time step及total time,time step根据CFL 数预先估计。
一般选择mass flow inlet p-static outlet给定入口流量及出口静压值。
图3.1.3计算定义进入interface definition 定义interface,其中转动部分和静止部分之间的类型是frozen rotor,静止部件之间是none。
进入boundary definition,定义边界条件,一般系统会自动的根据前面的设置自动定义边界条件,也可根据需要自行定义。
以上部分参见图3.1.3。
进入general mode,根据需要设置迭代次数及残差。
残差推荐是1e-4。
可以在Solve control 中可以选择时间步(timescal control和physical timescal,其中physical timescal为转速分之一)。
见图3.1.4。
保存设置。
图3.1.4 solver control 定义3.2进入计算在菜单栏中选择write solver file 按钮,弹出write solver file对话框,见图3.1.5,选择ok,进入CFX-solver界面(见图3.1.6),选择start run 进入计算。
图3.1.5 write solver file对话框图3.1.6 CFX-solver界面CFX设置中存在问题的部分:1)在定义interface时,对于动静之间的interface其transformation type有时需设置为none 才可以计算,否则会提示出错。
2)入口湍流强度该如何定义。
3)出口边界为opening和outlet对计算结果的影响。
附件一:CFX中设置的具体操作如下:a.创建区域1modmodmod...location jsdomaintype fluid domaingeneral optionsfluid list waterdomaindomainmotionoption stationaryheattrancefer el option nonefluid elturbulence el option kdcreat adomainε⎧→⎧⎪⎪→⎪⎪⎨⎪→⎪⎪→⎨⎪→⎩⎪⎪→⎧⎪⎨→-⎪⎩⎩2()modmodmod...location impellerdomaintype fluid domaingeneral optionsfluid list wateromaindomainmotionoption rotationaryheattrancefer el option nonefluid elturbulence el option kε⎧→⎧⎪⎪→⎪⎪⎨⎪→⎪⎪→⎨⎪→⎩⎪⎪→⎧⎪⎨→-⎩⎩注意转速方向3modmodmod...4location guidedomaintype fluid domaingeneral optionsfluid list waterdomaindomainmotionoption stationaryheattrancefer el option nonefluid elturbulence el option kgdomainε⎪⎧→⎧⎪⎪→⎪⎪⎨⎪→⎪⎪→⎨⎪→⎩⎪⎪→⎧⎪⎨→-⎪⎩⎩→modmodmod...location csdomaintype fluid domaineneral optionsfluid list waterdomainmotionoption stationaryheattrancefer el option nonefluid elturbulence el option kε⎧⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧→⎧⎪⎪→⎪⎪⎨⎪→⎪⎪⎨⎪→⎩⎪⎪→⎧⎪⎨→-⎪⎩⎩⎩⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪b.设置边界条件1det boundarytype inlet basic settings location name inlet creat boundary domain domain mass and momentum boundary ail ⎧⎧→⎧⎪⎨⎪⎩⎪⎪→⎫⎪⎪→⎧⎧⎨⎬⎨→⎭⎪⎨⎪⎪⎨⎩⎪⎪⎪⎪⎪⎩⎩⎩质量流量速度入口wall 1noslip det boundarytype basic settings location name wall creatboundary domain domain boundary ail smoothwall ⎧⎧⎧⎪⎪⎪⎨⎪⎫⎪⎪⎪⎩⎨⎬⎨⎧⎭⎪⎪⎪⎨⎪⎪⎪⎩⎩⎪⎩→→→→opening out 1det boundarytype outlet basicsettings location name let creatboundary domain domain massand momentum boundary ail ⎧⎧⎧⎪⎪⎪⎨⎪⎪⎪⎩⎪⎪⎫⎪⎪⎧⎧⎨⎬⎨⎪⎪⎭⎪⎪⎪⎨⎨⎪⎪⎪⎩⎪⎪⎪⎪⎩⎪⎪⎩⎩→→→→()质量流量相对压力C.建立交界面int erface domain domain1int erface side1region list interface11creat int erface name interface domain domain2int erface side2region list interface12int mod type fluid fluid erface els option general connecti →→⎧⎨→⎩→→⎧⎨→⎩→(1)/on framechange mixing option frozenrotor ⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪→⎩int erface domain domain2int er faceside1region list interface21creat int erface name interface domain domain3int er faceside2region list interface22int mod type fluid fluid erface els option general connecti →→⎧⎨→⎩→→⎧⎨→⎩→(2)/on framechange mixing option frozenrotor⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪→⎩ int erface domain domain3int erface side1region list interface31creat int erface name interface domain domain4int erface side2region list interface32int mod type fluid fluid erface els option general connecti →→⎧⎨→⎩→→⎧⎨→⎩→(3)/none on framechange mixing option ⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪→⎩d.设置时间步、收敛精度及迭代次数(,..)max int solve ax convegencecriteria arg ()advection schemeoption high resolution up wind erations timescalecontrol fluid timescalecontrol control physical timescale residual type residual t et ⎧→⎨→⎩⎧⎨⎩()转速分之一(RMS,M )收敛精度⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩。