Macroscopic representation of structural geometry for simulating water and solute movement in dual-porosity media

The structure of macroporous or aggregated soils and fractured rocks is generally so complex that it is impractical to measure the geometry at the microscale (i.e., the size and the shape of soil aggregates or rock matrix blocks, and the myriad of fissures or fractures), and use such data in geometry-dependent macroscale flow and transport models. This paper analyzes a first-order type dual-porosity model which contains a geometry-dependent coefficient, β, in the mass transfer term to macroscopically represent the size and shape of soil or rock matrix blocks. As a reference, one- and two-dimensional geometry-based diffusion models were used to simulate mass transport into and out of porous blocks of defined shapes. Estimates for β were obtained analytically for four different matrix block geometries. Values for β were also calculated by directly matching analytical solutions of the diffusion models for a number of selected matrix block geometries to results obtained with the first-order model assuming standard boundary conditions. Direct matching improved previous results for cylindrical macropore geometries, especially when relatively small ratios between the outer soil mantle and the radius of the inner cylinder were used. Results of our analysis show that β is closely related to the ratio of the effective surface area available for mass transfer, and the soil matrix volume normalized by the effective characteristic length of the matrix system. Using values of β obtained by direct matching, an empirical function is derived to estimate macroscopic geometry coefficients from medium properties which in principle are measurable. The method permits independent estimates of β, thus allowing the dual-porosity approach eventually to be applied to media with complex and mixed types of structural geometry.