Praha-Prosek, ČR

Multidimensional Exact Solution of Two-Phase Flow in Porous Media

This is the online implementation of the d-dimensional integral solution of two-phase flow in porous media by Fučík et al. (submitted). 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
ddimension (can be any real number $\geq 2$)
$A$[$m^{d}s^{-\frac{d}{2}}$]injection rate constant
$S_i$[1]initial saturation
$t$[s]the solution will be plotted and exported at this time
nodeslength of the discrete vector
$psilon_{g}$stopping criterion of the bisection method
max_glob_itermaximum number of iterations for the bisection method
$psilon$stopping criterion of local iterations
max_itermaximum number of local iterations


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