RESUMEN ]]> El flujo en medios porosos a través de contornos paralelos ha sido ampliamente investigado en las últimas cinco décadas. Pero son pocos los estudios con permeámetros convergentes, menos aún cuando el flujo no sigue la ley de Darcy. Este trabajo utiliza resultados experimentales publicados recientemente por otros autores donde afirmaron que los parámetros de la ecuación binómica de Forchheimer cambian en el flujo convergente. El objetivo consiste en estudiar esta ecuación en el flujo convergente y demostrar su validez bajo estas condiciones. Se obtiene un excelente ajuste gráfico de una serie experimental cuando se emplean los valores de todas las secciones del modelo plano convergente. Este resultado demuestra que los parámetros lineal y no lineal de un medio poroso homogéneo tienen un valor que no depende del tipo de flujo.

ABSTRACT

Flow behavior through porous media with parallel boundaries has been widely studied in the last five decades. Nevertheless, only few studies have been reported about flow behavior through porous media with converging boundaries, even less when flow does not obey Darcy's law. This work is based on experimental results through a convergent flow permeameter published by other authors recently. They concluded that the Forchheimer equation parameters change in converging flow. The objective here is to study non-Darcy steady converging flow and to show the validity of the Forchheimer equation under these flow conditions. An excellent graphical adjustment is obtained for different width ratios when all values of sections on the plane convergent model are used. This result shows that linear and non-linear parameters of a homogeneous porous media have one value not depending on flow conditions.

Keywords: porous media, converging flow, non-Darcy flow, Forchheimer law.

INTRODUCCIÓN

El estudio del flujo en medios porosos es importante en las ingenierías hidráulica y civil, así como en ramas relacionadas, como la geología, el petróleo y otras. Parece a simple vista un estudio conceptualmente simple, pero al entrar en detalles se observa una complejidad tal que no facilita su tratamiento (Kovacs 1981). Una revisión de las diferentes ecuaciones en estudios de flujo de filtración indica que la tradicional ley lineal de Darcy es ampliamente usada para flujo a baja velocidad (Porges y Hammer 2001) expresada como:

Donde V es la velocidad aparente o macroscópica; K es la conductividad hidráulica del medio poroso e I el gradiente hidráulico del flujo.

Frecuentemente la ecuación de Forchheimer (Domenico y Schwartz 1998) se emplea para relacionar la variación no lineal de la resistencia del flujo con su naturaleza, expresada como:

Donde a y b son coeficientes determinados por las propiedades del medio poroso y del fluido, conocidos como el parámetro lineal o de Darcy y el parámetro no lineal o no darciano, respectivamente.

La ecuación anterior de Forchheimer fue redefinida adicionándole un tercer término para considerar una condición transicional y un cuarto término para condición no estacionaria (Polubarinova-Kochina 1962).La expresión más general de Forchheimer es ligeramente mejor que la de los dos únicos términos, Nasser (1970), (citado por Porges y Hammer 2001) ya que según este autor resultan los mismos valores de gradiente por ambas expresiones.

Las investigaciones llevadas a cabo por Ward (1964) y Ahmed y Sunada (1969) requieren especial atención. Ward expresó dimensionalmente una ecuación válida tanto para flujo lineal como no lineal, expresada como:

donde: µ, ρ, g: viscosidad, densidad y aceleración gravitacional, respectivamente.
C_w y k son constantes del medio poroso.

Pérez Franco (1995) reconoce que estas constantes del medio poroso son verdaderas propiedades hidrogeológicas. Comparando la ecuación binómica de Forchheimer con su ecuación, Ward obtuvo las siguientes expresiones para los parámetros a y b:

Por su parte, Ahmed y Sunada (1969) partiendo de las ecuaciones generales de Navier-Stokes obtuvieron una ecuación idéntica a la propuesta por Ward. El flujo en medios porosos a través de contornos paralelos ha sido materia de amplia investigación en las ultimas cinco décadas (Ward 1964; Ahmed y Sunada 1969). En este caso de flujo para un caudal dado, a través de un medio poroso con estructura uniforme conocida, el gradiente hidráulico es constante así como la velocidad del flujo.

Cuando el flujo es convergente el gradiente es una función puntual, por lo que se requiere una expresión separada para estimar el gradiente hidráulico teniendo en cuenta el efecto de la variación de la velocidad (Nasser 1970). Este tipo de flujo resulta de interés ya que muchas situaciones de campo son convergentes por naturaleza, por ejemplo, el flujo a través del área adyacente a un pozo de bombeo, el flujo en presas de enrocamiento, el flujo a través de un banco de depósitos de rocas, el flujo a través de los filtros empleados en las plantas de purificación de agua, etc. Sin embargo, son pocos los estudios reportados en la literatura con permeámetros convergentes, menos aún cuando el flujo no obedece a la ley de Darcy.

Nasser (1970) supone que el parámetro de Darcy es el mismo en ambos flujos, sea paralelo o convergente, y la convergencia sólo afecta al parámetro no lineal de la ecuación binómica de Forchheimer (Porges y Hammer 2001). Thiruvengadam y Kumar (1997), basados en experimentos con roca y esferas, han obtenido ecuaciones empíricas para a y b, válidas solamente en el rango limitado de su experimentación.

Reddy y Rao (2006), sostienen que los valores de a y b varían cuando el flujo es convergente, no solamente a lo largo de la dirección del flujo sino también en la dirección radial o transversal al mismo. El tratamiento dado en la literatura a los parámetros de la ecuación binómica de Forchheimer cuando el flujo es convergente motiva el objetivo de este trabajo, el cual consiste en la evaluación del comportamiento del flujo no lineal a través de un medio poroso con contornos* convergentes y demostrar la validez de la ecuación de Forchheimer bajo estas condiciones de flujo, donde sus parámetros a y b son funciones únicas del medio poroso, y no de las condiciones del flujo que ocurre a través del mismo, sea éste paralelo o convergente.

GRADIENTE HIDRÁULICO EN FLUJO CONVERGENTE

Thiruvengadam y Kumar (1997), basados en la experimentación en cuatro muestras de roca y tres de esferas de vidrio y en el principio de conservación de la energía, obtuvieron una expresión exponencial para el gradiente hidráulico, dada por:

Donde I es el gradiente, x es la distancia a lo largo del flujo. Pi y m son determinados por ajuste empírico, lo que hace que la ecuación (5) sea válida para el rango limitado de condiciones en las cuales se realizaron sus experimentos.

Una fundamentación para determinar el gradiente hidráulico entre dos puntos cualesquiera en la dirección principal de un flujo no lineal, a través de un medio poroso con contornos convergentes se basa en adoptar las tres siguientes hipótesis:

· Medio poroso homogéneo y uniforme.

· Las condiciones de entrada garantizan un flujo radial con contornos convergentes.

Según el esquema mostrado en la figura 1, P1 y P2 son puntos localizados en la línea de corriente central del permeámetro convergente, en los radios R1 y R2 desde O, respectivamente. Se designan: hL la pérdida de carga entre P1 y P2 para un caudal Q, W el ancho de flujo y θ el ángulo formado en O por el permeámetro. Sea dh_i la pérdida de carga en una franja elemental de ancho dR.

El caudal a través de un radio Ri será: Q = V_i R_i θ

donde R_i θ W es el área normal al flujo. La ecuación (2) de Forchheimer aplicada en la franja elemental ABCD es:

La ecuación de continuidad expresa la velocidad en la sección V_i como:

La expresión (8) es similar en forma a la de Nasser (1970). Los parámetros a y b de la ecuación binómica pueden obtenerse como parámetros de ajuste de la función dada en (8), a partir de los valores medidos en la experimentación bajo flujo convergente, vale decir, de la variación de las cargas medidas entre la entrada R1 y la sección R_i, o sea, h₁ h_i , y del caudal de circulación Q.

Es necesario vincular la expresión (8) con la ecuación de continuidad (7) para expresarla en función de la velocidad en la entrada del permeámetro V₁ cuando ésta es la información experimental dada o procesada, obteniéndose al final la ecuación (9a):

o también,

Donde a y b serían los parámetros de ajuste de la función (9b), los cuales dependen de la velocidad en la entrada del permeámetro convergente V₁ ; de la variación de las cargas medidas en la sección h_i y en la entrada h₁ y de la relación F₁ entre los radios de la sección evaluada y la de entrada. F_i es llamado factor de la sección:

Obtenidos los parámetros a y b, se calcula ahora el gradiente puntual I_i en la sección R_i con la velocidad de la sección V_i a través de la ecuación binómica (10):

De esta forma se evita emplear el llamado gradiente hidráulico medio, y como consecuencia llevar a cabo un manejo inadecuado de la información experimental, el cual ha sido realizado por algunos autores en sus experiencias con flujo convergente en medios porosos (Venkataraman y Rao 2000; Reddy y Rao 2006).

El gradiente hidráulico medio I_m, antes referido, se define como:

Una excepción de lo anterior es la desarrollada por Thiruvengadam y Kumar (1997), que obtienen una ecuación experimental para el gradiente hidráulico en función de la velocidad, expresión (5) ya comentada.

El flujo no lineal radial convergente utilizando la ecuación de Dupuit-Forcheimer ha sido ampliamente estudiado en la hidráulica de pozos en los trabajos de Pérez Franco (1995) y otros. Por lo que el caso bajo estudio puede ser analizado por el mismo procedimiento. El gradiente en la sección que interviene en la ecuación binómica (10) no es posible medirlo en la experimentación, por lo que deberá ser calculado a través de ella una vez conocidos sus parámetros a y b.

Por otra parte, si se divide la expresión (8) entre (R₁- R_i) se obtiene el gradiente medio definido en (11). Sustituyendo en ella el caudal Q en función de la velocidad de entrada, auxiliado con la ecuación de continuidad (7), se obtiene la ecuación de Forchheimer en función del gradiente medio y de la velocidad de entrada, expresada en la ecuación (12):

Los términos en corchetes de (12) son los llamados «factores de convergencia» por Venkataraman y Rao (2000), los cuales reconocen que para un caudal dado estos factores son puramente geométricos ya que dependen sólo de la relación entre los radios de los puntos piezométricos.

Por su parte, Reddy y Rao (2006) afirman que la ecuación de Forchheimer cambia sus parámetros a y b en el flujo convergente y sin embargo no hacen mención de estos «factores de convergencia».

MANEJO DE LA ECUACIÓN DE FORCHHEIMER EN FLUJO CONVERGENTE

La tabla 1 muestra los valores puntuales observados en la experimentación realizada por Reddy y Rao (2006) del suelo de 3.25 mm, donde los valores de velocidad se muestran en la columna B y la relación de gradiente y velocidad en la columna C, tomados de los gráficos publicados por los autores para cada sección F_i. Ellos sostienen que los valores de a y b varían cuando el flujo es convergente, no solamente a lo largo de la dirección del flujo sino también en la dirección radial o transversal al mismo.

El gradiente medio mostrado en la columna D de la tabla 1 es obtenido a partir de las columnas B y C. Los autores antes referenciados realizan sus experimentos en un modelo plano vertical. La figura 1 ilustra un esquema del permeámetro empleado.

El flujo comienza en un radio R₁ bajo una carga hidráulica constante H₁ y va atravesando ocho secciones inferiores de radio Ri donde el factor F₁ varía en dichas secciones de la siguiente manera: 1.44; 1.55; 1.83; 2.02 y 2.25.

Obtención de los parámetros a y b según las cargas piezométricas y las velocidades medidas en la entrada.

Como las cargas hidráulicas medidas en los piezómetros no son dadas (no se publican los datos originales básicos de los experimentos), se estima la diferencia de carga entre la entrada y la sección i, Δh_i = (h₁ h_i), como sigue: Δh_i = I_m. (R1Ri), donde I_m es el gradiente medio dado en la columna D de la tabla 1, obtenido de los valores puntuales medidos.

Se considera un valor de R₁ de 1282 mm y un ángulo è de convergencia de 40.95 grados. Los valores así estimados de Δh_i se muestran en la Columna E de la tabla 1.

Resultados obtenidos con los valores de la experimentación en todas las secciones del perméametro

La expresión (9b) dividida entre el término de a: (V₁*R₁* ln ( Fi )) da:

La función (13) resulta una recta en un gráfico donde el eje Y sea el término a su izquierda y el eje X el que acompaña al parámetro b.

En la figura 2 se presenta este gráfico con todos los valores de la experimentación de Reddy y Rao (2006), correspondiente con las cinco secciones y ocho velocidades analizadas. Los valores obtenidos de la recta de ajuste confirman la justeza de lo tratado: a = 0.0259, b = 0.0089 con una correlación de 0.983.

Cálculo del Gradiente hidráulico

El gradiente hidráulico en cualquier sección se halla por la ecuación binómica (10), una vez determinados sus parámetros a y b, empleando la velocidad correspondiente a la sección. En la tabla 2, columna D, se presentan los valores calculados del gradiente hidráulico puntual en todas las secciones en dependencia de la velocidad en la sección, calculada por la ecuación de continuidad (7) con la velocidad de entrada de dato, mostradas en las columna C y B, respectivamente.

]]>

Variación de ( Ii / Vi ) versus Vi

Es evidente que el ajuste en la sección, con los gradientes y velocidades puntuales, dará siempre la misma recta de ajuste en todas ellas, donde como se sabe su intercepto con el eje vertical es el parámetro a y la pendiente es b, ya determinados antes.

Se demuestra en este trabajo que los parámetros a y b son únicos en un medio poroso homogéneo, relacionados con sus verdaderas propiedades hidrogeológicas (ecuación 4), y no deben por tanto variar atendiendo a que el flujo sea paralelo o convergente.

Resultados obtenidos por secciones

Se muestra a manera de ejemplo en la figura 3 la correlación en una sección.

CONCLUSIONES ]]> · El trabajo presenta un procedimiento para el manejo de la ecuación binómica de Forchheimer cuando el flujo es convergente en un medio poroso. El empleo de esta ecuación con el gradiente medido en la experimentación, llamado gradiente medio, y a base de la velocidad en la entrada del permeámetro, independiente de la sección analizada, ha conducido a poner en duda su cumplimiento cuando el flujo es convergente por algunos autores en la literatura científica.

· Se obtiene un excelente ajuste en el gráfico que permite determinar los parámetros de la ecuación binómica, cuando se emplean los datos experimentales de todas las secciones de un perméametro plano convergente. Este resultado demuestra que el medio poroso homogéneo tiene un único valor de sus parámetros a y b, por lo que sus valores no dependen del tipo de flujo que ocurra a través del mismo.

· Los resultados obtenidos de los parámetros a y b, con todos los valores de una experimentación bajo flujo convergente, demuestran la validez de la ecuación de Forchheimer donde sus parámetros son únicos en un medio poroso homogéneo. Además, alertan la necesidad de manejar con cuidado la ecuación binómica en flujo convergente cuando se opera con gradientes medios y empleando sólo la velocidad de entrada, lo cual es lo que obliga a introducir en ella los llamados factores de convergencia.

· Se puede apreciar que a medida que aumenta la separación entre las secciones o lo que es lo mismo, aumente la proporción entre ellas, también aumentará la diferencia entre el gradiente medio entre dichas secciones, dado en la columna D de la tabla 1, y el gradiente teórico correspondiente a la sección, en la columna D de la tabla 2, lo que es la causa evidente de la dispersión obtenida en los coeficientes a y b de la ecuación binómica y que los autores Venkataraman y Rao (2000), Reddy y Rao(2006), entre otros, atribuyen al flujo convergente.

AGRADECIMIENTOS

Los autores agradecen a los ingenieros Eric Cabrera Estupiñán y Yoel Martínez González, profesores de nuestro centro laboral, gracias a los cuales fue posible llevar a cabo la manipulación de los datos experimentales obtenidos por otros autores.

Ahmed, N. y Sunada, D.K. «Non linear flow in porous media». Journal of Hydraulics Division, ASCE, 1969, Vol.95, No. HY3, pp.1847-1857, (1969).

Domenico, P.A. y Schwartz, W. «Physical and chemical hydrogeology». Second Edition, John Wiley, (1998).

Kovacs, G. «Seepage hydraulics». Developments in Water Science, ESPC, New York, (1981).

Pérez Franco, D. «La explotación del agua subterránea. Un nuevo enfoque». Editorial Científico-Técnica. La Habana, (1995).

Polubarinova-Kochina, P.Y. «Theory of the motion of ground water», J. M. De Weist, translator, Princeton University Press, Princeton, N.J, (1962).

Porges, R.E. y Hammer, M.J. «The compendium of hydrogeology». National Ground Water Association, John Wiley & Sons, USA, (2001).

Reddy,N.B.P. y Rao, P.R.M. . «Effect of convergence on non linear flow in porous media». Journal of Hydraulic Engineering, ASCE, vol.132, no.HY4, pp.420- 427, (2006).

Thiruvengadam, M. y Pradip Kumar,G.N. «Validity of Forchheimer equation in radial flow through coarse granular media» Journal of Engineering Mechanics, ASCE, vol.123, no. HY 7, pp. 696-705, (1997).

Venkataraman, P. y Rao, P. R. M. «Validation of Forchheimer's law for flow through porous media with converging boundaries». Journal of Hydraulic Engineering, ASCE, vol.126, no. HY1, pp. 6371, (2000).

Ward,J.C. «Turbulent flow in porous media». Journal of Hydraulics Division, ASCE, 1964, Vol.90, Núm. HY5, pp.1-13, (1964).