PHOTON USE
p
GR OU X 21 Y 1 3 Z 1 1
GR OU X 1 Y 1 3 Z 1 1
GR OU Y 4 X 1 20 Z 1 1
GR OU Y 1 X 1 20 Z 1 1
GR OU X 6 Y 4 10 Z 1 1
GR OU X 1 Y 4 10 Z 1 1
GR OU Y 11 X 1 5 Z 1 1
GR OU Y 4 X 1 5 Z 1 1
GR OU X 4 Y 11 14 Z 1 1
GR OU X 1 Y 11 14 Z 1 1
GR OU Y 15 X 1 3 Z 1 1
GR OU Y 11 X 1 3 Z 1 1
GR OU X 6 Y 15 22 Z 1 1
GR OU X 1 Y 15 22 Z 1 1
GR OU Y 23 X 1 5 Z 1 1
GR OU Y 15 X 1 5 Z 1 1
GR OU X 21 Y 15 22 Z 1 1
GR OU X 14 Y 15 22 Z 1 1
GR OU Y 23 X 14 20 Z 1 1
GR OU Y 15 X 14 20 Z 1 1
GR OU X 21 Y 11 14 Z 1 1
GR OU X 16 Y 11 14 Z 1 1
GR OU Y 15 X 16 20 Z 1 1
GR OU Y 11 X 16 20 Z 1 1
GR OU X 21 Y 9 10 Z 1 1
GR OU X 14 Y 9 10 Z 1 1
GR OU Y 11 X 14 20 Z 1 1
GR OU Y 9 X 14 20 Z 1 1
GR OU X 16 Y 11 14 Z 1 1 COL 14
GR OU X 4 Y 11 14 Z 1 1 COL 14
GR OU Y 15 X 4 15 Z 1 1 COL 14
GR OU Y 11 X 4 15 Z 1 1 COL 14
MSG COOLER GEOMETRY
msg Press Enter to continue
PAUSE;red
VEC Z 1 Y 1 28 SH
MSG Velocity vectors
msg Press Enter to continue
pause;vec off;red
con tmp1 z 1 Y 1 28 fil;0.01
msg Temperature contours
msg Press Enter to continue
pause;con off;red
con metl z 1 Y 1 28 fil;0.01
msg Metal temperature
msg Press Enter to continue
pause;con off;red
con p1 z 1 Y 1 28 fil;0.01
msg Pressure contours
msg Press Enter to continue
pause;con off;red
con ke z 1 Y 1 28 fil;0.01
msg Turbulent energy contours
msg Press Enter to continue
pause;con off;red
con enut z 1 Y 1 28 fil;0.01
msg Turbulent viscosity contours
msg Press Enter to end
ENDUSE
************************************************************
* GROUP 1. Run identifiers and other preliminaries.
TEXT(Natural cooler in a still atmosphere:134
DISPLAY
The packed bed of heat generating materials is cooled by air
flow generated by buoyancy. The pressure drop across a packed
bed is obtained from Ergun equation. The material temperatures
are also calculated. PLANT is used to introduce the
corresponding sources.
ENDDIS
==================
REAL(TENV,PORI,QDVOL,VVOL); TENV=-40.; PORI=0.944
VVOL=5.2*4.62*2.345
QDVOL=20.E6/VVOL*(1-PORI)
REAL(DIAM,LENGTH);DIAM=0.02
LENGTH=DIAM*PORI/(1.-PORI)
************************************************************
* GROUP 2. Time-dependence and related parameters.
************************************************************
* GROUP 3. x-direction grid specification.
NX = 25
XULAST = 14.4
XFRAC ( 1) = 4.630E-02 ;XFRAC ( 2) = 9.259E-02
XFRAC ( 3) = 1.389E-01 ;XFRAC ( 4) = 1.649E-01
XFRAC ( 5) = 1.910E-01 ;XFRAC ( 6) = 2.179E-01
XFRAC ( 7) = 2.448E-01 ;XFRAC ( 8) = 2.717E-01
XFRAC ( 9) = 2.986E-01 ;XFRAC ( 10) = 3.255E-01
XFRAC ( 11) = 3.524E-01 ;XFRAC ( 12) = 3.793E-01
XFRAC ( 13) = 4.063E-01 ;XFRAC ( 14) = 4.323E-01
XFRAC ( 15) = 4.583E-01 ;XFRAC ( 16) = 4.972E-01
XFRAC ( 17) = 5.361E-01 ;XFRAC ( 18) = 5.750E-01
XFRAC ( 19) = 6.139E-01 ;XFRAC ( 20) = 6.528E-01
XFRAC ( 21) = 6.838E-01 ;XFRAC ( 22) = 7.406E-01
XFRAC ( 23) = 8.142E-01 ;XFRAC ( 24) = 9.012E-01
XFRAC ( 25) = 1.000E+00
************************************************************
* GROUP 4. y-direction grid specification.
NY = 30
YVLAST = 27.0
YFRAC ( 1) = 6.173E-02 ;YFRAC ( 2) = 1.235E-01
YFRAC ( 3) = 1.852E-01 ;YFRAC ( 4) = 2.022E-01
YFRAC ( 5) = 2.193E-01 ;YFRAC ( 6) = 2.363E-01
YFRAC ( 7) = 2.533E-01 ;YFRAC ( 8) = 2.704E-01
YFRAC ( 9) = 2.778E-01 ;YFRAC ( 10) = 2.852E-01
YFRAC ( 11) = 3.065E-01 ;YFRAC ( 12) = 3.278E-01
YFRAC ( 13) = 3.491E-01 ;YFRAC ( 14) = 3.704E-01
YFRAC ( 15) = 4.028E-01 ;YFRAC ( 16) = 4.352E-01
YFRAC ( 17) = 4.676E-01 ;YFRAC ( 18) = 5.000E-01
YFRAC ( 19) = 5.324E-01 ;YFRAC ( 20) = 5.648E-01
YFRAC ( 21) = 5.972E-01 ;YFRAC ( 22) = 6.296E-01
YFRAC ( 23) = 6.759E-01 ;YFRAC ( 24) = 7.222E-01
YFRAC ( 25) = 7.685E-01 ;YFRAC ( 26) = 8.148E-01
YFRAC ( 27) = 8.611E-01 ;YFRAC ( 28) = 9.074E-01
YFRAC ( 29) = 9.537E-01 ;YFRAC ( 30) = 1.000E+00
* GROUP 7. Variables (including porosities) named,
* stored & solved.
SOLVE(P1,U1,V1,H1,METL); SOLUTN(P1,Y,Y,Y,N,N,N)
STORE(TMP1)
STORE(ENUT)
************************************************************
* GROUP 8. Terms (in differential equations) and devices.
TERMS(METL,N,N,N,N,Y,N)
************************************************************
* GROUP 9. Properties of the medium (or media).
PRESS0=1.0000E+05
* CP1 = Specific heat (J/kg K)
REAL(CAP);CAP=996.
* TREF1 = Reference temperature (K)
REAL(TREF1);TREF1=273.
* HREF1 = Reference enthalpy (J/kg)
REAL(HREF1);HREF1=CAP*(TREF1+TENV)
* VEXP1 = Volume expansivity coefficient (1/K)
REAL(VEXP1);VEXP1=3.6720E-03
* Density (kg/m^3)
RHO1=1.2200E+00
* Set TMP1 = TREF + H1/CP1
TMP1=LINH;TMP1A=0.0;TMP1B=1./CAP;CP1=CAP
* Prandtl number of the fluid/media
PRNDTL(H1)=7.3600E-01
* Laminar kinematic viscosity (m^2/s)
ENUL=1.4650E-05
TURMOD(KEMODL)
* Turbulent Prandtl number of the fluid/media
PRT(H1)=1
************************************************************
* GROUP 10. Interphase-transfer processes and properties.
************************************************************
* GROUP 11. Initialization of fields of variables,
* porosities, etc.
CONPOR(0,CELL,1,-20,1,-3,1,1)
CONPOR(0,CELL,1,-5,4,10,1,1)
CONPOR(0,CELL,1,-3,11,14,1,1)
CONPOR(0,CELL,1,-5,15,-22,1,1)
CONPOR(0,CELL,-14,-20,15,-22,1,1)
CONPOR(0,CELL,-16,-20,11,14,1,1)
CONPOR(0,CELL,-14,-20,-9,10,1,1)
CONPOR(PORI,VOLUME,4,15,11,14,1,1)
CONPOR(PORI,EAST,4,15,11,14,1,1)
CONPOR(PORI,NORTH,4,15,11,14,1,1)
FIINIT(U1)=0.0
FIINIT(V1)=0.1
FIINIT(H1)=HREF1
FIINIT(KE)=5.0E-02
FIINIT(EP)=2.4573E-01
*
************************************************************
* GROUP 12. Convection and diffusion adjustments
************************************************************
* GROUP 13. Boundary conditions and special sources
OUTLET(OUTR,EAST,NX,NX,1,NY,1,NZ,1,1)
VALUE(OUTR,U1,SAME)
VALUE(OUTR,V1,SAME)
VALUE(OUTR,P1,0.)
VALUE(OUTR,H1,CAP*(TREF1+TENV))
VALUE(OUTR,KE,SAME)
VALUE(OUTR,EP,SAME)
OUTLET(OUTL,WEST,1,1,23,30,1,NZ,1,1)
VALUE(OUTL,U1,SAME)
VALUE(OUTL,V1,SAME)
VALUE(OUTL,P1,0.)
VALUE(OUTL,H1,CAP*(TREF1+TENV))
VALUE(OUTL,KE,SAME)
VALUE(OUTL,EP,SAME)
* OUTLET boundary condition, name OUTUP
OUTLET(OUTUP,NORTH,1,NX,NY,NY,1,NZ,1,1)
VALUE(OUTUP,U1,SAME)
VALUE(OUTUP,V1,SAME)
VALUE(OUTUP,P1,0.)
VALUE(OUTUP,H1,CAP*(TREF1+TENV))
VALUE(OUTUP,KE,SAME)
VALUE(OUTUP,EP,SAME)
* HEAT-SOURCE boundary condition, name HOTHEAT
PATCH(HEAT,VOLUME,4,15,11,14,1,1,1,1)
COVAL(HEAT,H1,FIXFLU,QDVOL)
COVAL(HEAT,METL,FIXFLU,QDVOL)
** Ergun's resitance law
* Linear resistance term
PATCH(ERLIN,PHASEM,4,15,11,14,1,1,1,1)
COVAL(ERLIN,U1,0.1*ENUL/LENGTH**2,0.0)
COVAL(ERLIN,V1,0.1*ENUL/LENGTH**2,0.0)
PLANTBEGIN
* Quadratic resistance term
PATCH(ERQUA,PHASEM,4,15,11,14,1,1,1,1)
CO=RG(2)*ABS(U1)
COVAL(ERQUA,U1,GRND,0.0)
CO=RG(2)*ABS(V1)
COVAL(ERQUA,V1,GRND,0.0)
* Metal temperature
PATCH(METEMP,VOLUME,4,15,11,14,1,1,1,1)
VAL=TMP1
COVAL(METEMP,METL,100.,GRND)
PLANTEND
* GRAVITY boundary condition, name BUOYANCY
PATCH(BUOYANCY,PHASEM,1,NX,1,NY,1,1,1,1)
COVAL(BUOYANCY,U1,FIXFLU,BOUSS)
COVAL(BUOYANCY,V1,FIXFLU,BOUSS)
* Set gravity resolutes
BUOYA=0.;BUOYB=-9.8100E+00;BUOYC=0.
* Set contants for Boussinesq approximation
BUOYD=VEXP1/CAP;BUOYE=HREF1
* GROUP 15. Termination criteria for sweeps and
LSWEEP=100
* GROUP 17. Under-relaxation and related devices.
REAL(MAXV,MINL,RELX)
MAXV=4.
MINL=5.7600E-01
RELX=1
RELAX(P1,LINRLX,0.8)
RELAX(U1,FALSDT,MINL/MAXV*RELX)
RELAX(V1,FALSDT,MINL/MAXV*RELX)
RELAX(H1,FALSDT,MINL/MAXV*1.E3)
RELAX(KE,FALSDT,MINL/MAXV*RELX/10)
RELAX(EP,FALSDT,MINL/MAXV*RELX/10)
* GROUP 19. Data communicated by SATELLITE to GROUND
NAMSAT=MOSG
RG(1)=150.*ENUL/LENGTH**2
RG(2)=0.0175/LENGTH
* GROUP 22. Location of spot-value & frequency of
* residual printout.
IXMON=14; IYMON=13
* GROUP 23. Variable-by-variable field printout and plot
ITABL=3
tstswp=-1
dmpstk=t
DISTIL=T
STORE(PRPS,NPOR,EPOR)
EX(P1)=1.865E+00; EX(U1)=4.052E-01; EX(V1)=7.230E-01
EX(KE)=7.989E-03; EX(EP)=4.578E-03; EX(H1)=1.641E+05
EX(PRPS)=6.325E+01; EX(NPOR)=9.964E-01; EX(EPOR)=9.964E-01
EX(VPOR)=6.804E-01; EX(ENUT)=1.459E-03
EX(TMP1)=1.641E+05; EX(METL)=1.665E+04
LIBREF=134
STOP