Mass transfer process in a two‐region medium

[1] Two-region mass transfer models have been widely considered for describing flow and transport in highly heterogeneous porous media. Previous studies have indicated that in two-region systems the effective mass transfer coefficient is not a constant but depends upon a time convolution of the concentration history. Although such a representation does capture the mass transfer dynamics correctly, it can be difficult to apply because the convolution operation couples the effective mass transfer rate to the average concentration. In this work, we present the upscaling of a two-region system complete with a method for predicting the time-varying effective mass transfer coefficient that arises. We consider in detail the particular case that occurs when the convective flux in the low conductivity region can be neglected and the inclusions are spherical. This puts the problem in a classical context, and we develop explicit closed-form expressions for the time-varying effective mass transfer coefficient within the context of volume averaging. Our goals in this work are to (1) determine the general (convolution) form of the macroscale two-equation model for transport in such a medium using volume averaging, (2) assess the influence of several simplifications to the convolution form of the two-equation model, and (3) compare this result with existing classical results for the simplified geometry of spherical inclusions. A comparison of the results predicted by theory and experimental data from two macroscale, two-region systems is made. Additional discussion about how the method can be extended to other geometries and to cases where the convective flux is not negligible is also presented.

[1]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[2]  W. G. Gray,et al.  A derivation of the equations for multi-phase transport , 1975 .

[3]  William G. Gray,et al.  On the definition and derivatives of macroscale energy for the description of multiphase systems , 2002 .

[4]  E. Glueckauf,et al.  Theory of chromatography. Part 10.—Formulæ for diffusion into spheres and their application to chromatography , 1955 .

[5]  A. Valocchi,et al.  Constraints on the Validity of Equilibrium and First-Order Kinetic Transport Models in Structured Soils , 1986 .

[6]  R. E. Jessup,et al.  Experimental and mathematical description of nonadsorbed solute transfer by diffusion in spherical aggregates. , 1980 .

[7]  A. G. Greenhill,et al.  Handbook of Mathematical Functions with Formulas, Graphs, , 1971 .

[8]  John C. Slattery,et al.  Flow of viscoelastic fluids through porous media , 1967 .

[9]  Michel Quintard,et al.  Transport in chemically and mechanically heterogeneous porous media V. Two-equation model for solute transport with adsorption , 1998 .

[10]  J. Maxwell.  A Treatise on Electricity and Magnetism , 1873, Nature.

[11]  Charles F. Harvey,et al.  Temporal Moment‐Generating Equations: Modeling Transport and Mass Transfer in Heterogeneous Aquifers , 1995 .

[12]  Michel Quintard,et al.  Comparison of theory and experiment for solute transport in highly heterogeneous porous medium , 2007 .

[13]  S. Whitaker,et al.  Transport in chemically and mechanically heterogeneous porous media—III. Large-scale mechanical equilibrium and the regional form of Darcy's law , 1998 .

[14]  S. Whitaker.  The method of volume averaging , 1998 .

[15]  J. Crank.  XX. A diffusion problem in which the amount of diffusing substance is finite.—II. Diffusion with non-linear adsorption , 1948 .

[16]  C. Moyne.  Two-equation model for a diffusive process in porous media using the volume averaging method with an unsteady-state closure , 1997 .

[17]  S. Whitaker,et al.  Transport in chemically and mechanically heterogeneous porous media IV : large-scale mass equilibrium for solute transport with adsorption , 1998 .

[18]  Michel Quintard,et al.  Volume averaging for determining the effective dispersion tensor: Closure using periodic unit cells and comparison with ensemble averaging , 2003 .

[19]  Mark N. Goltz,et al.  Three‐Dimensional Solutions for Solute Transport in an Infinite Medium With Mobile and Immobile Zones , 1986 .

[20]  Sean Andrew McKenna,et al.  On the late‐time behavior of tracer test breakthrough curves , 2000 .

[21]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[22]  L. Rayleigh,et al.  LVI. On the influence of obstacles arranged in rectangular order upon the properties of a medium , 1892 .

[23]  X. Sanchez‐Vila,et al.  On matrix diffusion: formulations, solution methods and qualitative effects , 1998 .

[24]  Michel Quintard,et al.  Two‐medium description of dispersion in heterogeneous porous media: Calculation of macroscopic properties , 2003 .

[25]  Michel Quintard,et al.  Quasi-steady two-equation models for diffusive transport in fractured porous media: large-scale properties for densely fractured systems , 2001 .

[26]  P. Roberts,et al.  Adsorption of Normal Paraffins from Binary Liquid Solutions by 5A Molecular Sieve Adsorbent , 1967 .

[27]  Michel Quintard,et al.  Two-domain description of solute transport in heterogeneous porous media: Comparison between theoretical predictions and numerical experiments , 2007 .

[28]  P. Carman,et al.  Measurement of diffusion coefficients , 1954, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.