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 radek.fucik@fjfi.cvut.cz.

Parameters

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

Model selection
primary variable is the wetting phase saturation $S_w$
primary variable is the non-wetting phase saturation $S_n$

Model selection
Brooks and Corey model
van Genuchten model

Material properties
$\phi$	[1]	porosity
$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
nodes		length of the discrete vector
max_iter		maximum number of iterations
$\epsilon$		stopping criterion of the iterations

Computational method selection
Original method (McWhorter and Sunada 1990)
Modified method - variant A (Fučík et al. 2007)
Modified method - variant B (Fučík et al. 2007)

References

link R. Fučík, T. H. Illangasekare, and M. Beneš Multidimensional self-similar analytical solutions of two-phase flow in porous media, Advances in Water Resources, Volume 90, April 2016, Pages 51–56 doi.org/10.1016/j.advwatres.2016.02.007
pdf R.Fučík, J. Mikyška, T. H. Illangasekare and M. Beneš: Semianalytical Solution for Two-Phase flow in Porous Media with a Discontinuity Vadose Zone Journal 2008 vol. 7, no. 3: pages 1001-1009
pdf R.Fučík, J. Mikyška, T. H. Illangasekare and M. Beneš: An Improved Semi-Analytical Solution for Verification of Numerical Models of Two-Phase Flow in Porous Media Vadose Zone Journal 2007 no. 6: pages 93-104
pdf (1,01 MB) corrected version