Capítulo 3

TEORÍA DE LA CAPA LIMITE

El movimiento de un fluido con relación a una superficie sólida o líquida determina la

formación de una zona cercana a la superficie en la que se dan gradientes fuertes de velocidad.
Esto ha sido comprobado experimentalmente e indica que en esa zona el esfuerzo cortante es
elevado, así como las consecuentes pérdidas de energía por fricción. Esa zona es la capa límite
hidrodinámica.

Figura 3.1. Capa límite hidrodinámica

Figura 3.2. Capa límite turbulenta y subcapa laminar

3.1. DESCRIPCIÓN DEL FENÓMENO

Si un fluido se aproxima laminarmente con una velocidad o a una superficie sólida, en

dirección paralela a la misma, la capa de fluido en contacto con la superficie se adhiere a ella.
(Fig. 3.1). En cualquier posición "x" a lo largo de la superficie de la placa habrá capas del fluido
colocadas a cierta distancia de la superficie, que aún "no se habrán enterado" de la existencia de
ésta y conservarán su velocidad original o. Entre ambas capas fluidas, la de velocidad nula
adherida a la superficie sólida y la de velocidad o se darán todos los valores intermedios de la
misma. Se define la capa límite hidrodinámica laminar como la zona comprendida entre x =0 y x
=0.99o . Así puede definirse el "borde" de la capa límite como el lugar geométrico de los puntos
en los cuales x =0.99o. Al aumentar el valor de "x" aumenta el espesor de la capa límite
hidrodinámica laminar hasta un punto en que el número de Reynolds local, definido como Re x
 3-2

=xo / tiene un valor de 1.5-2.0x10 5 . En este punto comienzan a aparecer turbulencias que se
establecen completamente al llegar el valor local del reynolds a 4-5x10 5. La zona laminar, no
obstante, nunca desaparece completamente sino que se reduce a una subcapa laminar,
adyacente a la superficie sólida y separada de la turbulenta por una zona de transición o "buffer"
(Fig. 3.2). La subcapa laminar es importante pues en ella se dan casi totalmente las caídas de
velocidad, temperatura y concentración. Fué Prandtl en 1904 quien estableció el concepto de
capa límite como una zona del fluido en movimiento, adyacente a la superficie sólida, en la cual
se manifiestan los efectos del esfuerzo cortante. El valor máximo de este último se localiza en
contacto con la superficie sólida. En cuanto a la presión, ésta no varía prácticamente con la
posición a lo largo de la placa plana y es igual a la que existe fuera de ella. La hipótesis de la
capa límite simplificó el tratamiento matemático del movimiento de los fluidos viscosos. Una
extensa obra sobre el tema es la ref(1).

Si un observador se moviera en la dirección "x" a una distancia fija de la placa(es decir a

un valor dado constante de "y") vería que la velocidad x disminuye, es decir, que el fluido va
"frenando" según puede apreciarse en la misma Fig. 3.1. Si la densidad es constante esto
implica que parte del fluido se desvía hacia arriba en la figura, generándose así un componente
y de la velocidad. Se trata pues de un flujo bidimensional en régimen permanente. El objetivo
aquí es obtener una expresión para la fuerza de arrastre. Esta expresión es importante pues de
ella se partirá para definir el factor de fricción y establecer posteriormente la analogía entre los
tres fenómenos de transporte. Una solución exacta del problema parte de la solución de Blasius a
la ecuación de movimiento. Otra solución aproximada es debida a Von Karman.

Figura 3.3. Solución de Blasius

3.2. SOLUCIÓN EXACTA

Las ecuaciones de continuidad y las de movimiento, obtenidas de las tablas del cap.1 son
respectivamente:

 X  y
 0 (3.1)
x y
 3-3

Que el esfuerzo cortante tenga su máximo valor en la superficie sólida se traduce en que
el valor máximo de la derivada   x /  y se ubique en y=0 mientras que el valor de esta derivada
es prácticamente cero en el “borde” de la capa límite. Este hecho implica que la primera derivada
  x /  y varíe grandemente en el pequeño espesor de la capa límite. En cambio en la dirección
"x" puede decirse que la variación de la velocidad x es pequeña y casi lineal, de manera que
 2  x /  x 2 tiene un valor cercano a cero y por consiguiente despreciable en relación con
 2  x /  y 2 pudiendo eliminarse en la última ecuación quedando:

(3.2)

Las condiciones límite para esta ecuación son:

x (x=0, y)= o
x (x ,y=0)=0
x (x ,y=)=o

Las Ecs. (3.1, 3.2) pueden transformarse en ordinarias por medio de un cambio de
variable apropiado. Blasius introdujo una función “de corriente” que sustituye a las dos variables
espaciales llegando mediante un proceso matemático a una ecuación diferencial ordinaria no-
lineal que lleva su nombre:

d 2 d 3
  2 0 (3.3)
d 2 d 3

siendo  es la velocidad adimensional definida por

 = x /o (3.4)

La variable  es la variable espacial adimensional, definida como

y
 Re x (3.5)
x

La Ec. (3.3) fue resuelta primeramente por Blasius en forma numérica obteniendo para 
una solución asintótica con =1. En la Tabla 3.1 aparecen algunos valores numéricos de las
variables , . La pendiente de la curva en el origen es igual a 0.332 según se ve en la tabla.(Fig.
3.3)

La validez de la solución de Blasius fue comprobada experimentalmente por Nikuradze en

1942. No existiendo contribución de la presión a la fuerza de arrastre, ésta última es de origen
netamente viscoso y por lo tanto puede ser calculada partiendo de los componentes del tensor
esfuerzo:
 3-4

  
 yx y 0   x  y0 (3.6)
 y 
La derivada en la Ec. (3.6) puede ser obtenida más fácilmente a partir de la solución de
Blasius, por medio de la regla de la cadena

(3.7)
y por lo tanto,

Re x
 yx y 0  0.332  0 (3.8)
x

El signo(-) en la última ecuación obedece al hecho de que la transferencia de momentum

se da de zonas de mayor velocidad a otras de menor velocidad y por consiguiente en sentido
contrario al positivo de "y"; puede ser eliminado para tener una ecuación general. El valor
obtenido de yx es local, al igual que el del número de Reynolds. Un valor promedio  para toda la
longitud de la placa L es de mayor utilidad:

Tabla [Link]ón de Blasius

= (y/x)Rex  =x/o

0 0
0.4 0.1328
0.8 0.2647
1.6 0.5168
2.4 0.729
3.2 0.8761
4.0 0.9555
4.8 0.9878
5.6 0.9962
6.4 0.9996
7.2 1.

de donde:

(3.9)
 3-5

siendo ReL= Lo /. Para la fuerza de arrastre la ecuación resultante es:

FARR=0.664Wo (3.10)

De acuerdo con la Ec. (3.9) el valor promedio del esfuerzo cortante es igual al doble del
valor local en x=L:

= 2 L (3.11)

Definición del factor de fricción

Cuando el flujo es turbulento o la geometría del sistema es complicada no es posible
llegar a ecuaciones teóricas como las que se han obtenido aquí para placa plana pues los
mecanismos de flujo no son conocidos. Es así que en la mayor parte de los casos prácticos no es
posible obtenerlas. Para estos casos es necesario recurrir a un factor de fricción que para placa
plana se define mediante la ecuación:

= fP(K) = fp(o2/2-s2/2) (3.12)

K es la concentración de energía cinética(energía cinética por unidad de volumen =2/2.

En la superficie sólida la velocidad s es nula y por lo tanto la energía cinética también y por lo
tanto =fpo2/2. fp es el factor de fricción en placa, o sea el coeficiente de transferencia de
momentum. En flujo turbulento es empírico. En flujo laminar puede obtenerse combinando las
ecs(3.9, 3.12):

(3.13)

donde fp es valor medio para la toda la placa.

Espesor de la capa límite hidrodinámica

De la solución numérica de la ecuación de Blasius se obtiene que para =5, = 0.99.
Esto indica, de acuerdo con la definición de capa límite, que en el "borde" de la misma se cumple
 = 5. Denotando el espesor de la capa límite hidrodinámica como , se tiene entonces que

 H Re x
 5 (3.14)
x

de donde

(3.15)

La Ec. (3.15) da el espesor de la capa límite hidrodinámica en función de la posición a lo

largo de la placa.

Ejemplo 3.1. Fluye aire a presión de 1 atmósfera, T=38ºC, a lo largo de una placa plana
de 1.2 m de largo por 0.6 m de ancho. La velocidad de aproximación es de 1.5 m/s. Calcular el
 3-6

espesor máximo de la capa límite, la fuerza de arrastre y el factor de fricción. Calcular la longitud
que debería tener la placa para que la capa límite se hiciera turbulenta? =1.91x10 -5 Pas.

Solución:
=29/(0.082x311)=1.137 kg/m3.
ReL= 107,152 (flujo laminar)

a) =0.018 m
b) FARR =0.664x0.6x1.91x10-5x1.5 =3.74x10 -3 N

c) =4.057x10 -3

d)L= = =4.48 m

3.3. METODO APROXIMADO

La solución aproximada del problema de la capa límite es debida a Von Karman y no

toma en cuenta la componente del movimiento perpendicular a la superficie de la placa, y. Un
balance de masa en el elemento de volumen de la Fig. 3.4 es:

(3.16)

Figura 3.4. Balance de masa en un elemento de volumen

W es el ancho de la capa límite. o es el gasto másico que entra al elemento de volumen

con velocidad o. El espesor  depende del valor de x. El balance de momentum en el mismo
elemento de volumen es:
(3.17)
 3-7

s= -(dx/dt)y=0 es el esfuerzo cortante en la superficie de la placa. Despejando o de la

Ec. (3.16) e introduciendo en la Ec. (3.17) se obtiene

(3.18)

Dividiendo entre x y haciendo tender x a cero se obtiene:

(3.19)

que también puede escribirse:

(3.20)

La ecuación anterior es llamada ecuación integral de movimiento. Es válida para cualquier

fluido, en flujo laminar o turbulento. Conociendo el perfil x=f(y) puede integrarse la Ec. (3.20).
Von Karman postuló un polinomio de tercer grado:

(3.21)

Para evaluar las constantes se requiere cuatro condiciones frontera; tres son evidentes:
CF1. en y=0, x=0; CF2.- en y=, x=o; CF3.- en y=, x/y=0. La cuarta condición es: CF4.-en
y=0, 2x/y2 =0; se obtiene imponiendo a la Ec. (3.2) las restricciones x=0, y=0, válidas en y=0.

Por la CF1.- se obtiene a=0. Por CF4, c=0. Imponiendo las otras condiciones se obtiene:
b=3/2; d=-1/2. El perfil de velocidad queda:

(3.22)

Con la Ec. (3.22) es posible resolver la integral de la Ec. (3.20). Introduciéndola en ella se
llega a:

=0.1393 (3.23)

Por otro lado, s=-(x/y)y=0. Introduciendo x de la Ec. (3.22) se obtiene s=-3o/2.

Introduciendo ahora esta última ecuación y la Ec. (3.23) en (3.20) queda:

(3.24)
Integrando la ecuación anterior se obtiene: 2/2=10.768x/o+C1 y puesto que en x=0,
=0, C1=0 la ecuación anterior queda:
 3-8

(3.25)

Comparando la Ec. (3.25) con la (3.15) se observa que el resultado obtenido mediante el
método integral tiene -7% de error.

Si en la expresión s=-3o/2 se sustituye  por su expresión en Ec. (3.25) se obtiene: s=

y puesto que por definición fp=s/(o2), combinando estas dos últimas
ecuaciones, eliminando s, se llega a: . Por último, el valor medio para toda la
placa es dado por fp= . El resultado es:

(3.26)

Este resultado es comparable con el obtenido por el método exacto, con error <3%. El
método aproximado es útil para tratar matemáticamente la capa límite de fluidos no-newtonianos.

Otros perfiles de velocidad han sido evaluados. Así por ejemplo si se postula un perfil de segundo

grado, ; se requieren tres condiciones frontera para evaluar los coeficientes. Usando

las tres primeras ya empleadas se obtiene: a=0; b=2; c=-1. Siguiendo el mismo procedimiento se obtiene:
=5.47x/ , fp=1.464/ . De todas estas soluciones se desprende que  incrementa con y
disminuye con

3.4. LIQUIDOS NO NEWTONIANOS

Dos casos diferentes constituyen los líquidos que siguen la ley de potencia y los plásticos
de Bingham. En el primer caso la capa límite está limitada por el plano en el cual x/o=0.99 al
igual que en los newtonianos. Más allá de este plano no existe prácticamente gradiente de
velocidad. En el caso de los plásticos Bingham, en cambio, en el plano que delimita la capa límite
el valor del esfuerzo cortante es igual casi siempre al de cedencia y por lo tanto más allá de esa
superficie no existirá gradiente alguno de velocidad. Es esta condición yx=o la que establece la
separación entre la capa límite y el resto del fluido. En caso contrario el esfuerzo de cedencia se
alcanzará fuera de la capa límite y dentro de ella el comportamiento será newtoniano. El
procedimiento ya empleado en el método integral puede ser aplicado a los líquidos no-
newtonianos.

Líquidos con modelo de potencia

En la Ec. (3.20) no hay implícita consideración alguna en relación con la naturaleza del
fluido, de manera que es aplicable para no-newtonianos. El perfil de velocidad de tercer grado ya
obtenido para el caso newtoniano es igualmente el mismo, pues las condiciones frontera son las
mismas. No obstante, la tercera condición no es exacta para suspensiones pues no se tiene en
ese caso un medio continuo. Considerando válida la Ec. (3.22) la expresión para el esfuerzo
cortante en la superficie sólida es: s=-K(/y)ny=0=-K(3o/2)n. Introduciendo esta última
ecuación, y la Ec. (3.23) en la (3.20) queda:
 3-9

(3.27)

que integrada queda:

(3.28)

C1=0 pues =0 en x=0. De la Ec. (3.28) se obtiene la siguiente expresión para el espesor
:

(3.29)

que se puede escribir

(3.30)

siendo f(n)= , Rex=o2-nxn/K

Introduciendo en la ecuación =yxy=0=-K(3o/2)n el valor de  dado por la Ec. (3.30) se

obtiene:

(3.31)

y para el valor medio,

(3.32)

siendo ReL=o2-nLn/K. Por último, fp=s/(½o2). De la Ec. (3.32):

(3.33)

Ejemplo 3.2. Calcular la fuerza de arrastre que ejerce un líquido sobre una placa plana
sobre la cual fluye tangencialmente con velocidad de aproximación o=20 m/s. El modelo
reológico del líquido es =0.8(-d/dx)0.72; su densidad es '=1.05, las medidas de la placa son
L=1.2 m, W=0.6 m

Solución:
 3-10

f(n)=

ReL=

Por la Ec. (3.32): Pa

FARR=458.5x1.2x0.6=330 N

No hay datos experimentales que permitan establecer el valor del Reynolds en que la
capa se hace turbulenta. Según una ecuación propuesta por Skelland ( 4), el valor del Reynolds
crítico Re'c=o2-nxcn/K es menor para los líquidos con ley de potencia que para los newtonianos:
para n=0.5 el Re'c es del orden de 12,000

Plásticos de Bingham
En el balance de momentum y fuerzas (Ec. 3.17) es necesario considerar el esfuerzo
cortante o que actúa en la parte superior del elemento de volumen de la Fig. 3.4. La ecuación
integral de momentum queda:

(3.34)

Asumiendo el mismo perfil de velocidad de la Ec. (3.22) queda para este caso:
(3.35)
s-o=-(dx/dy)y=0=-3oo/2. Integrando la Ec. (3.35) y reordenándola se llega a la
misma ecuación obtenida para newtonianos(Ec. 3.25) siendo ahora Re x=xo/o. Por el mismo
procedimiento seguido para los newtonianos se obtiene(3):

(3.36)

(3.37)

El grupo adimensional 2o/o2 es el número de Bingham.

Ejercicios (ver respuestas en el Apéndice D)

1. Si se definiera que en el borde de la capa límite la velocidad adimensional vale 0.9¿cuál sería la
ecuación para el espesor ?
 3-11

2. Obtener las ecuaciones para el espesor de la capa límite =(x) y el coeficiente de fricción
mediante el método integral. Postular perfil de velocidad de primer grado: x/o=a+b(y/); usar las
condiciones 1,2 ya empleadas.
3. A lo largo de una placa plana de 1.2 m de largo pasa agua con velocidad de aproximación de 0.5 m/s. Se
quiere saber a qué distancia "y" de la superficie de la placa la velocidad del agua es de 0.3645 m/s en una posición
ubicada a x=40 cm del borde de ataque.
4.-Se ha propuesto para la capa límite turbulenta la ecuación f x=0.0576Rex-0.2 ¿Cuál es la ecuación para el
valor medio de fp en la zona turbulenta de longitud L?
5. Considerando que el Reynolds para el cual la capa límite se hace turbulenta es Re c=500,000,
obtener una ecuación para el factor de fricción en una placa plana de longitud L suficientemente grande para
que se alcance el flujo turbulento. Para capa límite turbulenta en placa plana Blasius postuló la ecuación
fx=0.0576(Rex)-0.2

Bibliografía

1. Schlichting,H.,"Boundary-Layer Theory", McGraw-Hill, Nueva York .1955

2. Blasius,H.,[Link].1(1908)
3. Chabra, R.P.,Richardson,J.F. "Non-Newtonian Flow in the Process Industries".Butterworth-Heinemann. 1999
4. Skeland, H.P., Diffusional Mass Transfer. John Wiley and Sons.1974

Notación
A área, m2
fp factor de fricción en placa plana, adimensional (valor medio)
fx ídem, valor local
ReX,ReL reynolds local, y evaluado en x=L, adimensional
W ancho considerado de la capa límite, m

 espesor de la capa límite hidrodinámica, m

 variable espacial adimensional,
 viscosidad absoluta, kgm-1s-1
 viscosidad cinemática, m2s-1
 densidad, kgm-3
yx densidad local de flujo de cantidad de movimiento, Pa
L valor de tyx, evaluado en x=L, Pa
 valor medio de tyx, evaluado en y=0, Pa
o velocidad de aproximación, ms -1
 función velocidad normalizada, adimensional

Subíndices:
x valor local
L valor en x=L