Abstract: | This paper presents a numerical study of infiltration processes on stratified porous media. The study is carried out to examine the performance of a finite volume method on problems with discontinuous solutions due to the transmission conditions in the interfaces. To discretize the problem, a curvilinear co‐ordinate system is used. This permits matching the interface with the boundary of the control volumes that interchange fluxes between layers. The use of the multigrid algorithm for the resulting systems of equations allows problems involving a large number of nodes with low computational cost to be solved. Finally, some numerical experiments, which show the capillary barrier behaviour depending on the material used for the different layers and the geometric design of the interface, are presented. Copyright © 2001 John Wiley & Sons, Ltd. |