Added simple DiffusionSolver added a bit of code main simulation loop