```
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

(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.
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

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
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.)

* 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
```