MECÁNICA DE FLUIDOS

CAPA LÍMITE LAMINAR Y TURBULENTA

Docente:

Alumnos:
RESUMEN:
Primeramente, se plantearán las ecuaciones de la capa limite laminar, limitándonos al
flujo bidimensional e incomprensible. De dichas ecuaciones extraeremos, mediante
estimación de órdenes de magnitud, una serie de conclusiones prácticas. Se resolverá
exactamente el caso del movimiento alrededor de una placa plana en el caso laminar.
Seguidamente, se planteará la integral de Karman, que da un método aproximado de
resolución para el caso turbulento, y también para capas límites alrededor de paredes de
formas arbitrarias. Esta ecuación se utilizará para obtener un resultado aproximado de la
capa límite turbulento alrededor de una placa plana. Seguidamente, se estudiará el
fenómeno de desprendimiento de la capa límite. Por último, de forma breve, y
basándonos fundamentalmente en la estimación de órdenes de magnitud, se estudiarán
las capas limites térmicas, asociadas a la transmisión de calor. (Crespo M., 2006)
 :
En mecánica de fluidos, la capa límite de un fluido se entiende como aquella en la que la
velocidad del fluido respecto al sólido en movimiento varía desde cero hasta el 99% de la
velocidad de la corriente no perturbada.

La capa límite puede ser laminar o turbulenta; aunque también pueden coexistir en ella zonas de
flujo laminar y de flujo turbulento. En ocasiones es de utilidad que la capa límite sea turbulenta.
En aeronáutica aplicada a la aviación comercial, se suele optar por perfiles alares que generan
una capa límite turbulenta, ya que ésta permanece adherida al perfil a mayores ángulos de
ataque que la capa límite laminar, evitando así que el perfil entre en pérdida, es decir, deje de
generar sustentación aerodinámica de manera brusca por el desprendimiento de la capa límite.

El espesor de la capa límite en la zona del borde de ataque o de llegada es pequeño, pero
aumenta a lo largo de la superficie. Todas estas características varían en función de la forma del
objeto (menor espesor de capa límite cuanta menor resistencia aerodinámica presente la
superficie: ej. forma fusiforme de un perfil alar).
 Objetivo general:

 Aprender la teoría sobre capa limite, utilizando diferentes textos sobre mecánica de
fluidos e hidráulica.

Objetivos específicos:

 Conocer e identificar la naturaleza del flujo de un líquido, ya sea de tipo laminar o

turbulenta.
 Aplicar la teoría resolviendo ejercicios de aplicación.
Una mejor explicación del concepto de capa limite se consigue al suponer la existencia de un
flujo uniforme horizontal, de velocidad constante v0, al cual se interpone un cuerpo agudo,
como se muestra en la figura 1.
El campo de velocidades del flujo se adapta a la forma del cuerpo y modifica su carácter de
uniformidad. Además, las fuerzas viscosas que retardan el movimiento del fluido en una capa
muy delgada (de espesor 0) en el punto en el que la corriente toca al cuerpo, se manifiesta
sobre una capa de espesor creciente al transmitir progresivamente el efecto de filete en filete.
(Sotelo A, 2012)

Figura 1: Concepto de Capa Limite

Fuente: (Sotelo A, 2012)
 La región interior de la capa limite presenta las siguientes características:

a) la velocidad del fluido, en el punto de contacto con la frontera vale 0.

b) el gradiente transversal de velocidades y, por lo mismo, el esfuerzo cortante, tiene valores
máximos en la frontera; disminuye a medida que el punto se aleja de esta.
c) fuera de la capa limite, el gradiente de velocidades es prácticamente cero, porque también lo
es el esfuerzo cortante y los efectos viscosos.
d) fuera de la capa limite las líneas de corriente se conforman, según un flujo con potencial,
sufriendo un ligero desplazamiento hacia el exterior de la pared.

Figura 2: Desarrollo de la capa limite

en la pared interna del conducto.
Fuente: (Crespo M., 2006)
Se va a limitar al movimiento bidimensional, plano, incomprensible, laminar, de viscosidad
constante y con fuerzas máximas que derivan de un potencial. Las ecuaciones que rigen el
movimiento son:

: ….. (1)

(2)

(3)

Dónde: 𝜹 es el espesor típico de la capa limite, L el tamaño típico del cuerpo y U una velocidad
típica de deslizamiento del fluido
 Ec. (12) Ec. (13)

En la superficie de la pared las dos componentes de la velocidad son cero:

Ec. (16)

Y lejos de la pared, la velocidad Vx debe tender al valor correspondiente al deslizamiento del

fluido ideal U:
Ec. (17)

Dado que la velocidad de deslizamiento U correspondientes a:

 (17)

Hay una aparente contradicción entre estas ecuaciones, podemos resolverla de la

forma siguiente: el infinito al que se refiere la Ecuación (*) corresponde a la
coordenada y medida con la escala pequeña 𝛿, mientras que el cero al que se
refiere la solución ideal corresponde a la coordenada y medida con el tamaño del
cuerpo L. Pero al ser L >> 𝛿, hay margen suficiente para poder considerar que un
valor muy grande de 𝑦/𝛿 es un valor muy pequeño de 𝑦/𝐿.
 :

La anterior estimación de orden de magnitud nos permite estimar la fuerza de fricción

ejercida por el fluido sobre un cuerpo perfilado, como el mostrado en la Figura 6.

Figura 6: Esquema de un cuerpo

de forma aerodinámica, su
longitud L y velocidad incidente
U. Fuente: (Crespo M., 2006)

El esfuerzo cortante seria de Y la fuerza:

un orden de magnitud:
…..Ec. (21)
…..Ec. (20)
Donde: B es la envergadura del cuerpo.

Sustituyendo en (21) el valor

de 𝝉 de (20) y el valor de δ …..Ec. (22)
dado en (9), resulta:
Esta fuerza, adimensionalizada como se indica en la Ecuación (30)
resultaría:

…….. Ec. (23)

Es decir, este valor adimensional decrece al aumentar el número de Reynolds. En este

caso al coeficiente CD le llamamos también Cf, pues al no tener en cuenta el
desprendimiento de capa limite, toda la fuerza es debida a la fricción. (Crespo M., 2006)
 Figura 3: Efecto de
desplazamiento en la
capa límite.
Fuente: (White, 2004)

En este caso el campo fluido tendrá una

velocidad uniforme, 𝑈𝑖Ԧ , para el caso de
movimiento ideal. La velocidad de …….. Ec. (24)
deslizamiento será pues constante, y de la
ecuación (15) se tiene:
Debemos pues resolver las Ecuaciones (12) y (13) con la condición (24) y
las condiciones de contorno (16) y (17):
…….. Ec. (25)

…….. Ec. (26)

Y la condición (19):

…….. Ec. (27)

Que expresa que en el borde delantero de la placa la velocidad es constante, o

según veremos, que la capa limite tiene espesor nulo. El problema que acabamos
de plantear es algo similar al de la placa con movimiento impulsivo, de Rayleigh. En
aquel problema el efecto de la aceleración por arrastre del flujo por la placa se iba
difundiendo temporalmente a capas de fluido cada vez más anchas.
Ahora, según se vería, el papel del tiempo lo desempeñara la coordenada x, y el efecto de
frenado de la corriente exterior se va experimentando aguas abajo a capas cada vez más
anchas. En aquel problema encontramos que había una variable de semejanza de la

forma dada en la ecuación:

Por analogía con esa ecuación, y haciendo una equivalencia entre el tiempo del problema de
Rayleigh y el tiempo que se tarda en recorrer la coordenada X, nuestra variable de semejanza
seria:

…….. Ec. (28)

También, por analogía con la …….. Ec. (29)

Ecuación (35), ponemos:
Es nuestro caso, para satisfacer la Ecuación, se define, de acuerdo con la

Ecuación, una función de corriente Ѱ:

…….. Ec. (30)

Para satisfacer la primera de esta ecuación se debe cumplir:

…….. Ec. (31)

Donde f es una función de n. De (31) y (30) se obtiene:

…….. Ec. (32)

…….. Ec. (33)

Como se puede comprobar a partir de las Ecuaciones anteriores.
La Ecuación (34) junto con las condiciones (35) y (36) se puede
integrar numéricamente. Una vez obtenido f, de puede calcular f´, y
obtener la distribución de velocidad según X, que se representa en
la figura 8.

Figura 8: Diagrama de la distribución de

velocidades para la solución de
semejanza de la capa limite. Fuente:
(Crespo M., 2006)
A un valor de n constante una velocidad constante. Esto quiere decir que, a medida que
nos movemos aguas abajo, el espesor de la capa limite crece. Por ejemplo el valor de n
que da que Vx es 0.99U, es 5; si suponemos que prácticamente para ese n la capa
límite ha terminado, su espesor seria (a partir de [28]):

... Ec. (37)

Resultado que se debe comparar con la Ecuación (9), que ahora nos diría (se llama la
atención sobre que la desigualdad δ<<x no se cumple muy cerca del borde de ataque,
donde:
 Y por tanto allí no vale la teoría de la capa límite, lo que
puede introducir dudas sobre la validez de la condición
(27).) δ<<x.

Figura 9: Diagrama que muestra

la interpretación física del
espesor de desplazamiento.
Fuente: (Crespo M., 2006)

Esta forma de definir el espesor de la capa límite es algo

arbitrario, ya que el valor U para la velocidad solo se alcanza en el
infinito. Una forma menos arbitraria de definir el espesor es:

... Ec. (38)

De las Ecuaciones (28), (32) y (38) y de una integración
numérica, a partir de los valores del grafico anterior, se
puede obtener:

... Ec. (39)

O también:

... Ec. (40)

Derivando en (40) respecto a x, respecto a x, y utilizamos las

Ecuaciones (33), (36) y (18) se obtiene:

... Ec. (41)

Al igual que para el fluido másico, se puede definir un espesor
asociado a la pérdida de impulso o cantidad de movimiento a
través de la capa límite:

…... Ec. (42)

De (28) y (32) y una integración numérica, se obtiene:

…... Ec. (43)

Y este espesor, según veremos, está asociado a la fuerza de

fricción sobre la placa. Para calcular la fuerza de fricción, se va a
calcular lo que vale el esfuerzo cortante en la pared:

…... Ec. (44)

Donde se han utilizado las ecuaciones (28) y (32) y que f’’(0) =
0.664/2 = 0.332. Integrando τ a lo largo de la placa, se obtiene la
fuerza sobre una sola cara de la misma por unidad de
envergadura:

…... Ec. (45)

Ecuación que se debe con la (22). La fuerza de fricción

adimensionalizada seria:

…... Ec. (46)

ECUACIÓN INTEGRAL DE KARMAN:

Observamos que la ecuación de la continuidad (1) y de la

contorno (16), la velocidad normal a la placa se puede
expresar como:

…... Ec. (53)

Y por tanto el segundo sumando de (13), dividido por la

densidad e integrado a través de la capa limite, se puede escribir
en la forma:

…... Ec. (54)

Donde hemos efectuado una integración por partes. La
ecuación (13), dividida por la densidad, integrada a través
de la capa límite, habiendo sustituido el segundo sumando
por la anterior expresión, resulta:

…... Ec. (55)

Donde el último termino es el esfuerzo cortante en la pared.

Reagrupando términos en (55), con la ecuación (15) para la presión,
y sumando y restando dentro del integrando el término:

Tenemos (cambiando de signo los 2 miembros):

…... Ec. (56)

Si ahora utilizamos las definiciones (38) y (42) para los espesores
de desplazamiento e impulso, tenemos:
…... Ec. (57)

APLICACIÓN A UNA PLACA PLANA

Supongamos que para la placa plana, antes estudiada, tomamos un
perfil de velocidad lineal como el de la figura 13:
…... Ec. (58)

Que termina a una distancia δ. Esto, si se compara con la

solución exacta, figura 5, es una aproximación muy mala.

Figura 13: Diagrama que muestra la distribución lineal de

velocidades en la capa límite para la aplicación de la
ecuación integral de Karman.
Fuente: (Crespo M., 2006)
La capa límite se hace inestable cuando el número de Reynolds basado en su espesor
y en la velocidad exterior se hace mayor de un cierto valor.
En general el número de Reynolds que marca la
transición de laminar a turbulento depende de:
1. Gradiente de presión reducida, dP/dx. Cuanto menor sea, la capa limite es más estable y
el número de Reynolds de transición se hace mayor.
Los valores negativos de dP/dx hacen estable la capa limite, y los positivos inestable.
Según la ecuación (15), cuanto mayor sea la aceleración de la corriente exterior, dU/dx,
más estable es la capa límite.
2. Nivel de turbulencia que tiene la corriente exterior; si este es grande, se favorece la
transición.
3. Rugosidad de la pared; cuanto mayor es, más se favorece la transición a la turbulencia.
4. Otras: succión de capa limite, transición de calor, etc.
Para una placa plana lisa, con un nivel de turbulencia ambiente del 1% la transición
corresponde aproximadamente a la condición:

…... Ec. (62)

Siendo δ1 el espesor de desplazamiento; lo que
corresponde, de acuerdo con la ecuación (39), a una
distancia critica al borde de ataque de la placa tal que:

…... Ec. (63)

La ecuación (62) se puede interpretar,

aunque solo sea como regla nemotécnica,
como si la capa limite fuese la mitad de la
corriente en un tubo; en aquel caso la
transición correspondía a un numero de
Reynolds 2.000.
Así pues, en general, en la capa limite hay dos
zonas una primera que es laminar y llega
hasta la distancia dada por la ecuación (36) FIGURA 14: Deposito que descarga agua
para una plaza plana, y otra a continuación,
por un tubo, formándose una capa limite
que es turbulenta. La capa limite turbulenta
se caracteriza por ser más gruesa y dar en sus paredes. Fuente: (Crespo M.,
esfuerzos cortantes mayores que la laminar. 2006)
ECUACIONES DE LA CAPA LÍMITE TURBULENTA
Se va a seguir un procedimiento parecido y se van a hacer las
mismas hipótesis que para la capa limite laminar, apartado 1.
Se va a suponer que el flujo es estacionario, bidimensional y
con fuerzas másicas que derivan de una potencial. Las
ecuaciones de la capa limite turbulenta se obtienen de las
ecuaciones de Reynolds,

Con lo que aparecerán unos nuevos términos debidos a los

esfuerzos de agitación turbulenta cuyo orden de magnitud es:

𝜌𝑣𝑖 𝑣𝑗 ~ 𝜌𝑣 ∗2 = 𝜎𝑝 = 𝜏0 ……. (64)

Donde, σρ es el esfuerzo cortante en la pared, ahora denominado τ0, y v ∗
La teoría de la capa limite turbulenta se basa en que los esfuerzos
cortantes turbulentos, al igual que sucedía con los viscosos, son
depreciables en todo el campo fluido, excepto en la capa limite.
Para ello, supondremos que la velocidad de agitación turbulenta
es pequeña comparada con la velocidad media de deslizamiento
del fluido:

V∗ ≪ U ……. (65)

Con la condición (65), al ser además Re grande, podemos

asegurar que fuera de la capa limite los esfuerzos viscosos y
los turbulentos son despreciables, y que por tanto el fluido
se comporta como ideal.
Sin embargo, dentro de la capa limite, que es muy estrecha,
los términos de transporte más grande deben ser del orden
de los convectivos. Si se supone que todas las velocidades
de agitación turbulenta son de igual orden, resultaría que la
contribución más importante seria de los términos de
transporte con derivadas según el eje y. En particular:
∗2 ∗2
𝑣 𝜕 2
𝜕 ′ ′
𝑣
~ 𝑣′𝑥 ≪ 𝑣𝑥 𝑣𝑦 ~
𝐿 𝜕𝑥 𝜕𝑦 𝛿
Con esta condición y las utilizadas para la capa limite laminar,
apartado 1, la ecuación
Proyectada según x nos daría:

𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑃 𝜕𝜏 ……. (66)

𝜌 𝑣𝑥 + 𝑣𝑦 + =
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦

Donde el esfuerzo cortante es:

𝜕𝑣𝑥 ′ ′
𝜏=𝑢 − 𝜌𝑣𝑥 𝑣𝑦 ……. (67)
𝜕𝑦
Comparando en la ecuación (66) el término correspondiente al
esfuerzo cortante turbulento con el término convectivo, es fácil
ver que si el término viscoso es despreciable y los otros dos son
del mismo orden, se debe cumplir que:

……. (68)
Lo que confirma que el espesor de la capa limite es
pequeño. La ecuación de cantidad de movimiento,
proyectada según y, nos daría el orden de magnitud de la
variación transversal de presión que resulta muy pequeña
comparada con la longitudinal, dada por la ecuación (10):

……. (69)

De esta relación se deduce que la ecuación (14) sigue valiendo. Por lo tanto, el
sistema de ecuaciones estará ahora formado por las ecuaciones de continuidad
(12), y conservación de cantidad de movimiento, (66), que sustituye a la (13) de
caso laminar, y la (14), que expresa que la presión no cambia a través de la capa
limite. Las condiciones de contorno del caso laminar siguen valiendo, y en
particular el flujo fuera de la capa limite sigue siendo ideal. La principal diferencia
con el caso laminar estriba en el término turbulento de (67), que no aparecía
antes y que además no se sabe cómo calcularlo, aunque se den algunos criterios
para su estimación turbulentas, que consisten esencialmente en procedimiento
para modelar ese término. (Crespo M., 2006)
CAPA LÍMITE TURBULETA PARA PLACA PLANA LISA

Está fuera de los objetivos de este texto deducir o intentar resolver las
ecuaciones de capa límite de flujo turbulento. Las expresiones para la forma del
perfil de capa límite y otras propiedades de la capa límite turbulenta se
obtienen empíricamente (o, en el mejor de los casos, semi-empíricamente),
porque no se pueden resolver analíticamente las ecuaciones de capa límite para
flujo turbulento. Note también que los flujos turbulentos son inherentemente
no-estacionarios y que la forma del perfil de velocidad instantánea varía con el
tiempo (Fig. 15). (Cengel & Cimbala , 2006)

Figura 15: Ilustración de la condición

no estacionaria de una capa límite
turbulenta; las delgadas líneas
negras son perfiles instantáneos, y
la línea gruesa azul es un perfil
promediado en el tiempo largo.
Fuente: (Cengel & Cimbala , 2006)
Con la ecuación (66) se puede hacer exactamente
el mismo tratamiento que con la (13), y siguiendo
los pasos (53) a (57) se llegaría a la ecuación
integral de Karman (57). Particularizando para
placa plana:
U= constante
……. (70)
Esa ecuación nos daría:
En principio se podrían utilizar las leyes logarítmicas de la pared,
para calcular el espesor de cantidad de movimiento e introducirlo
en la ecuación (70), y así obtener la evolución de la capa limite. Sin
embargo, los cálculos necesarios son largos y los resultados difíciles
de aplicar. Aunque sea menos riguroso, vamos a utilizar el perfil de
velocidades en forma de potencias dado en

……. (71)
Según se indicó allí mismo, esta ley no tiene ningún
fundamento teórico y no se puede usar en forma directa
para calcular τ_o. Para ello, se debe usar la ecuación

Junto con la :

(Teniendo en cuenta que lo que allí llamábamos σ_p ahora

es τ_0):

……. (72)
El espesor de la capa limite es equivalente al radio del tubo,
y la velocidad máxima en el centro del tubo equivale a la
velocidad exterior a la capa limite. De la ecuación

Con n = 7:
……. (73)

Y de (72) y (73) se tiene:

……. (74)
Llevando la ecuación (71) a la (42) e integrando, se obtiene
una relación entre los espesores δ y de cantidad de
movimiento δ_2:
𝛿
𝛿2 = 7 ……. (75)
72
Y llevando (75) y (74) a la ecuación integral de Karman, (70), se
obtiene una ecuación para δ:

……. (76)

Que podemos integrar con la condición inicial:

x=0 ,δ=0 ………… (77)

Resultado:

……. (79)
……. (78)
Que nos confirma que si el Re basado en la longitud de la
plaza es suficientemente grande, el espesor de la capa
limite es pequeño comparado con esa longitud.
Repitiendo los cálculos que llevan a las ecuaciones (45) y
(46), se llega a las siguientes expresiones para la fuerza
sobre la placa:

……. (80)

……. (81)c

Este mismo resultado se puede obtener utilizando el espesor de cantidad

de movimiento mediante las ecuaciones (50), (75) y (78). El margen de
validez de este resultado es:
 La siguiente formula aproximada de Schlichting:

……. (82)

Vale hasta un valor de 𝑅𝑒 de 109 . Se pueden comparar las ecuaciones (81) y (82)
para diferentes 𝑅𝑒 y ver que dan resultados parecidos. Si la capa limite laminar de la
zona inicial ocupase una parte importante de placa, habría que hacer la corrección
correspondiente. (Crespo M., 2006)
Ejercicio 05
Supóngase que una capa limite sobre una placa lisa, plana, es laminar al principio y luego se
hace turbulenta a un número critico de Reynolds de 5 x 105 . Si tenemos una placa de 3m de
largo y 1m de ancho, y si circula aire a 20°C y presión atmosférica normal alrededor de esta
placa con una velocidad de 30 m/s, ¿Cuál será el coeficiente de resistencia promedio Cf para
la placa? También, ¿Cuál será la resistencia al corte total de un lado de la placa y cual la
resistencia debida a la parte turbulenta y la parte laminar de la capa límite?
SOLUCIÓN
𝑈02
La resistencia total es 𝐹𝑠 = 𝐶𝑓 𝐵𝐿𝜌 , donde
2

de la ecuación:

0.523 1520
𝐶𝑓 = 2 −
𝑙𝑛 0.06𝑅𝑒𝐿 𝑅𝑒𝐿
𝑈𝐿
Además, 𝑅𝑒𝐿 = 𝑉
, o sea

30𝑚/𝑠 × 3𝑚
𝑅𝑒𝐿 = −5 = 5.96 × 106
1.51 10
Entonces, resolviendo la ecuación:

0.523 1520
𝐶𝑓 = −
𝑙𝑛2 0.06𝑅𝑒𝐿 𝑅𝑒𝐿
Tenemos:
𝐶𝑓 = 0.00320 − 0.00026 = 0.00294
La resistencia total se calcula como:

𝑈02 302
𝐹𝑠 = 𝐶𝑓 𝐵𝐿𝜌 = 0.00294 × 1 × 3 × 1.2 × = 4.76𝑁
2 2

Entonces 𝑥𝑐𝑟 se determina como:

𝑈𝑥𝑐𝑟
= 500000
𝑣
500000×1.51×10−5
O bien, 𝑥𝑐𝑟 = = 0.252𝑚
30
Por lo tanto, la resistencia laminar será

1.33 302
𝐹𝑠, 𝑙𝑎𝑚 = × 1 × 0.252 × 1.2 × = 0.256𝑁
5 × 105 1/2 2

Entonces, 𝐹𝑆, 𝑡𝑢𝑟𝑏 = 4.76𝑁 − 0.26𝑁 = 4.50𝑁