DISPLAY

  TWO TUBES IN A SEMI-RECTANGULAR CHAMBER 

  The natural convection around two differently heated tubes
  placed in an adiabatic rectangular chamber with circular roof
  is solved here by PARSOL.

  The case demonstrates the PARSOL 'cut-cell' technique of
  representing curvilinear shapes in a Cartesian grid.

  The case considers the natural convection arising
  in a 2D adiabatic rectangular enclosure. It is covered by the
  roof of cylindrical shape. Two cylindrical tubes are placed at
  the bottom part of the chamber. They have different temperatures
  provoking the gravity-induced circulation in the surrounding air.

  ENDDIS

 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 24/07/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =       1 ;LIBREF =     400
 ************************************************************
  Group 1. Run Title
 TEXT(Two tubes in a semi-rectangular chamber  )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,25,25,1)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 3.200000E-02; YSI= 3.200000E-02; ZSI= 1.400000E+00
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(150) =TEM1
 SOLVE(P1  ,U1  ,V1  ,TEM1)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)

 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  0)
 RHO1    = 1.189000E+00
 RHO1A   = 0.000000E+00 ;RHO1B  = 3.496503E-03
 RHO1C   = 7.142857E-01
 PRESS0  = 1.000000E+05
 TEMP0   = 2.730000E+02
 CP1     = 1.005000E+03
 ENUL    = 1.544000E-05 ;ENUT   = 0.000000E+00
 ENULA   = 1.460000E-06 ;ENULB  = 1.100000E+02
 ENULC   = 0.000000E+00
 DVO1DT  = 3.410000E-03
 PRNDTL(TEM1) = -2.580000E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
   No PATCHes used for this Group
 FIINIT(TEM1)=0.5 * (20.0 + 150.0)

 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      , GRND3       )
 COVAL (BUOYANCY,V1  , FIXFLU      , GRND3       )

 BUOYA   = 0.000000E+00 ; BUOYB =-9.810000E+00
 BUOYC   = 0.000000E+00
 BUOYD   = 3.410000E-03
 BUOYE   = 2.000000E+01
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =    350
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 1.000000E-01)
 RELAX(V1  ,FALSDT, 1.000000E-01)
 RELAX(TEM1,FALSDT, 1.000000E+01)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(TEM1) = 1.000000E+10 ;VARMIN(TEM1) = 2.000000E+01
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 LSG3    =    T
 CSG3    ='FGEM'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(P1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(U1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(V1  ,Y,Y,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      13 ;IYMON  =       6 ;IZMON  =       1
 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
 NOWIPE  =    T

 GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        3.200000E-02, 3.200000E-02, 1.400000E+00
> DOM,    MONIT,       1.550000E-02, 7.500000E-03, 7.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 3.000000E-02
> DOM,    SNAPSIZE,    1.000000E-03
> DOM,    RELAX,       5.000000E-01

> OBJ1,   NAME,        B2
> OBJ1,   POSITION,    1.600000E-02, 1.600000E-02, 0.000000E+00
> OBJ1,   SIZE,        1.600000E-02, 1.600000E-02, 1.400000E+00
> OBJ1,   CLIPART,     hh1
> OBJ1,   ROTATION24,        1
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        BLOCKAGE
> OBJ1,   MATERIAL,      198

> OBJ2,   NAME,        B3
> OBJ2,   POSITION,    0.000000E+00, 1.600000E-02, 0.000000E+00
> OBJ2,   SIZE,        1.600000E-02, 1.600000E-02, 1.400000E+00
> OBJ2,   CLIPART,     hh1
> OBJ2,   ROTATION24,        2
> OBJ2,   GRID,        2.000000E+00
> OBJ2,   TYPE,        BLOCKAGE
> OBJ2,   MATERIAL,      198

> OBJ3,   NAME,        B6
> OBJ3,   POSITION,    4.999997E-03, 7.300000E-03, 0.000000E+00
> OBJ3,   SIZE,        8.000000E-03, 8.000000E-03, 1.400000E+00
> OBJ3,   CLIPART,     cylinder
> OBJ3,   ROTATION24,        1
> OBJ3,   GRID,        2.000000E+00
> OBJ3,   TYPE,        BLOCKAGE
> OBJ3,   MATERIAL,      100
> OBJ3,   FIXED_TMP,     0.000000E+00, 2.000000E+01

> OBJ4,   NAME,        B7
> OBJ4,   POSITION,    3.200001E-03, 1.440001E-02, 0.000000E+00
> OBJ4,   SIZE,        1.600000E-03, 1.600000E-03, 1.400000E+00
> OBJ4,   CLIPART,     cubet
> OBJ4,   ROTATION24,        1
> OBJ4,   VISIBLE,    -1.000000E+00
> OBJ4,   TYPE,        PRESSURE_RELIEF
> OBJ4,   PRES_RELIEF,   1.000000E+03, 0.000000E+00

> OBJ5,   NAME,        B8
> OBJ5,   POSITION,    1.900001E-02, 7.300000E-03, 0.000000E+00
> OBJ5,   SIZE,        8.000000E-03, 8.000000E-03, 1.400000E+00
> OBJ5,   CLIPART,     cylinder
> OBJ5,   ROTATION24,        1
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      100
> OBJ5,   FIXED_TMP,     0.000000E+00, 1.500000E+02
STOP
distil=t
EX(P1  )=   2.133E-02 
EX(U1  )=   1.225E-02 
EX(V1  )=   1.384E-02 
EX(PRPS)=   2.152E+01 
EX(TEM1)=   8.837E+01