On a singular diffusion equation originating from porous media flow

Article type: Research Article

Authors: van Duijn, C.J. | Floris, F.J.T.

Affiliations: Department of Mathematics, Delft University of Technology, P.O. Box 5031, 2600 GA Delft, The Netherlands | TNO Institute of Applied Geoscience, P.O. Box 6012, 2600 JA Delft, The Netherlands

Abstract: In this paper we consider the singular boundary value problem, −z1/pz″=D(f)/(1+p) for 0<f<1 with z(0)=z(1)=0. Here p>0 and D: [0,1]→R is a given function. This problem arises in a model for two-phase capillary induced flow in porous media. Considering the special case D(f)=fαw(1−f)α0, with α0,αw>−2, we investigate the singular behaviour of the solution z(f) as α0,αw↓−2. We show that the solution then becomes unbounded. We investigate the behaviour of z and z′ in this limit process. The results are incorporated in an algorithm which we use to solve the problem numerically. The numerical results show significant improvement over standard discretisation techniques near the limit. Non-existence arises for α0 or αw≤−2.

DOI: 10.3233/ASY-1995-10102

Journal: Asymptotic Analysis, vol. 10, no. 1, pp. 29-48, 1995