#$r002
AUTOPLOT USE
file
phi 5
cl
msg Green curve----Temperature at Top Plate
msg Red curve ----Temperature at Lower Plate
da 1 tem1 y 15;da 1 tem1 y 1
col7 1;colf 2
ENDUSE
VRV USE
SCALE 1.000000E+00 1.000000E+01 1.000000E+00
SLICE Z
VIEW 0.000000E+00 0.000000E+00 1.000000E+00
UP 0.000000E+00 1.000000E+00 0.000000E+00
VIEW CENTRE 1.000000E+00 5.000000E-03 5.000000E-01
VIEW SIZE 1.073323E+00; VIEW DEPTH 1.656186E+02
VARIABLE Temperature; CONTOUR ON
MSG Temperature contours
MSG (Y scale increased by 10)
PAUSE
ENDUSE
************************************************************
* GROUP 1. Run identifiers and other preliminaries.
TEXT(LAMINAR FLOW BETWEEN PARALLEL PLATE:R102
TITLE
DISPLAY
This Q1 is setup to simulate the convective-radiative interaction
in a parallel plate channel. Such heat transfer is present in air
operated solar collectors.
The channel considered is 2m long and 0.1m high. The flow
considered is laminar and developing. A constant heat flux of
375W/m2 is applied uniformly to the top plate, whilst the lower
plate is considered to be adiabatic. The entrance and exit of
the channel are treated as apertures at radiative temperatures
of 285K.
This Q1 is setup to enable the PHOENICS predictions to be
compared with the base case results presented by C H Liu
and E M Sparrow in 'Convective-Radiative Interaction in a
Parallel Plate Channel - Application to Air-Operated Solar
Collectors',J Heat Mass Trans,Vol23 pp 1137-1146.
The Q1 contains AUTOPLOT and Viewer macros.
ENDDIS
*
REAL(REYNO,GPRL,GCHNL,GCHNH,GCOND,GCCP)
REAL(GTIN,GQDOT,GUIN,GENUL)
*---------------------------------------------------------------
* REYNO - Reynolds number for flow
* GPRL - Laminar Prandtl number
* GCGHNL - Channel length (m)
* GCHNH - Channel width (m)
* GCOND - Thermal conductivity (W/m/K)
* SIGMA - Stefan-Boltzmann constant (W/m2/K4)
* GCCP - Specific Heat (J/kg/K)
* GTIN - Inlet gas temperature
* GQDOT - Applied heat flux (W/m2)
* GUIN - Mean inflow velocity
* GENUL - Laminar kinematic viscosity (m2/s)
*----------------------------------------------------------------
REYNO=2500.;GPRL=0.7;GCHNL=2.0;GCHNH=GCHNL/200.
GCOND=0.0263;GCCP=1007.
GTIN=0.5*GCOND/(SIGMA*GCHNH);GTIN=GTIN**0.33
GQDOT=SIGMA*GTIN**4.
GENUL=GPRL*GCOND/GCCP;GUIN=REYNO*GENUL/(2.*GCHNH)
* There are in total 42 thermal zones: 2 are apertures;
* 20 are fixed-flux thermal zones at the top plate, and
* 20 are zero-flux zones at the lower plate
INTEGER(NTZ);NTZ=42
* GROUP 3. x-direction grid specification.
* GDX - length of thermal zone
REAL(GDX);GDX=GCHNL/(NTZ/2-1);NREGX=20
IREGX=1;GRDPWR(X,4,GDX,1);IREGX=2;GRDPWR(X,2,GDX,1)
IREGX=3;GRDPWR(X,2,GDX,1);IREGX=4;GRDPWR(X,2,GDX,1)
IREGX=5;GRDPWR(X,2,GDX,1);IREGX=6;GRDPWR(X,2,GDX,1)
IREGX=7;GRDPWR(X,2,GDX,1);IREGX=8;GRDPWR(X,2,GDX,1)
IREGX=9;GRDPWR(X,2,GDX,1);IREGX=10;GRDPWR(X,2,GDX,1)
IREGX=11;GRDPWR(X,2,GDX,1);IREGX=12;GRDPWR(X,2,GDX,1)
IREGX=13;GRDPWR(X,2,GDX,1);IREGX=14;GRDPWR(X,2,GDX,1)
IREGX=15;GRDPWR(X,2,GDX,1);IREGX=16;GRDPWR(X,2,GDX,1)
IREGX=17;GRDPWR(X,2,GDX,1);IREGX=18;GRDPWR(X,2,GDX,1)
IREGX=19;GRDPWR(X,2,GDX,1);IREGX=20;GRDPWR(X,4,GDX,-1.5)
* GROUP 4. y-direction grid specification.
NREGY=1;IREGY=1;GRDPWR(Y,-15,GCHNH,1.5000E+00)
* GROUP 7. Variables (including porosities) named,
* stored & solved.
SOLVE(U1,V1,P1);SOLUTN(P1,Y,Y,Y,N,N,N)
SOLVE(TEM1);SOLUTN(TEM1,Y,Y,Y,N,N,Y)
STORE(DEN1,PRPS,ENUL)
* GROUP 9. Properties of the medium (or media).
RHO1=1.0;ENUL=GENUL; PRNDTL(TEM1)=-GCOND; CP1=GCCP
* GROUP 11. Initialization of fields of variables,
* porosities, etc.
* Initialize Properties Field and block-correction
FIINIT(PRPS)=-1;FIINIT(U1)=GUIN;FIINIT(V1)=0.0;FIINIT(P1)=1.
FIINIT(TEM1)=GTIN
* GROUP 13. Boundary conditions and special sources
INLET
INLET(IN,WEST,#1,#1,#1,#1,#1,#1,#1,#1)
VALUE(IN,P1,GUIN*RHO1);VALUE(IN,U1,GUIN)
VALUE(IN,V1,0.0);VALUE(IN,TEM1,GTIN)
OUTLET
OUTLET(OUT,EAST,#20,#20,#1,#1,#1,#1,#1,#1)
VALUE(OUT,P1,0.0);VALUE(OUT,U1,0.0)
VALUE(OUT,V1,0.0)
* Define the plates as walls for the velocities
WALL(TOP,NORTH,#1,#20,#1,#1,#1,#1,#1,#1)
WALL(BOT,SOUTH,#1,#20,#1,#1,#1,#1,#1,#1)
* RADIATIVE THERMAL ZONES *
* define inlet as an aperture
PATCH(@RI001,WEST,#1,#1,#1,#1,#1,#1,#1,#1)
COVAL(@RI001,TEM1,0.0,GTIN)
* heated wall thermal zones *
PATCH(@RI002,NORTH,#1,#1,#1,#1,#1,#1,#1,#1)
COVAL(@RI002,TEM1,GRND1,GQDOT)
PATCH(@RI003,NORTH,#2,#2,#1,#1,#1,#1,#1,#1)
COVAL(@RI003,TEM1,GRND1,GQDOT)
PATCH(@RI004,NORTH,#3,#3,#1,#1,#1,#1,#1,#1)
COVAL(@RI004,TEM1,GRND1,GQDOT)
PATCH(@RI005,NORTH,#4,#4,#1,#1,#1,#1,#1,#1)
COVAL(@RI005,TEM1,GRND1,GQDOT)
PATCH(@RI006,NORTH,#5,#5,#1,#1,#1,#1,#1,#1)
COVAL(@RI006,TEM1,GRND1,GQDOT)
PATCH(@RI007,NORTH,#6,#6,#1,#1,#1,#1,#1,#1)
COVAL(@RI007,TEM1,GRND1,GQDOT)
PATCH(@RI008,NORTH,#7,#7,#1,#1,#1,#1,#1,#1)
COVAL(@RI008,TEM1,GRND1,GQDOT)
PATCH(@RI009,NORTH,#8,#8,#1,#1,#1,#1,#1,#1)
COVAL(@RI009,TEM1,GRND1,GQDOT)
PATCH(@RI010,NORTH,#9,#9,#1,#1,#1,#1,#1,#1)
COVAL(@RI010,TEM1,GRND1,GQDOT)
PATCH(@RI011,NORTH,#10,#10,#1,#1,#1,#1,#1,#1)
COVAL(@RI011,TEM1,GRND1,GQDOT)
PATCH(@RI012,NORTH,#11,#11,#1,#1,#1,#1,#1,#1)
COVAL(@RI012,TEM1,GRND1,GQDOT)
PATCH(@RI013,NORTH,#12,#12,#1,#1,#1,#1,#1,#1)
COVAL(@RI013,TEM1,GRND1,GQDOT)
PATCH(@RI014,NORTH,#13,#13,#1,#1,#1,#1,#1,#1)
COVAL(@RI014,TEM1,GRND1,GQDOT)
PATCH(@RI015,NORTH,#14,#14,#1,#1,#1,#1,#1,#1)
COVAL(@RI015,TEM1,GRND1,GQDOT)
PATCH(@RI016,NORTH,#15,#15,#1,#1,#1,#1,#1,#1)
COVAL(@RI016,TEM1,GRND1,GQDOT)
PATCH(@RI017,NORTH,#16,#16,#1,#1,#1,#1,#1,#1)
COVAL(@RI017,TEM1,GRND1,GQDOT)
PATCH(@RI018,NORTH,#17,#17,#1,#1,#1,#1,#1,#1)
COVAL(@RI018,TEM1,GRND1,GQDOT)
PATCH(@RI019,NORTH,#18,#18,#1,#1,#1,#1,#1,#1)
COVAL(@RI019,TEM1,GRND1,GQDOT)
PATCH(@RI020,NORTH,#19,#19,#1,#1,#1,#1,#1,#1)
COVAL(@RI020,TEM1,GRND1,GQDOT)
PATCH(@RI021,NORTH,#20,#20,#1,#1,#1,#1,#1,#1)
COVAL(@RI021,TEM1,GRND1,GQDOT)
* adabatic wall *
PATCH(@RI022,SOUTH,#1,#1,#1,#1,#1,#1,#1,#1)
COVAL(@RI022,TEM1,GRND1,0.0)
PATCH(@RI023,SOUTH,#2,#2,#1,#1,#1,#1,#1,#1)
COVAL(@RI023,TEM1,GRND1,0.0)
PATCH(@RI024,SOUTH,#3,#3,#1,#1,#1,#1,#1,#1)
COVAL(@RI024,TEM1,GRND1,0.0)
PATCH(@RI025,SOUTH,#4,#4,#1,#1,#1,#1,#1,#1)
COVAL(@RI025,TEM1,GRND1,0.0)
PATCH(@RI026,SOUTH,#5,#5,#1,#1,#1,#1,#1,#1)
COVAL(@RI026,TEM1,GRND1,0.0)
PATCH(@RI027,SOUTH,#6,#6,#1,#1,#1,#1,#1,#1)
COVAL(@RI027,TEM1,GRND1,0.0)
PATCH(@RI028,SOUTH,#7,#7,#1,#1,#1,#1,#1,#1)
COVAL(@RI028,TEM1,GRND1,0.0)
PATCH(@RI029,SOUTH,#8,#8,#1,#1,#1,#1,#1,#1)
COVAL(@RI029,TEM1,GRND1,0.0)
PATCH(@RI030,SOUTH,#9,#9,#1,#1,#1,#1,#1,#1)
COVAL(@RI030,TEM1,GRND1,0.0)
PATCH(@RI031,SOUTH,#10,#10,#1,#1,#1,#1,#1,#1)
COVAL(@RI031,TEM1,GRND1,0.0)
PATCH(@RI032,SOUTH,#11,#11,#1,#1,#1,#1,#1,#1)
COVAL(@RI032,TEM1,GRND1,0.0)
PATCH(@RI033,SOUTH,#12,#12,#1,#1,#1,#1,#1,#1)
COVAL(@RI033,TEM1,GRND1,0.0)
PATCH(@RI034,SOUTH,#13,#13,#1,#1,#1,#1,#1,#1)
COVAL(@RI034,TEM1,GRND1,0.0)
PATCH(@RI035,SOUTH,#14,#14,#1,#1,#1,#1,#1,#1)
COVAL(@RI035,TEM1,GRND1,0.0)
PATCH(@RI036,SOUTH,#15,#15,#1,#1,#1,#1,#1,#1)
COVAL(@RI036,TEM1,GRND1,0.0)
PATCH(@RI037,SOUTH,#16,#16,#1,#1,#1,#1,#1,#1)
COVAL(@RI037,TEM1,GRND1,0.0)
PATCH(@RI038,SOUTH,#17,#17,#1,#1,#1,#1,#1,#1)
COVAL(@RI038,TEM1,GRND1,0.0)
PATCH(@RI039,SOUTH,#18,#18,#1,#1,#1,#1,#1,#1)
COVAL(@RI039,TEM1,GRND1,0.0)
PATCH(@RI040,SOUTH,#19,#19,#1,#1,#1,#1,#1,#1)
COVAL(@RI040,TEM1,GRND1,0.0)
PATCH(@RI041,SOUTH,#20,#20,#1,#1,#1,#1,#1,#1)
COVAL(@RI041,TEM1,GRND1,0.0)
* outlet radiation patch
PATCH(@RI042,EAST,#20,#20,#1,#1,#1,#1,#1,#1)
COVAL(@RI042,TEM1,0.0,GTIN)
* GROUP 15. Termination criteria for sweeps and
* outer iterations.
LSWEEP=300;TSTSWP=-1
* GROUP 16. Termination criteria for inner iterations.
RESREF(P1)=1.E-12*GUIN*GCHNH*1.0
RESREF(U1)=RESREF(P1)*GUIN;RESREF(V1)=RESREF(U1)
RESREF(TEM1)=RESREF(P1)*GTIN
* GROUP 17. Under-relaxation and related devices.
* Variable declarations
REAL(MAXV,MINL,RELX)
* Estimate of the maximum velocity within domain
MAXV=2
* Estimate of the minimum cell dimension
MINL=2.0000E-03
RELX=100.
RELAX(P1,LINRLX,0.8);RELAX(U1,FALSDT,MINL/MAXV*RELX)
RELAX(V1,FALSDT,MINL/MAXV*RELX)
* GROUP 19. Data communicated by SATELLITE to GROUND
S2SR=t
* GROUP 23. Variable-by-variable field printout and plot
LIBREF = 102
spedat(set,cvd,radcvd,l,t)
STOP