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