Matlab Problem -

Assume you have a 20m*50m rectangular pond. A pollutant enters a 1m*1m rectangular section at the center of pond and initially has the solute concentration of 1000[ML^-3], while the solute concentration at other sections of this pond is zero.

At the boundaries the concentration of the pollutant will remain zero. Following the Advection dispersion formula stated below:

∂C/∂t = D_x ∂²C/∂x² + D_y ∂²C/∂y²

where:

Dx = 1.0, Dy = 0.2

Evaluate and demonstrate using appropriate mathematical tools and graphs how the solute concentration varies in time and how the pollutant is advected and dispersed in this pond.

For solving the above equation use the following discretization methods:

(1) Explicit technique

(2) Crank Nicholson method.

And use iterative techniques (three iterative methods were explained in the class, you may provide results only using Gauss seidel method) for solving your system of linear equations.

And compare your results based on the discretization method you use.

Note that for explicit technique the stability requirement is the following:

Δt (1/Δx² + 1/Δy²) ≤ 1/2 Δt(1/Δx² + 1/Δy²) ≤ 1/2

and the Crank Nicholson method is unconditionally stable.

In preparing the report you need to clearly explain your solution algorithm. You also need to specify sections in your code corresponding to each step.

Your graphs need to be analyzed and your results need to be explained clearly.

In the manual report you have to give the mentioned method's (Explicit and Crank Nicholson) at least one trial by explaining procedure and steps solving the problem. And then in the matlab you have to give explanation of each step and in the ouput you have to provide all the graphs tables etc.