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