The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. Response to Reviewers: My comments to the reviewers are in the attached word-document, "comments reviewers final.docx".
Powered by Editorial Manager® and Preprint Manager® from Aries Systems CorporationComputational Geoscience manuscript No. The former ignores capillary pressure and can be solved analytically. The latter has traditionally been formulated as an iterative integral solution, which suffers from convergence problems as the injection saturation approaches unity. Here an alternative approach is presented that solves the MSE using a pseudospectral Chebyshev differentiation matrix. The resulting pseudospectral solution is compared to results obtained from the original integral implementation and the Buckley-Leverett limit, when the capillary pressure becomes negligible. A self-contained MATLAB code to implement the new solution is provided within the manuscript.The new approach offers a robust and accurate method for verification of numerical codes solving two-phase flow with capillary pressure.