1.
Fundamento y solución de la capa límite de VON KARMAN
Del perfil de la velocidad en la capa límite de la figura, se observa que la velocidad pasa de
u=0.99U ∞ en y=δ a u=0 en y=0 a lo largo de una distancia muy corta (el espesor de la capa
límite).
Figura 1. Capa límite en aire con ℜcrít =3 ×105 .
Fuente: [CITATION Pot15 \p 387 \n \l 10250 ]
Si el perfil de la velocidad puede considerarse conocido, las ecuaciones integrales de continuidad
y de la cantidad de movimiento harán posible predecir el espesor de la capa límite y el cortante
en la pared y, por tanto, la resistencia al avance.
Figura 2. Volumen de control.
Fuente: [CITATION Pot15 \p 388 \n \l
10250 ]
Considerando un volumen de control infinitesimal,
mostrado en la figura. La ecuación integral de continuidad nos permite hallar:
δ
∂
ṁ parte superior =ṁsalida −ṁentrada =
∂x (∫ )
0
ρu dy dx (1)
Figura 3. Flujo másico.
Fuente: [CITATION Pot15 \p 388 \n \l 10250 ]
La ecuación integral de la cantidad de movimiento toma la forma:
∑ F x =mom
˙ salida −mom
˙ entrada −mom
˙ parte superior (2)
donde mom
˙ representa el flujo de la cantidad de movimiento en la dirección x. Teniendo en cuenta las
siguientes figuras:
Figura 4. Fuerzas.
Fuente:[CITATION Pot15 \p 388 \n \l 10250 ]
Figura 5. Flujo de la cantidad de movimiento.
Fuente: [CITATION Pot15 \p 388 \n \l 10250 ]
Y despreciando los términos de orden superior, esto se convierte en:
δ δ
−δdp−τ 0 dx=
∂
∂x ( 0
) ( (∫ ) )
∫ ρ u2 dy dx− ∂
∂x 0
ρu dy dx U ( x )
Si se divide todo entre −dx se obtiene:
δ δ
dp d
τ 0+ δ
dx
=U ( x )∫ ρu dy−
0
dx 0
(
∫ ρu 2 dy ) (3)
Donde se utilizan derivadas ordinarias porque las integrales solo son funciones de x. Esta
ecuación se conoce como "ecuación integral de Von Karman .
Solución aproximada de la capa límite laminar
Es posible usar la ecuación integral de Von Karman y obtener una aproximación bastante precisa
de la capa limite laminar en una placa plana con gradiente de presión cero. Tenemos cuatro
condiciones que un perfil de velocidad propuesto debe satisfacer:
u=0 en y =0
u=U ∞ en y =δ
∂u
=0 en y =δ
∂y
∂2 u
=0 en y =0
∂ y2
Las tres primeras condiciones son obvias si se observa un dibujo de perfil de velocidad, mientras
que la cuarta condición proviene de la ecuación de Navier Stokes para la componente x, ya que
u=v=0 en la pared y dp /dx=0 para el flujo permanente sobre la placa plana considerada.
Supongamos que el perfil de la velocidad se aproximado a uno de tipo parabólico
Si se aproxima el perfil de velocidades a uno de tipo cuadrático:
u
=A + By+C y 2
U
Y se aplican las condiciones de contorno anteriormente descritas, se genera el siguiente sistema:
0=A
1= A+ Bδ +C δ 2
0=B+ 2C δ
Una solución simultanea da
A=0
2
B=
δ
−1
C=
δ2
El perfil de la velocidad es entonces:
u y y2
=2 − 2
U∞ δ δ
Esto se sustituye en la ecuación integral de Von Karman para obtener
δ δ
dp d
τ 0+ δ =U ( x )∫ ρu dy−
dx 0
dx (
∫ ρu 2 dy
0
)
δ
d y y2 2 y y2
τ 0= ρU
2
(
∫ ∞ δ δ2
dx 0
2 − )(
1− −
δ δ2
dy )
2 dδ
τ 0= ρ U ∞2
15 dx
Se sabe que en la pared τ 0=μ ∂u /∂ y ¿ y=0, y si se utiliza el perfil cuadrático:
2
τ 0=ρ U ∞
δ
Igualando las dos expresiones anteriores, obtenemos
v
δ dδ=15 dx
U∞
Usando δ =0 en x=0, esto se integra a:
vx x
δ =5,48
U∞
=5,48
√
√ Rnx
(4)
Donde Rnx es el número de Reynolds local. Se sustituye esto de nuevo en la ecuación (5.7) para
obtener el esfuerzo cortante en la pared, sabiendo que v=u/ ρ:
2 μ U∞
τ 0=
δ
v
τ 0=0,365 ρ U ∞2
√ xU ∞
0,365 ρU ∞2
τ 0= (5)
√ Rn x
Si se hace adimensional este esfuerzo cortante dividiéndolo entre 1/2 ρU 2, es resultado es el
coeficiente de fricción superficial C f :
τ0
Cf =
1
ρ U2
2
730 730
Cf = =
x
√ Rn x (6)
√U v
Si se integra el esfuerzo cortante en la pared a lo largo de la longitud L, se obtiene, por unidad de
anchura, la expresión siguiente para la fuerza de arrastre:
L
D F=∫ τ 0 ( x ) dx
0
D F=0,730 ρU √ ULv
0,730 ρ U 2 L
D F=
√ RnL
O, en términos del coeficiente de fricción superficial C f :
DF
Cf =
1
ρ U2 L
2
1,460 1,460
Cf = =
UL √ Rn L
√ v
Supongamos que el perfil de la velocidad se aproximado a uno de tipo cúbico
Si se aproxima el perfil de velocidades a uno de tipo cuadrático:
u
= A+ By+C y 2+ D y 3
U∞
Y se aplican las condiciones de contorno anteriormente descritas, se genera el siguiente sistema:
0=A
1= A+ Bδ +C δ 2
0=B+ 2C δ
Utilizando las cuatro condiciones de contorno, quedan los valores:
A=0
3
B=
2δ
C=0
−1
D=
2 δ3
El perfil de la velocidad es entonces:
3
u 3y 1 y
= −
U∞ 2 δ 2 δ ( )
Esto se sustituye en la ecuación integral de Von Karman para obtener
δ δ
dp d
τ 0+ δ =U ( x )∫ ρu dy−
dx 0
dx (
∫ ρu 2 dy
0
)
δ
d 3 y 1 y3 3 y 1 y3
τ 0= ρ (
∫ 2 δ 2 δ3
dx 0
− 1− − )(
2δ 2 δ3 ∞
2
U dy )
dδ
τ 0=0,139 ρ U ∞2
dx
Se sabe que en la pared τ 0=μ ∂u /∂ y ¿ y=0, y si se utiliza el perfil cúbico:
3
(
τ 0=μ U ∞
2δ )
Igualando las dos expresiones anteriores para τ 0 ( x ), obtenemos
3
μU
2 v
δ dδ= dx=10,8 dx
0,139 ρ U 2
U
Usando δ =0 en x=0 (el borde de entrada), se integra la ecuación anterior para obtener
vx x
δ =4,65
U∞
=4,65
√√ Rnx
Donde Rnx es el número de Reynolds local. Se sustituye esto de nuevo en la ecuación para
(7)
obtener el esfuerzo cortante en la pared, sabiendo que v=u/ ρ:
2 μ U∞
τ 0=
δ
v
τ 0=0,323 ρ U ∞2
√ xU ∞
0,323 ρU ∞2
τ 0= (8)
√ Rn x
El esfuerzo cortante se hace adimensional al dividir entre 1/2 ρU 2. El coeficiente de fricción
superficial local C f resultante es:
τ0
Cf =
1
ρ U ∞2
2
0,646 0,646
Cf = =
x √ Rn x (9)
√U∞
v
Si se integra el esfuerzo cortante en la pared a lo largo de la longitud L, se obtiene, por unidad de
anchura, la expresión siguiente para la fuerza de arrastre:
L
D F=∫ τ 0 ( x ) dx
0
D F=0,646 ρU ∞ √ U ∞ Lv
0,646 ρ U 2 L
D F=
√ RnL
O, en términos del coeficiente de fricción superficial C f :
DF
Cf =
1
ρ U2 L
2
1,29 1,29
Cf = =
√ U ∞ L/v √ Rn L
