Mesh Adaptation

This brings in mesh adaptivity capabilities from libMesh. Not all available error estimators are there, but at least 1 adjoint-based one is as well as a couple of global error estimators. Currently only have a solver for the steady case. cavity_benchmark example now uses adaptivity.