PHOTON USE
p
p1;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 1 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p2;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 2 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p3;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 3 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p4;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 4 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p5;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 5 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p6;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 6 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p7;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 7 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p8;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 8 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p9;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 9 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p10;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 10 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p11;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 11 sec.
msg Velocity vectors and concentration field
msg Hit Enter to continue
pause
p
p12;
con C1 z 1 fil;.001
surf mark z .99
ve z 1 sh
msg Time = 12 sec.
msg Velocity vectors and concentration field
msg Hit Enter to finish
ENDUSE
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
DISPLAY
The geometry of this case is a rectangular vessel with
two-padlle impellor inside. There are therefore 2
velocity components.
The vessel is filled with the liquid initial state of
which is characterized by difference in concentration of
inert contaminant: the upper part of the vessel is pure
liquid while lower part concentartion is equal to unity.
The paddle is supposed to be suddenly set in rotation.
The computational task is to predict the mixing,
represented by the subsequent distributions of velocity,
pressure and time-average concentration.
PHOTON use commands are supplied.
ENDDIS
PLANT information :
* Data input groups used: 13, 19
* Ground groups planted : 13, 19-3, 19-7
* Headings used : SORC??, SC03??, SC07??
* Functions used : SPHERE, XYELLP
* Commands used : IF, REGION
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
GROUP 1. Run title and other preliminaries
TEXT( Unsteady mixing in two paddle-stirred reactor
** Impellor cells marker = 198
REAL(mrk);mrk=198.;RG(1)=mrk
REAL(PI)
+ PI =3.14159
** Number of revolutions, 1/s.
+ RG(2)=PI/6.
** Impellor is represented as two paddles of virtually
ellipsoidal shapes fixed at the cylindrical shaft.
* Paddle half axis,
in X direction in Y direction
REAL(AA , BB )
+ AA =8.0 ; BB =1.0
+ RG(3)=AA ; RG(4)=BB
* Rotation axis coordinates,
X coordinate Y coordinate
REAL(XIC , YIC )
+ XIC =10.0 ; YIC =10.0
+ RG(5)=XIC ; RG(6)=YIC
GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,12,12.,1.0)
GROUP 3. X-direction grid specification
XULAST=20.
GRDPWR(X,40,XULAST,1.0)
GROUP 4. Y-direction grid specification
YVLAST=20.
GRDPWR(Y,40,YVLAST,1.)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,C1)
STORE(CONC,MARK)
GROUP 9. Properties of the medium (or media)
RHO1=1000.; ENUL=0.01
GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0;FIINIT(U1)=0.0
FIINIT(MARK)=0.0
** Concentration initialization
PATCH(WC1,INIVAL,1,NX,1,NY/2,1,NZ,1,LSTEP)
INIT(WC1,C1,0.0,0.0)
PATCH(EC1,INIVAL,1,NX,1,NY/2+1,1,NZ,1,LSTEP)
INIT(EC1,C1,0.0,1.0)
GROUP 13. Boundary conditions and special sources
** North wall
PATCH(WALLN,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
COVAL(WALLN,U1,1.,0.0)
** South wall
PATCH(WALLS,SWALL,1,NX,1,1,1,NZ,1,LSTEP)
COVAL(WALLS,U1,1.,0.0)
** East wall
PATCH(WALLE,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
COVAL(WALLE,V1,1.,0.0)
** West wall
PATCH(WALLW,WWALL,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(WALLW,V1,1.,0.0)
** Pressure relief
PATCH(FIXPRESS,CELL,1,1,1,1,1,1,1,1)
COVAL(FIXPRESS,P1,1000.*FIXP,0.0)
COVAL(FIXPRESS,U1,ONLYMS,0.0)
COVAL(FIXPRESS,V1,ONLYMS,0.0)
COVAL(FIXPRESS,C1,ONLYMS,SAME)
GROUP 15. Termination of sweeps
LSWEEP=50
NAMSAT=MOSG
PLANTBEGIN
** Impellor cell cartesian components
PATCH(SS198VEL,CELL,1,NX,1,NY,1,NZ,1,lSTEP)
CO=1.e5
VAL=-RG(2)*(YG2D-RG(6))
COVAL(SS198VEL,U1,GRND,GRND)
CO=1.e5
VAL=RG(2)*(XG2D-RG(5))
COVAL(SS198VEL,V1,GRND,GRND)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The cartesian velocity components for each cell of the
impellor are set above by reference to their markers
indicated by PATCH names.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Impellor geometry
MARK =0.
IF(ISWEEP.EQ.1)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
First, nulify the markers values over the whole domain
at the start of first sweep.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
MARK =XYELLP(RG(1),RG(5),RG(6),RG(3),RG(4),RG(2)*TIM,0.)
IF(ISWEEP.EQ.1)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
Then, specify the ellipse of MARK=198 with the center at
the axis of rotation, having the half axis as above and
rotation angle equal to product of number of revolution
per second and current time.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
MARK =SPHERE(RG(1),RG(5),RG(6),0.5,3.0)
IF(ISWEEP.EQ.1)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The above statement represents the cylindrical shaft as
circular object centered with paddles at the middle of
Z-slab with radius equal 3. m.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
C1 = 0
IF(ISWEEP.EQ.LSWEEP.AND.MARK.EQ.198.)
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
The above statement provides the nulification of
concentration at the cells occupied by impellor at the
end of last sweep.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,0.3)
RELAX(V1,Falsdt,0.3)
relax(c1,falsdt,1.)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
GROUP 22. Spot-value print-out
TSTSWP=-1;IXMON=NX/4; IYMON=NY/4
GROUP 23. Print-out & plot control
idispa=1; csg1=p; selref=t; resfac=1.e-2
tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=6.974E+08; EX(U1)=1.161E+00; EX(V1)=1.135E+00
EX(C1)=3.418E-01; EX(MARK)=2.079E+01; EX(CONC)=1.000E-10
LIBREF=610
STOP