Physics">
Informe Dittus Boelter
Informe Dittus Boelter
Informe Dittus Boelter
COEFICIENTE GLOBAL DE
TRANSMISIÓN DE CALOR
2. FUNDAMENTO TEÓRICO
Un intercambiador de calor es un dispositivo donde se produce una transferencia de
calor entre dos fluidos que se encuentran a distintas temperaturas. Dichos fluidos en
ningún momento entran en contacto. El calor se transfiere del fluido de mayor
temperatura (pierde calor y se enfría) al de menor temperatura (gana calor).
Hay diferentes tipos de intercambiadores de calor: de circulación en corrientes
paralelas, en contracorriente o de dos pasos en tubos. Nosotros haremos el
experimento en un intercambiador de flujo paralelo:
T1
t1 t2
T2
• Ecuación de Dittus-Boelter:
𝑁𝑢 = 0,023𝑅𝑒 0,8 𝑃𝑟 𝑛
Con n igual a 0,3 para un enfriamiento y n igual a 0,4 si es un calentamiento. Como
vamos a hacer un calentamiento n=0,4.
Por tanto, como está expresado en el apartado “objetivo” el fin de esta práctica
consiste en demostrar que los coeficientes son correctos. De modo que podría decirse
que partimos de la siguiente expresión:
𝑁𝑢 = 𝑎 ∗ 𝑅𝑒 𝑏 ∗ 𝑃𝑟 0,4
Siendo a y b las variables a calcular.
A continuación, muestro paso a paso el desarrollo teórico de las expresiones y
ecuaciones empleadas en la práctica:
A) Velocidad del aire en la placa perforada
De la ecuación de Bernouilli entre la placa perforada y el exterior:
2 ∗ 𝑔(𝑔𝑟𝑎𝑣𝑒𝑑𝑎𝑑) ∗ Δ𝐻𝑓𝑙𝑢𝑖𝑑𝑜
𝑣𝑜 = 𝐶𝐷 (𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑑𝑒𝑠𝑐𝑎𝑟𝑔𝑎) ∗ √
𝜙
1 − ( 𝜙𝑜 )4
𝐴𝑜 𝑣𝑜 (𝜋𝑟𝑜2 ) 2 ∗ 𝑔 ∗ Δ𝐻 2 ∗ 𝑔 ∗ Δ𝐻 𝑟𝑜4
𝑣𝑡 = = ∗ 𝐶 = 𝐶 ∗ ⇒
𝐴𝑡 (𝜋𝑟𝑡2 ) 𝐷 √ 1 − (𝜙𝑜 )4 𝐷√
𝜙𝑜 4 𝑟𝑡4
1 − (𝜙 )
𝜙𝑡 𝑡
2 ∗ 𝑔 ∗ Δ𝐻𝑎𝑖𝑟𝑒
⇒ 𝑣𝑡 = 𝐶𝐷 √
𝜙
(𝜙𝑜 )4 − 1
𝑡
B) Altura manométrica del aire
Como se puede ver en la ecuación de velocidad del aire en la tubería, para poder
calcular 𝑣𝑡 debemos saber la altura manométrica del aire, y como nosotros medimos
experimentalmente la altura manométrica del agua es necesario hacer una
transformación.
Partimos de la ecuación de diferencia de presión para el agua en un manómetro y la
multiplicaremos por el seno de 60º que es la inclinación del manómetro de agua en el
laboratorio:
C) Reynolds
𝜌∗𝑣∗𝜙
𝑅𝑒 =
𝜇
La densidad se toma a la temperatura de entrada del aire y la viscosidad a la
temperatura media.
D) Nusselt
Como bien dije antes:
ℎ∗𝜙
𝑁𝑢 =
𝑘
1 1 𝑒𝑠𝑝𝑒𝑠𝑜𝑟 𝑡𝑢𝑏𝑒𝑟í𝑎 1
= + +
𝑈 ∗ 𝐴 ℎ(𝑒𝑥𝑡) ∗ 𝐴(𝑒𝑥𝑡) 𝑘∗𝐴 ℎ(𝑖𝑛𝑡) ∗ 𝐴(𝑖𝑛𝑡)
𝑁𝑢 = 𝑎 ∗ 𝑅𝑒 𝑏 𝑃𝑟 0,4
E igualamos a*Pr0,4= a’
𝑁𝑢 = 𝑎′ ∗ 𝑅𝑒 𝑏
Linealizamos:
𝐿𝑛(𝑁𝑢) = 𝐿𝑛(𝑎′ ) + 𝑏 ∗ 𝐿𝑛(𝑅𝑒)
𝑦 = 𝑎 + 𝑏𝑥
De ahí sacamos el valor de “a” y “b”.
Esta representación se hace con los valores de Nusselt experimental y teórico.
G) Error relativo
Obtenemos el error relativo al aplicar su fórmula:
(𝑁𝑢)𝑡𝑒ó𝑟𝑖𝑐𝑜 − (𝑁𝑢)𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙
% 𝑒𝑟𝑟𝑜𝑟 = 100 ∗ (
(𝑁𝑢)𝑡𝑒ó𝑟𝑖𝑐𝑜
- Equipo:
1. Placa perforada
2. Intercambiador de calor
3. Soplante
4. Depósito de agua
5. Manómetro
6. Termómetro de salida y entrada
7. Válvula
- Procedimiento experimental:
Abrimos la instalación poniendo el caudal máximo y esperamos aproximadamente una
hora hasta que se alcance el estado estacionario. Sabremos que estamos en el estado
estacionario cuando la temperatura del aire a la entrada y salida no sufre variaciones.
Después, con la válvula variamos el caudal y con él la altura manométrica del agua. Es
importante saber la altura manométrica mínima a partir de la cual tenemos un flujo
turbulento y podemos emplear la ecuación de Dittus-Boelter. También debemos
conocer la altura manométrica máxima.
Los cálculos de dichas alturas se encuentran en el apartado de cálculos.
El flujo de caudal lo variamos fijándonos en el manómetro. Las medidas del
manómetro a tomar deben encontrarse entre el intervalo de alturas calculado
previamente.
Los datos que recogemos, además de ∆𝐻: temperatura de entrada y salida del aire,
temperatura del agua.
La temperatura del agua es la temperatura del vapor a la presión del laboratorio.
Buscamos las tablas de propiedades termodinámicas del agua e interpolamos con la
presión.
4. DATOS Y CÁLCULOS
- Datos obtenidos:
- Cálculos:
En este apartado se encuentra explicado detalladamente el procedimiento de
cálculos para el primero de los datos experimentales
𝜇 ∗ 10.000
𝑣𝑚í𝑛 =
𝜙∗𝜌
2 ∗ 𝑔 ∗ Δ𝐻𝑎𝑖𝑟𝑒
𝑣𝑡 = 𝐶𝐷 √
𝜙
(𝜙𝑡 )4 − 1
𝑜
𝑚 2 ∗ 9,81 ∗ Δ𝐻𝑎𝑖𝑟𝑒
7,06 = 0,61√ ⇒ Δ𝐻𝑎𝑖𝑟𝑒 = 109,4 𝑚
𝑠 0,022 4
(0,011) − 1
1. Tm = (35+59)/2 = 47 ºC
2. Altura manométrica aire:
B) Coeficiente de convección
2 ∗ 9,81 ∗ 201 𝑚
𝑣𝑡 = 0,61√ ⇒ 𝑣𝑡 = 9,90
0,022 𝑠
(0,011)4 − 1
𝐾𝑔 𝑚 0,022 2 𝐾𝑔
𝑚̇ = 1,145 ( 3
) ∗ 9,90 ( ) ∗ 𝜋 ∗ ( ) = 0,0043
𝑚 𝑠 2 𝑠
2. Flujo de calor
̇ ∗ ∆𝑇
𝑞̇ = 𝑚𝑎𝑖𝑟𝑒 ∗ 𝐶𝑝
Cp a 47 ºC = 1007 J/Kg*K
Diferencia de T = (59-35,5)ºC = 24 ºC
𝐾𝑔 𝐽 ̇ 𝐽 𝐾𝑐𝑎𝑙
𝑞̇ = 0,0043 ( ) ∗ 1007 ( ) ∗ 24 𝐾 = 103,92 = 89,5
𝑠 𝐾𝑔 ∗ 𝐾 𝑠 ℎ
4. Coeficiente de convección:
𝐾𝑐𝑎𝑙 𝑲𝒄𝒂𝒍
89,5 ( ) = ℎ(𝑖𝑛𝑡) ∗ 0,069𝑚2 ∗ 49,44º𝐶 ⇒ 𝒉 ( ) = 𝟐𝟔, 𝟑𝟎
ℎ 𝒉 ∗ º𝑪 ∗ 𝒎𝟐
C) Nusselt
ℎ∗𝜙
𝑁𝑢 =
𝑘
Coef. conducción térmica del aire a la temperatura media (47ºC)= 0,0234(Kcal/m*h*ºC)
26,30 ∗ 0,022
𝑁𝑢 = = 24,73
0,0234
D) Reynolds
𝜌 ∗ 𝑣 ∗ 𝜙 1,145 ∗ 9,90 ∗ 0,022
𝑅𝑒 = = = 12775
𝜇 1,952 ∗ 10−5
E) Nusselt teórico
F) Nusselt experimental
ΔH vt m aire ΔT h int
ΔH (aire)/cm (aire)/m (m/s) Re (kg/s) q (Kca/h) ΔT1 ΔT2 Ln(ΔT1/ΔT2) útil (Kcal/m2*h*ºc)
20135,38 201,35 9,90 12775,03 0,0043 87,82 61,9 38,4 0,48 49,23 25,81
18400,87 184,01 9,46 12212,41 0,0041 85,38 62,3 38,4 0,48 49,40 25,01
15459,75 154,60 8,67 11193,95 0,0038 80,55 63 38,4 0,49 49,70 23,45
13800,66 138,01 8,20 10576,26 0,0036 76,72 63,2 38,4 0,50 49,78 22,30
11312,01 113,12 7,42 9575,29 0,0032 70,30 63,5 38,4 0,50 49,91 20,38
Nu teórico y experimental
4
3,5
2,5
y = 0,8183x - 4,5432
y = 0,8x - 3,8998
Ln(Nu)
2 R² = 0,9967
R² = 1
1,5
0,5
0
9,15 9,2 9,25 9,3 9,35 9,4 9,45 9,5
Ln(Re) experimental
teórico
H) Constantes de Dittus-Boelter
Ahora que sabemos las ecuaciones de la recta despejamos el valor de las constantes.
El valor del Prandlt se toma a la temperatura de entrada del aire al intercambiador y se
considera constante. Hacemos la media de las temperaturas de entrada del aire
(34,62ºC) y el Prandlt es igual a 0,7269.
El Prandlt se obtiene interpolando en las tablas termodinámicas de propiedades del
aire.
- CONSTANTES EXPERIMENTALES
𝑦 = 𝑛 + 𝑚𝑥
ln(𝑁𝑢) = ln(𝑎′ ) + 𝑏 ∗ ln(𝑅𝑒)
ln(𝑎′ ) = −4,5432 = 𝑎 ∗ 𝑃𝑟 0,4 = 𝑎 ∗ 0,72690,4
0,01064
𝑎′ = 𝑒 −4,5432 = 0,01064 = 𝑎 ∗ 0,72690,4 ⟹ 𝑎 = = 0,012
0,72690,4
𝑏 = 0,8183
- CONSTANTES TEÓRICAS
𝑦 = 𝑛 + 𝑚𝑥
ln(𝑁𝑢) = ln(𝑎′ ) + 𝑏 ∗ ln(𝑅𝑒)
I) Error relativo
Con el fin de saber cuánto se desvía el Nusselt experimental del teórico se calcula el
error relativo:
5. CONCLUSIONES
Ecuación de Dittus-Boelter:
𝑁𝑢 = 0,023𝑅𝑒 0,8 𝑃𝑟 0,4
Ecuación de Dittus-Boelter experimental:
𝑁𝑢 = 0,012𝑅𝑒 0,8183 𝑃𝑟 0,4
Ecuación de Dittus-Boelter teórica:
0,023 − 0,012
∗ 100 = 47,8%
0,023
En el caso del segundo coeficiente:
0,8 − 0,8183
∗ 100 = −2,287 %
0,8
Resultado bastante acertado.
En resumen, se ha obtenido los resultados esperados pues ya se sabía de partida que el
error no iba a ser menor al 30-40%