向版主请教

rainbowyoung

[这个贴子最后由rainbowyoung在 2004/08/30 10:36am 第 1 次编辑]

版主你好，我是做火灾的，
最近在处理边界条件的时候，
发现将壁面留出，不予任何设定
和设置一个plate，之间有一定的差别，
设置了plate以后，紧靠壁面的速度不真实地增大，
我的第一个网格点是布置在湍流旺盛区的。
请问默认的区域边界是不是就相当于一个“CONPOR”呢，
什么情况下应该再加入一个“WALL”(Q1中）或是“PLATE”(VR中）？
这些问题我一直都不太明白，还请版主不吝赐教，先谢过了！

rainbowyoung

[这个贴子最后由rainbowyoung在 2004/08/30 10:41am 第 1 次编辑]

还有一个问题
我在用Q1编制程序的时候，
将CONPOR 的六个网格起迄之都设为相应的负值，
但是却并没有自动生成应有的壁面函数，
请问这是为什么呢
我用的是3.5版本

rainbowyoung

我学习PHOENICS是从Q1入门的，
现在逐渐向VR过渡，但是还是有一些Q1命令和VR对应不上
所以
麻烦您多费心了！

phoenics

那你应该是高手.
为什么用CONPOR?

rainbowyoung

版主这么说
我实在是不敢当的
使用CONPOR建立各种房间的墙体和挡块
是我的学习环境如此
所以习惯也就是这样了
我不适应长时间的盯着电脑在VR中建立各种所需的模块
所以更加喜欢动手写写Q1文件
这种方式虽然比较怀旧
但是有时却可以更完全的表达出我的意思
我知道现在这样的Q族很少
更多的是VR一族
但是这个问题实在是困扰了我一阵子了
所以
还请版主指点一二
或者你愿意推荐一种更加行之有效的方法呢？

rainbowyoung

还有啊
我的第一个问题
您有什么见解呢
——俺小心翼翼地问……

phoenics

您也太客气了.俺主要是吃不准,所以没有敢回答.PLANT加入以后和缺剩设置实际上差别了一个壁面函数的作用.设置了PLATE壁面函数也自动调用了.至于为什么速度会增大,我觉得还是计算误差.
你做的是烟气流动还是燃烧?
如果是烟气流动我可以一起讨论.如果是燃烧,本人确实是外行,需要请教高人.
如果能提供Q1文件就比较容易讨论了.

sfsm

下面引用由phoenics在 2004/09/04 10:18am 发表的内容：
您也太客气了.俺主要是吃不准,所以没有敢回答.PLANT加入以后和缺剩设置实际上差别了一个壁面函数的作用.设置了PLATE壁面函数也自动调用了.至于为什么速度会增大,我觉得还是计算误差.
你做的是烟气流动还是燃烧?
...

按照phoenics斑竹所说的，如果不设置plate的话，那么在group13里的EGWF应该是默认值F啊，但是，我以前所建立的冷态模型中并没有设置PLATE，但在Q1文件中的GROUP13里的EGWF怎么会是T呢？

rainbowyoung

多谢版主的答复
那么在壁面加入plate以后
预报的准确性是否应该得到提高呢
我是指在计算正确的前提下
我主要是做烟气流动的，正在向燃烧迈进，
但是对燃烧模型的理解一点都不够，所以步履维艰。这次我的准备不够充分
下次把q1一起带过来，共同探讨一下吧
顺便说一句，原来你们周末不休息啊，真的很敬业！

rainbowyoung

EGWF的默认值是F,
sfsm你可能是设置了其它的命令，自动激活了EGWF
也许是这样吧

rainbowyoung

版主你还过来吗？
这是我的q1文件
俺现在潜心修炼烟气流动的模拟
可以一起探讨一下吗？
另外
你知道怎样做开口处流入和流出气体的流量吗？
我想大概要用积分来做
也听公司的代表人员说过
但是没有详细的深入了解

rainbowyoung

TALK=T;RUN( 1, 1)

************************************************************
   Q1 created by VDI menu, Version 3.5, Date 08/10/02
CPVNAM=VDI;SPPNAM=FLAIR
************************************************************
IRUNN   =       1 ;LIBREF =       0
************************************************************
  Group 1. Run Title
TEXT(tanaka350,transient                     )
************************************************************
  Group 2. Transience
STEADY  =    T
************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,64,22,29)
************************************************************
  Group 6. Body-Fitted coordinates
************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
    * Non-default variable names
NAME( 38) =HPOR ; NAME( 39) =NPOR
NAME( 42) =ABSV ; NAME( 43) =BLOK
NAME( 44) =PRPS ; NAME( 46) =DEN1
NAME( 47) =TEM1 ; NAME(148) =EPOR
NAME(149) =VPOR ; NAME(150) =SPH1
    * Solved variables list
SOLVE(P1  ,U1  ,V1  ,W1  ,C1  ,TEM1)
    * Stored variables list
STORE(SPH1,VPOR,EPOR,DEN1,PRPS,BLOK,ABSV,NPOR)
STORE(HPOR)
    * Additional solver options
SOLUTN(P1  ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
IVARBK  =      -1 ;ISOLBK =       1
TURMOD(KEMODL)

************************************************************
  Group 8. Terms & Devices
ISOLX   =       0 ;ISOLY  =       0 ;ISOLZ  =       0
DENPCO  =    T
************************************************************
  Group 9. Properties
PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
RHO1    = FILE
ENUL    = FILE
CP1     = FILE
DVO1DT  = 3.330000E-03
PRNDTL(TEM1) =  CONDFILE
    * Non-default SOLUTN setting for KE
SOLUTN(KE  ,Y,Y,N,N,N,Y)
    * Non-default SOLUTN setting for EP
SOLUTN(EP  ,Y,Y,N,N,N,Y)
PRT   (EP  ) =  1.314000E+00
EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
************************************************************
  Group 10.Inter-Phase Transfer Processes
************************************************************
  Group 11.Initialise Var/Porosity Fields
FIINIT(KE  ) =  7.225000E-03 ;FIINIT(EP  ) =  1.007000E-02
FIINIT(C1  ) =  1.001000E-10 ;FIINIT(HPOR) =  1.000000E+00
FIINIT(NPOR) =  1.000000E+00 ;FIINIT(ABSV) =  1.001000E-10
FIINIT(BLOK) =  1.000000E+00 ;FIINIT(PRPS) =  2.000000E+00
FIINIT(DEN1) =  1.161000E+00 ;FIINIT(TEM1) =  2.000000E+01
FIINIT(EPOR) =  1.000000E+00 ;FIINIT(VPOR) =  1.000000E+00
FIINIT(SPH1) =  1.001000E-10
   No PATCHes used for this Group


INIADD  =    F
************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
************************************************************
  Group 13. Boundary & Special Sources

PATCH (FIRE1   ,SOUTH ,11,0,0,0,0,0,1,10)
COVAL (FIRE1   ,C1  , FIXVAL      , 1.000000E+00)
COVAL (FIRE1   ,TEM1, FIXFLU      , 5.900000E+05)

PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,U1  , FIXFLU      , GRND2       )
COVAL (BUOYANCY,V1  , FIXFLU      , GRND2       )
COVAL (BUOYANCY,W1  , FIXFLU      , GRND2       )

BUOYA   = 0.000000E+00 ; BUOYB =-9.810000E+00
BUOYC   = 0.000000E+00
BUOYD   = 1.189000E+00
EGWF    =    T
************************************************************
  Group 14. Downstream Pressure For PARAB
************************************************************
  Group 15. Terminate Sweeps
LSWEEP  =    5000
RESFAC  = 1.000000E-02
************************************************************
  Group 16. Terminate Iterations
LITER (TEM1) =       15
ENDIT (TEM1) =  GRND1
************************************************************
  Group 17. Relaxation
RELAX(P1  ,LINRLX, 1.000000E+00)
RELAX(U1  ,FALSDT, 1.000000E-02)
RELAX(V1  ,FALSDT, 1.000000E-02)
RELAX(W1  ,FALSDT, 1.000000E-02)
RELAX(KE  ,FALSDT, 1.000000E-02)
RELAX(EP  ,FALSDT, 1.000000E-02)
RELAX(DEN1,LINRLX, 7.000000E-01)
RELAX(TEM1,FALSDT, 1.000000E+02)
KELIN   =       1
************************************************************
  Group 18. Limits
VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
************************************************************
  Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
NAMGRD  =FLAR
GENK    =    T
ASAP    =    T
IDISPB  =       1 ;IDISPC =      10
CSG1    ='m'
CSG2    ='m'
************************************************************
  Group 20. Preliminary Printout
ECHO    =    T
************************************************************
  Group 21. Print-out of Variables
OUTPUT(SPH1,N,N,Y,Y,N,N)
************************************************************
  Group 22. Monitor Print-Out
IXMON   =       7 ;IYMON  =       7 ;IZMON  =      20
NPRMON  =  100000
NPRMNT  =       1
TSTSWP  =      -1
************************************************************
  Group 23.Field Print-Out & Plot Control
NPRINT  =  100000
NXPRIN  =       1
NYPRIN  =       1
NZPRIN  =       4
NPLT    =       1
ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
************************************************************
  Group 24. Dumps For Restarts
NOWIPE  =    T
IDISPA  =       1 ;IDISPB =       1 ;IDISPC =      10
CSG1    ='m'
CSG2    ='m'

GVIEW(P,2.388419E-01,5.271721E-01,8.155023E-01)
GVIEW(UP,-1.296399E-01,8.496006E-01,-5.112460E-01)

> DOM,    SIZE,        8.649999E+00, 2.170000E+00, 3.451000E+00
> DOM,    MONIT,       7.875000E-01, 6.500000E-01, 2.235000E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,      1, 1.000000E+00
> GRID,   RSET_X_2,     10, 1.000000E+00
> GRID,   RSET_X_3,      5, 1.000000E+00
> GRID,   RSET_X_4,     12, 1.000000E+00
> GRID,   RSET_X_5,      1, 1.000000E+00
> GRID,   RSET_X_6,      1, 1.000000E+00
> GRID,   RSET_X_7,     23, 1.000000E+00
> GRID,   RSET_X_8,      1, 1.000000E+00
> GRID,   RSET_X_9,     10, 1.000000E+00
> GRID,   RSET_Y_1,      1, 1.000000E+00
> GRID,   RSET_Y_2,      1, 1.000000E+00
> GRID,   RSET_Y_3,     15, 1.000000E+00
> GRID,   RSET_Y_4,      1, 1.000000E+00
> GRID,   RSET_Y_5,      3, 1.000000E+00
> GRID,   RSET_Y_6,      1, 1.000000E+00
> GRID,   RSET_Z_1,      6, 1.000000E+00
> GRID,   RSET_Z_2,      4, 1.000000E+00
> GRID,   RSET_Z_3,      1, 1.000000E+00
> GRID,   RSET_Z_4,      6, 1.000000E+00
> GRID,   RSET_Z_5,      2, 1.000000E+00
> GRID,   RSET_Z_6,      4, 1.000000E+00
> GRID,   RSET_Z_7,      6, 1.000000E+00

> OBJ,    NAME,        X1
> OBJ,    POSITION,    1.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.170000E+00, 3.451000E+00
> OBJ,    CLIPART,     wirexyz
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        NULL

> OBJ,    NAME,        X4
> OBJ,    POSITION,    3.475000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.170000E+00, 3.451000E+00
> OBJ,    CLIPART,     wirexyz
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        NULL

> OBJ,    NAME,        Y5
> OBJ,    POSITION,    0.000000E+00, 2.100000E+00, 0.000000E+00
> OBJ,    SIZE,        8.650001E+00, 0.000000E+00, 3.451000E+00
> OBJ,    CLIPART,     wirexyz
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        NULL

> OBJ,    NAME,        DOOR11
> OBJ,    POSITION,    3.600000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        4.999995E-02, 1.700000E+00, 1.280000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        DOOR12
> OBJ,    POSITION,    3.600000E+00, 0.000000E+00, 2.171000E+00
> OBJ,    SIZE,        4.999995E-02, 1.700000E+00, 1.280000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        DOOR13
> OBJ,    POSITION,    3.600000E+00, 1.700000E+00, 0.000000E+00
> OBJ,    SIZE,        4.999995E-02, 4.700003E-01, 3.451000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        DOOR21
> OBJ,    POSITION,    7.099999E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000067E-02, 1.800000E+00, 7.680000E-01
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        DOOR22
> OBJ,    POSITION,    7.099999E+00, 0.000000E+00, 2.683000E+00
> OBJ,    SIZE,        5.000067E-02, 1.800000E+00, 7.680001E-01
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        DOOR23
> OBJ,    POSITION,    7.099999E+00, 1.800000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000067E-02, 3.699995E-01, 3.451000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        FLOOR1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        7.150000E+00, 1.000000E-01, 3.451000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199

> OBJ,    NAME,        FIRE1
> OBJ,    POSITION,    1.350000E+00, 2.000000E-01, 1.379000E+00
> OBJ,    SIZE,        6.250000E-01, 0.000000E+00, 5.939999E-01
> OBJ,    CLIPART,     default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED

> OBJ,    NAME,        DOOR2
> OBJ,    POSITION,    7.150000E+00, 0.000000E+00, 7.679999E-01
> OBJ,    SIZE,        0.000000E+00, 1.800000E+00, 1.915000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   2.000000E+01
> OBJ,    COEFFICIENT,   1.000000E+00
> OBJ,    TURBULENCE,   -1.026000E+04,-1.026000E+04

> OBJ,    NAME,        WFN1
> OBJ,    POSITION,    3.600000E+00, 1.700000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 4.700003E-01, 3.451000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    POROSITY,      0.000000E+00
> OBJ,    SIDE,        HIGH

> OBJ,    NAME,        WFN2
> OBJ,    POSITION,    7.150000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E+00, 0.000000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube13
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ,    NAME,        OT1
> OBJ,    POSITION,    8.650000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.170000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   2.000000E+01
> OBJ,    COEFFICIENT,   2.000000E+10
> OBJ,    TURBULENCE,    0.000000E+00, 0.000000E+00

> OBJ,    NAME,        OT2
> OBJ,    POSITION,    7.150000E+00, 2.170000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E+00, 0.000000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   2.000000E+01
> OBJ,    COEFFICIENT,   2.000000E+10
> OBJ,    TURBULENCE,    0.000000E+00, 0.000000E+00

> OBJ,    NAME,        OT3
> OBJ,    POSITION,    7.150000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E+00, 2.170000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   2.000000E+01
> OBJ,    COEFFICIENT,   2.000000E+10
> OBJ,    TURBULENCE,    0.000000E+00, 0.000000E+00

> OBJ,    NAME,        OT4
> OBJ,    POSITION,    7.150000E+00, 0.000000E+00, 3.451000E+00
> OBJ,    SIZE,        1.500000E+00, 2.170000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   2.000000E+01
> OBJ,    COEFFICIENT,   2.000000E+10
> OBJ,    TURBULENCE,    0.000000E+00, 0.000000E+00

> OBJ,    NAME,        WALL1
> OBJ,    POSITION,    0.000000E+00, 2.170000E+00, 0.000000E+00
> OBJ,    SIZE,        7.150000E+00, 0.000000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE

> OBJ,    NAME,        WALL2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        8.650000E+00, 0.000000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE

> OBJ,    NAME,        WALL3
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.170000E+00, 3.451000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE

> OBJ,    NAME,        WALL11
> OBJ,    POSITION,    3.600000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.700000E+00, 1.280000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    POROSITY,      0.000000E+00
> OBJ,    SIDE,        BOTH

> OBJ,    NAME,        WALL12
> OBJ,    POSITION,    3.600000E+00, 0.000000E+00, 2.171000E+00
> OBJ,    SIZE,        0.000000E+00, 1.700000E+00, 1.280000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    POROSITY,      0.000000E+00
> OBJ,    SIDE,        BOTH
STOP