*add a checkOnePhaseSolution function, but can only check L2 norm for velocity now...