PHOENICS的原理与应用
- 格式:pdf
- 大小:608.48 KB
- 文档页数:8
第1卷 第4期 数值计算与工程仿真 Vol.1 No.4 2005年 3月 Numerical Calculation and Engineering Simulation March,2005
——
清洁能源技术论坛版权所有,www.myCFD.cn. ─ 73 ─计算流体力学软件PHOENICS及其工程应用
中国石油大学机电学院 马培勇,仇性启,唐鹏
摘 要: 计算流体力学是用计算机来分析流动的尖端技术,用计算模型描述所建立和研
究的物理系统。流体流动物理学运用于有效模型时,CFD软件可模拟出流体流动、传热、
传质、相变、化学反应、机械运动以及固体结构的压力变形。PHOENICS是流体力学中通
用性很强的商业软件之一,能够有效对上述现象进行数量上的预测。本文对该软件的基本
结构、物理模型及算法、边界条件及源项的处理方法以及其功能特点与应用作简要介绍。
Application of CFD Code PHOENICS in Engineering
MA Pei-yong, QIU Xing-qi, TANG Peng
(College of Machinery and Engineering, University of Petroleum, Dongying
257061,China)
Abstract: Computational fluid dynamics (CFD) is a sophisticated analysis technique that uses
computers to simulate fluid flow and enables a computational model representing a physical
system to be built and studied. When fluid flow physics is applied to this virtual prototype, the
CFD application outputs a prediction of the fluid dynamics. The software predicts not only the
fluid flow behavior, but also the transfer of heat, mass, phase change, chemical reaction,
mechanical movement, and stress or deformation of related solid structures. PHOENICS is one
of the general-purpose software packages and it can predict the above phenomena quantitatively.
The basic structure , physics model, arithmetic, management of boundary condition and source ,
abilities and application of PHOENICS is introduced simply in the paper.
1前言
计算流体力学(Computational Fluid Dynamics,简称CFD,以下同)是流
体力学、偏微分方程数值理论、数值分析和计算机科学等多种领域的交叉学
科,主要用于解决工程中的流体和传热问题。随着计算机技术的推广普及和计第1卷 第4期 数值计算与工程仿真 Vol.1 No.4 2005年 3月 Numerical Calculation and Engineering Simulation March,2005
——
清洁能源技术论坛版权所有,www.myCFD.cn. ─ 74 ─算方法的发展,取得了蓬勃的发展。由于数值模拟相对实验研究有成本低、周
期短的优点,可模拟出实际运行过程中各种状态的数据,对于设计、改造等技
术方案有重要的指导作用,因此CFD技术得到了越来越多的应用。目前应用
比较好的CFD软件有Phoenics、Fluent、CFX、Star-CD等。
PHOENICS是英国CHAM公司开发的模拟流动、传热、反应、燃烧过程
的通用CFD软件,有30多年的历史。PHOENICS是Parabolic Hyperbolic Or
Elliptic Numerical Integration Code Series起始字母的缩写词。该软件的基本功
能可计算质量、动量和能量守恒的有限体积方程,可以对三维稳态或非稳态的
可压缩流或不可压缩流进行模拟,包括非牛顿流、多孔介质中的流动,并且可
以考虑粘度、密度、温度变化的影响。对于特定方向的流动现象有相应的“抛
物线”计算模块,并且还有一些附件,如用于电子冷却的HOTBOX模块,用
于导热和通风的FLAIR,以及用于化学汽相淀积(Chemical Vapor Deposition)
模拟的CVD。PHOENICS内置22种适合于各种Re数场合的湍流模型,包括
雷诺应力模型、多流体湍流模型和通量模型及k-ε模型的各种变异,共计21
个湍流模型,8个多相流模型,10多个差分格式。网格系统包括直角、圆柱、
曲面、多重网格、精密网格。可以直接读入CAD软件建立的模型,使复杂几
何体的生成更为方便,在边界条件的定义方面也极为简单,并且网格自动生
成,另外,自带了1000多个例题与验证题,附有完整的可读可改的输入文
件。开放性好,提供对软件现有模型进行修改、增加新模型的功能和接口,可
以用FORTRAN语言进行二次开发[1]。
2基本结构
PHOENICS程序以“行星”状布置其子程序(图1),主要由三个工作模
块组成:前处理器Pre-processor(问题的定义与表述)、求解器Solver(数值模
拟计算)和后处理器Post-processor(计算结果的显示)。
Q1为PHOENICS的输入文件,用于建立数值模拟模型,如激活数值模拟
中所要用到的控制方程,输入材料的物理性能参数,设置迭代次数,控制方程
的求解方法,以及实现网格的划分等。这些任务均由Q1的24组语句在相应的
位置用相应的命令中实现。
SATELLITE程序的功能是接受数据,包括三个FORTRAN子程序:
MAIN、SAT和SATLIT,主程序MAIN将Q1文件翻译成FORTRAN程序可第1卷 第4期 数值计算与工程仿真 Vol.1 No.4 2005年 3月 Numerical Calculation and Engineering Simulation March,2005
——
清洁能源技术论坛版权所有,www.myCFD.cn. ─ 75 ─读取的数据文件,用户可以在SATLIT子程序中加入自己的FORTRAN程序,
以补充Q1文件无法表示的数据项;SATLIT子程序的格式与Q1文件相同也分为24组;SAT负责调用用户子程序。
图1 PHOENICS程序结构图
EARTH程序是模拟的核心部分,也包括三个FORTRAN子程序:
MAIN、GROSTA和GROUND,主程序MAIN负责读取SATELLITE生产的数
据文件,并把方程离散化和解方程组,用于用户扩展求解器EARTH功能的工
具GROUND子程序具有与Q1相似的结构,在PHOENICS程序本身不能满足
用户需要时,可以在GROUND中相应的位置加入用户自己的程序,从而引进
复杂的非线性边界条件,提供物性公式等。图中所表示的工作过程是:在Q1
文件中写入数据来建立模型;SATELLITE读取Q1并且编写中介数据文件
EARDAT;求解器EARTH读取EARDAT之后,对问题进行求解,并且写出
PHI、RESULT和PATGEO文件(其中PHI包括所有输入输出数据;对问题进
行求解,并且写出PHI、RESULT和PATGEO文件(其中PHI包括所有输入输
出数据;RESULT则包括可读的文本输入和计算结果;而PATGEO则用于
PHOTON进行绘图);计算完成后,进行后处理工作,可以采用VR(Virtual-
Reality,虚拟现实)视图器显示各种计算变量的等值线图、矢量图或等值面
图,或者应用PHOTON绘制矢量图,而AUTOPLOT也是一种绘图工具,最第1卷 第4期 数值计算与工程仿真 Vol.1 No.4 2005年 3月 Numerical Calculation and Engineering Simulation March,2005
——
清洁能源技术论坛版权所有,www.myCFD.cn. ─ 76 ─后的RESULT文件则用来存储计算结果的数据等。
3 物理模型及算法
一个软件计算问题的准确性,主要取决于他所建立的模型的准确性、真实
性。对于PHOENICS软件的控制方程是由连续性方程、动量方程、能量方程
和组分方程所组成,可统一写成[2]:
ΦΦ=Φ∇Γ−Φ+Φ∂∂Sudivt)()(ρρ (1)
式中,Φ表示所求解的变量,例如焓、温度、速度、电势和质量等;ρ表示密
度;u表示对应方向的速度;ΓΦ表示对流项;SΦ表示源项;t表示时间。
PHOENICS控制方程组具有封闭性、形式一致性、耦合性和非线性的特
点。该软件采用有限容积法来实现控制方程的离散化,可以选择一阶迎风、混
和格式及QUICK格式等;流速和压力耦合采用用SIMPLEST算法[2];对两相流
纳入IPSA算法(适用于两相介质互相穿透时)、PSI-Cell(粒子跟踪法)
GENTRA算法(颗粒轨迹);代数方程可以采用整场求解或点迭代、块迭代方
法,同时纳入块修正加快收敛;用交错网格解决一阶导数项离散化降阶问题,
用标准k-ε湍流模型封闭时均方程;对流扩散项离散化用混合格式,界面粘
度系数与密度取相邻结点的算术平均值,压力场全域求解,速度场逐行
TDMA。
SIMPLEST算法与SIMPLE算法基本相同,不同点主要在求解动量方程的
方法上。SIMPLE算法在求解动量方程时收敛性差,主要原因是动量方程离散
化后的对流项有问题。因此SIMPLEST对SIMPLE作如下改进:1)采用迎风
格式。因为迎风格式是一个绝对稳定的格式,且扩散项与对流项的影响系数可
以分离。至于由迎风差分所引起的假扩散问题,则采取逐步加密网格以获得与
网格稀密程度无关的解这种做法加以克服。2)动量方程的系数只包含扩散项
的作用,而将对流项的作用都放入源项中;采用Jacobi逐点求解动量方程,免
去原来求解速度场时所需要的很强的低松弛,促进强烈非线性问题的迭代过程
收敛,SIMPLEST的计算步骤则与SIMPLE基本相同,算法流程图如图2。