Problem after discretizing a PDE to ODE and solve with DifferentialEquations.jl

[ChrisRackauckas]

Are you sure that’s a stable discretization of the PDE? Do a Von Neumann analysis.

For reference, look through the following:

PSA: How to help yourself debug differential equation solving issues

Debugging differential equation solver issues almost always boils down to doing the same thing, so this is a summary to help you out. For more information on specific issues, check out the FAQ section of the DifferentialEquations.jl documentation which highlights common issues and questions: https://diffeq.sciml.ai/dev/basics/faq/ How do I debug why the differential equation solver is diverging? dt <= dtmin. Aborting. There is either an error in your model specification or the true solution i…