Encyclopaedia Index
TALK=F;RUN( 1, 1);VDU=VGAMOUSE
TEXT( Library case Y617: Unsteady mixing in two paddle-stirred reactor
>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
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 <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
** 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 • direction in Y direction
REAL(AA , BB )
+ AA =8.0 ; BB =1.0
+ RG(3)=AA ; RG(4)=BB
* Rotation axis coordinates, ¬
• 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
** 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
REGION() 198 /ISWEEP.EQ.LSWEEP
>>>>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>
The above statement provides the nulification of concentration at
the cells occupied by impellor at the end of last sweep.
<<<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,0.3)
RELAX(V1,Falsdt,0.3)
relax(c1,falsdt,1.)
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
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
STOP