Acessibilidade / Reportar erro

Sedimentación de partículas con distribución de tamaño fractal

Sedimentation of fractal size distribution particles

Resúmenes

Desde hace varios años, el modelo de fragmentación fractal ha atraido la atención de los investigadores, como un camino lógico para describir e interpretar distribuciones de partículas observadas. El análisis textural de un suelo ha mostrado ser muy importante, pues se utiliza para diagnosticar y predecir el funcionamiento y uso del mismo. Los métodos más populares para determinar la textura han sido los de sedimentación en agua utilizando el hidrómetro o la pipeta. Ambos tienen como objetivo encontrar la fracción de masa de partículas que se encuentran en suspensión a tiempos prefijados y relacionarla con los diámetros de las mismas. En este trabajo se ha desarrollado una nueva función potencial que relaciona la fracción de masa en suspensión con el tiempo de sedimentación. Utilizando la misma se puede determinar la dimensión fractal de fragmentación de una distribución de partículas en sedimentación. La nueva ecuación ha sido chequeada con datos propios obtenidos por el hidrómetro de Bouyoucos y otros publicados en la literatura, obtenidos mediante la pipeta de Robinson. El acuerdo logrado entre la teoría y los datos experimentales, mediante la técnica de regresión no lineal, ha sido excelente. Los valores de la dimensión fractal de fragmentación resultaron entre 2,404 y 2,512, para muestras de La Plata, Argentina, y entre 2,434 y 2,819 para los suelos de California, USA. El coeficiente de determinación, R², fue en todos los casos mayor que 0,9.

suelos; textura; fragmentación


Since several years the fractal fragmentation model has attracted the attention of researchers, as a logic way to describe and interprete observed particle size distributions. Textural analysis has shown to be very important because of its usefulness in the dignosis and inferences about soil functioning and use. Most popular methods of textural analysis employ sedimentation of particles in water using the hydrometer or the pipet. Both have the objective of determining the particle fraction remaining in suspension at predetermined time and to relate them with particle diameters. In the present work a new power law relationship between the mass fraction in suspension and the time was developed. Using this relationship it was possible to determine the fragmentation fractal dimension of a set of particles in sedimentation. The new equation has been checked with data obtained in this research by the Bouyoucos's hydrometer, and others published in the literature, using the Robinson's pipet method. The agreement between the model and the experimental data, using non linear regression, was excellent. Resulting fractal fragmentation dimensions ranged from 2.404 to 2.512, for samples from La Plata, Argentina, and between 2.434 and 2.819 for soils from California, USA. Determination coefficients, R², were always higher than 0.9.

soils; texture; fragmentation


MANEJO DE ÁGUA E SOLO

Sedimentación de partículas con distribución de tamaño fractal

Sedimentation of fractal size distribution particles

Roberto R. FilgueiraI; Yakov A. PachepskyII; Lidia L. FournierI; Mirta G. GarcíaIII; Guillermo O. SarliI

IFacultad de Ciencias Agrarias y Forestales, UNLP, Calles 60 y 119, CC 31, 1900 La Plata, Argentina. Fone: +54 221 4236758. Email: filgueira@infovia.com.ar; llfournier@infovia.com.ar; gsarli@ceres.agro.unlp.edu.ar

IIUS Dept. of Agriculture, Animal Waste Pathogen Laboratory, Bldg. 173, Rm. 203, BARC-EAST, Powder Mill Road, Beltsville, MD 20705. USA. Fone: +1 301 504 7468. Email: ypachepsky@anri.barc.usda.gov

IIIFacultad de Ciencias Agrarias y Forestales, UNLP, Calles 60 y 119, CC 31, 1900 La Plata, Argentina. Fone: +54 221 4236758. Email: mycsuelos@ceres.agro.unlp.edu.ar

RESUMEN

Desde hace varios años, el modelo de fragmentación fractal ha atraido la atención de los investigadores, como un camino lógico para describir e interpretar distribuciones de partículas observadas. El análisis textural de un suelo ha mostrado ser muy importante, pues se utiliza para diagnosticar y predecir el funcionamiento y uso del mismo. Los métodos más populares para determinar la textura han sido los de sedimentación en agua utilizando el hidrómetro o la pipeta. Ambos tienen como objetivo encontrar la fracción de masa de partículas que se encuentran en suspensión a tiempos prefijados y relacionarla con los diámetros de las mismas. En este trabajo se ha desarrollado una nueva función potencial que relaciona la fracción de masa en suspensión con el tiempo de sedimentación. Utilizando la misma se puede determinar la dimensión fractal de fragmentación de una distribución de partículas en sedimentación. La nueva ecuación ha sido chequeada con datos propios obtenidos por el hidrómetro de Bouyoucos y otros publicados en la literatura, obtenidos mediante la pipeta de Robinson. El acuerdo logrado entre la teoría y los datos experimentales, mediante la técnica de regresión no lineal, ha sido excelente. Los valores de la dimensión fractal de fragmentación resultaron entre 2,404 y 2,512, para muestras de La Plata, Argentina, y entre 2,434 y 2,819 para los suelos de California, USA. El coeficiente de determinación, R2, fue en todos los casos mayor que 0,9.

Palabras-clave: suelos, textura, fragmentación

ABSTRACT

Since several years the fractal fragmentation model has attracted the attention of researchers, as a logic way to describe and interprete observed particle size distributions. Textural analysis has shown to be very important because of its usefulness in the dignosis and inferences about soil functioning and use. Most popular methods of textural analysis employ sedimentation of particles in water using the hydrometer or the pipet. Both have the objective of determining the particle fraction remaining in suspension at predetermined time and to relate them with particle diameters. In the present work a new power law relationship between the mass fraction in suspension and the time was developed. Using this relationship it was possible to determine the fragmentation fractal dimension of a set of particles in sedimentation. The new equation has been checked with data obtained in this research by the Bouyoucos's hydrometer, and others published in the literature, using the Robinson's pipet method. The agreement between the model and the experimental data, using non linear regression, was excellent. Resulting fractal fragmentation dimensions ranged from 2.404 to 2.512, for samples from La Plata, Argentina, and between 2.434 and 2.819 for soils from California, USA. Determination coefficients, R2, were always higher than 0.9.

Key words: soils, texture, fragmentation

INTRODUCCION

Los suelos se han formado a partir de material primario fragmentado por meteorización, dando origen a una distribución de tamaño de partículas y agregados, que frecuentemente es interpretada como una función acumulativa, ya sea del número de partículas mayores que un cierto diámetro, o como la masa menor que un cierto diámetro. Una de las relaciones matemáticas más conocidas para este proceso es la siguiente (Mandelbrot, 1982; Turcotte, 1986):

N (r > R) = k R-D

(1)

donde N (r > R) es el número de objetos por unidad de volumen que tienen un radio mayor que un valor cualquiera R, mientras que k y D son constantes. Esta relación potencial define lo que se denomina un fractal probabilístico, también llamado fractal de fragmentación, siendo D la dimensión fractal de fragmentación (Baveye & Boast, 1998). En este caso, es esperable que el modelo sea válido dentro de un rango limitado de radios, dado que no es un fractal matemático. Turcotte (1986) propuso un modelo para el proceso de fragmentación de un cubo, que dividido hipotéticamente en ocho partes iguales puede fragmentarse en ocho cubos más pequeños. Cada cubo resultante puede a su vez dividirse en otros ocho, y así sucesivamente. La fragmentación del cubo tiene una probabilidad p de ocurrir, y se considera constante para todos los órdenes de tamaño. La probabilidad máxima sería p = 1, que implica una división en ocho cubos más chicos, y la mínima sería p = 1/8, para el desprendimiento de un solo cubo de los ocho posibles. En estas condiciones se dedujo la siguiente relación.

(2)

donde D es la dimensión fractal, que puede variar entre 0 y 3, y p la probabilidad de fragmentación.

Para partículas del suelo, Tyler & Wheatcraft (1992) sugirieron otra forma de expresar la Ec. 1, introduciendo la fracción de masa en vez del número de partículas. La ventaja consiste en que, experimentalmente, es más sencillo y preciso obtener fracciones de masa que número de partículas. Con este modelo, la masa acumulada de las partículas con un radio mayor que un valor R es:

(3)

donde Cm y lm son constantes.

La masa total de la distribución de estas partículas es:

(4)

La Ec. 3 puede ser normalizada, dividiendo por la Ec. 4, con lo que se obtiene la fracción de masa de partículas mayores que un valor R:

(5)

Para evaluar la constante lm suponemos un límite superior para el comportamiento fractal, r = RLsup. Para el caso de partículas del suelo, este límite sería el de la arena muy gruesa. Cuando r = RLsup se cumplirá que M(r > RLsup)/MT = 0, es decir, ya no habrá partículas por encima de este límite. De la Ec. 5 resulta entonces RLsup = lm.

En el caso del experimento de sedimentación es común obtener la fracción de masa de partículas que aún permanecen en la suspensión, es decir, M(r < RLsup)/MT. Podemos transformar fácilmente la Ec. 5, de la siguiente forma:

(6)

Este modelo ha sido aplicado también a muestras de suelos por otros autores (Bitelli et al., 1999; Filgueira et al., 2002a,b). Bitelli et al. (1999) estudiaron la distribución de tamaño de partículas de varios suelos mediante la técnica de tamizado en húmedo, la pipeta, y la difracción de la luz, en el rango 0,1 - 2000 mm. Estos autores identificaron en las distribuciones de tamaño de partículas, tres rangos separados por diferentes dimensiones fractales, que llamaron dominios arcilla, limo y arena, respectivamente. Los resultados mostraron que las dimensiones fractales obtenidas eran muy diferentes entre los dominios y aún dentro de cada uno de éstos. Filgueira et al. (2002b) investigaron el comportamiento de la distribución de tamaño de partículas de un suelo, a profundidades de 0 a 400 mm, mediante la técnica de Bouyoucos y con la aproximación fractal basada en la fracción de masa acumulada respecto del diámetro de las partículas (Ec. 6). Se corroboró que la distribución de masa de partículas acumulada versus diámetro sigue una ley potencial, en todos los casos. La dimensión fractal de fragmentación del dominio limo creció a medida que aumentó la profundidad del suelo.

Cuando un cuerpo cualquiera se mueve dentro de un fluido viscoso, el medio ejerce una fuerza de oposición al avance. De igual manera, si mantenemos el cuerpo quieto y es el fluido el que se mueve alrededor del objeto, éste siente una fuerza debida a esta corriente.

Stokes (1851) investigó el caso particular del movimiento de una esfera en un medio viscoso y encontró una expresión para la fuerza viscosa que se opone al movimiento, para el caso en que este movimiento (esfera respecto al fluido) sea tal que se conserven las condiciones de flujo laminar.

La expresión matemática de esta fuerza es:

Fv = 6 p h R v

(7)

donde h es el coeficiente absoluto de viscosidad del líquido, R es el radio de la esfera, y v la velocidad de la esfera respecto al fluido. La Ec. 7 es conocida como Ley de Stokes. El estudio de la caida de un cuerpo esférico en el interior de un líquido viscoso, en reposo, nos muestra que éste alcanza una velocidad de caída constante llamada velocidad límite. Planteando el sistema de fuerzas que actúan sobre el cuerpo en la caída, se deduce la siguiente expresión:

(8)

donde vlim es la velocidad límite, R el radio de la esfera, g la aceleración de la gravedad, h el coeficiente de viscosidad absoluta del líquido, rc y rL son las densidades de la esfera y del líquido, respectivamente. La velocidad límite se puede obtener fácilmente, midiendo una cierta distancia de caída h, tomando el tiempo empleado t. Luego, reemplazando vlim por h/t en la Ec. (8) y reagrupando nos queda:

(9)

Una aplicación muy difundida de la Ec. 9 es su uso en sedimentación, para separar partículas de diferentes radios, que están moviéndose en un líquido. Un ejemplo es el caso del estudio de la granulometría de un suelo; siendo los métodos más conocidos y aceptados los del hidrómetro de Bouyoucos y el de la pipeta de Robinson (Gee & Bauder, 1986). En estos casos se supone que las partículas del suelo son esferas con un radio equivalente al resultante de la Ec. 9. Kaddah (1974) realizó una investigación comparativa detallada de 24 suelos de California, utilizando las técnicas del hidrómetro y la pipeta, llegando a la conclusión de que tomando los recaudos necesarios ambos son equivalentes. Este autor publicó sus resultados en forma de tabla donde se daban los porcentajes de partículas que permanecían en suspensión acuosa, para diámetros de las mismas menores que 20, 10, 5, 2, y 1 mm. También desarrolló una metodología gráfica para convertir los datos del experimento de sedimentación a distintas temperaturas, en diámetros de partículas y porcentajes, y viceversa.

No hay antecedentes de que el modelo de fragmentación haya sido relacionado con la distribución de masa de partículas en función del tiempo de sedimentación en un líquido.

El objetivo de este trabajo fue encontrar una ecuación que relacione la distribución fractal de partículas esféricas con el tiempo, a una profundidad determinada, cuando las mismas están sedimentando.

MATERIAL Y METODOS

Modelo teórico

El análisis de la distribución de partículas del suelo se hace, generalmente, por sedimentación en agua. Esta técnica se basa, como ya mencionáramos, en la ley de Stokes, siendo los dos métodos más difundidos el del hidrómetro de Bouyoucos y el de la pipeta de Robinson (Gee & Bauder, 1986). En ambos casos lo que se busca medir es la concentración de partículas dispersas del suelo por unidad de volumen de la suspensión. Si suponemos que los tamaños de estas partículas esféricas están distribuidos según el modelo de fragmentación probabilística de Turcotte (1986), podemos combinar las ecuaciones 6 y 9 de la siguiente forma:

(10)

donde RLsup es el límite superior de radio de la esfera equivalente, h es la viscosidad absoluta del líquido de la suspensión, h es la distancia recorrida por las esferas en el tiempo t, g es la aceleración gravitatoria y rp y rL son las densidades de la esfera y del líquido, respectivamente. Fijado h, a es una constante.

La Ec. 10 es una función potencial y por ajuste con los datos experimentales obtenidos mediante la técnica de Bouyoucos o de la pipeta de Robinson se podría verificar si la distribución de partículas en sedimentación tiene un comportamiento fractal y obtener también la dimensión fractal de fragmentación.

Verificación experimental

Para la verificación del correcto funcionamiento del modelo se usaron dos tipos de datos experimentales: a) datos propios y b) datos tomados de la literatura.

Para obtener los datos propios, las muestras de suelos estudiadas se extrajeron de un Argiudol Típico de la Estación Experimental J. J. Hirschhorn dependiente de la Facultad de Ciencias Agrarias y Forestales de la UNLP, situada en La Plata, provincia de Buenos Aires (Latitud: 34º 55', longitud: 57º 50'), Argentina. La textura del horizonte superficial (Ap) es franco limosa. Se tomaron muestras con cilindros de 100 mm de diámetro y 200 mm de altura. Los sitios muestreados fueron cuatro, unos 100 m alejados entre sí. En cada sitio se tomaron dos muestras, las cuáles se llevaron al laboratorio, se dividieron a las profundidades de 0-100 mm y 100-200 mm y se dejaron secar a la temperatura ambiente. Posteriomente, se fragmentaron y mezclaron las dos muestras de cada sitio, para la misma profundidad, y se pasaron por un tamiz de 2 mm, para realizar el análisis textural. Este se realizó utilizando el método de sedimentación en agua, con el hidrómetro calibrado según Bouyoucos (1962). De las mezclas de muestras de suelos de las mismas características se utilizaron 50 g para las medidas. El tratamiento previo de las muestras, para la dispersión de las partículas, fue el convencional (Gee & Bauder, 1986). Los tiempos utilizados para medir la densidad de la suspensión acuosa fueron 40 s, y 2, 8, 15, 30, 60, 120, 240 y 1440 min. La determinación de las densidades según Bouyoucos fue hecha por duplicado, y se tomó el promedio de las medidas para la determinación de la dimensión fractal mediante el ajuste de una curva potencial que relacionó el diámetro de las partículas y el porcentaje acumulado que permanecía en la suspensión acuosa.

Los datos de la literatura fueron tomados de Kaddah (1974). Este autor estudió muestras superficiales de veinticuatro suelos con diferentes texturas, mediante los métodos de la pipeta de Robinson y del hidrómetro de Bouyoucos. De los resultados publicados por este autor se eligieron los 24 suelos cuya textura fuera medida por el primer método. La dependencia del porcentaje de partículas variando con el tiempo fue obtenida por interpolación usando la Fig. 1 de la publicación, de acuerdo a lo sugerido por el autor.


RESULTADOS Y DISCUSIÓN

En la Figura 1, se puede observar una representación típica de la fracción de masa de partículas acumulada en suspensión, obtenida en este estudio, versus el tiempo de sedimentación. Estos datos experimentales fueron ajustados mediante regresión no lineal con la Ec. 10 desarrollada en este trabajo. En la Figura 2 se muestra el mismo conjunto de datos experimentales, pero la fracción de masa se ha representado en función del diámetro equivalente de las partículas. También en este caso se hizo un ajuste de los datos con una función potencial, la Ec. 6, obtenida de Tyler & Wheatcraft (1992). Las dimensión fractal obtenida de ambas ecuaciones es conceptualmente la misma. Los valores de D obtenidos mediante la Ec. 6 dio, en efecto, resultados similares (no reportados en este trabajo), dentro de un 1,5% respecto a los obtenidos mediante la Ec. 10. En la Tabla 1 se pueden ver los valores de D obtenidos en este trabajo a partir de nueve datos medidos, para los suelos de Argentina, utilizando la Ec. 10, con los correspondientes coeficientes de determinación, R2. La variación de estos coeficientes, entre 0,927 y 0,992, da una idea del buen ajuste logrado por el modelo utilizado.


En las Figuras 3 y 4 se pueden ver representaciones gráficas de la fracción de masa versus tiempo y versus diámetro, para uno de los suelos de California (Kaddah, 1974). En este caso se usó la fracción limo, rango de 20 a 2 mm, con el fin de tener partículas pertenecientes a un sólo dominio de diámetros. En la Tabla 2 se pueden observar los valores de D estimados con la Ec. 10, y los coeficientes de determinación R2. El ajuste logrado en este caso con cuatro datos experimentales, fue también muy bueno. También en este caso los resultados del ajuste de los datos experimentales mediante la Ec. 6 se encuentran dentro un 1,5% de diferencia respecto a los obtenidos mediante la Ec. 10. Los valores de los coeficientes de ajuste, R2, variaron entre 0,900 y 0,996.



Se observa, en el caso del suelo de Argentina, que la dimensión fractal en la superficie, de 0 a 200 mm, toma valores entre 2,404 y 2,512. Para los suelos estudiados por Kaddah (1974) el rango de dimensiones fractales obtenidas varió entre 2,434 y 2,819. Tanto Wu et al. (1993) como Borkovec et al. (1993) sugirieron un valor universal D = 2,8 ± 0,1, lo cual está fuera del rango obtenido en este trabajo. Bitelli et al. (1999) reportaron, para lo que ellos denominaron el dominio limo, valores de D en el rango 1,728 - 2,792, siendo los valores más altos correspondientes a suelos con mayor contenido de arcilla, es decir suelos que presentan una distribución de partículas más finamente particionada. Filgueira et al. (2002) trabajando en el dominio limo, encontraron que la dimensión fractal de fragmentación variaba entre 2,410 y 2,679, siendo los valores correspondientes a muestras de mayor profundidad.

CONCLUSIONES

1. El modelo de sedimentación propuesto para un conjunto de partículas que tienen una distribución fractal de masa en función del diámetro, que obedecen la ley de Stokes, da también como resultado una función potencial de la fracción acumulada de partículas que permanecen en la suspensión versus el tiempo de sedimentación.

2. La ecuación resultante ajusta adecuadamente los datos experimentales obtenidos por el método de Bouyoucos o de la pipeta de Robinson.

3. La dimensión fractal obtenida mediante la nueva ecuación es la misma que resulta del modelo de fragmentación de Turcotte (1986). Los valores que toma esta dimensión es variable y dependiente del tipo de suelo.

LITERATURA CITADA

Protocolo 082.03 - 06/05/2003 - Aprovado em 23/10/2005

  • Baveye, P.; Boast, C. W. Concepts of "fractal" in soil science: demixing apples and oranges. Soil Science Society of America Journal, Madison, v.62, n.5, p.1469-1470, 1998.
  • Bitelli, M.; Campbell, G. S.; Flury, M.. Characterization of particle-size distribution in soil with a fragmentation model. Soil Science. Society of America Journal, Madison, v.63, n.4, p.782-788, 1999.
  • Borkovec, M.; Wu, Q.; Degovics, P.; Laggner, P.; Sticher, H.. Surface area and size distributions of soil particles. Colloids and surfaces A: Physicochemical and Engineering Aspects, Amsterdam, v.73, n.1, p.65-76, 1993.
  • Bouyoucos G. J. Hydrometer method improved for making particle size analysis of soils. Agronomy Journal, Madison, v.54, n.3, p.464-465, 1962.
  • Filgueira, R. R.; García, M. G.; Roggiero, M. F.; Cerisola, C.; Aragón, A.; Sarli, G. O; Soracco G. Caracterización de la distribución de tamaño de partículas en suelos con un modelo fractal. En: Congreso Argentino de la Ciencia del Suelo, 18, 2002, Resúmenes... Puerto Madryn: AACS, 2002a, v.1, p.13
  • Filgueira, R. R.; García, M. G.; Roggiero, M. F.; Cerisola, C. I.; Aragón, A.; Sarli, G. O. Uso del modelo fractal para caracterizar la distribución de tamaño de partículas en suelos. Ciencia del Suelo, Santa Rosa, v.20, n.2, p.114-117, 2002b.
  • Gee, G. W.; Bauder J.W. Particle-size analysis, in A. Klute (ed.) Methods of soil analysis. Part 1. Physical and mineralogical methods, 2.ed., Madison: ASA-SSSA, 1986, p.383-411. Agronomy Monograph No.9
  • Kaddah, M. T. The hydrometer method for detailed particle-size analysis: 1. Graphical interpretation of hydrometer readings and test of method. Soil Science, Baltimore, v.118, n.2, p.102-108, 1974.
  • Mandelbrot, B. B. The fractal geometry of nature. New York: Freeman & Co.:1982. 468p.
  • Stokes, G. G. On the effect of the lateral friction of fluids on the motion of pendulums. Transactions of Cambridge Philosophical Society, Cambridge, v.9, n.1, p.8-106, 1851.
  • Turcotte, D. L. Fractals and fragmentation. Journal of Geophysical Research, Washington, v.91, n.2, p.1921-1926, 1986.
  • Tyler, S. W.; Whitecraft, S. W. Fractal scaling of soil particle size distributions: analysis and limitations. Soil Science Society of America Journal, Madison, v.56, n.2, p.362-369, 1992.
  • Wu, Q.; Borkovec, M.; Sticher, H. On particle-size distributions in soils. Soil Science Society of America Journal, Madison, v.57, n.4, p.883-890, 1993.

Fechas de Publicación

  • Publicación en esta colección
    20 Oct 2006
  • Fecha del número
    Jun 2006

Histórico

  • Recibido
    06 Mayo 2003
  • Acepto
    23 Oct 2005
Departamento de Engenharia Agrícola - UFCG Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com