Asymptotical computations for a model of flow in saturated porous media

We discuss an initial value problem for an implicit second order ordinary differential equation which arises in models of flow in saturated porous media such as concrete. Depending on the initial condition, the solution features a sharp interface with derivatives which become numerically unbounded. By using an integrator based on finite difference methods and equipped with adaptive step size selection, it is possible to compute the solution on highly irregular meshes. In this way it is possible to verify and predict asymptotical theory near the interface with remarkable accuracy.

[1]  C. Budd,et al.  ASYMPTOTIC BEHAVIOUR OF WETTING FRONTS IN POROUS MEDIA WITH EXPONENTIAL MOISTURE DIFFUSIVITY , 2008 .

[2]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[3]  Ewa Weinmüller,et al.  Numerical Solution of Singular Two Point BVPs 1 , 2009 .

[4]  J. Parlange,et al.  A Note on A Three-Parameter Soil-Water Diffusivity Function-Application to the Horizontal Infiltration of Water , 1973 .

[5]  Peter F. Dux,et al.  Unsaturated diffusivity functions for concrete derived from NMR images , 2003 .

[6]  D. K. Babu.  Infiltration analysis and perturbation methods: 1. Absorption with exponential diffusivity , 1976 .

[7]  Mathias J. M. Römkens,et al.  Extension of the Heaslet‐Alksne Technique to arbitrary soil water diffusivities , 1992 .

[8]  Wilfried Brutsaert,et al.  Universal constants for scaling the exponential soil water diffusivity , 1979 .

[9]  L. Brugnano,et al.  Solving differential problems by multistep initial and boundary value methods , 1998 .

[10]  D. Lockington,et al.  A new perturbation expansion for horizontal infiltration and sorptivity estimates , 1988 .

[11]  John M. Stockie,et al.  Multilayer Asymptotic Solution for Wetting Fronts in Porous Media with Exponential Moisture Diffusivity , 2015, 1505.07197.

[12]  John Crank,et al.  The Mathematics Of Diffusion , 1956 .

[13]  M. E. Rose,et al.  Numerical methods for flows through porous media. I , 1983 .

[14]  Pierluigi Amodio,et al.  High order generalized upwind schemes and numerical solution of singular perturbation problems , 2007 .

[15]  D. K. Babu.  Infiltration analysis and perturbation methods: 3. Vertical infiltration , 1976 .

[16]  Pierluigi Amodio,et al.  A finite differences MATLAB code for the numerical solution of second order singular perturbation problems , 2012, J. Comput. Appl. Math..

[17]  Giuseppina Settanni,et al.  Variable Step/Order Generalized Upwind Methods for the Numerical Solution of Second Order Singular Perturbation Problems 1 2 , 2009 .

[18]  Pierluigi Amodio,et al.  High-order finite difference schemes for the solution of second-order BVPs , 2005 .

[19]  Michael Frankfurter,et al.  Numerical Recipes In C The Art Of Scientific Computing , 2016 .

[20]  Yuanle Ma,et al.  Computational methods for multiphase flows in porous media , 2007, Math. Comput..

[21]  Othmar Koch,et al.  Efficient Numerical Solution of the Density Profile Equation in Hydrodynamics , 2007, J. Sci. Comput..

[22]  Brent Clothier,et al.  Measurement of Sorptivity and Soil Water Diffusivity in the Field , 1981 .

[23]  Pierluigi Amodio,et al.  High Order Finite Difference Schemes for the Solution of Second Order Initial Value Problems , 2010 .