photon use
p
gr z 1
msg velocity vectors
vec z 1 sh
pause
msg pressure contours
con p1 z 1 fi;0.001
pause
con off
red
msg temperature contours
con temp z 1 fi;0.001
vec z 1
msg
msg Press e to END
enduse
DISPLAY
This case tests the "link" procedure in the x-y plane.
The geometry and boundary conditions are as shown:
----------------------------------
| | |
| outflow |link |
<-------- Isection I
I I I
I barrier I open I
I--------------------------------I
I I I
I inflow I I
-------> I I
I I I
I I I
I----------------I---------------I
ENDDIS
delay(500)
GROUP 1. Run title and other preliminaries
TEXT(East-West Link Test In X-Y Plane
TITLE
mesg(PC486/50 time last reported as 1.5 min
REAL(XLENGTH,YLENGTH,ZLENGTH)
INTEGER(NXNOM,IYSHFT,ISHIFT,IYPLUS,IYBOT)
XLENGTH=1.0;YLENGTH=1.0;ZLENGTH=1.0
NXNOM=10;NY=10;NZ=1
NPHI=20
NX=NXNOM+2
ISHIFT=2*NY-IYSHFT
GROUP 3. X-direction grid specification
**Domain is XLENGTH m long in x-direction, with equal intervals
XULAST=XLENGTH
XFRAC(1)=-NXNOM/2
XFRAC(2)=1/NXNOM
XFRAC(5)=-XFRAC(1)
XFRAC(6)=XFRAC(2)
dummy link cells are made thin for print-out purposes
XFRAC(3)=2
XFRAC(4)=0.01*XFRAC(2)
GROUP 4. Y-direction grid specification
**Domain is YLENGTH m long in y-direction, with equal intervals
GRDPWR(Y,NY,YLENGTH,1.0)
GROUP 5. Z-direction grid specification
**Domain is ZLENGTH m long in z-direction, with equal intervals
GRDPWR(Z,NZ,ZLENGTH,1.0)
GROUP 7. Variables stored, solved & named
**Choose first-phase enthalpy (H1) as dependent variable
and activate the whole-field elliptic solver
SOLUTN(H1,Y,Y,Y,N,N,N);NAME(H1)=TEMP
SOLVE(P1,U1,V1);SOLUTN(P1,Y,Y,Y,N,N,N)
GROUP 8. Terms (in differential equations) & devices
**For pure conduction, cut out built-in source and convection
terms
TERMS(TEMP,N,N,Y,N,Y,Y)
GROUP 9. Properties of the medium (or media)
**Thermal conductivity will be ENUL*RHO1/PRNDTL(TEMP), so :
ENUL=1.0e-3;PRNDTL(TEMP)=1.0
GROUP 11. Initialization of variable or porosity fields
INIADD=F
FIINIT(U1)=1.0
PATCH(INIT,INIVAL,1,NX,NY/2+1,NY,1,1,1,1)
COVAL(INIT,U1,0.0,-1.0)
GROUP 13. Boundary conditions and special sources
** Cold inflow boundary on the bottom left
PATCH(COLD,WEST,1,1,1,NY/2,1,1,1,1)
COVAL(COLD,U1,ONLYMS,1.0)
COVAL(COLD,TEMP,FIXVAL,-0.9)
COVAL(COLD,P1,FIXFLU,1.0)
** hot outflow boundary on the top left
PATCH(HOT,CELL,1,1,NY/2+1,NY,NZ,NZ,1,1)
COVAL(HOT,TEMP,FIXVAL,0.9)
COVAL(HOT,P1,1.E-2,0.0)
** the link patches
PATCH(+1,EAST,NXNOM/2+1,NXNOM/2+1,1,NY,1,1,1,1)
COVAL(+1,TEMP,fixval,ISHIFT)
COVAL(+1,P1,FIXVAL,ISHIFT)
COVAL(+1,U1,FIXVAL,ISHIFT)
COVAL(+1,V1,FIXVAL,ISHIFT)
PATCH(+2, WEST,NXNOM/2+2, NXNOM/2+2, 1,NY,1,1,1,1)
COVAL(+2,TEMP,FIXVAL,-ISHIFT)
COVAL(+2,P1,FIXVAL,-ISHIFT)
COVAL(+2,V1,FIXVAL,-ISHIFT)
** halving the diiffusion coefficients because the dummy
cells are very small
PATCH(GP12DFE1,EAST,NXNOM/2,NXNOM/2,1,NY,1,1,1,1)
COVAL(GP12DFE1,TEMP,0.5,0.0)
PATCH(GP12DFE2,EAST,NXNOM/2+2,NXNOM/2+2,1,NY,1,1,1,1)
COVAL(GP12DFE2,TEMP,0.5,0.0)
** the barrier between the top and bottom left-hand halves
PATCH(FENCE,NORTH,1,NX/2,NY/2,NY/2,1,1,1,1)
COVAL(FENCE,V1,FIXVAL,0.0)
PATCH(GP12DFN1,NORTH,1,nxnom/2,NY/2,NY/2,1,1,1,1)
COVAL(GP12DFN1,TEMP,0.0,0.0)
GROUP 15. Termination of sweeps
LSWEEP=20
GROUP 16. Termination of iterations
** Set the frequencies of application of the one-dimensional
correction features in the linear-equation solver to once
per iteration for each direction.
ISOLX=1;ISOLY=1;ISOLZ=1
LITER(TEMP)=10;LITER(P1)=20
GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.5);RELAX(V1,FALSDT,0.01);RELAX(U1,FALSDT,0.01)
GROUP 21. Print-out of variables
**Print fields of temperature
OUTPUT(TEMP,Y,Y,Y,Y,Y,Y)
GROUP 22. Spot-value print-out
IYMON=NY/2+1;IZMON=NZ/2+1;IXMON=NX-1;ITABL=1
GROUP 23. Field print-out and plot control
IXPRF=NXNOM/2;IXPRL=NXNOM/2+4
GROUP 24. Dumps for restarts
NPLT=1;TSTSWP=-1
LSWEEP=200
SELREF=T;RESFAC=1.E-5
IF(NY.EQ.1) THEN
SOLUTN(V1,N,N,N,N,N,N);OUTPUT(V1,N,N,N,N,N,N)
ENDIF
lsg57=t