Kyoto, Japan

McWhorter and Sunada Exact Solution [testing]

This is the online implementation of the 1D integral solution of two-phase flow in porous media by McWhorter and Sunada (1990) and Fučík et al. (2007). Please, cite our work (see references below).
Warning, this is software is experimental: the accuracy of the solution is not guaranteed. Please, report errors to


If you see this, please, leave this form field blank and invest in CSS support.

Model selection

Model selection

Material properties
$K$[m$^2$]intrinsic permeability
$S_{wr}$[1]residual saturation of the wetting phase
$S_{nr}$[1]residual saturation of the non-wetting phase
Fluid properties
$\mu_w$[kg/m/s]dynamic viscosity of the wetting phase
$\mu_n$[kg/m/s]dynamic viscosity of the non-wetting phase
Brooks and Corey model parameters
$p_d$[Pa]Brooks and Corey model parameter: entry pressure
$\lambda$[1]Brooks and Corey model parameter: pore size distribution index
van Genuchten mode parameters
$\alpha$[1/Pa]van Genuchten model parameter
$m$[1]van Genuchten model parameter
$n$[1]van Genuchten model parameter
Problem parameters
$S_0$[1]injection boundary saturation
$S_i$[1]initial saturation
$R$[1]flux parameter
Computation parameters
$t$[s]the solution will be plotted and exported at this time
nodeslength of the discrete vector
max_itermaximum number of iterations
$psilon$stopping criterion of the iterations

Computational method selection


Copyright © Radek Fučík. All rights reserved. Version 2015.12