Artículos técnicos

 

Longitud de mezcla hiperbólica, perfil de las velocidades y ley de resistencia al flujo

 

Hyperbolic mixture length, velocity profile, and flow resistance law

 

Carlos Fuentes

 

Universidad Autónoma de Querétaro, México.

 

Dirección institucional del autor

Dr. Carlos Fuentes

Facultad de Ingeniería
Universidad Autónoma de Querétaro
C.U., Cerro de las Campanas
76010 Santiago de Querétaro, Querétaro, México
Teléfono: +52 (442) 1921 200, extensión 6036
Fax: +52 (442) 1921 200, extensión 6006
cfuentes@uaq.mx

 

Recibido: 25/02/09
Aprobado: 08/06/10

 

Resumen

Considerando que el esfuerzo cortante total en el flujo del agua resulta de la suma de un esfuerzo viscoso y uno turbulento, se introduce una longitud de mezcla hiperbólica en la hipótesis de Newton para describir la relación entre el esfuerzo cortante total y la derivada del perfil medio de las velocidades. La longitud hiperbólica se reduce a la longitud de mezcla de Prandtl, que es proporcional a la distancia a la pared cuando ésta es grande. La existencia de un punto de inflexión en el perfil de las velocidades induce un valor teórico de la constante de von Kármán muy cercano a su valor experimental. Con este perfil se deducen fórmulas para las velocidades medias del fluido en un tubo y en un canal rectangular; de éstas se deducen las leyes de resistencia al flujo. La ley de resistencia en el tubo llevada a la forma de Darcy-Weisbach indica que el factor de fricción presenta los comportamientos extremos de Poiseuille y de Prandtl. La rugosidad de la pared a altos números de Reynolds se introduce siguiendo la idea de Colebrook y White, al adicionar a la escala de la turbulencia otra proporcional al tamaño de los granos de arena de Nikuradse; el factor de fricción tiene el mismo comportamiento que el mostrado en el diagrama de Moody. La escala dependerá del número de Reynolds en régimen transitorio de acuerdo con los resultados experimentales de Nikuradse, y su dependencia a pequeños y moderados valores del número de Reynolds y de la rugosidad es un sujeto de estudio.

Palabras clave: un valor teórico de la constante de von Kármán, efectos de la rugosidad de la pared.

 

Abstract

Considering that in water flow total shear stress is the sum of a viscous shear stress and a turbulent shear stress, a hyperbolic mixture length is introduced in the Newton hypothesis to describe the relationship between total shear stress and the mean velocity profile derivative. The hyperbolic length is reduced to the Prandtl length which is proportional to the wall distance when the latter is large. The existence of an inflexion point in the velocity profile induces a theoretical value of the von Kármán constant very close to its experimental value. From this, profile formulas for average fluid velocities in a pipe and rectangular channel are deduced, and from these latter, resistance to flow laws are established. The pipe resistance law, expressed in a Darcy-Weisbach form, indicates that the friction factor shows the extreme behaviors of Poiseuille and Prandtl laws. The wall roughness for high Reynolds numbers is introduced following the Colebrook and White idea, adding to the turbulence scale another one proportional to the Nikuradse sand grains size; the friction factor has the same behavior as that of the Moody diagram. The last scale will depend on the Reynolds number in transient regime in agreement with the experimental results of Nikuradse and its dependency with respect to small and moderate Reynolds number values, and the wall roughness is a study subject.

Keywords: theoretical value of the von Kármán constant, wall roughness effects.

 

Introducción

En hidráulica de canales y tuberías, hidrología superficial, ingeniería de riego y drenaje, entre otras disciplinas, es de fundamental importancia la utilización de una ley de resistencia al flujo para estimar la fuerza de fricción. La ley de resistencia relaciona la pendiente de fricción (J) con la velocidad media (U) y el radio hidráulico (R_h) (Du Buat, 1786; Trueba-Coronel, 1954).

En flujo a superficie libre es común utilizar la ley en potencia U = cR^b_hJ^d, que incluye las leyes extremas en flujo turbulento de Chézy (1776) (b = 1/2, d = 1/2) y en flujo laminar de Poiseuille (1840, 1841) (b = 2, d = 1) originalmente obtenida en tubos capilares. En un estudio sobre el acoplamiento de las ecuaciones de Saint-Venant (1871) y de Richards (1931) en el riego por gravedad en un canal rectangular muy ancho, donde R_h = h, siendo h el tirante de agua, Fuentes et al. (2004) establecen que la singularidad presente en el tiempo inicial requiere que los exponentes estén relacionados por b = 3d - 1. Esta relación es exactamente satisfecha por las leyes de Chézy y Poiseuille, aproximadamente por la ley de Hazen-Williams (b = 0.63, d = 0.54) utilizada en el diseño de redes de tuberías (King et al., 1952) y claramente la ley de Manning (1891) (b = 2/3, d = 1/2) es excluida por la singularidad.

La ley de resistencia en potencia que satisface la singularidad, haciendo intervenir el coeficiente de viscosidad cinemática (v) y la aceleración gravitacional (g), puede ser escrita como U = kν^1-2dh^3d-1(gJ)^d, donde 1/2 ≤ d ≤ 1 y k es un coeficiente empírico; la expresión para el caudal unitario es q = Uh = kν(h³gJ/v²)^d. Una posible limitación en la aplicación de esta ley es la dependencia de d con respecto al número de Reynolds (1883) (Re), d → 1 cuando Re → 0 y d → 1/2 cuando Re → ∞. Por ejemplo, en el riego por gravedad es usual aportar un caudal constante, la presencia de un sumidero (infiltración del agua en el suelo) induce que el caudal sea nulo en el frente de onda, y puesto que Re ∝ q, es claro que Re varía desde cero hasta un valor máximo proporcional al caudal unitario de aporte.

Existen investigaciones experimentales y teóricas sobre la dependencia de este exponente con respecto al número de Reynolds. En efecto, introduciendo el esfuerzo cortante en la pared τ_p = ρgR_hJ se definen las escalas de velocidad (u_*) y longitud (λ_c) siguientes:

donde ρ es la densidad del agua y u_* es denominada la velocidad de fricción.

Con las escalas de la ecuación (1), la ley en potencia adquiere la forma U/u_* = k(h/λ_c)^δ, donde 0 ≤ δ = 2d - 1 ≤ 1. Ésta puede ser deducida aceptando que el perfil de las velocidades medias en potencia u(y) = k_δu_*(y/λ_c)^δ es válido en todo el dominio 0 ≤ y ≤ h, a partir de la expresión que define el caudal unitario:

con k_δ= k(1 + δ) y en donde se ha agregado el subíndice C, para denotar que se trata de la velocidad media en un canal rectangular muy ancho. Para el caudal en un tubo se tiene:

De acuerdo con las experiencias en tubos llevadas a cabo por Nikuradse, los parámetros k_δ y δ varían con el número de Reynolds (Re). En tubos lisos y para 5 × 10³ < Re < 10⁵, Blasius propone k_δ = 8.74 y δ = 1/7 (Sotelo-Ávila, 1974). Barenblatt (1993) y Barenblatt y Prostokishin (1993) argumentan que el comportamiento asintótico (Re → ∞) de k_δ y δ puede ser k_δ = k₁/ε + k₂ + k₃ε + ... y δ = δ₁ε + δ₂ε² + δ₃ε³ + ..., donde ε = 1/h(Re) y sugieren, a partir de los datos de Nikuradse, los valores empíricos k₁ = 0.577, k₂ = 2.50 y δ₁ = 1.5, y nulos el resto de los coeficientes. Zagarola y Smits (1997), con nuevos datos sugieren k₁ = 0.7053, k₂ = 0.3055, δ₁ = 1.085, δ₂ = 6.535, y nulos el resto de los coeficientes, para el rango del número de Reynolds de 3.1 × 10⁴ a 3.5 × 10⁷. La dependencia de δ con respecto al número Reynolds pequeños no es bien conocida.

La aplicabilidad de la ley de resistencia en potencia es posiblemente limitada, ya que la variación del exponente en la vecindad del régimen laminar no es plenamente conocida. Ante esta dificultad se investiga la posibilidad de extender la ley de resistencia que resulta del perfil logarítmico de las velocidades de Prandtl a números de Reynolds pequeños.

En el flujo turbulento cercano a la pared se distinguen tres regiones principales (e.g. Dally y Harleman, 1969; White, 1979): a) la región interior, donde los esfuerzos viscosos son dominantes, también denominada subcapa viscosa; b) la región exterior, donde los esfuerzos turbulentos son dominantes, y c) la región intermedia o capa de solape, donde ambos tipos de esfuerzos son importantes. En la supcapa viscosa, la velocidad es una función básicamente lineal de la distancia a la pared: u/u_* ≃ y/λ_c cuando 0 < y/λ_c < 4. La región exterior sigue una ley de deflexión (u_máx - u)/u_* = G(y/δ), donde u_máx es la velocidad máxima ubicada en el espesor de la capa límite δ, generalmente asimilada al radio R en un tubo. En la capa intermedia se ha argumentado que la velocidad sigue una ley logarítmica (Sotelo-Ávila, 1974):

donde κ y γ son dos parámetros empíricos adimensionales. El parámetro κ es conocido como la constante de von Kármán.

Los valores experimentales reportados en la literatura se encuentran alrededor de los valores medios k ≅ 0.4 y γ ≅ 0.13; algunos autores señala que en particular la constante de von Kármán puede depender del número de Reynolds (e.g. Daily y Harleman, 1966; Landau y Lifchitz, 1989). Un análisis sobre el sujeto se hace ulteriormente.

La ley logarítmica es válida experimentalmente cuando y/λ_c > 30 a 70, y/δ > 0.15. En cuanto a la función G, en la ley de deflexión se ha argumentado experimentalmente que es una función logarítmica: (u_máx - u)/u_* = -κ^-1ln(y/δ) + B, donde B es una constante.

Considerando que la ley de deflexión puede ser incluida aproximadamente por la ecuación haciendo B ≈ 0, y que la contribución al caudal por la subcapa viscosa es despreciable, la velocidad media se obtiene, a pesar de que la velocidad calculada en la pared es u(0) = -∞ en lugar de cero, como el valor principal en el sentido de Cauchy de las integrales que definen el caudal en las ecuaciones (2) y (3). La forma general es la siguiente:

donde β = γe ≅ 0.35 en el canal; β = γe^3/2 ≅ 0.58 en un tubo de radio R completamente lleno, h se reemplaza por R.

Es claro que la ecuación (5) no es aplicable cuando h = 0 o R = 0, ya que proporciona U(0) = -∞ en lugar de cero. Este comportamiento de la velocidad media y del perfil de las velocidades se debe a que en la deducción de la ecuación (4) se desprecia la fricción debida a la viscosidad frente a la fricción por la turbulencia. No obstante, la ley de resistencia al flujo representa adecuadamente los datos experimentales para altos números de Reynolds.

El objetivo de este trabajo es la extensión del perfil logarítmico de las velocidades, a fin de obtener el comportamiento correcto en la pared y su vecindad, para establecer la ley de resistencia, que incluya la ley de Poiseuille cuando el número de Reynolds sea pequeño.

 

Longitud de mezcla hiperbólica

Para alcanzar el objetivo es necesario recordar que la ley de Newton, válida en régimen laminar, establece que el esfuerzo cortante (τ) y el gradiente vertical de las velocidades horizontales son proporcionales, es decir:

donde ρν = μ es el coeficiente de viscosidad dinámica.

El perfil parabólico se obtiene tomando en cuenta que el esfuerzo cortante varía linealmente, siendo nulo en la superficie libre del agua o en el centro del tubo, y adquiriendo su valor máximo en la pared (τ_p), es decir:

el perfil de velocidades resultante es el siguiente:

Las velocidades medias en el canal y el tubo se obtienen introduciendo la ecuación en las ecuaciones (2) y (3):

En régimen turbulento, el cortante, según las hipótesis de Reynolds, es definido por la media del producto de las fluctuaciones alrededor de las medias de las velocidades en el sentido horizontal (u') y en el sentido vertical (v'), es decir (Landau y Lifchitz, 1989):

Para la modelación del cortante turbulento citemos textualmente a Levi (1989): "En la teoría cinética de los gases se había utilizado un concepto importante: el de 'camino libre medio'. Análogamente, Prandtl introduce una longitud característica capaz de definir el estado de turbulencia, longitud que se podría interpretar como la distancia que una partícula fluida recorre sin perder su individualidad al mezclarse con sus vecinas, o bien como el diámetro de los agregados de partículas que avanzan juntas; distancia y diámetro que -cuando el número de Reynolds es bastante grande- tendrían que ser proporcionales entre sí, pues la energía cinética de la masa debe ser igual al trabajo realizado por las partículas para vencer la resistencia a su avance que las vecinas oponen. A esta longitud característica, Prandtl la llama 'longitud de mezcla' (Mischungsweg) y la representa con la letra , basándose en la hipótesis de que, en un movimiento dirigido en promedio según el eje x (o sea = 0), las partículas fluidas tienden a abandonar la capa a la cual pertenecen y -luego de cruzar, paralelamente al eje y, la distancia - alcanzan otra posición, donde la velocidad media es distinta de la que la partícula poseía originalmente. Ahora, por ser muy pequeña, la diferencia de velocidades entre la capa de salida y la de llegada es ∂/∂y; y esta expresión representaría -para Prandtl- la fluctuación de velocidad u' = ∂/∂y. Como otras partículas se mueven de manera análoga en dirección opuesta, habrá choque de masas fluidas con valores contrarios de u'; por tanto nacerán fluctuaciones v' (en el sentido transversal) del mismo orden de magnitud que u'. En conclusión, el esfuerzo aparente -ρ❬u'υ'❭ será proporcional a ρ²(∂/∂y)²''. Para establecer la igualdad, Prandtl (1925) propone que el coeficiente de proporcionalidad sea absorbido por la longitud de mezcla, y escribir:

en la cual se ha introducido el subíndice P a la longitud de mezcla de Prandtl; se ha eliminado la barra que denota el perfil de las velocidades medias; se han reemplazado las derivadas parciales por las totales, y se ha introducido la función valor absoluto para que el signo quede determinado por du/dy.

El perfil logarítmico de las velocidades, definido por la ecuación (4), se obtiene de la ecuación (12), bajo el supuesto de que el esfuerzo cortante es una constante e igual a su valor en la pared (τ = τ_p = ρu_*²) y que la longitud de mezcla en un punto dado es proporcional a la distancia a la pared: _P = κy. Como se ha visto, la ecuación (12) no es aplicable en la pared debido a que se han despreciado los efectos viscosos. Lo anterior puede ser corregido, considerando que el esfuerzo cortante es el resultado de los efectos viscosos y turbulentos, es decir:

La sustitución de las ecuaciones (6) y (12) en la ecuación (13) conduce a:

La ecuación (14) se escribe bajo la forma de la hipótesis de Newton, reemplazando el coeficiente de viscosidad cinemática en la ecuación (6) por un coeficiente de 'viscosidad cinemática total' (ν_t), a saber (Newton, 1985):

en la cual ν_t = ν + ν_turb y ν_turb = _P²du/dy, es decir la 'viscosidad cinemática turbulenta' es proporcional a u'_P.

De acuerdo con la ecuación (1), se define una longitud de mezcla total como λ = ν_t/u_*, de modo que λ = λ_c, cuando ν_t = ν. En consecuencia, las ecuaciones (6), (12) y (15) se escriben de manera adimensional, introduciendo las variables adimensionales y₊ = y/λ_c, u₊ = u/u_*, τ₊ = τ/τ_p, λ₊ = λ/λ_c = ν_t/ν, a saber:

donde τ_visc+ = τ_visc/τ_p, τ_turb+ = τ_turb/τ_p, _P+ = _P/λ_c.

Puesto que ν_t = νλ/λ_c = u_*λ, se tendrá u_*λ = ν + _P²du/dy, o en escritura adimensional λ₊ = 1 + ²_{P +}du₊/dy₊ , y considerando la ecuación (16) se tiene λ₊ = 1 + ²_P+τ₊/λ₊. De esta última se deduce:

donde ₊ = /λ_c, con l una longitud de mezcla efectiva, que coincide con la longitud de mezcla de Prandtl cuando τ = τ_p.

Las relaciones entre los diferentes esfuerzos cortantes se deducen de las ecuaciones (16) y (17), a saber:

Para τ₊ = 1 y cuando y >>λ_c según Prandtl se tiene ~ κy, y de la ecuación (17) se infiere también λ ~ κy (λ₊ = λ/λ_c >> 1); en la pared (y = 0) se tiene λ = λ_c (λ₊ = 1). Se propone unir estos comportamientos a través de la hipérbola λ₊² = 1 + by₊ + (κy₊)². Para que existan las inestabilidades que originan la turbulencia se acepta que el perfil de velocidades presente un punto de inflexión a una distancia y₊ > 0 de la pared (Landau y Lifchitz, 1989).

De la ecuación (16) se deduce d²u₊/dy₊² = -(1/λ²₊ )dλ₊/dy₊ para τ₊ = 1, y en consecuencia dλ₊/dy₊ = 0 en un punto cercano a la pared. La hipérbola en cuestión presenta un mínimo en y₊ = -b/2κ², implicando b < 0. Haciendo b = -2κα con α > 0, la ecuación de la hipérbola se escribe como sigue:

El punto mínimo de la hipérbola queda definido por y_m+ = α/κ, λ_m+ = , y se tiene que λ₊ = 1 tanto en la pared como a una distancia y_ƒ+ = 2y_m+ de la misma, por consiguiente λ_m+ ≤ λ₊ ≤ 1 en el intervalo 0 ≤ y₊ ≤ y_ƒ+ .

De acuerdo con la ecuación (17), en el intervalo 0 ≤ y₊ ≤ y_ƒ+ , se tiene ²₊ = -λ₊ (1 -λ₊) ≤ 0, es decir que ₊ es imaginaria, ya que λ_m+ ≤λ₊ ≤ 1. Esto se traduce en que el esfuerzo cortante turbulento definido en la ecuación (18) cambia de signo en esta región vecina a la pared. El hecho de que ₊ es imaginaria puede ser interpretado considerando que partículas de fluido en el plano xy, donde x es la dirección principal del flujo, son transferidas a otro plano x'y paralelo al anterior; otras partículas llegan al plano xy procedentes de planos paralelos, de modo que la masa se conserva en cada plano. La igualdad τ_visc+ = τ_turb+ se da cuando λ₊ = 2 y corresponde a y_νt+ = y_m+.

 

Un valor teórico de la constante de von Kármán

Los valores de los parámetros α y κ pueden ser estimados bajo algunas consideraciones. Derivando la ecuación (17) se obtiene:

Ahora bien, el punto mínimo de la longitud de mezcla hiperbólica dλ₊/dy₊|_ym+ = 0 implica que d₊/dy₊|_ym+ = 0. De acuerdo con la ecuación (20), este resultado se obtiene también si λ_m+ = ½, lo que implica α = . De la ecuaciones (16) y (17) se deduce que en el punto crítico se tiene ²_m+ = -¼ o sea _m+ = ±½i, τ_visc+ = 2τ₊ y τ_turb+ = -τ₊.

En la región 0 < y < y_ƒ, las escalas importantes son precisamente las que forman el punto crítico (y_m,λ_m); la primera caracteriza una masa crítica y la segunda es una longitud de mezcla crítica que caracteriza las inestabilidades que originan la turbulencia. Haciendo intervenir la viscosidad cinemática, el análisis dimensional proporciona las velocidades correspondientes u_m = ν/y_m y v_m = ν/λ_m, así como el cortante en la pared τ_p = ρu_mv_m; este cortante tiene la forma de aquel definido por la ecuación (11). Dado que τ_p = ρu²_*, ecuación (1), se tiene u²_* = u_mv_m, y λ²_c = y_mλ_m, es decir, en escritura adimensional u_m+v_m+ = 1 y y_m+λ_m+ = 1. De modo que λ_m+ = ½ implica que y_m+ = α/κ = 2, esto es κ = y, ver ecuación (23), γ = , valores muy cercanos a los reportados en la ecuación (4). Por la importancia de estos parámetros se establecen a continuación los valores resultantes del análisis precedente:

Es importante señalar que este valor teórico de la constante de von Kármán es muy cercano al valor obtenido por Zaragola (1996), a partir de ensayos realizados en el rango 3.1 × 10⁴ a 3.5 × 10⁷ del número de Reynolds más amplio que el rango estudiado por Nikuradse (1932, 1933) (Zaragola y Smits, 1997; McKeon et al., 2004), a saber: κ ≅ 0.436.

Con los valores de α y κ definidos por la ecuación (21) se deduce la frontera superior de la región crítica cercana a la pared y_ƒ+ = 4 y la distancia a la pared donde los esfuerzos viscoso y turbulento son iguales y_νy+= . Es de notar que el valor de y_ƒ+ = 4 coincide con el valor experimental del espesor de la subcapa viscosa.

Las gráficas de la longitud de mezcla hiperbólica definida por la ecuación (19) y la asíntota κy₊ - α, con los valores de los parámetros proporcionados por la ecuación (21), se muestran en la figura 1.

 

Perfil de las velocidades

El perfil de velocidades adimensional correspondiente a τ₊ = 1 se obtiene de las ecuaciones (16) y (19), con la condición u₊(0) = 0, a saber:

en donde se ha introducido la notación u_∞+(y₊) para indicar que τ₊ = 1 se obtiene de la ecuación (8) cuando h₊ → ∞.

Cuando y₊ >>, 0 la ecuación (22) tiene el comportamiento de la ecuación (4):

y cuando y₊ → 0, su desarrollo es el siguiente:

el cual muestra el comportamiento lineal en la subcapa viscosa cerca de la pared (Panton, 1984).

El perfil adimensional de las velocidades correspondiente a la ecuación (7), τ₊ = 1 - y₊/h₊, se obtiene de las ecuaciones (16) y (19), con la condición u₊(0) = 0, a saber:

Es importante señalar que cuando y₊ es muy grande y en consecuencia también h₊, ya que 0 < y₊ < h₊, de la ecuación (25) se deduce el comportamiento u₊(y₊,h₊) = κ^-1(1 - α/κh₊) ln(y₊/γ). Escribiendo la ecuación (4) como u₊ = κ ^-1_vk ln(y₊/γ), en donde se ha cambiado la notación para el parámetro original de von Kármán (κ_vk); es claro que κ ^-1_vk = κ ^{- 1}(1 - α/κh₊), la cual puede justificar la dependencia de dicho parámetro con respecto al número de Reynolds definido por Re^* = h₊ = u_*h/ν, argumentada por algunos autores cuando el perfil logarítmico es aplicado en flujos donde el esfuerzo cortante no puede ser considerado constante. Claramente, cuando h₊ → ∞, la ecuación (25) se reduce a la ecuación (22).

El comportamiento de u₊ (y₊,h₊) cuando y₊ → 0 es el siguiente:

en donde de nuevo aparece el comportamiento lineal en la subcapa viscosa.

El perfil parabólico correspondiente al régimen laminar y proporcionado por la ecuación (8) en escritura adimensional es u ₊ = y₊ - y²₊ /2h₊, el cual se deduce de la ecuación (26) cuando h ₊ → 0.

La velocidad máxima ocurre a la distancia y₊ = h₊ de la pared; la diferencia entre esta velocidad y aquella en y₊ cercana a h₊ se calcula con las ecuaciones (22) y (25). Los resultados son las leyes de deflexión siguientes:

donde la primera se deduce de la segunda haciendo h₊ → ∞ y h₊ - y₊ → 0.

Los perfiles de las velocidades definidos por la ecuación (22) y (25) se muestran en la figura 2, en los cuales se debe observar la existencia del punto de inflexión.

 

Velocidades medias

Se obtienen a continuación las expresiones de la velocidad media del agua en un tubo completamente lleno y en un canal rectangular muy ancho.

 

Tubo. La velocidad media adimensional para τ₊ = 1 se obtiene introduciendo la ecuación (22) en la ecuación (3):

Los comportamientos de U_T∞+(R₊), cuando R₊ → 0 y R₊ → ∞ son, respectivamente, los siguientes:

la ecuación (31) corresponde a la ecuación (5).

La velocidad media adimensional para τ_* = 1 - y_*/R_* se deduce de las ecuaciones (3) y (25):

Los comportamientos de U_T+(R₊) cuando R₊ → 0 y R₊ → ∞ son, respectivamente, los siguientes:

En escritura adimensional, la ley de Poiseuille de los tubos capilares definida por la ecuación (10) es U_T+(R₊) = ¼R₊, la cual se deduce precisamente de la ecuación (33).

 

Canal rectangular. La velocidad media adimensional para τ₊ = 1 se obtiene introduciendo la ecuación (22) en la ecuación (2):

Los comportamientos de U_C∞+(h₊) cuando h₊ → 0 y h₊ → ∞ son, respectivamente, los siguientes:

la ecuación (37) corresponde a la ecuación (5).

La velocidad media adimensional para τ₊ = 1 - y₊/h₊ se obtiene introduciendo la ecuación (25) en la ecuación (2):

Los comportamientos de U_C+(h₊) cuando h₊ → 0 y h₊ → ∞ son, respectivamente, los siguientes:

En escritura dimensional, la ecuación (39) corresponde al régimen de Poiseuille sobre una placa plana, ecuación (9).

 

El factor de fricción

La teoría anterior es válida en paredes lisas. En paredes rugosas y a grandes números de Reynolds, la escala λ_c se reemplaza por la escala λ_ε, que caracteriza la rugosidad de la pared; siguiendo a Colebrook y White (Colebrook, 1939), en régimen intermedio, la escala λ_c se reemplaza por la escala λ_t, definida como λ_t = λ_c + λ_ε.

A altos números de Reynolds, la escala λ_ε es proporcional a la altura media de las rugosidades de la pared, tomada como la altura equivalente de granos de arena (ε) o escala de Nikuradse, es decir:

donde c es un coeficiente de proporcionalidad.

En un tubo, la ley de resistencia se aplica utilizando el coeficiente de fricción de Fanning (ƒ) definido en la ecuación de Darcy-Weisbach J = ƒU²/2gD, donde D es el diámetro interno del tubo (Darcy, 1857; Darcy y Bazin, 1865; Weisbach, 1850). Considerando la velocidad de fricción, este coeficiente queda definido como ƒ = 8/U₊ ². Reemplazando la escala λ_c por la escala λ_t = λ_c + λ_ε, la relación entre el coeficiente de fricción, el número de Reynolds definido por Re = UD/ν, y la escala relativa definida por λ_εr = λ_ε/D, se obtiene de la ecuación (29) o (32):

donde λ_εr = cε_r cuando Re → ∞ y ε_r = ε/D es la rugosidad relativa de Nikuradse.

El parámetro c se estima a partir de la información obtenida por Nikuradse en sus experiencias realizadas en tubos con seis valores de rugosidad relativa a grandes números de Reynolds (Sotelo-Ávila, 1974). Los valores de (ε/D,ƒ) se presentan en el cuadro 1; asimismo, se reportan los valores correspondientes de λ_ε/D obtenidos utilizando en la ecuación (42) las ecuaciones (29) y (32) con Re → ∞. La relación entre λ_ε/D y ε/D se muestra en la figura 3. Para la ecuación (29), c ≅ 0.125 con R² = 0.9996, mientras que para la ecuación (32), c ≅ 0.09 con R² = 0.9997.

En la zona de transición, la función λ_ε(Re, ε) es desconocida y podrá ser aproximada empíricamente a partir de los datos experimentales. Con la finalidad de que la función ƒ[Re, λ_ε(Re, ε)], que define al factor de fricción implícitamente a través de la ecuación (42), manifieste algunas de sus características, se tomará el comportamiento asintótico definido por la ecuación (41).

La función ƒ(Re, cε), mostrada en la figura 4 para los seis valores de la rugosidad relativa de Nikuradse, utiliza la ecuación (32). La función así construida no proporciona el comportamiento observado del factor de fricción en la zona de transición; su comportamiento asintótico es similar al del factor de fricción mostrado en el clásico diagrama de Moody (1944, 1947). Sin embargo, el diagrama de la figura 4 tiene la principal ventaja de que contiene como caso límite, a bajos números de Reynolds, al régimen laminar de Poiseuille para el cual ƒ = 64/Re.

En el caso de tubos o canales rugosos, si la rugosidad de la pared se obtiene de manera inversa en un problema de flujo en régimen laminar o de transición, es posible que el valor no tenga ningún significado físico, por ejemplo, que sea del mismo orden que el del diámetro del tubo, debido a que ƒ(Re, cε) es aproximada en esta región de flujo. En esta situación, tal vez es mejor retener el valor obtenido de la escala λ_ε(Re, ε), la cual estará asociada con un número de Reynolds del régimen de flujo y con una rugosidad de la pared, y ser utilizada como un valor medio en un pequeño rango de variación de este número.

 

Conclusiones

Se ha considerado que el esfuerzo cortante en el flujo del agua es el resultado de la suma de un esfuerzo viscoso y de un esfuerzo turbulento. El primero es proporcional, de acuerdo con Newton, a la derivada del perfil de las velocidades, mientras que el segundo, de acuerdo con Prandtl, al cuadrado del producto de la derivada de este perfil y una longitud de mezcla. Se ha introducido una longitud de mezcla hiperbólica para escribir la relación entre el esfuerzo cortante total y la derivada del perfil en la forma de Newton.

La longitud hiperbólica ha sido construida de tal manera que tenga el comportamiento lineal con respecto a la distancia a la pared de la de Prandtl, en la cual el coeficiente de proporcionalidad es la constante de von Kármán, y que tenga un valor constante en la pared para incluir el régimen laminar. Esto induce que la hipérbola presente un punto mínimo cerca de la pared, hecho que a su vez induce un punto de inflexión en el perfil de las velocidades. La existencia de este punto crítico permite deducir un valor teórico de la constante de von Kármán bastante cercano a su valor experimental.

Con el perfil de las velocidades obtenido se han deducido las fórmulas para calcular las velocidades medias del fluido en un tubo completamente lleno de agua y en un canal rectangular de ancho muy grande; a partir de estas fórmulas se obtienen las leyes de resistencia al flujo. Se ha estudiado en particular la ley de resistencia en el tubo llevándola a la forma de la ecuación de Darcy-Weisbach; el factor de fricción resultante presenta los comportamientos extremos clásicos, en función del número de Reynolds, es decir el régimen laminar de Poiseuille a bajos valores y el de Prandtl a altos valores de este número.

El efecto de la rugosidad de la pared a altos números de Reynolds se ha introducido en la ley de resistencia, siguiendo la idea de Colebrook y White de sumar a la escala de la turbulencia una escala proporcional al tamaño de los granos de arena de Nikuradse. El factor de fricción resultante tiene el mismo comportamiento que el mostrado en el diagrama clásico de Moody. Sin embargo, esta escala depende del número de Reynolds en régimen transitorio de acuerdo con los resultados experimentales de Nikuradse. La dependencia de esta escala a pequeños y moderados valores del número de Reynolds y de la rugosidad es un sujeto de investigación.

 

Agradecimientos

El autor manifiesta su profundo agradecimiento a don Samuel Trueba Coronel, ingeniero agrónomo especialista en irrigación de la Escuela Nacional de Agricultura de Chapingo y doctor Honoris causa por el Colegio de Postgraduados, por sus enseñanzas en hidráulica.

 

Referencias

BARENBLATT, G.I. Scaling laws for fully developed turbulent shear flows. Part 1. Basic hypotheses and analysis. Journal of Fluid Mechanics. Vol. 248, 1993, pp. 513-520.         [ Links ]

BARENBLATT, G.I. and PROSTOKISHIN, V.M. Scaling laws for fully developed turbulent shear flows. Part 2. Processing of experimental data. Journal of Fluid Mechanics. Vol. 248, 1993, pp. 521-529.         [ Links ]

CHÉZY, A. Rapport sur le canal de l'Yvette. En Herschel, C. On the origin of the Chézy formula. J. Assoc. Eng. Soc. 1897. Vol. 18, 1776, pp. 363-365.         [ Links ]

COLEBROOK, C.F. Turbulent flow in pipes with particular reference to the transition region between the smooth and rough pipe laws. J. Inst. Civ. Eng. Vol. 11, 1939, pp. 133-156.         [ Links ]

DAILY, J.W. and HARLEMAN, D.R.F Fluid dynamics. Reading, USA: Addison-Wesley Publishing Co., 1966, 511 pp.         [ Links ]

DARCY, H.P.G. Recherches expérimentales relatives au mouvement de l'eau dans les tuyaux. Paris: Imprimerie Nationale, 1857.         [ Links ]

DARCY, H. et BAZIN, H. Recherches hydrauliques; première partie, recherches expérimentales sur l'écoulement de l'eau dans les canaux découverts. Paris: Imprimerie Nationale, 1865.         [ Links ]

DU BUAT, P.L.G. Principes d'hydraulique, vérifiés par un grand nombre d'expériences faites par ordre du gouvernement. Secund. ed. Paris: Imprimerie de Monsieur, 1786.         [ Links ]

FUENTES, C., DE LEÓN, B., SAUCEDO, H., PARLANGE, J.-Y. y ANTONINO, A.C.D. El sistema de ecuaciones de Saint-Venant y Richards del riego por gravedad: 1. La ley de resistencia hidráulica en potencia. Ingeniería hidráulica en México. Vol. XIX, núm. 2, abril-junio de 2004, pp. 65-75.         [ Links ]

KING, H.W., WISLER, C.O. and WOODBURN, J.G. Hydraulics. New York: John Wiley & Sons, Inc., 1952, 351 pp.         [ Links ]

LANDAU, L. et LIFCHITZ, E. Physique théorique. Tome 6: Mécanique des fluides. Seconde édition. Moscou: Éditions Mir, 1989, 748 pp.         [ Links ]

LEVI, E. El agua según la ciencia. Conacyt-Castell Mexicana, México, 1989, 677 pp.         [ Links ]

MANNING, R. On the flow of water in open channels and pipes. Trans. Inst. Civil Engr. Ireland. Vol. 20, 1891, pp. 161-207.         [ Links ]

McKEON, B.J., SWANSON, C.J., ZARAGOLA, M.V., DONNELLY, R.J. and SMITS, A.J. Friction factors for smooth pipe flow. J. Fluid Mechanics. Vol. 541, 2004, pp. 41-44.         [ Links ]

MOODY, L.F., Friction factors for pipe flows. Trans. ASME. Vol. 66, 1944, pp. 641.         [ Links ]

MOODY, L.F. An approximate formula for pipe friction factors. Trans. ASME. Vol. 69, 1947, pp. 1005.         [ Links ]

NEWTON, I. De philosophiae naturalis principia mathematica. Préface de S. Hawking. Postface de F. Biarnais. Paris: Christian Bourgois Éditeur, 1985.         [ Links ]

NIKURADSE, J. Gesetzmäßigkeit der turbulenten Stromüng in glatten Rohren. Forsch. Arb. Ing.-Wes. No. 356. English translation NACA TT F-10, 1932, 359 pp.         [ Links ]

NIKURADSE, J. Stroemungsgesetze in rauhen Rohren. Ver. Dtsch. Ing. Forsch. 1933, 361 pp.         [ Links ]

PANTON, R.L. Incompressible flow. New York: John Wiley & Sons, Inc., 1984, 780 pp.         [ Links ]

POISEUILLE, J.L.M. Recherches expérimentales sur le mouvement des liquides dans les tubes de très petits diamètres. Comptes rendus de l'Académie des Sciences. Vol. 11, 1840, pp. 961-967 y 1041-1049; vol. 12, 1841, pp. 112-115.         [ Links ]

PRANDTL, L. Über die ausgebildete Turbulenz. ZAMM (Zeitschrift für Angewandte Mathematik und Mechanik). Vol. 5, 1925, pp. 136-139.         [ Links ]

REYNOLDS, O. An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous and of the Law of Resistance in parallel channels. Philos Trans. Royal Soc. London. Vol. 174, 1883, pp. 935-982.         [ Links ]

RICHARDS, L.A. Capillary conduction of liquids through porous mediums. Physics. Vol. 1, 1931, pp. 318-333.         [ Links ]

SAINT-VENANT, A.J.C. BARRÉ DE. Théorie du mouvement non permanent des eaux avec application aux crues des rivières, et à l'introduction des marées dans leurs lits. Comptes rendus des Séances de l'Académie des Sciences. Vol. 73, 1871 pp. 147-154 y 237-240.         [ Links ]

SOTELO-ÁVILA, G. Hidráulica General. México, D.F.: Editorial Trillas, 1974.         [ Links ]

TRUEBA-CORONEL, S. Hidráulica. México, D.F.: CECSA, 1954, 454 pp.         [ Links ]

WEISBACH, J. Lehrbuch der Ingenieur-und Maschinen-Mechanik. 2nd ed. Braunschweig: Vieweg und Sohn, 1850, 535 pp.         [ Links ]

WHITE, F.M. Mecánica de fluidos. México, D.F.: McGraw Hill, 1979, 757 pp.         [ Links ]

ZARAGOLA, M.V. Mean-flow Scaling of Turbulent Pipe Flow. Ph.D. thesis. Princeton: Princeton University, 1996.         [ Links ]

ZARAGOLA, M.V. and SMITS, A.J. Scaling of the mean velocity profile for turbulent pipe flow. Physics Review Letters. Vol. 78, núm. 2, 1997, pp. 239-242.         [ Links ]