Has anybody programmed a successful Crank-nicholson algorithm for integrating systems of PDEs? I want to implement a set of 6 simultaneous advective-diffusive PDEs - I am going to have a go myself but would appreciate any help
Rob Low? Really--Probably and I mean probably the spatio-temporal part of your spatio can be assumed to go to zero. Then use the same time constant on each equation reducing it to a 1st ODE and then use MathCads solver. I make this assumption because the advection is usually uniform with the solution. If this is not the case you have a non-linear dynamic model and you must isolate the dampening terms.