Add the option of advecting theta. You can register an AdvDiffHierarchyIntegrator object that knows how to advect theta. Also adds in an initial four roll mill problem with the make target four_roll_mill.
This runs, but if theta is small or large, i.e. thn > 0.8 or thn < 0.2, we need to crank up the preconditioner sweeps considerable. I don't quite understand why yet. We should investigate the order of accuracy of the time integrator before the advection diffusion integrator is merged.
Add the option of advecting theta. You can register an
AdvDiffHierarchyIntegratorobject that knows how to advect theta. Also adds in an initial four roll mill problem with the make targetfour_roll_mill.This runs, but if theta is small or large, i.e. thn > 0.8 or thn < 0.2, we need to crank up the preconditioner sweeps considerable. I don't quite understand why yet. We should investigate the order of accuracy of the time integrator before the advection diffusion integrator is merged.