INSTITUTO POLITÉCNICO NACIONAL
E.S.I.M.E. UNIDAD PROFESIONAL TICOMÁN
INGENIERÍA AERONÁUTICA
Aerodinámica
Practica No.3:
Distribución de presión alrededor de un cilindro.
Integrantes del Equipo:
Del Valle Salgado Mauricio
Profesor
Moyado Camargo Pablo Rodrigo
Peñalva Martínez Diana Medrano Mejía Antonio
Piña Ramos Cecilia Vianel
Grupo: 4AM1
Fecha de elaboración de la práctica: 16 de octubre de 2015
Fecha de Entrega: 26 de octubre de 2015
Calificación: __________
1 CÁLCULO DE LA DENSIDAD EN EL LABORATORIO
Tabla 1 - Condiciones ambientales.
Condiciones Ambientales
Iniciales Finales Promedio
Temperatura Ambiental 24 °C 26 °C 25 °C
Presión Barométrica 577.5 mmHg 577.5 mmHg 577.5 mmHg
Humedad Relativa 68% 67% 67.5%
Cálculo de la presión barométrica corregida:
𝟏
𝟏 + 𝟎. 𝟎𝟎𝟎𝟎𝟏𝟖𝟒 (𝑻𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒂 ℃)
𝑷𝒄𝒐𝒓𝒓 = 𝑷𝒃𝒂𝒓𝒐𝒎𝒆𝒕𝒓𝒊𝒄𝒂 [ ℃ ]
𝟏
𝟏 + 𝟎. 𝟎𝟎𝟎𝟏𝟖𝟏𝟖 (𝑻𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒂 ℃)
℃
1
1 + 0.0000184 (25℃)
𝑃𝑐𝑜𝑟𝑟 = 577.5𝑚𝑚𝐻𝑔 [ ℃ ]
1
1 + 0.0001818 (25℃)
℃
𝑷𝒄𝒐𝒓𝒓 = 𝟓𝟕𝟓. 𝟏𝟓𝟏𝟓𝟖𝟔𝒎𝒎𝑯𝒈
Conversión de la presión barométrica corregida:
𝑘𝑔𝑓
[1 𝑚𝑚𝐻𝑔] = [13.595101 ]
𝑚2
𝑘𝑔𝑓
13.595101
𝑚2
𝑃𝑐𝑜𝑟𝑟 = 575.151586 𝑚𝑚𝐻𝑔 [ ]
1 𝑚𝑚𝐻𝑔
𝑘𝑔𝑓
𝑃𝑐𝑜𝑟𝑟 = 7819.24218
𝑚2
𝒌𝒈𝒇
𝑷𝒁 = 7819.24218
𝒎𝟐
Conversión de la temperatura:
𝟗
℉= (𝑻𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒂 ℃) + 𝟑𝟐
𝟓
9
℉ = (25 ℃) + 32
5
77℉ = 25℃
Cálculo de la presión de saturación (𝑷𝒔 ) :
𝑷𝒔 = 𝟐. 𝟔𝟖𝟓 + 𝟑. 𝟓𝟑𝟕𝒙𝟏𝟎−𝟑 (𝑻𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒂 ℉)𝟐.𝟐𝟒𝟓
𝑃𝑠 = 2.685 + 3.537𝑥10−3 (77)2.245
𝑙𝑏𝑓
𝑃𝑠 = 63.471444
𝑓𝑡 2
Pero:
𝑙𝑏𝑓 𝑘𝑔𝑓
[1 2
] = [4.8824276364 2 ]
𝑓𝑡 𝑚
𝑙𝑏𝑓
𝑃𝑠 = 63.471444 ⁄𝑓𝑡 2
𝑲𝒈𝒇
𝑷𝒔 = 𝟑𝟎𝟗. 𝟖𝟗𝟒𝟏𝟖𝟒 ⁄ 𝟐
𝒎
Cálculo de la presión de vapor:
𝑃𝑉 = 𝐻𝑟 𝑃𝑠
𝐻𝑟 = 𝐻𝑢𝑚𝑒𝑑𝑎𝑑 𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑎
𝑘𝑔𝑓
𝑃𝑉 = (0.67.5) (296.0600972 )
𝑚2
𝒌𝒈𝒇
𝑷𝑽 = 𝟐𝟎𝟗. 𝟏𝟕𝟖𝟓𝟕𝟒
𝒎𝟐
Cálculo de la densidad:
𝑷𝒛 − 𝟎. 𝟑𝟕𝟕𝟗𝑷𝑽
𝝆𝒛 =
𝒈𝑹𝑻
Donde:
𝑘𝑔𝑓
𝑃𝑧 = 𝑃𝑟𝑒𝑠𝑖ò𝑛 𝑏𝑎𝑟𝑜𝑚𝑒𝑡𝑟𝑖𝑐𝑎 𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎 [ ]
𝑚2
𝑘𝑔𝑓
𝑃𝑉 = 𝑃𝑟𝑒𝑠𝑖ò𝑛 𝑑𝑒 𝑣𝑎𝑝𝑜𝑟 [ ]
𝑚2
𝑚
𝑔 = 𝐴𝑐𝑒𝑙𝑒𝑟𝑎𝑐𝑖ò𝑛 𝑑𝑒 𝑙𝑎 𝑔𝑟𝑎𝑣𝑒𝑑𝑎𝑑 (9.81 2)
𝑠
𝑚
𝑅 = 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒𝑙 𝐴𝑖𝑟𝑒 (29.256 )
°𝑘
𝑢𝑡𝑚
𝜌𝑧 = 𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝐴𝑖𝑟𝑒 [ 3 ]
𝑚
𝑇 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 [𝐾]
𝑘𝑔𝑓 𝑘𝑔𝑓
7819.24218 2 − 0.3779(209.178574 2 )
𝜌𝑧 = 𝑚 𝑚
𝑚 𝑚
(9.81 2 )(29.256 )(25 + 273.15)𝐾
𝑠 °𝐾
𝑢𝑡𝑚
𝜌𝑧 = 0.09050063 [ 3 ]
𝑚
[𝟏 𝒖𝒕𝒎] = [𝟗. 𝟖𝟎𝟔𝟔𝟓 𝒌𝒈]
𝑢𝑡𝑚 9.80665 𝑘𝑔
𝜌𝑧 = 0.09050063 [ 3
]×( )
𝑚 1 𝑢𝑡𝑚
𝑘𝑔𝑚
𝜌𝑧 = 0.88781118
𝑚3
Densidad dentro del laboratorio:
𝒌𝒈𝒎
𝝆𝒛 = 𝟎. 𝟖𝟖𝟕𝟖𝟏𝟏𝟏𝟖
𝒎𝟑
2 MEDICIÓN DE LA PRESIÓN LOCAL ALREDEDOR DEL CILINDRO Y DETERMINACIÓN
DEL COEFICIENTE DE RESISTENCIA AL AVANCE POR DISTRIBUCIÓN DE PRESIONES.
Para obtener la viscosidad cinemática se interpola el valor de la
temperatura promedio en la tabla de atmósfera estándar. Una vez obtenido
𝑚2
el valor de la viscosidad cinemática en , se multiplica por la densidad en
𝑠
𝑈𝑇𝑀 𝑈𝑇𝑀
para obtener la viscosidad absoluta o dinámica en
𝑚3 𝑚∙𝑠
Cálculo del número de Reynolds
Velocidad del viento m/s 35.3115254
Longitud de referencia (diámetro del
cilindro) m
0.0762
Viscosidad del aire a temperatura
ambiente promedio UTM/ms
1.30343× 10−5
Densidad dentro del laboratorio UTM/m3 0.09050063
Re=(ρVL/μ) 183,275.9394
mmH2O Sin mmH2O en 0° mmH2O Real
Flujo
PT 254 134 60.00012724
PE 255 248 3.500007423
q = PT - PE 56.50011982
1-ET 0.999
q corr. 56.4436197
𝐶𝑝𝐶𝑂𝑆𝜃
Determinar el coeficiente de presión total y obtener el producto 2
PL
° (mmH2O) PL Corregida CP 1/2CPCos θ
Sin Flujo 253.5 126.7502688 - -
0 130 65.00013785 1.030973451 0.515486726
10 145 72.50015375 0.898230088 0.442291946
20 184 92.00019511 0.553097345 0.25987067
30 242 121.0002566 0.039823009 0.017243856
40 300 150.0003181 -0.473451327 -0.181342131
50 366 183.0003881 -1.057522124 -0.339880233
60 423 211.5004485 -1.561946903 -0.390485069
70 458 229.0004856 -1.871681416 -0.320073861
80 473 236.5005016 -2.004424779 -0.174029132
90 463 231.5004909 -1.915929204 3.5188E-06
100 432 216.0004581 -1.64159292 0.142533109
110 340 170.0003605 -0.827433628 0.141501229
120 301 150.5003192 -0.482300885 0.120576244
130 292 146.0003096 -0.402654867 0.129411598
140 290 145.0003075 -0.384955752 0.147447314
150 290 145.0003075 -0.384955752 0.16669132
160 289 144.5003064 -0.376106195 0.176712528
170 288 144.0003054 -0.367256637 0.180838813
180 287 143.5003043 -0.35840708 0.17920354
190 288 144.0003054 -0.367256637 0.180838345
200 290 145.0003075 -0.384955752 0.180869502
210 290 145.0003075 -0.384955752 0.166689906
220 290 145.0003075 -0.384955752 0.147445497
230 292 146.0003096 -0.402654867 0.129409332
240 300 150.0003181 -0.473451327 0.118360824
250 343 171.5003637 -0.853982301 0.14603548
260 436 218.0004623 -1.67699115 0.145594466
270 464 232.000492 -1.924778761 -1.06052E-05
280 474 237.0005026 -2.013274336 -0.174812039
290 465 232.5004931 -1.933628319 -0.33068067
300 423 211.5004485 -1.561946903 -0.390495007
310 364 182.000386 -1.039823009 -0.334197712
320 294 147.0003117 -0.420353982 -0.161006681
330 229 114.5002428 0.154867257 0.067060011
340 177 88.50018768 0.615044248 0.28897773
350 141 70.50014951 0.933628319 0.459723361
360 131 65.50013891 1.022123894 0.511061947
Con los datos obtenidos realizar la gráfica polar de la distribución de
presiones alrededor del cilindro en donde la magnitud del radio vector será
el coeficiente de presión para cada grado indicando el número de Reynolds
para el que es válida.
Gráfica Polar
0
360 10
350 1.5 20
340 30
1
330 40
0.5
320 50
0
310 60
-0.5
300 -1 70
-1.5
290 80
-2
280 90
-2.5
270 100
260 110
250 120
240 130
230 140
220 150
210 160
200 170
190 180
Determinar el coeficiente de resistencia al avance por distribución de
presiones utilizando un método de integración aproximada.
𝐶𝑝 cos 𝜃
𝑓(𝜃) =
2
Método Simpson
𝑎−𝑏 0.174533333
ΔX = h = ( )
𝑛
Impares 3.611582539
Pares 0.870842729
CDPR 0.320499862
Coeficiente de resistencia al avance por distribución de presiones
𝐶𝐷𝑃𝑅 = 0.320499862
Para un número de Reynolds de 𝑅𝐸 = 183,275.9394
3 OBTENCIÓN DEL COEFICIENTE DE RESISTENCIA AL AVANCE TOTAL Y DEL
COEFICIENTE DE RESISTENCIA AL AVANCE POR FRICCIÓN DEL CILINDRO
𝐷
𝑪𝑫𝑻𝑶𝑻 =
0.5𝜌𝑉 2 𝑆
𝐷 = 8.2 𝑁
𝐷 = 0.835881753 𝑘𝑔
8.2 𝑁
𝑪𝑫𝑻𝑶𝑻 =
𝑘𝑔𝑚 𝑚2
0.5 × 0.88781118 3 × 35.31152547 2 × 0.03483864 𝑚2
𝑚 𝑠
𝑪𝑫𝑻𝑶𝑻 = 0.42523508
El coeficiente de resistencia al avance por fricción se obtiene de la
diferencia entre el coeficiente de resistencia al avance total y el coeficiente
de resistencia al avance por distribución de presiones.
𝑪𝑫𝑭𝑹𝑰𝑪𝑪 = 𝑪𝑫𝑻𝑶𝑻 − 𝑪𝑫𝑷𝑹
𝑪𝑫𝑭𝑹𝑰𝑪𝑪 = 0.42523508 − 0.320499862
𝑪𝑫𝑭𝑹𝑰𝑪𝑪 = 0.104153789
4 CUESTIONARIO
1. Describa el comportamiento del coeficiente de resistencia al avance del cilindro en
función del número de Reynolds
2. Si se realiza la integración exacta de la ecuación (5) se obtiene 𝑪𝑫 = 𝟎 ¿Cómo explica
usted esto?
3. Explique cómo influye la rugosidad de la superficie del cilindro en su coeficiente de
resistencia al avance por distribución de presiones.
La rugosidad se define como cualquier irregularidad o aspereza sobre una superficie, estas
irregularidades afectan o perturban la capa límite (que es la región adyacente a una superficie sólida
donde se considera que se manifiestan los efectos de la viscosidad) y por lo tanto afectan el flujo;
decimos que el flujo turbulento podría estar en función de la rugosidad de la superficie. Esa región
turbulenta tiende a adherirse mejor a las asperezas, lo que permite que el flujo se desprenda
después y genera entonces una estela menor, reduciendo así su coeficiente de resistencia de
manera notable.
4.- En un experimento se determinó un valor de Cp= 0.55 en un punto situado en la superficie de
𝒎
un cilindro. La velocidad de la corriente libre es igual a 44 𝒔 . Calcular el valor de la velocidad local.
A partir de lo visto en clase tenemos:
Pe∞ + q∞ = Pex + qx
q∞ − qx = Pex − Pe∞
Dividiendo entre q∞
q∞ − qx Pex − Pe∞
=
q∞ q∞
qx
1 − = Cpx
q∞
1⁄ ρV 2
Cpx = 1 − 2 x
1⁄ ρV 2
2 ∞
𝑉𝑥 2
𝐶𝑝𝑥 = 1 − [ ]
𝑉∞
Despejando 𝑉𝑥
𝑉𝑥 = √1 − 𝐶𝑝𝑥 × 𝑉∞
Vx = √1 − 0.55 × 44 m⁄s
𝑽𝒙 = 𝟐𝟗. 𝟓𝟏𝟔𝟎 𝒎⁄𝒔
