一、 在PHOENICS中导入gridgen生成的贴体网格
虽然gridgen号称能够生成PHOENICS的网格文件,但实际上是有一定问题的。两者的网格文件格式存在一些差别。
PHOENICS的网格文件格式:
NX+1 NY+1 NZ+1
XC(1,1,1) XC(1,2,1) XC(1,3,1) XC(1,4,1) XC(1,5,1)
XC(1,6,1) .....
............. XC(1,NY+1,1) XC(2,1,1) XC(2,2,1)
............. XC(2,NY+1,1) XC(3,1,1) ........
.............XC(NX+1,NY+1,1)
YC(1,1,1) YC(1,2,1) ............................
.............YC(NX+1,NY+1,1)
ZC(1,1,1) ZC(1,2,1) ............................
.............ZC(NX+1,NY+1,1)
XC(1,1,2) XC(1,2,2).............................
..............ZC(NX+1,NY+1,2)
XC(1,1,3).......................................
..............ZC(NX+1,NY+1,NZ+1)
Gridgen生成的PHOENICS网格文件格式:
NX+1 NY+1 NZ+1
XC(1,1,1) XC(2,1,1) XC(3,1,1) XC(4,1,1) XC(5,1,1)
XC(6,1,1) ...XC(NX+1,1,1) XC(1,2,1) XC(2,2,1)
……XC(NX+1,2,1)……XC(NX+1,NY+1,1)………
XC(1,1,2)…………XC(NX+1,NY+1,NZ+1)
YC(1,1,1) YC(2,1,1) YC(3,1,1) YC(4,1,1) YC(5,1,1)
YC(6,1,1) ...YC(NX+1,1,1) YC(1,2,1) YC(2,2,1)
……YC(NX+1,2,1)……YC(NX+1,NY+1,1)………
YC(1,1,2)…………YC(NX+1,NY+1,NZ+1)
ZC(1,1,1) ZC(2,1,1) ZC(3,1,1) ZC(4,1,1) ZC(5,1,1)
ZC(6,1,1) ...ZC(NX+1,1,1) ZC(1,2,1) ZC(2,2,1)
……ZC(NX+1,2,1)……ZC(NX+1,NY+1,1)………
ZC(1,1,2)…………ZC(NX+1,NY+1,NZ+1)
以上文件中,NX,NY,NZ分别为X,Y,Z三个方向的网格数,XC,YC,ZC分别为X,Y,Z三个方向的网格节点坐标。
解决办法:可以自己编写一个的程序,读入gridgen生成的网格文件,然后进行转换,生成适合PHOENICS所需格式的网格文件。
例子:如何将gridgen生成的贴体网格导入PHOENICS
第一步:gridgen生成了一个网格文件:BLK1。现在用自己编写的转换程序进行转换,生成新的网格文件BLK1-N。
第二步:进入PHOENICS,将坐标改为贴体坐标系。点击 FileOpen file for editingQ1,进入Q1文件,在其中的GROUP6中添加语句:
BFC=T
READCO(BLK1-N)
即可。生成的贴体坐标如下图:
相应的Q1文件为:Q1.Q1 |