The permeability of two-dimentional porous medium of square fibers (cell model)

E. V. Mosina ; I. V. Chernyshev

Matematičeskaâ fizika i kompʹûternoe modelirovanie, no. 2 (2017), pp. 56-64.

Résumé

In the present study we have solved numerically the problem of a creeping two-dimensional viscous incompressible flow through a fibrous porous medium. The medium is modeled by a system of square fibers (rods) arranged across the flow. The viscous flow around the fibers is solved under the Stokes approximation using the cell model. The Happel boundary conditions (zero normal velocity and tangential shear stresses) have been posted on the boundary of the fluid cell. The numerical realization is based on finite differences of the second order accuracy with successive over-relaxation method for solving a system of linear equations. For regular lattice of square cylinders the hydrodynamic microscopic velocity fields were determined also by the direct numerical simulations. The permeability coefficient of the model porous medium is evaluated by macroscopic averaging and Darcy low assumption. It was found that the cell model with the Happel boundary conditions gives higher values of permeability in comparison with direct numerical calculations. A comparison between these results obtained by cell model and other direct numerical calculations, empirical and analytical relationships in a wide range of the solid concentration $\phi$ made it possible to make the following inferences. The permeability of a rarefied system of square cylinders can be described by an explicit expression that was obtained for circular cylinders [11] but with form factor $\sigma^2$ that reflects the difference between the square and circular cross sections of the rods $$\frac{k}{a^2}=\frac{1}{\sigma^2}\frac{1}{8\phi}\left(\ln\frac{1}{\phi}-1.476+2\phi\right),\quad \phi\leq 0.2.$$ For moderately rarefied porous media the semi-empirical Kozeny — Carman formula is appropriate. In this case corrective factor $1/\sigma^2$ gives error of the permeability more than 25 %. To decrease the disagreement of the numerical values with the Kozeny — Carman formula this relation is better to use with the factor 8/9. The specified expression for permeability is $$\frac{k}{a^2}=\frac{1}{9\phi}\left(\ln\frac{1}{\phi}-\frac{1-\phi^2}{1+\phi^2}\right),\quad 0.1\phi0.6.$$ At higher fiber volume fraction the Kozeny coefficient is $c_k=6.0$, which corresponds to a large number of experimental data on the real porous beds consisting of particles of various forms [7] $$\frac{k}{a^2}=\frac{(1-\phi)^3}{24\phi^2},\quad 0.6\leqslant \phi0.8.$$ These formulas are good approximations for the permeability coefficient for the system of square rods. Obtained relations for the permeability can be useful in the solutions of the macroscopic filtration problems with fibrous porous media.

Détail

Keywords: model porous medium, square fibers, Stokes approximation, permeability.
Mots-clés : Happel cell model

Bibliographie

[1] D. Anderson, J. Tannehill, R. Pletcher, Computational Fluid Mechanics and Heat Transfer, v. 2, Mir Publ., Moscow, 1990, 392 pp. | MR

[2] P.\;S. Kirdan, The Coefficient of Permeability of Two Dimentional Array of Rigit Particles in Slow Parallel Flow: Master Thesis for Degree on Applied Mathematics and Computer Science, Volgograd, 2012, 27 pp.

[3] L.\;S. Leybenzon, The Motion of Nature Fluid and Gas Through Porous Medium, Gostekhizdat Publ., Moscow, Leningrad, 1947, 245 pp.

[4] E.\;V. Mosina, I.\;V. Chernyshev, “Fluid Flow Near the Porous Boundary”, Vestnik Nizhegorodskogo universiteta im. N.I. Lobachevskogo, 4:3 (2011), 999–1001

[5] E.\;V. Mosina, I.\;V. Chernyshev, “Slip Condition on the Surface of a Model Fibrous Porous Medium”, Technical Physics Letters, 35:5 (2009), 103–110

[6] E.\;V. Mosina, “Numerical Study of Flow at a Liquid — Porous Medium Interface”, Theoretical Foundations of Chemical Engineering, 44:5 (2010), 536–542

[7] J. Happel, H. Brenner, Low Reynolds Numbers Hydrodynamics, Mir Publ., Moscow, 1976, 631 pp.

[8] A.\;C. Baytaş, D. Erdemş, H. Acar, O. Çetiner, H. Başci, “Analytical determination of the permeability for slow flow past periodic arrays of cylinders with different cross sections”, J. Porous Med., 15 (2012), 1009–1018 | DOI

[9] J.\;E. Drummond, M.\;I. Tahir, “Laminar viscous flow through regular arrays of parallel solid cylinders”, Int. J. Multiphase Flow, 10 (1984), 515–540 | DOI | Zbl

[10] M. Hellou, J. Martinez, M.\;El Yazidi, “Stokes flow through microstructural model of fibrous media”, Mech. Res. Commun., 31 (2004), 97–103 | DOI | Zbl

[11] D.\;F. James, A.\;M.\;J. Davis, “Flow at the interface of a model fibrous porous medium”, J. Fluid Mech., 426 (2001), 47–72 | DOI | MR | Zbl

[12] M.\;A. Tahir, H.\;V. Tafreshi, “Influence of fiber orientation on the transverse permeability of porous media”, Phys. Fluids, 21 (2009) | DOI | Zbl

[13] A. Tamayol, M. Bahrami, “Transverse permeability of fibrous porous medi”, Phys. Review, 83 (E) (2011), 046314-1-046314-9 | DOI

[14] C.\;Y. Wang, “Stokes flow through an array of rectangular fibers”, Int. J. Multiphase Flow, 22 (1996), 185–194 | DOI | Zbl

[15] P.\;K. Yada, “Slow motion of a porous cylindrical shell in a concentric cylindrical cavity”, Meccanica, 48 (2013), 1607–1622 | DOI | MR | Zbl