Dual-Permeability Module
The DualPerm Module was developed as a supplemental module of the HYDRUS (2D/3D) software package (Version 2), to model the two-dimensional variably-saturated water movement and solute transport in dual-permeability porous media.
Warning: This module is usually significantly less numerically stable than the regular HYDRUS code due to much larger nonlinearity usually associated with the fracture domain and it therefore requires much finer spatial discretization then the standard HYDRUS. It is thus intended mainly for more experienced HYDRUS users!
Preferential flow in structured media (both macroporous soils and fractured rocks) can be described using a variety of dual-porosity, dual-permeability, multi-porosity, and/or multi-permeability models. Dual-porosity and Dual-Permeability models both assume that the porous medium consists of two interacting regions, one associated with the inter-aggregate, macropore, or fracture system, and one comprising micropores (or intra-aggregate pores) inside soil aggregates or the rock matrix. While dual-porosity models assume that water in the matrix is stagnant, dual-permeability models allow for water flow in the matrix as well. While dual-porosity models are available in the standard version of HYDRUS (2D/3D), dual-permeability models are not.
Dual-porosity models have long been applied to solute transport studies. Especially popular early on were dual-porosity models in which distinct mobile and immobile flow regions are assumed to be present. Dual-permeability models in which water can move in both the inter- and intra-aggregate pore regions are now also becoming more popular. Available dual-permeability models differ mainly in how they implement water flow in and between the two pore regions, especially with respect to the degree of simplification and empiricism. Approaches to calculating water flow in macropores or inter-aggregate pores range from those invoking Poiseuille’s equation, the Green and Ampt or Philip infiltration models, the kinematic wave equation, and the Richards equation (Gerke and van Genuchten, 1993a).
We have implemented the dual-permeability model based on the approach suggested by Gerke and van Genuchten (1993a). The dual-permeability formulation for water flow is based on a mixed form of the Richards equation, describing water flow in both the fractures (macropores) and the matrix (micropores) domains [Gerke and van Genuchten, 1993a]. The dual-permeability formulation for solute transport is based on a convection-dispersion equation, describing solute transport in both the fractures (macropores) and the matrix (micropores) domains [Gerke and van Genuchten, 1993a]. The mass transfer of water between the two domains is driven by the gradient of pressure heads. The mass transfer for solute includes both convective mass transfer with water mass transfer, as well as diffusive mass transfer driven by the concentration gradient. See the detailed description of this model either in references given below or in the DualPerm Module manual.
|
|
Download the DualPerm Module manual (pdf).
Šimůnek, J., M. Šejna, and M. Th. van Genuchten, The DualPerm Module for HYDRUS (2D/3D) Simulating Two-Dimensional Water Movement and Solute Transport in Dual-Permeability Porous Media, Version 1.0, PC Progress, Prague, Czech Republic, 32 pp., 2012.
- Project Group: Dual-Permeability Module
- Description: Examples demonstrating the use of the special module of HYDRUS that considers dual-permeability model
- Availability: Download HYDRUS projects now: zip (5.11 MB)
- Note: These projects were created with version 2, and users using higher Hydrus versions need to convert them to their particular version.
Project
|
Description
|
Test1a
|
Dual-Permeability Model, Infiltration into fracture, a=1.0 cm.
|
Test1a1
|
Dual-Permeability Model, Infiltration into matrix, a=1.0 cm.
|
Test1a2
|
Dual-Permeability Model, Infiltration into matrix, a=1.0 cm. Higher flux.
|
Test1a3
|
Same as Test1a1 + solute pulse.
|
Test1a4
|
Dual-Permeability Model, Atmospheric BC, Root Water Uptake, a=1.0 cm
|
Test1b
|
Dual-Permeability Model, Infiltration into fracture, a=0.1 cm.
|
Test1c
|
Dual-Permeability Model, Infiltration into fracture, a=0. 0316 cm.
|
Test1d
|
Same as Test1a + solute transport.
|
Test1e
|
Same as Test1b + solute transport.
|
Test1f
|
Same as Test1c + solute transport.
|
Test2a
|
Dual-Permeability Model, Surface ponding, a=1.0 cm.
|
Test2a
|
Dual-Permeability Model, Surface ponding, a=1.0 cm.
|
Test2b
|
Dual-Permeability Model, Surface ponding, a=0.2 cm.
|
Test2c
|
Dual-Permeability Model, Surface ponding, a=3.0 cm.
|
Test2d
|
Dual-Permeability Model, Surface ponding, a=1.0 cm, Solute Transport.
|
Test2e
|
Dual-Permeability Model, Surface ponding, a=0.2 cm, Solute Transport.
|
Test2f
|
Dual-Permeability Model, Surface ponding, a=3.0 cm, Solute Transport.
|
Test3
|
Dual-Permeability Model, Infiltration from a Surface Ring into isotropic layered system, Solute Transport.
|
Test3a
|
Dual-Permeability Model, Infiltration from a Surface Ring into anisotropic layered system (different anisotropy in Fracture (10:1) and Matrix (1:1) domains), Solute Transport.
|
Test3b
|
Dual-Permeability Model, Infiltration from a Surface Ring into anisotropic layered system (different anisotropy in Fracture (1:10) and Matrix (1:1) domains), Solute Transport.
|
References:
Gärdenäs, A., J. Šimunek, N. Jarvis, and M. Th. van Genuchten, Two-dimensional modelling of preferential water flow and pesticide transport from a tile-drained field, J. Hydrology, 329, 647-660, 2006.
Gerke, H. H., & van Genuchten, M. Th., 1993. A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media. Water Resour. Res. 29, 305-319.
Šimůnek, J., N. J. Jarvis, M. Th. van Genuchten, and A. Gärdenäs, Review and comparison of models for describing non-equilibrium and preferential flow and transport in the vadose zone, Journal of Hydrology, 272, 14-35, 2003.
Šimůnek, J., M. Th. van Genuchten, and M. Šejna, The HYDRUS Software Package for Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media, Technical Manual, Version 1.0, PC Progress, Prague, Czech Republic, pp. 241, 2006.
van Genuchten, M. Th., and J. Šimůnek, Integrated modeling of vadose zone flow and transport processes, Proc. Unsaturated Zone Modelling: Progress, Challenges and Applications, Eds. R. A. Feddes, G. H. de Rooij, and J. C. van Dam, Wageningen UR Frontis Series, Vol. 6, Chapter 2, pp. 37- 69, x-xi, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2004.
Šimůnek, J., M. Th. van Genuchten, and M. Šejna, Modeling Subsurface Water Flow and Solute Transport with HYDRUS and Related Numerical Software Packages, In: Garcia-Navarro & Playán (eds.), Numerical Modelling of Hydrodynamics for Water Resources, An International Workshop, Centro Politecnico Superior, University of Zaragoza Spain, June 18-21 2007. Taylor & Francis Group, London, ISBN 978-0-415-44056-1, 95-114, 2007.
Köhne, J. M., S. Köhne, and J. Šimůnek, A review of model applications for structured soils: a) Water flow and tracer transport, J. Contam. Hydrology, Special Issue “Flow Domains”, doi:10.1016/j.jconhyd.2008.10.002, 104(1-4), 4-35, 2009.
Köhne, J. M., S. Köhne, and J. Šimůnek, A review of model applications for structured soils: b) Pesticide transport, J. Contam. Hydrology, Special Issue “Flow Domains”, doi:10.1016/j.jconhyd.2008.10.003, 104(1-4), 36-60, 2009.
Šimůnek, J., M. Šejna, and M. Th. van Genuchten, The DualPerm Module for HYDRUS (2D/3D) Simulating Two-Dimensional Water Movement and Solute Transport in Dual-Permeability Porous Media, Version 1.0, PC Progress, Prague, Czech Republic, 32 pp., 2012.
Check out the corresponding 1D Problems.