Module 5: GROUNDWATER TRANSPORT PROCESS

Lecture 37: Solution of Advection-Diffusion equation using MATLAB

Comparing (37.16) with (37.2), we have

co = 0
(37.23)

Comparing (37.17) and (37.18) with (37.2), we have

pl = cl -1

ql = 0

pr = cr

qr = 0

(37.24)

(37.25)

(37.26)

(37.27)

Now the Matlab functions can be written as

Function pdefun

Function icfun

Function bcfun