TALK=T;RUN(1,1)
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
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
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)
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)
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)
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)
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.)
The above statement provides the nulification of
concentration at the cells occupied by impellor at the
end of last sweep.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT( Unsteady mixing in two paddle-stirred r)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 610
************************************************************
Group 2. Time dependence
STEADY = F
* Set overall time and no. of steps
TFIRST =0. ;TLAST =12.
FSTEP = 1 ;LSTEP = 12
TFRAC(1)=0.083333 ;TFRAC(2)=0.166667
TFRAC(3)=0.25 ;TFRAC(4)=0.333333
TFRAC(5)=0.416667 ;TFRAC(6)=0.5
TFRAC(7)=0.583333 ;TFRAC(8)=0.666667
TFRAC(9)=0.75 ;TFRAC(10)=0.833333
TFRAC(11)=0.916667 ;TFRAC(12)=1.
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 40
XULAST =20.
XFRAC(1)=0.025 ;XFRAC(2)=0.05
XFRAC(3)=0.075 ;XFRAC(4)=0.1
XFRAC(5)=0.125 ;XFRAC(6)=0.15
XFRAC(7)=0.175 ;XFRAC(8)=0.2
XFRAC(9)=0.225 ;XFRAC(10)=0.25
XFRAC(11)=0.275 ;XFRAC(12)=0.3
XFRAC(13)=0.325 ;XFRAC(14)=0.35
XFRAC(15)=0.375 ;XFRAC(16)=0.4
XFRAC(17)=0.425 ;XFRAC(18)=0.45
XFRAC(19)=0.475 ;XFRAC(20)=0.5
XFRAC(21)=0.525 ;XFRAC(22)=0.55
XFRAC(23)=0.575 ;XFRAC(24)=0.6
XFRAC(25)=0.625 ;XFRAC(26)=0.65
XFRAC(27)=0.675 ;XFRAC(28)=0.7
XFRAC(29)=0.725 ;XFRAC(30)=0.75
XFRAC(31)=0.775 ;XFRAC(32)=0.8
XFRAC(33)=0.825 ;XFRAC(34)=0.85
XFRAC(35)=0.875 ;XFRAC(36)=0.9
XFRAC(37)=0.925 ;XFRAC(38)=0.95
XFRAC(39)=0.975 ;XFRAC(40)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =20.
YFRAC(1)=0.025 ;YFRAC(2)=0.05
YFRAC(3)=0.075 ;YFRAC(4)=0.1
YFRAC(5)=0.125 ;YFRAC(6)=0.15
YFRAC(7)=0.175 ;YFRAC(8)=0.2
YFRAC(9)=0.225 ;YFRAC(10)=0.25
YFRAC(11)=0.275 ;YFRAC(12)=0.3
YFRAC(13)=0.325 ;YFRAC(14)=0.35
YFRAC(15)=0.375 ;YFRAC(16)=0.4
YFRAC(17)=0.425 ;YFRAC(18)=0.45
YFRAC(19)=0.475 ;YFRAC(20)=0.5
YFRAC(21)=0.525 ;YFRAC(22)=0.55
YFRAC(23)=0.575 ;YFRAC(24)=0.6
YFRAC(25)=0.625 ;YFRAC(26)=0.65
YFRAC(27)=0.675 ;YFRAC(28)=0.7
YFRAC(29)=0.725 ;YFRAC(30)=0.75
YFRAC(31)=0.775 ;YFRAC(32)=0.8
YFRAC(33)=0.825 ;YFRAC(34)=0.85
YFRAC(35)=0.875 ;YFRAC(36)=0.9
YFRAC(37)=0.925 ;YFRAC(38)=0.95
YFRAC(39)=0.975 ;YFRAC(40)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(16)=C1
NAME(149)=MARK ;NAME(150)=CONC
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,N,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(C1,Y,Y,N,N,N,Y)
SOLUTN(MARK,Y,N,N,N,N,Y)
SOLUTN(CONC,Y,N,N,N,N,Y)
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(C1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1000. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-02 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(C1)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(C1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(U1)=0.
FIINIT(V1)=1.0E-10 ;FIINIT(C1)=1.0E-10
FIINIT(MARK)=0. ;FIINIT(CONC)=1.0E-10
PATCH(WC1 ,INIVAL, 1, 40, 1, 20, 1, 1, 1, 12)
INIT(WC1 ,C1 ,0. ,0. )
PATCH(EC1 ,INIVAL, 1, 40, 1, 21, 1, 1, 1, 12)
INIT(EC1 ,C1 ,0. ,1. )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(WALLN ,NWALL , 1, 40, 40, 40, 1, 1, 1, 12)
COVAL(WALLN ,U1 ,1. ,0. )
PATCH(WALLS ,SWALL , 1, 40, 1, 1, 1, 1, 1, 12)
COVAL(WALLS ,U1 ,1. ,0. )
PATCH(WALLE ,EWALL , 40, 40, 1, 40, 1, 1, 1, 12)
COVAL(WALLE ,V1 ,1. ,0. )
PATCH(WALLW ,WWALL , 1, 1, 1, 40, 1, 1, 1, 12)
COVAL(WALLW ,V1 ,1. ,0. )
PATCH(FIXPRESS,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(FIXPRESS,P1 ,1000. ,0. )
COVAL(FIXPRESS,U1 ,0. ,0. )
COVAL(FIXPRESS,V1 ,0. ,0. )
COVAL(FIXPRESS,C1 ,0. , SAME )
PATCH(SS198VEL,CELL , 1, 40, 1, 40, 1, 1, 1, 12)
COVAL(SS198VEL,U1 , GRND , GRND )
COVAL(SS198VEL,V1 , GRND , GRND )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 50 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10 ;LITER(C1)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(C1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,0.3)
RELAX(V1,FALSDT,0.3)
RELAX(C1,FALSDT,1.)
RELAX(MARK,LINRLX,1.)
RELAX(CONC,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(C1)=1.0E+10 ;VARMIN(C1)=-1.0E+10
VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10
VARMAX(CONC)=1.0E+10 ;VARMIN(CONC)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
RG( 1) =198.
RG( 2) =0.523598
RG( 3) =8.
RG( 4) =1.
RG( 5) =10.
RG( 6) =10.
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=6.974E+08 ;EX(U1)=1.161
EX(V1)=1.135 ;EX(C1)=0.3418
EX(MARK)=20.790001 ;EX(CONC)=1.0E-10
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(C1,Y,N,Y,Y,Y,Y)
OUTPUT(MARK,Y,N,Y,N,N,N)
OUTPUT(CONC,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;IYMON = 10 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='P'
STOP