hot121

Library case H121:


talk=t;run(1,1)
  PHOTON USE
  p




  use patgeo
  up .54 .84 0
  red
  scal
  msg
  msg Temperature contours
  msg
  cont tem1 z 1 fil;.01
  msg Press  to continue
  pause
  msg Velocity field
  vec z 1 sh
  msg
  MSG Type E to end
  ENDUSE
  DISPLAY

  Library case H121:

  THE EFFECT OF ENCLOSURE ORIENTATION

  This case models simple  two-dimensional  flow  through  a
  tilted box (angle is 1 rad).  Air has  two  apertures  for
  flow  and  the  buoyancy  force  is  generated  from   the
  temperature difference;  single solid in the middle of the
  domain generates heat.

  ENDDIS
 ************************************************************
 CPVNAM=VDI;SPPNAM=HOTBOX
 ************************************************************
 IRUNN   =       1 ;LIBREF =     121
 ************************************************************
  Group 1. Run Title
 TEXT(The Effect Of Enclosure Orientation
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,10,10,1)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
 RSET(D,BOX     )
    * Set objects: x0     y0     z0
    *              dx     dy     dz
                   name
 XPO= 4.000000E-01; YPO= 4.000000E-01; ZPO= 0.000000E+00
 XSI= 2.000000E-01; YSI= 2.000000E-01; ZSI= 1.000000E+00
 RSET(B,BLK     )
 XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
 XSI= 0.000000E+00; YSI= 2.000000E-01; ZSI= 1.000000E+00
 RSET(B,APT1    )
 XPO= 8.000000E-01; YPO= 1.000000E+00; ZPO= 0.000000E+00
 XSI= 2.000000E-01; YSI= 0.000000E+00; ZSI= 1.000000E+00
 RSET(B,APT2    )
    * Modify default grid
 RSET(X,1,4,1.000000E+00,G)
 RSET(X,2,2,1.000000E+00)
 RSET(X,3,2,1.000000E+00)
 RSET(X,4,2,1.000000E+00)
 RSET(Y,1,2,1.000000E+00,G)
 RSET(Y,2,2,1.000000E+00)
 RSET(Y,3,2,1.000000E+00)
 RSET(Y,4,4,1.000000E+00)
 ************************************************************
  Group 6. Body-Fitted coordinates
    * X-cyclic boundaries switched
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME( 46) =BLOK ; NAME( 47) =PRPS
 NAME( 49) =DEN1 ; NAME( 50) =TEM1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,TEM1)
    * Stored variables list
 STORE(DEN1,PRPS,BLOK)
    * Additional solver options
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 IVARBK  =      -1 ;ISOLBK =       1
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T
 NEWENL  =    T
 ISOLX   =       0 ;ISOLY  =       0 ;ISOLZ  =       0
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  2)
 RHO1    = GRND5
 RHO1A   = 0.000000E+00 ;RHO1B  = 3.484300E-03
 RHO1C   = 7.142900E-01
 PRESS0  = 1.000000E+05
 TEMP0   = 2.730000E+02
 CP1     = 1.004000E+03
 ENUL    = GRND2        ;ENUT   = 0.000000E+00
 ENULA   =-4.946790E-06 ;ENULB  = 4.583940E-08
 ENULC   = 8.097400E-11
 DRH1DP  = GRND5
 DVO1DT  = 3.330000E-03
 PRNDTL(TEM1) = -2.630000E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(BLOK) =  1.000000E+00 ;FIINIT(PRPS) =  2.000000E+00
 FIINIT(TEM1) =  1.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      , GRND2       )
 COVAL (BUOYANCY,V1  , FIXFLU      , GRND2       )

 BUOYA   =-5.300000E+00 ; BUOYB =-8.255000E+00
 BUOYC   = 0.000000E+00
 BUOYD   = 1.231200E+00
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     100
 RESFAC  = 1.000000E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER (TEM1) =   15
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 5.000000E-03)
 RELAX(V1  ,FALSDT, 5.000000E-03)
 RELAX(TEM1,FALSDT, 1.000000E+04)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 NAMGRD  =HTBX
 ASAP    =    T
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(BLOK,N,N,N,N,N,N)
 OUTPUT(PRPS,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       4 ;IYMON  =       5 ;IZMON  =       1
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T
 GVIEW(P,-2.999549E-02,-1.079130E-02,9.994916E-01)
 GVIEW(UP,4.879574E-01,8.725354E-01,2.406454E-02)
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       3.500000E-01, 4.500000E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       1.000000E+00

> OBJ1,   NAME,        BLK
> OBJ1,   POSITION,    4.000000E-01, 4.000000E-01, 0.000000E+00
> OBJ1,   SIZE,        2.000000E-01, 2.000000E-01, 1.000000E+00
> OBJ1,   CLIPART,     cube4
> OBJ1,   ROTATION24,        1
> OBJ1,   TYPE,        BLOCKAGE
> OBJ1,   MATERIAL,      112
> OBJ1,   HEAT_FLVX,     0.000000E+00, 3.000000E+02

> OBJ2,   NAME,        APT1
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 2.000000E-01, 1.000000E+00
> OBJ2,   CLIPART,     cube12
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        APERTURE
> OBJ2,   PRESSURE,      0.000000E+00
> OBJ2,   TEMPERATURE,   1.000000E+01
> OBJ2,   COEFFICIENT,   1.000000E+00

> OBJ3,   NAME,        APT2
> OBJ3,   POSITION,    8.000000E-01, 1.000000E+00, 0.000000E+00
> OBJ3,   SIZE,        2.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ3,   CLIPART,     cube12
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        APERTURE
> OBJ3,   PRESSURE,      0.000000E+00
> OBJ3,   TEMPERATURE,   1.000000E+01
> OBJ3,   COEFFICIENT,   1.000000E+00
M
stop