Efficient reordered nonlinear Gauss–Seidel solvers with higher order for black-oil models