A Fast Explicit Operator Splitting Method for Modified Buckley Leverett Equations
In this paper, we propose a fast explicit operator splitting method to solve the modified Buckley Leverett equations which include a third order mixed derivatives term resulting from the dynamic effects in the pressure difference between the two phases. The method splits the original equation into two equations, one with a nonlinear convective term and the other one with high order linear terms so that appropriate numerical methods can be applied to each of the split equations: the high order linear equation is numerically solved using a pseudo spectral method, while the nonlinear convective equation is integrated using the Godunov type central upwind scheme. The spatial order of the central upwind scheme depends on the order of the piecewise polynomial reconstruction: we test both the second order minmod based reconstruction and fifth order WENO5 one to demonstrate that using higher order spatial reconstruction leads to more accurate approximation of solutions. A variety of numerical examples in both one and two space dimensions show that the solutions may have many different saturation profiles depending on the initial conditions, diffusion parameter, and the third order mixed derivatives parameter. The results are consistent with the study of traveling wave solutions and their bifurcation diagrams.
J. Sci. Comput.