1. As we discussed in class, solve the diffusion equation in 1D using finite differences. Use the dike injection scenario with cooling as a function of time. Inject a second dike at x=10 m some time after the inital dike injection. Can you solve this problem analytically?
2. Solve the diffusion advection equation in 2D for a given velocity field. Consider the animation above. A point source of some solute (say Ar) is released in a fixed velocity field - shown by the streamlines (which might not look continuous on your browser).
a. modify the program (given below) to calculate the downstream concentration at several (3) points as a function of time for an instantaneous release at the point source. Plot your results as three curves on a single graph.
b. modify the program again to have a continuous point source and plot the same three curves.
c. Discuss your results
d. What happens to the results when you change the diffusivity and / or the grid cell spacing? Be careful and make small changes.
You need the following files
Turn in everything in one MS document. Due as soon as you can get it done (and absolutely before April 27).
Chuck Connor