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
msg Press return to load results after third time step.
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
msg Press return to load results after the last time step.
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
STEADY=F
* 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)
* Additional solver options
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)
INIADD = F
************************************************************
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