Capa Límite PDF
Capa Límite PDF
Capa Límite PDF
Captulo 10
Capa Lmite
Un cuerpo que este inmerso en un flujo experimenta una fuerza resultante debido a la accion
entre el flujo y el cuerpo. Esta es la fuerza resultante de los esfuerzos de corte en la pared del
cuerpo (w ) y de los esfuerzos normales a la superficie. La resultante de las fuerzas en direccion
horizontal se denomina Arrastre (D) y las fuerza con direccion vertical se denomina Empuje o
sustentacion (L).
La magnitud de estas fuerzas dependera de la forma que tome el flujo alrededor del cuerpo y por
lo tanto de la forma del cuerpo, de las condiciones del flujo y de la posicion relativa del cuerpo
con respecto al flujo. La figura 10.2 muestra el flujo alrededor de 3 cuerpos con distinta forma.
Separacin
Capa lmite Capa lmite
(a) Placa plana delgada. (b) Cuerpo aerodinamico. (c) Cuerpo obstructor.
La influencia sobre el flujo para el caso de la placa plana delgada paralela al flujo es mnima y
las lneas de corriente tenderan a ser paralelas a la placa. Alrededor de un cuerpo aerodinamico
el flujo que se establece es tal que las lneas de corriente se cierran detras del cuerpo. Alrededor
de un cuerpo obstructor por el contrario, las lneas de corriente no son capaces de cerrarse detras
del cuerpo generando detras de este lo que se conoce como estela.
10.1 Desarrollo de la capa lmite 117
La figura 10.3 muestra el flujo alrededor de una placa plana paralela al flujo para distintos
numeros de Reynolds. Se puede ver de esta figura que a medida que aumenta Re disminuye la
region donde los efectos viscosos son importantes y por lo tanto tambien su influencia en el flujo
externo.
Efectos viscosos
despreciables
Re=Ul/n=0.1 u<0.99U Efectos viscosos
Efectos viscosos despreciables Efectos viscosos
importantes Re=10
importantes
y y
l l
Capa Lmite x Capa Lmite x
u<0.99U
Lneas de corriente
desviadas levemente
Lneas de corriente desviadas
U U
considerablemente U U
y l
Efectos viscosos
despreciables Efectos viscosos
7 Capa Lmite importantes
Re=10
x
d<<l u<0.99U
Lneas de corriente
inalteradas
U
U
Figura 10.3: Flujo alrededor de una placa plana para diferentes Re.
Se denomina Capa lmite a la region alrededor de un cuerpo en la cual los efectos viscosos ( ) no
son despreciables. Como se vio anteriormente los esfuerzos de corte estan asociados a gradientes
de velocidad por lo que en la capa lmite los gradientes de velocidad no son despreciables. Se
debe recordar que las partculas fluidas en contacto con un cuerpo tienen la misma velocidad del
cuerpo. Lo anterior indica que existe una diferencia de velocidades entre el contorno del cuerpo
y el flujo libre lejos del cuerpo. Fuera de la capa lmite se puede considerar el flujo como ideal.
En este captulo se vera, en una primera parte, el desarrollo teorico-emprico de las ecuaciones
que gobiernan el flujo dentro de la capa lmite sobre una placa plana delgada paralela al flujo. En
una segunda parte se analizaran los efectos que la curvatura de los cuerpos produce sobre el flujo
y su influencia en la interaccion dinamica entre el flujo y el cuerpo (lease sustentacion y arrastre).
Se considerara, para efectos de estos apuntes, que el flujo es incompresible y permanente.
La figura 10.4 muestra el desarrollo de la capa lmite sobre una placa plana inmersa en un flujo
para un numero de Reynolds elevado. es el espesor de la capa lmite, es decir, es el lmite de
la region donde los esfuerzos de corte no son despreciables. El numero de Reynolds se define
como:
U x
Re = . (10.1)
C. Gherardelli U. de Chile
10.1 Desarrollo de la capa lmite 118
En una primera parte se desarrolla la capa lmite laminar (x pequeno Re pequeno). En esta
region el flujo es laminar por lo que las partculas se encontraran sometidas a esfuerzos de corte
laminar y no existira mezcla entre las capas. El espesor de la capa lmite aumenta con x
debido al flujo que entra en esta region desde la corriente libre. Como Re es una funcion de la
posicion x sobre la placa, este aumenta con x. Lo anterior indica que para una placa dada y
una velocidad de la corriente libre U dada siempre se alcanzara el regimen turbulento siempre y
cuando la placa sea lo suficientemente larga. Por lo tanto, si la placa es lo suficientemente larga,
existira un punto de transicion (en realidad existe una zona de transicion) donde el regimen se
torna turbulento. La aparicion de un regimen turbulento esta asociado a un aumento notable
en el espesor de la capa lmite. En esta region las partculas estaran sometidas a deformaciones
en todas direcciones y existira mezcla o difusion entre las distintas capas del fluido.
Partcula
fluida
En la zona cercana a la placa las velocidades relativas entre el flujo y la placa son pequenas
generando una zona donde el flujo es laminar. Esta region se denomina sub-capa laminar. Se
puede ver tambien de la figura 10.4 que las partculas que pasan fuera de la capa lmite no se
encuentran sometidas a deformaciones y, por lo tanto, tampoco a esfuerzos de corte.
En cuerpos curvados existe una aceleracion de las partculas fluidas en la parte anterior del
cuerpo hasta el punto de mayor espesor (lneas de corriente se juntan). Pasado este punto el
flujo se desacelera (lneas de corriente se abren) y aumenta la presion. Si la curvatura del cuerpo
es muy grande, y dependiendo de las condiciones del flujo, se puede generar un flujo inverso y lo
que se denomina desprendimiento de la capa lmite (seccion 10.6). Detras del cuerpo se generara
un zona denominada estela (Fig. 10.5).
Capa lmite
Cuerpo
Pto.de Zona de
desprendimiento torbellinos o estela
C. Gherardelli U. de Chile
10.2 Espesor de capa lmite 119
Hasta aqu se ha definido la capa lmite como la region del flujo donde los esfuerzos viscosos, y
por lo tanto los gradientes de velocidad, son importantes. Ya que la aproximacion de la velocidad
del flujo entre el cuerpo y el flujo libre es asintotica, no existe en realidad un lmite determinado
para definir el espesor de la capa lmite. Debido a lo anterior se ha definido (en forma arbitraria)
el espesor como la distancia a la pared donde la velocidad es igual a un 99% la velocidad de
la corriente libre (Fig. 10.6).
u(y)
d=d(x)
Otra forma de definir la capa lmite es considerar el hecho de que la existencia de un gradiente
de velocidades en la region cercana al cuerpo tiene como consecuencia una disminucion tanto
del caudal masico como de la cantidad de movimiento del flujo que pasa por esta region si se le
compara con un flujo no viscoso. Este tipo de definicion elimina la arbitrariedad de la definicion
anterior. Se define el espesor de desplazamiento (ec. 10.3) como la distancia a la que habra
que desplazar el contorno del cuerpo si se supone que no existe roce pero se mantiene el caudal
masico real.
Z
q = (U u)dy = U (10.2)
0
Z
u
= 1 dy . (10.3)
U
0
y y
U U
u(y)
d
A
d*
x x
C. Gherardelli U. de Chile
10.3 Capa lmite laminar 120
Se define el espesor de cantidad de movimiento (ec. 10.5) como la distancia a la que habra
que desplazar el contorno del cuerpo suponiendo que no hay roce pero se mantiene el flujo de
cantidad de movimiento del flujo viscoso (Fig. 10.7), es decir:
Z
2
U = u(U u)dy (10.4)
0
Z
u u
= 1 dy . (10.5)
U U
0
y
U
x
Figura 10.8: Espesor de cantidad de movimiento .
= (x)
= (x)
= (x)
Las ecuaciones que gobiernan el flujo viscoso en la capa lmite son las ecuaciones de Navier
Stokes. Consideraremos en el siguiente desarrollo un flujo bidimensional, permanente y laminar.
Las ecuaciones de Navier Stokes para este caso son
!
u u 1 p 2u 2u
u +v = + + 2 , (10.6)
x y x x2 y
!
v v 1 p 2v 2v
u +v = + + . (10.7)
x y y x2 y 2
C. Gherardelli U. de Chile
10.3 Capa lmite laminar 121
u v
+ = 0. (10.8)
x y
2. ( /x) 1.
3. (/x) 1.
4. v u.
5. x y .
Alguna de estas aproximaciones son solo validas para numeros de Reynolds altos (Re > 1000).
Introduciendo las hipotesis anteriores el sistema de ecuaciones se reduce a
u u 1 p 2u
u +v = + 2 , (10.9)
x y x y
u v
+ = 0. (10.10)
x y
Para la variacion de la presion con x se puede suponer que el flujo de la corriente principal
p
impone su distribucion de presiones dentro de la capa lmite, es decir, x dentro de la capa
p
lmite es igual al x existente fuera de la capa lmite, el cual se puede determinar mediante un
analisis de flujo potencial. Lo anterior indica que para una placa delgada, donde las lneas de
corriente son paralelas a la placa y por lo tanto no existe una variacion de la velocidad y por lo
p
tanto tampoco de la presion (ec. 4.2) el termino x = 0.
Las condiciones de borde, necesarias para resolver el sistema de ecuaciones, son u = v = 0 para
y = 0 y u U cuando y .
De un analisis dimensional se puede determinar que el numero de Reynolds (Rex ) es del orden
de magnitud de (1/ 2 ), es decir
1/2 1/2
x
=x
U Ux
1
Rex 2 . (10.11)
x
Blasius resolvio las ecuaciones simplificadas de Navier Stokes para el caso de la placa delgada
p
(donde x = 0) encontrando que la constante de proporcionalidad de la ecuacion anterior es
4.96
1
= 4.96 Rex 2 . (10.12)
x
C. Gherardelli U. de Chile
10.3 Capa lmite laminar 122
1
= 1.73 Rex 2 . (10.13)
x
1
= 0.664 Rex 2 . (10.14)
x
u
w =
y y=0
Se puede apreciar que el esfuerzo de corte decrece con x debido al aumento del espesor de la
capa lmite .
Como resolver las ecuaciones de Navier Stokes analticamente es practicamente imposible, salvo
en casos muy simplificados, se debe tener una forma alternativa de resolver el problema en forma
aproximada. Una forma de resolver este tipo de problemas es aplicando la ecuacion de cantidad
de movimiento sobre un elemento diferencial de capa lmite como se el que muestra la figura 10.9
y realizando aproximaciones sobre el perfil de velocidades tanto dentro como fuera de la capa
lmite.
y
d d+dd
dx
x
x x+dx
Las fuerzas externas segun x que actuan sobre el elemento diferencial de la figura 10.9 son
(fig. 10.10):
dp 1 dp
dFx = p p + dx ( + d) + p + dx d w dx . (10.16)
dx 2 dx
Reordenando se obtiene
dp
dFx = + w dx . (10.17)
dx
C. Gherardelli U. de Chile
10.3 Capa lmite laminar 123
(p+ 21 dx
dp
dx)dd
dd
pd
(p+ dx
dp
dx)(d+dd)
dx
twdx
Esta fuerza neta debe ser igual al flujo neto de cantidad de movimiento que sale del volumen de
control (fig. 10.11), es decir,
Z
dp d
+ w dx = u2 dy dx + U (V dx U d) . (10.18)
dx dx
0
V
U n
u n
n
dx
Z
d
Z Z
u dy V dx + U d u dy + u dy dx = 0, (10.19)
dx
0 0 0
La ecuacion 10.21 se denomina ec. integral de la cantidad de movimiento de von Karman. Para
el caso de la placa placa delgada, donde (dp/dx) = 0, se cumple que U = cte., es decir, U 6= U (x)
C. Gherardelli U. de Chile
10.3 Capa lmite laminar 124
por lo que se puede introducir dentro de la derivada y la integral en el segundo termino del lado
derecho de la ecuacion 10.21.
Z
d
w = (u2 U u)dy . (10.22)
dx
0
w d 2 dU d dU
= U + U = U2 + ( + ) U . (10.23)
dx dx dx dx
w d
= U2 . (10.24)
dx
Integrando w sobre toda la superficie se obtiene el arrastre que el fluido ejerce sobre el cuerpo.
En forma general, para calcular el espesor de la capa lmite sobre una placa plana que se
encuentra inmersa en un flujo laminar, incompresible y permanente se debe:
1. utilizar el gradiente de presiones (dp/dx) determinado para el flujo exterior a la capa lmite.
2. Definir o aceptar una forma razonable del perfil de velocidades dentro de la capa lmite.
Ejemplo
Determinar el espesor de la capa lmite que se genera sobre una placa plana delgada suponiendo
que el perfil de velocidades dentro de la capa lmite es cuadratico.
El perfil de velocidades debe tener la siguiente forma
u = y + y 2
2
y y
u = 2U U .
C. Gherardelli U. de Chile
10.4 Transicion a la turbulencia 125
1
= 5.48 Rex2 .
x
Este resultado es del orden de un 10% mayor al obtenido en forma analtica por Blasius.
Analogamente se obtiene para el espesor de desplazamiento la siguiente relacion
1
= 1.835 Rex2 ,
x
resultado que es del orden de un 6% superior al resultado de Blasius. Las diferencias exis-
tentes entre ambas soluciones son atribuibles a que la solucion de la ecuacion de von Karman
supone un perfil de velocidades de forma conocida y por lo tanto es solo una solucion aproxi-
mada. La figura 10.12 muestra distintos perfiles tpicos de velocidad utilizados para resolver la
ecuacion 10.22 as como la solucion obtenida por Blasius.
Lineal
Cbico
Sinusoidal
Parablico
Los resultados presentados hasta aqu son validos para la capa lmite laminar. Estos resulta-
dos concuerdan con mediciones experimentales hasta el punto donde la capa lmite se torna
turbulenta.
Dado que el numero de Reynolds, parametro principal que gobierna la transicion a la turbulencia,
es una funcion de la posicion relativa sobre el cuerpo, siempre es posible alcanzar el regimen
turbulento en la capa lmite, independiente de la velocidad del flujo libre U y del fluido. Para
el caso particular de la placa plana basta con que esta sea lo suficientemente larga.
La transicion a la turbulencia depende de varios factores como
C. Gherardelli U. de Chile
10.4 Transicion a la turbulencia 126
2. rugosidad de la superficie,
3. curvatura de la superficie,
5. vibracion de la placa.
Para una placa plana y delgada la transicion ocurre para un rango de la coordenada x tal que
2 105 Rex 3 106 . La figura 10.13 muestra el efecto del numero de Reynolds en el espesor de
la capa lmite sobre una placa lisa. Se ve que para Re = 3.2105 el espesor aumenta rapidamente
y se puede considerar este valor de Reynolds como crtico para la transicion a la turbulencia.
Como Re es funcion de x se puede determinar la posicion x a lo largo de la placa en que mas
probablemente aparecera la transicion a la turbulencia.
Figura 10.13: Variacion del espesor de la capa lmite con el numero de Reynolds local para el
flujo sobre una placa plana.
La figura 10.14 muestra el efecto de el grado de turbulencia de la corriente libre sobre la transicion
para un flujo sobre una placa plana. Las lneas de este grafico representan los lmites entre los
distintos tipos de flujo y por lo tanto la curva inferior representa la variacion del numero de
Reynolds crtico con el grado de turbulencia de la corriente libre U 0 /U .
La rugosidad de la superficie de la placa producira un adelantamiento del punto de transicion
si se le compara con una superficie lisa (pelotas de golf). La transicion a la turbulencia conlleva
tambien una variacion notable del perfil de velocidades dentro de la capa lmite como se puede
observar en la figura 10.15.
C. Gherardelli U. de Chile
10.5 Capa lmite turbulenta 127
Figura 10.14: Efecto del grado de turbulencia de la corriente libre sobre la transicion para el
flujo sobre una placa plana.
Transicin
Turbulento
Laminar
Figura 10.15: Variacion tpicas del perfil de velocidades para distintos tipos de regmenes.
La estructura de la capa lmite turbulenta es muy compleja, irregular y aleatoria. No existe, por
lo tanto, una solucion exacta para el flujo en esta zona por lo que se recurre a aproximaciones
C. Gherardelli U. de Chile
10.5 Capa lmite turbulenta 128
y validacion experimental.
Utilizaremos la ecuacion de von Karman deducida anteriormente, donde u representara la media
temporal de la velocidad, es decir, u = u.
Blasius encontro que para Re 107 y una superficies lisa se cumple
1/4
w = 0.0225U 2 . (10.25)
U
1/7
y
u=U . (10.26)
Reemplazando las ecuaciones anteriores en la ec. de von Karman y desarrollando para una placa
plana delgada se obtiene
1
= 0.376 Rex 5 . (10.27)
x
1
5 4
= 0.0463 x5 , (10.28)
U
U 2
w = 0.0288 1 . (10.29)
Rex5
Integrando sobre toda la placa es posible determinar la fuerza de arrastre debida a la friccion
Df
ZL
Df = bw dx , (10.30)
0
area A
z}|{
bL
Df = 0.036 U 2 1 . (10.31)
Rel5
1
Df = CDf U 2 A (10.32)
2
C. Gherardelli U. de Chile
10.6 Separacion / Efecto del gradiente de presion 129
Df
CDf = 1 2
. (10.33)
2 U A
15
CDf = 0.072 Rel . (10.34)
Turbulento
desarrollado
Turbulento
Transicin
Figura 10.16: Coeficiente de arrastre CD para una placa plana paralela al flujo (Recrit = 5 105 ).
Las ecuaciones que aproximan las distintas zonas del grafico son
1.328
CDf = Re0.5 flujo laminar
l
0.455 1700
CDf = (log Rel )2.58
Rel transicion con Rex,c = 5 105
Hasta ahora se ha analizado el flujo en la capa lmite sobre una placa plana delgada, donde
p
consideramos que no existe variacion de la presion en la direccion del flujo, es decir x = 0. Lo
C. Gherardelli U. de Chile
10.6 Separacion / Efecto del gradiente de presion 130
anterior debido a que se supone que el gradiente de presion en la capa lmite es igual al gradiente
de presion existente fuera de la capa lmite y que como la placa es delgada las lneas de corriente
son paralelas a la placa por lo que no existe variacion de la velocidad y por ende de la presion
fuera de la capa lmite.
Para cuerpos distintos a una placa plana delgada existe una variacion (magnitud y sentido)
de la velocidad de la corriente libre alrededor del cuerpo y por lo tanto una variacion de la
p
presion x 6= 0. Como el flujo tiene que abrirse para que pase el cuerpo en la parte anterior
de estos se producira un aumento de la velocidad y una disminucion de la presion. Pasado el
punto de maximo espesor del cuerpo el flujo se expande nuevamente disminuyendo la velocidad
y aumentando la presion. Una consecuencia importante de la existencia de un gradiente de
presion es lo que se conoce como separacion de la capa lmite.
La disminucion de la presion en la direccion del flujo se denomina gradiente de presion favorable.
El aumento de presion en la direccion del flujo se denomina gradiente de presion adverso ya que
tiende a frenar el flujo. Por otro lado, los esfuerzos de corte existentes dentro de la capa lmite
tienden tambien a frenar el flujo. Las partculas fluidas tienden a fluir a regiones del flujo
donde el gradiente de presion sea menos adverso, generando de esta manera la separacion de
la capa lmite. Detras de la zona de separacion se produce una reestructuracion del flujo que
comunmente se conoce como estela. Como efecto de la separacion se produce un aumento notable
en el arrastre y una baja en la sustentacion.
El punto de separacion se caracteriza por tener una velocidad y un gradiente de velocidad nulo
en la pared como muestra la figura 10.17.
Capa lmite
u=0
tw= u =0
y y=0
La ecuaciones de Navier Stokes para la capa lmite en la superficie del cuerpo, donde u = v = 0,
es
!
1 dp 2u
= . (10.35)
dx y 2 y=0
dp
Para el caso de un gradiente de presion favorable, es decir, dx < 0 se cumple que
!
2u
< 0. (10.36)
y 2 y=0
Al alejarse de la pared del cuerpo y a medida que y se cumple que u U por lo que
u
y disminuye con y. Lo anterior indica que no habra un cambio en la curvatura del perfil de
C. Gherardelli U. de Chile
10.6 Separacion / Efecto del gradiente de presion 131
dp
Para el caso de un gradiente de presion adverso, es decir, dx > 0 se cumple que
!
2u
> 0. (10.38)
y 2 y=0
2
Cerca del borde la la capa lmite sin embargo, se sigue cumpliendo que yu2 < 0. Esto indica
que existe un cambio de la curvatura (punto de inflexion) del perfil de velocidades entre y = 0
e y = . De lo anterior se concluye que para que exista la separacion debe existir un gradiente
de presion adverso. Esta es una condicion necesaria pero no suficiente. La separacion produce
ademas de un aumento en el arrastre una disminucion de la sustentacion pudiendo llegar a
situaciones extremas como la perdida total de la sustentacion (stall).
En la capa lmite turbulenta las partculas fluidas tienen mayor cantidad de movimiento que
en una capa lmite laminar (esto se puede ver en la diferencia entre los perfiles de velocidad
existentes en ambas situaciones) y por lo tanto pueden resistir mejor un gradiente de presion
adverso y, por lo tanto, resiste mejor la separacion. Analizado desde otro punto, si se adelanta
la transicion a la turbulencia de la capa lmite, se retardara la separacion de esta influyendo
positivamente en el arrastre (menor) y la sustentacion (mayor). Lo anterior puede ser visto en
la figura 10.18 donde, para una misma velocidad de la corriente libre, se ha producido mediante
un alambre en la parte anterior de la esfera la transicion adelantada a la turbulencia en forma
artificial (fig. 10.18(a)) produciendo un retraso en la separacion de la capa lmite.
Figura 10.18: Desprendimiento de la capa lmite sobre una esfera con y sin generadores de
vortices.
Este fenomeno se utiliza en la practica por ejemplo en pelotas de golf y en algunos perfiles alares.
C. Gherardelli U. de Chile