Tortuosity

Last updated
A tortuous river (meander of Nowitna River, Alaska) Nowitna river.jpg
A tortuous river (meander of Nowitna River, Alaska)

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]

Contents


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):

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:

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:

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]

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

Tortuosity calculation from an x-ray tomography reconstruction of a porous sandstone (the pores are shown): the color represents the shortest distance within the pore space from the left limit of the image to any point in the pores. Comparing this distance to the straight-line distance shows that the tortuosity is about 1.5 for this sample. It has been demonstrated that the tortuosity increases when the porosity decreases. Tortuosity via geodesic distance calculation.tif
Tortuosity calculation from an x-ray tomography reconstruction of a porous sandstone (the pores are shown): the color represents the shortest distance within the pore space from the left limit of the image to any point in the pores. Comparing this distance to the straight-line distance shows that the tortuosity is about 1.5 for this sample. It has been demonstrated that the tortuosity increases when the porosity decreases.

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. Intuitively, the curvature is the amount by which a curve deviates from being a straight line, or a surface deviates from being a plane.

<span class="mw-page-title-main">Helix</span> Space curve that winds around a line

A helix is a shape like a corkscrew. It is a type of smooth space curve with tangent lines at a constant angle to a fixed axis. Helices are important in biology, as the DNA molecule is formed as two intertwined helices, and many proteins have helical substructures, known as alpha helices. The word helix comes from the Greek word ἕλιξ, "twisted, curved". A "filled-in" helix – for example, a "spiral" (helical) ramp – is a surface called a helicoid.

<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.

Permeability in fluid mechanics and the Earth sciences 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.

<span class="mw-page-title-main">Soil mechanics</span> Branch of soil physics and applied mechanics that describes the behavior of soils

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, hydrology and soil physics.

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.

<span class="mw-page-title-main">Contact angle</span> The angle between a liquid–vapor interface and a solid surface

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.

<span class="mw-page-title-main">Curved space</span>

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 space. The Friedmann–Lemaître–Robertson–Walker metric is a curved metric which forms the current foundation for the description of the expansion of space and shape of the universe.

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 . 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.

<span class="mw-page-title-main">Horocycle</span> Curve whose normals converge asymptotically

In hyperbolic geometry, a horocycle, sometimes called an oricycle or limit circle, is a curve of constant curvature which converges asymptotically in both directions to a single ideal point, called the centre of the horocycle. The perpendicular geodesics through every point on a horocycle are limiting parallels, and also all converge asymptotically to the centre. It is the two-dimensional case of a horosphere.

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.

<span class="mw-page-title-main">Curve-shortening flow</span> Motion of a curve based on its curvature

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.

<span class="mw-page-title-main">Pore structure</span>

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

  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.
  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