phoenics3.5中对气液两相流的自由表面如何设定？？

飞龙在天

请教专家
phoenics3.5中对气液两相流的自由表面如何设定？？
比如说，在一个上面开口的容器底部吹气，气体从上口逸出。其自由表面条件如何设定？
请回复ht-hu@mail.china.com
附我的q1文件
TALK=T;RUN( 1, 1)
************************************************************
  Q1 created by VDI menu, Version 3.5, Date 07/03/03
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN   =         1 ;LIBREF =         0
************************************************************
Group 1. Run Title
TEXT(No title has been set for this run.     )
************************************************************
Group 2. Transience
STEADY  =    T
************************************************************
Groups 3, 4, 5  Grid Information
   * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,20,20,20)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    F
   * Non-default variable names
NAME(146) =ENUT ; NAME(147) =EL1
NAME(148) =GENK ; NAME(149) =CFIP
NAME(150) =EPKE
   * Solved variables list
SOLVE(P1  ,U1  ,U2  ,V1  ,V2  ,W1  ,W2  ,R1  )
SOLVE(R2  )
   * Stored variables list
STORE(EPKE,CFIP,GENK,EL1 ,ENUT)
   * Additional solver options
SOLUTN(P1  ,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
SETPRPS(1,  0)
SETPRPS(2, 67)
DVO1DT  = 3.330000E-03
DVO2DT  = 1.180000E-04
PRT   (EP  ) =  1.314000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
CFIPS   = GRND7
RLOLIM  = 1.000000E-03 ;CMDOT  = 0.000000E+00
CFIPA   = 1.000000E-02 ;CFIPB  = 5.000000E-03
CFIPC   = 2.000000E-02 ;CFIPD  = 4.000000E+00
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(R1  ) =  5.000000E-01 ;FIINIT(R2  ) =  5.000000E-01
  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 (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,U1  , FIXFLU      , GRND1       )
COVAL (BUOYANCY,V1  , FIXFLU      , GRND1       )
COVAL (BUOYANCY,W1  , FIXFLU      , GRND1       )
BUOYA   = 0.000000E+00 ; BUOYB = 0.000000E+00
BUOYC   =-9.800000E+00
EGWF    =    T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =         1
RESFAC  = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (R1  ) =        1 ;LITER (R2  ) =        1
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 1.000000E+00)
RELAX(KE  ,LINRLX, 5.000000E-01)
RELAX(EP  ,LINRLX, 5.000000E-01)
KELIN   =         3
************************************************************
Group 18. Limits
VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
VARMAX(U2  ) = 1.000000E+06 ;VARMIN(U2  ) =-1.000000E+06
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(V2  ) = 1.000000E+06 ;VARMIN(V2  ) =-1.000000E+06
VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
VARMAX(W2  ) = 1.000000E+06 ;VARMIN(W2  ) =-1.000000E+06
VARMAX(R1  ) = 1.000000E+00 ;VARMIN(R1  ) = 1.000000E-06
VARMAX(R2  ) = 1.000000E+00 ;VARMIN(R2  ) = 1.000000E-06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
GENK    =    T
ASAP    =    T
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =        10 ;IYMON  =        10 ;IZMON  =        10
NPRMON  =    100000
NPRMNT  =         1
TSTSWP  =        -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =    100000
ISWPRF  =         1 ;ISWPRL =    100000
  No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.845450E-01, 4.750000E-01, 4.750000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> OBJ,    NAME,        B1
> OBJ,    POSITION,    3.600000E-01, 4.500000E-01, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 1.000000E-01, 0.000000E+00
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       0.000000E+00
> OBJ,    DENSITY-2,     9.982300E+02
> OBJ,    VELOCITY,      0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    TEMPERATURE,   0.000000E+00
> OBJ,    TEMPERATURE-2, 0.000000E+00
> OBJ,    TURB-INTENS,   5.000000E+00
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   SAME
> OBJ,    COEFFICIENT,   1.000000E+03
> OBJ,    VELOCITY,      SAME        , SAME        , 0.000000E+00
> OBJ,    TEMPERATURE-2, SAME
> OBJ,    COEFFICIENT-2, 1.000000E+06
> OBJ,    VELOCITY-2,    SAME        , SAME        , SAME
> OBJ,    TURBULENCE,    SAME        , SAME
STOP

merryfish

需要由计算自己自动决定