Update `Extrapolate` function inside incompressible NS solver to use new boundary functions
As of !539 (merged) we have additional functions to deal with boundary elements that could be used to radically simplify the construction of the high-order pressure boundary conditions. The Extrapolate
classes should therefore be updated to use these new functions.