The Permeability of Boolean Sets of Cylinders
La perméabilité de modèles Booléens de cylindres
Mines ParisTech, PSL Research University, Center for mathematical morphology, 35 rue St Honoré, 77300
Fontainebleau – France
email: francois.willot@ensmp.fr
^{*} Corresponding author
Numerical and analytical results on the permeability of Boolean models of randomlyoriented cylinders with circular crosssection are reported. The present work investigates cylinders of prolate (highlyelongated) and oblate (nearly flat) types. The fluid flows either inside or outside of the cylinders. The Stokes flow is solved using fullfields Fourierbased computations on 3D binarized microstructures. The permeability is given for varying volume fractions of pores. A new upperbound is derived for the permeability of the Boolean model of oblate cylinders. The behavior of the permeability in the dilute limit is discussed.
Résumé
Ce travail numérique et analytique porte sur la perméabilité de milieux Booléens de cylindres à section circulaire. On considère les deux cas de cylindres très allongés ou très aplatis. L’écoulement du fluide se fait à l’extérieur ou encore à l’intérieur des cylindres. L’équation de Stokes est résolue par méthode de Fourier sur des microstructures binarisées et les perméabilités calculées. Le comportement de la perméabilité d’un modèle Booléen de cylindre au voisinage de la limite diluée est étudié au moyen de bornes analytiques et les résultats discutés.
© F. Willot et al., published by IFP Energies nouvelles, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Introduction
In the recent years, several socalled “Fourierbased” algorithms have been introduced to solve the problem of Stokes flow in heterogeneous media [1–3]. These promising methods, which rely on digitized images of the microstructures and do not require meshing, are especially useful for predicting the permeability of materials based on microstructure models or segmented microtomography images. They have been applied to compute the permeability of textiles [4, 5], ceramic foams [6] or anode layers in fuel cells [7], making use of 3D models.
Numerical computations on largescale 3D microstructures are also useful for investigating theoretical and empirical homogenization estimates, like the classical CarmanKozeny [8, 9] formula, Doi’s upperbound [10, 11] based on correlation functions or dilute limit expansions [12]. The Boolean model of spheres has been considered in the literature [13] and, as expected, follows the CarmanKozeny approximation [14], whereas the upperbound of Doi is valid in the dilute limit only, i.e. for a small volume fraction of obstacles.
The permeability of ordered arrays of cylinders has received considerable interest in experimental, theoretical and numerical works [15–19]. Attention has also been devoted in the literature to various random fibrous media, where fluid flows around cylindrical obstacles. Many of the theoretical homogenization methods used for regular arrays, however, such as the construction of trial fields, cannot be employed for random media. To study random media, numerical methods, such as the LatticeBoltzmann method, and comparisons with analytical formula for regular arrays have been used. In one of the first such study [20], an analytical estimate, based on numerical predictions, is proposed for disordered fibrous media. The estimate has been found to be accurate for porosities less than about 25% [21]. A fiberweb material, with a structure similar to that found in fibrous sheets or paper, has been investigated numerically in [22]. Extensive numerical computations have been performed in [23] on models of prolate spheroids with varying aspect ratio. More recently, the efficient Lattice Boltzmann method has been used [24] to compute the permeability fibrous media, where the fluid flows outside of the cylinders. In one such study [25], several models made of random cylinders are investigated. An empirical law is proposed for the permeability, based on the theoretical analysis of the permeability of regular arrays of cylinders [26]. The effect of the cylinders curvature is found to be small compared to straight cylinders. Computations have also been carried out with varying cylinders aspect ratios. The combined numerical and experimental study in [5] considers fibrous media. The authors investigate in particular the domain of validity of the relevant analytical estimates. In another experimental investigation, Darcy and DarcyForchheimer laws have been identified on packings of nonspherical particles, including prisms and cylinders [27].
Much less work has been devoted to the “reverse” model of fluid flow through fibers, despite some early works [28]. Systems of random flat cylinders have been investigated in hydrogeological literature on fractured rocks. The permeability of 3D random models of fracture networks has been investigated numerically and theoretically in various works [29, 30]. Among important geometrical factors that influence the permeability, authors have highlighted the role of contact regions [31] or that of a wide size distribution in fractures [32].
In this work, a systematic numerical investigation of the Boolean model of cylinders, is carried out with emphasis on effective permeability, fluid velocity distributions, and analytical bounds. The study is organized as follows. The microstructure model is introduced in Section 1.1, the equations to solve for Stokes flow are given in Section 1.2. Fast Fourier Transform (FFT) numerical results are presented in Section 2 with permeability of the various models and numerical maps. Analytical upperbounds are derived in Section 3 for flows in Boolean models of cylinders in the dilute limit. Concluding remarks are given in Conclusion.
1 Permeability of random media
1.1 Boolean Model of Cylinders
The microstructures considered hereafter consist in Boolean sets of cylinders. Boolean models [33] are a classical type of random models which depend on a primary grain and on a Poisson point process [34]. In the present work, the primary grain is defined by a distribution of cylinders, and the homogeneous Poisson point process by a spatiallyconstant intensity n, the mean number of points per unit volume. A cylinder is attached to each point of the Poisson point process, and the Boolean set is the union of cylinders. In the following, we choose to focus on cylinders with high aspect ratio that significantly depart from sphere packings. For such ideal materials, we expect in particular a strong influence of the microstructure on the permeability.
The cylinders have circular crosssection of radius r and height h but their main axis is randomlyoriented. The direction of the axis follows a uniform distribution on the sphere, resulting in (macroscopically) isotropic Boolean sets. Cylinder intersections are allowed so that the volume fraction of cylinder p varies from 0 to 1 and is given by [33]:(1)
where V_{0} = πr^{2}h is the volume of a cylinder. The Boolean sets accordingly depend on n, h and r, or, equivalently q, h and r. In the following, cylinders with very high or small aspect ratio h/r are considered, with p fixed, as detailed below. All microstructures accordingly depend on two parameters, the cylinders volume fraction q and a characteristic size for the cylinder which is either r or h.
Let first r ≪ h and consider a model of prolate cylinders. Two configurations are possible: “porous cylinders” (model A): the fluid flows inside cylinders; or “cylindrical obstacles” (model B): the fluid flows in the complementary set of the cylinders. Let f be the volume of the phase where the fluid flows, so that f = p for model A and f = 1 − p for model B. Macroscopic flow is permitted if the cylinders percolate (A) or if the complementary set of the cylinders percolates (B). The range of useful values for f therefore reads f_{c}^{A} ≤ f ≤ 1 and f_{c}^{B} ≤ f ≤ 1 where f_{c}^{A}, f^{B} are the two percolation thresholds for the Boolean model of cylinders. Note however that f_{c}^{A} → 0 when r/h → 0 as demonstrated for spheroids [35] and ellipsoids [36]. Likewise, the cylinders become infinitesimal lines as r/h → 0 and one expects f_{c}^{B} → 0.
Secondly, let h ≪ r and consider oblate cylinders. The configuration of interest is that of “porous cylinders” (C) where the fluid flows inside cylinders. In the case h ≪ r, the cylinders have asymptotically the same shape as sections enclosed by two (randomlyoriented) parallel planes. The percolation threshold f_{c}^{C} again tends to 0 (we refer to [37] for a study on the connectivity of such systems). Accordingly, f is defined in the entire interval [f_{c}^{C}; 1] ≈ [0; 1]. Note that, for r = ∞, the “reverse” model where fluids flow outside the cylinders presents no interest as the complementary set does not percolate.
Hereafter, all microstructures for models A, B and C are generated on grids of L^{3} voxels with L = 512. The cylinders highest dimension is set to 512 voxels in all cases so that h = L for models A, B and r = L/2 for model C. In configurations A and B, we take r = 3 voxels, to keep a large aspect ratio with a crosssection of the cylinders that resembles disks. By comparison, radii lengths of 6 to 12 lattice units were chosen in media containing 100^{3} to 300^{3} elements in [21].
Likewise, the value of h for model C cannot be as small as 1 due to discretization effects in the numerical method. We take h = 4 voxels for model C. The aspect ratio is accordingly h/r = 171 (A, B) and r/h = 64 (C). Sections of the microstructures used for models A, B and C are shown in Figure 1.
Figure 1. 2D sections of Boolean models of cylinders (black) of the prolate a) and oblate b) types. a) Both phases percolate. Fluid flow occurs inside the black phase (model A) or inside the white phase (model B). b) The black phase only percolates. Fluid flow occurs inside the black phase (model C). 
1.2 Stokes Flow
We consider a viscous fluid of velocity u satisfying the Stokes equation in the pores:(2)
where p(x) is the pressure at point x and μ the fluid viscosity. The fluid velocity is zero (u = 0) along the interface with the solid. Periodic boundary conditions are applied on the computational domain, of the form [38]:(3)
here # denotes periodicity and ΔP is the macroscopic pressure drop, oriented along the first axis e_{1} of a Cartesian coordinate system. The permeability k is a scalar for isotropic media. It is given by Darcy’s law:(4)
The above problem is solved using fullfield computations with the FFTStokes algorithm [3]. The number of iterations used to solve each problem ranges from 15,000 to about 35,000, for an error criterion of η = 10^{−6} (we refer to [3] for details on the FFT algorithm). The convergence of η toward 0 is roughly a powerlaw of the number of iterations (Fig. 2). In many instances, however, convergence is either faster or slower during the computations. Each computation took about 2 days on a 16cores machine with a CPUparallelized code.
Figure 2. Error criterion of the MINRES algorithm used in the FFT scheme for Stokes flow computation in model C, as a function of the number of iterations, for various values of the porosity. 
2 FFT results: permeability
The results are compared to the CarmanKozeny estimate:(5)
where c = 5 is an empirical factor and γ is the specific surface area. For Boolean models A, B and C, we have [39]:(6)
where S_{0}/V_{0} is the surface/volume ratio of a cylinder.
FFT results for the permeability of models A and B are plotted in Figure 3. The data pertaining to model B (flow outside cylinders) is qualitatively similar to the results of LatticeBoltzmann computations [21] performed for cylinders with aspect ratio h/r ≈ 29, i.e. 6 times smaller than the one used here. Specifically, the authors distinguish three “concentrated” (f ≤ 0.5), “intermediate” (0.5 ≤ f ≤ 0.8) and “semidilute” (0.8 ≤ f) regimes. The intermediate regime corresponds to an exponential scaling law of k with respect to the porosity k ~ − log f. Our data presents a similar scaling law in the region 0.4 ≤ f ≤ 0.8 with a slower decrease of k as a function of f in the concentrated regime (f ≤ 0.3). The FFT results for model B also compare well to the CarmanKozeny estimates (Eq. 5) and (Eq. 6) in a quite wide region 0.3 ≤ f ≤ 0.9 (solid black curve). As expected, the CarmanKozeny estimate is not valid in the dilute and concentrated regimes, in agreement with [21]. For model B, we also make use of the JacksonJames estimate [20], valid for fibrous media in the dilute region f ≈ 1 (solid green curve).
Similar results hold for models A and C (Fig. 3 and 4). Again, FFT data points are close to the CarmanKozeny estimates (Eq. 5) and (Eq. 6). This is especially true for model C, except in two narrow regions near the concentrated and highlydilute regimes.
Figure 3. Effective permeability k vs porosity for a Stokes flow inside and outside a Boolean model of highlyelongated cylinders. Symbols: FFT results, solid lines: CarmanKozeny estimates. 
Figure 4. Effective permeability k vs porosity for a Stokes flow inside a Boolean model of highlyflate cylinders. Symbols: FFT results, solid line: CarmanKozeny estimate. 
Histograms of the fluid velocity component u_{1} are represented in Figure 5 for the various models. Remarkably, the shape of the distribution takes different forms with respect to the models. This is especially true of models A vs B. The field maps show a very high concentration of the fluid velocity field around selected points, for all models and all values of the porosity, except in dilute regimes. Maps of the velocity field at selected porosity are shown in Figure 6.
Figure 5. Histograms of the fluid velocity component u_{1} for models A a), B b) and C c). The curves displayed in each plot correspond to different porosity values with wide histograms centered on the right for high porosity and narrow histograms at the lowest porosity values. 
Figure 6. 2D sections of the fluid velocity field u_{1} along a plane normal to e_{3} for models A a), B b) and C c). The porosity is respectively f = 50%, 43% and 14%. Highest field velocity in yellow, lowest in black, with intermediary values in blue and green. 
3 Upperbound for the permeability of Boolean sets of oblate cylinders
This section is devoted to upperbounds for the permeability of a Boolean model of cylinders of the “flat” type, relevant to model C. Recall that flow occurs inside the cylinders in model C. The results in this section are based on exact results previously obtained for the covariogram of cylinders [40]. Microstructures with elongatedtype cylinders (h ≫ r), relevant to models A and B, have much more complex covariance functions and are not treated here.
The Doi [10] upperbound for permeability reads:(7)
where s is the specific surface area of the voidsolid interface, and F_{vv}, F_{sv} and F_{ss} are the “volumevolume”, “surfacevolume” and “surfacesurface” correlation functions:(8) (9) (10)
In the above, Z is the indicator function for the complementary B^{c} of set B (the obstacles), ∇ is the gradient operator and u is a unit vector (u = 1). The quantity F_{vv} is given by the probability that two points separated by a distance t lie in B^{c}. The quantities F_{sv} and F_{ss} can also be interpreted in terms of probability, replacing the interface of B by an interphase of vanishingly small width δ. The quantity F_{sv} is proportional to the probability (divided by δ) that one point lies in B^{c} and the other in the interphase. The quantity F_{ss} is related to the probability that the two points lie in the interphase. Equation (7) is a rigorous bound which provides the exact lowestorder correction to the permeability for dilute concentration of spherical obstacles q → 0 [12].
The correlation function F_{vv} of a Boolean model of cylinders is given by [39]:(11)
where K(t; h, r) (not to be confounded with k) is the normalized mean cylinder covariogram. The latter is defined as the volume intersected by two cylinders as follows:(12)
where C is the volume of the cylinder C and C_{t} is C translated by a vector of length t. The quantity is averaged over orientations uniformly distributed on the sphere and normalized by the cylinder volume, so that k(t = 0; h, r) ≡ 1.
3.1 Case r = ∞
We first consider the case of infinitely flat cylinders with r = ∞ [40]. Note that h and q are kept finite in this limit. Accordingly, the Poisson intensity θ (mean number of cylinders per unit of volume) should tend to 0 in a manner inversely proportional to the cylinder volume. Indeed θ is related to q by [39]:(13)
and so θ ∼ 1/r^{2} → 0. Accordingly, the microstructure should be understood as a set of randomlyoriented, very large cylinders, with a very low density (in number per unit volume), so that the cylinders volume fraction remains finite. This limit microstructure could be usefuly compared to Poisson tessellations (or dilated versions thereof). We refer to [41, 42] for theoretical results on the latter.
The equation for F_{sv} and F_{ss} in Boolean media are similar to that of F_{vv} [43, 44] but quite more complex. The formula involve first and secondorder derivatives at h′ = h of the following functional:(15)
where, C′ is a thick plane of width h′. Straightforward computations lead to:(16)
The three expressions for the correlation functions allow one to compute the integrand in Equation (7):(17)
and so the integral in Equation (7) diverges and the upperbound is infinite. The bound accordingly provides no information on the permeability of materials with “infinitely” flatcylinders. This negative result will nevertheless prove useful for interpreting the results of Section 3.2.
3.2 Case r Finite
We now seek for the computation of upperbounds on the permeability for cylinders of the oblate type but with finite radius r. The isotropized covariogram K of general cylinders with finite r and h was studied in [40, 45]. The exact form provided in [40] is a rather complicated analytical expression involving incomplete Elliptic functions. In the present work, we use the following approximate expression [40]:(18a) (18b)
In the above, cos^{−1} is the inverse cosine (also denoted arccos) function. For r > h, (Eq. 18) is a very good approximation of K(t) which becomes asymptotically exact when r/h → ∞. More precisely, the maximal difference ϵ = sup_{t} K(t) − K_{app}(t) is attained when t ⪆ h for all r > h. For fixed h, the maximal error ϵ is a decreasing function of r, and ϵ → 0 when r → ∞. As an indication, we have ϵ ≈ 1.1% for r = h and ϵ ≈ 0.1% for r = 4h (see [40] for details).
The derivation of Doi’s bounds for general cylinders requires one to compute the volume of the intersection of two cylinders of different heights and radii averaged over uniformlydistributed directions, which we do not have. Accordingly, hereafter we use the simpler BerrymanMilton bound [12, 46] which reads:(19)
For cylinders with r > h, replacing F_{vv} with Equation (11) and using Equation (18) yields, after integrating:(20)
The upperbound is only useful in the dilute limit. In fact, the latter largely overestimates the FFT data for the permeability and represented in Figure 4 (not shown on the plot). Its interest lies in the dilute regime q → 0. Fix r and expand the above to zeroorder correction O(1) in the dilute limit q → 0 and let r → ∞ afterwards. We obtain:(21)
At the lowestorder correction, one recovers an asymptotic regime ∼ 1/[q(log q)^{2}], a higher correction than that predicted by the JamesJackson model [20] which gives ∼ −1/[q log q] as firstorder correction. The latter is itself higher than the dilute expansion for spheres ∼1/q.
Interestingly, the second leadingorder term in Equation (21) scales as . A similar nonanalytic dependence appears in the dilute limit expansion of the permeability of a bed of totally impenetrable spheres, obtained by selfconsistent methods [47]:(22)
where a is the sphere radius. The term is related to hydrodynamic screening effects [48, 49]. For a dilute packing of spheres, as shown in [12], no screening term is predicted by the bound Equation (19) or by Doi’s bound [10]. As shown in expansion in Equation (21), a screening term is indeed predicted by the variational bound, for the particular highlyelongated shapes considered here. Note that in the dilute limit considered here, the obstacles are not cylinders, but have a variety of shapes with very large surface/volume ratio.
Although the bound in Equation (19)is in general less accurate than Doi’s bound which uses the surfacesurface correlation information [12] both are relevant in the dilute limit. Doi’s bound predicts the correct leadingorder term in the dilute limit for Boolean spherical obstacles. The bound in Equation (19) predicts the correct scaling law ∼ a^{2}/q, with a slightly overestimated prefactor 4/15 ≈ 0.27 instead of 2/9 ≈ 0.22 for the Boolean model of spheres.
Compute now the limit r → ∞ and afterwards q → 0 in Equation (20):(23)
The expansion above depends on r unlike Equation (22) and tends to +∞ which explains the result obtained previously for Doi’s bound with r = ∞. When taking the limits r → ∞ and q → 0, different regimes are obtained depending on the order one takes the limits. The two regimes appear when the permeability is plotted as a function of q in loglog plot (Fig. 7). The regime change occurs at some points r ≈ r_{c}(q), q ≈ q_{c}(r) which are obtained by equaling the two highestorder corrections Equation (21) and Equation (23):(24)
Figure 7. BerrymanMilton upperbound (Eq. 19) for the permeability k of a Boolean model of flat cylinders with height h = 1 and radius r = 10^{2}, 10^{3}, 10^{4}vs concentration of obstacles q, in loglog plot (black). Solid lines in green: r → ∞ expansion (Eq. 23). Solid lines in pink: dilute limit expansion q → 0 (Eq. 22). 
When r ≫ r_{c}(q), or equivalently q =≪ q_{c}(r), the permeability follows expansion (Eq. 23). Expansion (Eq. 22) holds in the domain r ≪ r_{c}(q) or q =≫ q_{c}(r).
Comparison with FFT data would be needed to confirm the existence of the two regimes, however, accurate computations in the dilute limit requires high computational costs [14]. This is left to future work.
Conclusion
In this work, the permeability of various Boolean models of cylinders with very high aspect ratio have been computed, with fluid flow inside or outside cylinders. The flat and elongated cylinders are tantamount to models of planes or lines of finite width. Good agreement is found with the CarmanKozeny approximation in the intermediate (nondilute, nonconcentrated) regime. The estimate is particularly relevant for the model of plane or flat cylinders.
The behavior of certain upperbound in the dilute limit has been computed for the model of flat cylinders in the dilute limit, making use of exact formula for the covariogram of cylinders. The upperbound predicts two distinct regimes depending on the diameter of the cylinders. In the stronglydilute regime, a scaling law ∼ 1/[q(log q)^{2}] is predicted, whereas a much weaker regime ∼ log q arises when r ≫ h/q. Numerical computations, which are quite challenging with the methods currently available, should be carried out to confirm or infirm the predicted behavior.
Acknowledgments
This project has received funding from the European Union’s Seventh Framework Programme (FP7/20072013) for the Fuel Cells and Hydrogen Joint Technology Initiative under grant agreement No. 303429.
References
 Bignonnet F., Dormieux L. (2014) FFTbased bounds on the permeability of complex microstructures, International Journal of Numerical and Analytical Methods in Geomechanics 38, 16, 1707–1723. [Google Scholar]
 Monchiet V., Bonnet G., Lauriat G. (2009) A FFTbased method to compute the permeability induced by a Stokes slip flow through a porous medium, Comptes Rendus Mécanique 337, 4, 192–197. [Google Scholar]
 Wiegman A. (2007) Computation of the permeability of porous materials from their microstructure by FFTStokes, Report of the Fraunhofer ITWM 129, http://kluedo.ub.unikl.de/files/1984/bericht129.pdf, accessed 22 July 2015. [Google Scholar]
 Gervais P.C., BardinMonnier N., Thomas D. (2012) Permeability modeling of fibrous media with bimodal fiber size distribution, Chemical Engineering Science 73, 239–248. [CrossRef] [Google Scholar]
 Rief S., Latz A., Wiegmann A. (2006) Computer simulation of air filtration including electric surface charges in threedimensional fibrous micro structures, Filtration 6, 2, 169–172. [Google Scholar]
 Redenbach C., Wirjadi O., Rief S., Wiegmann A. (2011) Modeling a ceramic foam for filtration simulation, Advanced Engineering Materials 13, 3, 171–177. [CrossRef] [Google Scholar]
 Masson D., Abdallah B., Willot F., Jeulin D., Mercadelli E., Sanson A., Chesnaud A., Thorel A. (2015) Morphological modeling of a metal foam SOFC configuration, ECS Transactions 68, 1, 2951–2960. [CrossRef] [Google Scholar]
 Carman P.C. (1937) Fluid flow through granular beds, TransactionsInstitution of Chemical Engineers 15, 150–166. [Google Scholar]
 Kozeny J. (1927) Ueber kapillare leitung des wassers im boden, Sitzungsber Akad Wiss Wien 136, 271–306. [Google Scholar]
 Doi M. (1976) A new variational approach to the diffusion and the flow problem in porous media, Journal of the Physical Society of Japan 40, 2, 567–572. [CrossRef] [Google Scholar]
 Torquato S., Lu B. (1990) Rigorous bounds on the fluid permeability: Effect of polydispersivity in grain size, Physics of Fluids A: Fluid Dynamics 2, 4, 487–490. [CrossRef] [Google Scholar]
 Rubinstein J., Torquato S. (1989) Flow in random porous media: mathematical formulation, variational principles, and rigorous bounds, Journal of Fluid Mechanics 206, 25–46. [CrossRef] [Google Scholar]
 Martys N.S., Torquato S., Bentz D.P. (1994) Universal scaling of fluid permeability for sphere packings, Physical Review E 50, 1, 403. [CrossRef] [Google Scholar]
 Abdallah B., Willot F., Jeulin D. (2015) Stokes flow through a Boolean model of spheres: Representative volume element, Transport in Porous Media 109, 3, 711–726. [CrossRef] [Google Scholar]
 Martin J.J., McCabe L.W., Monrad C.C. (1951) Pressure drop through stacked spheres. Effect of orientation, Chemical Engineering Progress 47, 2, 91–94. [Google Scholar]
 Richardson F.J., Zaki W.N. (1954) The sedimentation of a suspension of uniform spheres under conditions of viscous flow, Chemical Engineering Science 3, 2, 65–73. [CrossRef] [Google Scholar]
 Shahmardan M.M., Nazari M., Khaksar M., Khatib M. (2013) A new mathematical model for permeability of composites, Journal of Solid Mechanics 5, 4, 371–379. [Google Scholar]
 Tamayol A., Bahrami M. (2009) Analytical determination of viscous permeability of fibrous porous media, International Journal of Heat and Mass Transfer 52, 9, 2407–2414. [CrossRef] [Google Scholar]
 Yazdchi K., Srivastava S., Luding S. (2011) Microstructural effects on the permeability of periodic fibrous porous media, International Journal of Multiphase Flow 37, 8, 956–966. [CrossRef] [Google Scholar]
 Jackson G.W., James D.F. (1986) The permeability of fibrous porous media, The Canadian Journal of Chemical Engineering 64, 3, 364–374. [CrossRef] [Google Scholar]
 Clague D.S., Kandhai B.D., Zhang R., Sloot P.M.A. (2000) Hydraulic permeability of (un)bounded fibrous media using the lattice Boltzmann method, Physical Review E 61, 1, 616. [CrossRef] [Google Scholar]
 Koponen A., Kandhai D., Hellen E., Alava M., Hoekstra A., Kataja M., Niskanen K., Sloot P., Timonen J. (1998) Permeability of threedimensional random fiber webs, Physical Review Letters 80, 4, 716. [CrossRef] [Google Scholar]
 Claeys I.L., Brady J.F. (1993) Suspensions of prolate spheroids in Stokes flow. Part 2. Statistically homogeneous dispersions, Journal of Fluid Mechanics 251, 443–477. [CrossRef] [Google Scholar]
 Thomas M.L.R., Ingham D.B., Pourkashanian M. (2010) Prediction of the permeability of fibrous porous media using the lattice Boltzmann method in conjuction with coarse numerical lattices, Open Transport Phenomena Journal 2, 80–89. [Google Scholar]
 Nabovati A., Llewellin E.W., Sousa A.C.M. (2009) A general model for the permeability of fibrous porous media based on fluid flow simulations using the lattice Boltzmann method, Composites Part A: Applied Science and Manufacturing 40, 6, 860–869. [CrossRef] [Google Scholar]
 Gebart B.R. (1992) Permeability of unidirectional reinforcements for RTM, Journal of Composite Materials 26, 8, 1100–1133. [Google Scholar]
 Clavier R., Chikhi N., Fichot F., Quintard M. (2015) Experimental investigation on singlephase pressure losses in nuclear debris beds: Identification of flow regimes and effective diameter, Nuclear Engineering and Design 292, 222–236. [CrossRef] [Google Scholar]
 Matheron G. (1979) L’émergence de la loi de Darcy, École Nationale Supérieure des Mines, Note Géostatistique No. 592, www.cg.ensmp.fr/bibliotheque/public/MATHERON_Rapport_00216.pdf, accessed August 21, 2015. [Google Scholar]
 Koudina N., Garcia R.G., Thovert J.F., Adler P.M. (1998) Permeability of threedimensional fracture networks, Physical Review E 57, 4, 4466. [CrossRef] [Google Scholar]
 Lang P.S., Paluszny A., Zimmerman R.W. (2014) Permeability tensor of threedimensional fractured porous rock and a comparison to trace map predictions, Journal of Geophysical Research: Solid Earth 119, 8, 6288–6307. [CrossRef] [Google Scholar]
 Berkowitz B. (2002) Characterizing flow and transport in fractured geological media: A review, Advances in Water Resources 25, 8–12, 861–884. [Google Scholar]
 Mourzenko V.V., Thovert J.F., Adler P.M. (2004) Macroscopic permeability of threedimensional fracture networks with powerlaw size distribution, Physical Review E 69, 6, 066307. [CrossRef] [Google Scholar]
 Matheron G. (1972) Random sets theory and its applications to stereology, Journal of Microscopy 95, 1, 15–23. [CrossRef] [Google Scholar]
 Miles R.E. (1970) On the homogeneous planar Poisson point process, Mathematical Biosciences 6, 85–127. [CrossRef] [MathSciNet] [Google Scholar]
 Jeulin D., Moreaud M. (2011) Percolation of random cylinder aggregates, Image Analysis & Stereology 26, 3, 121–127. [CrossRef] [Google Scholar]
 Garboczi E.J., Snyder K.A., Douglas J.F., Thorpe M.F. (1995) Geometrical percolation threshold of overlapping ellipsoids, Physical Review E 52, 1, 819. [CrossRef] [Google Scholar]
 Huseby O., Thovert J.F., Adler P.M. (1997) Geometry and topology of fracture systems, Journal of Physics A: Mathematical and General 30, 5, 1415. [CrossRef] [MathSciNet] [Google Scholar]
 Ene H.I., SanchezPalencia E. (1975) Équations et phénomènes de surface pour l’écoulement dans un modèle de milieu poreux, Journal de Mécanique 14, 1, 73–108. [Google Scholar]
 Matheron G. (1975) Random sets and integral geometry, Wiley, NewYork. [Google Scholar]
 Willot F. (2016) Mean covariogram of a cylinder and applications to Boolean random sets, Journal of Contemporary Mathematical Analysis, In Press. [Google Scholar]
 Jeulin D. (1991) Modèles de fonctions aléatoires multivariables, Sciences de la Terre 30, 225–256. [Google Scholar]
 Jeulin D. (1991) Modèles Morphologiques de Structures Aléatoires et de Changement d’Échelle, PhD Thesis, Université de Caen. [Google Scholar]
 Rubinstein J., Torquato S. (1988) Diffusioncontrolled reactions: Mathematical formulation, variational principles, and rigorous bounds, The Journal of Chemical Physics 88, 10, 6372–6380. [CrossRef] [Google Scholar]
 Torquato S. (1986) Microstructure characterization and bulk properties of disordered twophase media, Journal of Statistical Physics 45, 5–6, 843–873. [CrossRef] [Google Scholar]
 Gille W. (1987) The intercept length distribution density of a cylinder of revolution, Experimentelle Technik der Physik 35, 2, 93–98. [Google Scholar]
 Berryman J.G., Milton G.W. (1985) Normalization constraint for variational bounds on fluid permeability, The Journal of Chemical Physics 83, 2, 754–760. [CrossRef] [Google Scholar]
 Childress S. (1972) Viscous flow past a random array of spheres, The Journal of Chemical Physics 56, 6, 2527–2539. [CrossRef] [Google Scholar]
 Brady J.F., Bossis G. (1988) Stokesian dynamics, Annual Review of Fluid Mechanics 20, 111–157. [Google Scholar]
 Tam C.K.W. (1969) The drag on a cloud of spherical particles in low Reynolds number flow, Journal of Fluid Mechanics 38, 3, 537–546. [CrossRef] [Google Scholar]
Cite this article as: F. Willot, B. Abdallah and D. Jeulin (2016). The Permeability of Boolean Sets of Cylinders, Oil Gas Sci. Technol 71, 52.
All Figures
Figure 1. 2D sections of Boolean models of cylinders (black) of the prolate a) and oblate b) types. a) Both phases percolate. Fluid flow occurs inside the black phase (model A) or inside the white phase (model B). b) The black phase only percolates. Fluid flow occurs inside the black phase (model C). 

In the text 
Figure 2. Error criterion of the MINRES algorithm used in the FFT scheme for Stokes flow computation in model C, as a function of the number of iterations, for various values of the porosity. 

In the text 
Figure 3. Effective permeability k vs porosity for a Stokes flow inside and outside a Boolean model of highlyelongated cylinders. Symbols: FFT results, solid lines: CarmanKozeny estimates. 

In the text 
Figure 4. Effective permeability k vs porosity for a Stokes flow inside a Boolean model of highlyflate cylinders. Symbols: FFT results, solid line: CarmanKozeny estimate. 

In the text 
Figure 5. Histograms of the fluid velocity component u_{1} for models A a), B b) and C c). The curves displayed in each plot correspond to different porosity values with wide histograms centered on the right for high porosity and narrow histograms at the lowest porosity values. 

In the text 
Figure 6. 2D sections of the fluid velocity field u_{1} along a plane normal to e_{3} for models A a), B b) and C c). The porosity is respectively f = 50%, 43% and 14%. Highest field velocity in yellow, lowest in black, with intermediary values in blue and green. 

In the text 
Figure 7. BerrymanMilton upperbound (Eq. 19) for the permeability k of a Boolean model of flat cylinders with height h = 1 and radius r = 10^{2}, 10^{3}, 10^{4}vs concentration of obstacles q, in loglog plot (black). Solid lines in green: r → ∞ expansion (Eq. 23). Solid lines in pink: dilute limit expansion q → 0 (Eq. 22). 

In the text 