phoenics命令.docxVIP

ID=:YAMOS(JBY;71G1AS7.SNGRSFPHOENICS 输入库中第1部分件case91，92，93为在线教学。PRESSO为参考压力（通常为1atm）P1 为相对压力绝对压力=PRESSO+P1RHO1为第一相密度RHO1A 为燃料分子量RHO1B为氧化剂分子量RHO1C为产物分子量RHO2为第二相密度U1为第一相的X方向速度V1为第一相的Y方向速度W1为第一相的Z方向速度R1为第一相的体积分数H1为第一相的焓TEMP0为参考温度TMP1为第一相的温度TMP2为第二相的温度EL1为第一相的混合长度尺度EL2为第二相的混合长度尺度PRNDTL（PHI）为该变量的层流普郎特数PRT（PHI）为该变量的湍流普郎特数ENUL为层流粘性系数ENUT为湍流粘性系数KE 为湍流动能EP 为湍流耗散率EL1为长度尺度PRNDTL 普郎特数TEMP1 第一相温度TEMP2 第二相温度若计算网格为柱坐标情况下，周向定义为I（X），径向定义为J（Y），轴向定义为K（Z）；网格数为I ×J ×K。计算区域特殊情况处理，对燃气喷孔，按其面积相当的当量直径处理成正方形，即凡是流体所到达不了的区域，均设空度为零，这样有助于网格生成与编号处理湍流模型采用模型，激活湍流模型的命令为TURMOD（KEMODL）燃烧模型：三气体燃烧模型，采用综合辐射模型，应符合：简单化学反应假设；快速反映假设；综合辐射模型（六通量模型）对入口边界，首先有空气入口（助燃风）和燃料入口（水煤气）两部分，这时速度的各个分量取第一类边界条件，即：注：对空气（助燃风）而言，不同的分速度显然反映的是旋流数不同，由此可研究旋流数的不同所造成的影响对入口，取，其中为来流的平均速度对入口，取，其中为入口当量直径X方向循环的定义：流体通过X=0的边界与通过X=HIGH的边界相同。循环的使用情况为：带有涡旋的轴对称流动以及三维有涡旋的柱坐标情况EARTH 的计算结果为两个文件：RESULT：可以直接读；PHIDA：要通过PHOTON或AUTOPLOT处理，转换成可以输出的文件PHOENICS的六通量模型只使用于3维笛卡尔坐标系抛物问题不适宜BFC坐标（有待验证） PATCH 和COVAL的应用不适于压力P的边界，主要有三种边界：固定流边界；固定值边界；壁面边界固定流边界：设边界的固定质量流为，设C为很小值C=1E-10，V=Q/C，由于，则方程中的原项为：在PHOENICS中，设置系数为FIXFLU（默认很小值1E-10）28．固定值边界：C给大值（1E10），此时，，由于C很大，所以，即。在PHOENICS中，设置系数为FIXVAL（默认很小值1E10）壁面边界：（边界层）通过剪切力（粘性行为）输运着动量壁面到网格点的距离为d，原项为，线性化为系数由此壁面边界的PATCH类型为：NWALL，SWALL等初始值的输入有几种方法：FIINIT的应用或PATCH和INIT的组合应用。PIL语言最少应从第三列开始注意高低：X方向，低为西，高为东；在Y方向，底为南，高为北；网格生成的三种方法：让SATELLITE读一个在PHOENICS运行前已编好的网格生成程序（应采用此方法）在与SATELLITE交互方式中，逐点输入网格坐标值在与SATELLITE交互方式中，用户确定角坐标，其余网格坐标由给顶的代数方程计算得到注释行必须在第一列或第二列PIL命令中，键入STORED可以看到存储的变量，键入SOLVED可以看到要求解的变量；键入SEELIB 252可以看到252库文件的内容。假定网格角坐标数据被放在一个文件GRID名中，用命令读此文件进入SATELLITE，这种网格可由VDU来显示（重要）42. 程序的第一行应输入：TALK=T；RUN（1，1）；VDU=VGAMOUSE（交互模式）TALK=F；RUN（1，1）；VDU=VGAMOUSE（编写Q1文件）satellite中的帮助：commands?可获得PIL需要的命令EDITOR？可获得怎样运用内置编辑器ABORT结束不存储文件END结束存储文件GROUP学习PHOENICS数据输入的组织GROUP N学习本组的命令和变量LOAD？学习卸载库文件要卸载核心库，用命令LOAD（CASE-NUMBER）要卸载对话库，用命令LOAD（LIBKEYCASE-NUMBER）LIBKEY是单字母要卸载用户自己库文件，用命令：LOAD（UCASE-NUMBER）(8)SEELIB244可以看到库文件244的内容；也可先用命令LOAD（244）把文件卸载下来，然后用SEE命令浏览，或用SEE N浏览某一组的具体内容。PHOENICS命令应在第一列或第二列开始，从第三列开始的均属于注释。72．命令XULAST指X方向的总长度，单位是米，在笛卡坐标下指X方