PHOTON USE
p
red
msg u1 contours. Press return for w1
con u1 z 1 fi;0.001;con u1 z 5 fi;0.001;con u1 z 10 fi;0.001
gr ou z m;gr ou y m;gr ou x m
pause
con off;red
msg w1 contours, and velocity vectors on the top.
msg Press return for h1
con w1 y 1 fi;0.001;con w1 y 5 fi;0.001;con w1 y m fi;0.001
vec y m
pause
vec off;con off;red
msg temperature contours. Press e to end
con h1 z 1 fi;0.001;con h1 z 5 fi;0.001;con h1 z 10 fi;0.001
enduse
xprtbegin ... start of the inputs to expert via Q1
expsol f .... expert is not used in the solver
expdtf t .... expert is used to adjust false time step
indtf 2 ... adjust the false time step of all velocities
inres 3 ... monitor the residuals of variable 3
ifrsts 10 ... first adjustment sweep
ifrequ 10 ... adjustments will be made every ifrequ sweeps
facdec 0.5 .. decreases will be by a factor of facdec
facinc 2.0 .. increases will be by a factor of facinc
exprin t .... print adjusted false time steps in RESULT file
expend ..... end of the inputs to expert via Q1
xprtend
EXPERT=t
GROUP 1. Run title and other preliminaries
TEXT(Cubical Cavity With Moving Lid
TITLE
DISPLAY
This simulation is similar to that of case 249; but the
cavity is limited in the third direction by adiabatic
end walls, at rest. Economy could be affected by noting
that symmetry must prevail about the NZ/2 plane; but
this is not done here. The symmetry of the solution
can therefore be checked..
_ _ _ _ _ _ _
/ /|
/ north wall / |
/ at H1=1. --/--|--->moving
/_ _ _ _ _ _ _/ | y
| | | ^ ^
| |east wall | /z
| low wall | /H1=0. | /
| adiabatic | / |/----> x
|_ _ _ _ _ _ _|/
EXPERT is in use
The Q1 contains PHOTON USE commands
ENDDIS
GROUP 3. X-direction grid specification
** Set a symmetrical grid, consisting of power-law
grids (varying with IX**2.0), which start from
each edge and meet in the middle.
GRDPWR(X,-10,0.1,2.0)
GROUP 4. Y-direction grid specification
** Set a symmetrical grid as in GROUP 3.
GRDPWR(Y,-10,0.1,2.0)
GROUP 5. Z-direction grid specification
** Set a symmetrical grid as in GROUP 3.
GRDPWR(Z,-10,0.1,2.0)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1,H1)
**Select whole-field solution procedure for pressure,
SOLUTN(P1,Y,Y,Y,N,N,N)
GROUP 8. Terms (in differential equations) & devices
** De-activate 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
GROUP 13. Boundary conditions and special sources
** South wall, at rest
WALL (SOUTH,SOUTH,1,NX,1,1,1,NZ,1,1)
COVAL(SOUTH,U1,1.0,0.0);COVAL(SOUTH,H1,1.0,0.0)
COVAL(SOUTH,W1,1.0,0.0)
** North wall, moving
WALL (MOVING,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(MOVING,U1,1.0,1.0);COVAL(MOVING,H1,1.0,1.0)
COVAL(MOVING,W1,1.0,0.0)
** West wall, at rest
WALL (WEST,WEST,1,1,1,NY,1,NZ,1,1)
COVAL(WEST,H1,1.0,0.0);COVAL(WEST,V1,1.0,0.0)
COVAL(WEST,W1,1.0,0.0)
** East wall, at rest
WALL (EAST,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(EAST,V1,1.0,0.0);COVAL(EAST,H1,1.0/PRNDTL(H1),0.0)
COVAL(EAST,W1,1.0,0.0)
** Adiabatic high wall, at rest
WALL (HIGH,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(HIGH,U1,1.0,0.0);COVAL(HIGH,V1,1.0,0.0)
** Adiabatic low wall, at rest
WALL (LOW,LOW,1,NX,1,NY,1,1,1,1)
COVAL(LOW,U1,1.0,0.0);COVAL(LOW,V1,1.0,0.0)
** Pressure relief
PATCH(RELIEF,CELL,NX/2,NX/2,NY/2,NY/2,1,1,1,1)
COVAL(RELIEF,P1,FIXP,0.0);COVAL(RELIEF,H1,ONLYMS,SAME)
GROUP 15. Termination of sweeps
LSWEEP=200;ISWC1=LSWEEP/2
SELREF=T;RESFAC=1.0
GROUP 16. Termination of iterations
LITER(P1)=100;CSG3=CNGR
LITER(U1)=20;LITER(V1)=20;LITER(W1)=20
GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,0.001);RELAX(V1,FALSDT,0.001)
RELAX(W1,FALSDT,0.001)
GROUP 21. Print-out of variables
GROUP 22. Spot-value print-out
IXMON=NX/2;IYMON=NY-1;IZMON=NZ/2
TSTSWP=-1
GROUP 23. Field print-out and plot control
XZPR=T
PATCH(MIDX,PROFIL,NX/2,NX/2,1,NY,NZ/2,NZ/2,1,1)
PLOT(MIDX,U1,0.0,0.0);PLOT(MIDX,H1,0.0,0.0)
PATCH(MAP,CONTUR,1,NX,1,NY,NZ/2,NZ/2,1,1)
PLOT(MAP,U1,0.0,10);PLOT(MAP,V1,0.0,10)
PLOT(MAP,H1,0.0,10);PLOT(MAP,W1,0.0,10)
#conprom
#maxabs
#endpause