Advection-diffusion equation q_t + u*q_x = D*q_{xx} solved with a fractional step method. The diffusion equation is solved in src1.f using Crank-Nicolson. The tridiagonal system is solved using the LAPACK routine in lapacktridiag.f.

The advection velocity u and diffusion coefficient dcoef are set in setprob.data.

afterframe.m computes and plots the exact solution along with the computed solution.

Animation of computed and exact solution.

See also:

- claw/applications/burgers/1d/viscous Viscous Burgers' equation.

Download directory in compressed tar file format (Note: Please unpack the tar file in your CLAW home directory.)