Thread Subject: 1-D advection-diffusion: how to benchmark.

Hello everyone:

I am learning Matlab and have written a program for solving the 1-D advection-diffusion equation by finite differences. I would like to check it.

If anyone has/knows of any existing code I could use to do this, I would very much appreciate knowing about it. I don't know much about the existence of analytical solutions to this problem; if anyone can advise me about this/provide a solution and relevant conditions, that would be very welcome.

I have used central differencing, with Crank-Nicholson in time, but I can change my scheme if necessary.

Thank you for any help/suggestions.

> Hello everyone:
>
> I am learning Matlab and have written a program for
> solving the 1-D advection-diffusion equation by
> finite differences. I would like to check it.
>
> If anyone has/knows of any existing code I could use
> to do this, I would very much appreciate knowing
> about it. I don't know much about the existence of
> analytical solutions to this problem; if anyone can
> advise me about this/provide a solution and relevant
> conditions, that would be very welcome.
>
> I have used central differencing, with
> Crank-Nicholson in time, but I can change my scheme
> if necessary.
>
> Thank you for any help/suggestions.

Take a look at MATLAB's pdepe if it can be used to check
the correctness of your coding.

Best wishes
Torsten.