[转帖]PHOENICS Q&A

frogfish

来自各个论坛的讨论

存储空间不够的问题


Q；我的程序3.4网格设置加密以后，运行出现了这个错误
大家帮忙亚
result文件如下
************************************************************


---------------------------------------------------------
CCCC HHH PHOENICS Version 3.4.0 - EARTH
CCCCCCCC H (C) Copyright 2001
CCCCCCC See H Concentration Heat and Momentum Ltd
CCCCCCC our new H All rights reserved.
CCCCCC Web-site H Address: Bakery House, 40 High St
CCCCCCC www.cham. H Wimbledon, London, SW19 5AU
CCCCCCC co.uk H Tel: 0208-947-7651
CCCCCCCC H Fax : 0208-879-3497
CCCC HHH E-mail: phoenics@cham.co.uk
---------------------------------------------------------
This program forms part of the PHOENICS installation for:
CHAM
The code expiry date is the end of : Sep 2011
---------------------------------------------------------

************************************************************

prndtl of temperature was 1.0, so EARTH
re-set it to -grnd10.
OK? If not, please choose another value
material properties used are...
denst1
vistrb
visclm
tempr1
mixln1
thrme1
speht1
also, other related settings are ...
usegrx = T
usegrd = T
Number of F-array locations available is 30000000
Number used before BFC allowance is 33317715

* 3D velocity coeffs stored out of core
Number of F-array locations available is 30000000
Number used before BFC allowance is 30611667

* Dependent variables stored out of core
Number of F-array locations available is 30000000
Number used before BFC allowance is 22719027
Number used after BFC allowance is 22719027
F-array too small. Please re-dimension
to at least the size shown below
nftot =34521435 nfmax =30000000 needed =45158739
A：你计算所需要的储存空间 =34521435 而phoenics的最大空间 nfmax =30000000 需要册除needed =45158739 或增大此空间
Q；如何增大？
如何删除？
A；册除就是把计算网格\计算变量减少\储存变量减少\
增大就是把phoenics的数组增大

A；main.for里面的nfdim数组增大,然后编译连接,run private eart
Q；按照您说的做了之后出现这些问题

C:\phoenics\d_priv1>call \phoenics\d_satell\..\d_utils\d_windf\..\phoepath.bat
Adding DIGITAL PHOENICS to path
系统找不到指定的路径。
'df' 不是内部或外部命令，也不是可运行的程序
或批处理文件。
'df' 不是内部或外部命令，也不是可运行的程序
或批处理文件。
'link' 不是内部或外部命令，也不是可运行的程序
或批处理文件。
找不到 C:\phoenics\d_priv1\temp.exe
找不到 C:\phoenics\d_priv1\main.obj
找不到 C:\phoenics\d_priv1\ground.obj
请按任意键继续 . . .


怎么办呀？ df是什么意思？
A：DF means the Compaq FORTRAN Compiler.
If you have not installed the compiler,you can not use ground
to recompile the earexe.
Q：我安装了Compaq FORTRAN Compiler
可是里面没有编辑器呀

不知道是不是得另外安装编辑器？

A：I think you can get the answer from the POLIS with the title,
Documentation /Hard-copy documentation
/Installation of PHOENICS-3.5(or.3.4); TR 110
/Appendix B. Frequently Asked Questions and Troubleshooting

frogfish

Q:已设置入口条件。
出口条件为pressure 0Pa;
Coefficient 1000;（是什么意思？）
Velocity X 设为in-cell;
Velocity Y 设为in-cell;
Velocity Z 设为in-cell;
是不是出口的速度由软件自动计算给出？自己设置怎么办？
A: 我的理解是:在phoenics中出口是使用相对压力作为边界条件的,在这里的系数是压力源相的系数,即是SOURCE=T*COEFFICIENT*(value-P0),系数通常是10-1000
我的理解是在PHOENICS中很少设置在出口处的速度的, 不知道这样的理解对不对??
可以讨论一下哦
Q:请问T是什么？P0是不是大气压？SOURCE是什么？
A: T：源项类型是几何因子，由patch的type定,比如入口条件，你的选择的是east那么T就是单元EAST面的面积。
在需要固定一个值的时候，需要一个很大的系数FIXVAL=1E10
但是当需要把压强固定为0的时候，从数值方面考虑建议使用小于1E10
这里用1000是满足需要的
如果你设置体积力的时候TYPE 你设置的VOLUM 那么T就是单元体积。
p0:标准大气压 1e5
source应该是源项

Q:出口设置Coefficient 1000;（是什么意思？）那么这个系数是反应什么的呢？？
Velocity X 设为in-cell;与Velocity X=0
Velocity Y 设为in-cell;与Velocity Y=0
Velocity Z 设为in-cell;与Velocity Z=0
对计算结果有影响吗？一般应该如何设置？
A:我也不太明白，但凭我的经验是：出口也可以当入口，入口也可以当出口，只是一旦选择了类型为INLET，就作为入口处理，实际上给入口的速度可以为负值，就成为实际中的出口了。OUTLET也有这个问题。这些问题不要去管它，但只遵循一个原则就行了。这个原则是：INLET是固定速度的，OUTLET是固定压力的。所以，设入口还是出口完全是根据你要给的初始边界条件来选择的。
对于出口的系数，单相计算默认就行了，没有什么影响。但对于两相计算，就会有两个相的系数，这时就不能随便设了，而是两个相的系数与密度是成比例的才行，这一点我刚刚搞懂，现在正在算着呢

编者注：This command is used to declare a fixed-pressure boundary, where mass is allowed to leave (or enter) the solution domain

frogfish

Q：在温度设置中有个static/total的选项，请问有什么区别？
A：好像很多人都在问这个问题，
你可以参见TR326，TR326里讲到了VR里面所有按钮的具体含义。比如：

Models - Energy Equation, Total/Static
By default, the Temperature form is set to 'Total', the enthalpy form to 'Static'. The static form includes the substantial derivative of the pressure and the kinetic heating terms in the energy conservation equation, the Total form does not.

If the flow is highly compressible (high Mach number) the Temperature form should be switched to 'Static' otherwise incorrect solutions will be obtained. The Enthalpy form can be used in 'Total' form as long as a suitable temperature derivation is selected in the properties panel.
A：我的理解是一般的不可压缩流体的温度形式的能量方程就是用TOTAL，高马赫数的可压缩流体的温度形式的能量方程必须采用STATIC
焓形式的能量方程一般采用STATIC形式，因为其中包含压力以及热力加热源项的直接起源，只要能在属性面板上选择适当的温度起源，同样焓形式的能量方程也能使用TOTAL类型
我想知道有关于STATIC形式和TOTAL形式的能量方程究竟是什么？？请教各位高手！！
呵呵

frogfish

Q：请问图中所说的coefficient(系数）对气体常设为1000，
对液体设为1；在实际中都是这样做吗？
为什么？
有没有应相应的原则？
在设置出口时XYZ的速度都是0,还是设成in-cell?
A： 压力是相对压力，它相对于下面给的压力值的压力，就是说，如果下面的外部压力为1e05（1个大气压）时，它就是相对于这个压力的值，它的值就代表是表压了。所以要看下面的基准是多少。
系数的设定我也不太懂，但我的经验是，单相流默认就可以了，两相流时将有两个相的系数，这时要注意两个系数比等于两相密度比，这个问题我搞了好几天才试出来的。


A：所有边界条件都是以源项形式加入到PHOENICS里面的，所以outlet里面的系数就是源项的系数。
关于两相流出口的设置，在POLIS的一个星期教程里面说的很清楚。

- if Pext is not zero.

If only one phase is allowed to pass, then a PATCH must be used, with a COVAL for P1 or P2.

PATCH(OUT, HIGH, 1, NX, 1, NY, NZ, NZ, 1, LSTEP)
COVAL(OUT, P1,1E3, Pext)

- if only Phase 1 is allowed to pass

or

COVAL(OUT, P2,1E3, Pext)
- if only Phase 2 is allowed to pass.

frogfish

Q： 急问：RElAX(LTLS,FALSDT,1.00000E+00)是什么意思？是如何设置的？
A：设置松弛因子啊，LTLS是离开壁面的距离，FALSDT是虚拟时间步长法
Q：在VR的numbers什么地方设，才能有效如何设？才能反应到q1文件中？
A：RELAX(variable name,FALSDT,false time step), adds a source in each cell equal to the mass in the cell divided by the size of the time step.
在第17组中设置这个命令的,正的松弛因子是采用虚拟时间步长方法的,最基本的是(可以查看polis中的相关介绍)：
Its effect is to add to the phi-balance for each cell a source equal t
(latest_phi - next_phi) * (mass in cell) / DTFALS(phi)

where latest_phi is the current value of the variable being solved for, and next_phi is that which it will take after adjustment.

This has the effect of making next_phi very close to latest_phi when DTFALS is large, so that changes can occur only after many repetitions of the adjustment process.

A value of DTFALS=1.E-20 would prevent change entirely on many machines, a fact which it is sometimes convenient to exploit, for example when seeking to understand the causes of divergence

frogfish

Q：各位大侠，我想问的是：在一个计算空间内用隔板（导热的blockage或者thinplate）隔开两种流体进行耦合的换热器传热计算问题。两种不同的流体分别从两个进口进入计算空间。我想问的是一种流体可以在main menu 中的property设为一种流体，那么另外一种流体该如何设置其属性？在什么地方可以对其设置？例如：比热系数，导热系数，密度等。
A：
这其实就是一个管壳式换热器的问题,你可以把其中一中流体当作源项来处理
我当初就是在plant里做的,应该可以做到的

frogfish

Q：
在设置流场的时候，如果要设置固体，似乎有两种方法
一种是认为全场都是流体，但是固体位置粘性无限大，这样就体现了固体的存在
另一种是单独设置固体，但是要设置边界，边界要增加

是不是有这么两种情况？
Phoenics了有这两种情况的选择么？还是采用的是哪一种方法？
A：
是有两种方法，一种是常规的模型，即流流道是空的，这样设置整个区域为流体，物体设置为实际的物体。另一种方法是正好相反，流道做成模型物体，而整个区域设置为流体。两种方法的计算结果是一样的。只是对于两相流来说，只能用常规模型来做。
Q：
那怎么样体现使用那种方法呢

有什么选择么？
还是该怎样设置？
A：
没有什么选择，两种方法设置只是区域性质与物体性质设置相反。比如，物体代表流道或流动区域的话，就要把整个区域性质设成实际的固体，如钢什么的，而代表流道的物体就要设置成流场中的流体。反之则刚好相反。
这两种情况在PHOENICS例子中都有，大多数是常规设置，即物体就是物体，流道是空的。相反的例子是一个在弯管中的计算，可以查一查。

frogfish

Q：用proE画的三维模型导入phoenics计算完成后，计算文件用save as a case 保存，将保存的文件烤到别的计算机上，打开后导入的模型显示不出来，而是一个方形体代替，请问如何解决这个问题？ 如果要重新导入的话太麻烦了吧！
A：CAD导入PHOENICS以后的文件保存在\phoenics\d_satlit\d_vrgeom\fromstl\
目录下面,所以你COPY Q1文件的时候也需要把相应的 .DAT文件COPY到新的计算机里面.

frogfish

Q:对于大空间建筑，左右两侧完全对称，如何保证流场的对称？
发现几个实例计算结果不对称，估计是何原因？
如只设一半建筑，请问对称面如何设置边界条件？
A:缺省边界不就是对称边界吗?设一般就行了

frogfish

Q：我将phoenics装在了D盘，其他使用一切正常，只是帮助功能在使用时不能用，在VR环境下，如果用帮助，显示找不到文件C/phoenics/...，可有那位知道是否可有改路径等方法，使可以从使用界面上使用帮助文件？
A：Please do as follows:
VR Menu =>File=> Open file for Editing =>Cham.ini

in the file cham.ini you can modify the

[POLIS]
drive = c ==> drive = * ,here * is your drive in which the POLIS is.

frogfish

Q：phoenic3.5安装在winxp下可以吗　

A：可以
安装的时侯选WIN NT就行
因为win2000 和win xp 都是NT内核

frogfish

Q：我在进行流场计算的时候，遇到了一个问题
为什么流场计算的最后一个网格的速度U总是0
好多例子里也出现了这个问题
谁能帮我解释一下

如何改变这种设置？



A：You are right . But you should know that is the result from staggered grid.
With the grid there is no velocity at the boundary, for example,for 1D problem(N grid):

| | | ...... | | |
IX 1 2 N-1 N <==cell(or grid)

U1 1 2 N-1 N <==U1 position

Here U1 at N(i.e boundary) is no meaning in staggred grid,so
the U1 is 0 at IX=N.


Q：那是不是说第n-1各网格以前的数值都是可以用的正确的数值呢？
A：Yes ,you are right.
Here,the velocity is located at the face of cell.
The value at cell-centre can be got with 0.5*(U1(I-1)+U1(I)) for I cell
if you hope.