Numerical Simulation of the Water Infiltration into a Layered Soil

摘要

The Richards' equation represents the mathematical model of water flow into porous media. In the case of layered soils the different hydraulic properties of the layers lead to the discontinuity of the equation's coefficients. We use a numerical method which consists in the finite volume scheme with imposed flux continuity for the space discretization and in the implicit backward differentiation formulae (BDF) for the time integration. The numerical results are also presented in oreder to illustrate the performance of the method.