Contaminant and sediment transport by advection and diffusion 219
diff=sqrt(6*dif/Dt);
% Main program;
for k=1:nt
% Movement of particles by advection and diffusion;
for i=1:ipp
X=fix(x(i,1)/Dd)+1;
Y=fix(y(i,1)/Dd)+1;
if X>2 && X2 && Y u=fm(X,Y,fCoU);
ur=fm(X+1,Y,fCoU);
v=fm(X,Y,fCoV);
vo=fm(X,Y+1,fCoV);
uu=(u+ur)/2;
vv=(v+vo)/2;
tempx=x(i,1);
tempy=y(i,1);
ax=rand;
x(i,1)=x(i,1)+(2*(ax-0.5)*diff+uu)*Dt;
ay=rand;
y(i,1)=y(i,1)+(2*(ay-0.5)*diff+vv)*Dt;
X=fix(x(i,1)/Dd)+1;
Y=fix(y(i,1)/Dd)+1;
% Deflection on solid boundarues;
h=fm(X,Y,fDep);
if h<0.1
x(i,1)=tempx;
y(i,1)=tempy;
end
end
end
% Decaying of particles;
if dec>0
for i=1:depp
ad=rand;
X=fix(ad*ipp)+1;
x(X,1)=10;
y(X,1)=50;
end
end
% Computation of particle concentration in grid cells;
for i=1:nx
for j=1:ny
c(i,j)=0;
end
end
for i=1:ipp
X=fix(x(i,1)/Dd)+1;
Y=fix(y(i,1)/Dd)+1;
c(X,Y)=c(X,Y)+1;
end