声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3280|回复: 4

[基础理论] 前处理部分-Q1文件中文讲解

[复制链接]
发表于 2005-9-26 10:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
前处理部分-Q1文件

Q1文件的组成

GROUP 1
Run title and other preliminaries
定义 作业的名称,定义变量
(标准变量35个,主要是指PIL中没有的变量)
变量类型可以有:数组(array)
逻辑变量(boolean)
字符型变量 (char)
实型 (real)
整型 (integer)

GROUP 2
定义非定常计算相关的参数
Steady=T(定常,与时间无关)
F (非定常,与时间有关)
TFIRST 开始时间
TLAST 结束时间
TSTEP 时间步长

GROUP 3,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 表示求解的为抛物型方程

GROUP 6
贴体坐标以及相关的网格参数的设定
BFC = T (采用贴体坐标系统)
Domain 规定计算区域的网格数
Grdchk = T 网格的正交性的检查
Readco 读入外部的网格文件
注意:网格文件的格式必须满足PHOENICS
网格文件的格式,才能正确的读入。

GROUP 7
变量定义
定义需要求解的变量,和存储的变量(包括因变量和
各种辅助变量)。
激活模型(如激活K-epsilon湍流模型: Turbmod(kemod))
Store 存储变量
Solutn 定义求解变量
Onephs = T 表示单相流动
F 表示两相流动

定义需要求解的变量,和存储的变量(包括因变量和
各种辅助变量)。
激活模型(如激活K-epsilon湍流模型: Turbmod(kemod))
Store 存储变量
Solutn 定义求解变量
Onephs = T 表示单相流动
F 表示两相流动

GROUP 8
选择方程中需要保留的项
其中包括:瞬变项(时间项),对流项,扩散项,源相
这些项的选择可以根据方程求解的需要来定。
Terms(变量,Y/N,Y/N,Y/N,Y/N,Y/N,Y/N)
表示 1:激活内置的源相
2:激活对流相
3:激活扩散相
4:激活瞬变相
5: 变量属于第一相流体
6: 激活相间的输运


GROUP 9
定义介质的物理属性
密度 Rho1
粘性 Enul(层流),Enut (湍流)
热导系数
比热
激活湍流模型(Turmod)
注意:标量方程的交换系数是通过Prndtl数来规定其介质的输运系数的,Prndtl为负数的时候,表示扩散系数。

GROUP 10
定义多相流的能量交换,动量交换等特性
Cfips: 相间的输运系数(Grnd1…Grnd10)
Cmdot:相间的质量输运系数
GROUP 11
给定变量以及多孔度的初场,自启动功能
Fiinit(phi)给定变量的初场
初场也可以通过patch和coval命令来给定
Restrt:变量的自启动功能
对于强的非线性问题,自启动功能有利于迭代
的收敛

GROUP 12
用于调节对流和扩散的特性
主要用Patch命令来实现,不同的Patchname实现不同
的功能:
如:Gp12con? 调节所有的对流项
Gp12sor? 调节所有的内置源相
Gp12cne? 调节网格东面的对流项
. .
. .
. .
GROUP 13
边界条件与特殊源相的设置
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: 值
具体的几种常用的边界条件的给定
固定值
回复
分享到:

使用道具 举报

 楼主| 发表于 2005-9-26 10:12 | 显示全部楼层

回复:(hyacinth)前处理部分-Q1文件中文讲解

PHOENICS中可以用FIXVAL来表示,给定边界的固值。
比如:COVAL(Patchname,Phi,Fixval,value)
给定通量 GROUP 14
对于抛物型问题(PARA=T)给定下游的压力条件

GROUP 15
外层迭代(Sweep)的次数

GROUP 16
规定内层迭代(Iteration)的次数
即用迭代法求解线形方程组的迭代次数

GROUP 17
规定松弛因子
有两种方法给出松弛因子
线形松弛因子(Linrlx)
虚拟时间步松弛因子(Falsdt)
湍流源相的线型化方法(2种)

GROUP 18
对迭代变量的幅度进行限制
Varmax 给出变量幅度的最大值
Varmin 给出变量幅度的最小值

GROUP 19
用于Q1和Ground之间的数据传递
Usegrd=T Earth 访问Ground.for
Usegrx=T Earth 访问Grex3.for

GROUP 20-24
控制数据的输出,图形的输出,以及程序的调试功能
输出到Result 文件数据的计算结果可以通过Output 命令来控制。
程序的调试功能。
PIL命令的编程功能
PIL命令语言就象是一门高级编程语言,如果你熟练掌握的化,可以编非常精彩的前处理程序。在PHOENICS提供的一些列子中,可以看出PIL命令非常精彩的一面。
PIL中可以实现循环控制,数据的交互输入,条件控制,以及子程序等。
THE END
发表于 2006-3-16 09:10 | 显示全部楼层
谢谢了!!!!!!!
发表于 2013-10-14 14:07 | 显示全部楼层
发表于 2014-10-16 12:15 | 显示全部楼层
学习了学习了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-29 23:04 , Processed in 0.075514 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表