Environmental Engineering Reference
In-Depth Information
8
<
H
2
þ
'ð
x
Þ
for the confined aquifer
KH
hðxÞ¼
r
2
K
'ðxÞ
(13.15)
:
for the unconfined aquifer
The entire M-file is given by:
Qy0 = -1.e-6; Qx0 = 0.0;
K = 0.001;
h0 = 9.99;
H = 10;
[x,y] = meshgrid ([0:2:200],[0:2:200]);
phi = -Qx0*x - Qy0*y;
if
(h0>=H)
phi0 = -phi(1,1) + K*H*h0 - 0.5*K*H*H;
% confined
else
phi0 = -phi(1,1) + 0.5*K*h0*h0;
% unconfined
end
phi = phi + phi0;
phicrit = 0.5*K*H*H + phi0;
% margin between confined and unconfined
confined = (phi>=phicrit);
h = confined.*(0.5*H+(1/K/H)*phi)+~confined.*sqrt((2/K)*phi);
surf (x,y,h);
In the final part of the M-file we use a matrix containing the values 0 and 1,
depending on the hydraulic situation of the aquifer at the corresponding mesh
position. The array
con
ned
contains a 1 if the aquifer is confined, and a 0 if it is
Fig. 13.3 Head distribution for the reference input data
