PHOENICS 软件培训

PHOENICS软件培训前处理部分-Q1文件PHOENICS的基本结构Q1文件的组成TALK=T;RUN(1,1)GROUP1.RuntitleandotherpreliminariesGROUP2.Transience;time-stepspecificationGROUP3.X-directiongridspecificationGROUP4.Y-directiongridspecificationGROUP5.Z-directiongridspecificationGROUP6.Body-fittedcoordinatesorgriddistortionGROUP7.Variablesstored,solved&namedGROUP8.Terms(indifferentialequations)&devicesGROUP9.Propertiesofthemedium(ormedia)GROUP10.Inter-phase-transferprocessesandpropertiesGROUP11.InitializationofvariableorporosityfieldsGROUP12.ConvectionanddiffusionadjustmentsGROUP13.BoundaryconditionsandspecialsourcesGROUP14.DownstreampressureforPARAB=.TRUE.GROUP15.TerminationofsweepsGROUP16.TerminationofiterationsGROUP17.Under-relaxationdevicesGROUP18.LimitsonvariablesorincrementstothemGROUP19.DatacommunicatedbysatellitetoGROUNDGROUP20.Preliminaryprint-outGROUP21.Print-outofvariablesGROUP22.Spot-valueprint-outGROUP23.Fieldprint-outandplotcontrolGROUP24.DumpsforrestartsSTOPGROUP1Runtitleandotherpreliminaries定义作业的名称，定义变量（标准变量35个，主要是指PIL中没有的变量）变量类型可以有：数组（array）逻辑变量（boolean）字符型变量（char）实型（real）整型（integer）GROUP2定义非定常计算相关的参数Steady=T（定常，与时间无关）F（非定常，与时间有关）TFIRST开始时间TLAST结束时间TSTEP时间步长GROUP3，4，5X，Y，Z方向网格划分的相关参数坐标系的选取Cartes=T（直角坐标）F（柱坐标系统，X角度，Y半径）NX，NY，NZ表示网格的数目Xulast，Yulast，Zulast表示X，Y，Z方向相应的尺度Nregx，Nregy，Nregz表示X，Y，Z方向区域的个数Para=T表示求解的为抛物型方程GROUP6贴体坐标以及相关的网格参数的设定BFC=T（采用贴体坐标系统）Domain规定计算区域的网格数Grdchk=T网格的正交性的检查Readco读入外部的网格文件注意：网格文件的格式必须满足PHOENICS网格文件的格式，才能正确的读入。GROUP7变量定义定义需要求解的变量，和存储的变量（包括因变量和各种辅助变量）。激活模型（如激活K-epsilon湍流模型：Turbmod（kemod））Store存储变量Solutn定义求解变量Onephs=T表示单相流动F表示两相流动GROUP8选择方程中需要保留的项其中包括：瞬变项（时间项），对流项，扩散项，源相这些项的选择可以根据方程求解的需要来定。Terms（变量，Y/N，Y/N，Y/N，Y/N，Y/N，Y/N）表示1：激活内置的源相2：激活对流相3：激活扩散相4：激活瞬变相5：变量属于第一相流体6：激活相间的输运GROUP9定义介质的物理属性密度Rho1粘性Enul（层流），Enut（湍流）热导系数比热激活湍流模型（Turmod）注意：标量方程的交换系数是通过Prndtl数来规定其介质的输运系数的，Prndtl为负数的时候，表示扩散系数。GROUP10定义多相流的能量交换，动量交换等特性Cfips：相间的输运系数（Grnd1…Grnd10）Cmdot：相间的质量输运系数GROUP11给定变量以及多孔度的初场，自启动功能Fiinit（phi）给定变量的初场初场也可以通过patch和coval命令来给定Restrt：变量的自启动功能对于强的非线性问题，自启动功能有利于迭代的收敛GROUP12用于调节对流和扩散的特性主要用Patch命令来实现，不同的Patchname实现不同的功能：如：Gp12con？调节所有的对流项Gp12sor？调节所有的内置源相Gp12cne？调节网格东面的对流项......GROUP13边界条件与特殊源相的设置PHOENICS中边界条件的处理和其他的一些商用程序不同，它有自己的一套处理边界条件的方法。主要分两步来实现：1：什么地方是边界，什么时候？2：什么变量规定什么样的边界？PIL命令主要是：Patch（什么地方，什么时候）Coval（什么变量什么样的边界）PHOENICS中的边界条件PHOENICS中的边界条件主要是以源相的形式给出的对流扩散源相面积因子系数值经过线形化后的源相离散化后方程可以写成如下形式关于边界条件的PIL命令边界的位置PATCH（Name，Type，Ixf，Ixl，Iyf，Iyl，Izf，Izl，Itf，Itl）Name：规定边界的名称Type：边界的类型给出变量的边界Coval（Name，Variable，Coefficient，Value）Coefficient：系数Value：值具体的几种常用的边界条件的给定固定值PHOENICS中可以用FIXVAL来表示，给定边界的固值。比如：COVAL（Patchname，Phi，Fixval，Value）给定通量PHOENICS中可以用FIXVAL来表示，给定边界的固值。比如：COVAL（Patchname，Phi，Fixval，Value）GROUP14对于抛物型问题（PARA=T）给定下游的压力条件GROUP15外层迭代（Sweep）的次数DOISTEP=1,LSTEP!Transienceloop|VisitSec1:startoftime-step|DOISWEEP=1,LSWEEP!Outer-iterationloop||VisitSec2:startofsweep||DOIZ=1,NZ!Single-sweeploop|||VisitSec3:startofIZslab|||DOITHYD=1,LITHYD!Hydrodynamicloop||||VisitSec4:startofiteration||||Solveequations||||VisitSec5:finishofiteration|||ENDDO|||VisitSec6:finishofIZslab||ENDDO||VisitSec7:finishofsweep|ENDDO|VisitofSec8:finishoftimestepENDDOGROUP16规定内层迭代（Iteration）的次数即用迭代法求解线形方程组的迭代次数GROUP17规定松弛因子有两种方法给出松弛因子线形松弛因子（Linrlx）虚拟时间步松弛因子（Falsdt）湍流源相的线型化方法（2种）GROUP18对迭代变量的幅度进行限制Varmax给出变量幅度的最大值Varmin给出变量幅度的最小值GROUP19用于Q1和Ground之间的数据传递Usegrd=TEarth访问Ground.forUsegrx=TEarth访问Grex3.forGROUP20-24控制数据的输出，图形的输出，以及程序的调试功能输出到Result文件数据的计算结果可以通过Output命令来控制。程序的调试功能。PIL命令的编程功能PIL命令语言就象是一门高级编程语言，如果你熟练掌握的化，可以编非常精彩的前处理程序。在PHOENICS提供的一些列子中，可以看出PIL命令非常精彩的一面。PIL中可以实现循环控制，数据的交互输入，条件控制，以及子程序等。THEEND