PHOENICS菜鸟版相关资料（转贴）

linda

56.关于迭代次数的一个问题
Total number of iteration   与 Iteration contral 里面的LITER有何区别
？？

答：Total number of iteration   迭带次数,  就是SIMPLE里面的迭带次数.
另外一个是内迭带,就是每个SLAB里迭带次数


57.为什么流场计算的最后一个网格的速度U总是0

我在进行流场计算的时候，遇到了一个问题
为什么流场计算的最后一个网格的速度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.


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

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


59.将phoenics装在其他盘下的情况
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.


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

linda

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


62.我用PHOENICS的体会与收获
我是用PHOENICS来做毕业论文，方向是两相射流方面的，由于资料少、例题不详细，加之一外是日语，所以学起来非常难，到现在两相流还没有搞定。不过，经过这一段的应用与练习，还是有一点收获，想说出来供大家参考，也许对你能有一点帮助。
1.坐标系的选择
我觉得，如果你的模型很简单的话，最好用贴体坐标，因为贴体坐标计算的精度还是好些，收敛也快，运算量小。
如果你的模型是轴对称，如管道，轴对称射流等，最好用极坐标，这样运算量大大降低。
如果不是上述两种情况，那只能用直角坐标系了。
2.前处理问题
对于贴体坐标系，当然只有建贴体网格了，可以在PHOENICS提供的环境中直接建立，也可以利用专门的软件生成后调入。有一种软件叫AC3D不错，其界面如同3DMAX，只有不到2M大小，只是下载后不能输出，如果有注册版的话，还是不错的选择。
对于极坐标系，也可以用AC3D，但仍然存在不能输出的问题；在极坐标系下，要注意Ｘ方向代表角度。最好用3DMAX软件做一个二维图形，然后挤压出一定的厚度，最后导出STL格式，就可以调用。这里需要注意两个问题：一是导出时要选择ＡＳＣＩＩ码，另一个是在挤压成三维实体时，一定要在挤压方向给一定的段数。这是我昨天刚刚实现的，现在感觉很爽。
对于直角坐标系来说，生成三维实体的方法很多，可以用PHOENICS模型文件，也可以用autoCAD，还可以用proE，当然也可以用3DMAX。这几中方法我都试过，但是对比来看，还是用proE和3DMAX更好些。PHOENICS的模型文件有限，有时不一定能给合出合适的模型，另外，它的模型很粗糙，如果要求模型较为精细的话，就不行了；autoCAD我试过，但没有成功，因为我的软件不能导出需要的格式。据说很多人用这个；3DMAX制作非常简单，也很方便，只是相对于proE来说，3DMAX在准确度方面有些差。尢其是进行布尔运算时，它常常出错。还用，用3DMAX做的模型，常常不被承认。所以，用proE更好些。
3. 关于后处理
后处理软件中，tecplot和origin不错，tecplot我用的不多，我最的最熟的是origin，这要看每个人的喜好了。
值得提出来的是，运算结果中往往数据量很大，如果用记事本打开一行一列地找就很费事了，有人贴出来俄罗斯人编的处理小程序，我自从用了它之后，效率大增。用它可以把你想要的曲线保存成数据，再用origin调出来，几条曲线画在一起，这样就可以看出对比的结果了。
就这些，是我的收获吧，都告诉大家了，可是我还苦呀，我的两相流和多相流还没有搞定，曾经有一段时间为此而头大大的，差点就把脑子搞坏掉了。后来我一想，不能让脑子坏掉的，因为一坏掉的话，我就不能把我学到的东西告诉大家了，所以脑子就又好了。
相信以后我会把两相流和多相流搞定的，也希望大家相信我。
请记住，只要有信念，就没有学不会的。
好了，我要非常非常吃力地看我的英文资料了


63.如何画等值线（利用photon)
以library case 100为例
load 100
run earth
进入PHOTON，按如下操作：
phi        读入PHI文件
con
temp
x 5         (假如需要画出X=5上的温度等值线）
再回车       表示缺省
val          表示需要画温度等于某个值的曲线
1            画 1 个温度值
0.5          温度等于0.5

结束

如果需要写入0.5
M      进入菜单操作
TEXT
0。5
鼠标点击需要放置0.5的位置。

64.有关phoenics中流场初始化的问题
问：phoenics采用的是simplest算法，simplest的计算步骤与simple相同，只是在simplest中对于对流－扩散项的离散格式做了明确的规定。
                        ――――――――《数值传热学》第二版，陶文铨，224页
   simple算法的计算步骤如下：
   1、假定一个速度分布，u0、v0，以此计算动量立三方程中的系数及常数项。
   2、假定一个压力场p*;
   3、4、5、6、省略
――――――――《数值传热学》第二版，陶文铨，207页
   问题：既然simplest的计算步骤与simple相同，那么在phoenics中，开始计算的假定  值：u0、v0、p*在什么地方设置？？
   是在主“menu”－Initialisation－initial value for each variable中设置吗？
   流场的初值就是SIMPLE算法中的假定值：u0、v0、p*吗？？
答：是的


65.有关斜面设置的问题讨论
问：我用phoenics3.4计算一导热问题，边界有一斜面，为对流换热边界条件
   我尝试直接设定，没有成功
   只好利用很小的阶梯代替，很繁琐
   问问各位老师，这种问题可以不可以直接设定
   如果没有，那么小阶梯代替斜面，有没有标准，多小才可以
答：你把平面转一定角度不就是斜面了吗?
答：使用PARSOL看是否有效果.3.5.1的PARSOL估计没有问题.3.4够戗
答：是不是要设置一个非线形边界条件啊，如果是那样，Q1里面不也是很容易

66.有关多孔物体的设置
对于多孔物体，设置气物性的时候，把有孔的地方的物性设置为domain material即可

linda

67.有关出口流量的设置问题
问：我得模型有很多出口很多入口，
   只设出口其中一个的固定流量。这样可以么？
答：都要设
(d) Fixed values and fixed fluxes
A large C makes the source term the dominant term in the balance equation for PHI, which for large C becomes,
T*C*(V-phiP)=0.0 ,

ie. a large C has the effect of fixing phiP to V. The commands COVAL(name,PHI,2.E10,V), and COVAL(name,PHI,FIXVAL,V) both have this effect, as FIXVAL is numerically equal to 2.0E10.

A very small C and a very large V set the source equal to T*C*V independently of the value of phiP, for example thus: COVAL(name,PHI,1.1E-10,1.1E10*flux) .

COVAL(name,PHI,FIXFLU,flux) has the same effect, because the variable FIXFLU, as well as being a large number (viz 2.0 E10), is also taken by EARTH as a signal multiply flux by 1./FIXFLU when it adds the COVAL term to the source.


(e) Pressure and mass flow boundary conditions
Mass-flow conditions are supplied by way of pressure boundary conditions, eg. COVAL(name,P1,Cp,Vp), sets a mass source ( inflow is positive ),

T * Cp * (Vp-pressP).

A large coefficient fixes the in-cell pressure, pressP, to Vp; a small Cp and a large Vp set a fixed mass flux equal to T*Cp*Vp.

The value, Vphi, of variable PHI convected in to the domain when this mass source is positive is specified by, COVAL(name,PHI,Cphi,Vphi).

This gives a total source for PHI equal to,


  1----------------------------------------------------------1
  1   T*Cphi*(Vphi-phiP) + [T*Cp*(Vp-pressP)]*(Vphi-phiP).   1
  1----------------------------------------------------------1

The first term represents diffusive inflow; the second represents convective inflow. Usually, the diffusive inflow is neglected, so
Cphi is set to 0.0, but,
COVAL(name,PHI,ONLYMS,Vphi),
signifies the same thing, ie. ONLY MaSs flow.

For cells in which mass outflow occurs ( ie. in which pressP is greater than Vp ), the second term above is absent, in agreement with the upwind convention used for the cell-face fluxes generally. Thus, for outflow, the source is,

T * Cphi * (Vphi-phiP).

The magnitude of Cp needed to fix the pressure may be estimated from the following expression, 1.E3 * ( expected flow rate ) / ( T*Vp ). This gives a relative difference between Vp and the in-cell pressure of order 1.E-3.

When the pressure is to be fixed to zero, as is often done in incompressible flows in which the pressure level is immaterial, Vp is omitted from the above expression, which then gives an estimate of Cp needed to give an absolute difference between Vp (which is 0.0) and pressP of order 1.E-3.

A value of Cp 1000 times bigger than that given by the above expression, would result in an absolute pressure difference between zero and press"P" barely representable on a 32-bit machine, and should hence be avoided. The coefficient FIXP has a numerical value of 1.0; it is suitable for fixing pressures subject to its conformity to the above expression. FIXVAL ( =1.E10 ) is usually far too big.

When the pressure is "fixed", small changes resulting from dynamic- head variations can change inflow to outflow and vice versa, with consequent convergence difficulties.

When COVAL is used to prescribe pressure boundary conditions, it should be called with P1 or P2 as the variable argument. In this regard, it is helpful to note that the relationship between pressure and mass flux is of the same kind as that between temperature and heat flux. Thus, for fixed pressure, enter:

COVAL(patch name, P1, FIXP, value of pressure). For a fixed first-phase mass flux, enter: COVAL(patch name,P1,FIXFLU, value of mass flow rate) ; and for fixed second-phase mass flux, enter: COVAL(patch name,P2,FIXFLU, value of mass flow rate) .

IMPORTANT NOTE: The above statement IS correct even though, in the current version of PHOENICS, the second phase shares the same pressure as the first, so that no storage space is ordinarily allocated to it unless specifically requested by the user introducing SOLUTN(P2,Y,N,N,N,N,N) .

In two-phase flows, the formula used for mass outflows (ie. when pressP is greater than Vp ) is multiplied by the phase volume fraction: R1 for phase-1 mass outflow, and R2 for phase-2 mass outflow.

linda

68.关PATCH(name,PROFIL,NX/2,NX/2,NY/2,NY/2,NZ/2,NZ/2,1,LSTEP)
问：PATCH(name,PROFIL,NX/2,NX/2,NY/2,NY/2,NZ/2,NZ/2,1,LSTEP)
   PLOT(name,P1,0.0,0.0)
   PLOT(name,W1,0.0,0.0)
   In a parabolic calculation, the following setting,
   PATCH(name,PROFIL,1,1,1,1,1,NZ,1,1)
   PLOT(name,W1,0.0,0.0)

   这几个语句的作用是什么啊,谢谢了
答：PLOT
---- Autoplot Help ----
PL[OT] [j]
Plot data elements i - j using a solid line to join the data points. If i & j are unspecified, all data elements in memory but not on the screen will be plotted.

See also HELP on : BLOB, BLB, DOT
PLOT
------ Command; group 23 -------------- -

PLOT....command used for plotting field values as CONTURs or PROFILes, according to the PATCH-type set and over the PATCH prescribed. It has four arguments.

(a) CONTUR plots
Contour plots for one or more variables may be elicited over any portion of a plane of constant IX, IY or IZ indicated by the arguments of PATCH, when it has CONTUR as its second argument. In this case, the significance of the 3rd and 4th arguments of PLOT are as follows. The third argument of the corresponding PLOT command is significant only when the CARTES is F ( ie the polar-coordinate grid is employed) and the contours are being plotted for a plane of constant z; then a zero gives a rectangular plot (angle is horizontal, radius vertical ); whereas insertion of 1.0 will produce a polar plot. The fourth argument indicates how many contour intervals are required.

For example, the following commands give contour plots for H1 ( in 20 intervals ) and for C1 ( in 10 intervals ) in one quadrant of the plane at IX=3:

PATCH(name,CONTUR,3,3,NY/2,NY,NZ/2,NZ,1,1) PLOT(name,H1,0.0,20.0) PLOT(name,C1,0.0,10.0)

See CONTUR, NCOLCO and NROWCO for related information.

(b) PROFILe plots
Profile plots for one or more variables may be elicited over any row ( or portion of row ) of cells at constant IX, IY, constant IX, IZ, constant IY,IZ and constant IX,IY,IZ in transient calculations. The row of cells is indicated by the last 8 arguments of PATCH when it has PROFIL as its second argument.

The profiles corresponding to the dependent variables mentioned in the PLOT statements may be plotted either singly or together. In the former case, the values are plotted horizontally, and in the latter they are plotted vertically. IPROF=0 gives the former arrangement, and IPROF=1, 2 or 3 gives the latter. See IPROF, for further information, including the various meanings of the PLOT arguments.

For the default of IPROF=1, the following commands plot the variation with time of P1 and W1 at the cell NX/2, NY/2, NZ/2, on the same plot:


PATCH(name,PROFIL,NX/2,NX/2,NY/2,NY/2,NZ/2,NZ/2,1,LSTEP)
PLOT(name,P1,0.0,0.0)
PLOT(name,W1,0.0,0.0)

In a parabolic calculation, the following setting,
PATCH(name,PROFIL,1,1,1,1,1,NZ,1,1)
PLOT(name,W1,0.0,0.0)
plots the profile of the axial velocity as a function of IZ.

See also PROFIL and IPROF for related information, and ABSIZ and ORSIZ for control of the box size printed.

Many examples of the use of PROFIL and CONTUR line-printer plots are to be found in the PHOENICS Input Library.

linda

69.设置进口压力的方法
如果要设置进口的压力，可有以下几个方法：
1.把入口定义成outlet
2.用PATCH和COAVL来实现：
PATCH(name,type,ixf,ixl,iyf,iyl,izf,izl,timef,timel)
COVAL(name,P1,FIXVAL,P1in)
3.VR中的user define 中通过设置合适的COEF 和value 就可以也可实

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

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

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

71.有关湿度问题的问答
问：温度场和流场已经算出来了，不知道怎么加湿度？
答：1首先设置流体是混合的理想气体，即干空气和水蒸气的混合气体；
   2添加湿度变量，设置水份蒸发边界条件；
   3还可以计算水蒸气分压力和相对湿度等变量。
问：要设置的流体再哪儿？
答：在理想气体中选择混合气体，再修改密度，
   输入各组分的气体常数即可。
问：湿度变量怎么添加？
答：湿度变量在MODEL里面设置.
   建议看一下www.phoenics.cn里面的FLAIR新版本有关内容.


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

73.有关源相设置的问题
在phoenics中源项和边界条件都是通过PATH,COVAL命令添加的。
比如说要模擬一個長方形腔体,在長方形的腔体上方,有一個孔,再由孔噴入气体,
要观察气体的浓度变化，此时源相就可以通过如下的方法来设置：
在设定喷入气体的入口处加上一个coval(name,c1,fixval,1.0）就可以了

74.RElAX(LTLS,FALSDT,1.00000E+00)是什么意思？是如何设置的？
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 to:
(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

linda

75.有关object attribute里设置系数的选项


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.

linda

76.温度设置中static/tota选项的区别
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形式的能量方程究竟是什么？？请教各位高手！！
呵呵


77.关于出口条件的设置
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

78.在Q1文件中的某些项，在VR中对应什么设置
Q:1）group9中，EMISS=0.000000E+00;SCATT=0.000000E+00是在VR中什么地方设置的，是什么意思？
2）group15中，LSWEEP=100是在VR中什么地方设置的，是什么意思？
3）group17中，RELAX（V1，FALSDT，1。00E-0。1）是在VR中什么地方设置的，是什么意思？
A：1)Models/Radiation model settings.
  2)Numerics/Total number of iterations
  3)Numerics/Relaxation control settings


79.PHOENICS文件夹下各子文件夹的内容
What is the PHOENICS directory structure ?

A brief description of the first level subdirectories under /phoenics is given below.

d_manage is a global interface enabling users to access and display HTML files, activate PHOENICS executables and activate PHOENICS utilities (Windows only).

d_allpro contains the binary code for the PHOENICS service routines shared by SATELLITE and EARTH, the program libraries for all modules and the graphics drivers used by SATELLITE, EARTH and PHOTON;

d_chmkin, containing files for the CHEMKIN interface.

d_earth contains the binary code and some source files of the PHOENICS main solver program EARTH, together with the PHOENICS case libraries;

d_enviro contains program and data files for the PHOENICS COMMANDER;

d_includ contains the Fortran common blocks required for recompilation;

d_intfac contains the interfaces between PHOENICS and other codes;

d_photon contains the binary code of the PHOENICS post-processor program PHOTON, which now incorporates AUTOPLOT;

d_polis contains all files for the PHOENICS On-Line Information System;

d_priv1 is referred to as a "private" working directory from which the user may run all the PHOENICS programs, or use the PHOENICS-VR Environment;

d_satell contains the binary code of the PHOENICS pre-processor program SATELLITE, and the menu libraries;

d_shapem contains the executable of Shapemaker, a geometry generator (Windows only);

d_utils contains utility batch files, such as run scripts, link scripts, and compiler option scripts. Also included are directories, which contain many utility programs.

Please note that the files used only for Compaq Visual Fortran compiler will be in the subdirectory d_windf and those only for Salford Fortran compiler will be in the subdirectory d_winsd in each of the above-mentioned directories. The 'df' is a reminder of "Digital Fortran (now Compaq Fortran); and the 'sd' is a reminder of "Salford DBOS".

In addition to phoenics, the installation process will create the directory website, containing the additional files relevant to the PHOENICS On-Line Information System.
polis 的缺点就是要讲的东西太多，链接太多，看看就转晕了啊～

80.如何设置自己想要的材料
在VR环境中编辑事，如果自己想要的材料没有怎么办呢？？
用写字板打开这个文件phoenics\d_earth\props
然后自己编辑吧~~~~
气体，液体，固体
注意编号。

0--30  气体， 30---40  空的（留给用户添加气体物性）

其他的也类似。

然后重启PHOENICS
OK了~~~~~

81.ProE与PHOENICS
将ProE图形调入PHOENICS的步骤
(1). 主菜单上File-Export-Model, 选STL格式，注意：STL菜单中format选ASCII,不能选Binary格式。最后取一个文件名。还需要注意的是：单位问题，必须设置好单位，与Phoenics统一。
2. Phoenics中调入STL文件
   在Object中的CAD interface 一栏选STL file , 调入刚存的Pro/E文件。需要主意的是Domin 的比列问题。

注意：只能传入实体,在vr下用

linda

82.Y400的运行
Y400这个例子是专用模块中ESTER的，应该先运行satlite,再运行earth, 在commander下可一运行,在vr下不能运行.


存储空间不够的问题
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

linda

83.phoenics3.6
The many changes made in PHOENICS since the release of V3.5.1 include:
&#8226; Extension of the PHOENICS object concept to permit non-geometric and geometric attributes to be stored together
&#8226; Provision of a ‘guaranteed convergence’ feature
&#8226; Further generalization of the solid-stress feature
&#8226; Provision of numerous new features in the graphical user interface
&#8226; Many improvements to FLAIR, the specialpurpose program for HVAC and fires in buildings
&#8226; Extensions to the PHOENICS Commander
&#8226; A re-working of Shapemaker which can now add and store non-geometric objects, as well as creating their shapes
&#8226; Use of dynamic memory allocation in all modules
&#8226; Restoration and extension of parallel-PHOENICS capabilities.
The following object types have been added to the
existing list:
&#8226; “General_plate” (zero thickness, flat or curved)
&#8226; “Drag_lift” object (drag and lift calculation, 3D)
&#8226; “Assembly” object (master and related
components structure)
[The new “Assembly” object has been introduced to act
as a ‘container’ for any number of other objects.]
PHOENICS-3.6 is available in Beta form, with more
widespread release expected during Autumn 2004.
Flair Enhancements
Many new and improved features have been added to the
special-purpose FLAIR (building services) module of
PHOENICS:
&#8226; A “Person” object was created for a single
person acting as source of heat, humidity and
any solved scalars.
&#8226; Similarly, the “People” object was created to
represent large numbers of people.
&#8226; Diffuser (mass an heat source/sink, 3D)
&#8226; Fire object (piece-wise linear heat, mass and
smoke sources for transient calculations)
The Comfort-index calculation: PPD, PMV, TRES, Mean
Age of Air.
FLAIR provides for three comfort index calculations:
&#8226; Dry resultant temperature (TRES) (CIBSE Guide)
&#8226; Predicted mean vote (PMV) (ISO 7730)
&#8226; Predicted percentage dissatisfied (PPD)
(ISO 7730)
What’s New in PHOENICS-3.6
PHOENICS NEWS Autumn 2004 6
Standing - and sitting - person objects
Help entry on comfort indices
FLAIR provides a humidity calculation, including output of
humidity ratio and relative humidity. The solved-for
specific humidity equation, MH2O, has units of kg water
vapour/kg mixture. It is a mass fraction of water vapour.
Several derived quantities can be displayed. These are:
&#8226; Humidity ratio, which has units of (g/kg)
&#8226; Relative humidity (%)
[In this case, the water-vapour saturation pressure, partial
pressure and mole fraction are also made available for
storage.
The units used to specify boundary sources (inlets,
openings, volume sources) can be set on the dialogs for
the individual objects.
The options are:
&#8226; Specific Humidity (Mass fraction) (kg/kg, kg/s)
&#8226; Humidity ratio (g/kg, g/s)
&#8226; Relative humidity (%)
[The default is Specific Humidity]
The default Up direction is set to Z instead of Y as in
earlier versions.
The default domain size is set to 10 x 10 x 3m to be more
representative of an average room or office shape and
size. The automatic convergence control is set as the
default option.
A more comprehensive description of these, and other
improvements, can be found by clicking on:
www.cham.co.uk/phoenics\d_polis\d_docs\tr006\item6.htm
What’s New in PHOENICS-3.6 - continued
Melbourne Cricket Ground
Pictures featured in plot of the month - TECPLOT Inc - wwwtecplot.com
www.cham.co.uk/phoenics/d_polis/d_docs/tr006/tr006.htm

linda

84.PHOENICS读取AUTOCAD 文件步骤：
1、在AUTOCAD下生成“三维实体”，X、Y、Z必需为正值。
2、在命令状态输入：STLOUT
3、SELECT OBJECT
4、不要保存成二进制
5、保存到\PHOENICS\D_PRIV1目录
6、进入PHOENICS程序，读取AUTOCAD文件
7、结束

我在用CAD画图之后，导出图形的时候，采用的是菜单里的EXPORT命令，结果就出现了我上面的错误。
后来经高人指点， 采用命令形式
“ command状态输入命令:
stlout

AutoCAD文件存为STL格式，注意选ASCII码，不要选BIN格式，同时注意单位转换问题，可以传入Phoenics中。”
经过这样的处理就好用了。后来想想，可能是EXPORE命令导出的时候就是ASCII码。


特点：

PHOENICS软件是世界上第一套计算流体与计算传热学商用软件，它是Parabolic Hyperbolic Or Elliptic
Numerical Integration Code Series 几个字母的缩写，这意味着只要有流动和传热都可以使用PHOENICS程序来模
拟计算。除了通用计算流体/计算传热学软件应该拥有的功能外PHOENICS软件有自己独特的功能：
　　 1、开放性：PHOENICS最大限度地向用户开放了程序，用户可以根据需要任意修改添加用户程序、用户模型。
PLANT及INform功能的引入使用户不再需要编写FORTRAN源程序，GROUND程序功能使用户修改添加模型更加任意、方
便。
　　 2、CAD接口：PHOENICS可以读入任何CAD软件的图形文件。
　　 3、MOVOBJ：运动物体功能可以定义物体运动，避免了使用相对运动方法的局限性。
　　 4、大量的模型选择：20多种湍流模型，多种多相流模型，多流体模型，燃烧模型，辐射模型。
　　 5、提供了欧拉算法也提供了基于粒子运动轨迹的拉格朗日算法。
6、计算流动与传热时能同时计算浸入流体中的固体的机械和热应力。
7、VR（虚拟现实）用户界面引入了一种崭新的CFD建模思路。
8、PARSOL（CUT CELL) ：部分固体处理。
9、软件自带1000多个例题，附有完整的可读可改的原始输入文件。
　　 10、PHOENICS专用模块

PHOENICS软件以低速热流输送现象为主要模拟对象，由于长期积累以及Spalding在建立模拟模型上非凡的创造力，
PHOENICS包含的湍流模型、多相流模型、燃烧与化学反应模型等相当丰富，其中有不少原创性的成分，如将湍流与
层流成分假设为两种流体的双流体湍流模型、专门为组件杂阵的狭小空间（如计算机箱体）内的流通和传热计算而
设计的代数湍流模型LVEL等都是Spalding 与其它作者共同提出的。
PHOENICS的边界条件设置也有一定的特点，是以源项的方式给定的。这个软件附带了从简到繁的大量算例，一般的
工程应用问题几乎都可以从中找到相近的范例，在做一些修改就可以用以计算用户的课题，所以能给用户带来极大
方便。


PHOENICS的暖通空调计算模块FLAIR被广泛应用，也被一些别的软件包采纳，如英国集成环境公司（IES）的虚拟环
境软件，一用它来模拟局部空间的热流现象。由于PHOENICS以压力校正法为基本的解法，从而不大适应高速可压流
计算。PHOENICS采用非正交贴体网格，但网格畸变较大时可能会发生困难，估计是在算法中采用交错网格离散的缘
故。


85.Phoenics
简介：

PHOENICS程序是世界著名的计算流体与计算传热学(CFD/NHT)软件，它是英国皇家学会D.B.SPALDING教授及40多位
博士20多年心血的典范之作。PHOENICS已广泛应用于航空航天、船舶、汽车、暖通空调、环境、能源动力、化工等
各个领域。在核电方面，利用PHOENICS不仅可节约大量经费，更为核电的安全可靠运行提供了可靠保证。

PHOENICS是Parbolic,Hyperbolic or Ellicpic Numerical Integration Code Series的缩写。它可以用来模拟流
体流动、传热、化学反应及相关现象。程序有前处理、求解器、后处理模块构成. PHOENICS程序语言是标准ANSI
FORTRAAN77语言，与机器无关，程序总共大约110,000条语句，2000个子程序。

Joy726

谢谢^_^

mengjiale

真的很棒。

ChaChing

像这种占25层的东西, 既然是转载, 是否直接给个连接即可!?

砸牛顿的苹果

谢谢楼主的慷慨分享