Original suggestion here #676, with #872 we could have
fields = (u = Lagrange{RefTriangle, 2}()^2, p = Lagrange{RefTriangle, 1}())
qr = QuadratureRule{RefTriangle}(2)
ipg = geometric_interpolation(grid)
dh = DofHandler(grid)
add!(dh, fields)
close!(dh)
cmv = MultiFieldCellValues(qr, fields, ipg)
(And for convenience we could also support dh = DofHandler(grid, fields) which closes automatically which is nice for quick examples :))