SciELO - Scientific Electronic Library Online

 
vol.88 número217PV system interconnected to the electricity grid with hourly control of energy injectionQuantification of water erosion and characterization of surface sealing in Ultisols in semiarid areas in Brazil índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google

Compartilhar


DYNA

versão impressa ISSN 0012-7353versão On-line ISSN 2346-2183

Dyna rev.fac.nac.minas vol.88 no.217 Medellín abr./jun. 2021  Epub 11-Nov-2021

https://doi.org/10.15446/dyna.v88n217.88713 

Articles

Numerical approximation to the scaling law describing the geometrical tortuosity in porous media

Aproximación numérica a la ley de escalamiento que describe la tortuosidad geométrica en medios porosos desordenados

Alejandro Ramírez-Véleza 
http://orcid.org/0000-0002-9240-3477

a Instituto de Química, Universidad de Antioquia, Medellín, Colombia. alejandro.ramirez@udea.edu.co


Abstract

When studying porous media transport properties, it is crucial to ascertain tortuosity (() and its variation with porosity (ϕ). In this work, numerical methods were used to investigate this relationship. First, a digital representation of media was derived, and thereby implement an algorithm for calculating tortuosity. The program allows deriving several statistics of the paths present within the pores. The results complement the theoretical studies that suggest the existence of a scaling law in disordered media. However, this paper proposes that the relationship between ( and 𝜙 depends on the average fractal dimension instead of the fractal dimensionality of the optimal path. It was also confirmed that the geometry of the latter can be considered in the same universality class as those described by loopless compressible invasion percolation.

Keywords: tortuosity; porous media; percolation; numerical method; scaling law; fractal dimension

Resumen

Cuando se estudian las propiedades de transporte de un medio poroso, es importante conocer la tortuosidad (() y su variación con la porosidad (𝜙). En este trabajo se utilizan métodos numéricos para buscar dicha relación. Primero se obtiene una representación digital de los medios y luego se implementa un algoritmo para el cálculo de la tortuosidad. El programa permite conocer varios estadísticos de los caminos presentes al interior de los poros. Los resultados sirven de complemento a los estudios teóricos que sugieren la existencia de una ley de escalamiento en medios desordenados. Sin embargo, se propone que la relación entre ( y 𝜙 depende de la dimensión fractal promedio en lugar de la dimensión fractal del camino óptimo. También se verifica que la geometría de este último, puede considerarse dentro de la misma clase de universalidad que las descritas por la percolación por invasión compresible sin bucles.

Palabras clave: tortuosidad; medios porosos; percolación; método numérico; ley de escalamiento; dimensión fractal

1. Introduction

When studying the behavior of a porous medium, it is crucial to gain information about its effective transport properties (i.e., conductivity, permeability, permittivity, and diffusivity). For example, knowledge of the soil hydraulic permeability helps predict productivity in an oil well [1] or the pollutant migration toward aquifers or other water bodies [2]. The thermal conductivity of soil allows estimate the geothermal energy production efficiency [3]. Insight into the diffusivity of soot particles is key to designing metallic or ceramic porous media useful for the treatment of exhaust gases [4]. It is also relevant to know the hydrogen effective diffusivity when synthesizing materials for use in energy storage devices . A main subject of study in materials science is the relationship between transport properties and the microstructure descriptors of a porous medium (porosity (𝜙), pore diameter, connectivity, and tortuosity ((), among others). Through numerical, experimental, and analytical approaches, such as percolation theory, critical path analysis, and effective medium approximation; different mathematical expressions for these relations have been obtained (see references [6-8] for some examples). Porosity and tortuosity are among the most widely used descriptors. Conceptually, ( can be interpreted as a measurement of the geometrical complexity of a porous structure. Despite the simplicity of this concept, some authors state that tortuosity is neither understood nor consistently defined [9-11]. Carman [12] used it for the first time. He recoursed to this parameter to match his experimental permeability data with the results predicted by the analytical solution to the capillary tubes problem. In literature, four classes of ( can be found. Three of them are related to the process under study (hydraulic, electric, or diffusive tortuosity) and the fourth class is associated with the morphological properties of porous media (geometrical tortuosity). In 2013, Ghanbarian et al. [9] published a critical review of the deduction, meaning, and application field of each of these definitions. The geometrical tortuosity (𝜏𝑔)defines the relationship between the average length of the paths existing within the porous channels (〈Lg 〉) and the linear dimension of the system (L):

Some authors propose the use of the shortest pathway, or optimal path (Lopt), instead of 〈Lg 〉 (9(. However, this definition is not the most widely adopted. Experimentally, τg can be directly obtained from tomographic or micrographic image registration [13-15]; or indirectly by using some function that relates it to another microstructural descriptor. The last one is the most widely used option, and porosity is the preferable descriptor. In literature, many equations that relate τ and 𝜙 are found. Some of them are compiled by Sobieski and Lipinski [16]. The abundance of mathematical expressions is attributed to the fact that each study was carried out for particular cases of porous media, thus their application ranges are limited. Some of them are useful for 3D beds of spherical particles but cannot be used for cubic particles. Other studies apply to 2D sections of overlapping square particles but not to the non-overlapping ones. Many of these equations have a common deficiency; they predict a value of τg=0 when 𝜙 →1. Nevertheless, the geometrical tortuosity in this limit must be 1.

The lack of a predictive unified model has motivated investigations on more universal formulations. Based on the percolation theory and the finite-size scaling approach, Ghanbarian et al. [17]proposed a model that allows predict τg as a function of 𝜙 in a relatively wide range of two- and three-dimensional porous media. In this analytical model, which is described in more detail in section 2, the authors use the concepts of fractal length and dimension as well as the percolation theory and its critical exponents. They found that tortuosity follows a scaling law that gradually increases, while porosity decreases and diverges as 𝜙 approaches the percolation threshold (or critical porosity). The authors found that many of the equations relating τ and 𝜙 matched their proposal. Consequently, they could be considered within the percolation theory framework. This model has been used for predicting different transport parameters, such as effective permeability and electrical conductivity in concrete [6]; and thermal conductivity in porous alumina ceramics [6] and in materials with polyhedral porosity [7].

Using numerical methods, this work found that media with saturated, disordered, and randomly distributed pores follow a scaling law that relates these two microstructure descriptors: geometrical tortuosity and porosity. The percolation model is used for obtaining the digital representations of media with different porosities. An advanced algorithm was developed to derive tortuosity from the percolating paths that exist within the porous channels. Therefore, in addition to the average length, it is possible to obtain statistics, such as the median as well as the shortest and largest paths, which are not possible to derive using state-of-the art algorithms. The program validity is verified by comparing its results with recently published results. Moreover, the universality class to which the shortest path belongs was found. The paper is organized as follows: the analytical model presenting tortuosity as a scaling law is described in section 2. The procedure used for generating the digital images of porous media and the algorithm developed for obtaining the percolating paths are described in section 3. Section 4 presents the scaling law obtained numerically and other results, whereas section 5 reports the conclusion.

2. Percolation-based tortuosity model [17]

According to the percolation theory, the mean distance between two points in a same finite cluster of pores, or correlation length (χ), is given by the scaling law:

𝜙 c is the critical porosity, a minimum porosity value needed for the existence of a cluster of pores that spans the network from one side to another (percolating cluster). 𝜙 c depends on the type of percolation studied (bound or site, with or without trapping) and the lattice (square, cubic) [18]. For example, for site percolation in a square lattice, 𝜙 c has an accepted value of 0.5927. Similar to other percolation properties, this is an accepted value since an analytical result does not exist. The accepted values are numerically obtained. C is a fitting coefficient with an accepted magnitude of 0.85±0.4 for a square lattice [19]. ν is a critical exponent which only depends on the Euclidean dimension. Unlike 𝜙 c, ν neither depends on the type of percolation nor the type of lattice and has an exact value. For this reason, it is known as a universal exponent. Its exact value in 2D is ν=4/3.

The porous channels have a characteristic fractal configuration [17]. The total length of a fractal path (LF), constructed of steps of length ε, can be obtained from eq. (3):

Where Df is the fractal dimension

Close to the critical porosity, ε is inversely proportional to the correlation length and eq. (3) can be written as [17]:

Taking LF (χ) as the average length of the paths, the geometrical tortuosity can be expressed as [17]:

The correlation length predicted by percolation theory is only independent of the system size if χ≤L. When χ>L, eq. (2) and consequently eq. (5) must be modified to include this dependency. Ghanbarian et al. [17] proposed to use the scaling factor (C/L)1/ν.

The fact that τg=1 when 𝜙 →1 also suggests the inclusion of the factor |1- 𝜙 c+(C/L)^1/ν | in eq. (6) for ultimately obtaining the percolation-based tortuosity model

3. Numerical approximation to the tortuosity-porosity relationship

3.1 Porous media generation

The digital representations of the microstructures are obtained by randomly selecting a fraction 𝜙 of the nodes constituting a square lattice of linear dimension L (in this work an L=120 is arbitrarily selected). These nodes are turned into the mass center of a square pore with an edge equal to a lattice unit. The phase function (z(x,y)) for this system is:

Where 𝜙 is the porosity. A cluster of pores is defined as two or more mass centers that have z(x,y)=1 and occupy consecutive sites in the lattice. The control over the distribution and number of clusters is carried out by using the labeling algorithm proposed by Ramirez et al. [20]. As in real solids, there are pores in the microstructures that do not belong to the percolating cluster (occluded pores) and therefore cannot be taken into account when calculating τ.There are also blind pores that unjustifiably increase the percolating path length. Both types of pores are eliminated by doing the following transformation z(x,y)=1→0 (Fig. 1 a and b).

3.2 Algorithm for calculating tortuosity (Camper)

The algorithm, named Camper from caminos percolantes, has the objective to find as many paths as possible, but not all. The general diagram is presented in Scheme 1. Paths are found to depart from the left side of the digital representation to arrive to the right side and also in the opposite direction.

All the nodes in the lattice edges [i.e., the nodes (0,0) (0,L) and (L,0)⋯(L,L)] with z(x,y)=1 are the starting point for the search. These paths are formed by steps with an equivalent length of a lattice unit. In this case, there are four possibilities. When departing from the left side, we will have steps from (x,y) to (x,y±1), which will be cited as right (-) or left (+), or as (x±1,y), cited as forward (+) or backward (-). When departing from the right side, the signs change. The steps are preferentially taken forward provided that there are not obstacles. Obstacles are the pores walls. There are two types of pores walls in the digital representations: those that can only be evaded in one direction (to the left or right, they will be named single obstacles) and those that can be evaded in both directions (double obstacles) (Fig. 1 d). In the first case, the obstacle is bordered in the direction, in which z=1, until there is the possibility of continuing the path forward. In the second case, before bordering the pore wall a direction is randomly selected. Finally, 〈Lg 〉 is obtained after averaging the length of the different paths. Any configuration of steps with loops is rejected for the τg calculation.

Source: elaborated by the author

Figure 1 a) Digital representation of porous media with 𝜙 =0.65; the occluded pores are shown in gray. b) Lattice used for calculating 𝜏, after eliminating occluded and blind pores. c) Zoomed view of the framed region, whereby arrows indicate some blind pores. d) Arrows point out a single obstacle (only evadable to the left) and a double one (Obstacles are described in section 3.2.) 

Source: elaborated by the author

Scheme 1 Flowchart of the algorithm used to find the percolating paths. 

4. Results and discussion

4.1 Algorithm validation

The algorithm validity is corroborated by comparing results with those reported by Roque and Costa in 2020 [14] and Al-raoush and Madhoun in 2017 [15]. They also developed algorithms for calculating τg from image registration. From reference [14], the synthetic image called curved line is taken (see Fig. 2 a). From reference [15], the micrography of the natural media called M1, a mixed sand of silica and quartz, is taken (see Fig. 2 b). As stated in section 3.1, Camper only reads square images; for this reason, the rectangular image of M1 is cut as presented in Fig. 2 c. Both digital representations were adjusted to have a resolution of 240×240 pixels (L=240). Table 1 presents the results.

Source: adapted from references [14] and [15] with the author’s permission

Figure 2 Digital representations used for corroborating the algorithm. a) Curved line, b) Micrography of the mixed sand (M1 sample), the frame encloses the area used for the calculation. c) Zoomed view of the used area. As an example, in c) two of the paths found by the program are presented. Note how they evade the blind pores. 

Table 1 Results used for validating the algorithm 

a Length of the shortest percolating path or optimal path

b Median for paths

c Largest path length

d Tortuosities published in [14] and [15] and the one obtained by Camper

Source: elaborated by the author

One of the advantages of this program is the descriptive statistical analysis of the detected paths, which is not possible with the other proposals. Camper allows establish the number of paths used for calculating τg as well as the maximum and minimum lengths, and the median. In general, there is a good agreement between the tortuosities reported in literature and those calculated in this work. Particularly for M1, where the same value of tortuosity was obtained.

4.2 Detection of a scaling law

With the algorithm, τg is calculated in media with porosity from 0.63 to 0.95. Fig. 3 presents the results. This range of porosity was used because χ≤L; hence, the correlation length is independent of the size of the system and the scaling factor can be eliminated in eq. (7). However, while we are closer to 𝜙=0.63, the finite-size effects appear as greater standard deviations, which lead to greater uncertainty in the calculation of τg. Ghanbarian et al. [17] proposed to use the fractal dimension of the optimal path [21, 22] (Dopt=1.22 in 2D) as Df in eq. (7). By replacing the accepted values for 𝜙c and Dopt, the universal exponent in (7), and by solving the denominator such equation can be rewritten as:

Source: elaborated by the author

Figure 3 Relationship between the geometrical tortuosity and porosity. The circles indicate the numerically obtained results. The vertical bar represents the standard deviation. The solid line is the equation that better fits the data. The dashed line is the scaling law analytically obtained by Ghanbarian et al. (eq. 8). 

This equation is presented as a dashed line in Fig. 3. The solid line is the scaling law that better fits the simulated tortuosity values and is

There is a good correlation between both equations. After solving Df from ν-νDf=-0.35, a fractal dimension value of 1.26 was found. As expected, Df>Dopt since all percolating paths, and not only the shortest one, were used for its calculation. Consequently, Df can be interpreted as an average fractal dimension. When Dopt is used in eq. (7), as Ghanbarian et al. do so, the geometrical tortuosity is obtained in terms of the Lopt/L relationship and not as the relationship between the average length of the paths and the linear dimension of the system. To adjust to the most accepted definition of tortuosity (eq. 1), the use of the average fractal dimension found in this work, instead of Dopt, is proposed. Eq. (9) fulfills τg→1 when 𝜙 →1.

4.3 Optimal path and fractal dimension

In this section, the capacity of Camper for assigning a value to the optimal path length was used. As mentioned before, when χ>L the percolation properties depend on the system dimension. It has been possible to establish that under these conditions and in the percolation threshold, the length of such path scales with its Euclidean distance, meaning that Lopt~Le Dopt. [22,23] Le is the distance between the starting point and the arrival point, Le≥L. It is proposed to obtain Dopt from the slop of the relation log10 (Lopt ) vs.log10 (Le) constructed from twelve media with𝜙=0.6 and L=10,20,…,110,120. The critical porosity(𝜙c=0.5927) could not be used because, in some cases, the program does not find a percolating path. Results are presented in Fig. 4. Despite not using 𝜙c ,the value obtained for Dopt was 1.21, which agrees adequately with the accepted value of 1.22 [21,22]. This suggests that the optimal path existing within the porous channels can be considered as in the same universality class of the geometries described by loopless compressible invasion percolation [21,22,24].

Source: elaborated by the author

Figure 4 Logarithmic relationship between the optimal path length and the linear dimension of the system. The slope gives the value of the fractal dimension, Dopt 

5. Conclusions

The sequence followed by the developed algorithm, i.e., removing occluded and blind pores, paths starting in different points of each media edge, steps forward, which border obstacles, and avoiding loops allows to derive a solid statistical description of the percolating paths existing within porous media. For media with saturated, disordered, and randomly distributed pores, there is a power law that relates geometrical tortuosity and porosity. The exponent of such law can be obtained from the critical exponent of the correlation length and the average fractal dimension instead of the fractal dimension of the optimal path. The fractal dimension of the optimal path that exists within this type of microstructures is in the same universality class that the Eden growth observed in tissue or bacterial cultures or the geometry of the undirected polymers.

Acknowledgments

The author acknowledges Universidad de Antioquia for the time given for the development of this work. The author would like to thank Enago (www.enago.com) for the English language review.

References

[1] Nguyen, S.-T., Pham, D.-C., Vu, M.-N. and To, Q.-D., On the effective transport properties of heterogeneous materials. International Journal of Engineering Science, 104, pp. 75-86, 2016. DOI: 10.1016/j.ijengsci.2016.04.001 [ Links ]

[2] Di Palma, P.R., et al., Impact of synthetic porous medium geometric properties on solute transport using direct 3D Pore-Scale Simulations. Geofluids, (2019), pp. 1-13, 2019. DOI: 10.1155/2019/6810467 [ Links ]

[3] Jia, G.S., et al., Review of effective thermal conductivity models of rock-soil for geothermal energy applications. Geothermics, 77, pp. 1-11, 2019. DOI: 10.1016/j.geothermics.2018.08.001 [ Links ]

[4] Ou, Q., Maricq, M.M. and Pui, D.Y.H., Evaluation of metallic filter media for sub-micrometer soot particle removal at elevated temperature. Aerosol Science and Technology, 51(11), pp. 1303-1312, 2017. DOI: 10.1080/02786826.2017.1349871 [ Links ]

[5] Koizumi, K., Nobusada, K. and Boero, M., Hydrogen storage mechanism and diffusion in metal-organic frameworks. Physical Chemistry Chemical Physics, 21(15), pp. 7756-7764, 2019. DOI: 10.1039/c8cp07467d [ Links ]

[6] Xu, W. and Jiao, Y., Theoretical framework for percolation threshold, tortuosity and transport properties of porous materials containing 3D non-spherical pores. International Journal of Engineering Science, 134, pp. 31-46, 2019. DOI: 10.1016/j.ijengsci.2018.10.004 [ Links ]

[7] Ghanbarian, B. and Daigle, H., Thermal conductivity in porous media: Percolation-based effective-medium approximation. Water Resources Research, 52(1), pp. 295-314, 2016. DOI: 10.1002/2015wr017236 [ Links ]

[8] Berg, C.F. and Held, R., Fundamental transport property relations in porous media incorporating detailed pore structure description. Transport in Porous Media, 112(2), pp. 467-487, 2016. DOI: 10.1007/s11242-016-0661-7 [ Links ]

[9] Ghanbarian, B., Hunt, A.G., Ewing, R.P. and Sahimi, M., Tortuosity in Porous media: a critical review. Soil Science Society of America Journal, 77(5), pp. 1461-1477, 2013. DOI: 10.2136/sssaj2012.0435 [ Links ]

[10] Hunt, A.G. and Sahimi, M., Flow, Transport, and reaction in porous media: percolation scaling, critical-path analysis, and effective medium approximation. Reviews of Geophysics, 55(4), pp. 993-1078, 2017. DOI: 10.1002/2017rg000558 [ Links ]

[11] Niya, S.M.R. and Selvadurai, A.P.S., Statistical Correlation between permeability, porosity, tortuosity and conductance. Transport in Porous Media, 121(3), pp. 741-752, 2018. DOI: 10.1007/s11242-017-0983-0 [ Links ]

[12] Carman, P.C., Fluid flow through granular beds. Chemical Engineering Research and Design, 75, pp. S32-S48, 1997. DOI: 10.1016/S0263-8762(97)80003-2 [ Links ]

[13] Amien, M.N., Pantouw, G.T., Juliust, H. and Eljabbar Latief, F.D., Geometric tortuosity analysis of porous medium using simple neurite tracer. IOP Conference Series: Earth and Environmental Science, 311, art. 012041, 2019. DOI: 10.1088/1755-1315/311/1/012041 [ Links ]

[14] Roque, W.L. and Costa, R.R.A., A plugin for computing the pore/grain network tortuosity of a porous medium from 2D/3D MicroCT image. Applied Computing and Geosciences, 5, art. 100019, 2020. DOI: 10.1016/j.acags.2020.100019 [ Links ]

[15] Al-Raoush, R.I. and Madhoun, I.T., TORT3D: a MATLAB code to compute geometric tortuosity from 3D images of unconsolidated porous media. Powder Technology, 320, pp. 99-107, 2017. DOI: 10.1016/j.powtec.2017.06.066 [ Links ]

[16] Sobieski, W. and Lipiński, S., The analysis of the relations between porosity and tortuosity in granular beds. Technical Sciences, 1(20), pp. 75-85, 2017. DOI: 10.31648/ts.2912 [ Links ]

[17] Ghanbarian, B., et al., Percolation theory generates a physically based description of tortuosity in saturated and unsaturated porous media. Soil Science Society of America Journal, 77(6), pp. 1920-1929, 2013. DOI: 10.2136/sssaj2013.01.0089 [ Links ]

[18] Isichenko, M.B., Percolation, statistical topography, and transport in random media. Reviews of Modern Physics, 64(4), pp. 961-1043, 1992. DOI: 10.1103/RevModPhys.64.961 [ Links ]

[19] Kapitulnik, A., Aharony, A., Deutscher, G. and Stauffer, D., Self similarity and correlations in percolation. Journal of Physics A: Mathematical and General, 16(8), pp. L269-L274, 1983. DOI: 10.1088/0305-4470/16/8/003 [ Links ]

[20] Ramírez, A., Sierra, L., Mesa, M. and Restrepo, J., Simulation of nitrogen adsorption-desorption isotherms. Hysteresis as an effect of pore connectivity. Chemical Engineering Science, 60(17), pp. 4702-4708, 2005. DOI: 10.1016/j.ces.2005.03.004 [ Links ]

[21] Cieplak, M., Maritan, A. and Banavar, J.R., Invasion percolation and eden growth: geometry and universality. Physical Review Letters, 76(20), pp. 3754-3757, 1996. DOI: 10.1103/PhysRevLett.76.3754 [ Links ]

[22] Porto, M., Havlin, S., Schwarzer, S. and Bunde, A., Optimal path in strong disorder and shortest path in invasion percolation with trapping. Physical Review Letters, 79(21), pp. 4060-4062, 1997. DOI: 10.1103/PhysRevLett.79.4060 [ Links ]

[23] Sahimi, M., Applications of Percolation Theory, Taylor & Francis, USA, 1994, 17 P. [ Links ]

[24] Cieplak, M., Maritan, A. and Banavar, J.R., Optimal paths and domain walls in the strong disorder limit. Physical Review Letters, 72(15), pp. 2320-2323, 1994. DOI: 10.1103/PhysRevLett.72.2320 [ Links ]

A. Ramirez-Velez, received his BSc. in Chemistry in 1998, and his PhD. in Chemistry Sciences in 2004, all of them from the Universidad de Antioquia, Medellin, Colombia. From 2005 to 2008, he was assistant professor in the School of Chemistry, Universidad Nacional de Colombia, Medellin. Currently, he is associate professor in the Institute of Chemistry, Universidad de Antioquia, Medellin. His research interests include: the use of computational methods for characterizing the structure of porous media and studying their transport properties, and the synthesis and application of silicon-based materials.ORCID: 0000-0002-9240-3477

How to cite: Ramírez-Vélez, A., Numerical approximation to the scaling law describing the geometrical tortuosity in porous media.. DYNA, 88(217), pp. 91-96, April - June, 2021.

Received: June 30, 2020; Revised: December 12, 2020; Accepted: February 26, 2021

Creative Commons License The author; licensee Universidad Nacional de Colombia