```

talk=f;run(1,1)
DISPLAY

Library case H107:

LINEAR HEAT DISSIPATION

This case is geometricaly identical to case 104.  However,
the problem is modelled  as  transient  and  heat  sources
increase linearly with  time,  from  0  at  the  start  of
calculation to maximum value at the end.

(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
****************************************************
PHOTON USE
p
h1

use patgeo
msg Temperature contours after first time step.
msg
cont tem1 y 6 fil;.001
msg Velocity field after first time step.
vec x 5 sh
vec z m sh
msg
pause
p
h3

use patgeo
msg Temperature contours after third time step.
msg Note increase in the maximum temperature.
msg
cont tem1 y 6 fil;.01
msg Velocity field after third time step.
vec x 5 sh
vec z m sh
msg
pause
p
h5

use patgeo
msg Temperature contours at the end of time period.
msg Note further increase in maximum temperature.
msg
use hotspot
cont tem1 y 6 fil;.001
msg Velocity field at the end of time period.
vec x 5 sh
vec z m sh
msg Type E to end
ENDUSE
************************************************************
Group 1. Run Title
TEXT(Linear Heat Dissipation
************************************************************
Group 2. Transience
* Set overall time and no. of steps
RSET(U,0.000E+00,3.000E+02,5)
* Cut regions
RSET(T,L,0,6.000E+01)
RSET(T,L,0,1.200E+02)
RSET(T,L,0,1.800E+02)
RSET(T,L,0,2.400E+02)
* Modify regions
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,10,10)
* Set overall domain extent:
*        xulast  yvlast  zwlast    name
XSI= 5.000E-01;YSI= 7.000E-01;ZSI= 5.000E-01;RSET(D,BOX     )
* Set objects: x0     y0     z0
*              dx     dy     dz    name
XPO= 1.000E-01;YPO= 7.000E-01;ZPO= 1.500E-01
XSI= 3.000E-01;YSI= 0.000E+00;ZSI= 2.000E-01;RSET(B,VENT    )
XPO= 1.500E-01;YPO= 2.000E-01;ZPO= 5.000E-01
XSI= 2.000E-01;YSI= 3.000E-01;ZSI= 0.000E+00;RSET(B,FAN     )
XPO= 1.000E-01;YPO= 3.500E-01;ZPO= 1.000E-01
XSI= 3.000E-01;YSI= 0.000E+00;ZSI= 3.000E-01;RSET(B,PCB     )
XPO= 1.500E-01;YPO= 3.500E-01;ZPO= 1.500E-01
XSI= 5.000E-02;YSI= 2.000E-02;ZSI= 7.000E-02;RSET(B,CHIP1   )
XPO= 3.000E-01;YPO= 3.500E-01;ZPO= 1.500E-01
XSI= 5.000E-02;YSI= 2.000E-02;ZSI= 7.000E-02;RSET(B,CHIP2   )
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME( 43) =NPOR ; NAME( 44) =BLOK
NAME( 46) =PRPS ; NAME( 48) =DEN1
NAME( 49) =SPH1 ; NAME( 50) =TEM1
* Solved variables list
SOLVE(P1  ,U1  ,V1  ,W1  ,TEM1)
* Stored variables list
STORE(SPH1,DEN1,WDIS,PRPS,BLOK,NPOR)
SOLUTN(P1  ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
IVARBK  =      -1 ;ISOLBK =       1
TURMOD(LVEL)
************************************************************
Group 8. Terms & Devices
NEWRH1  =    T; NEWENL  =    T; NEWENT  =    T
ISOLX   =       0 ;ISOLY  =       0 ;ISOLZ  =       0
************************************************************
Group 9. Properties
RHO1    = FILE
PRESS0  = 1.000E+05; TEMP0   = 2.730E+02
CP1     = FILE; ENUL    = FILE; DRH1DP  = FILE
PRNDTL(LTLS) =  GRND2     ;PRNDTL(TEM1) = CONDFILE
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(NPOR) =  1.000E+00 ;FIINIT(BLOK) =  1.000E+00
FIINIT(PRPS) =  2.000E+00 ;FIINIT(WDIS) =  5.000E-02
FIINIT(TEM1) =  2.000E+01
CONPOR(PCB     ,  0.00,NORTH ,-#2,-#6,-#2,-#2,-#2,-#5)
CONPOR(CHIP1   , -1.00,CELL  ,-#3,-#3,-#3,-#3,-#3,-#3)
INIT(CHIP1   ,BLOK, 0.000E+00, 2.000E+00)
INIT(CHIP1   ,PRPS, 0.000E+00, 1.030E+02)
CONPOR(CHIP2   , -1.00,CELL  ,-#5,-#5,-#3,-#3,-#3,-#3)
INIT(CHIP2   ,BLOK, 0.000E+00, 3.000E+00)
INIT(CHIP2   ,PRPS, 0.000E+00, 1.030E+02)
************************************************************
Group 13. Boundary & Special Sources
INLET (FAN     ,HIGH  ,#3,#5,#2,#4,#6,#6,#1,#5)
VALUE (FAN     ,P1  , 1.783E+00)
VALUE (FAN     ,W1  ,-1.500E+00)
VALUE (FAN     ,TEM1, 2.000E+01)

PATCH (VENT    ,NORTH ,#2,#6,#5,#5,#3,#4,#1,#5)
COVAL (VENT    ,P1  , 1.000E+00, 0.000E+00)
COVAL (VENT    ,TEM1, 0.000E+00, 2.000E+01)

PATCH (TM1CHIP1,VOLUME,#3,#3,#3,#3,#3,#3,#1,#5)
COVAL (TM1CHIP1,TEM1, FIXFLU   , GRND     )

PATCH (TM2CHIP2,VOLUME,#5,#5,#3,#3,#3,#3,#1,#5)
COVAL (TM2CHIP2,TEM1, FIXFLU   , GRND     )

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

BUOYA   = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
BUOYD   = 1.189E+00
EGWF    =    T
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =      20 ;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, 4.667E-01)
RELAX(V1  ,FALSDT, 4.667E-01); RELAX(W1  ,FALSDT, 4.667E-01)
RELAX(LTLS,LINRLX, 1.000E+00); RELAX(TEM1,FALSDT, 4.667E+03)
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD  =HTBX
RG( 41) = 7.143E+04; RG( 42) = 7.143E+04
************************************************************
Group 21. Print-out of Variables
OUTPUT(BLOK,N,N,N,N,N,N); OUTPUT(PRPS,N,N,N,N,N,N)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON   =       5 ;IYMON  =       7 ;IZMON  =       4
TSTSWP  =      -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN  =       1; NYPRIN  =       1; NZPRIN  =       1
NPLT    =       1; ITABL   =       3
************************************************************
Group 24. Dumps For Restarts
IDISPA  =       1; CSG1    =H
STOP

```