We developed faster methods for simulating Earth's interior dynamics. Standard iterative algorithms struggle to solve these equations efficiently. We introduced two improvements. First, we reformulated how calculation errors are measured. Second, we added a mass conservation correction. The first reduces the number of expensive solves, and the second improves accuracy over long simulations. We tested the method on several benchmarks, and both improvements come at little extra cost.