simulation_type 2 seems to be working. A time dependent mesh is used and optimal convergence rates are still seen (using backwards differences).