资源描述
2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-1,ANSYS/FLOTRAN流体动力学(CFD)分析,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-2,目录,第一章FLOTRAN计算流体动力学(CFD)分析概述一、FLOTRANCFD分析的概念二、FLOTRAN分析的种类1、层流分析2、紊流分析3、热分析4、可压缩流分析5、非牛顿流分析6、多组份传输分析,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-3,目录,第二章FLOTRAN分析基础一、FLOTRAN单元的特点1、FLUID141单元2、FLUID142单元3、FLUID单元的其他特征二、使用FLOTRAN单元的一些限制及注意事项三、FLOTRAN分析的主要步骤四、FLOTRAN分析中产生的一些文件五、提高收敛性和稳定性的常用的工具六、FLOTRAN分析过程中应处理的问题七、对一个FLOTRAN分析进行评价八、验证结果,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-4,目录,第三章FLOTRAN设置命令一、FLOTRAN求解控制命令二、FLOTRAN执行及输出控制命令稳态控制参数设置三、FLOTRAN执行及输出控制命令瞬态控制参数设置四、FLOTRAN输出及保存文件控制五、FLOTRAN输出收敛监测量的控制六、FLOTRAN流体类型及其特性的可变性控制七、FLOTRAN流体性质参数设定八、FLOTRAN分析的坐标系统的指定九、FLOTRAN分析参考条件的设置十、指定FLOTRAN分析的旋转坐标系统十一、指定FLOTRAN分析的重力加速度十二、设置有助于FLOTRAN求解稳定的参数,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-5,目录,第三章FLOTRAN设置命令(续)十三、设定FLOTRAN自由度松弛系数十四、设定FLOTRAN流体性质松弛因子十五、设置FLOTRAN分析的自由度限值十六、选择FLOTRAN各自由度相应的求解器十七、对FLOTRAN各求解器的控制十八、设置FLOTRAN湍流模型的一些常数十九、重新设定FLOTRAN各分析参数的值二十、控制FLOTRAN面积积分的阶次二十一、FLOTRAN多组份疏运分析的设置及控制二十二、定义FLOTRAN的重启动(续算)控制二十三、设置并执行一个零迭代FLOTRAN分析,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-6,目录,第四章FLOTRAN边界条件一、边界条件的类型二、进口/出口边界条件三、对称边界条件四、固定壁面五、移动壁面六、内流七、外流八、周期边界九、热边界条件十、边界条件总结十一、边界条件的菜单形式,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-7,目录,第五章FLOTRAN层流和湍流分析算例一、问题描述二、分析方法及假定三、几何尺寸及流体性质四、分析过程第1步:进入ANSYS第2步:设置分析选择第3步:定义单元类型第4步:生成分析区域的几何面第5步:定义单元形状第6步:划分有限元网格第7步:生成并应用新的工具栏按钮第8步:施加边界条件第9步:求解层流,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-8,目录,第五章FLOTRAN层流和湍流分析算例(续)四、分析过程第10步:观察层流分析的结果第11步:确定流体粘性如何影响流场特性第12步:进行湍流分析第13步:对新的出口区划分网格第14步:施加湍流分析的载荷第15步:改变FLOTRAN分析选项和流体性第16步:进行求解第17步:将流体速度结果以向量图和路径图的方式进行显示第18步:绘制压力等值线图,Guidelines,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-9,第一章FLOTRAN计算流体动力学(CFD)分析概述,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-10,在完成本章学习后,我们应该对流体动力学分析的基本概念有所了解,并知道它的基本分类。,ModuleObjective,第一讲、FLOTRANCFD分析的概念第二讲、FLOTRAN分析的种类第三讲、层流分析第四讲、紊流分析第五讲、热分析第六讲、可压缩流分析第七讲、非牛顿流分析第八讲、多组份传输分析,LessonObjectives,目标,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-11,一、FLOTRANCFD分析的概念,T-1.FLOTRANCFD分析的概念ANSYS程序中的FLOTRANCFD分析功能是一个用于分析二维及三维流体流动场的先进的工具,使用ANSYS中用于FLOTRANCFD分析的FLUID141和FLUID142单元,可解决如下问题:作用于气动翼(叶)型上的升力和阻力超音速喷管中的流场弯管中流体的复杂的三维流动同时,FLOTRAN还具有如下功能:计算发动机排气系统中气体的压力及温度分布研究管路系统中热的层化及分离使用混合流研究来估计热冲击的可能性用自然对流分析来估计电子封装芯片的热性能对含有多种流体的(由固体隔开)热交换器进行研究,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-12,二、FLOTRAN分析的种类,T-2.FLOTRAN分析的种类FLOTRAN可执行如下分析:层流或紊流传热或绝热可压缩或不可压缩牛顿流或非牛顿流多组份传输这些分析类型并不相互排斥,例如,一个层流分析可以是传热的或者是绝热的,一个紊流分析可以是可压缩的或者是不可压缩的。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-13,层流分析,T-2.FLOTRAN分析的种类层流中的速度场都是平滑而有序的,高粘性流体(如石油等)的低速流动就通常是层流。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-14,紊流分析,T-2.FLOTRAN分析的种类紊流分析用于处理那些由于流速足够高和粘性足够低从而引起紊流波动的流体流动情况,ANSYS中的二方程紊流模型可计及在平均流动下的紊流速度波动的影响。如果流体的密度在流动过程中保持不变或者当流体压缩时只消耗很少的能量,该流体就可认为是不可压缩的,不可压缩流的温度方程将忽略流体动能的变化和粘性耗散。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-15,热分析,T-2.FLOTRAN分析的种类流体分析中通常还会求解流场中的温度分布情况。如果流体性质不随温度而变,就可不解温度方程。在共轭传热问题中,要在同时包含流体区域和非流体区域(即固体区域)的整个区域上求解温度方程。在自然对流传热问题中,流体由于温度分布的不均匀性而导致流体密度分布的不均匀性,从而引起流体的流动,与强迫对流问题不同的是,自然对流通常都没有外部的流动源。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-16,可压缩流分析,T-2.FLOTRAN分析的种类对于高速气流,由很强的压力梯度引起的流体密度的变化将显著地影响流场的性质,ANSYS对于这种流动情况会使用不同的解算方法。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-17,非牛顿流分析,T-2.FLOTRAN分析的种类应力与应变率之间成线性关系的这种理论并不能足以解释很多流体的流动,对于这种非牛顿流体,ANSYS程序提供了三中粘性模式和一个用户自定义子程序。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-18,多组份传输分析,T-2.FLOTRAN分析的种类这种分析通常是用于研究有毒流体物质的稀释或大气中污染气体的传播情况,同时,它也可用于研究有多种流体同时存在(但被固体相互隔开)的热交换分析。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-19,第二章FLOTRAN分析基础,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-20,本章学习,我们简单介绍FLOTRAN分析的基础知识,以助于大家对以后的FLOTRAN分析有个大体的概念。,ModuleObjective,第一讲、FLOTRAN单元的特点第二讲、FLUID141单元第三讲、FLUID142单元第四讲、FLUID单元的其他特征第五讲、使用FLOTRAN单元的一些限制及注意事项第六讲、FLOTRAN单元使用中的一些限制第七讲、FLOTRAN分析的主要步骤第八讲、FLOTRAN分析中产生的一些文件第九讲、提高收敛性和稳定性的常用的工具第十讲、FLOTRAN分析过程中应处理的问题,LessonObjectives,目标,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-21,一、FLOTRAN单元的特点,T2-1.FLOTRAN单元的特点ANSYS中的FLOTRAN单元,即FLUID141和FLUID142,用于解算单相粘性流体的二维和三维流动、压力和温度分布。对于这些单元,ANSYS通过质量、动量和能量三个守恒性质来计算流体的速度分量、压力、以及温度。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-22,FLUID141单元,T2-1.FLOTRAN单元的特点FLUID141单元具有下列特征:维数:二维形状:四节点四边形或三节点三角形自由度:速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额,Objective,FLUID141单元,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-23,FLUID142单元,T2-1.FLOTRAN单元的特点FLUID142单元具有下列特征:维数:三维形状:四节点四面体或八节点六面体自由度:速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额,Objective,FLUID142单元,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-24,FLUID单元的其他特征,T2-1.FLOTRAN单元的特点FLUID单元的其他特征包括:用于模拟紊流的二方程紊流模式有很多推导结果,诸如:流场分析中的马赫数、压力系数、总压、剪应力、壁面处的y-plus、以及流线函数;热分析中的热流、热交换(膜)系数等。流体边界条件,包括:速度、压力、紊流动能以及紊流能量耗散率。用户无需提供流场进口处紊流项的边界条件,因FLOTRAN对此提供的缺省值适用于绝大多数分析。热边界条件,包括:温度、热流、体积热源、热交换(膜)系数。用户可使用的坐标系有:的卡尔坐标系、柱坐标系、极坐标系和轴对称坐标系。如果所计算的问题是轴对称的,激活旋转(swirl)选项即可算出垂直于对称平面的速度分量。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-25,二、使用FLOTRAN单元的一些限制及注意事项,T2-2.FLOTRAN单元的特点FLOTRAN单元的一些局限性:在同一次分析中不能改变求解的区域单元不支持自由流面边界条件ANSYS程序的某些特征不能同FLOTRAN单元一起使用使用FLOTRAN单元时不能使用某些命令或菜单当使用ANSYS的图形用户界面时,程序将只能显示那些在菜单和对话框中的FLOTRANSetUp部分要求了的特征和选项。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-26,FLOTRAN单元使用的一些限制,T2-2.FLOTRAN单元的特点当使用FLOTRAN单元时,要避免使用ANSYS的某些特征和命令,至少,要注意到在使用FLOTRAN单元时与别的分析稍微有些不同,当使用了无效的命令时,程序会给出相应的警告或错误信息。使用FLOTRAN单元要注意如下几点(下面所列命令相应的菜单路径请参见ANSYS命令手册或联机帮助中的“CommandsandTheirLocationintheGUI”FLOTRAN单元不能和其他单元联合使用。节点坐标系必须与总体坐标系一致。/CLEAR命令并不破坏业已存在的FLOTRAN结果文件(Jobname.RFL),这有助于防止用户不小心破坏那些花了很多时间和精力才求得的结果,用户必须在操作系统里才能删除那些无用的结果文件。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-27,FLOTRAN单元使用的一些限制(续),T2-2.FLOTRAN单元的特点CP命令通过对自由度进行耦合来形成周期边界条件,ANSYS命令手册对CP命令的描述是可以只对某些自由度进行耦合,但作FLOTRAN分析时,周期边界的所有自由度都将被耦合。用户不能对同一个单元中的节点进行耦合,相邻单元间节点耦合也很困难。ADAPT命令不适用于FLOTRAN分析。不能用ANTYPE命令来引入FLOTRAN的瞬态分析。FLOTRAN分析不支持自动时间步长功能,详见“FLOTRAN瞬态分析”。如果用户通过BFCUM、BFDELE或BFUNIF来定义节点热源,则ANSYS会在内部用BFE命令来代替。不能使用LDREAD,FORC命令来电磁载荷转换到FLOTRAN分析中,而必须使用相应的宏来进行转换。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-28,FLOTRAN单元使用的一些限制(续),T2-2.FLOTRAN单元的特点FLOTRAN分析不能使用CE、CECMOD、CEDELE和DEINTF命令。FLOTRAN分析不能使用CNVTOL命令来设置收敛容差。不能用DSYM命令来定义FLOTRAN的对称和反对称边界条件。FLOTRAN不支持旋转坐标系中的角加速度向量。在FLOTRAN分析中,使用FLDATA4,TIME命令而不是DELTIM命令来定义一个载荷步的时间。对于FLOTRAN分析,不能使用DESOL命令或PRESOL命令来修改节点的热(HEAT)、流(FLOW)、或流密(FLUX)结果。FLOTRAN分析不允许将积分点结果外推到节点上(ERESX命令)。FLOTRAN分析不允许通过镜象操作来形成单元。FLOTRAN分析不能用KBC命令来施加渐变载荷,而必须用多个载荷步来逐渐改变载荷。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-29,FLOTRAN单元使用的一些限制(续),T2-2.FLOTRAN单元的特点FLOTRAN分析不允许用LCCALC、LCDEF、LCFA、LCFI等命令来作不同载荷状况之间的运算操作。NCNV命令中的收敛工具不能用于FLOTRAN中相互独立的求解器。FLOTRAN分析不允许用NEQUIT命令来定义非线性分析的平衡迭代数。FLOTRAN分析要求节点坐标系必须是总体的卡尔坐标系,故不能用N命令、NMODIF命令、和NROTAT命令的旋转(rotational)域。FLOTRAN分析中,FLDATA2,ITER命令和FLDATA4,TIME命令用来控制一个载荷步中的总体迭代数。FLOTRAN分析不允许用户自定义单元。FLOTRAN分析中,FLDATA2,ITER命令和FLDATA4,TIME命令用来对写入数据库中的结果进行控制。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-30,FLOTRAN单元使用的一些限制(续),T2-2.FLOTRAN单元的特点PRNLD命令不能用于FLOTRAN分析中,因其不能将边界条件作为可打印的单元节点载荷来保存。FLOTRAN分析中没有节点反力解。部分和预定义求解选项(PSOLVE命令定义)不适用于FLOTRAN各自独立的求解器。TIME命令不能用在FLOTRAN分析中。FLOTRAN用FLDATA1,SOLU命令而不是TIMINT命令来定义瞬态载荷步。FLOTRAN用FLDATA4,TIME命令而非TRNOPT命令来定义瞬态分析选项。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-31,三、FLOTRAN分析的主要步骤,T2-3.FLOTRAN分析的主要步骤一个典型的FLOTRAN分析有如下七个主要步骤:1.确定问题的区域。2.确定流体的状态。3.生成有限元网格。4.施加边界条件。5.设置FLOTRAN分析参数。6.求解。7.检查结果。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-32,第一步:确定问题的区域,T2-3.FLOTRAN分析的主要步骤用户必须确定所分析问题的明确的范围,将问题的边界设置在条件已知的地方,如果并不知道精确的边界条件而必须作假定时,就不要将分析的边界设在靠近感兴趣区域的地方,也不要将边界设在求解变量变化梯度大的地方。有时,也许用户并不知道自己的问题中哪个地方梯度变化最大,这就要先作一个试探性的分析,然后再根据结果来修改分析区域。这些在后面章节中都有详述。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-33,第二步:确定流体的状态,T2-3.FLOTRAN分析的主要步骤用户在此需要估计流体的特征,流体的特征是流体性质、几何边界以及流场的速度幅值的函数。FLOTRAN能求解的流体包括气流和液流,其性质可随温度而发生显著变化,FLOTRAN中的气流只能是理想气体。用户须自己确定温度对流体的密度、粘性、和热传导系数的影响是否是很重要,在大多数情况下,近似认为流体性质是常数,即不随温度而变化,都可以得到足够精确的解。通常用雷诺数来判别流体是层流或紊流,雷诺数反映了惯性力和粘性力的相对强度。通常用马赫数来判别流体是否可压缩。流场中任意一点的马赫数是该点流体速度与该点音速之比值,当马赫数大于0.3时,就应考虑用可压缩算法来进行求解;当马赫数大于0.7时,可压缩算法与不可压缩算法之间就会有极其明显的差异。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-34,第三步:生成有限元网格,T2-3.FLOTRAN分析的主要步骤用户必须事先确定流场中哪个地方流体的梯度变化较大,在这些地方,网格必须作适当的调整。例如:如果用了紊流模型,靠近壁面的区域的网格密度必须比层流模型密得多,如果太粗,该网格就不能在求解中捕捉到由于巨大的变化梯度对流动造成的显著影响,相反,那些长边与低梯度方向一致的单元可以有很大的长宽比。为了得到精确的结果,应使用映射网格划分,因其能在边界上更好地保持恒定的网格特性,映射网格划分可由命令MSHKEY,1或其相应的菜单MainMenuPreprocessor-Meshing-Mesh-entity-Mapped来实现。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-35,第四步:施加边界条件,T2-3.FLOTRAN分析的主要步骤可在划分网格之前或之后对模型施加边界条件,此时要将模型所有的边界条件都考虑进去,如果与某个相关变量的条件没有加上去,则该变量沿边界的法向值的梯度将被假定为零。求解中,可在重启动之间改变边界条件的值,如果需改变边界条件的值或不小心忽略了某边界条件,可无须作重启动,除非该改变引起了分析的不稳定。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-36,第五步:设置FLOTRAN分析参数,T2-3.FLOTRAN分析的主要步骤为了使用诸如紊流模型或求解温度方程等选项,用户必须激活它们。诸如流体性质等特定项目的设置,是与所求解的流体问题的类型相关的,该手册的其他部分详细描述了各种流体类型的所建议的参数设置。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-37,第六步:求解,T2-3.FLOTRAN分析的主要步骤通过在观察求解过程中相关变量的改变率,可以监视求解的收敛性及稳定性。这些变量包括速度、压力、温度、动能(ENKE自由度)和动能耗散率(ENDS自由度)等紊流量以及有效粘性(EVIS)。一个分析通常需要多次重启动。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-38,第七步:检查结果,T2-3.FLOTRAN分析的主要步骤可对输出结果进行后处理,也可在打印输出文件里对结果进行检查,此时用户应使用自己的工程经验来估计所用的求解手段、所定义的流体性质、以及所加的边界条件的可信程度。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-39,四、FLOTRAN分析中产生的文件,T2-4.FLOTRAN分析中产生的文件在ANSYS中进行的大多数流体分析都是通过多次中断和重启动来完成的,通常,分析人员需要在各个重启动之间改变诸如松弛系数等参数或开关某些项(如求解温度方程的开关)。每当用户继续一个分析时,ANSYS程序会自动将数据附加在所有的由FLOTRAN单元产生的文件中。下面将对FLOTRAN单元产生的所有文件进行说明:结果文件,Jobname.RFL,包含节点结果。打印文件,Jobname.PFL,包含各量的收敛记录及进/出口状态(如流量等)。壁面文件,Jobname.RSW,包含壁面剪切应力以及Y-Plus信息。残差文件,Jobname.RDF,包含节点残差。调试文件,Jobname.DBG,包含数学求解器的有关信息。结果备份文件,Jobname.RFO,包含结果文件数据的一个拷贝。重启动文件,Jobname.CFD,包含FLOTRAN的数据结构。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-40,结果文件,T2-4.FLOTRAN分析中产生的文件FLOTRAN分析的结果不自动保存,在每次求解之后,程序会将一个结果集附加在结果文件Jobname.RLF中。用户可对结果文件的更新频率进行控制,ANSYS命令手册中对FLDATA5,OUTP命令的介绍就详细说明了结果文件会基于用户的选择而保存些什么内容。在稳态分析中,在求解的初期多保存几个结果可以比较各结果集之间的变化、可以使用不同的选项或松弛系数来从一个分析的较早状态重新开始分析。当开始一个新分析时(第一次迭代之前),ANSYS程序会保存一个结果,然后在当中断发生时保存再保存结果,在这些事件之间,用户还可通过设置将一些中间结果附在结果文件里,这样就可以从较早的分析状态开时,通过激活一些不同的选项和特征来重新分析。设置覆盖频率的方式如下:命令:FLDATA2,ITER,OVER,value菜单:MainMenuSolutionFLOTRANSetUpExecutionCtrl设置附加频率的方式如下:命令:FLDATA2,ITER,APPE,value菜单:MainMenuSolutionFLOTRANSetUpExecutionCtrl,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-41,打印文件,T2-4.FLOTRAN分析中产生的文件打印文件Jobname.PFL文件包含了所有FLOTRAN输入参数的完整记录,该信息每在发出一个求解命令时就保存一次以完整地记录整个分析历程。同时,所有激活了的变量的收敛过程也记录了下来,还有一个对结果的总结,即每个性质和自由度的最大最小值,这些记录的频率都由用户自己设定。所记录的其他量还有:各记录量的平均值、质量流的边界、质量平衡的计算、所有热传导和热源的相关信息。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-42,节点残差文件,T2-4.FLOTRAN分析中产生的文件节点残差文件Jobname.RDF,显示了当前解的收敛好坏程度。在求解过程的每一个阶段,流场、性质场、温度场都用于对每个自由度计算系数矩阵和强迫函数,如果解完全收敛,这些矩阵和强迫函数将会生成一个与产生它们的速度场一样的速度场,同时,矩阵方程的残差也会变得很小。要得到一个残差文件,必须至少执行一次迭代。当求解过程发生振荡时,残差的幅值将显示分析的错误所在。(矩阵的主对角元素对残差作归一化处理)这种归一化使用户可对自由度的值及其残差作比较。对每一个激活了的自由度计算残差并将其存入残差文件的方式如下:命令:FLDATA5,OUTP,RESI,TRUE菜单:MainMenuSolutionFLOTRANSetUpAdditionalOutResidualFile要读取残差文件,可通过菜单MainMenuGeneralPostprocFLOTRAN2.1A或命令FLREAD来实现。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-43,重启动文件,T2-4.FLOTRAN分析中产生的文件通常,FLOTRAN在一个重启动的起始处计算数据结构,对于大模型,这种计算将消耗大量的时间。为了避免它,可要求FLOTRAN将数据结构保存在重启动文件Jobname.CFD中,FLOTRAN从ANSYS的数据库中产生该文件。对Jobname.CFD文件的读和写的方式如下:命令:FLDATA32,REST,RFIL,T菜单:MainMenuPreprocessorFLOTRANSetUpRestartOptionsCFDRestartFile可将RFIL状态设置为ON或OFF,若设为ON,则FLOTRAN开始执行分析时将读入重启动文件,若文件不存在,将产生一个重启动文件。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-44,重启动文件(续),T2-4.FLOTRAN分析中产生的文件如果在改变了边界条件之后再进行重启动分析,则必须覆盖掉已存在的.CFD文件,以使得ANSYS能用新的边界条件进行重新分析,覆盖.CFD文件的方式如下:命令:FLDATA32,REST,WFIL,T菜单:MainMenuPreprocessorFLOTRANSetUpRestartOptionsCFDRestartFile这就使FLOTRAN在下一载荷步产生一个新的重启动文件,并自动将RFIL状态设置为关闭。当新的重启动文件产生之后,用FLDATA32,REST,RFIL,T命令使随后的重启动能使用新的重启动文件。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-45,FLOTRAN重启动分析(续算),T2-4.FLOTRAN分析中产生的文件用户可在结果文件Jobname.RFL中任意一个解集的基础上开始一个重启动分析,重启动位置的设置可基于解集号(NSET)、迭代数(ITER)、载荷步/子步号(LSTP)或瞬态分析的时间(TIME),方式如下:命令:FLDATA32,REST,lable,value其中,lable为上面的NSET、ITER、LSTP、TIME等菜单:MainMenuPreprocessorFLOTRANSetUpRestartOptionsRestart/Iteratio(或Restart/Loadstep,Restart/Set,等)当重启动一个分析时,ANSYS将原始的结果文件拷贝到Jobname.RSO中并将重启动点、所有在重启动点之前的结果集、所有的后续结果集放在新的结果文件中。如果在FLDATA32,REST命令中的value值是一个负值,则将不产生Jobname.RSO文件,而重启动的点将由value的绝对值来指定。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-46,五、提高收敛和稳定性的常用工具,T2-5.提高收敛和稳定性的常用工具ANSYS程序提供几个有助于收敛和求解稳定的工具,理论手册对其机理有详述。松弛系数惯性松弛人工粘性速度限值面积积分阶次,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-47,松弛系数,T2-5.提高收敛和稳定性的常用工具松弛系数是一个其值介于0和1之间的小数,它表示旧结果与附加在旧结果上以形成新结果的最近一次计算量之间的变化量。设置松弛系数的方式如下:命令:FLDATA25,RELX,lable,value菜单:MainMenuPreprocessorFLOTRANSetUpRelax/Stab/CapDOFRelaxationMainMenuPreprocessorFLOTRANSetUpRelax/Stab/CapPropRelaxationMainMenuSolutionFLOTRANSetUpRelax/Stab/CapDOFRelaxationMainMenuSolutionFLOTRANSetUpRelax/Stab/CapPropRelaxation注:命令手册中对该命令的自由度和性质量有详述。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-48,惯性松弛,T2-5.提高收敛和稳定性的常用工具对某个自由度的方程组的惯性松弛就是使其矩阵的主对角占优以保持求解的稳定性。如果当一个解在收敛过程中没有发生舍入误差,则惯性释放的值不会影响到求解的最终结果。但是通常的求解过程都会发生舍入误差,故惯性松弛可能对结果产生影响。用户可对动量方程(MOME)、紊流方程(TURB)。压力方程(PRES)和温度方程(TEMP)施加惯性松弛,其方式如下:命令:FLDATA26,STAB,lable,value菜单:MainMenuPreprocessorFLOTRANSetUpRelax/Stab/CapStabilityParmsMainMenuSolutionFLOTRANSetUpRelax/Stab/CapStabilityParms惯性松弛系数是以所加项的分母的形式出现的,故其值越小,所起作用越大,其典型值介于1.0(作用中等)到1.0107(作用很大)之间。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-49,人工粘性,T2-5.提高收敛和稳定性的常用工具人工粘性用于在梯度较大的区域平抑速度解。它有助于可压缩问题的收敛,也有助于对有分布阻力的不可压缩问题的速度解进行平抑。对于不可压缩问题,应使人工粘性的幅值与有效粘性的幅值处于相同的数量级。施加人工粘性的方式如下:命令:FLDATA26,STAB,VISC,value菜单:MainMenuPreprocessorFLOTRANSetUpRelax/Stab/CapStabilityParmsMainMenuSolutionFLOTRANSetUpRelax/Stab/CapStabilityParms,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-50,速度限值,T2-5.提高收敛和稳定性的常用工具速度限值使所求解量不能超出用户所定义的值,可对速度、压力和温度自由度进行限制(VX、VY、VZ、PRES、TEMP),方式如下:命令:FLDATA31,CAPP菜单:MainMenuPreprocessorFLOTRANSetUpRelax/Stab/CapResultsCappingMainMenuSolutionFLOTRANSetUpRelax/Stab/CapResultsCapping速度限值可消除速度尖峰的不利影响,这种速度尖峰通常发生在收敛过程中的较早阶段。它还特别适合用于可压缩流分析,因这类分析中速度尖峰通常使动能项大到产生负的静温。当对压力进行限值时,所限的值是由压力方程解算出来的压力而不是松弛后的压力,故当限值后作重启动时,压力值仍有可能超出限值。注意:当有速度限值时,质量有可能不守恒。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-51,面积积分阶次,T2-5.提高收敛和稳定性的常用工具面积积分阶次(QuadratureOrder)缺省的用于计算单元面积积分的阶次是单点积分,用户可对其进行控制。对于轴对称问题,求解时,该值自动设为2,因为当面积积分阶次为2时,可使含有异常形状单元的问题收敛到更精确的解。用下面的方式改变动量、压力、热或紊流项的面积积分阶次:命令:FLDATA30,QUAD,lable,value其中,lable为要改变的单元积分,value为积分点的数目。菜单:MainMenuPreprocessorFLOTRANSetUpModRes/QuadOrdCFDQuadOrdersMainMenuSolutionFLOTRANSetUpModRes/QuadOrdCFDQuadOrders,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-52,六、分析过程中应处理的问题,T2-6.FLOTRAN分析过程中应处理的问题ANSYS程序的FLOTRAN分析过程中应处理的问题主要有以下几个选项,理论手册对其机理有详述。确定总体迭代的数目收敛监测中断一个FLOTRAN分析,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-53,确定总体迭代的数目,T2-6.FLOTRAN分析过程中应处理的问题FLOTRAN分析是一个非线性的序列求解过程,故每次分析首先得确定要让程序执行多少次迭代。一次总体迭代就是对所有相关的控制方程按序列进行求解,并且在求解过程中流体性质会随时更新。在瞬态分析中,时间步循包含了总体迭代循环。在一个总体迭代中,程序首先获得动量方程的近似解,再在质量守恒的基础上将动量方程的解作为强迫函数来求解压力方程,然后用压力解来更新速度,以使速度场保持质量守恒。如果要求了程序求解温度,则程序会同时求解温度方程并更新与温度相关的流体性质。最后,如果激活了紊流模型,则程序将求解紊流方程并用紊流动能及其耗散率来计算有效粘性和热传导系数,有效粘性和热传导系数将分别代替层流粘性和热传导系数以在平均流上模拟紊流的影响。用下面的方式定义总体迭代的数目:命令:FLDATA2,ITER,EXEC,value(value即为迭代数)菜单:MainMenuPreprocessorFLOTRANSetUpExecutionCtrlMainMenuSolutionFLOTRANSetUpExecutionCtrl,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-54,收敛监测,T2-6.FLOTRAN分析过程中应处理的问题在FLOTRAN求解过程中,程序在每一个总体迭代里对每一个自由度计算出一个收敛监测量,这些自由度包括:速度(VX、VY、VZ)、压力(PRES)、温度(TEMP)、紊流动能(ENKE)、动能耗散率(ENDS)、以及激活了的多组份传输方程(SP01SP06)。收敛监测量就是两次迭代之间结果改变量的归一化值,若以表示任一自由度,则该自由度的收敛监测量可由下式表示:,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-55,收敛监测(续),T2-6.FLOTRAN分析过程中应处理的问题收敛监测量表示变量在当前迭代(kth)的结果和前一次迭代((k-1)th)的结果之间差值的总和除以当前值的总和,这种求和是在所有节点上进行的,并且使用的是差值的绝对值。在批处理或交互式运行过程中,当求解进行时,程序的“图形求解跟踪(GST)”功能将实时显示出所计算的收敛监测量,GST的缺省值在交互运行时是开(ON),而在批处理运行时是关(OFF)。用户可用下面的方式定义其开关:命令:/GST菜单:MainMenuSolutionOutputCtrlsGrphSoluTrack图21是两个典型的GST图形。图21(b)是一个FLOTRAN的瞬态分析过程,图中的每一个尖峰表示了一个新时间步的开始。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-56,收敛监测(续),T2-6.FLOTRAN分析过程中应处理的问题,Objective,由GST显示的收敛监测量(a)稳态求解,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-57,收敛监测(续),T2-6.FLOTRAN分析过程中应处理的问题,Objective,由GST显示的收敛监测量(a)瞬态求解,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-58,收敛监测(续),T2-6.FLOTRAN分析过程中应处理的问题在初始阶段可能出现的一些振荡之后,收敛监测量的大小将随着分析过程的收敛而逐渐减小,但其减小的程度将依赖于几个因素,诸如:几何边界的复杂程度高梯度区域有限元网格的精度紊流的严重程度(由雷诺数确定)出口边界处流场的发展是否充分当使用图形求解跟踪(GST)功能时,还应注意:不单是FLOTRAN分析有GST功能,非线性的结构分析、非线性的热分析和非线性的电磁场分析都有GST功能。详见各自的分析指导手册。GST可同时显示多达10条的跟踪曲线,如果用户的模型有多于10个的自由度,则GST将只显示前10个自由度的收敛跟踪曲线。当GST开始显示时,程序会弹出一个带STOP按钮的对话框,用户可在任意时刻通过点取该STOP按钮来中断求解过程,而后要进行重启动分析时,可通过执行命令SOLVE或其相应的菜单MainMenuSolutionRunFlotran来实现。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-59,中断一个FLOTRAN分析,T2-6.FLOTRAN分析过程中应处理的问题用户可以定义一个基于压力和温度收敛监测量的目标值来中断一个FLOTRAN分析,定义方式如下:命令:FLDATA3,TERM,PRES,valueFLDATA3,TERM,TEMP,value菜单:MainMenuPreprocessorFLOTRANSetUpExecutionCtrlMainMenuSolutionFLOTRANSetUpExecutionCtrl压力和温度的收敛缺省值都是1.010-8,如果没有激活温度方程的求解,则程序只检测压力的收敛值是否满足要求。若同时激活了流体和温度方程的求解,则二者的收敛标准都必须同时满足。在满足了压力和温度的收敛条件或总体平衡迭代数达到了所要求的值后,FLOTRAN求解过程就自动中断。要中断一个FLOTRAN分析,需在当前工作目录下生成一个Jobname.ABT文件,该文件的第一行应含有terminate字样。在每一次总体迭代之前,FLOTRAN都会在当前目录下搜寻Jobname.ABT文件,如果程序找到该文件并发现含有terminate,则立即完成该次总体迭代并正常中断程序的执行,且将结果写入结果文件。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-60,七、对一个分析进行评价,T2-7.对一个FLOTRAN分析进行评价分析员必须回答的两个问题是:1.所作的分析是何时结束的?2.所作的分析是否是正确的?这两个问题是相互关联的,因为,如果没有正确地设置和正确地分析一个流体问题,它一般都是不会收敛的。如果所输入的初始参数和所有的边界条件都是正确的,则当所有变量的收敛监测量都停止增长,以及所有求解量的平均、最大、最小值都不再升降时,求解过程就算是完成了。然而,这并不能保证所求解的结果是唯一正确的,因为自然界本身并不保证存在唯一解。振荡问题(例如:柱体绕流的旋涡脱落问题)用稳态或瞬态求解技术都不能得到一个稳定的解。要验证一个分析是稳定的或是振荡的,可以通过对它执行大量的迭代求解来实现。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-61,七、对一个分析进行评价,T2-7.对一个FLOTRAN分析进行评价ANSYS将求解变量的平均、最小、最大值保存在文件Jobname.RFL中,该文件同时还保存了FLOTRAN的输入数据和计算出的收敛监测量、所有自由度的结果总结、层流特性和有效特性。可用下面的方式来规定ANSYS进行结果总结的频率:命令:FLDATA5,OUTP,SUMF,value菜单:MainMenuPreprocessorFLOTRANSetUpAdditionalOutRFLOutDerivedMainMenuSolutionFLOTRANSetUpAdditionalOutRFLOutDerived,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-62,八、验证结果,T2-7.验证结果验证求解结果的可靠性是所有分析人员的责任,如果一个FLOTRAN分析得到了非预期的结果,则应进行下列所示的一些操作,这些操作的大部分都可以在开始一个分析前完成。即使只进行了零次迭代,ANSYS也会生成一个Jobname.RFL文件并检查所有的输入数据。1.检查作为结果总结的一部分而打印出来的质量平衡情况。内部检查将确定是否有任何的可能会通过模型的质量流,允许质量流的边界条件是:确定的速度边界条件确定的压力边界条件未确定的边界(这有可能是由于用户忘了施加边界条件而致)ANSYS会将进口和出口边界编号列表,而这些应与所希望的条件相对应。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-63,八、验证结果,T2-7.验证结果2.在ANSYS里检查边界条件,以保证其正确性。3.检查所定义的流体性质及其随温度的可变动性正确与否,这可在.RFL文件中方便地检查。4.检查用以建立模型的单位制与用以定义流体性质的单位制是否一致。5.有时,还需确认与所选选项相联系的方程的求解是否正确(例如:可压缩流中的压力方程)。,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-64,八、验证结果,T2-7.验证结果6.如果求解发散,可能的原因还有:有限元网格不够精细、或者邻近出口处流场梯度太大,要解决这些问题,可以使用一些诸如惯性松弛等有助于收敛的手段,本手册的后面将详述各种松弛技术。7.如果仅仅只有某个特定的量产生发散,则可将该量重新初始化到一个单值,并作重启动分析,方式如下:命令:FLDATA29,MODV菜单:MainMenuPreprocessorFLOTRANSetUpModRes/QuadOrdModifyResultsMainMenuSolutionFLOTRANSetUpModRes/QuadOrdModifyResults,Objective,2001年10月1日,ANSYS培训教程版本5.5XJTUMSSVBy:HaichGao(011001),Intro-65,第四章FLOTRAN边界条件,2001年10月1日,ANSYS培训教程版本5.5XJTUMS
展开阅读全文