Improvement of fully-implicit two-phase pore-network models by employing generalized flux functions with additional throat variables

In fully-implicit two-phase pore-network models, developing a well-converged scheme remains a major challenge, primarily due to the discontinuities in the phase conductivities. This paper addresses these numerical issues by proposing a generalized flux function that establishes a continuous flux expression for two-phase flows by introducing an additional throat variable $Θ$. Two approaches for expressing this additional throat variable are introduced: the first applies regularization strategies, while the second constructs an additional residual constraint equation. It is shown that this approach significantly improves accuracy and ensures the temporal convergence, as demonstrated through various numerical examples.