Library case H101:

THE EFFECT OF ENCLOSURE ORIENTATION



talk=f;run(1,1)
  
#cls
  DISPLAY

  Library case H101:

  
  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.

  (Use PHOTON to display  results.  Enter  USE Q1  in PHOTON
  command mode,  or create a file called U,  which  contains
  the one line: USE Q1   before typing  RUNPHO.)
  ENDDIS
#pause
  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
 ************************************************************
  Group 1. Run Title
 TEXT(Tilted Box                       
 ************************************************************
  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.000E+00;YSI= 1.000E+00;ZSI= 1.000E+00;RSET(D,BOX     )
    * Set objects: x0     y0     z0
    *              dx     dy     dz    name
 XPO= 4.000E-01;YPO= 4.000E-01;ZPO= 0.000E+00
 XSI= 2.000E-01;YSI= 2.000E-01;ZSI= 1.000E+00;RSET(B,BLK     )
 XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00
 XSI= 0.000E+00;YSI= 2.000E-01;ZSI= 1.000E+00;RSET(B,APT1    )
 XPO= 8.000E-01;YPO= 1.000E+00;ZPO= 0.000E+00
 XSI= 2.000E-01;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,APT2    )
 ************************************************************
  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
 RHO1    = FILE
 PRESS0  = 1.000E+05; TEMP0   = 2.730E+02
 CP1     = FILE
 ENUL    = FILE    ;ENUT   = 0.000E+00
 PRNDTL(TEM1) = CONDFILE
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(BLOK) =  1.000E+00 ;FIINIT(PRPS) =  2.000E+00
 FIINIT(TEM1) =  1.000E+01
 CONPOR(BLK     , -1.00,CELL  ,-#2,-#2,-#3,-#3,-#1,-#1)
 INIT(BLK     ,BLOK, 0.000E+00, 2.000E+00)
 INIT(BLK     ,PRPS, 0.000E+00, 1.120E+02)
 INIADD  =    F
 ************************************************************
  Group 13. Boundary & Special Sources
 PATCH (APT1    ,WEST  ,#1,#1,#1,#1,#1,#1,1,1)
 COVAL (APT1    ,P1  , 1.000E+00, 0.000E+00)
 COVAL (APT1    ,TEM1, 0.000E+00, 1.000E+01)

 PATCH (APT2    ,NORTH ,#4,#4,#4,#4,#1,#1,1,1)
 COVAL (APT2    ,P1  , 1.000E+00, 0.000E+00)
 COVAL (APT2    ,TEM1, 0.000E+00, 1.000E+01)

 PATCH (HOTBLK  ,VOLUME,#2,#2,#3,#3,#1,#1,#1,#1)
 COVAL (HOTBLK  ,TEM1, FIXFLU   , 3.000E+02)

 PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT)
 COVAL (BUOYANCY,U1  , FIXFLU   , DENSDIFF    )
 COVAL (BUOYANCY,V1  , FIXFLU   , DENSDIFF    )

 BUOYA   =-5.300E+00 ; BUOYB =-8.255E+00 ; BUOYC = 0.000E+00
 BUOYD   = 1.2312E+00
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     100; SELREF  =    T; RESFAC  = 1.000E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER (TEM1) =   15
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000E+00); RELAX(U1  ,FALSDT, 5.000E-03)
 RELAX(V1  ,FALSDT, 5.000E-03); RELAX(TEM1,FALSDT, 1.000E+04)
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 NAMGRD  =HTBX
 ************************************************************
  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
 TSTSWP=-1
 ************************************************************
  Group 24. Dumps For Restarts
STOP