ANSYS - CFD之Flotran中文讲解说明（全+重点标注）11

目录

第一章 FLOTRAN 计算流体动力学(CFD)分析概述 ............................................................ 2 FLOTRAN CFD 分析的概念 ........................................................................................................ 2 FLOTRAN 分析的种类 ................................................................................................................. 2 FLOTRAN单元的特点 .................................................................................................................. 4 FLUID141单元 ............................................................................................................................... 4 FLUID142单元 ............................................................................................................................... 4 FLUID单元的其他特征 ................................................................................................................. 4 FLOTRAN单元使用中的一些限制 .............................................................................................. 5 FLOTRAN分析的主要步骤 .......................................................................................................... 6 FLOTRAN分析中产生的一些文件 .............................................................................................. 8 FLOTRAN重启动分析（续算） .................................................................................................. 9 Restart/Iteratio（或Restart/Load step, Restart/Set, 等） .................................................. 10 ANSYS程序提供几个有助于收敛和求解稳定的工具，理论手册对其机理有详述。 ........... 10 Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Prop Relaxation ................ 10 Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms ......................... 10 Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms ......................... 11 Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Results Capping ........................ 11 FLOTRAN分析过程中应处理的问题 ........................................................................................ 11 Main Menu>Solution>FLOTRAN SetUp>Execution Ctrl ........................................................ 12 FLDATA3,TERM,TEMP,value ....................................................................................................... 14 Main Menu>Solution>FLOTRAN SetUp>Execution Ctrl ........................................................ 14 FLOTRAN设置命令 ..................................................................................................................... 15 FLOTRAN边界条件 ..................................................................................................................... 56 第五章 FLOTRAN层流和湍流分析算例 ............................................................................. 61

第一章 FLOTRAN 计算流体动力学(CFD)分析概述

FLOTRAN CFD 分析的概念

ANSYS程序中的FLOTRAN CFD分析功能是一个用于分析二维及三维流体流动场的先进的工具，使用ANSYS中用于FLOTRAN CFD分析的FLUID 141和FLUID 142 单元，可解决如下问题：

? 作用于气动翼(叶)型上的升力和阻力 ? 超音速喷管中的流场

? 弯管中流体的复杂的三维流动 同时，FLOTRAN还具有如下功能：

? 计算发动机排气系统中气体的压力及温度分布 ? 研究管路系统中热的层化及分离

? 使用混合流研究来估计热冲击的可能性

? 用自然对流分析来估计电子封装芯片的热性能 ? 对含有多种流体的(由固体隔开)热交换器进行研究

FLOTRAN 分析的种类

FLOTRAN可执行如下分析： ? 层流或紊流 ? 传热或绝热

? 可压缩或不可压缩 ? 牛顿流或非牛顿流 ? 多组份传输

这些分析类型并不相互排斥，例如，一个层流分析可以是传热的或者是绝热的，一个紊流分析可以是可压缩的或者是不可压缩的。

层流分析

层流中的速度场都是平滑而有序的，高粘性流体（如石油等）的低速流动就通常是层流。

紊流分析

紊流分析用于处理那些由于流速足够高和粘性足够低从而引起紊流波动的流体流动情况，ANSYS中的二方程紊流模型可计及在平均流动下的紊流速度波动的影响。如果流体的密度在流动过程中保持不变或者当流体压缩时只消耗很少的能量，该流体可认为是不可压缩的，不可压缩流的温度方程将忽略流体动能的变化和粘性耗散。

热分析

流体分析中通常还会求解流场中的温度分布情况。如果流体性质不随温度而变，就可不解温度方程。在共轭传热问题中，要在同时包含流体区域和非流体区域（即固体区域）的整个区域上求解温度方程。在自然对流传热问题中，流体由于温度分布的不均匀性而导致流体密度分布的不均匀性，从而引起流体的流动，与强迫对流问题不同的是，自然对流通常都没有外部的流动源。

可压缩流分析

对于高速气流，由很强的压力梯度引起的流体密度的变化将显著地影响流场的性质，ANSYS对于这种流动情况会使用不同的解算方法。

非牛顿流分析

应力与应变率之间成线性关系的这种理论并不能足以解释很多流体的流动，对于这种非牛顿流体，ANSYS程序提供了三中粘性模式和一个用户自定义子程序。

多组份传输分析

这种分析通常是用于研究有毒流体物质的稀释或大气中污染气体的传播情况，同时，它也可用于研究有多种流体同时存在（但被固体相互隔开）的热交换器分析。

第二章 FLOTRAN分析基础

FLOTRAN单元的特点

ANSYS中的FLOTRAN单元，即FLUID141和FLUID142，用于解算单相粘性流体的二维和三维流动、压力和温度分布。对于这些单元，ANSYS通过质量、动量和能量三个守恒性质来计算流体的速度分量、压力、以及温度。

FLUID141单元

FLUID141单元具有下列特征： 维数：二维 形状：四节点四边形或三节点三角形 自由度：速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额

FLUID142单元

FLUID142单元具有下列特征： 维数：三维 形状：四节点四面体或八节点六面体 自由度：速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额

FLUID141单元

FLUID142单元

FLUID单元的其他特征

FLUID单元的其他特征包括：

? 用于模拟紊流的二方程紊流模式

? 有很多推导结果，诸如：流场分析中的马赫数、压力系数、总压、剪应

力、壁面处的y-plus、以及流线函数；热分析中的热流、热交换（膜）系数等。

? 流体边界条件，包括：速度、压力、紊流动能以及紊流能量耗散率。用

户无需提供流场进口处紊流项的边界条件，因FLOTRAN对此提供的缺省值适用于绝大多数分析。

? 热边界条件，包括：温度、热流、体积热源、热交换（膜）系数。 用户可使用的坐标系有：的卡尔坐标系、柱坐标系、极坐标系和轴对称坐标系。如果所计算的问题是轴对称的，激活旋转（swirl）选项即可算出垂直于对称平面的速度分量。

使用FLOTRAN单元的一些限制及注意事项

FLOTRAN单元的一些局限性：

? 在同一次分析中不能改变求解的区域 ? 单元不支持自由流面边界条件

? ANSYS程序的某些特征不能同FLOTRAN单元一起使用 ? 使用FLOTRAN单元时不能使用某些命令或菜单

? 当使用ANSYS的图形用户界面时，程序将只能显示那些在菜单和对话

框中的FLOTRAN SetUp部分要求了的特征和选项。

FLOTRAN单元使用中的一些限制

当使用FLOTRAN单元时，要避免使用ANSYS的某些特征和命令，至少，要注意到在使用FLOTRAN单元时与别的分析稍微有些不同，当使用了无效的命令时，程序会给出相应的警告或错误信息。使用FLOTRAN单元要注意如下几点（下面所列命令相应的菜单路径请参见ANSYS命令手册或联机帮助中的“Commands and Their Location in the GUI”

? FLOTRAN单元不能和其他单元联合使用。 ? 节点坐标系必须与总体坐标系一致。

? /CLEAR命令并不破坏业已存在的FLOTRAN结果文件（Jobname.RFL），

这有助于防止用户不小心破坏那些花了很多时间和精力才求得的结果，用户必须在操作系统里才能删除那些无用的结果文件。

? CP命令通过对自由度进行耦合来形成周期边界条件，ANSYS命令手册

对CP 命令的描述是可以只对某些自由度进行耦合，但作FLOTRAN分析时，周期边界的所有自由度都将被耦合。

用户不能对同一个单元中的节点进行耦合，相邻单元间节点耦合也很困难。

? ADAPT命令不适用于FLOTRAN分析。

? 不能用ANTYPE命令来引入FLOTRAN的瞬态分析。 ? FLOTRAN分析不支持自动时间步长功能，详见“FLOTRAN瞬态分析”。 ? 如果用户通过BFCUM、BFDELE或BFUNIF来定义节点热源，则

ANSYS会在内部用BFE命令来代替。

? 不能使用LDREAD,FORC命令来电磁载荷转换到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命令来施加渐变载荷，而必须用多个载荷步

? ? ? ? ? ? ? ? ? ? ? ? ?

来逐渐改变载荷。

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命令用来对写入数据库中的结果进行控制。

PRNLD命令不能用于FLOTRAN分析中，因其不能将边界条件作为可打印的单元节点载荷来保存。

FLOTRAN分析中没有节点反力解。

部分和预定义求解选项（PSOLVE命令定义）不适用于FLOTRAN各自独立的求解器。

TIME命令不能用在FLOTRAN分析中。

FLOTRAN用FLDATA1,SOLU命令而不是TIMINT命令来定义瞬态载荷步。

FLOTRAN用FLDATA4,TIME命令而非TRNOPT命令来定义瞬态分析选项。

FLOTRAN分析的主要步骤

一个典型的FLOTRAN分析有如下七个主要步骤： 1. 确定问题的区域。 2. 确定流体的状态。 3. 生成有限元网格。 4. 施加边界条件。

5. 设置FLOTRAN分析参数。 6. 求解。 7. 检查结果。

第一步：确定问题的区域

用户必须确定所分析问题的明确的范围，将问题的边界设置在条件已知的地方，如果并不知道精确的边界条件而必须作假定时，就不要将分析的边界设在靠近感兴趣区域的地方，也不要将边界设在求解变量变化梯度大的地方。有时，也许用户并不知道自己的问题中哪个地方梯度变化最大，这就要先作一个试探性的分析，然后再根据结果来修改分析区域。这些在后面章节中都有详述。

第二步：确定流体的状态

用户在此需要估计流体的特征，流体的特征是流体性质、几何边界以及流场的速度幅值的函数。FLOTRAN能求解的流体包括气流和液流，其性质可随温度

而发生显著变化，FLOTRAN中的气流只能是理想气体。用户须自己确定温度对流体的密度、粘性、和热传导系数的影响是否是很重要，在大多数情况下，近似认为流体性质是常数，即不随温度而变化，都可以得到足够精确的解。 通常用雷诺数来判别流体是层流或紊流，雷诺数反映了惯性力和粘性力的相对强度，详见第四章。

通常用马赫数来判别流体是否可压缩，详见第七章。流场中任意一点的马赫数是该点流体速度与该点音速之比值，当马赫数大于0.3时，就应考虑用可压缩算法来进行求解；当马赫数大于0.7时，可压缩算法与不可压缩算法之间就会有极其明显的差异。

第三步：生成有限元网格

用户必须事先确定流场中哪个地方流体的梯度变化较大，在这些地方，网格必须作适当的调整。例如：如果用了紊流模型，靠近壁面的区域的网格密度必须比层流模型密得多，如果太粗，该网格就不能在求解中捕捉到由于巨大的变化梯度对流动造成的显著影响，相反，那些长边与低梯度方向一致的单元可以有很大的长宽比。 为了得到精确的结果，应使用映射网格划分，因其能在边界上更好地保持恒定的网格特性，映射网格划分可由命令MSHKEY,1或其相应的菜单Main Menu>Preprocessor > -Meshing-Mesh>-entity-Mapped来实现。

第四步：施加边界条件

可在划分网格之前或之后对模型施加边界条件，此时要将模型所有的边界条件都考虑进去，如果与某个相关变量的条件没有加上去，则该变量沿边界的法向值的梯度将被假定为零。求解中，可在重启动之间改变边界条件的值，如果需改变边界条件的值或不小心忽略了某边界条件，可无须作重启动，除非该改变引起了分析的不稳定。

第五步：设置FLOTRAN分析参数

为了使用诸如紊流模型或求解温度方程等选项，用户必须激活它们。诸如流体性质等特定项目的设置，是与所求解的流体问题的类型相关的，该手册的其他部分详细描述了各种流体类型的所建议的参数设置。

第六步：求解

通过在观察求解过程中相关变量的改变率，可以监视求解的收敛性及稳定性。这些变量包括速度、压力、温度、动能 (ENKE自由度) 和动能耗散率 (ENDS自由度)等紊流量以及有效粘性（EVIS）。一个分析通常需要多次重启动。

第七步：检查结果

可对输出结果进行后处理，也可在打印输出文件里对结果进行检查，此时用户应使用自己的工程经验来估计所用的求解手段、所定义的流体性质、以及所加的边界条件的可信程度。

FLOTRAN分析中产生的一些文件

在ANSYS中进行的大多数流体分析都是通过多次中断和重启动来完成的，通常,分析人员需要在各个重启动之间改变诸如松弛系数等参数或开关某些项（如求解温度方程的开关）。每当用户继续一个分析时，ANSYS程序会自动将数据附加在所有的由FLOTRAN单元产生的文件中。下面将对FLOTRAN单元产生的所有文件进行说明:

? 结果文件， Jobname. RFL，包含节点结果。 ? 打印文件， Jobname.PFL，包含各量收敛记录及进/出口状态(如流量等)。 ? 壁面文件， Jobname.RSW，包含壁面剪切应力以及Y-Plus信息。 ? 残差文件， Jobname.RDF，包含节点残差。

? 调试文件， Jobname.DBG，包含数学求解器的有关信息。

? 结果备份文件， Jobname.RFO，包含结果文件数据的一个拷贝。 ? 重启动文件， Jobname.CFD，包含FLOTRAN的数据结构。

结果文件

FLOTRAN分析的结果并不自动保存在ANSYS的数据库中，在每次求解之后，程序会将一个结果集附加在结果文件Jobname.RLF中。用户可对结果文件的内容及程序对结果文件的更新频率进行控制，ANSYS命令手册中对FLDATA5,OUTP命令的介绍就详细说明了结果文件会基于用户的选择而保存些什么内容。

在一个稳态FLOTRAN分析中，结果文件能保存多少个结果集是没有限制的，在求解的初期多保存几个结果有很多好处：可以比较各结果集之间的变化、可以使用不同的选项或松弛系数来从一个分析的较早状态重新开始分析。 当开始一个新分析时（在其第一次迭代之前），ANSYS程序会保存一个结果, 然后在当中断发生时再保存结果，在这些事件之间，用户还可通过设置将一些中间结果附在结果文件里，这样就可以从较早的分析状态开时，通过激活一些不同的选项和特征来重新分析，例如，可以通过这种方式来提高分析的稳定性。 使用 ANSYS 的覆盖频率选项是一个明智的方法，它就可以周期性地保存和更新一个临时的结果集，这样，当由于断电或其他系统原因而发生求解中断时，总可以有一个可用的结果集用于重新开始分析。设置覆盖频率的方式如下： 命令：FLDATA2,ITER,OVER,value 菜单：Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl 设置附加频率的方式如下： 命令：FLDATA2,ITER,APPE,value 菜单：Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl

打印文件（ Jobname.PFL）

Jobname.PFL文件包含了所有FLOTRAN输入参数的完整记录，该信息每在发出一个求解命令时就保存一次以完整地记录整个分析历程。同时，所有激活了的变量的收敛过程也记录了下来，还有一个对结果的总结，即每个性质和自由度的最大最小值，这些记录的频率都由用户自己设定。所记录的其他量还有：各记录量的平均值、质量流的边界、质量平衡的计算、所有热传导和热源的相关信息。

节点残差文件

节点残差文件，即Jobname.RDF，显示了当前解的收敛好坏程度。在求解过程的每一个阶段，流场、性质场、温度场都用于对每个自由度计算系数矩阵和强迫函数，如果解完全收敛，这些矩阵和强迫函数将会生成一个与产生它们的速度场一样的速度场，同时，矩阵方程的残差也会变得很小。要得到一个残差文件，必须至少执行一次迭代。 当求解过程发生振荡时，残差的幅值将显示分析的错误所在。（矩阵的主对角元素对残差作归一化处理）这种归一化使用户可对自由度的值及其残差作比较。 对每一个激活了的自由度计算残差并将其存入残差文件的方式如下： 命令：FLDATA5,OUTP,RESI,TRUE 菜单：Main Menu>Solution>FLOTRAN Set Up>Additional Out>Residual File

要读取残差文件，可通过菜单Main Menu>General Postproc>FLOTRAN 2.1A或命令FLREAD来实现。

重启动文件

通常，FLOTRAN在一个重启动的起始处计算数据结构，对于一个大模型，这种计算将消耗大量的时间，为了避免这种重新计算，可要求FLOTRAN将数据结构保存在重启动文件 Jobname.CFD中，FLOTRAN从ANSYS的数据库中产生该文件。 对 Jobname.CFD文件的读和写的方式如下： 命令： FLDATA32,REST,RFIL,T 菜单： Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options>CFD Restart File 可将RFIL状态设置为开(ON)或关(OFF)，若设为开，则FLOTRAN开始执行分析时将读入重启动文件，若此时重启动文件不存在，则将产生一个重启动文件。 如果在改变了边界条件之后再进行重启动分析，则必须覆盖掉业已存在的.CFD 文件以使得ANSYS能用新的边界条件进行重新分析，覆盖.CFD文件的方式如下： 命令： FLDATA32,REST, WFIL,T 菜单： Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options>CFD Restart File

这就使FLOTRAN在下一载荷步产生一个新的重启动文件，并自动将RFIL状态设置为关闭。当新的重启动文件产生之后，用FLDATA32,REST,RFIL,T命令使随后的重启动能使用新的重启动文件。

FLOTRAN重启动分析（续算）

用户可在结果文件 Jobname.RFL中任意一个解集的基础上开始一个重启动分析，重启动位置的设置可基于解集号(NSET)、迭代数(ITER)、载荷步/子步号(LSTP)或瞬态分析的时间(TIME)，方式如下： 命令： FLDATA32,REST,lable,value 其中，lable为上面的NSET、ITER、LSTP、TIME等

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options> Restart/Iteratio（或Restart/Load step, Restart/Set, 等） 当重启动一个分析时，ANSYS将原始的结果文件拷贝到Jobname.RSO中并将重启动点、所有在重启动点之前的结果集、所有的后续结果集放在新的结果文件中。如果在FLDATA32,REST命令中的value值是一个负值，则将不产生Jobname.RSO文件，而重启动的点将由value的绝对值来指定。

提高收敛性和稳定性的常用的工具

ANSYS程序提供几个有助于收敛和求解稳定的工具，理论手册对其机理有详述。

松弛系数

松弛系数是一个其值介于0和1之间的小数，它表示旧结果与附加在旧结果上以形成新结果的最近一次计算量之间的变化量。设置松弛系数的方式如下： 命令：FLDATA25,RELX,lable,value 菜单： Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>DOF

Relaxation

Main Relaxation

Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Prop SetUp>Relax/Stab/Cap>DOF SetUp>Relax/Stab/Cap>Prop

MainMenu>Solution>FLOTRAN Relaxation

Main Menu>Solution>FLOTRAN Relaxation

注：命令手册中对该命令的自由度和性质量有详述。

惯性松弛

对某个自由度的方程组的惯性松弛就是使其矩阵的主对角占优以保持求解的稳定性。如果当一个解在收敛过程中没有发生舍入误差，则惯性释放的值不会影响到求解的最终结果。但是通常的求解过程都会发生舍入误差，故惯性松弛可能对结果产生影响。用户可对动量方程（MOME）、紊流方程（TURB）。压力方程（PRES) 和温度方程（TEMP）施加惯性松弛，其方式如下： 命令：FLDATA26,STAB,lable,value 菜单：Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Stability

Parms Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms

惯性松弛系数是以所加项的分母的形式出现的，故其值越小，所起作用越大，其典型值介于1.0（作用中等）到1.0×10－7（作用很大）之间。

人工粘性

人工粘性用于在梯度较大的区域平抑速度解。它有助于可压缩问题的收敛，也有助于对有分布阻力的不可压缩问题的速度解进行平抑。对于不可压缩问题，应使人工粘性的幅值与有效粘性的幅值处于相同的数量级。施加人工粘性的方式如下： 命令：FLDATA26,STAB,VISC,value

12 设置有助于FLOTRAN求解稳定的参数 13 设定FLOTRAN自由度松弛系数 58

33 34 35 36 38 39 42 45 48 49 55 14 设定FLOTRAN流体性质松弛因子 15 设置FLOTRAN分析的自由度限值 16 选择FLOTRAN各自由度相应的求解器 17 对FLOTRAN各求解器的控制 18 设置FLOTRAN湍流模型的一些常数 19 重新设定FLOTRAN各分析参数的值 20 控制FLOTRAN面积积分的阶次 21 FLOTRAN多组份疏运分析的设置及控制 22 定义FLOTRAN的重启动（续算）控制 23 设置并执行一个零迭代FLOTRAN分析 1 FLOTRAN求解控制命令

命令：FLDATA1, SOLU, Lable, Value 功能：设定求解控制

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Solution Options Main Menu>Solution>FLOTRAN Set Up>Solution Option 其中，Lable的选项及其各自含义如下：

TRAN：求解稳态或非稳态选项，缺省为稳态

FLOW：是否求解流动方程选项，缺省为求解流动方程

TEMP：是否求解温度方程选项，缺省为绝热(不求解温度方程) TURB：层流或湍流选项，缺省为层流

COMP：不可压缩或可压缩选项，缺省为不可压缩流 SWRL：求解轴对称旋流选项，缺省为否 SPEC：求解多组份疏运选项，缺省为否 Value的值为“是”或“否”二者之一。 该命令的菜单形式如下：

2 FLOTRAN执行及输出控制命令?稳态控制参数设置

命令：FLDATA2,ITER, Lable, Value 功能：设置稳态流的迭代及输出控制

其中， Lable的选项及其各自含义如下： EXEC：总体迭代次数，缺省为10次

OVER：结果文件覆盖频率(每隔Value次迭代)，缺省为0 APPE：结果文件附加频率(每隔Value次迭代)，缺省为0 Value的值以上Lable各自的迭代次数 命令：FLDATA3,TERM, Lable, Value 功能：设置压力和温度的收敛准则

其中， Lable的选项及其各自含义如下： PRES：设置压力收敛准则，缺省为10-8 TERM：设置温度收敛准则，缺省为10-8 Value为以上两项各自的收敛值

命令：FLDATA5, OUTP, Lable, Value 功能：设置输出控制

其中， Lable的选项及其各自含义如下：

SUMF：输出总结的频率(每隔Value次迭代)，缺省为10 Value为上面SUMF的迭代数

以上三条命令的菜单路径如下：

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Execution Ctrl Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl 以上三条命令的菜单形式如下：

3 FLOTRAN执行及输出控制命令?瞬态控制参数设置

命令：FLDATA4, TIME, Lable, Value

功能：基于瞬态时间来设置瞬态分析的求解及输出控制 命令：FLDATA4A, STEP, Lable, Value

功能：基于时间步来设置瞬态分析的输出控制

菜单：Main Menu>Preprocessor>FLOTRAN Set Up> Execution Ctrl Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl 在点取上面的菜单路径时，首先会弹出一个下图所示的对话框，该对话框的第二项即是激活上面的第二条命令，即用户需要选择到底是基于瞬态分析的时间值、或是时间步、或是二者兼顾来对求解及输出进行控制。同时，该对话框的第一项是用以定义到底是由程序还是由用户自己来对瞬态分析的时间步长进行控制，也就是下面将要讲到的STEP项。这两个选项的值直接影响了随后的弹出菜单的内容，本手册在此是合起来讲的，谨请注意。

第一条命令的 Lable选项及其各自含义如下：

STEP：定义时间步长，当在上面的对话框中设置为用户自定义时，该选项的 Value值应为一正值，表示真正的时间步长，命令输入就在Value处输一真正时间即可；当在上面的对话框中设置为Advection时，即为由程序定义时间步长，命令输入时可在Value处输-1、-2、-3或-4四值之一，其含义分别如下：

-1：时间步长会小到在单一的任何一个时间步长内，流场中任意一点的 运动距离都不会大于一个单元的长度

-2：只用于可压缩流分析，它使时间步长会小到在单一的任何一个时间 步长内，流场中压力信号的传输距离都不会大于一个单元的长度

-3：只用于可压缩流分析，它使时间步长取上面-1和-2项中的最小值 -4：只用于纯传热分析（不计算流场方程），它使时间步长会小到在单

一的任何一个时间步长内，任何一个“热点”的传导或对流距离不会大于一个单元的长度。

ISTEP：指定初始时间步长，仅在时间步长由程序定义时有效，缺省为0 NUMB：指定时间步的数目，缺省为10

TEND：指定瞬态分析的终止时间，缺省为1×106

GLOB：指定每一时间步的总体迭代数，缺省为20 PRES：指定每一时间步的压力收敛准则，缺省为1×10-6 TEMP：指定每一时间步的温度收敛准则，缺省为1×10-6

OVER：指定结果文件的覆盖频率 (每隔Value时间覆盖一次)，缺省为0，即不覆盖

APPE：指定结果的输出频率(每隔Value时间输出一次结果到结果文件中)， 缺省为1×106

SUMF：指定结果总结的输出频率(每隔Value时间输出一次结果总结)，缺省为1×106

BC：指定瞬态边界条件的变化形式，缺省为阶跃变化 其Value的含义如上面各项所述

第二条命令的 Lable选项及其各自含义如下： OVER：指定结果文件的覆盖频率(每隔Value个时间步覆盖一次)，缺省为0， 即不覆盖

APPE：指定结果的输出频率(每隔Value个时间步输出一次结果到结果文件中)，缺省为10

SUMF：指定结果总结的输出频率(每隔Value个时间步输出一次结果总结)，缺省为10

下页所示为该二命令的菜单形式，第一幅图为用户自定义时间步长，第二幅图为程序定义时间步长（上述对话框的第一项设定），同时，该二图所显示的都是设定由时间值和时间步同时控制求解及输出选项时（上述对话框的第二项设为“both”）的菜单情况：

ENKE：重新将整个模型的紊流动能设为一均值

ENDS：重新将整个模型的紊流动能耗散率设为一均值 DENS：重新将整个模型的密度设为一均值 VISC：重新将整个模型的粘性设为一均值

COND：重新将整个模型的导热系数设为一均值 EVIS：重新将整个模型的有效粘性设为一均值

ECON：重新将整个模型的有效导热系数设为一均值 TTOT：重新将整个模型的总温设为一均值 SPHT：重新将整个模型的比热设为一均值

另外， FLDATA29命令还有其它三个用于多组份疏运分析的Lable，其名称及含义如下：

SP0n：重新设定第n（n＝1～6）组份的质量份额

LMDn：重新设定第n（n＝1～6）组份的层流质量扩散系数 EMDn：重新设定第n（n＝1～6）组份的有效质量扩散系数

第一条命令（FLDATA28）的Value值为“是(T)”或“否(F)”之一，缺省为否。 第二条命令（FLDATA29）的Value值为各变量重新设定后的值，缺省为 0

注意：只有当在FLDATA28的相应菜单中将某变量的Value设为“是”后，才会

在FLDATA29命令的菜单中显示出该变量的Lable，以允许修改该变量的值。 该两条命令相应的菜单形式如下： 第一条命令（FLDATA28）：

第二条命令（FLDATA29）：

20、控制FLOTRAN面积积分的阶次(Quadrature Order)

命令：FLDATA30, QUAD, Lable, Value

功能：控制面积积分的阶次，亦即积分点的个数 菜单：

Main Menu>Solution>FLOTRAN Set Up>Mod Res/Quad Ord>CFD Quad Orders

其中，Lable的选项及其各自的含义如下：

MOMD：控制动量扩散项的面积积分阶次，缺省为0点积分

MMOS：控制动量源项的面积积分阶次，缺省为0点积分 PRSD：控制压力扩散项的面积积分阶次，缺省为1点积分 PRSS：控制压力源项的面积积分阶次，缺省为1点积分 THRD：控制热扩散项的面积积分阶次，缺省为0点积分 THRS：控制热源项的面积积分阶次，缺省为0点积分

TRBD：控制湍流扩散项的面积积分阶次，缺省为0点积分 TRBS：控制湍流源项的面积积分阶次，缺省为2点积分 Value即为上述积分点的数目，各值含义如下：

0：单点面积积分，但积分所用的是扩散系数的平均值 1：单点面积积分 2：两点面积积分

注意：程序总会将积分点数目设为最优值，例如，对于轴对称分析以及极坐标

下的分析，积分点数目会自动设为2。另，当单元畸形时，使用2点积分有助于提高精度。该命令的菜单形式如下：

21、FLOTRAN多组份疏运分析的设置及控制

下面所述所有命令的菜单路径都是：

Main Menu>Preprocessor>FLOTRAN Set Up>Multiple Species Main Menu>Solution>FLOTRAN Set Up>Multiple Species 命令：KEYOPT, 1, 1, Value

功能：设定组份数目，Value的值介于1到6之间

命令：MSDATA, ALGEB, UGAS

功能：设定施加在所有组份上的多组份数据

其中：ALGEB：设定一特定组份号，该组份的质量份额等于1.0减去其它所有组份质量份额之和，如此以保证总的质量份额为1.0。其缺省值为2 UGAS：设定大气常数，缺省为8314.3（国际单位制） 以上两条命令的菜单形式如下：

命令：MSSPEC, SPNUM, Name, MOLWT, SCHMIDT 功能：设定各组份的名字、分子量以及Schmidt数 其中：SPNUM：组份号，介于1到6之间

Name： 该组份的名字(1到4个字符, 不能与已存在的自由度的符号相同),缺省为SP0n，n为组份号

MOLWT： 该组份的分子量，只适用于气体(由MSPROP命令定义)

SCHMIDT：该组份的Schmidt数，只适用于气体(由MSPROP命令定义) 命令：MSQUAD, QDIF, QSRC

功能：定义多组份单元的面积积分阶次

其中：QDIF：扩散项面积积分的阶次，其值及含义如下： 0：单点积分（此为缺省值） 1：与0相似，只是在计算与温度相关的流体性质时，所使用的温度是分布温度（而不是平均值）

2：两点积分（作为轴对称分析时的缺省值）

QSRC：源项面积积分的阶次，其值及含义如下： 0：单点积分（此为缺省值） 1：与0相似，只是在计算与温度相关的流体性质时，所使用的温度是分布温度（而不是平均值）

2：两点积分（作为轴对称分析时的缺省值） 命令：MSNOMF, SPNUM, FRACTION 功能：定义各组份名义质量份额的初始值 其中：SPNUM：组份号，介于1到6之间

FRACTION：该组份在整个分析区域内的初始质量份额，缺省值为1/n。所有组份的该值之和应为1.0 上面三条命令的菜单形式如下：

命令：MSMETH, SPNUM, KEY

功能：指定多组份疏运方程的求解方法

其中：SPNUM：组份号，介于1到6之间（必须定义）

KEY：定义多组份疏运方程的求解方法，其值及其含义下： 0：不对该组份的方程进行求解

1：三对角矩阵法（TDMA法，此为缺省值） 2：共轭残差法（Conj residual）

3：预条件共轭残差法（Precond conj res）

4：预条件广义最小残差法（GMRES/precond） 该命令的菜单形式如下：

命令：MSSOLU, SPNUM, NSWEEP, MAXI, NSRCH, CONV, DELMAX 功能：定义所组份疏运分析的求解选项

5、FLOTRAN输出收敛监测量的控制

命令：FLDATA6, CONV, Lable, Value 功能：控制收敛监测量的输出

菜单：Main Menu>Solution>FLOTRAN Set Up>Additional Out>Print Controls 其中，Lable的选项及其各自含义如下： OUTP：控制收敛监测量的显示模式，可用Value值有LAND、BLOC或BNOW 详见Value之解释，缺省为BNOW模式 ITER：当OUTP为LAND或BLOC时，控制每隔 n 次迭代打印收敛监测参数，缺省 n 为1

Value的个值含义如下：

LAND：所有自由度都以水平模式打印输出(Lable=OUTP时有效) BLOC：所有自由度都以条状模式打印输出(Lable=OUTP时有效) BNOW：当一个迭代完成时，所有自由度都立即以条状模式打印输出(Lable=

OUTP时有效)

n ：迭代数(Lable=ITER时有效) 该命令的菜单形式如下：

6、FLOTRAN流体类型及其特性的可变性控制

命令：FLDATA7, PROT, Lable, Value 或 FLDATA12, PROT, Lable, Value 功能：定义流体类型

命令：FLDATA13, VARY, Lable, Value 功能：控制流体特性是否可变 以上二命令的菜单路径如下：

Main Menu>Preprocessor>FLOTRAN Set Up>Fluid Properties Main Menu>Solution>FLOTRAN Set Up>Fluid Properties 第一条命令的 Lable 的选项及其各自含义如下： DENS：指定流体密度类型，缺省为常值 VISC：指定流体粘性类型，缺省为常值

COND：指定流体导热系数类型，缺省为常值 SPHT：指定流体比热类型，缺省为常值

第一条命令的Value为流体性质的类型，有如下值： CONSTANT：性质为常值 GAS：气体性质 LIQUID：流体性质

TABLE：由材料性质表输入(MPTEMP和MPDATA命令) POWL：非牛顿流的 Power Law 粘性类型 CARR：非牛顿流的 Carreau 粘性类型 BING： 非牛顿流的 Bingham 粘性类型

USERV：用户自定义粘性类型(通过用户子程序USERVISLAW实现) AIR：国际单位制的空气性质

AIR_B：国际单位制的空气性质，其计算密度时的压力为参考压力 AIR-SI：国际单位制的空气性质

AIR-SI_B： 国际单位制的空气性质，其计算密度时的压力为参考压力 AIR-CM： 厘米-克-秒制的空气性质

AIR-CM_B：厘米-克-秒制的空气性质，其计算密度时的压力为参考压力 AIR-MM： 毫米-克-秒制的空气性质

AIR-MM_B：毫米-克-秒制的空气性质，其计算密度时的压力为参考压力 AIR-FT：英尺-斯-秒值的空气性质

AIR-FT_B：英尺-斯-秒值的空气性质，其计算密度时的压力为参考压力 AIR-IN：英寸-(磅-秒2/英寸)-秒制空气性质

AIR-IN_B：英寸-(磅-秒2/英寸)-秒制空气性质，其计算密度时的压力为

参考压力

CMIX： 多组份流体质量比 USER： 用户自定义流体性质 各性质的详细描述请参见理论手册

第二条命令的 Lable 的选项及其各自含义如下： DENS：指定流体密度是否变化，缺省为否 VISC：指定流体粘性是否变化，缺省为否

COND：指定流体导热系数是否变化，缺省为否 SPHT：指定流体比热是否变化，缺省为否

第二条命令的 Value 取“是”或“否”二者之一 这两条命令的菜单形式如下：

7、FLOTRAN流体性质参数设定

上面6种介绍的各种不同的流体类型，其特性的定义是不一样的，如，CONSTANT、 GAS、LIQUID等流体性质就分别需要定义其性质的名义项、第一参数项、第二参数项、第三参数项等，然后利用这些项通过各自特定的公式来拟合流体性质，定义这些项目的命令如下：

命令：FLDATA8, NOMI, Lable, Value 功能：定义流体性质方程的名义项参数

命令：FLDATA9, COF1, Lable, Value 功能：定义流体性质方程的第一参数项

命令：FLDATA10, COF2, Lable, Value 功能：定义流体性质方程的第二参数项

命令：FLDATA11, COF3, Lable, Value 功能：定义流体性质方程的第三参数项 这四条命令的菜单路径都是：

Main Menu>Preprocessor>FLOTRAN Set Up>Fluid Properties Main Menu>Solution>FLOTRAN Set Up>Fluid Properties

这四条命令的 Lable 的选项及其各自含义如下： DENS：指定密度特性的各参数项 VISC：指定粘性特性的各参数项

COND：指定导热系数特性的各参数项 SPHT：指定比热特性的各参数项

这些命令的菜单形式如下(该图只显示了LIQUID性质系数的定义，其它性质流体的性质系数的定义与此类似)：

本文来源：https://www.bwwdw.com/article/6mrh.html