Ley de la pared

Ley de la pared, velocidad horizontal cerca de la pared con modelo de longitud de mezcla

En dinámica de fluidos, la ley de la pared (también conocida como ley logarítmica de la pared) establece que la velocidad media de un flujo turbulento en un punto determinado es proporcional al logaritmo de la distancia desde ese punto a la "pared", o límite de la región de fluido. Esta ley de la pared fue publicada por primera vez en 1930 por el matemático húngaro-estadounidense, ingeniero aeroespacial y físico Theodore von Kármán.^[1]​ Sólo es técnicamente aplicable a las partes del flujo que están cerca de la pared (<20% de la altura del flujo), aunque es una buena aproximación para todo el perfil de velocidad de las corrientes naturales.^[2]​

Formulación logarítmica general

La ley logarítmica de la pared es una solución autosimilar para la velocidad media paralela a la pared, y es válida para flujos a número de Reynoldss altos -en una región de solapamiento con esfuerzo cortante aproximadamente constante y lo suficientemente lejos de la pared para que los efectos (directos) viscosos sean despreciables:^[3]​

${\displaystyle u^{+}={\frac {1}{\kappa }}\ln \,y^{+}+C^{+},}$   with   ${\displaystyle y^{+}={\frac {y\,u_{\tau }}{\nu }},}$   ${\displaystyle u_{\tau }={\sqrt {\frac {\tau _{w}}{\rho }}}}$   and   ${\displaystyle u^{+}={\frac {u}{u_{\tau }}}}$

donde

${\displaystyle y^{+}}$	es la coordenada de la pared: la distancia y a la pared, hecha adimensional con la velocidad de fricción uτ y la viscosidad cinética ν,
${\displaystyle u^{+}}$	es la velocidad adimensional: la velocidad u paralela a la pared en función de y (distancia a la pared), dividida por la velocidad de fricción uτ,
${\displaystyle \tau _{w}}$	es el esfuerzo cortante de la pared,
${\displaystyle \rho }$	es la densidad del fluido,
${\displaystyle u_{\tau }}$	se llama velocidad de fricción o velocidad de corte,
${\displaystyle \kappa }$	es la constante de von Kármán,
${\displaystyle C^{+}}$	es una constante y
${\displaystyle \ln }$	es el logaritmo natural.

A partir de experimentos, la constante de von Kármán resulta ser ${\displaystyle \kappa \approx 0.41}$ y ${\displaystyle C^{+}\approx 5.0}$ para una pared lisa.^[3]​

Con dimensiones, la ley logarítmica de la pared puede escribirse como:^[4]​

${\displaystyle {u}={\frac {u_{\tau }}{\kappa }}\ln \,{\frac {y}{y_{0}}}\ }$

donde y₀ es la distancia desde la frontera a la que la velocidad idealizada dada por la ley de la pared llega a cero. Esto es necesariamente distinto de cero porque el perfil de velocidad turbulenta definido por la ley de la pared no se aplica a la subcapa laminar. La distancia de la pared a la que llega a cero se determina comparando el espesor de la subcapa laminar con la rugosidad de la superficie sobre la que fluye. Para una subcapa laminar próxima a la pared de espesor ${\displaystyle \delta _{\nu }}$ y una rugosidad característica de escala de longitud ${\displaystyle k_{s}}$,^[2]​

${\displaystyle k_{s}<\delta _{\nu }\,}$	: flujo hidráulicamente suave ,
${\displaystyle k_{s}\approx \delta _{\nu }\,}$	: flujo de transición,
${\displaystyle k_{s}>\delta _{\nu }\,}$	: flujo hidráulicamente áspero .

Intuitivamente, esto significa que si los elementos de rugosidad están ocultos dentro de la subcapa laminar, tienen un efecto muy diferente en la ley turbulenta del perfil de velocidad de la pared que si sobresalen en la parte principal del flujo.

Esto también se formula a menudo más formalmente en términos de un número de Reynolds límite, ${\displaystyle Re_{w}}$, donde

${\displaystyle Re_{w}={\frac {u_{\tau }k_{s}}{\nu }}.\ }$

El flujo es hidráulicamente suave para ${\displaystyle Re_{w}<3}$, hidráulicamente rugoso para ${\displaystyle Re_{w}>100}$, y transicional para valores intermedios.^[2]​

Los valores de ${\displaystyle y_{0}}$ vienen dados por:^[2]​^[5]​

${\displaystyle y_{0}={\frac {\nu }{9u_{\tau }}}}$	para un flujo hidráulicamente suave
${\displaystyle y_{0}={\frac {k_{s}}{30}}}$	para flujos hidráulicamente rugosos.

Los valores intermedios suelen venir dados por el diagrama de Nikuradse derivado empíricamente,^[2]​ aunque también se han propuesto métodos analíticos para resolver este rango.^[6]​

Para canales con una frontera granular, como los sistemas fluviales naturales,

${\displaystyle k_{s}\approx 3.5D_{84},\ }$

donde ${\displaystyle D_{84}}$ es el diámetro medio del percentil 84 más grande de los granos del material del lecho.^[7]​

Soluciones de ley de potencia

Los trabajos de Barenblatt y otros han demostrado que además de la ley logarítmica de la pared - el límite para infinitos números de Reynolds - existen soluciones de ley de potencia, que dependen del número de Reynolds.^[8]​^[9]​ En 1996, Cipra presentó pruebas experimentales en apoyo de estas descripciones de ley de potencia.^[10]​ Esta evidencia en sí no ha sido totalmente aceptada por otros expertos.^[11]​ En 2001, Oberlack afirmó haber derivado tanto la ley logarítmica de la pared, como las leyes de potencia, directamente de las ecuaciones de Navier-Stokes promediadas por Reynolds, explotando las simetrías en un enfoque de grupo de Lie.^[3]​^[12]​ Sin embargo, en 2014, Frewer et al.^[13]​ refutaron estos resultados.

Para escalares

Para escalares (más notablemente la temperatura), la ley logarítmica autosimilar de la pared ha sido teorizada (formulada por primera vez por B. A. Kader^[14]​) y observado en estudios experimentales y computacionales.^[15]​^[16]​^[17]​^[18]​ En muchos casos, las extensiones a la ley original de la formulación de la pared (generalmente a través de transformaciones integrales) son generalmente necesarias para tener en cuenta la compresibilidad, la propiedad variable y los efectos reales del fluido.

Cerca de la pared

Por debajo de la región donde es aplicable la ley de la pared, existen otras estimaciones para la velocidad de rozamiento.^[19]​

Subcapa viscosa

En la región conocida como subcapa viscosa, por debajo de 5 unidades de pared, la variación de ${\displaystyle u^{+}}$ a ${\displaystyle y^{+}}$ es aproximadamente 1:1, tal que:

Para  ${\displaystyle y^{+}<5}$

${\displaystyle u^{+}=y^{+}}$

donde

${\displaystyle y^{+}}$	es la coordenada de la pared: la distancia y a la pared, hecha adimensional con la velocidad de fricción ${\displaystyle u_{\tau }}$ y viscosidad cinemática ${\displaystyle \nu }$,
${\displaystyle u^{+}}$	es la velocidad adimensional: la velocidad u paralela a la pared en función de y (distancia a la pared), dividida por la velocidad de fricción ${\displaystyle u_{\tau }}$,

Esta aproximación puede usarse a más de 5 unidades de pared, pero para ${\displaystyle y^{+}=12}$ el error es superior al 25%.

Capa intermedia

En la capa intermedia, entre 5 unidades de pared y 30 unidades de pared, no se cumple ninguna de las dos leyes, de modo que:

Para  ${\displaystyle 5<y^{+}<30}$

${\displaystyle u^{+}\neq y^{+}}$

${\displaystyle u^{+}\neq {\frac {1}{\kappa }}\ln \,y^{+}+C^{+}}$

con la mayor variación de cualquiera de las leyes ocurriendo aproximadamente donde las dos ecuaciones se interceptan, en ${\displaystyle y^{+}=11}$. Es decir, antes de 11 unidades de pared la aproximación lineal es más precisa y después de 11 unidades de pared se debe utilizar la aproximación logarítmica, aunque ninguna de las dos es relativamente precisa a 11 unidades de pared.

El perfil de velocidad media en sentido de la corriente ${\displaystyle u^{+}}$ se mejora para ${\displaystyle y^{+}<20}$ con una formulación de viscosidad de remolino basada en una energía cinética turbulenta cercana a la pared. ${\displaystyle \kappa ^{+}}$ y la ecuación de longitud de mezcla de van Driest. Las comparaciones con datos de simulación numérica directa de flujos de canales turbulentos completamente desarrollados para ${\displaystyle 109<Re_{\tau }<2003}$ mostraron una buena concordancia.^[20]​

Referencias

1. von Kármán, Th. (1930), «Mechanische Ähnlichkeit und Turbulenz», Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Fachgruppe 1 (Mathematik) 5: 58-76 . (también como: "Similitud mecánica y turbulencia", Tech. Mem. NACA, nº 611, 1931).
2. Mohrig, David (2004). edu/courses/earth-atmospheric-and-planetary-sciences/12-110-sedimentary-geology-fall-2004/lecture-notes/l2_mass_mom_cons.pdf «Conservación de la masa y el momento». 12.110: Geología sedimentaria, otoño de 2004. MIT OCW. Consultado el 27 de marzo de 2009. 
3. Schlichting & Gersten (2000) pp. 522-524.
4. Schlichting & Gersten (2000) p. 530.
5. Whipple, Kelin (2004). «Rugosidad hidráulica». 12.163: Procesos superficiales y evolución del paisaje. MIT OCW. Consultado el 27 de marzo de 2009. 
6. Le Roux, J.P. (2004), «Una ley integrada de la pared para flujos hidrodinámicamente transitorios sobre lechos planos», Sedimentary Geology 163 (3-4): 311-321, Bibcode:2004SedG..163..311L .
7. Haws, Benjamin. edu/~fange/Research/Articles/Thesis/Equivalent%20sand%20roughness%20of%20Nikuradse.doc «Rugosidad equivalente de la arena de Nikuradse (ks)». Consultado el 27 de marzo de 2009. 
8. Lynn Yarris. «Un fallo en la ley». Berkeley Lab: Highlights 97-98. Lawrence Berkeley National Laboratory, U.S. Department of Energy. 
9. Barenblatt, G.I. (1993), «Leyes de escala para flujos de cizalla turbulentos completamente desarrollados. Parte 1. Hipótesis básicas y análisis. Parte 1. Hipótesis básicas y análisis», Journal of Fluid Mechanics 248: 513-520, Bibcode:1993JFM...248..513B, doi:10.1017/S0022112093000874 .
   Barenblatt, G.I.; Prostokishin, V.M. (1993), «Leyes de escala para flujos turbulentos de cizalla completamente desarrollados. Part 2. Procesamiento de datos experimentales», Journal of Fluid Mechanics 248: 521-529, Bibcode:1993JFM...248..521B, doi:10.1017/S0022112093000886 .
   Barenblatt, G.I.; Goldenfeld, N. (1995), «¿Existe la turbulencia completamente desarrollada? Independencia del número de Reynolds frente a covarianza asintótica», Physics of Fluids 7 (12): 3078-3084, Bibcode:1995PhFl....7.3078B, arXiv:cond-mat/9507132, doi:10.1063/1.868685 .
   Barenblatt, G. I.; Chorin, A. J. (1998), «Leyes de escala y límites de viscosidad decreciente para flujos de cizalla limitados por paredes y para estructura local en turbulencia desarrollada», Communications on Pure and Applied Mathematics 50 (4): 381-398, doi:10.1002/(SICI)1097-0312(199704)50:4<381::AID-CPA5>3.0.CO;2-6 .
10. Cipra, Barry Arthur (May 1996), «Una nueva teoría de la turbulencia causa revuelo entre los expertos», Science 272 (5264): 951, Bibcode:1996Sci...272..951C .
11. Zagarola, M.V.; Perry, A.E.; Smits, A.J. (1997), «Leyes Log o Leyes de Potencia: La escala en la región de solapamiento», Physics of Fluids 9 (7): 2094-2100, Bibcode:1997PhFl....9.2094Z, doi:10.1063/1.869328 .
12. Oberlack, Martin (2001), «Un enfoque unificado para simetrías en flujos de cizalla turbulentos planos paralelos», Journal of Fluid Mechanics 427: 299-328, Bibcode:2001JFM...427..299O .
13. Frewer, Michael; Khujadze, George; Foysi, Holger (2014), ¿Es la ley logarítmica un resultado de primer principio del análisis de invariancia de grupos de Lie?, pp. 1-32, Bibcode:2014arXiv1412.3069F, arXiv:1412.3069 .
14. Kader, B. A. (1 de septiembre de 1981). «Perfiles de temperatura y concentración en capas límite totalmente turbulentas». International Journal of Heat and Mass Transfer 24 (9): 1541-1544. ISSN 0017-9310. 
15. Simonich, J. C.; Bradshaw, P. (1 de noviembre de 1978). 3450875 «Efecto de la turbulencia de chorro libre en la transferencia de calor a través de una capa límite turbulenta». Journal of Heat Transfer 100 (4): 671-677. ISSN 0022-1481. 
16. Patel, Ashish; Boersma, Bendiks J.; Pecnik, Rene (21 de agosto de 2017). «Estadística escalar en flujos de canal turbulento de propiedad variable». Physical Review Fluids 2 (8): 084604. 
17. Toki, Takahiko; Teramoto, Susumu; Okamoto, Koji (1 de enero de 2020). aiaa.org/doi/10.2514/1.B37381 «Perfiles de velocidad y temperatura en flujo de canal turbulento a presión supercrítica». Journal of Propulsion and Power 36 (1): 3-13. 
18. Guo, J.; Yang, X. I. A.; Ihme, M. (Marzo 2022). «Estructura de la capa límite térmica en flujos de canales turbulentos en condiciones transcríticas». Journal of Fluid Mechanics 934. ISSN 0022-1120. 
19. Flujos turbulentos (2000) pp. 273-274.Pope, Stephen (2000), Flujos turbulentos (1ª revisada edición), Cambridge University Press, ISBN 0-521-59125-2 .
20. Absi, Rafik (2009), «A simple eddy viscosity formulation for turbulent boundary layers near smooth walls», Comptes Rendus Mécanique 337 (3): 158-165, Bibcode:2009CRMec.337..158A, arXiv:1106.0985, doi:10.1016/j.crme.2009.03.010 .

Bibliografía

• Chanson, H. (2009), Applied Hydrodynamics: An Introduction to Ideal and Real Fluid Flows, CRC Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages, ISBN 978-0-415-49271-3 .
• Schlichting, Hermann; Gersten, K. (2000), Boundary-layer Theory (8th revised edición), Springer, ISBN 3-540-66270-7 .

Enlaces externos

• Definition from ScienceWorld
• Formula on CFD Online
• Y+ estimator