188 Computational Modelling in Hydraulic and Coastal Engineering
dl=fl-zb(i-1,j);
do=fo-zb(i,j+1);
du=fu-zb(i,j-1);
dc=h(i,j)-zb(i,j);
fc=h(i,j);
Qwell=0;
if n(i,j)<0
Qwell=n(i,j)/Dd^2/p;
end
hn(i,j)=h(i,j)+rr*((dr+dc)*
(fr-fc)-(dc+dl)*(fc-fl)+(do+dc)*(fo-fc)(dc+du)*(fc-fu));
hn(i,j)=hn(i,j)+Dt*Qwell;
end
chk=0;
end
end
for i=1:nx
for j=1:ny
h(i,j)=hn(i,j);
end
end
test=k
WellDepth(k)=h(nxm,nym);
end
m=1:nt;
WD=WellDepth(m);
figure, plot(m,WD,'b','Linewidth',1.5);
xlabel('Number of time steps');
ylabel('Water depth in the well [m]');
i=1:nx;
j=1:ny;
fpp=hn(i,j);
figure, surf(i,j,fpp)
figure
contourf(i,j,fpp)
PROBLEM 7.2
Modify the application by making the following suggested changes while
keeping the rest of the data unchanged. Run the simulations, and comment
and explain the results.
1. Change the hydraulic conductivity to 0.05, 0.1, 0.25 and 0.75 m/hr.
2. Change the porosity to 0.1, 0.3 and 0.4.
3. Along with the pumping well, add a 100 m 3 /hr recharging well at
node (i = 25, j = 25).