186 Computational Modelling in Hydraulic and Coastal Engineering
Computer code 7.2
% Example 7.2 Unconfined Aquifer with Pumping Well
% h = Water table elevation [m];
% zb(j,j) = Elevation of the bed from the reference datum
%[m];
% h(i,j) = Water surface elevation [m];
% K = Permeability coefficient [m/hr];
% p = Porosity;
% n = Node parameter (interior, imperveous, constant head,
% well site);
% Dd = Spatial step (the same in both x and y directions
%[m];
% Dt = Time step [hr];
% nx = Number of steps in the x-direction;
% ny = Number of steps in the y-direction;
clc; clear all; close all;
% Input data;
Dd=500;
Dt=24;
nx=30;
ny=30;
nt=7200;
K=0.5;
p=0.2;
rr=Dt/Dd^2/2*K/p;
% Initialization of variables;
for i=1:nx
for j=1:ny
h(i,j)=50;
zb(i,j)=0;
% n=0 indicates the interior nodes to be
% calculated;
n(i,j)=0;
end
end
nxm=nx/2;
nym=ny/2;
% Definition of the boundary conditions;
for j=1:ny
% n=1 indicates impervious boundaries;
n(1,j)=1;
n(nx,j)=1;
end
for i=1:nx
n(i,1)=1;
% n different than 0 and 1 indicates Piezometric head in
% meters;
n(i,ny)=50;
end