Kamakura, Japan

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


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
$\epsilon_{g}$stopping criterion of the bisection method
max_glob_itermaximum number of iterations for the bisection method
$\epsilon$stopping criterion of local iterations
max_itermaximum number of local iterations


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