Demonstrate approximate Benzi preconditioner for grad-div perturbations to the velocity block of Navier-Stokes

This example is based on https://epubs.siam.org/doi/abs/10.1137/22M1505529. It tests one new API which provides the ability to use a matrix other than the identity in between applications of the two special composite preconditioners. We will be attempting to use this in MOOSE for high Reynolds number flows as it provides an entirely algebraic alternative to some of the excellent work of Farrell et. al. using geometric multigrid