Transient Diffusion and Conduction in Heterogeneous Media: Beyond the Classical Effective-Medium Approximation

Local details versus effective medium approximation: A study of diffusion in microfluidic random networks made from Voronoi tessellations

We measured the effective diffusion coefficient in regions of microfluidic networks of controlled geometry using the fluorescence recovery after photobleaching (FRAP) technique. The geometry of the networks was based on Voronoi tessellations, and had varying characteristic length scale and porosity. For a fixed network, FRAP experiments were performed in regions of increasing size. Our results indicate that the boundary of the bleached region, and in particular the cumulative area of the channels that connect the bleached region to the rest of the network, are important in the measured value of the effective diffusion coefficient. We found that the statistical geometrical variations between different regions of the network decrease with the size of the bleached region as a power law, meaning that the statistical error of effective medium approximations decrease with the size of the studied medium with no characteristic length scale.

Evaluating transport in irregular pore networks

A general approach for investigating transport phenomena in porous media is presented. This approach has the capacity to represent various kinds of irregularity in porous media without the need for excessive detail or computational effort. The overall method combines a generalized effective medium approximation (EMA) with a macroscopic continuum model in order to derive a transport equation with explicit analytical expressions for the transport coefficients. The proposed form of the EMA is an anisotropic and heterogeneous extension of Kirkpatrick's EMA [Rev. Mod. Phys. 45, 574 (1973)] which allows the overall model to account for microscopic alterations in connectivity (with the locations of the pores and the orientation and length of the throat) as well as macroscopic variations in transport properties. A comparison to numerical results for randomly generated networks with different properties is given, indicating the potential for this methodology to handle cases that would pose significant difficulties to many other analytical models.

Multicomponent effective medium-correlated random walk theory for the diffusion of fluid mixtures through porous media

Molecular transport in nanoconfined spaces plays a key role in many emerging technologies for gas separation and storage, as well as in nanofluidics. The infiltration of fluid mixtures into the voids of porous frameworks having complex topologies is common place to these technologies, and optimizing their performance entails developing a deeper understanding of how the flow of these mixtures is affected by the morphology of the pore space, particularly its pore size distribution and pore connectivity. Although several techniques have been developed for the estimation of the effective diffusivity characterizing the transport of single fluids through porous materials, this is not the case for fluid mixtures, where the only alternatives rely on a time-consuming solution of the pore network equations or adaptations of the single fluid theories which are useful for a limited type of systems. In this paper, a hybrid multicomponent effective medium-correlated random walk theory for the calculation of the effective transport coefficients matrix of fluid mixtures diffusing through porous materials is developed. The theory is suitable for those systems in which component fluxes at the single pore level can be related to the potential gradients of the different species through linear flux laws and corresponds to a generalization of the classical single fluid effective medium theory for the analysis of random resistor networks. Comparison with simulation of the diffusion of binary CO(2)/H(2)S and ternary CO(2)/H(2)S/C(3)H(8) gas mixtures in membranes modeled as large networks of randomly oriented pores with both continuous and discrete pore size distributions demonstrates the power of the theory, which was tested using the well-known generalized Maxwell-Stefan model for surface diffusion at the single pore level.

Upscaling the diffusion equations in particulate media made of highly conductive particles. II. Application to fibrous materials

In paper I [Vassal, Phys. Rev. E77, 011302 (2008)] of this contribution, the effective diffusion properties of particulate media with highly conductive particles and particle-particle interfacial barriers have been investigated with the homogenization method with multiple scale asymptotic expansions. Three different macroscopic models have been proposed depending on the quality of contacts between particles. However, depending on the nature and the geometry of particles contained in representative elementary volumes of the considered media, localization problems to be solved to compute the effective conductivity of the two first models can rapidly become cumbersome, time and memory consuming. In this second paper, the above problem is simplified and applied to networks made of slender, wavy and entangled fibers. For these types of media, discrete formulations of localization problems for all macroscopic models can be obtained leading to very efficient numerical calculations. Semianalytical expressions of the effective conductivity tensors are also proposed under simplifying assumptions. The case of straight monodisperse and homogeneously distributed slender fibers with a circular cross section is further explored. Compact semianalytical and analytical estimations are obtained when fiber-fiber contacts are perfect or very poor. Moreover, two discrete element codes have been developed and used to solve localization problems on representative elementary volumes for the same types of contacts. Numerical results underline the significant roles of the fiber content, the orientation of fibers as well as the relative position and orientation of contacting fibers on the effective conductivity tensors. Semianalytical and analytical predictions are discussed and compared with numerical results.

Upscaling the diffusion equations in particulate media made of highly conductive particles. I. Theoretical aspects

Many analytical and numerical works have been devoted to the prediction of macroscopic effective transport properties in particulate media. Usually, structure and properties of macroscopic balance and constitutive equations are stated a priori. In this paper, the upscaling of the transient diffusion equations in concentrated particulate media with possible particle-particle interfacial barriers, highly conductive particles, poorly conductive matrix, and temperature-dependent physical properties is revisited using the homogenization method based on multiple scale asymptotic expansions. This method uses no a priori assumptions on the physics at the macroscale. For the considered physics and microstructures and depending on the order of magnitude of dimensionless Biot and Fourier numbers, it is shown that some situations cannot be homogenized. For other situations, three different macroscopic models are identified, depending on the quality of particle-particle contacts. They are one-phase media, following the standard heat equation and Fourier's law. Calculations of the effective conductivity tensor and heat capacity are proved to be uncoupled. Linear and steady state continuous localization problems must be solved on representative elementary volumes to compute the effective conductivity tensors for the two first models. For the third model, i.e., for highly resistive contacts, the localization problem becomes simpler and discrete whatever the shape of particles. In paper II [Vassal, Phys. Rev. E 77, 011303 (2008)], diffusion through networks of slender, wavy, entangled, and oriented fibers is considered. Discrete localization problems can then be obtained for all models, as well as semianalytical or fully analytical expressions of the corresponding effective conductivity tensors.

Voronoi tessellation of the packing of fine uniform spheres

The packing of uniform fine spherical particles ranging from 1 to 1000 microm has been simulated by means of discrete particle simulation. The packing structure is analyzed, facilitated by the well established Voronoi tessellation. The topological and metric properties of Voronoi polyhedra are quantified as a function of particle size and packing density. The results show that as particle size or packing density decreases, (i) the average face number of Voronoi polyhedra decreases, and the distributions of face number and edge number become broader and more asymmetric; (ii) the average perimeter and area of polyhedra increase, and the distributions of polyhedron surface area and volume become more flat and can be described by the log-normal distribution. The topological and metric properties depicted for the packing of fine particles differ either quantitatively or qualitatively from those reported in the literature although they all can be related to packing density. In particular, our results show that the average sphericity coefficient of Voronoi polyhedra varies with packing density, and although Aboav-Weaire's law is generally applicable, Lewis's law is not valid when packing density is low, which are contrary to the previous findings for other packing systems.