文档介绍：Xi’an Jiaotong University
数值传热学
第11章求解二维椭圆型流动与换热问题通用程序介绍
(续)
主讲陶文铨
西安交通大学能源与动力工程学院
热流科学与工程教育部重点实验室
2012年12月5日, 西安 1/63
Xi’an Jiaotong University
11. 5 程序使用方法及主程序说明
使用方法
主程序说明
MAIN PROGRAM
SUBROUTINE DIFLOW
SUBROUTINE SOLVE
SUBROUTINE SETUP
SUBROUTINE SUPPLY
2/63
Xi’an Jiaotong University
11. 5 程序使用方法及主程序说明
使用方法
形式相比较:
()
div() u div ( grad ) S
t
确定 S, 及.
,mon块,修正
其余部分,并以新的USER名称存储.
3/63
Xi’an Jiaotong University
(每一方向5~7个节点),较
少迭代 LAST=3-5次,先通过语法调试. 通过后逐
步增加问题的复杂性,例如要计算紊流的先计算层流,
计算高Re 数的先计算低Re 。
,按每个ENTRY的功能及具
体问题对有关变量及参数赋值。特别注意以下几个变量
的赋值:
(1) LSOLVE(NF)—.
(2) LPRINT(NF)—.
4/63
Xi’an Jiaotong University
(3) TITLE(NF)—要打印的NF变量给出打印标题
(限七个字符)。
(4) LBLK(NF)—对要采用块修正的NF变量赋
.TRUE.,.,. 。
(5) LAST—最大迭代次数,预置值为5。
(6) NTIMES(NF)—预置值为1,稳态非线性问题取
12, 非稳态线性问题取 5~6。
(7) DT—非稳态问题选定时间步长,预置值为 1010 。
0 V
全隐格式的特点:离散方程中多了项aP ,
0 t
而,当ta,0P ,导致稳态计算。
5/63
Xi’an Jiaotong University
(8) RELAX(NF)—
弛因子之值。
(9) 确定参考压力点的ij, 值, IPREF及JPREF的预
置值为1,1。
(10)如果定义了新的变量,如T(NI,NJ),则应用
EQUIVALENCE 与 F(NI,NJ,NFMAX) 对应起来,
如EQUIVALENCE (F(1,1,4),T(1,1)).
6/63
Xi’an Jiaotong University
11-5-2 主程序说明
C puter program was copied from the graduate student course
C program of the University of Minnesota. Part of it was re-formulated
C to meet the puter environment. Some inappropriate
C expressions were also corrected. The program is used only for the
C teaching purpose. No part of it may be published. You may use it
C as a frame to re-develop your own code for research purpose.
C --------Instructor of Numerical Heat Transfer, XJTU,-----------------
C***************************************************************
7/63
Xi’an Jiaotong University
C***********************************************************
---------------------