Capa Límite

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 9

UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

(Universidad del Perú, Decana de América)


FACULTAD DE CIENCIAS FISICAS
ESCUELA ACADEMICO PROFESIONAL DE INGENIERIA MECANICA DE FLUIDOS

PRACTICA CALIFICADA TEMA 2

CAPA LIMITE

JULIO 2017

PROFESOR DEL CURSO: Ing. Alvarado Torres Wilson Emilio


ESTUDIANTE: Marca Andía Anthony Michael
CÓDIGO: 14130158
PROBLEMA

1. Una placa plana se sumerge en un fluido de agua. En el fluido se ubica una sonda en miniatura
que mide la presión de estancamiento, esta sonda se puede desplazar por el eje normal a la
pared, en este mismo lugar en la pared se tiene ubicada una apertura que permite obtener el
valor de la presión estática. La sonda y el manómetro están unidas por un manómetro
diferencial de agua que muestra la diferencia entre la presión de estancamiento y la presión
estática de la pared. El resultado de la medición es:

Y (mm) 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Po-P (Pa) 1039 1356 1573 1743 1885 2008 2117 2215 2305 2388 2433 2533 2606 2700
Para y>14 mm la diferencia Po-P no cambia

La distribución de velocidades no corresponde a una función exponencial:

a) ¿Calcular es espesor de la capa límite?


b) ¿Calcular los espesores convencionales?
c) ¿Encontrar la distribución de velocidad?
d) ¿Calcular la fuerza de rozamiento sobre ambos lados de la placa?

Si 𝓁=1.2 m y tiene un ancho b=1.0 m


SOLUCIÓN:

Consideraciones:
 El flujo de la Capa Límite permanece constante 𝑈 = 𝑈∞
 Para valores superiores a y=14 mm la diferencia de presión 𝑃0 − 𝑃 no varía
 Para determinar si la C. L. es laminar o turbulento se debe efectuar el siguiente
algoritmo: 𝑢 = 𝑓(𝑥 )

 Dentro de la Capa limite se cumple:

𝑢2 2(𝑃0 − 𝑃)
𝑃0 = 𝑃 + 𝜌 → 𝑢=√ … (1)
2 𝜌
 Fuera de la C. L.

𝑈2 2(𝑃0 − 𝑃´)
𝑃0 = 𝑃´ + 𝜌 → 𝑈=√ … (2)
2 𝜌

Luego, para y>14 mm al dividirse (1) y (2) se tiene:


𝑢 𝑃0 − 𝑃
=√ ; 𝑐𝑜𝑚𝑜 𝑃0 − 𝑃´ = 2700 𝑃𝑎
𝑈 𝑃0 − 𝑃´
𝑢
Luego tomamos en la siguiente tabla de cada punto:
𝑈
Y/δ 0.071 0.143 0.214 0.286 0.357 0.429 0.500 0.571 0.643 0.714 0.786 0.857 0.929 1

u/U 0.620 0.709 0.763 0.803 0.836 0.862 0.885 0.906 0.924 0.940 0.949 0.969 0.982 1

Y/δ (u/U) vs (Y/δ)


1.200

1.000

0.800

0.600

0.400

0.200

0.000 u/U
0.000 0.200 0.400 0.600 0.800 1.000 1.200

𝑢 y
Por la forma de la gráfica: = 𝑓( )
𝑈 δ
Se observa que corresponde a una distribución potencial:
𝑢 y
= ( )𝑛
𝑈 δ
𝑢 𝑦
Aplicando logaritmo: 𝐿𝑜𝑔( ) = 𝑛
𝑈 δ

Log(Y/14) -1.146 - -0.669 -0.544 -0.447 -0.368 -0.301 -0.243 -0.192 -0.1461 -0.1047 -0.0669 -0.032 0
0.845
Log(u/U) -0.207 -0.149 -0.117 -0.095 -0.078 -0.064 -0.053 -0.043 -0.034 -0.0267 -0.0226 -0.0139 -0.008 0

Podemos representar en la gráfica y=Kx

(Y/δ) vs (u/U)
0
-1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0

-0.05

-0.1

-0.15

-0.2

-0.25

Encontramos el valor de “n”, para una C. L. Turbulenta


K=n donde:
𝑙𝑜𝑔(0.982) − 𝑙𝑜𝑔(0.620)
𝑇𝑎𝑛(𝐾) = = 0.1788
log(0.924) − log(0.0071)

Luego tomando: arctan(𝐾) = 0.177 rad = n


Empezamos a hallar los espesores para n= 0.177.
0.177
𝑛 0.177
𝛿 ∗ = 𝛿 ∫ (1 − 𝑛𝑛 )𝑑𝑛 = 14 × [ ]
𝑛+1 0
0

𝛿 ∗ = 2.105 𝑚𝑚

𝑛
𝛿 ∗∗ = ∗𝛿
(𝑛 + 1)(2𝑛 + 1)
𝛿 ∗∗ = 1.55 𝑚𝑚

0.177
𝑢 𝑢
𝛿 ∗∗∗ = ∫ (1 − ) 𝑑𝑦 … (∗)
𝑈 𝑈
0
𝑢 y
Como: = ( )0.177 , reemplazamos (*)
𝑈 δ
14
∗∗∗
𝑦 0.177 𝑦
𝛿 = ∫( ) (1 − ( )2×0.177 ) 𝑑𝑦
14 14
0
∗∗∗
𝛿 = 2.75 𝑚𝑚

Además:
2(𝑃0 −𝑃´) 𝑘𝑔
𝑈=√ 𝜌
; 𝜌𝑎𝑔𝑢𝑎 = 998.2 𝑚2 𝑎 𝑇 = 20 0 𝐶 , 𝑃0 − 𝑃 = 2700

2(2700´)
𝑈=√ → 𝑈 = 2.32 𝑚/𝑠
998.2

𝑚2
Sabemos que 𝜈𝑎𝑔𝑢𝑎 = 1.007 ∗ 10−6 𝑠
𝑎 𝑇 = 20 0 𝐶 𝑦 𝑈 = 𝑈∞ fuera de la C. L.

𝜈𝑎𝑔𝑢𝑎 0.25
𝜏𝜔 = 0.0232 ∗ 𝜌 ∗ 𝑈∞ 2 ( )
𝑈∞ ∗ 𝛿

−6 0.25
1.007 ∗ 10
𝜏𝜔 = 0.0232 ∗ 998.2 ∗ 2.322 ( )
2.32 ∗ 14 ∗ 10−3

𝜁𝜔 = 9.301 𝑁/𝑚2
PROBLEMA

2. Un tanque cilíndrico con radio de 8 m y altura 8 m que se muestra en la Figura 1, el viento que
impacta contra el cilindro tiene una velocidad de 50 m/s el cual corre en una zona plana de
500 m sin obstáculos, se tuene un P=101.3 KPa, ν =1.57*10-5 m2/s y T=20 0C.
Encontrar mediante el método de Von Karman los parámetros de la capa límite.
𝑈(𝑥) = 2𝑈∞ ∗ 𝑆𝑒𝑛(𝜑)

SOLUCIÓN:

Tenemos los siguientes datos:


R=8 m
V=50 m/s
P=101.3 KPa
ν=1.57*10-5 m2/s
T=20 0C

A partir de la ecuación de Von Karman, se


obtiene la ecuación de Twartez, para el
flujo sobre un cuerpo cilíndrico.
𝜇𝑢𝜉(𝑓)
𝜏𝜔 =
𝛿 ∗∗
Donde:
𝜏𝜔 : Tensión de Rozamiento
Μ: Viscosidad dinámica del fluido
𝛿 ∗∗ : Espesor de la C. L. de pérdida de impulso

Como: 𝜉(𝑓) = (𝑓 + 0.09)0.62


∗∗ 2
𝛿 𝑑𝑢
𝑓= Conocido como gradiente de presiones adimensional … (α)
υ 𝑑𝑥
2 0.45𝜈 𝑥 5
𝛿 ∗∗ 2 = 𝛿𝑜 ∗∗ + ∫ 𝑢 𝑑𝑥 … (𝛽)
𝑢6 0

Para simplificar:
Elegimos una distribución de velocidad en un sistema de coordenadas polares
donde la velocidad tiene componentes Ur y Uθ.

 De acuerdo al flujo potencial sobre un


cilindros e tiene:
𝑅2
𝑈𝜃 = −𝑈∞ (1 + 𝑟1 ) 𝑆𝑒𝑛𝜃

Dónde: r = R y θ = 180 - ϕ
𝑈𝜃 = −𝑈∞ 𝑆𝑒𝑛𝜃
La velocidad que rodea al cuerpo cilíndrico será:

𝑈(𝑥) = −𝑈𝜃 = 2𝑈∞ 𝑆𝑒𝑛𝜃 ← Se observa que 𝑈(𝑥) Crece en sentido opuesto a 𝑈𝜃

𝑈(𝑥) = 2𝑈∞ 𝑆𝑒𝑛(𝛱 − 𝜑)


𝑈(𝑥) = 𝑈𝜃 = 2𝑈∞ 𝑆𝑒𝑛(𝜑)

Ahora evaluamos la integral de Twarzen en (β).

𝜑
0.45𝜈
𝛿 ∗∗ 2
= 𝛿0∗∗ 2 + ∫(2𝑈∞ 𝑆𝑒𝑛(𝜑) )5 𝑅𝑑𝑦 … (1)
(2𝑈∞ 𝑆𝑒𝑛(𝜑) )6
0

Como: 𝑥 = 𝑅𝜑 → 𝑑𝑥 = 𝑅𝑑𝜑
Para: x=0 → 𝛿𝑜 = 0 → 𝑃𝑜𝑟 𝑙𝑜 𝑡𝑎𝑛𝑡𝑜: 𝛿 ∗∗ 2 = 0

Resolviendo la Integral:
𝜑
2 1 8
∫(2𝑈∞ 𝑆𝑒𝑛(𝜑) )5 𝑑𝑦 = − cos(𝜑) + 𝐶𝑜𝑠(𝜑)3 − 𝐶𝑜𝑠(𝜑)5 + … (2)
3 5 15
0

Resolviendo (2) en (1):

𝜑
0.45 ∗ 𝜈 ∗ (2𝑈∞ )5 𝑅
𝛿 ∗∗ =√ ∫(𝑆𝑒𝑛(𝜑) )5 𝑑𝑦
(2𝑈∞ 𝑆𝑒𝑛(𝜑) )6
0

0.45 ∗ 𝜈 ∗ (2𝑈∞ )5 𝑅 2 1 8
𝛿 ∗∗ = √ 3 5
6 (− cos(𝜑) + 3 𝐶𝑜𝑠(𝜑) − 5 𝐶𝑜𝑠(𝜑) + 15)
(2𝑈∞ 𝑆𝑒𝑛(𝜑))

0.45 ∗ 𝜈 ∗ 1.57 ∗ 10−5 ∗ 8 2 1 8


𝛿 ∗∗ = √ 6 (− cos(𝜑) + 𝐶𝑜𝑠(𝜑)3 − 𝐶𝑜𝑠(𝜑)5 + )
(2 ∗ 50 ∗ 𝑆𝑒𝑛(𝜑)) 3 5 15

2 1 8
− cos(𝜑) + 3 𝐶𝑜𝑠(𝜑)3 − 𝐶𝑜𝑠(𝜑)5 +
𝛿 ∗∗ √
= 0.000752 ( 5 15 )
𝑆𝑒𝑛(𝜑)6

Para hallar el espesor de desplazamiento 𝛿 ∗ usaremos:

2.61 − 3.75𝑓 + 5.24𝑓 2 ; 𝑓 > 0


∗ ∗∗
𝛿 = 𝐻(𝑓) ∗ 𝛿 𝑑𝑜𝑛𝑑𝑒: 𝐻(𝑓) = { 0.0731 …(3)
0.15+𝑓
+ 2.008 ;𝑓 < 0
Primero calculamos f de (α)
∗∗ 2
𝛿 𝑑𝑢
𝑓=
υ 𝑑𝑥

𝑑𝑢 1 𝑑𝑢 1 𝑑(2𝑈∞ 𝑆𝑒𝑛(𝜑) ) 2
Pero: = = = 𝑈∞ 𝐶𝑜𝑠(𝜑)
𝑑𝑥 𝑅 𝑑𝜑 𝑅 𝑑𝜑 𝑅

Entonces:

𝛿 ∗∗ 2 2
𝑓= 𝑈 𝐶𝑜𝑠(𝜑)
υ 𝑅 ∞
𝛿 ∗∗ 2 2 ∗ 50
𝑓= 𝐶𝑜𝑠(𝜑)
1.57 ∗ 10−5 𝑅

𝑓 = 796178.34𝛿 ∗∗ 2 𝐶𝑜𝑠(𝜑) … (4)

Reemplazando (4) en (3):

2
2.61𝛿 ∗∗ − 29.8 ∗ 105 𝛿 ∗∗ 3 𝐶𝑜𝑠(𝜑) + 5.24 ∗ (79.6 ∗ 104 𝛿 ∗∗ 2 𝐶𝑜𝑠(𝜑)) ;𝑓 > 0
𝐻(𝑓) = { 0.0731
+ 2.008 ;𝑓 < 0
0.15 + 796178.34𝛿 ∗∗ 2 𝐶𝑜𝑠(𝜑)

Finalmente: 𝛿 ∗ sería en (θ):

2.61 − 3.75 ∗ 79.6 ∗ 104 𝛿 ∗∗ 2 𝐶𝑜𝑠(𝜑) + 41.7 ∗ 105 𝛿 ∗∗ 2 𝐶𝑜𝑠(𝜑)2 ;𝑓 > 0


𝛿∗ = { 0.0731𝛿 ∗∗
2 + 2.008𝛿 ∗∗ ;𝑓 < 0
∗∗
0.15 + 796178.34𝛿 𝐶𝑜𝑠(𝜑)

(Espesor de Desplazamiento de la C. L. para el cuerpo cilíndrico)


Al reemplazar 𝛿 ∗∗

Ahora tendremos la Tensión de Rozamiento 𝜏𝜔 :


Necesitamos de:

𝜉(𝑓) = (𝑓 + 0.09)0.62
𝛿 ∗∗ 2 𝑑𝑢
𝜉(𝑓) = ( + 0.09)0.62
υ 𝑑𝑥
𝛿 ∗∗ 2 2
𝜉(𝑓) = ( 𝑈 𝐶𝑜𝑠(𝜑) + 0.09)0.62
υ 𝑅 ∞

De la Ecuación de Twartez:
𝜇 ∗ 𝑢 ∗ 𝜉(𝑓)
𝜏𝜔 =
𝛿 ∗∗
2
𝜇 ∗ 2𝑈∞ 𝑆𝑒𝑛(𝜑) 𝛿∗∗ 2
𝜏𝜔 = ∗∗
[ 𝑈 𝐶𝑜𝑠(𝜑) + 0.09]0.62 … (5)
𝛿 υ 𝑅 ∞

De la Ecuación de Gases Ideales, hallaremos la densidad del aire en las condiciones dadas:

𝑃 101.3 ∗ 103 𝑃𝑎
𝜌= → 𝜌=
𝑅𝑇 𝑃𝑎. 𝑚3
287 𝐾𝑠. 𝐾 ∗ 273.15𝐾

𝜌 = 1.204 𝐾𝑔/𝑚3

1.204𝐾𝑔
𝜇 =𝜌∗𝜈 = ∗ 1.57 ∗ 10−5 𝑚2 /𝑠
𝑚3
𝑁∗𝑠
𝜇 = 1.82 ∗ 10−6
𝑚2

Reemplazando datos en (5):

2
1.82 ∗ 10−6 ∗ 2 ∗ 50𝑆𝑒𝑛(𝜑) 𝛿∗∗ 2
𝜏𝜔 = ∗∗
[ −5
∗ ∗ 50𝐶𝑜𝑠(𝜑) + 0.09]0.62
𝛿 1.57 ∗ 10 8

1.82 ∗ 10−4 ∗ 𝑆𝑒𝑛(𝜑) ∗∗ 2


𝜏𝜔 = [7.961783.44 ∗ 𝛿 ∗ 𝐶𝑜𝑠(𝜑) + 0.09]0.62
𝛿 ∗∗
(Tensión de Rozamiento de la C. L. del Cilindro)

También podría gustarte