El coeficiente de descarga y la densidad beta

The Discharge Coefficient and the Beta Density

José Roberto Mercado*, Mauro Íñiguez, Pedro Guido, Javier Ramírez-Luna, Arturo González-Casillas
Instituto Mexicano de Tecnología del Agua,
*Autor de correspondencia.

Dr. José Roberto Mercado
Dr. Mauro Íñiguez
Dr. Pedro Guido
Dr. Javier Ramírez Luna
Dr. Arturo González Casillas

Instituto Mexicano de Tecnología del Agua (IMTA)
Paseo Cuauhnáhuac 8532, Colonia Progreso
62550 Jiutepec, Morelos, México
Teléfono: +52 (777) 3293 600 ]]> rmercado@tlaloc.imta.mx
mic@tlaloc.imta.mx
pedroguido@tlaloc.imta.mx
jramirez@tlaloc.imta.mx
arturo_gonzalez@tlaloc.imta.mx

Recibido: 08/12/2009
Aceptado: 26/06/2013

Resumen

Palabras clave: ecuaciones de descarga, autoafinidad, modelos de turbulencia, funciones de distribución, función de Kummer, multifractales, energía cinética.

Abstract

Discharge coefficient and turbulence intensity distribution are studied. With Torricelli's theorem and the approach of probability theory, flow discharge and discharge coefficient equation are derivate, following an unimodal Beta density function, renormalized, with two shape parameters. A multifractals model for the kinetic energy cascade in the turbulence was build, starting from the methods of Pearson and Kolmogorov. For turbulence intensity, with the first method, a Beta distribution was created; with the second, a power law. The multifractals model is completed, recognizing the structure function as a Kummer function. The compatibility between the two models are searched and so the identification of its parameters. It is found that the two shape parameters determine the cascade model resolution. Local dimension and dimension spectra are determine for the two states that produces Torricelli theorem. Redefining the structure function, resolution is defined by the water depth for the regime change. Analogously, different prototypes could be define, which we have call: the four experimentals, three channels, Kolmogorov, Kármán, Taylor, Verhulst (logistic), Cauchy-Manning, and Euclides (golden proportion). We conclude that the discharge coefficient is a renormalized Beta; turbulence intensities distribution is a Beta; Torricelli prototype results representative for the four experimentals and the Euclides, far away from the Gaussian distribution that is contained in von Karman model, meanwhile the Taylor's model yield the Dirac function.

Keywords: Discharge equations, self-similarity, turbulence models, density function, Kummer function, multifractals, kinetic energy.

Introducción

Existen muchos trabajos sobre la determinación del coeficiente de descarga para flujos a superficie libre y en diferentes documentos se proponen ecuaciones para flujo libre y ahogado. Existe además la deducción clásica de ecuaciones de descarga para vertedores de cresta ancha en canal rectangular (Henderson, 1966). En Sotelo (1999) se citan diferentes ecuaciones de descarga para vertedores por métodos empíricos. En dichas ecuaciones se observa la dificultad al pasar de un régimen a otro, con coeficientes empíricos que conserven su validez cuando ocurre dicho cambio de régimen. En general, para las estructuras de control y extracción no existe una ecuación de descarga que trabaje de forma adecuada en todos los regímenes de funcionamiento impuestos por la dinámica del flujo en los canales de riego. Una ecuación de descarga con invarianza para los diferentes regímenes es útil para su acoplamiento con la dinámica del flujo en el canal. De acuerdo con Baume (1992), la derivación de ecuaciones de descarga invariantes, bajo los diferentes regímenes de funcionamiento, ha sido raramente tratada.

La derivación de ecuaciones de descarga se ha abordado con diferentes enfoques empíricos y visión determinista, razón por la cual esas ecuaciones contienen coeficientes de descarga empíricos que dependen de varios parámetros (Henderson, 1966; Chow, 1959; Swamee et al., 1993; Baume, 1992).

El propósito de este trabajo es estudiar el coeficiente de descarga, desde un punto de vista experimental y teórico, y abordar su invarianza para diferentes vertedores y en diferentes regímenes de funcionamiento, posibilidad que se abre a través del enfoque estadístico. Kármán estudió la función de correlación a partir de la ecuación de Navier-Stokes; aproximando hasta la doble correlación y suprimiendo desde la triple correlación (correlación entre tres componentes) en adelante y bajo la hipótesis de la turbulencia isotrópica llega a una ecuación de Kummer, cuya transformada de Laplace inversa produce una distribución beta (De Kármán y Howarth, 1938). Nosotros invocamos las simetrías de la ecuación de Navier-Stokes para encontrar algunas propiedades de la función de correlación y completamos su descripción por medio de un proceso de Cantor generalizado, (Mercado, 2008). Ahora, con el coeficiente de descarga descrito a través de una distribución beta, y una transformada de Laplace directa, comprobamos la función de correlación, la cual, por el trabajo de Kármán, refleja la ecuación de Navier-Stokes como una aproximación hasta la doble correlación. Además, las simetrías de la ecuación de Navier-Stokes son independientes de la viscosidad, por lo que la función de correlación vale también para la ecuación de Navier-Stokes fraccional.

El esquema del trabajo se presenta en cinco apartados. En el primero se define el concepto de coeficiente de descarga y se desarrolla el método esbozado en esta introducción, la cual se sustenta en la ecuación de Navier-Stokes, pasa por la ecuación de Bernoulli y llega al teorema de Torricelli. En el segundo se aplica el teorema anterior y se elabora una forma covariante para el caudal adimensional generalizado; se reconoce la densidad beta y se muestra el cambio de la energía específica de la sección, en términos del tirante, hasta su forma beta. En el tercero se estudia la relación del coeficiente de descarga con la función Kummer, conocida como hipergeométrica confluente, y su descripción como función Fox; con el modelo (Mercado, 2008) se establece la relación entre los parámetros de forma de la beta y la resolución, para culminar la construcción del multifractal Torricelli, identificando la función de estructura como función Kummer. En la cuarta subsección se estudia la función de correlación a partir de las simetrías de la ecuación de Navier-Stokes hasta proponer la variable dual al tirante y la función de correlación con un exponente que contiene la información de los números característicos, y los números de Reynolds y Froude. En la última subsección se describen los prototipos y los números que los identifican, y se construyen tablas que recogen esa información; se estudian los números de Euler para diferentes vertedores y se identifican los valores de sus parámetros, después se enuncian las principales conclusiones.

Métodos

El coeficiente de descarga puede concebirse como un número de Euler que se expresa en función de las longitudes relativas características y de los números de Froude y Reynolds; es un objetivo primordial el esclarecer esa dependencia. Los números de Euler y de Froude se definen, respectivamente, por:

donde ρ, v, Δp, g y l son, respectivamente, densidad, velocidad, variación de presión, gravedad y una longitud característica.

En la referencia de Mercado et al. (2012) se expone una visión fraccional de la ecuación de Navier-Stokes y se considera su versión clásica para el presente trabajo, cuya forma vectorial es (2):

donde v, p, ρ, v, t y ɸ son velocidad, presión, densidad, viscosidad cinemática, tiempo y potencial escalar de la fuerza externa, como la gravedad. Integrando a lo largo de una línea de flujo entre dos cualesquiera de sus puntos, y bajo la hipótesis de flujo estacionario e inviscoso, surge la ecuación de D. Bernoulli (1738). Se trata de una ecuación de la energía y expresa la constancia de la misma. La expresión se identifica como carga total y se establece que la carga total es la misma a lo largo de cada línea de corriente (Rouse, 1946). La aplicación de la ecuación de D. Bernoulli para diferentes estados a lo largo de una línea de corriente puede expresarse como variaciones entre dos de estos estados, cualesquiera de ellos.

Desde un punto de vista histórico, la ecuación de D. Bernoulli es una generalización del teorema de Torricelli. Este teorema se obtiene al considerar que la superficie libre superior de un recipiente de agua está en el estado que se describe por la presión atmosférica, siendo constante y revalorada como p = 0; la carga piezométrica por h; y la velocidad v = 0; en la otra superficie, la inferior, el estado está dado por la presión atmosférica p = 0, la carga por h = 0, y la velocidad de salida por v. Por tanto, la variación de la energía es , entonces , luego . Este resultado en forma equivalente corresponde a la velocidad adquirida por una partícula fluida en el proceso de caída libre, que partiendo del reposo descienda la altura h bajo la acción única del campo gravitacional.

Coeficiente de descarga

Para las condiciones del teorema de Torricelli, con la caída del tirante H-h, y las diferencias de energía a lo largo de las líneas de corriente ΔE = ΔK + ΔU, con C como coeficiente de Coriolis, se tiene . Para pérdida de carga nula ΔE = 0, entonces , de tal manera que para un ancho unitario, A = h, el gasto es , la cual puede formularse como:

Dado que el coeficiente de Coriolis es mayor que 1, y puede llegar a acercarse al valor 2, su efecto puede interpretarse como un reforzamiento del número de Froude, visto como una intensificación de la inercia o un debilitamiento de la gravedad. En forma adimensional, el gasto, para un vertedor de cresta ancha en canal de sección rectangular, es (Ramírez, 1997). Que tiene la forma . Para definir el caudal máximo, según el cambio del régimen de funcionamiento, se aplica el principio de mínima energía en la ecuación adimensional y se obtiene la condición de cambio de régimen (la moda), con lo cual se puede definir:

Por lo que la ecuación debe formularse en términos de un gasto adimensional, adoptando una forma que generalice el gasto derivado del teorema de Torricell . Desde el punto de vista de la teoría de probabilidades, el factor (1–h/H)^b–1 se ve como proporcional al flujo (volumen) por unidad de tiempo y por unidad de área y, por tanto, representa la densidad de corriente de probabilidad; en tanto, el factor (h/H)^a–1, como proporcional al área, representa la función de riesgo, por lo que el producto de los dos designa la probabilidad por unidad de tiempo, y donde h/H es la cota superior del dominio de los valores que puede asumir la variable aleatoria en el intervalo (–∞, h/H). Aunque es igualmente probable definir al primero como el flujo y al segundo como el área, porque bastaría con hacer un cambio de variable y medir h/H de derecha a izquierda, y no al revés, como se acostumbra, lo que equivale al cambio ; en todo caso, uno de ellos representa al flujo de energía y el otro, al riesgo.

En Ramírez (1997) se formula un coeficiente de descarga de la siguiente forma:

Donde Q_*máx es el valor del caudal a descarga libre, y a y b son parámetros de forma, que dependen de las longitudes relativas características y de los números de Froude y Reynolds. Si L es la longitud de la cresta del vertedor, se observa que Q = C_dL Q_máx. El cambio de régimen (·)_cr viene dado por , la cual corresponde a la moda de la densidad de la distribución beta. Por tanto, se representa al coeficiente de descarga de forma general como:

Se trata de una densidad beta renormalizada , de tal manera que valga 1, por lo que queda normalizada por su máximo y no por el área debajo de la curva, como se acostumbra con la densidad beta. Ésta es la función de distribución más probable en el sentido del principio de mínima energía y queda determinada por la relación beta renormalizada (6), definida hasta el cambio de régimen.

En la sección de exposición de los prototipos se ilustra su aplicación al vertedor de cresta ancha, de cresta aguda, de compuerta, etcétera. Se conoce el comportamiento característico de la energía específica de una sección, que muestra un subrango dado por una función creciente de la energía con el tirante, el flujo subcrítico o lento; otro subrango con una función decreciente con el tirante es el flujo supercrítico o rápido; un punto que los separa y a la vez los une corresponde al tirante crítico, y al mínimo de la energía con respecto al tirante.

El tirante puede adimensionalizarse con el tirante crítico, definiendo h/h_c. Con el cambio de variable se contrae el intervalo de variación desde los reales positivos hasta el intervalo unitario y con el cambio de signo de la energía específica de la sección . La rama creciente de la función representa el subrango supercrítico y la rama decreciente el subrango subcrítico, separados por el punto máximo, representando el estado crítico. Así, el cambio de la densidad de probabilidad es .

Coeficiente de descarga y la función Kummer

La representación obtenida en la sección anterior para el coeficiente de descarga permite otra interpretación del mismo, que lo vincula con la ecuación diferencial de Kummer. El coeficiente de descarga es la transformada de Laplace inversa de la función Kummer, en tanto que es una solución de la ecuación diferencial homónima. O recíprocamente, la transformada de Laplace del coeficiente de descarga es solución de la ecuación diferencial de Kummer. Una sección transversal sobre el vertedor se cuadricula y sobre cada cuadro se toma un valor del tirante representativo, se adimensionaliza, se define la función completamente monótona ; con el proceso de Cantor generalizado y la discretización —entendiendo como s la variable dual al tirante, la cual también debe ser adimensional— se obtiene ₁F₁(–s;a,a+b) (Mercado, 2008). Pero ₁F₁(z;a,a+b) es una solución de la ecuación diferencial ordinaria zu"+(a+b–z)u'–au=0 y además , la cual es renormalizada y produce (7):

Por otra parte, esa función tiene también su representación como función Fox y, por tanto, como transformada inversa de Mellin (Metzler y Klafter, 2000):

en donde , y Γ(t) es la función gamma de Euler.

Multifractal Torricelli

En Mercado (2008) se presenta un modelo multifractal para la cascada de energía de la turbulencia partiendo de los métodos de Pearson y de Kolmogorov, simbolizando con q la llamada base de la resolución, siendo q > 1; además, se tienen, el espectro multifractal; , la dimensión fractal local; D_B, el máximo del espectro multifractal y medida del soporte; ε, la tasa de transferencia de la energía cinética; l₀, la longitud que corresponde a la dimensión de los vórtices iniciales; y n representando la enésima etapa del proceso de fractalización. Se encuentra que la intensidad de la turbulencia I_T, entendida como la media cuadrática del módulo de las variaciones de las velocidades (como sustituto del valor rms, el cual queda indeterminado para las distribuciones Levy, salvo la Gaussiana), puede describirse por , con u₀ = εl₀^1/3, , y D_T es la dimensión topológica, por tanto se obtiene (9):

Si se asume que la distribución de la intensidad de turbulencia sigue una distribución beta , se quiere calcular el sustituto del valor rms, para lo cual se requiere evaluar la integral ; pero su cuadrado es la transformada de Mellin de orden 2σ + 1 de la densidad beta, recordando que en tanto distribución se define como nula fuera del intervalo unitario, ; en términos de los símbolos de Pochhammer, el valor de la integral es , de tal manera que para una etapa bastante elevada del proceso de fractalización, ese valor se comporta como , por lo que la intensidad de la turbulencia resulta representada por su valor medio elevado a la potencia σ. Para que exista compatibilidad, la base de la resolución queda definida a partir de los dos parámetros de forma , y se denotará como k, más adelante. Además, el primer momento de la beta es la tangente en cero de la función Kummer, en tanto función generatriz de momentos; por tanto, si se construye un modelo multifractal, en donde la función de estructura sea la función Kummer, se puede afirmar que la transformada de Laplace inversa del coeficiente de descarga aporta la función de estructura. Para completar la descripción del modelo multifractal, se obtiene la dimensión local , por lo que el espectro resulta:

En particular, para los datos de Torricelli, el inverso de la resolución es y entonces la potencia debe ser .

Por tanto, del teorema de Torricelli y la cascada energética puede construirse la distribución de las intensidades de la turbulencia y renormalizarla para obtener el coeficiente de descarga, para que entonces su transformada de Laplace aporte la función generatriz de momentos, la cual es la función Kummer. Pero, además, para cada uno de los prototipos puede construirse el correspondiente multifractal con base en los parámetros de forma a y b, y la función de Kummer, de una manera similar al que hemos denominado de Torricelli. Una forma alternativa es definir la función de estructura como combinación lineal de la variable de reiteración s y la función de Kummer desplazada en la unidad g(s) = ₁F₁(–s;a,a+b)–1, así que g(s) satisface las tres condiciones g(0) = g'(∞) = 0, g(∞) = –1 y los dos coeficientes de la combinación lineal deben determinarse para cada caso (Liu et al., 2003).

Función de correlación

El presente camino constituye un recíproco al anterior, donde se describe la función de correlación de la velocidad y puede obtenerse la distribución de probabilidades asociada. En general, la función de correlación depende de la distancia entre los dos puntos y su orientación, pero como la ecuación de Navier-Stokes es invariante bajo rotaciones, la función de correlación depende sólo de la distancia entre los dos puntos y del tiempo (Olver, 1993).

Por otra parte, dentro de las simetrías de la ecuación de Navier-Stokes se incluyen las relativas a la escala, las que poseen como generadores a: G₁ = xδ_x + yδ_y + zδ_z + tδ_t y G₂ = tδ_t – uδ_u – vδ_v – wδ_w – 2pδ_p, siendo x, y, z, y t las coordenadas cartesianas usuales y el tiempo; y u, v, w, y p las coordenadas de la velocidad usuales y la presión. Para un canal con el eje x en la dirección de la corriente principal, las líneas del campo vectorial satisfacen o , así .

En general, como ya se dijo, la función de correlación será función de la distancia y del tiempo, pero a través de las transformaciones de escala dependerá de las longitudes características y de los números de Strouhal, Reynolds y Froude. Para el número de Strouhal resulta la dependencia de la velocidad proporcional, la cual, a través del factor de fricción, quedará en términos de los números de Reynolds y de Froude. De acuerdo con las longitudes características hay que tomar en cuenta los dos estados hidráulicos posibles: el liso y el rugoso. Cuando la rugosidad queda sumergida en la capa límite y se tiene el estado completamente liso, la escala de longitud de turbulencia está dada por la longitud intrínseca ; en tanto, en el estado rugoso está dada por la longitud asignada a la rugosidad k_s. Los números de Reynolds locales se dan por para el estado liso y para el estado rugoso. También para el estado rugoso, el número de Strouhal se da por y para el liso , así que para los dos estados ese número invariante se da por la velocidad proporcional, que a través del factor de fricción incluirá la dependencia de las longitudes características y los números de Reynolds y Froude, dependencia que finalmente será transferida a la potencia de la velocidad. Por tanto, la función de correlación se expresa a través de las variables adimensionales y , en la subcapa inercial; así (De Kármán, 1938; Chen, 2006).

Con base en los resultados experimentales, para el estado rugoso se admite representar la velocidad como distribución en potencia , con 0 < σ ≤ 1. Con los cambios de escala se tiene o bien repotenciados:

siendo m_σ > 1, el entero más grande, tal que el medido. Si en (11), con 0 < σ ≤ 1, se modifica la escala, así que y , entonces , por tanto en la ecuación (11) se tiene y entonces .

El proceso se basa en los criterios de Taylor y Kármán. Sea m(σ) = m, el mínimo σ admitido en la secuencia (o el máximo σ^–1). Bajo el criterio mk – 1 = k de Kármán, luego , con lo que , o mk – 1 > 0; pero , por tanto , y para Kármán, b = 1; por tanto, , con m > 2. Los valores propuestos por Kármán y también por Taylor son m_σ = 5 (De Kármán, 1938); sin embargo, esta normalización deja fuera importantes resultados reportados en Cheng-Lung (1991), pero puede escogerse m > 5 para que queden cobijados otros exponentes de rango rugoso, como el de Manning. Para ello, en la expresión de Kolmogorov , el término debe tender a cero cuando se suprima la turbulencia y se transforme en mono-fractal, con , por lo que , si D_B → 3. Sin embargo, el término es relativo, ya que puede ser absorbido por la resolución para sustituirlo por otro, así que la expresión de Kolmogorov sea y la condición D_B → 3, en cuyo caso m – 2 = 5, por tanto m = 7 es la cota inferior y el máximo σ^–1 admitido es 6 o , que corresponde al resultado experimental propuesto por Manning y está reportado en Cheng-Lung (1991). Está suficientemente documentado el vínculo entre la cascada energética de Kolmogorov y la distribución Gaussiana para las velocidades (Chen, 2006). Este vínculo se mantiene, porque si , entonces α = 2.0. Recíprocamente, si se toma b = 1, para mimetizar el factor (1 – h)^b–1, como se hace implícitamente en el modelo de cascada, y el inverso de la resolución es , se encuentra de , que , tal como se propuso en el modelo de cascada. Por otra parte, si el máximo aceptado para el inverso del exponente de la velocidad es y bajo el criterio de Kármán mk – 1 = k, luego , por tanto , y si m = 7, criterio de Taylor y Kármán, entonces , luego . Entonces, de , se infiere que , y como , entonces . Por tanto, bajo el criterio Kármám-Taylor: .

Por ejemplo, si se examina el exponente ligado a Lacey , conduce a 6 – 4 = 2; lo satisface. Por el contrario, el de Blasius , se tiene ; no satisface este criterio de clasificación. Por su parte, el de Manning produce y lo satisface en el límite. En conclusión, desde un punto de vista físico podría señalarse a la fracción como la fracción que separa los exponentes en los rangos lisos y rugosos. El intervalo experimental reportado por Agroskin , para el extremo , se tiene , y para es ; por tanto, todo el intervalo de Agroskin queda incluido y puede ser agregado al rango rugoso (Agroskin, 1980). Finalmente, la propuesta para la función de correlación es , siendo y , normalizada mediante: si α → 2^–, se tiene

σ(α) → 1; y para Taylor , que simboliza una delta de Dirac.

Resultados y discusión

En el presente apartado se describen los prototipos, los números que los identifican y se construyen los cuadros que muestran esa información. Además, se estudian los números de Euler para diferentes vertedores y se identifican los valores de sus parámetros.

Prototipos

Se consideran a continuación varios prototipos que provienen de muy diversas fuentes.

Para ellos, si se aportan los parámetros de forma, se obtiene el inverso de la resolución dada por . En seguida puede hallarse el exponente de la velocidad . Después se encuentra el índice Lévy resolviendo para α la ecuación .

Finalmente, el índice de ocupación espacial se obtiene de (Ramírez et al., 2009).

Al prototipo principal, los autores lo han llamado Torricelli, porque se obtiene del teorema homónimo. Este prototipo surge de los estados que dan origen al teorema de Torricelli y además se realiza como el vertedor de cresta ancha. Los exponentes que definen la densidad beta se determinan de a = 2, b = 3/2, luego , en seguida y . De la referencia Mercado (2008), reformulada como , se determina que α = 1.1843, siendo además .

Para un vertedor de cresta ancha, y de acuerdo con el teorema de Torricelli, si la carga total aguas arriba es H, para luego reducirse a h, la velocidad del flujo se determina como un punto material en caída libre, bajo la exclusiva acción de la gravedad, así que v = (2g(H – h))^1/2, y para el caudal en un ancho unitario es y en forma adimensional ; en particular, sobre la cresta del vertedor será , con . El inverso de la resolución es , por lo que coincide con el mismo valor para el de Torricelli; por tanto, el vertedor de cresta ancha está representado por el caso que se ha denominado Torricelli.

En el primer grupo se abordan los prototipos que provienen principalmente de resultados experimentales. Unos se sitúan dentro de los vertedores, en tanto que otros, como canales. Para los cuatro experimentales se determinan los parámetros de los vertedores: cresta ancha de arista viva, redondeado y de arista redondeada, triangular y Creager (cuadro 1), según la ecuación (6), y se anexa el valor de las variables al cambio de régimen de funcionamiento hidráulico. Los datos experimentales son presentados en Ramírez (1997).

En la figura 1 se muestran las gráficas del coeficiente de descarga en el eje vertical y de acuerdo con la ecuación (6); en tanto, en el eje horizontal se representa (1 – h/H)^1/2, y las gráficas tienen validez hasta el cambio de régimen. Se reconocen de izquierda a derecha de acuerdo con las columnas del cuadro 1; así, la primera línea continua corresponde a la primera columna o vertedor de cresta ancha y arista viva; la última, de cuadritos, al vertedor Creager. En correspondencia, en el cuadro 2 se muestran los valores de los parámetros de forma de la función beta unimodal, la resolución, el índice Lévy, el de ocupación espacial y el exponente de la velocidad, aunque se omite la columna 4 del vertedor de cresta ancha redondeada.

Los resultados experimentales citados en Agroskin (1980), que se ubican en el rango rugoso, producen lo que se tiene en el cuadro 4.

El resultado de Lacey (1930), citado por Cheng-Lung (1991), se originó de forma experimental, y es una de las potencias más altas dentro del rango rugoso. Se obtiene de la potencia de la velocidad , lo que produce , entonces α = 1.45, por lo que β(α,3) = 0.77011.

Orificio ahogado

Para el caso del flujo a través de un orificio, siendo C_c el coeficiente de contracción del chorro, un poco después de su salida y que está dado por el cociente del área del chorro sobre el área del orificio de salida (Rouse, 1946):

por tanto, con los valores de la variable aleatoria ligada al tirante y expresados por , se tiene .

En consecuencia, el inverso de la resolución es , luego produce α = 1.1262 con y β(α,3) = 0.70402.

Orificio horizontal ahogado

Para la descarga a través de un orificio horizontal ahogado se debe sustituir la proporción entre fronteras , siendo relativamente grande con respecto a , o sea la sustitución del ancho , entonces (Rouse, 1946):

por tanto, con se encuentra y, en consecuencia, se repiten los valores del caso anterior.

En el resultado anterior, si , se transforma en un vertedor de cresta aguda, y con , el caudal es (Rouse, 1946), desplazando es , y donde números F_d bajos implican grandes deflexiones gravitacionales:

así que , y y se ubica en el rango rugoso y cerca del de Manning. Además, se obtiene α = 1.1022; y β(α,3) = 0.69757.

Se recuerda que en este vertedor se presenta un decrecimiento de la presión piezométrica desde el dato, ubicado en la base plantilla, hasta el borde del vertedor de altura w, en donde se hace cero, acompañado con un crecimiento de la carga de velocidad; este borde representa un cambio brusco, un quiebre de la pendiente de la presión piezométrica. Luego dentro del chorro, un posterior crecimiento del tirante de presión hasta un máximo en el interior del chorro, para luego presentar un decrecimiento hasta la superficie exterior superior del chorro, que vuelve a ser cero, con el concomitante decrecimiento de la carga de velocidad, seguida de un crecimiento hasta la superficie exterior del borde superior del chorro (Rouse, 1946).

En este prototipo se presenta una diferencia con el de la abertura horizontal ahogada debido a que la presión en el exterior del chorro de aquel es la atmosférica, en tanto que para el presente se tiene una distribución hidrostática. Con la ecuación de continuidad y la de D. Bernoulli se encuentra (Rouse, 1946), luego:

Entonces y se reproduce el de Lacey. Se presenta un comportamiento del tirante de presión y de la velocidad similar al vertedor anterior con el cambio de la altura w por la abertura .

El segundo grupo de prototipos proviene de la teoría (algunos a partir de los parámetros de forma), otros del índice Lévy y algunos de la potencia.

Kolmogorov

Este prototipo surge de la cascada energética I_T (ecuación (9)), produciendo un modelo en potencia, que corresponde a una densidad beta B = h^a–1(1 – h)^b–1, con b = 1 y siendo el índice Lévy α = 2, que designa a la distribución Gaussiana, por lo que y de se obtiene . Además, . En tanto, el exponente de la velocidad es .

El tirante para cambio de régimen . El modelo de cascada se ilustra en la figura 2. En el eje vertical se representa la intensidad de la turbulencia I_T y en el horizontal h = q^–n. La curva de trazos es , ya que el coeficiente de la función beta es . Sobre esta curva de trazos, en la secuencia de pequeños círculos, se aprecia el proceso de fractalización que ocurre de derecha a izquierda, en donde la sucesión de la resolución se toma como de acuerdo con el planteamiento original de Kolmogorov, aunque también puede tomarse por ejemplo como . En tanto, la línea continua representa y sirve de referencia para el cambio del parámetro de forma hasta , y también para el cambio en la sucesión de la resolución desde hasta .

El de Kármán surge de la función de correlación, al establecer el inverso de la resolución como , así que en la distribución b = 5a. La potencia de la velocidad σ(α) = 1, por lo tanto coincide con la de Kolmogorov. El cociente de funciones gammas puede expresarse como una función Pochhammer y ésta se reduce a la identidad si , lo que equivale a α = 2, y entonces y . De forma recíproca, cerca del régimen laminar α → 2,

por lo que , lo cual significa la participación de la viscosidad.

Taylor

El modelo de Taylor brota de la función de correlación, siendo el inverso de la resolución. De ello se obtiene α = 2.1563, lo cual no es posible, porque α ≤ 2 y resulta σ → ∞.

Manning

Si el índice Lévy es α = 1, la distribución es la de Cauchy, luego el inverso de la resolución es k(1) = 1, entonces y b → 0; la potencia de la velocidad es y el índice de ocupación espacial es . De forma análoga, de acuerdo con los resultados de Manning, la potencia de la velocidad es , lo que produce k(α) = 1, entonces se obtiene α = 1 y . Este exponente de modo tradicional se ubica en el rango rugoso.

Euclides

Se origina en el arte y puede describirse a partir de la distribución, recordando la definición de la proporción de oro de Euclides: un segmento puede partirse en dos subsegmentos de longitud; el mayor es la media proporcional entre la totalidad y el menor , lo que da lugar a la solución , o en forma aproximada surge de dividir un segmento en cinco partes y destacar tres de ellas contra las dos restantes. Entonces, el inverso de la resolución es . El índice Lévy es tal que α = 1.1652. La potencia de la velocidad es . El índice de ocupación espacial es β = 0.71393.

Logística

Puede originarse desde cualquiera de dos aspectos: las funciones de distribución o de correlación. Desde la distribución, si a = 2 = b, entonces el inverso de la resolución es y el índice Lévy α = 1.2401. La potencia de la velocidad es y . En la figura 3 se muestra la gráfica de los inversos de la resolución con respecto al índice de estabilidad, donde se aprecian sus ubicaciones relativas y sus posiciones en referencia a los dos extremos marcados por α = 1 y α = 2. En el eje horizontal se representan los valores del índice Lévy α. La línea continua describe la gráfica de . En particular, se destaca el caso de Taylor, que se ubica fuera del intervalo señalado. Tanto los cuatro experimentales como el de Euclides de la proporción dorada se agrupan alrededor de la posición del de Torricelli y por tanto se asemejan al vertedor de cresta ancha; todos ellos, a su vez, se ubican más cerca del extremo 1 de la figura 3, por lo que se destaca su comportamiento turbulento.

La potencia de la velocidad se ilustra en la gráfica de la figura 4, en donde se señalan los dos extremos representados por α = 1, caso de Cauchy , y α = 2, caso de Gauss . En particular, para el de Taylor, la potencia diverge a infinito, así que está fuera del intervalo 0 < σ ≤ 1, como ya se señaló. La línea continua representa . En el eje horizontal se representa el índice Lévy.

La ecuación (6) permite representar de forma general las ecuaciones de descarga para diferentes vertedores. Se conjetura que dicha ecuación es válida para otros vertedores, lo que se ilustra encontrando la forma de los números de Euler para el orificio ahogado, el orificio horizontal ahogado y la compuerta. Los datos de los parámetros de los tres canales muestran la validez de las relaciones presentadas en el primer párrafo de la sección Prototipos. Los parámetros correspondientes al modelo teórico Torricelli son representativos de los cuatro experimentales, los que además se ubican claramente en régimen turbulento (figura 3), mientras que, por el contrario, el modelo teórico de Kolmogorov resulta más representativo para los tres canales, que se ubican más cerca del régimen laminar.

Conclusiones

1. El coeficiente de descarga adopta la forma covariante que se deriva del teorema de Torricelli o vertedor de cresta ancha.

]]> 2. La distribución de la intensidad de la turbulencia es la beta.

3. La función de correlación es la función Kummer.

4. Los exponentes de la velocidad del rango rugoso se extienden desde el de Manning hasta el de Kolmogorov o Kármán.

5. El coeficiente de descarga (6) refleja, a través de una aproximación, las ecuaciones de Navier-Stokes.

Referencias

AGROSKIN, I.I. Hidráulica. Tomo I. Torralba, V.J. (traductor). La Habana: Ministerio de Educación Superior, Instituto Superior de Ciencias Agropecuarias de La Habana, 1980, 521 pp.         [ Links ]

BAUME, J.P. Cemagref Mtp. Modélisation des ouvrages de type : déversoir, vanne, orifice, dans les modèles d'hydraulique à surface libre. Rapport interne. Montpellier: Cemagref, 1992.         [ Links ]

CHEN, W. A Speculative Study of 2/3-order Fraccional Laplacian Modeling of Turbulence: Some Thoughts and Conjectures. Chaos. Vol. 16, 2006, 023126.         [ Links ]

CHENG-LUNG, C. Unified Theory on Power Laws for Flow Resistance. J. of Hydraulic. Eng. Vol. 137, No. 12, 1991, pp. 1696-1699.         [ Links ]

CHOW, V.T. Open-channels hydraulics. New York: McGraw-Hill Book Company, 1959, pp. 33-37.         [ Links ]

DE KÁRMÁN, T. and HOWARTH, L. On the Statitical Theory of Isotrópic Turbulence. Proc. Roy. Soc. A, 164, London, 1938, pp. 192-215.         [ Links ]

HENDERSON, P.J. Open Channel Flow. New York: MacMillan Publishing Co., 1966.         [ Links ]

LIU, S.D., LIU, S.K., FU, Z.T., REN, K., and GUO, Y. The Most Intensive Fluctuation in Chaotic Time Series and Relativity Principle. Chaos, Solitons and Fractals. Vol. 15, 2003, pp. 627-630.         [ Links ]

MERCADO, J.R. La ecuación de Navier-Stokes y multifractales. Revista de Matemática: Teoría y Aplicaciones. Vol. 15, núm. 1, 2008, pp. 49-72.         [ Links ]

MERCADO, J.R., GUIDO, P., OJEDA, W., SÁNCHEZ, J., and OLVERA, E. Saint-Venant Fractional Equation and Hydraulic Gradient. Journal of Math. and System Science. Vol. 2, No. 8, 2012, pp. 494-503.         [ Links ]

METZLER, R., and KLAFTER, J. The Random Walk's Guide to Anomalous Diffusion: A Fractional Dynamics Approach. Phys. Rep. Vol. 339, 2000, pp. 1-77.         [ Links ]

OLVER, P.J. Applications Lie Groups to Differential Equations. New York: Springer-Verlag, 1993, 513 pp.         [ Links ]

RAMÍREZ, J. Modélisation des ouvrages frontaux et lateraux dans les canaux d'irrigation. These de Doctorat. Paris: L'École Nationale du Génie Rural, des Eaux et des Forêts (ENGREF), 1997.         [ Links ]

RAMÍREZ, J., MERCADO, J.R., PEREA, H., OLVERA, E., RUIZ, V., and ÍÑIGUEZ, M. Reynolds Equations and Water Distribution in Irrigation Canals. Hydraulic Engineering in México. Vol. XXIV, No. 3, July-September, 2009, pp. 121-130.         [ Links ]

ROUSE, H. Elementary Mechanics of Fluids. New York: Dover, Publ., 1946, pp. 376.         [ Links ]

SOTELO, G. Hidráulica General. Volumen I. México, D.F.: Editorial Limusa, 1999, pp. 126.         [ Links ]

SWAMEE, P.K., PATHAK, S.K., and ALI, M.S. Analysis of Rectangular Sluice Gate. Journal of Irrigation and Drainage Engineering. Volume 119, Issue 6, November, 1993, pp. 1026-1035.         [ Links ]