A Simulation Model of the Gravity Separation of Concentrated Solids
469
420
VEL(L)=PI*FI**3/(AK*(1 .-FI)**2)
430
PS(L)=PI
440
300 P(L)=PI
450
WRITE(6,*) AK.VELW.PI.FI.FUJLTH.CO.VIN.VO)
460
TIME=0.
470
NN=0
480
Kl=500
490
K2=501
500
K86=0
510
NN1=0
520
K3=l
530
K85=0
540
DO 807 J= 1,10000
550
J2=(J+3)/4
560
TIME=TIME+DT
570
NN=NN+1-NN1
580
DO 101=1,10
590
SUM=0.
600
DO201N=l,NN
610
L=K1-NN+N
620
UP=(F(L)-FU)/(FI-FU)
630
K=AK + (R-1 .)* AK*( 1 -UP)
640
P(L)=U(L)*PS(L-NN1)+(1.-U(L-1))*PS(L-1-NN1)*UP
650
IF(F(L).LT.FU) F(L)=FU
660
VEL(L)=P(L)*F(L)**3/(K*(1 .-F(L))**2)
670
201 CONTINUE
680
D O 92L=1,M
690
92 PS(L)=P(L)
700
DO200N=l,NN
710
L=K2-N
720
IF(J.EQ.1)SDT=1.
730
IF(J.EQ.1)TIME=1.
740
IF(J.GT.1)SDT=8.
750
VTR(L)=VEL(L-NN1)*SDT
760
VR(L)=VTR(L)-VTR(L+1)
770
VN(L)=V(L-NN1)-VR(L)
780
F(L)-(VN(L)-(1 .-F0LD(L-NN1))*V(L-NN1))/VN(L)
790
IF(F(L).GT.FI) F(L)=FI
800
CON(L)=V(L-NN 1 )/VN(L)*CON(L-NN 1)
810
V(L)=VN(L)
820
IF(FOLD(L).EQ.FU)GO T O 308
830
IF(F(L)-FU) 308,308,309
840
308 F(L)=FU
850
U(L)=0.
860
G O T O 310
870
309 U(L)=(F(L)-FU)/(F0LD(L-NN1)-FU)
880
310 CONTINUE
890
IF(U(L).GT. 1.) U(L)= 1.