Aproximación de capa límite
Como ya se ha indicado anteriormente, la capa límite es una zona
delgada junto a una pared en la que están confinados los efectos
viscosos en flujos con elevados números de Reynolds. En ella las
fuerzas de inercia y las debidas a la viscosidad son del mismo orden.
Cuando el número de Reynolds basado en el espesor de la capa límite y
en la velocidad exterior supera un cierto valor, la capa límite se hace
inestable y se produce la transición a la turbulencia. El número de
Reynolds de transición depende del gradiente de presión reducida (la
capa límite es más estable si el gradiente de presión reducida es
positivo y elevado); del nivel de turbulencia de la corriente exterior (un
nivel alto favorece la transición); de la rugosidad de la pared (una
rugosidad elevada favorece la transición), y de otros factores tales
como la succión de la capa límite, fenómenos de transmisión de calor,
etc. En una placa plana lisa, con un nivel de turbulencia en el flujo
exterior del 1%, la transición se produce para un valor de
Re ( es la velocidad de la corriente exterior y el
espesor de desplazamiento). La capa límite turbulenta se caracteriza
por ser más gruesa y dar lugar a tensiones cortantes mayores que la
capa límite laminar.
Cuando la capa límite permanece adherida a la pared, es posible tratar
el flujo utilizando ecuaciones apropiadas para describir el movimiento
en la capa límite y las ecuaciones de Euler en el exterior de ésta.
En esta sección se presentarán las ecuaciones simplificadas que
permiten describir el flujo en la capa límite, así como otros flujos en los
que se satisfagan las condiciones de aplicabilidad que se indicarán en lo
que sigue. Las ecuaciones de conservación de partida serán las
correspondientes a un flujo bidimensional, estacionario y laminar, de un
fluido incompresible con viscosidad constante (se supondrá que las
fuerzas másicas derivan de un potencial, y se utilizará la presión
reducida, en lugar de la presión, como variable dependiente):
(5.11
5)
(5.11
6)
(5.11
7)
La coordenada es perpendicular a la pared. Las ecuaciones para la
capa límite turbulenta pueden deducirse a partir de las ecuaciones
(5.96) y (5.97) utilizando consideraciones análogas a las que se
presentan a continuación para la capa límite laminar.
En la capa límite la velocidad varía desde un valor nulo en la pared (en
una pared en reposo) hasta un cierto valor característico de un flujo
ideal en el exterior de la capa límite. Si el espesor de la capa límite (a lo
largo del cual se produce dicha variación de velocidad) se supone muy
pequeño frente a la longitud característica paralela a la pared en la que
la velocidad del fluido varía apreciablemente, mediante estimación de
órdenes de magnitud en las ecuaciones se deduce fácilmente que las
ecuaciones (5.115) a (5.117) pueden simplificarse y adoptar la siguiente
forma aproximada (véase, por ejemplo, Liñán, 1967; Crespo,
1997; Batchelor, 1967):
(5.11
8)
(5.11
9)
(5.12
0)
Obsérvese que en la ecuación (5.119) se ha despreciado el término de
difusión en la dirección del flujo y que la ecuación de conservación de la
cantidad de movimiento según ha quedado reducida a la condición
de la que las variaciones transversales de presión reducida son nulas.
La presión es la presión reducida en el exterior de la capa límite,
que en principio puede determinarse a partir de la solución del
problema de flujo ideal en el exterior de la capa límite. La ecuación de
Bernoulli en el flujo exterior es
(5.12
1)
constante
siendo la componente de la velocidad en el exterior de la capa
límite, o la velocidad de deslizamiento de un fluido ideal sobre la pared.
Derivando la ecuación (5.121), resulta
(5.12
2)
Como consecuencia de la eliminación del término , el carácter
mixto elíptico/hiperbólico de las ecuaciones (5.115) a (5.117) pasa a ser
parabólico/hiperbólico en las ecuaciones (5.118) a (5.120) de la capa
límite, jugando en éstas la coordenada el papel que desempeña el
tiempo en problemas no estacionarios, lo que permite obtener la
solución siguiendo un procedimiento de avance según la coordenada ,
tal como se explicará en la Sección 6.1. Por tanto, será necesario
imponer una condición inicial
(5.12
3)
y condiciones de contorno,
(5.12
4)
siendo