SciELO - Scientific Electronic Library Online

 
vol.64 número3Wavelet characterization of hyper-chaotic time seriesRecovery of transit times and frequencies of multiple pulses via the short-time Fourier transform índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.64 no.3 México may./jun. 2018

 

Research

Research in Gravitation, Mathematical Physics and Field Theory

Fractal imbibition in Koch’s curve-like capillary tubes

D. Samayoaa 

L. Alvarez-Romeroa 

L.A. Ochoa-Ontiverosa 

L. Damián-Adamea 

E. Victoria-Tobona 

G. Romero-Paredesb 

a Instituto Politécnico Nacional, SEPI-ESIME, Unidad Poblacional Adolfo López Mateos, CDMEX, 07738, México.

b Department of Electrical Engineering, Solid State Electronic Section (SEES), Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, CDMEX, México.


ABSTRACT

Fractal dimension effects in capillary imbibition process are analytically studied. The fractal formulation of tortuous flow with the assumption of a fractal tortuous path introduced by Cai and Yu is used to analyse the capillary rise through the tubes with deterministic fractal geometry. Capillary rise in Koch’s curve-like tubes was investigated. A new permeability parameter that takes into account the tortuosity of the flow path is deduced, and a geometrical relationship for fractal dimension of flow tortuosity (dτ) in porous media is obtained. The equilibrium height time as a function of fractal dimension of the flow tortuosity in capillary tubes with tortuous path was derived.

Keywords: Spontaneous imbibition; linear fractal; Koch curve; fractal dimension of flow tortuosity

PACS: 47.53-+n; 47.56.+r

1. Introduction

Spontaneous capillary imbibition is a transport phenomenon present in a variety of technological applications such as oil recovery, building materials, soil science, textile and hydrology1. Due to this diversity of applications, an infinity of theoretical2-6 and experimental7-14 studies have been carried out based on the pioneering works of Lucas15 and Washburn16 in order to understand the imbibition mechanism and the related phenomena (for review, see1,17,18, and references therein).

It was experimentally demonstrated that the imbibition speed becomes slower than that of Lucas-Washburn Regime (xt) 7-14 as a consequence of the heterogeneity of flow in porous media due to the complex network of randomly distributed pores that connect to each other forming tortuous capillaries through which fluid flows. This tortuous path has a fractal character and the Lucas-Washburn equation is scarce to model imbibition process in permeable materials that possess fractal architecture.

Li and Zhao19 added a fractal parameter to the classical model of spontaneous imbibition to describe the heterogeneity of the porous medium as xtdf-2 where 2<df<3 is the fractal dimension of the medium. However, Cai and Yu5 presented a fractal formulation of tortuous flow through the assumption of a fractal tortuous path given by xt1/dτ, where dτ is the fractal dimension of flow tortuosity to depict the previous experimental results for a scaling exponent within the range of 0.25δ=1/dτ<0.5. The mentioned fractal approach is used in this work to study of liquid travel distance in the Euclidean case regarding the fractal case and find the time when the equilibrium height is achieved as a function of the fractal dimension of flow tortuosity in tubes with the shape of the Koch curve, Modified Koch curve, Minkowski curve and the “Carpintieri curve"20.

The objective of this work is to investigate the effects of the fractality in spontaneous imbibition processes on capillaries with paths similar to theoretical curves with exact fractal dimension.

The rest of the Letter is organized as follows. In Sec. 2 the mathematical tools needed in this paper were studied and defined. In Sec. 3 a new permeability parameter and an analytical model to describe the capillary rise by spontaneous imbibition on tubes with deterministic fractal geometry are established. In Sec. 4 an illustrative example in order to discuss some physical implications has been solved. In Sec. 5 the main findings and conclusions are outlined.

2. Basic tools in the fractal imbibition

In what follows, some basic concepts of fractal imbibition in highly tortuous capillaries are presented.

2.1. Governing equations in the Euclidean space

The governing equation of the liquid rise in a vertical capillary tube embedded in the Euclidean space E 1, when the capillary forces are dominant and the radius is very small, is given by a balance in the capillary, viscous and gravitational forces21,22, respectively represented as:

4σcosθΦ=32μΦ2xdxdt+ρgx, (1)

where σ is the surface tension of the liquid having viscosity μ, θ is the contact angle between the liquid and the capillary surface, Φ is the capillary diameter, x denotes the distance penetrated by the liquid, t is the imbibition time, ρ is the fluid density and, g is the gravitational acceleration. The process of the liquid rise in a vertically straight capillary tube is obtained from Eq. (1) as dx/dt=K(Pc-ρgx)/μx, where K=Φ2/32 is the intrinsic permeability of the tube1,23. When an equilibrium is reached, the force that drives the capillary rise Pc=4σcosθ/Φ equals the weight of the column of liquid, ρgx; therefore the net force on the liquid vanishes and the rising of the liquid stops.

Liquid-air interface equilibrium height is defined as xeq=Pc/ρg; and the liquid rise in the capillary tube between the time of initial contact and the final equilibrium is described by the relation1,23:

dxdt=Kμρgxeqx-1. (2)

When the gravitational forces can be considered very small, for example when the capillary tube is in the upright position and the penetration of the liquid x is very small, or the gravitational force is equal to zero (when the capillary tube is in the horizontal position), the above mentioned equation reduced to dx/dt=Φσcosθ/8μx and integrating it with the initial conditions x=0,t=0, the following Lucas-Washburn relation holds:

x(t)=2KPcμ12t12, (3)

when the gravitational forces are included, integrating Eq. (2) with the initial conditions x = 0, t = 0, we get:

x(t)xeq1-e-t/teq, (4)

where teq=-128σ cosθμ/Φ3ρ2g2.

2.2. Governing equations in fractal space

The actual tortuous path of a flow in a porous medium is defined as xτ=τx, where τ is the tortuosity. In the24 a scaling relationship for a flow through heterogeneous media was developed for the actual length xτ versus the scale of observation ϵ given by xτ=ϵ1-dτxdτ where dτ is the fractal dimension of the flow tortuosity. Expression mentioned implies the property of self-similarity, which means that the value of dτ is constant for a range of length scales εcxl, where εc>0 is the lower cutoff with value of the order of an average pore radius and l is the upper cutoff that describes the straight-line distance that the particle travels between the starting and the ending points of fractal (tortuous) path (see Fig. 1). Hence the diameters of capillaries are analogous to the length scale ϵ 25,26 and the fractal scaling relationship between the diameter Φ and the length of capillaries x and xτ in a porous medium can be written as:

xτ=xxΦdτ-1, (5)

where

1dτ=1+lnτln(x/Φ)<2. (6)

Figure 1 Fractal parameters of scale invariance for dτ . Note that lτ as c0

When dτ=1 the capillary tube is straight and theone-dimensional Euclidean case is presented, meanwhile for dτ=2 the capillary tube is a highly tortuous line so irregular that fills a two-dimensional space. The capillaries with tortuous path are embedded in the Euclidean space E2.

Differentiating Eq. (5) with respect to time t for a single capillary results in24:

vτ=dτΦ1-dτxdτ-1v, (7)

where vτ=dxτ/dt is the actual velocity of liquid through distance of a tortuous capillary; v=dx/dt is the straight-line imbibition velocity; and as for the Euclidian case dτ=1, it holds that vτ=v. The scaling ratio between both velocities can be rewritten as

dxτdt=dτΦ1-dτxdτ-1dxdt. (8)

When a wetting liquid is contacted with a tortuous capillary of any shape (Euclidean or Fractal), the capillary rise is described as5,16:

dxτdt=Kμxxτρgxeqx-1. (9)

Substituting Eqs. (5) and (8) in (9) the following differential equation that governs the capillary rise in a tortuous capillary is obtained :

dxdt=Kμρgdτx2dτ-2Φ2-2dτxeqx-1. (10)

If the gravitational force is negligible, Eq. (10) is reduced to:

dxdt=KμdτΦ2-2dτPcx2dτ-1, (11)

with the initial conditions x = 0, t = 0. Integration leads to Lucas-Washburn-Cai Equation5:

x(t)=2KPcΦ2-2dτμ12dτtδ=12dτ, (12)

and when dτ=1 the last expression is reduced to the classical of Lucas-Washburn equation.

3. Capillary rise on the linear fractals

Fractal propierties of capillary tubes with deterministic fractal geometries are obtained (see Fig. 2) for the first six iterations of each linear fractal including their fractal dimension, df, as it is shown in Table I.

Figure 2 Capillary tubes with deterministic fractal geometries for first 6 iterations. a) Modified Koch curve27; b) Classic Koch curve; c) Carpintieri curve; d) Cuadratic Koch curve or “Minkowski sausage”. 

Table I Properties of capillary tubes with deterministic fractal geometries, unitary porosity, constant diameter Φ = 0.001 cm and vertical height x = 64 cm. 

Parameter Fractal curve
Koch MKa Carpintieri Minkowski
df 001.261 001.161 001.340 0001.500
dτ0 001.000 001.000 001.000 0001.000
dτ1 001.025 001.020 001.036 0001.062
dτ2 001.051 001.040 001.073 0001.125
dτ3 001.077 001.060 001.109 0001.187
dτ4 001.103 001.090 001.146 0001.250
dτ5 001.129 001.100 001.183 0001.313
xτ0 064.000 064.000 064.000 0064.000
xτ1 085.330 080.000 096.000 0128.000
xτ2 113.770 100.000 144.000 0256.000
xτ3 151.703 125.000 216.000 0512.000
xτ4 202.271 156.250 324.000 1024.000
xτ5 269.695 195.310 486.000 2048.000

a Modified Koch curve.

The generalized permeability for a tortuous capillary tube is given by28-31:

K=Φ232ϕτ. (13)

However, the fractal approach to capillary cylinder permeability with fractal geometry is obtained inserting Eq. (5) in Eq. (13):

Kτ=Φdτ+132xdτ-1, (14)

and for capillary tubes with unitary porosity, ϕ=1. Where Kτ, is a permeability ratio for permeable media with scale invariance in materials where a pre-fractal pore network exists32.

It was shown that the equilibrium height xeq is the same for both, straight tubes and tubes with tortuous path. Substituting xeq, τ and Kτ in Eq. (10) and the following expression that describes the capillary rise by spontaneous imbibition in tubes with fractal geometry is obtained:

dxdt=Kτρgμdτxdτ-1Φ1-dτxeqx-1. (15)

When dτ=1, Eq. (15) reduces to Eq. (2). In the first stage of imbibition the liquid rise in the capillary tube is given by x(t)tδ where the time exponent δ=1/2dτ and the second stage of imbibition can be found solving numerically Eq. (15). The fractal dimension of flow tortuosity increases with the increase of the capillary tubes tortuosity, meanwhile in the case of tubes with a shape similar to linear fractals shown in Table I, the tortuosity increases with the increase of the fractal curve iteration (see Fig. 3) as described in Eq. (6), or it can be calculated by the following relation:

dτ=λilnτ+1, (16)

where λi=0.2327709525Φ0.1667743386 (see the insert of Fig. 3 where the slope λi depends on the capillary diameter by the power law).

Figure 3 Semi-log plot τ vs. dτ for different capillary diameters of linear fractals with logarithmic fittings, where dτ = slope ln(τ) + 1 represented by straight lines. Dashed line shows a graph obtained from linear fractals (that is in agreement with Fig. 1 of Ref. 5. Insert shows data of Fig. 3 versus capillary diameter, fitted curve for linear regression analysis. 

4. Example in Koch’s curve-like capillary tubes

In this section an illustrative example to clarify the physical implications of the introduced models is presented.

Consider a tortuous flow path, similar to a classic Koch curve, shown in Figure 1.b, with the following mechanical properties: θ=0, ρ=0.998 g/cm3, μ=0.01 dina sec/cm2, σ=725 dina/cm and x=64 cm for the first (i=0,1,,5) iterations. The results of the first stage of spontaneous imbibition, are in agreement with the time scaling exponent δ=0.5dτ-1 as it is shown in Fig. 4. The influence of the fractality is reflected only by the fractal dimension of flow tortuosity, meanwhile the fractal dimension of the deterministic fractal-like capillary tube does not influence on the behaviour of the flow.

Figure 4 Comparison between the Euclidean case and fractal case in the first stage of imbibition. Graph shows Washburn Regime in Euclidean case (dτ = 1) while the fractal case (dτ > 1) evolves slower then the Washburn regime. 

The flow behaviour in the second stage of imbibition, is obtained solving Eq. (15), where the equilibrium height is always the same and does not depend on the fractal dimension of the flow tortuosity (see Fig. 5). However, the time that fluid requires to reach the equilibrium height depends directly on dτ, as:

t(xeq)dτα, (17)

where α=18.44, as in Fig. 6.

Figure 5 Comparison between Euclidean case and fractal case in the second stage of imbibition. Travel distance as function of time, for different iterations of the basic Koch’s curve-like capillary tube. 

Figure 6 Time equilibrium height versus fractal dimension of flow tortuosity for 1 ≤ dτ < 2. 

The previous results are common for ideal liquids, when the contact angle is equal to zero. The capillary rise is directly affected by the hydrostatic effect depending on the contact angle. This effect is calculated with the information of Fig. 4 according to Lucas-Washburn Eq. (3) and such effects are shown in Fig. 7. In an upcoming report the experimental results of fractal imbibition in tortuous Koch’s curves-like capillary tubes will be given.

Figure 7 Contact angle as time function deduced from the Figure 4 data according to Eq.(3) and the parameters from Table I

5. Conclusions

In this paper the permeability relation for flow paths in capillary tubes with fractal geometry was generalized. This relation was used in Eq. (10) (Jian Chao Cai Equation) to describe the spontaneous fractal imbibition model given by Eq. (15). Also it was found that the fractal dimension of flow tortuosity increases, as the tortuosity and the capillary diameters of the cylinders increases, and does not depend on the fractal dimension of capillary tube (see Eq. (16)). An illustrative example of capillaries with the shape of linear fractal of Koch curve was presented. Results obtained are in agreement with the previously reported findings standard calculations and with the Lucas-Washburn-Cai equation, showing that as the fractality increases the penetration distance decreases. Finally, it was found that the necessary time to reach the equilibrium height of the tortuous capillary tube in the second stage of spontaneous imbibition is a function of the fractal dimension of flow tortuosity such that t(xeq)dτα.

This results provide a more detailed description of the physical phenomena of the spontaneous imbibition in linear fractal-like capillary tubes.

Acknowledgments

This work was supported by the Instituto Politécnico Nacional under the research SIP-IPN grant No. 20180008.

REFERENCES

1. M. Alava, M. Dube, and M. Rost, Advances in Physics 53 (2004) 83. [ Links ]

2. R. Laughlin and J. Davies, Textile Res. J., 31 (1961) 39. [ Links ]

3. A. Leventis, D. Verganelakis, M. Halse, and J. Webber, Transp. Porous Med. 39 (2000) 143. [ Links ]

4. C. Hammecker, L. Barbiro, P. Boivin, J. Maeght, and E. Diaw, Transp. Porous Med. 54 (2004) 193. [ Links ]

5. J. Cai and B. Yu, Chinese Physics Letters 27 (2010) 054701. [ Links ]

6. R. Masoodi, E. Languri, and A. Ostadhos- sein, journal of Colloid and Interface Science 389 (2013) 268. [ Links ]

7. J. D. Boer, Textile Research Journal 50 (1980) 624. [ Links ]

8. V. Horváth and H. Stanley, Physical Review E 52 (1995) 525166. [ Links ]

9. C. Lam and V. Horváth, Physical Review Letters 85 (2000) 1238. [ Links ]

10. A. Balankin and O. Susarrey, Physical Review Letters 96 (2006) 056101. [ Links ]

11. J. Milczareck, I. Fijal-Kirejczyc, and J. Zoladek, Acta Physica Polonica A 113 (2008) 1245. [ Links ]

12. T. Kwon and A. H. dan S.E. ODonnell, Physical Review E 54 (1996) 685. [ Links ]

13. A. Miranda, I. Menezes-Sobrinho, and M. S. Couto, Physical Review Letters 104 (2010) 086101. [ Links ]

14. A. Balankin, E. G. Otamendi, D. Samayoa, J. Patino, and M. Rodríguez, Physical Review E 85 (2012) 036313. [ Links ]

15. R. Lucas, Kolloid-Zeitschrift 23 (1918) 15. [ Links ]

16. E. Washburn, Physical Review 17 (1921) 273. [ Links ]

17. N. Morrow and G. Mason, Current Opinion in Colloid and Interface Science 6 (2001) 321. [ Links ]

18. J. Cai and B. Yu, Transport in Porous Media 89 (2011) 251. [ Links ]

19. K. Li and H. Zhao, Transport in Porous Media 91 (2012) 363. [ Links ]

20. A. Carpintieri, Mechanics of Materials 18 (1994) 89. [ Links ]

21. G. Martic, F. Gentner, D. Seveno, D. Coulon, and J. D. Coninck, Langmuir 18 (2002) 7971. [ Links ]

22. J. Polansky and T. Kaya, International Journal of Thermal Sciences 91 (2015) 25. [ Links ]

23. M. Dubé, M. Rost, and M. Alava, The European Physical Journal B-Condensed Matter and Complex Systems, 15 (2000) 691. [ Links ]

24. S. W. Wheatcraft and S. W. Tyler, Water Resources Research, 24 (1988) 566. [ Links ]

25. B. Yu and P. Cheng, International Journal of Heat and Mass Transfer, 45 (2002) 2983. [ Links ]

26. B. Yu, Chinesse Physics Letters 22 (2005) 158. [ Links ]

27. L. Damian, A. Kryvko, D. Samayoa, and A. R. Castellanos, Rev. Mex. Fis. 63 (2017) 12. [ Links ]

28. M. Panda and W. Lak, AAPG Bulletin 78 (1994) 1028. [ Links ]

29. A. Costa, Geophysical Research Letters 33 (2006) L02318. [ Links ]

30. N. Montersen, F. Okkels, and H. Bruus, Physical Review E 71 (2005) 057301. [ Links ]

31. F. Dullien, M. El-Sayed, and V. Batra, Journal of Colloid and Interface Science 60 (1977) 497. [ Links ]

32. A. Balankin, J. Valdivia, J. Márquez, O. Susarrey, and M. A. Solorio-Avila, Phyics Letters A 380 (2016) 2767. [ Links ]

Received: December 08, 2017; Accepted: February 16, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License