Introdución a la capa límite laminar bidimensional y

estacionaria

M. Rodríguez

1 Introducción

En los movimientos a altos números de Reynolds (basado en la longitud característica del mo-
vimiento), los efectos viscosos son despreciables. Los efectos de conducción de calor también lo
son si el producto del Reynolds por el Prandtl es grande. Esto lleva implícito el despreciar los
términos de mayor orden en las derivadas de la velocidad y de la temperatura, de modo que no
pueden imponerse todas las condiciones de contorno. Como consecuencia de ello, las condiciones
de contorno en el movimiento de un fluido ideal en presencia de una pared se reducen a decir
que la velocidad es tangente a la pared (si no hay paso de masa a través de dicha pared). Sin
embargo, dentro de la aproximación de un fluido como medio continuo, sabemos que la velocidad
de un fluido en contacto con una pared es igual a la velocidad de la pared y que la temperatura
del fluido debe coincidir con la de la pared (si a través de dicha pared no hay paso de masa ni
reacción química o evaporación en su superficie).

Para poder imponer todas las condiciones de contorno es necesario que los términos viscosos
y de conducción de calor sean tan importantes como los convectivos. Sin embargo, si se utiliza la
longitud característica ` del movimiento, el número de Reynolds es tan alto que estos términos
serían despreciables. Por lo tanto, es evidente que cerca de las paredes (donde se deben imponer
las condiciones de contorno) la velocidad (y temperatura) sufre variaciones del orden de ella
misma en distancias δ ¿ `. El orden de magnitud de δ se determina de la condición de que
los términos viscosos sean tan importantes como los convectivos, que son los más importantes
cuando el número de Reynolds es alto.

El primero en indicar la existencia de una zona en la que los efectos viscosos son importantes,
a pesar de que el número de Reynolds del movimiento sea alto, fue Prandtl en 19041 . La idea
de Prandtl de una región donde los efectos viscosos son importantes puso en claro multitud
de fenómenos que hasta entonces no habían obtenido explicación satisfactoria. En particular
explicó el porqué la teoría de fluidos ideales (altos números de Reynolds) proporciona buenos
resultados cuando se quiere determinar la sustentación o fuerza normal a la dirección de la
corriente incidente sobre un obstáculo, y sin embargo esta teoría ideal es incapaz de proporcionar
la resistencia (o componente de la fuerza en la dirección de la corriente incidente). También
explico el desprendimiento de la capa límite en cuerpos romos (en general con gradientes adversos
de presión), y como consecuencia de ello la existencia de una resistencia de forma (que no depende
de la viscosidad, pero es causada por ella).
1
Prandtl, L., ”Über Flüssigkeitsbewegung bei sehr kleiner Reibung”, Proc. III Intern. Math. Congress,
Heidelberg (1904).
La traducción al español puede encontrarse en ”Versión Crítica en Español del Trabajo de Ludwig Prandtl
sobre el Movimiento de Fluidos con Viscosidad muy Pequeña”, Ingeniería Aeronáutica y Astronáutica, No 328,
Julio 1992, por M. Rodríguez y R. Martínez-Val.

1
 En cuerpos fuselados la capa límite no se desprende más que en la parte final del cuerpo
(como es el caso del perfil de la figura 1), formando una estela muy delgada. En este caso
la resistencia es prácticamente toda ella debida al rozamiento (viscosidad). Sin embargo, en
un cuerpo romo (figura 2), la capa límite se desprende generando una estela amplia, donde la
presión toma valores próximos a la estática en el infinito, mientras que en el borde de ataque
la presión es del orden de la de remanso, generándose una fuerza de resistencia del orden de la
diferencia de presiones por el área frontal. Esta fuerza, aunque originada por el desprendimiento
de la capa límite, que a su vez tiene su origen en la viscosidad, no depende de dicha viscosidad
si el número de Reynolds es muy alto.

Figura 1.- Capalímite adherida al perfil.

Figura 2.- Estela detrás de un cilindro circular.

En el caso de un cuerpo fuselado, en el que la corriente no está lo suficientemente ”alineada”

con su geometría, puede desprenderse la corriente como en el caso de un cuerpo romo. Este es
el caso del perfil de la figura 1 cuando el ángulo de ataque es elevado (véase figura 3).

En el movimiento de los fluidos alrededor de cuerpos o en presencia de paredes, si el número

de Reynolds es grande, hay una región exterior donde los efectos viscosos y de conducción de
calor son despreciables. Las ecuaciones se reducen a las de Euler, a las que no se le pueden
imponer todas las condiciones de contorno. La solución de Euler proporciona una velocidad
tangente a la pared y que varía a lo largo de ella, que denominaremos corriente exterior. Esta
denominación se debe a que es la corriente que ”se ve” en el exterior de la capa límite.

2
 Figura 3.- Perfil con capa límite desprendida.

Dado que el espesor de la capa límite es pequeño comparado con la dimensión longitudinal, `,
característica del movimiento, consideraremos el eje de coordenadas x a lo largo de la superficie
del cuerpo y el eje y perpendicular y transversal a la capa. Dado que salvo algunas excepciones,
el radio de curvatura de la superficie es grande frente al espesor de la capa límite, consideraremos
los ejes x − y como un sistema de coordenadas cartesianas.

2 Ecuaciones de la capa límite bidimensional compresible y es-

tacionaria

Las ecuaciones de la capa límite bidimensional estacionaria de un fluido compresible, pueden

escribirse en la forma
∂ (ρu) ∂ (ρv)
+ = 0, (1)
∂x ∂y
µ ¶
∂u ∂u dpe ∂ ∂u
ρu + ρv =− + μ , (2)
∂x ∂y dx ∂y ∂y
µ ¶ µ ¶
∂h0 ∂h0 ∂ ∂T ∂ ∂u
ρu + ρv = k + μu . (3)
∂x ∂y ∂y ∂y ∂y ∂y

La ecuación (1) es la ecuación de la continuidad de la que se deduce que

δ
v ∼ u ¿ u,
`
por ser δ ¿ `. Esto último se deduce de la ecuación de cantidad de movimiento, en la que los
términos convectivos son del orden del viscoso si
r
δ μ 1
∼ ∼√ ¿ 1.
` ρue ` Re
Dado que δ ¿ `, la ecuación de cantidad de movimiento en la dirección del eje y se reduce a
decir que la presión no cambia con y, de modo que adquiere su valor en el exterior de la misma.
Por esa razón se ha sustituido ∂p/∂x = dpe /dx, siendo pe la presión en el exterior de la capa y
conocida de la solución exterior de Euler.

En la ecuación (3) de la energía h0 es la entalpía de remanso

1¡ 2 ¢ 1
h0 = h + u + v2 ≈ h + u2 .
2 2

3
Como la presión es constante a través de la capa límite, la densidad se determina mediante la
ecuación de estado ρ = pe /Rg T y la entalpía mediante la relación h = cp T en el caso de los
gases.

Para que el término de conducción sea comparable al convectivo en la ecuación (3) de la

energía, es necesario que
δT 1
∼√ ,
` ReP r
siendo δT el espesor de la capa donde los efectos de conducción de calor son importantes. Esta
capa se denomina capa límite térmica. El número de Prandtl P r = μcp /k es de orden unidad en
el caso de los gases, de modo que el espesor de ambas zonas, viscosa y térmica, es comparable.
Puede no ocurrir así en el caso de líquidos, dando lugar a una capa térmica mucho más gruesa
que la viscosa si P r ¿ 1 (pero P rRe À 1), o mucho más delgada si P r À 1. Para los líquidos,
la disipación viscosa es despreciable.

Las condiciones de contorno para integrar las ecuaciones anteriores son

∂T
En y = 0 : u = 0; v = vs ; h0 = hp ( = 0),
∂y
En y → ∞ : u = ue ; h0 = h0e , (4)
En x = 0 : u = ui (y) ; h0 = h0i (y).

En las condiciones anteriores se ha supuesto una velocidad de succión o soplado, vs , distinta de

cero en la pared. En la condición de contorno para la temperatura en la pared se han supuesto
dos posibilidades: la temperatura del gas coincide con la temperatura de la pared conocida
(hp = cp Tp ), o bien la pared está aislada térmicamente (∂T /∂y = 0).

Las ecuaciones para el caso incompresible toman la forma simplificada

∂u ∂v
+ = 0, (5)
∂x ∂y

∂u ∂u dpe ∂2u
u +v =− +ν 2, (6)
∂x ∂y ρdx ∂y
∂T ∂T ν ∂2T
u +v = , (7)
∂x ∂y P r ∂y 2
don P r = μc/k es el número de Prandtl. Las condiciones de contorno son

∂T
En y = 0 : u = 0; v = vs ; T = Tp ( = 0),
∂y
En y → ∞ : u = ue ; T = Te , (8)
En x = 0 : u = ui (y) ; T = Ti (y).

2.1 Desprendimiento de la capa límite

Como se ha indicado anteriormente, existe la posibilidad de que el fluido ”frenado” por la capa
límite pueda ser transportado a la corriente exterior y arrastrado por ella. Cuando este ocurre se
dice que la capa límite se ha desprendido una estela como la de las figuras 2 y 3. Cuando existe
una región con presiones crecientes a lo largo de la pared (dpe /dx > 0 ó gradiente adverso de
presiones) las partículas fluidas frenadas en la capa límite no pueden, en general, penetrar muy

4
lejos en la zona de gradiente adverso de presión a causa de su baja energía cinética, y la capa
límite es deflectada separándose de la pared. El punto de separación corresponde a esfuerzo nulo
en la pared, lo que equivale a decir que (∂u/∂y)y=0 = 0.

Figura 4.- Esquema de los perfiles de velocidad en la capa límite a lo largo de la pared. S es el
punto de separación. P I es el punto de inflexión en el perfil de velocidades.

En la figura 4 se muestra la distribución de velocidades cerca del punto de separación y, si el

esfuerzo se anula en la pared, hay un punto de inflexión en el perfil de velocidades.

Si particularizamos la ecuación (2) en la pared con vs = 0 y μ constante, se tiene

µ 2 ¶ µ 3 ¶
∂ u dpe ∂ u
μ = de modo que = 0. (9)
∂y2 y=0 dx ∂y3 y=0

En las proximidades de la pared la curvatura del perfil de velocidades depende sólo del gra-
¡diente de¢presiones externo. Si el gradiente de presiones es favorable (dpe /dx < 0) entonces
∂ 2 u/∂y2 y=0 < 0 y por lo tanto lo es en toda la capa límite, tal como se muestra en la figura
5.

Figura 5.- Distribución de velocidades, esfuerzos y curvatura en la capa límite cuando el

gradiente de presiones es favorable.

5
 ¡ ¢
Sin embargo, cuando el gradiente de presiones es adverso (dpe /dx > 0) se tiene ∂ 2 u/∂y2 y=0 >
0, pero a lejos de la pared se tiene ∂ 2 u/∂y2 < 0, de modo que en algún punto del interior debe
tenerse ∂ 2 u/∂y 2 = 0. Éste es el punto de inflexión del perfil de velocidades en la capa límite,
tal como se muestra en la figura 6.

Figura 6.- Distribución de velocidades, esfuerzos y curvatura en la capa límite cuando el

gradiente de presiones es adverso. P I es el punto de inflexión.

En la región de gradiente adverso de presiones, el perfil de velocidades en la capa límite

siempre presenta un punto de inflexión. Dado que cuando el esfuerzo es nulo en la pared, el
perfil de velocidades tiene un punto de inflexión (véase figura 4), la separación solo aparece
ocurre cuando el gradiente de presiones es adverso (a menos de que exista succión o soplado
como veremos a continuación).

2.1.1 Efecto de la succión o soplado

Si ahora suponemos que la velocidad vs es distinta de cero, la ecuación equivalente a la (5) sería
µ 2 ¶ µ ¶
∂ u dpe ∂u
μ = + ρv s . (10)
∂y 2 y=0 dx ∂y y=0

En la ecuación (6) puede observarse que la velocidad vs hace el mismo papel que un gradiente
adverso de presiones si es positiva (soplado), mientras que si es negativa (succión) hace el papel
de un gradiente favorable. Así, si se tiene un gradiente adverso de presiones, pero una velocidad
de succión adecuada, puede conseguirse que la capa límite no se desprenda.

En la figura 7 se muestra un canal convergente divergente. El chorro esta desprendido de las

paredes superior e inferior a causa del gradiente adverso de presiones en estas paredes, como se
observa en la foto superior. Cuando se introduce una succión en la pared superior, la corriente se
queda adherida a esta pared y se desprende de la inferior, como se observa en la foto intermedia.
Por último, en la foto inferior hay succión en ambas paredes, la corriente se queda adherida y
solo se desprende cuando desaparece la succión.

6
Figura 7.- Canal convergente divergente. En la foto superior se muestra el chorro desprendido.
En la foto intermedia hay una succión en la pared superior. En la foto inferior hay succión en
la pared superior e inferior.

2.2 Resistencia de forma

Si no hubiese desprendimiento de la capa límite, lo que ocurre en cuerpos fuselados como el

perfil de la figura 1, la resistencia del cuerpo seria debida a los esfuerzos de fricción en la pared.
Su orden de magnitud estaría dado por
µ ¶
∂u U μ ` ρU 2
τp = μ ∼ μ ∼ ρU 2 ∼√ ,
∂y y=0 δ ρU ` δ Re
donde U es el orden de magnitud de la velocidad exterior. Si Am es área del cuerpo bañada por
el fluido, la resistencia de fricción sería
ρU 2 Am
Dv ∼ τp Am ∼ √ .
Re

En el caso de un cuerpo romo como el de la figura 2, o un cuerpo fuselado a un ángulo de

ataque inadecuado como el de la figura 3, la capa límite se desprende generando una estela de
dimensión transversal comparable a la dimensión del cuerpo perpendicular a corriente incidente.
En la zona desprendida la presión es parecida a la presión aguas arriba no perturbada e incluso
menor, mientras que en la parte frontal la presión es del orden de la presión de remanso. Esta
diferencia de presiones es del orden de ρU 2 , de modo que la fuerza de resistencia debida a esta
diferencia de presiones es
Dp ∼ ρU 2 AF ,
siendo AF el área frontal del cuerpo (área proyectada en la dirección de la corriente incidente).
A esta resistencia se le denomina resistencia de presión o de forma.

Si se comparan ambas resistencias se obtiene

Dv Am /AF
∼ √ ,
Dp Re

7
que para Am /AF ∼ 1 (caso de una esfera, por ejemplo), la resistencia debida a la fricción es
muy pequeña comparada con la de forma ya que el número de Reynolds es muy alto.

2.3 Espesor de la capa límite

Debido a que la velocidad u tiende asintóticamente al valor exterior ue , el espesor de la capa

límite estaría definido de una forma ambigua, ya que dependería del grado de aproximación de
la velocidad u a ue . Por ejemplo se podría definir el espesor δ de modo que u (x, δ) = 0.99u (x),
pero también se podría definir de modo que u (x, δ) = 0.995u (x), con lo que obtendríamos un
valor mayor para el espesor.

Sin embargo se pueden definir espesores de la capa límite de una forma más física, y están
asociados a cada una de las ecuaciones del movimiento. Así, el espesor de desplazamiento, δ ∗ ,
esta asociado a la ecuación de la continuidad. La masa que atraviesa una sección de la capa
límite es Z ∞
ρudy,
0
y si no existiese capa límite, la masa sería
Z ∞
ρe ue dy,
0
y para que ambas magnitudes fuesen idénticas sería necesario desplazar la pared una cantidad
δ ∗ tal que
Z ∞ Z ∞ Z ∞ Z δ∗
ρudy = ρe ue dy = ρe ue dy − ρe ue dy,
0 δ∗ 0 0
de modo que Z ∞µ¶
∗ ρu
δ = 1− dy, (11)
0 ρe ue
ya que ρe y ue no dependen de y. Esta definición del espesor ya no depende del grado de
aproximación que se quiera con respecto a la velocidad exterior.

Asociado a la ecuación de cantidad de movimiento se tiene el espesor de cantidad de movi-

miento, δ ∗∗ , que es la distancia adicional a δ ∗ que hay que desplazar la pared para la igualdad
de cantidad de movimiento, esto es
Z ∞ Z ∞ Z ∞ Z δ∗ +δ∗∗
2 2 2
ρu dy = ρe ue dy = ρe ue dy − ρe u2e dy,
0 δ ∗ +δ ∗∗ 0 0

de modo que Z ∞µ ¶
∗ ∗∗ ρu2
δ +δ = 1− dy,
0 ρe u2e
y sustituyendo el valor de δ ∗ dado en (7) se tiene
Z ∞ µ ¶
∗∗ ρu u
δ = 1− dy, (12)
0 ρe ue ue
que representa el déficit de cantidad de movimiento.

Se pueden definir también espesores asociados a las distintas formas de la ecuación de la

energía. Por ejemplo, asociado a la entalpía de remanso sería
Z ∞ µ ¶
ρu h0 − hp
∆= 1− dy. (13)
0 ρe ue h0e − hp

8
 La evolución del espesor de cantidad de movimiento, δ ∗∗ , a lo largo de la capa límite está
relacionada con el coeficiente de fricción, mientras que la evolución de ∆ está relacionado con el
coeficiente de trasmisión de calor, como se verá en lecciones posteriores.

3 Temperatura de recuperación

La ecuación (3) de la energía puede escribirse en la forma

µ ¶ "µ ¶ ¡ ¢#
∂h0 ∂h0 ∂ k ∂h0 ∂ k ∂ u2 /2
ρu + ρv = + μ− ,
∂x ∂y ∂y cp ∂y ∂y cp ∂y

y, teniendo en cuenta que el número de Prandtl es P r = μcp /k, la ecuación anterior puede
escribirse en forma
µ ¶ ∙ µ ¶¸
∂h0 ∂h0 ∂ ∂h0 1 − Pr ∂ ∂ 1 2
ρu + ρv = μ + μ h0 − u . (14)
∂x ∂y ∂y ∂y P r ∂y ∂y 2
Cuando el número de Prandtl es igual a la unidad, la ecuación anterior toma la forma
µ ¶
∂h0 ∂h0 ∂ ∂h0
ρu + ρv = μ . (15)
∂x ∂y ∂y ∂y

En el caso de que la pared esté aislada térmicamente, la condición de flujo de calor nulo en la
pared se traduce en (∂h0 /∂y)y=0 = 0, ya que u (∂u/∂y) = 0 por ser u = 0 en y = 0. Por lo tanto,
la solución h0 = h0e satisface la ecuación (11) y las condiciones de contorno (∂h0 /∂y)y=0 = 0 y
h0 = h0e en y → ∞. Esto significa que la temperatura de la pared coincide con la de remanso
de la corriente exterior, esto es
µ ¶
γ−1 2
Tp = T0e = Te 1 + Me . (16)
2

Debido a que el número de Prandtl no es la unidad, si no un poco menor, la ecuación (10)

nos indica que la temperatura de la pared va a ser un poco menor que la de remanso
µ ¶
γ−1 2
Tp = Te 1 + < Me , (17)
2
donde < es el factor de recuperación, función del número de Prandtl, que es menor que la unidad
pero muy próximo a uno.

4 Analogía de Reynolds

En la capa limite de una placa plana la corriente exterior es uniforme, de modo que ue es
constante, lo mismo que la presión y la temperatura. Por lo tanto dpe /dx = 0 y h0e es también
constante. Si la placa está a temperatura constante Tp (hp = cp Tp ) y el número de Prandtl es
la unidad, la ecuación (11) de la energía puede escribirse en la forma
µ ¶
∂θ ∂θ ∂ ∂θ
ρu + ρv = μ , (18)
∂x ∂y ∂y ∂y

9
donde
h0 − hp
θ= .
h0e − hp
Esta ecuación ha de integrarse con las condiciones

θ = 0 en y = 0 y θ = 1 en y → ∞. (19)

Si ahora escribimos la ecuación (2), sin gradiente de presiones, utilizando la variable ϑ = u/ue ,
la ecuación toma la forma µ ¶
∂ϑ ∂ϑ ∂ ∂ϑ
ρu + ρv = μ , (20)
∂x ∂y ∂y ∂y
que ha de integrarse con las condiciones

ϑ = 0 en y = 0 y ϑ = 1 en y → ∞. (21)

Como puede verse, la ecuación (14) con las condiciones (15) es idéntica a la (16) con las
condiciones (17) . Por lo tanto θ (x, y) ≡ ϑ (x, y) 2 .

El flujo de calor en la placa esta dado por

µ ¶ µ ¶
∂T k ∂θ
qp = −k = − (h0e − hp ) , (22)
∂y y=0 cp ∂y y=0

mientras que el esfuerzo en la pared es

µ ¶ µ ¶
∂u ∂ϑ 1
τp = μ = μue = Cf ρe u2e , (23)
∂y y=0 ∂y y=0 2

donde Cf es el coeficiente de fricción.

De la ecuación (18) se obtiene

µ ¶
∂θ cp qp
= ,
∂y y=0 k (hp − h0e )

y de (19) µ ¶
∂ϑ 1 ρe ue
= Cf .
∂y y=0 2 μ
De la igualdad de ambas derivadas se tiene
1 qp
Cf = = St, (24)
2 ρe ue (hp − h0e )

ya que μcp /k = P r = 1. La ecuación (20) es la analogía de Reynolds, mediante la cual se obtiene

que el número de Stanton, St, es igual a la mitad del coeficiente de fricción.

El flujo de calor en la pared toma la forma

µ ¶
1 1 1 2
qp = Cf ρe ue (hp − h0e ) = Cf ρe ue hp − he − ue .
2 2 2

2
Tambien sería necesario que θi (y) ≡ ϑi (y).

10