TALK=F;RUN( 1, 1)
  When viewing via a browser. click
  here for
  explanation
  
#cls
TEXT(Flow In Stirred 3D Vessel
TITLE 
LIBREF=757
  DISPLAY
  
  This In-Form case uses "IMAT" In-Form option for simulation
  flow in stirred 3D vessel.

  The following statements make the vertical and horizontal
  velocities proportional to the corresponding distances from
  the paddle axis.

  (SOURCE of U1 at EI is 1.E10*(VEL*(YIC-YG)-U1) with IMAT>=100!LINE)
  (SOURCE of V1 at NI is 1.E10*(VEL*(XG-XIC)-V1) with IMAT>=100!LINE)

  These momentum sources are active only for cell having IMAT
  (the property index, as set by PRPS) values greater or equal
  to 100 .

  The following In-Form formulae describe density and viscosity:
    
  (PROPERTY of RHO1 is DENS)
  (PROPERTY VISL is VISC)

  The next statement sets PRPS to the IMAT value of water
  throughout the domain; but the value is changed below for
  the cells ocupied by the paddle
       
  (STORED of PRPS is 67.)

  Cells occupied by paddle material are marked by setting the
  value of PRPS in them to 100.

  (STORED of PRPS at TIME=1 is :MAT:)
  (STORED of PRPS at TIME=2 is :MAT:)
  (STORED of PRPS at TIME=3 is :MAT:)
  (STORED of PRPS at TIME=4 is :MAT:)
  (STORED of PRPS at TIME=5 is :MAT:)
  (STORED of PRPS at TIME=6 is :MAT:)
  (STORED of PRPS at TIME=7 is :MAT:)
  (STORED of PRPS at TIME=8 is :MAT:)

  The Q1 contains PHOTON USE commands
  ENDDIS  

    GROUP 1. Run title and other preliminaries
TEXT(Flow in stirred 3D vessel: 757
  INFORM1BEGIN
REAL(PROP);PROP=50.
    ** paddle cells material
REAL(MAT); MAT  =100.
     Angular velocity, 1/sec
REAL(VEL); VEL=2/TLAST
     ** Origin of rotation axes, m
REAL(XIC, YIC); XIC=7.5; YIC=7.5
     ** Density
REAL(DENS); DENS=1000.
     ** Viscosity
REAL(VISC); VISC=1.
  INFORM1END
 
  GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,8,8.,1.0)
 
    GROUP 3. X-direction grid specification
GRDPWR(X,15,15.,1.0)

    GROUP 4. Y-direction grid specification
GRDPWR(Y,15,15.,1.0)

    GROUP 5. Z-direction grid specification
GRDPWR(Z,15,15.,1.0)

    GROUP 7. Variables stored, solved , named
SOLVE(P1,U1,V1,W1); STORE(PRPS)

    GROUP 9. Properties of the medium (or media)
  INFORM9BEGIN
RHO1=DENS; ENUL=VISC
(PROPERTY of RHO1 is DENS)
(PROPERTY VISL is VISC)
  ** water 
FIINIT(PRPS)=67
  INFORM9END

    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0;FIINIT(U1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0
 
    GROUP 13. Boundary conditions and special sources
PATCH(PRESREL,CELL,NX,NX,NY,NY,NZ,NZ,1,LSTEP)
COVAL(PRESREL,P1,1000.,0.0)
  INFORM13BEGIN
  ** paddle cell cartesian components
PATCH(EI,EAST,1,NX,1,NY,1,NZ,1,lSTEP)
(SOURCE of U1 at EI is 1.E10*(VEL*(YIC-YG)-U1) with IMAT>=100!LINE)
PATCH(NI,NORTH,1,NX,1,NY,1,NZ,1,lSTEP)
(SOURCE of V1 at NI is 1.E10*(VEL*(XG-XIC)-V1) with IMAT>=100!LINE)
  INFORM13END
  INFORM7BEGIN
  ** Water 
(STORED of PRPS is 67.)
  
  ** paddle geometry
   * time: 1 sec.
PATCH(TIME=1,CELL,8,8,4,12,6,10,1,1)
(STORED of PRPS at TIME=1 is :MAT:)

   * time: 2 sec.
PATCH(TIME=2,CELL,4,12,8,8,6,10,2,2)
(STORED of PRPS at TIME=2 is :MAT:)

   * time: 3 sec.
PATCH(TIME=3,CELL,8,8,4,12,6,10,3,3)
(STORED of PRPS at TIME=3 is :MAT:)

   * time: 4 sec.
PATCH(TIME=4,CELL,4,12,8,8,6,10,4,4)
(STORED of PRPS at TIME=4 is :MAT:)

   * time: 5 sec.
PATCH(TIME=5,CELL,8,8,4,12,6,10,5,5)
(STORED of PRPS at TIME=5 is :MAT:)

   * time: 6 sec.
PATCH(TIME=6,CELL,4,12,8,8,6,10,6,6)
(STORED of PRPS at TIME=6 is :MAT:)

   * time: 7 sec.
PATCH(TIME=7,CELL,8,8,4,12,6,10,7,7)
(STORED of PRPS at TIME=7 is :MAT:)

   * time: 8 sec.
PATCH(TIME=8,CELL,4,12,8,8,6,10,8,8)
(STORED of PRPS at TIME=8 is :MAT:)
  INFORM7END
    GROUP 15. Termination of sweeps
LSWEEP=50
    GROUP 19. Data communicated by satellite to GROUND
SPEDAT(SET,MATERIAL,100,L,T) 
    GROUP 22. Spot-value print-out
NPRMON=LSWEEP;IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1
    GROUP 23. Print-out , plot control
TSTSWP=-1;idispa=1; csg1=p; selref=t; resfac=1.e-2
DISTIL=T
EX(P1)=4.077E+03; EX(U1)=1.238E+00
EX(V1)=1.336E+00; EX(W1)=8.445E-01
STORE(PRPS)
EX(PRPS)=6.744E+01
  PHOTON USE
  p
  p1;
 
  set prop off ; vi 1 1 5
 
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 1 sec.
  * con v1 y 5 fi;0.001;*dump picb1
  pause
  p
  p2;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 2 sec.
  * con v1 y 5 fi;0.001;*dump picb2
  pause
  p
  p3;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 3 sec.
  * con v1 y 5 fi;0.001;*dump picb3
  pause
  p
  p4;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 4 sec.
  * con v1 y 5 fi;0.001;*dump picb4
  pause
  p
  p5;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 5 sec.
  * con v1 y 5 fi;0.001;*dump picb5
  pause
  p
  p6;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 6 sec.
  * con v1 y 5 fi;0.001;*dump picb6
  pause
  p
  p7;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 7 sec.
  * con v1 y 5 fi;0.001;*dump picb7
  pause
  p
  p8;
 
  set prop off ; vi 1 1 5
  surf prps z 89 
  gr ou z 8;gr ou y 8;gr ou y m;gr ou z m;gr ou x m
  vec z 8 sh; vec y 1 sh; vec x 1 sh
  msg at 8 sec.
  * con v1 y 5 fi;0.001;*dump picb8
  ENDUSE
STOP