a recursive algorithm for connectivity analysis in a grid; application to 2d hydrodynamic modeling in heterogeneous soils：网格连通性分析的递推算法在非均匀土中二维水动力模型的应用.pdfVIP

Computers Geosciences 26 (2000) 121±135 A recursive algorithm for connectivity analysis in a grid; application to 2D hydrodynamic modeling in heterogeneous soils p B. Cappelaere*, J. Touma, C. Peugeot IRD±UMR HSM 5569, Montpellier, France Received 15 October 1998; received in revised form 26 April 1999; accepted 26 April 1999 Abstract A prerequisite for numerical simulation of water ¯ow in heterogeneous soils is to build a discrete model of the soil matrix that is a fair representation of the heterogeneities under study, while being compatible with the numerical equations used to compute unsaturated water ¯ow. When introducing signi®cant amounts of very coarse solid elements (gravels) in a discrete, multi-dimensional soil matrix model, in the form of internal boundaries that occupy a certain fraction of the grid nodes, the need arises to eliminate from the model the possible occurrence of isolated areas of soil, surrounded by continuous gravel barriers that keep them separate from other regular grid nodes. This situation is obviously an artefact resulting from the discrete representation of the physical system, since all nongravel areas ought to be considered as being submitted to at least some hydrodynamic linkage with each other and with the outer boundary conditions (rain and other water input at the top, gravity drainage or water table at the bottom). Computational nodes that do not connect in some way, through the computational grid, to these outer boundary conditions, are the source of computational problems due to system indetermination when solving the hydrodynamic equations. A method has thus been devised to