Tortuosity

Last updated December 10, 2024

Tortuosity is widely used as a critical parameter to predict transport properties of porous media, such as rocks and soils. But unlike other standard microstructural properties, the concept of tortuosity is vague with multiple definitions and various evaluation methods introduced in different contexts. Hydraulic, electrical, diffusional, and thermal tortuosities are defined to describe different transport processes in porous media, while geometrical tortuosity is introduced to characterize the morphological property of porous microstructures.^[1]

Tortuosity in 2-D

Subjective estimation (sometimes aided by optometric grading scales)^[2] is often used.

The simplest mathematical method to estimate tortuosity is the arc-chord ratio: the ratio of the length of the curve (C) to the distance between its ends (L):

\tau ={\frac {C}{L}}

Arc-chord ratio equals 1 for a straight line and is infinite for a circle.

Another method, proposed in 1999,^[3] is to estimate the tortuosity as the integral of the square (or module) of the curvature. Dividing the result by length of curve or chord has also been tried.

In 2002 several Italian scientists^[4] proposed one more method. At first, the curve is divided into several (N) parts with constant sign of curvature (using hysteresis to decrease sensitivity to noise). Then the arc-chord ratio for each part is found and the tortuosity is estimated by:

\tau ={\frac {N-1}{L}}\cdot \sum \limits _{i=1}^{N}{\left({{\frac {L_{i}}{S_{i}}}-1}\right)}

In this case tortuosity of both straight line and circle is estimated to be 0.

In 1993^[5] Swiss mathematician Martin Mächler proposed an analogy: it’s relatively easy to drive a bicycle or a car in a trajectory with a constant curvature (an arc of a circle), but it’s much harder to drive where curvature changes. This would imply that roughness (or tortuosity) could be measured by relative change of curvature. In this case the proposed "local" measure was derivative of logarithm of curvature:

{\frac {d}{dx}}\log \left(\kappa \right)={\frac {\kappa '}{\kappa }}

However, in this case tortuosity of a straight line is left undefined.

In 2005 it was proposed to measure tortuosity by an integral of square of derivative of curvature, divided by the length of a curve:^[6]

\tau ={\frac {\int \limits _{t_{1}}^{t_{2}}{\left({\kappa '\left(t\right)}\right)^{2}}dt}{L}}

In this case tortuosity of both straight line and circle is estimated to be 0.

Fractal dimension has been used to quantify tortuosity.^[7] The fractal dimension in 2D for a straight line is 1 (the minimal value), and ranges up to 2 for a plane-filling curve or Brownian motion.^[8]

In most of these methods digital filters and approximation by splines can be used to decrease sensitivity to noise.

Tortuosity in 3-D

Usually subjective estimation is used. However, several ways to adapt methods estimating tortuosity in 2-D have also been tried. The methods include arc-chord ratio, arc-chord ratio divided by number of inflection points and integral of square of curvature, divided by length of the curve (curvature is estimated assuming that small segments of curve are planar).^[11] Another method used for quantifying tortuosity in 3D has been applied in 3D reconstructions of solid oxide fuel cell cathodes where the Euclidean distance sums of the centroids of a pore were divided by the length of the pore.^[12]

Applications of tortuosity

Tortuosity of blood vessels (for example, retinal and cerebral blood vessels) is known to be used as a medical sign.

In mathematics, cubic splines minimize the functional, equivalent to integral of square of curvature (approximating the curvature as the second derivative).

In many engineering domains dealing with mass transfer in porous materials, such as hydrogeology or heterogeneous catalysis, the tortuosity refers to the ratio of the diffusivity in the free space to the diffusivity in the porous medium ^[13] (analogous to arc-chord ratio of path). Strictly speaking, however, the effective diffusivity is proportional to the reciprocal of the square of the geometrical tortuosity^[14]

Because of the porous materials found in several layers of fuel cells, the tortuosity is an important variable to be analyzed.^[15] It is important to notice that there are different kind of tortuosity, i.e., gas-phase, ionic and electronic tortuosity.

In acoustics and following initial works by Maurice Anthony Biot in 1956, the tortuosity is used to describe sound propagation in fluid-saturated porous media. In such media, when frequency of the sound wave is high enough, the effect of viscous drag force between the solid and the fluid can be ignored. In this case, velocity of sound propagation in the fluid in the pores is non-dispersive and compared with the value of the velocity of sound in the free fluid is reduced by a ratio equal to the square root of the tortuosity. This has been used for a number of applications including the study of materials for acoustic isolation, and for oil prospection using acoustics means.

In analytical chemistry applied to polymers and sometimes small molecules tortuosity is applied in gel permeation chromatography (GPC) also known as size exclusion chromatography (SEC). As with any chromatography it is used to separate mixtures. In the case of GPC the separation is based on molecular size and it works by the use of stationary media with an appropriate porous microstructure and adequate pore dimensions and distribution. The separation occurs because larger molecules cannot enter the smaller porosity for steric hindrance reasons (constrictivity of the narrow pores) and remain in the macropores, eluting more quickly, while smaller molecules can pass into smaller pores and take a longer, more tortuous path and elute later.

In pharmaceutical sciences, tortuosity is used in relation to diffusion-controlled release from solid dosage forms. Insoluble matrix formers, such as ethyl cellulose, certain vinyl polymers, starch acetate and others control the permeation of the drug from the preparation and into the surrounding liquid. The rate of mass transfer per area unit is, among other factors, related to the shape of polymeric chains within the dosage form. Higher tortuosity or curviness retards mass transfer as it acts obstructively on the drug particles within the formulation.

HVAC makes extensive use of tortuosity in evaporator and condenser coils for heat exchangers, whereas Ultra-high vacuum makes use of the inverse of tortuosity, which is conductivity, with short, straight, voluminous paths.

Tortuosity has been used in ecology to describe the movement paths of animals.^[8]

Related Research Articles

<span class="mw-page-title-main">Curvature</span> Mathematical measure of how much a curve or surface deviates from flatness

In mathematics, curvature is any of several strongly related concepts in geometry that intuitively measure the amount by which a curve deviates from being a straight line or by which a surface deviates from being a plane. If a curve or surface is contained in a larger space, curvature can be defined extrinsically relative to the ambient space. Curvature of Riemannian manifolds of dimension at least two can be defined intrinsically without reference to a larger space.

In trigonometry, the law of sines, sine law, sine formula, or sine rule is an equation relating the lengths of the sides of any triangle to the sines of its angles. According to the law, $where a, b, and c are the lengths of the sides of a triangle, and α, β, and γ are the opposite angles, while R is the radius of the triangle's circumcircle. When the last part of the equation is not used, the law is sometimes stated using the reciprocals; The law of sines can be used to compute the remaining sides of a triangle when two angles and a side are known—a technique known as triangulation. It can also be used when two sides and one of the non-enclosed angles are known. In some such cases, the triangle is not uniquely determined by this data and the technique gives two possible values for the enclosed angle.$

<span class="mw-page-title-main">Fermat's spiral</span> Spiral that surrounds equal area per turn

A Fermat's spiral or parabolic spiral is a plane curve with the property that the area between any two consecutive full turns around the spiral is invariant. As a result, the distance between turns grows in inverse proportion to their distance from the spiral center, contrasting with the Archimedean spiral and the logarithmic spiral. Fermat spirals are named after Pierre de Fermat.

Differential geometry of curves is the branch of geometry that deals with smooth curves in the plane and the Euclidean space by methods of differential and integral calculus.

<span class="mw-page-title-main">Frenet–Serret formulas</span> Formulas in differential geometry

In differential geometry, the Frenet–Serret formulas describe the kinematic properties of a particle moving along a differentiable curve in three-dimensional Euclidean space $or the geometric properties of the curve itself irrespective of any motion. More specifically, the formulas describe the derivatives of the so-called tangent, normal, and binormal unit vectors in terms of each other. The formulas are named after the two French mathematicians who independently discovered them: Jean Frédéric Frenet, in his thesis of 1847, and Joseph Alfred Serret, in 1851. Vector notation and linear algebra currently used to write these formulas were not yet available at the time of their discovery.$

In mathematics, the mean curvature $of a surface is an extrinsic measure of curvature that comes from differential geometry and that locally describes the curvature of an embedded surface in some ambient space such as Euclidean space.$

In fluid mechanics, materials science and Earth sciences, permeability is a measure of the ability of a porous material to allow fluids to pass through it.

<span class="mw-page-title-main">Porous medium</span> Material containing fluid-filled voids

In materials science, a porous medium or a porous material is a material containing pores (voids). The skeletal portion of the material is often called the "matrix" or "frame". The pores are typically filled with a fluid. The skeletal material is usually a solid, but structures like foams are often also usefully analyzed using concept of porous media.

Soil mechanics is a branch of soil physics and applied mechanics that describes the behavior of soils. It differs from fluid mechanics and solid mechanics in the sense that soils consist of a heterogeneous mixture of fluids and particles but soil may also contain organic solids and other matter. Along with rock mechanics, soil mechanics provides the theoretical basis for analysis in geotechnical engineering, a subdiscipline of civil engineering, and engineering geology, a subdiscipline of geology. Soil mechanics is used to analyze the deformations of and flow of fluids within natural and man-made structures that are supported on or made of soil, or structures that are buried in soils. Example applications are building and bridge foundations, retaining walls, dams, and buried pipeline systems. Principles of soil mechanics are also used in related disciplines such as geophysical engineering, coastal engineering, agricultural engineering, and hydrology.

In geometry, an intrinsic equation of a curve is an equation that defines the curve using a relation between the curve's intrinsic properties, that is, properties that do not depend on the location and possibly the orientation of the curve. Therefore an intrinsic equation defines the shape of the curve without specifying its position relative to an arbitrarily defined coordinate system.

The contact angle is the angle between a liquid surface and a solid surface where they meet. More specifically, it is the angle between the surface tangent on the liquid–vapor interface and the tangent on the solid–liquid interface at their intersection. It quantifies the wettability of a solid surface by a liquid via the Young equation.

Curved space often refers to a spatial geometry which is not "flat", where a flat space has zero curvature, as described by Euclidean geometry. Curved spaces can generally be described by Riemannian geometry, though some simple cases can be described in other ways. Curved spaces play an essential role in general relativity, where gravity is often visualized as curved spacetime. The Friedmann–Lemaître–Robertson–Walker metric is a curved metric which forms the current foundation for the description of the expansion of the universe and the shape of the universe. The fact that photons have no mass yet are distorted by gravity, means that the explanation would have to be something besides photonic mass. Hence, the belief that large bodies curve space and so light, traveling on the curved space will, appear as being subject to gravity. It is not, but it is subject to the curvature of space.

In the differential geometry of curves in three dimensions, the torsion of a curve measures how sharply it is twisting out of the osculating plane. Taken together, the curvature and the torsion of a space curve are analogous to the curvature of a plane curve. For example, they are coefficients in the system of differential equations for the Frenet frame given by the Frenet–Serret formulas.

A parametric surface is a surface in the Euclidean space $which is defined by a parametric equation with two parameters :\mathbb {R} ^{2}\to \mathbb {R} ^{3}} . Parametric representation is a very general way to specify a surface, as well as implicit representation. Surfaces that occur in two of the main theorems of vector calculus, Stokes' theorem and the divergence theorem, are frequently given in a parametric form. The curvature and arc length of curves on the surface, surface area, differential geometric invariants such as the first and second fundamental forms, Gaussian, mean, and principal curvatures can all be computed from a given parametrization.$

In hyperbolic geometry, a horocycle, sometimes called an oricycle or limit circle, is a curve of constant curvature where all the perpendicular geodesics ( normals) through a point on a horocycle are limiting parallel, and all converge asymptotically to a single ideal point called the centre of the horocycle. In some models of hyperbolic geometry it looks like the two "ends" of a horocycle get closer and closer to each other and closer to its centre, this is not true; the two "ends" of a horocycle get further and further away from each other and stay at an infinite distance off its centre. A horosphere is the 3-dimensional version of a horocycle.

Nuclear magnetic resonance (NMR) in porous materials covers the application of using NMR as a tool to study the structure of porous media and various processes occurring in them. This technique allows the determination of characteristics such as the porosity and pore size distribution, the permeability, the water saturation, the wettability, etc.

The table of chords, created by the Greek astronomer, geometer, and geographer Ptolemy in Egypt during the 2nd century AD, is a trigonometric table in Book I, chapter 11 of Ptolemy's Almagest, a treatise on mathematical astronomy. It is essentially equivalent to a table of values of the sine function. It was the earliest trigonometric table extensive enough for many practical purposes, including those of astronomy. Since the 8th and 9th centuries, the sine and other trigonometric functions have been used in Islamic mathematics and astronomy, reforming the production of sine tables. Khwarizmi and Habash al-Hasib later produced a set of trigonometric tables.

In mathematics, the curve-shortening flow is a process that modifies a smooth curve in the Euclidean plane by moving its points perpendicularly to the curve at a speed proportional to the curvature. The curve-shortening flow is an example of a geometric flow, and is the one-dimensional case of the mean curvature flow. Other names for the same process include the Euclidean shortening flow, geometric heat flow, and arc length evolution.

In Einstein's theory of general relativity, the interior Schwarzschild metric is an exact solution for the gravitational field in the interior of a non-rotating spherical body which consists of an incompressible fluid and has zero pressure at the surface. This is a static solution, meaning that it does not change over time. It was discovered by Karl Schwarzschild in 1916, who earlier had found the exterior Schwarzschild metric.

Pore structure is a common term employed to characterize the porosity, pore size, pore size distribution, and pore morphology of a porous medium. Pores are the openings in the surfaces impermeable porous matrix which gases, liquids, or even foreign microscopic particles can inhabit them. The pore structure and fluid flow in porous media are intimately related.

References

↑ Fu, Jinlong; Thomas, Hywel R.; Li, Chenfeng (January 2021). "Tortuosity of porous media: Image analysis and physical simulation" (PDF). Earth-Science Reviews. 212: 103439. doi:10.1016/j.earscirev.2020.103439. S2CID 229386129.
↑ Richard M. Pearson. Optometric Grading Scales for use in everyday practice. Optometry Today, Vol. 43, No. 20, 2003 Archived 2012-04-04 at the Wayback Machine ISSN 0268-5485
↑ Hart, William E.; Goldbaum, Michael; Cote, Brad; Kube, Paul; Nelson, Mark R. (1999). "Automated measurement of retinal vascular tortuosity". International Journal of Medical Informatics. 53 (2–3): 239–252. doi:10.1016/s1386-5056(98)00163-4. PMID 10193892. Archived from the original on 2009-01-09.
↑ Enrico Grisan, Marco Foracchia, Alfredo Ruggeri. A novel method for automatic evaluation of retinal vessel tortuosity. Proceedings of the 25th Annual International Conference of the IEEE EMBS, Cancun, Mexico, 2003
↑ M. Mächler, Very smooth nonparametric curve estimation by penalizing change of curvature, Technical Report 71, ETH Zurich, May 1993
↑ Patasius, M.; Marozas, V.; Lukosevicius, A.; Jegelevicius, D.. Evaluation of tortuosity of eye blood vessels using the integral of square of derivative of curvature // EMBEC'05: proceedings of the 3rd IFMBE European Medical and Biological Engineering Conference, November 20–25, 2005, Prague. - ISSN 1727-1983. - Prague. - 2005, Vol. 11, p. [1-4]
↑ Caldwell, I. R.; Nams, V. O. (2006). "A compass without a map: tortuosity and orientation of eastern painted turtles (Chrysemys picta picta) released in unfamiliar territory" (PDF). Canadian Journal of Zoology. 84 (8): 1129–1137. doi:10.1139/z06-102.
1 2 Benhamou, S (2004). "How to reliably estimate the tortuosity of an animal's path: straightness, sinuosity, or fractal dimension?". Journal of Theoretical Biology. 229 (2): 209–220. doi:10.1016/j.jtbi.2004.03.016. PMID 15207476.
↑ Gommes, C.J., Bons, A.-J., Blacher, S. Dunsmuir, J. and Tsou, A. (2009) Practical methods for measuring the tortuosity of porous materials from binary or gray-tone tomographic reconstructions. American Institute of Chemical Engineering Journal, 55, 2000-2012
↑ Espinoza-Andaluz, Mayken; Andersson, Martin; Sundén, Bengt (2017). "Computational time and domain size analysis of porous media flows using the lattice Boltzmann method". Computers & Mathematics with Applications. 74: 26–34. doi: 10.1016/j.camwa.2016.12.001 .
↑ E. Bullitt, G. Gerig, S. M. Pizer, Weili Lin, S. R. Aylward. Measuring tortuosity of the intracerebral vasculature from MRA images. IEEE Transactions on Medical Imaging, Volume 22, Issue 9, Sept. 2003, p. 1163 - 1171
↑ Gostovic, D., et al., Three-dimensional reconstruction of porous LSCF cathodes. Electrochemical and Solid State Letters, 2007. 10(12): p. B214-B217.
↑ Watanabe, Y.; Nakashima, Y. (2001). "Two-dimensional random walk program for the calculation of the tortuosity of porous media". Journal of Groundwater Hydrology. 43 (1): 13–22. Bibcode:2001JGHyd..43...13W. doi: 10.5917/jagh1987.43.13 .
↑ Gommes, C.J., Bons, A.-J., Blacher, S. Dunsmuir, J. and Tsou, A. (2009) Practical methods for measuring the tortuosity of porous materials from binary or gray-tone tomographic reconstructions. American Institute of Chemical Engineering Journal, 55, 2000-2012
↑ Espinoza Andaluz, M., Sundén, B., Andersson, M., & Yuan, J. (2014). Analysis of Porosity and Tortuosity in a 2D Selected Region of Solid Oxide Fuel Cell Cathode Using the Lattice Boltzmann Method. In Fuel Cell Seminar & Energy Exposition

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Fu, Jinlong; Thomas, Hywel R.; Li, Chenfeng (January 2021). "Tortuosity of porous media: Image analysis and physical simulation" (PDF). Earth-Science Reviews. 212: 103439. doi:10.1016/j.earscirev.2020.103439. S2CID 229386129.

[2] Richard M. Pearson. Optometric Grading Scales for use in everyday practice. Optometry Today, Vol. 43, No. 20, 2003 Archived 2012-04-04 at the Wayback Machine ISSN 0268-5485

[3] Hart, William E.; Goldbaum, Michael; Cote, Brad; Kube, Paul; Nelson, Mark R. (1999). "Automated measurement of retinal vascular tortuosity". International Journal of Medical Informatics. 53 (2–3): 239–252. doi:10.1016/s1386-5056(98)00163-4. PMID 10193892. Archived from the original on 2009-01-09.

[4] Enrico Grisan, Marco Foracchia, Alfredo Ruggeri. A novel method for automatic evaluation of retinal vessel tortuosity. Proceedings of the 25th Annual International Conference of the IEEE EMBS, Cancun, Mexico, 2003

[5] M. Mächler, Very smooth nonparametric curve estimation by penalizing change of curvature, Technical Report 71, ETH Zurich, May 1993

[6] Patasius, M.; Marozas, V.; Lukosevicius, A.; Jegelevicius, D.. Evaluation of tortuosity of eye blood vessels using the integral of square of derivative of curvature // EMBEC'05: proceedings of the 3rd IFMBE European Medical and Biological Engineering Conference, November 20–25, 2005, Prague. - ISSN 1727-1983. - Prague. - 2005, Vol. 11, p. [1-4]

[7] Caldwell, I. R.; Nams, V. O. (2006). "A compass without a map: tortuosity and orientation of eastern painted turtles (Chrysemys picta picta) released in unfamiliar territory" (PDF). Canadian Journal of Zoology. 84 (8): 1129–1137. doi:10.1139/z06-102.

[Benhamou_2004-8] 1 2 Benhamou, S (2004). "How to reliably estimate the tortuosity of an animal's path: straightness, sinuosity, or fractal dimension?". Journal of Theoretical Biology. 229 (2): 209–220. doi:10.1016/j.jtbi.2004.03.016. PMID 15207476.

[9] Gommes, C.J., Bons, A.-J., Blacher, S. Dunsmuir, J. and Tsou, A. (2009) Practical methods for measuring the tortuosity of porous materials from binary or gray-tone tomographic reconstructions. American Institute of Chemical Engineering Journal, 55, 2000-2012

[10] Espinoza-Andaluz, Mayken; Andersson, Martin; Sundén, Bengt (2017). "Computational time and domain size analysis of porous media flows using the lattice Boltzmann method". Computers & Mathematics with Applications. 74: 26–34. doi: 10.1016/j.camwa.2016.12.001 .

[11] E. Bullitt, G. Gerig, S. M. Pizer, Weili Lin, S. R. Aylward. Measuring tortuosity of the intracerebral vasculature from MRA images. IEEE Transactions on Medical Imaging, Volume 22, Issue 9, Sept. 2003, p. 1163 - 1171

[12] Gostovic, D., et al., Three-dimensional reconstruction of porous LSCF cathodes. Electrochemical and Solid State Letters, 2007. 10(12): p. B214-B217.

[13] Watanabe, Y.; Nakashima, Y. (2001). "Two-dimensional random walk program for the calculation of the tortuosity of porous media". Journal of Groundwater Hydrology. 43 (1): 13–22. Bibcode:2001JGHyd..43...13W. doi: 10.5917/jagh1987.43.13 .

[14] Gommes, C.J., Bons, A.-J., Blacher, S. Dunsmuir, J. and Tsou, A. (2009) Practical methods for measuring the tortuosity of porous materials from binary or gray-tone tomographic reconstructions. American Institute of Chemical Engineering Journal, 55, 2000-2012

[15] Espinoza Andaluz, M., Sundén, B., Andersson, M., & Yuan, J. (2014). Analysis of Porosity and Tortuosity in a 2D Selected Region of Solid Oxide Fuel Cell Cathode Using the Lattice Boltzmann Method. In Fuel Cell Seminar & Energy Exposition

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[11]

[12]

[13]

[14]

[15]