PHOTON USE
p
parphi
1 20 1
up z
vec x 1 sh
pause
con wcrt x 1 fi;0.001
pause
con h1 x 1 fi;0.001
ENDUSE
GROUP 1. Run title and other preliminaries
TEXT(Free Convec From Hot Vertical Wall
TITLE
DISPLAY
This run is similar to the run above,except that the flow is
driven by a buoyancy force rather than by a prescribed mass
inflow. As a consequence,the grid expands at a different rate.
| ^ ^ |/
^ ^ ^ | | ^ |/
Pressure, | ^ ^ ^ ^ | | | | | | ^ |/
velocity |/
and | |/
enthalpy |/
fixed to | | |/
zero |gravity |/
| v |/ Wall;
|/ heated with H1=1.0
| |/
|/
^ | |/
| z |/
<--| | ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ |/
y
ENDDIS
Interesting variants include changes to: the magnitude of the
added force; and to the wall boundary values.
GXBUOY is used for this case.
GROUP 4. Y-direction grid specification
GRDPWR(Y,20,4.0E-3,1.0)
*** Grid grows as (axial distance)**0.25
AZYV=0.25;ZWADD=1.0E-03
GROUP 5. Z-direction grid specification
*** The flow is parabolic
PARAB=T
*** Axial grid distance varies as IZ**4.0
GRDPWR(Z,100,0.1,4.0)
AZDZ=GRND
(ZGRID OF DZ IS 0.5*YVLAST)
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1,H1)
GROUP 8. Terms (in differential equations) & devices
*** Deactivate the built-in source term in H1 equation
TERMS(H1,N,Y,Y,Y,Y,Y)
GROUP 9. Properties of the medium (or media)
ENUL=1.E-5; PRNDTL(H1)=0.7
CP1 = 1.005000E+03
GROUP 11. Initialization of variable or porosity fields
FIINIT(V1)=0.0; FIINIT(W1)=0.0; FIINIT(H1)=0.0
GROUP 13. Boundary conditions and special sources
1. Wall boundary: constant temperature
WALL (CONST,SOUTH,1,1,1,1,1,NZ,1,1)
COVAL(CONST,W1,1.0,0.0)
COVAL(CONST,H1,1.0,1.0)
2. Free stream boundary
PATCH(HIGHY,NORTH,1,1,NY,NY,1,NZ,1,1)
COVAL(HIGHY,P1,FIXP,0.0);COVAL(HIGHY,W1,ONLYMS,0.0)
COVAL(HIGHY,H1,ONLYMS,0.0)
3. Buoyancy source
*** The added source in the Z-direction momentum equation
is RHO*Volume*grav.*dvo1dt/cp*H1 (See GREX3, Group 13,Sec.15,
and GXBUOY);
the following data need to be set for this purpose:
BUOYE=H1ref; BUOYC=grav.
BUOYE=0.0; BUOYC=-9.81
DVO1DT=1.005E+02
PATCH(BUOYANCY,PHASEM,1,1,1,NY,1,NZ,1,1)
COVAL(BUOYANCY,W1,FIXFLU,BOUSS)
GROUP 14. Downstream pressure for PARAB=T
*** The downstream pressure is set by the parameter PBAR
with its default value of 0.0.
IPARAB=1
GROUP 16. Termination of iterations
LITHYD=10
GROUP 17. Under-relaxation devices
RELAX(W1,FALSDT,0.1); RELAX(V1,FALSDT,0.1)
GROUP 22. Spot-value print-out
NPRMON=NZ/10; IYMON=3; NPLT=1; TSTSWP=2
GROUP 23. Field print-out and plot control
NZPRIN=NZ/4
PATCH(FIXEDZ,PROFIL,1,1,1,NY,1,NZ,1,1)
PLOT (FIXEDZ,W1,0.0,0.0); PLOT(FIXEDZ,H1,0.0,0.0)
IDISPA=1;IDISPB=10;IDISPC=NZ
#CONPROM
#MAXABS
#ENDPAUSE
TSTSWP=-1