“DETERMINACIÓN EXPERIMENTAL Y PREDICCIÓN DEL COEFICIENTE

DE TRANSFERENCIA DE CALOR ALREDEDOR DEL BULBO

DE UN TERMÓMETRO DE VIDRIO”

I.Q. MANUEL ALEJANDRO MAYORGA BETANCOURT

CÓDIGO 293714

Tesis de grado presentada para optar al título de Magíster en Ingeniería Química

UNIVERSIDAD NACIONAL DE COLOMBIA

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
BOGOTÁ, D.C.
2010
 “DETERMINACIÓN EXPERIMENTAL Y PREDICCIÓN DEL COEFICIENTE
DE TRANSFERENCIA DE CALOR ALREDEDOR DEL BULBO
DE UN TERMÓMETRO DE VIDRIO”

I.Q. MANUEL ALEJANDRO MAYORGA BETANCOURT

CÓDIGO 293714

Tesis de grado presentada para optar al título de Magíster en Ingeniería Química

DIRECTOR: I.Q. HÉCTOR ARMANDO DURÁN PERALTA, M.Sc.

CODIRECTOR: I.Q. LUIS FERNANDO CORDOBA CASTRILLON, PhD.

UNIVERSIDAD NACIONAL DE COLOMBIA

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
BOGOTÁ, D.C.
2010
 Nota de aceptación

_______________________________________
_______________________________________
_______________________________________
_______________________________________
_______________________________________
_______________________________________

Firma del Director de la Tesis

_______________________________________

Firma del Codirector de la Tesis

_______________________________________

Firma del jurado

_______________________________________

Firma del jurado

_______________________________________

Bogotá, D.C., 2010

 Dedico este trabajo al creador de este maravilloso universo, a toda la humanidad que
con su acumulado histórico de conocimientos permite que nos cuestionemos este tipo
de situaciones, a toda mi familia, en especial a mi madre Carmen por todo su amor y
dedicación, a mi Padre por su amor y compañía, a mis queridos hermanos Chuma y Yiyo,
a mi amada Luisa, a mis tías Hersilia y Fanny por su apoyo, al resto de mis ti@s, a mi
querida suegrita Dorita, a mis prim@s, a todos mis parientes, a Javier y Nana por su
amistad incondicional, a Roger y Angela por su sincera aunque distante amistad, a mis
colegas docentes: Alejo, Lucho, Paola y Leito por su compañía e interés –del bueno-, a
tod@s mis verdader@s amig@s, al Pueblo Colombiano que con sus impuestos mantiene
la Educación Pública y en fin a todos los que de una u otra manera se han preocupado
por mis asuntos y/o me han colaborado en el desarrollo de mis actividades cotidianas.
 AGRADECIMIENTOS

Al universo o Dios, por darme la sabiduría, el conocimiento, la fortaleza, el tiempo y los

recursos para desarrollar este trabajo.

Al Ingeniero Armando Duran por su paciencia, colaboración, ingenio y brillantez sin las cuales
este trabajo no se hubiera realizado satisfactoriamente.

Al Señor Edgar Alfonso Martínez Ramírez laboratorista de Lubricantes y Combustibles del

Laboratorio de Ingeniería Química de la Universidad Nacional de Colombia por su incondicional
colaboración, disposición y apoyo.

Al docente Germán Domingo Amat García del Instituto de Ciencias Naturales de la Universidad
Nacional de Colombia por las facilidades prestadas en el uso de equipos a su cargo.

Al Ingeniero Mario Enrique Velásquez, en este momento director del posgrado, por la
motivación y preocupación brindada durante todos estos años y que indudablemente aportó
para la culminación satisfactoriamente de mis estudios.

A Elvira Vargas de Mayorga por su aporte con el invaluable préstamo de la videograbadora.

A los estudiantes: Leonardo, Juan Carlos, Gonzalo y John Jairo por su invaluable colaboración.

A todos las personas que fueron mis docentes durante mis estudios desde el pre-escolar hasta
el posgrado ya que gracias a sus enseñanzas puedo alcanzar este logro.

A todas aquellas personas que estuvieron presentes en cada momento y no dudaron en

brindarme su apoyo y comprensión.

A los ingenieros Hermes Augusto Rangel Jara y Oscar J. Suárez Medina que como jurados
hicieron importantes aportes, correcciones y recomendaciones al presente trabajo.

A la Universidad Nacional donde me he formado y a la ECCI donde me he desempeñado.

 TABLA DE CONTENIDO

Pág.

INTRODUCCIÓN E IDENTIFICACIÓN DEL PROBLEMA ........................................................... 1

JUSTIFICACIÓN ............................................................................................................................ 2

1. GENERALIDADES ................................................................................................................. 4

1.1 TERMÓMETRO DE VIDRIO ................................................................................................... 4

1.1.1 Antecedentes ......................................................................................................................... 4
1.1.2 Descripción General .............................................................................................................. 6
1.1.3 Construcción .......................................................................................................................... 7
1.1.4 Clases .................................................................................................................................... 9
1.1.4.1 Profundidad de inmersión................................................................................................... 9
1.1.4.2 Escala de medición........................................................................................................... 10
1.1.5 Calibración ........................................................................................................................... 11
1.1.5.1 Examen visual .................................................................................................................. 11
1.1.5.2 Puntos de calibración ....................................................................................................... 11
1.1.5.3 Medición de columna emergente ..................................................................................... 12
1.1.5.4 Corrección de la separación de columna ......................................................................... 12
1.1.6 Termómetro para calorímetro.............................................................................................. 13

1.2 COEFICIENTES DE TRANSFERENCIA DE CALOR Y SU IMPORTANCIA ..................... 14

1.3 ANTECEDENTES SOBRE EL COEFICIENTE DE TRANSFERENCIA DE CALOR ............ 16

1.4 CONVECCIÓN LIBRE O NATURAL .................................................................................... 21

1.5 DESCRIPCIÓN DE ALGUNOS MÉTODOS EXPERIMENTALES PARA DETERMINAR

COEFICIENTES DE TRANSFERENCIA DE CALOR EN CONVECCIÓN LIBRE ..................... 23
1.5.1 Métodos desarrollados únicamente con base en la ley de enfriamiento de Newton
................................................................................................................... 24
1.5.2 Métodos basados en las mediciones del campo de velocidades y de temperaturas del
fluido y combinados con la ley de enfriamiento de Newton . ........................... 26
1.5.2 Métodos basados en las mediciones del campo de velocidades y de temperaturas del
fluido y combinados con la ley de enfriamiento de Newton ............................ 26

2. PREDICCIÓN DEL COEFICIENTE .................................................................................... 29

2.1 ESTADO DEL ARTE .......................................................................................................... 29

2.1.1 Análisis de la Convección.................................................................................................. 29
2.1.2 Correlaciones Empíricas ................................................................................................... 31
2.1.2.1 Cilindros Verticales ......................................................................................................... 32
2.1.2.2 Cilindros Horizontales..................................................................................................... 35

2.2 CONSIDERACIONES GENERALES.................................................................................. 36

2.3 MÉTODO DIFERENCIAL ................................................................................................... 40

2.3.1 Resolución Por Semejanza o Análisis Dimensional .......................................................... 46
2.3.1.1 Ecuaciones y Condiciones Límite Adimensionales ........................................................ 47
2.3.1.2 Solución .......................................................................................................................... 48
2.3.2 Resolución por Análisis de Escala .................................................................................... 48
2.3.2.1 Número Prandtl Grande ................................................................................................. 49
2.3.2.2 Número Prandtl Pequeño ............................................................................................... 50
2.3.2.3 Número de Prandt Unitario. ............................................................................................ 51
2.3.2.4 Resultado ....................................................................................................................... 52
2.3.3 Resolución por Teoría de la Capa Límite .......................................................................... 53

2.4 MÉTODO INTEGRAL ......................................................................................................... 56

2.4.1 Planteamiento.................................................................................................................... 56
2.4.1.1 Ecuación de Cantidad de Movimiento ............................................................................ 56
2.4.1.2 Ecuación de Energía ...................................................................................................... 59
2.4.2 Solución ............................................................................................................................. 61
2.4.2.1 Introducción de perfiles de velocidad y temperatura propuestos ................................... 61
2.4.2.2 Solución del sistema U y  .......................................................................................... 66
2.4.2.2.1 Resolución del sistema para n = 0 ............................................................................. 69
2.4.2.2.2 Resolución del sistema para n = 1 .............................................................................. 72
2.4.2.3 Número de Nusselt ......................................................................................................... 74
2.4.2.3.1 Local ............................................................................................................................ 75
2.4.2.3.2 Promedio ..................................................................................................................... 76
3. CONSTRUCCIÓN Y SOLUCIÓN DE LOS MODELOS MATEMÁTICOS .......................... 80

3.1 FUNCIONAMIENTO DEL TERMÓMETRO DE VIDRIO: FLUIDO TERMOMÉTRICO ....... 80

3.1.1 Descripción Fenomenológica ............................................................................................ 80
3.1.2 Construcción...................................................................................................................... 82
3.1.2.1 Termómetro .................................................................................................................... 85
3.1.2.2 Columna ......................................................................................................................... 92
3.1.2.3 Bulbo .............................................................................................................................. 93
3.1.3 Solución ............................................................................................................................. 94

3.2 CONVECCIÓN LIBRE ALREDEDOR DEL BULBO DEL TERMÓMETRO ....................... 98

3.2.1 Termómetro En Posición Vertical Con Flujo Transitorio .................................................... 98
3.2.1.1 Descripción Fenomenológica ......................................................................................... 98
3.2.1.2 Construcción................................................................................................................. 100
3.2.1.3 Solución ........................................................................................................................ 102
3.2.1.4 Parámetro  .................................................................................................................. 106
3.2.1.4.1 Despreciando la resistencia de pared ....................................................................... 106
3.2.1.4.2 Considerando la resistencia de la pared ................................................................... 106
3.2.1.5 Relación entre h y h* .................................................................................................... 107
3.2.1.6 Relación con el perfil de velocidad ............................................................................... 109
3.2.1.7 Número de Nusselt ....................................................................................................... 110
3.2.2 Termómetro En Posición Vertical Con Flujo Cuasiestacionario ...................................... 111
3.2.2.1 Planteamiento............................................................................................................... 111
3.2.2.2 Solución ........................................................................................................................ 111
3.2.2.3 Ecuación cantidad de movimiento ................................................................................ 113
3.2.2.3.1 Relación con el perfil de temperaturas ...................................................................... 113
3.2.2.3.2 Solución ..................................................................................................................... 114
3.2.2.3.3 Perfil de velocidad ..................................................................................................... 116
3.2.2.4 Número de Nusselt ....................................................................................................... 116
3.2.3 Termómetro En Posición Horizontal Con Flujo Cuasiestacionario .................................. 117
3.2.3.1 Consideraciones y Construcción .................................................................................. 117
3.2.3.2 Caso de gradiente angular de temperatura constante ................................................. 118
3.2.3.3 Caso sin gradiente angular de temperatura: ................................................................ 121
4. DETERMINACIÓN EXPERIMENTAL ............................................................................. 122

4.1 ENSAYOS Y SELECCIÓN DE EXPERIMENTOS .......................................................... 122

4.1.1 Ensayo-1: Calentamiento del bulbo con agua a temperatura máxima permisible .......... 122
4.1.2 Ensayo-2: Calentamiento del bulbo con agua a menor temperatura.............................. 123
4.1.3 Ensayo-3 calentamiento del bulbo con aceite ................................................................ 123
4.1.4 Ensayo-4: Enfriamiento del bulbo en agua y aceite a temperatura ambiente ................. 123
4.1.5 Ensayo-5: Enfriamiento del bulbo con aire a temperatura ambiente .............................. 124

4.2 DESCRIPCIÓN DEL EXPERIMENTO ............................................................................ 124

4.2.1 Experimento-1: Convección libre en aire alrededor del bulbo de un termómetro orientado
verticalmente (enfriamiento) ................................................................................................... 125
Figura 4.1 Orientación vertical del termómetro [El Autor ........................................................ 125
4.2.2 Experimento-2: Convección libre en aire alrededor del bulbo de un termómetro orientado
horizontalmente (enfriamiento) ............................................................................................... 126

4.3 MATERIALES E IMPLEMENTOS .................................................................................. 126

4.3.1 Termómetro de Precisión .............................................................................................. 127
4.3.2 Video grabadora y Software .......................................................................................... 129

4.4 MONTAJE ...................................................................................................................... 130

4.4.1 Cilindro Vertical ............................................................................................................. 130
4.4.2 Cilindro Horizontal ......................................................................................................... 130

4.5 PROCEDIMIENTO .......................................................................................................... 131

5. ANALISIS DE RESULTADOS ........................................................................................ 133

5.1 MEDICIONES ................................................................................................................. 133

5.1.1 Cilindro Vertical ............................................................................................................. 134
5.1.2 Cilindro Horizontal ......................................................................................................... 135

5.2 TRATAMIENTO DE DATOS........................................................................................... 136

5.2.1 Determinación del intervalo de temperatura a trabajar .................................................. 136
5.2.2 Parámetros ................................................................................................................... 137
5.2.2.1 Propiedades Físicas ................................................................................................... 137
5.2.2.2 Condiciones de Trabajo.............................................................................................. 138
5.2.2.3 Especificaciones del Bulbo ......................................................................................... 138
5.2.3 Regresión por mínimos cuadrados ................................................................................ 139
5.2.4 Derivación e integración numérica ................................................................................ 140
5.2.5 Corrección por Resistencia de Pared ............................................................................ 142

5.3 RESULTADOS ............................................................................................................... 142

5.3.1 Intervalos de temperatura.............................................................................................. 143
5.3.2 Regresión por Mínimos Cuadrados ............................................................................... 143
5.3.3 Derivación e integración numérica ................................................................................ 143
5.3.4 Valores corregidos por resistencia de pared ................................................................. 144
5.3.5 Valores esperados ........................................................................................................ 144
5.3.5.1 Cilindro Vertical .......................................................................................................... 144
5.3.5.2 Cilindro Horizontal ...................................................................................................... 145
5.3.6 Comparación ................................................................................................................. 145
5.3.6.1 Cilindro Vertical .......................................................................................................... 145
5.3.6.2 Cilindro Horizontal ...................................................................................................... 146

5.4 ANALISIS Y DISCUSIÓN ............................................................................................... 147

6. PRÁCTICA: COEFICIENTE DE TRANSFERENCIA DE CALOR

POR CONVECCIÓN LIBRE ALREDEDOR DE UN CILINDRO ............................................. 152

6.1 OBJETIVOS ................................................................................................................... 152

6.2 FUNDAMENTO TEÓRICO ............................................................................................. 152

6.2.1 Correlaciones Teóricas ................................................................................................. 153
6.2.2 Correlaciones Empíricas ............................................................................................... 154
6.2.3 Modelo Experimental..................................................................................................... 157

6.3 EQUIPO .......................................................................................................................... 159

6.4 PROCEDIMIENTO .......................................................................................................... 161

6.4.1 Ensayos preliminares cualitativos.................................................................................. 161
6.4.1.1 Ensayo-1: Calentamiento del bulbo con agua a temperatura máxima permisible ....... 161
6.4.1.2 Ensayo-2: Calentamiento del bulbo con agua a menor temperatura........................... 162
6.4.1.3 Ensayo-3 calentamiento del bulbo con aceite............................................................. 162
6.4.1.4 Ensayo-4: Enfriamiento del bulbo en agua y aceite a temperatura ambiente.............. 162
6.4.1.5 Ensayo-5: Enfriamiento del bulbo con aire a temperatura ambiente ........................... 162
6.4.2 Experimento Enfriamiento Cilindro Vertical ................................................................... 162
6.4.3 Ensayos Para Calentamiento ........................................................................................ 165

6.5 GRÁFICAS NECESARIAS PARA EL INFORME ........................................................... 165

6.6 CÁLCULOS Y REQUERIMIENTOS PARA EL INFORME .............................................. 166

6.6.1 Parámetros ................................................................................................................... 166
6.6.1.1 Propiedades Físicas ................................................................................................... 166
6.6.1.2 Condiciones de Trabajo.............................................................................................. 167
6.6.1.3 Especificaciones del Bulbo ......................................................................................... 167
6.6.2 Determinación del intervalo de temperatura a trabajar .................................................. 167
6.6.3 Determinación del Coeficiente de Transferencia de Calor y el Número de Nusselt ....... 168
6.6.3.1 Regresión por Mínimos Cuadrados ............................................................................ 168
6.6.3.2 Derivación e integración numérica ............................................................................. 169
6.6.3.4 Corrección por Resistencia de Pared ......................................................................... 171
6.6.3.5 Valores Esperados ..................................................................................................... 171
6.6.3.6 Errores Relativos Porcentuales .................................................................................. 171

6.7 RESULTADOS ............................................................................................................... 171

6.8 ANÁLISIS Y DISCUSIÓN DE RESULTADOS ................................................................ 171

6.9 ASPECTOS, ITEMS Y PREGUNTAS DE DISCUSIÓN .................................................. 172

6.10 CONCLUSIONES Y SUGERENCIAS ........................................................................... 172

6.11 REFERENCIAS ............................................................................................................. 172

7. CONCLUSIONES Y RECOMENDACIONES .................................................................. 173

REFERENCIAS ..................................................................................................................... 180

BIBLIOGRAFIA ...................................................................................................................... 180
INFOGRAFIA ......................................................................................................................... 187

APENDICE A
DEDUCCIÓN DEL NÚMERO DE NUSSELT PROMEDIO
PARA CONVECCIÓN NATURAL ALREDEDOR DE UN CILINDRO VERTICAL .................. 188
I. PRINCIPIO DE SEMEJANZA ......................................................................................... 188
II. ANÁLISIS DE ESCALA ................................................................................................... 190
III. RESOLUCIÓN POR TEORÍA DE LA CAPA LÍMITE ....................................................... 196
IV. MÉTODO INTEGRAL ..................................................................................................... 198

APENDICE B
CONSTRUCCIÓN DEL MODELO PARA EL TERMÓMETRO: FLUIDO TERMOMÉTRICO . 234

APENDICE C ......................................................................................................................... 254

CONVECCIÓN LIBRE ALREDEDOR DEL BULBO DEL TERMÓMETRO ............................ 254
I. TERMÓMETRO EN POSICIÓN VERTICAL CON FLUJO TRANSITORIO...................... 254
II. TERMÓMETRO EN POSICIÓN VERTICAL CON FLUJO CUASIESTACIONARIO ........ 265
III. TERMÓMETRO EN POSICIÓN HORIZONTAL CON FLUJO CUASIESTACIONARIO ... 270

APÉNDICE D
DATOS Y CALCULOS EXPERIMENTALES ......................................................................... 274
 LISTA DE TABLAS
Pág.

Tabla 1.1 Tipos de Termómetros ................................................................................................5

Tabla 1.2 Líquidos Termométricos más usuales .........................................................................8
Tabla 1.3 Rangos aproximados de valores medios del coeficiente de película .........................16
Tabla 1.4 Métodos Generales para evaluar coeficientes de trans. de calor...............................20

Tabla 2.1 Constantes para la ecuación de conv. libre en cilindros verticales ............................33
Tabla 2.2 Constantes para la ecuación de conv. libre en cilindros horizontales .......................35
Tabla 2.3 Valores del exponente para el número de Prandtl .....................................................53
Tabla 4.1 Dimensiones del termómetro usado ........................................................................ 127
Tabla 4.2 Especificaciones ASTM del termómetro usado ....................................................... 128

Tabla 5.1 Experimentos trabajados ......................................................................................... 134

Tabla 5.2 Especificaciones de la pared del bulbo ................................................................... 138
Tabla 5.3 Determinación de la temperatura de la constante de tiempo ................................... 143
Tabla 5.4 Evaluación del coeficiente por regresión con mínimos cuadrados ........................... 143
Tabla 5.5 Evaluación del coeficiente por derivación e integración numérica ........................... 143
Tabla 5.6 Valores corregidos por resistencia de pared ........................................................... 144
Tabla 5.7 Aumentos % debidos a la corrección por resistencia de pared................................ 144
Tabla 5.8 Valores esperados de coef. de trans. de calor para cilindros verticales .................. 144
Tabla 5.9 Condiciones y valores de números adimensionales para el aire usados en cilindros
verticales ................................................................................................................................ 145
Tabla 5.10 Valores esperados de coeficientes de transferencia de calor para cilindros
horizontales............................................................................................................................ 145
Tabla 5.11 Condiciones y valores de números adimensionales para aire usados en cilindros
horizontales............................................................................................................................ 145
Tabla 5.12 Coeficientes de transferencia de calor para cilindros verticales ............................. 146
Tabla 5.13 Errores relativos % de los coeficientes para cilindros verticales ............................ 146
Tabla 5.14 Errores relativos % promedio por ensayo para cilindros verticales ........................ 146
Tabla 5.15 Coeficientes de transferencia de calor para cilindros horizontales......................... 146
Tabla 5.16 Errores relativos % de los coeficientes para cilindros horizontales ........................ 147
Tabla 5.17 Errores relativos % promedio por ensayo para cilindros horizontales .................... 147

Tabla 6.1 Constantes para la ecuación de convección libre en cilindros verticales ................. 154
Tabla 6.2 Constantes para la ecuación de convección libre en cilindros horizontales ............. 157
Tabla 6.3 Especificaciones del termómetro ASTM 56C – 86 [ASTM E1. (2007)] .................... 160
Tabla 6.4 Especificaciones de la pared del bulbo ................................................................... 167

LISTA DE FIGURAS

Pág.

Figura 1.1 Termómetro de vidrio .................................................................................................6

Figura 1.2 Partes del termómetro de vidrio .................................................................................7
Figura 1.3 Tipos de termómetros por profundidad.....................................................................10
Figura 1.4 Termómetro en una bomba calorimétrica .................................................................14
Figura 1.5 Equipo para determinar coef. de película en convección forzada interna .................18
Figura 1.6 Exp. para medir el coef. promedio de trans. de calor por convección. ......................24
Figura 1.7 Experimento de Convección Natural ........................................................................25
Figura 1.8 Esquema de un interferómetro .................................................................................27
Figura 1.9 Interferogramas de isotermas de convección libre en una pared
vertical y en un cilindro inclinado ..............................................................................................28

Figura 2.1 Patrón de Flujo en superficies curvas ......................................................................37

Figura 2.2 Perfil de la capa limite sobre un cilindro vertical caso de caliente y frío ....................39
Figura 2.3 Interferogramas de convección libre externa en cilindros verticales ......................... 40
Figura 2.4 Capas laminares alrededor del cilindro vertical ........................................................42
Figura 2.5 Adimensionalización de los perfiles..........................................................................54
Figura 2.6 Esquema para tomar el elemento dif. de volumen sobre la capa limite ....................56
Figura 2.7 Balance de cantidad de movimiento sobre el elemento dif. de volumen. ..................57
Figura 2.8. Relación entre los espesores de capas límites y sus variaciones............................62
Figura 2.9 Perfil de la capa limite hidrodinámica .......................................................................63
Figura 2.10 Perfiles axiales de temperatura y velocidad y capas limites ...................................65

Figura 3.1 Subsistemas básicos del termómetro .......................................................................82

Figura 3.2. Termómetro ASTM 56C-86 usado en la experimentación .......................................86
Figura 3.3 Columna y capilar del termómetro usado ................................................................92
Figura 3.4 Bulbo del termómetro usado ....................................................................................93
Figura 3.5 Relación instantánea entre los perfiles de temperatura en la interfaz .......................99
Figura 3.6 Circuito de Resistencias Térmicas ......................................................................... 108

Figura 4.1 Orientación vertical del termómetro........................................................................ 125

Figura 4.2 Orientación horizontal del termómetro.................................................................... 126
Figura 4.3 Termómetro de la bomba calorimétrica .................................................................. 127
Figura 4.4 Termómetro ASTM 56C -86 ................................................................................... 128
Figura 4.5 Medición del espesor de la pared del Bulbo con el Estereoscopio ......................... 129
Figura 4.6 Elementos de grabación y edición.......................................................................... 129
Figura 4.7 Montaje para el termómetro vertical ....................................................................... 130
Figura 4.8 Montaje para termómetro horizontal....................................................................... 131

Figura 5.1 Fotogramas cilindro vertical ................................................................................... 134

Figura 5.2 Experimento V-1 ....................................................................................................134
Figura 5.3 Experimento V-2. ...................................................................................................135
Figura 5.4 Fotogramas cilindro horizontal (sentido inverso) .................................................... 135
Figura 5.5 Experimento H-1 ....................................................................................................135
Figura 5.6 Experimento H-2 ....................................................................................................136
Figura 5.7 Significado de la Constante de Tiempo .................................................................. 136

Figura 6.1 Termómetro ASTM 56C -86 ................................................................................... 159

Figura 6.2 Termómetro de la bomba calorimétrica .................................................................. 159
Figura 6.3 Elementos de grabación y edición.......................................................................... 160
Figura 6.4 Orientación vertical del termómetro........................................................................ 162
Figura 6.5 Montaje para el termómetro vertical ....................................................................... 163
Figura 6.6 Orientación horizontal del termómetro.................................................................... 164
Figura 6.7 Montaje para termómetro horizontal....................................................................... 164
 NOMENCLATURA

SIMBOLOS

Número de Biot Conductividad térmica

Número de Boussinesq Difusividad térmica
Número de Fourier Viscosidad cinemática
Número de Grashof Viscosidad dinámica
Número de Mach Calor específico a volumen constante
Número de Nusselt Calor específico a presión constante
Nusselt promedio
Número de Prandtl Potencia, flujo de trabajo
Número de Rayleigh Flujo de calor
Número de Reynolds Coeficiente de transferencia de calor (T.C.)
Coeficiente promedio de T.C.
Altura , Coeficiente global de T.C.
Radio Resistencia térmica
Diámetro
Longitud (Largo o alto) Cambio, paso. delta
Longitud característica , Espesor capa límite térmica
Área de transferencia de calor Espesor capa límite hidrodinámica
Volumen Relación entre capas límites
Distancia vertical adimensional
Masa , Distancia axial adimensional
Tiempo Temperatura adimensional
Dimensión espacial axial, altura Velocidad adimensional
Dimensión espacial radial Velocidad de escala
Dimensión espacial angular , Parámetros de perturbación o de curvatura
Temperatura
Masa del fluido termométrico
Velocidad Volumen del fluido termométrico
Flujo másico Velocidad fraccional
Tiempo adimensional
Fuerza(s) Número adimensional
Peso Número adimensional
Esfuerzo cortante Longitud fraccional
Κ Coeficiente adimensional de dilatación térmica
, Presión Velocidad de cambio de energía interna
Energía Volumen adimensional o fraccional
Energía interna específica
h Entalpia específica , Variable de asignación
, Coeficientes o constantes ley de potencias
Densidad , Funciones de la relación ⁄
,v Volumen específico , Funciones del Número de Prandtl
, Coeficiente de expansión térmica , Constante o Coeficiente
Coeficiente de comprensibilidad , ,… Integrales adimensionales
 Parámetro vinculante Masa molar
Vector unitario de escala Constante Universal de los Gases Ideales
Gradiente angular de temperatura
Índice de correlación
Aceleración gravitacional , Parámetros de linealización

SUBINDICES

Bulbo Pared Plana

Columna Radial
De película
, , Contadores Dimensión espacial angular
Inicial ∞ Del seno del fluido
Medio Adimensional
Radial, referencia Inicial
sonido
A Volumen Atmosférica
De pared Exterior, externa
Axial Flotación
Interior
Convectiv@ maxim@
Diametral Sistema
Longitudinal (Largo o alto) Viscosa

SUPERINDICES Y ÉNFASIS

, , , Exponentes de correlación potencias Exponente de correlación

* Adimensional Promedio
 INTRODUCCIÓN E IDENTIFICACIÓN DEL PROBLEMA

Convección es la forma de transporte de calor a través de un medio fluido, ya sea gas o líquido,
entre una temperatura en la superficie (correspondiente a la interfase fluido-sólido) y una
temperatura en el seno del mismo fluido. Tal manera de transporte de calor por un medio
material está gobernada por la ecuación llamada “Ley de Enfriamiento de Newton” que plantea
que el flujo de calor es directamente proporcional al área de la superficie que sirve de interfase
entre el fluido y el sólido, y a la diferencia entre estas dos temperaturas. A la constante de
proporcionalidad se le denomina coeficiente de transferencia de calor. En realidad esta ley no
es más que la definición del coeficiente de transferencia de calor [INCROPERA].

Dicha constante también recibe el nombre de coeficiente convectivo de transferencia de calor o

coeficiente de película y su valor depende de la geometría de la interfase, del régimen de flujo
y naturaleza del fluido, así como de las dos temperaturas. Esta multiplicidad de variables
independientes resulta porque la transferencia de calor por convección está determinada por
las capas limites que se generan sobre la superficie.

Para el diseño y apreciación de equipos donde tenga lugar fundamentalmente una

transferencia de calor por convección es necesario conocer el valor de dicho(s) coeficiente(s).

No se han desarrollado determinaciones experimentales de coeficientes de transferencia de

calor usando termómetro de vidrio. La determinación por este camino implica el modelamiento
a partir del balance de energía aplicado al respectivo sistema del termómetro.

El termómetro de vidrio es un instrumento que se encuentra disponible en el laboratorio, y a

parte de indicar la temperatura, puede también servir para la medición indirecta del coeficiente
de transferencia de calor a partir de los cambios de temperatura que experimenta el aparato en
el estado transitorio hasta alcanzar el equilibrio térmico. Esto puede a futuro conllevar al diseño
y construcción de un instrumento que reporte el valor local de un coeficiente de transferencia
de calor, aparato de medición que no se dispone.

La medida de coeficientes locales de transferencia de calor, y por analogía de masa, puede

tener notables aplicaciones en el diseño y optimización de equipos (como por ejemplo los

1
biorreactores ya que la heterogeneidad de los sistemas y la sensibilidad los hacen complejos
de modelar), para lo cual es deseable conocer los perfiles espaciales de ciertas propiedades.

Con el presente trabajo se busca principalmente determinar experimentalmente y predecir el

coeficiente de transferencia de calor alrededor del bulbo de un termómetro de vidrio. Para
lograr dicho objetivo se hace necesario: Plantear y solucionar modelos del funcionamiento del
termómetro de vidrio; determinar teóricamente el coeficiente de transferencia de calor alrededor
del bulbo del termómetro para determinadas situaciones; confrontar el grado de ajuste de los
modelos propuestos contra los datos experimentales; analizar la factibilidad del diseño del
instrumento de medición directo de coeficientes de transferencia de calor fundamentando en el
modelo propuesto; y diseñar por lo menos una práctica relacionada con fenómenos de
transporte (masa, calor y/o cantidad de movimiento) y la posibilidad de un diseño experimental
por analogía para transporte de masa.

JUSTIFICACIÓN

En la ingeniería química se realizan operaciones unitarias, algunas de las cuales involucran la

transferencia de calor por medio de fluidos siendo la convección el mecanismo predominante.

En el diseño y cálculo de equipos utilizados en ingeniería química normalmente se requiere

conocer el valor del coeficiente de transferencia de calor por convección.

Los coeficientes se pueden determinar por medición indirecta para un caso particular de flujo
en donde se requiere medir el área, las temperaturas, la velocidad y el flujo de calor,
generalmente usando medios potenciométricos.

No se dispone de una práctica experimental para la determinación local de dichos coeficientes,

razón por la que se propone un modelo que, a partir de la descripción del funcionamiento del
termómetro de vidrio en estado transitorio, permita lograr este objetivo en ciertas situaciones
físicas particulares. Además no se encuentra en la literatura un modelo para el termómetro de
vidrio que parece ser reserva de los fabricantes.

2
Hasta donde se tiene conocimiento no existe un aparato que realice el reporte directo del valor
del coeficiente, y tampoco hay registros del uso de termómetro de vidrio para lograr dicha
determinación, aunque si se reportan desarrollos experimentales que hacen uso de
termopares y termocuplas.

3
 1. GENERALIDADES

El presente trabajo plantea la determinación experimental del coeficiente de transferencia de

calor alrededor del bulbo de un termómetro. Este capítulo incluye elementos de carácter
general tanto de información teórica como experimental necesarios para la comprensión y
desarrollo del problema propuesto en este proyecto.

Como el problema involucra directamente un termómetro de vidrio, es preciso incluir aspectos

generales sobre el termómetro de vidrio.

1.1 TERMÓMETRO DE VIDRIO

1.1.1 Antecedentes

En la Antigüedad Clásica, periodo Greco-Romano, se realizaron los primeros registros de

medición de temperatura con instrumentos muy básicos sin escala y abiertos a la atmósfera,
influenciados así por la presión atmosférica y la temperatura.

Durante el Renacimiento, en la Era de los Descubrimientos, debido al auge por la investigación

en las ciencias naturales, se dió inicio al desarrollo de los termómetros de líquido en vidrio,
junto con su escala; de hecho fue el propio Galileo el que llamó la atención sobre la dilatación
de los líquidos y sugirió utilizar la expansión del alcohol como indicador de la temperatura,
entonces el antecesor del termómetro fue el termoscopio (1606) de Galileo Galilei (1564-1642);
hasta que el también italiano Sanctorius (1561-1636) le incorporó una escala con graduación
numérica, apareciendo así el termómetro de la escuela florentina que consistió en un bulbo
lleno de alcohol sellado por su parte superior.

Estos termómetros se usaron durante casi un siglo, pero presentan la desventaja que cuando el
alcohol se enfría, moja el vidrio y así dificulta una lectura correcta. Es Edmund Halley (1656-
1742) quien sugiere el uso del mercurio, pues aunque se dilata menos que el alcohol, no se
adhiere a las paredes de vidrio; el problema de la menor dilatación es resuelto haciendo más
estrecho el tubo del termómetro. Pero fue Fahrenheit (1686-1736), fabricante de termómetros,

4
quien difunde el uso del termómetro de mercurio, además de concebir la escala de temperatura
denominada con su nombre y que se impuso en ese momento [PEREZ].

Actualmente existen numerosos tipos de termómetros o sensores de temperatura que pueden

ser usados para diversas y específicas aplicaciones de medición; en la Tabla 1.1 se indican
algunos de los tipos más usuales, así como sus características más importantes.

RANGO CARACTERÍSTICAS
TIPO DE
NOMINAL COSTO LINEALIDAD NOTABLES
TERMÓMETRO
[ºC]
Termómetro de Simple, lento y de
-10 a 300 Bajo Buena
mercurio en vidrio lectura manual

Termorresistencia
(Pt, Ni) RTD -100 a 600 Medio Alta Exactitud

Bajo Alta Requiere referencia

Termocupla -150 a 1500
de temperatura

Termistor -15 a 115 Medio No lineal Muy sensible

Fácil conexión a
Integrado lineal Medio Muy alta sistemas de toma de
datos

Termómetro de gas -20 a 100 Medio Buena No muy versátil

Diodos -200 a 50 Bajo Alta Bajo costo

Tabla 1.1 Tipos de Termómetros

A pesar de la gran variedad de instrumentos para la medición de temperatura y de la

prohibición del mercurio, los termómetros de líquido en vidrio son frecuentemente utilizados por
su facilidad de uso, costo, portabilidad, estabilidad, amplio intervalo de trabajo y/o por
recomendaciones en normas. Pero aun así existe confusión al momento de su elección, uso y
calibración, ya que aún dentro de estos termómetros se fabrica una amplia gama de acuerdo a
determinadas aplicaciones específicas [ASTM].

5
1.1.2 Descripción General

El funcionamiento de los termómetros de vidrio con líquido se basa en la dilatación térmica de

un líquido termométrico introducido en un recipiente de vidrio de paredes delgadas (recipiente
termométrico, bulbo) comunicado con un capilar (capilar del termómetro), al cual se le ha
asignado una escala. El líquido actúa como un traductor, convierte el calor en un trabajo de
expansión mecánica. Con el incremento de la temperatura, el líquido y el vidrio del termómetro
se expanden con diferente coeficiente de expansión, causando que el líquido avance por el
tubo capilar. El más común de los termómetros es el de mercurio (Figura 1.1), que basa su
funcionamiento en la marcada diferencia entre las expansiones térmicas del mercurio .
y vidrio . .

Figura 1.1 Termómetro de vidrio [El Autor]

Al líquido termométrico en el capilar se le denomina columna y su menisco sirve como marca

de columna o indicador. Existen varios líquidos que pueden servir como fluidos termométricos
como Mercurio, Mercurio-Talio, Toluol, Alcohol, Mezcla de Pentano. La forma del menisco es:
para el mercurio, la parte superior de la curva o convexa, mientras que para líquidos orgánicos,
la parte inferior o cóncava; cuando las fuerzas adhesivas son mayores que las fuerzas
cohesivas, el menisco tiende a ser cóncavo mientras que cuando las fuerzas cohesivas son
superiores a las adhesivas, el menisco es convexo. El menisco del mercurio se debe al ángulo
de contacto del mercurio con el vidrio que varía dependiendo si está avanzando o
retrocediendo la columna y también a la variación de la tensión superficial. Esto combinado con
la pequeña compresibilidad del mercurio y la elasticidad del vidrio genera una pequeña
variación de la posición del menisco, especialmente en termómetros muy sensibles.

El capilar del termómetro consta del capilar de medida, el capilar de acoplamiento y posibles
ensanchamientos:

6
 Figura 1.2 Partes del termómetro de vidrio [www.metas.com]

 El capilar de medida es la parte del capilar del termómetro provista de una escala. El
capilar de acoplamiento es la sección del capilar que une el bulbo con el capilar de
medida.
 Los ensanchamientos de contracción son ampliaciones en el capilar, mediante los
cuales se suprime o se interrumpe la indicación de una determinada zona, que permiten
la instalación de una escala adicional, o que impiden el descenso de la columna al bulbo
(como en el caso de los termómetros clínicos).
 El ensanchamiento de expansión (burbuja de seguridad) es una ampliación, al final del
capilar, que hasta cierto punto protege al termómetro en caso de sobrecalentamiento.

1.1.3 Construcción

En la manufactura de tales termómetros, el capilar es marcado en dos puntos (0 °C y 100 °C) y

luego se hacen graduaciones uniformes entre las marcas, donde se asume que el volumen de
una masa fija de mercurio en el vidrio, es una función lineal de la temperatura. El error
resultante de esta consideración (+0.12 °C a 50 °C) es ordinariamente más pequeño que el
debido a las variaciones en el capilar.

7
Los materiales usados en la fabricación del termómetro son: vidrio, líquido termométrico y gas.

Todos los tipos de vidrios no son apropiados para el trabajo termométrico; actualmente se han
desarrollado vidrios termométricos especiales, con buena estabilidad y resistentes a las
temperaturas, generalmente boro-silicato, pyrex, schott, etc.

El líquido termométrico debe tener las siguientes propiedades físico-químicas:

• Tener un coeficiente de expansión lineal

• Tener color o ser opaco, para su fácil lectura.
• Tener un menisco bien definido, para fácil lectura.
• Mantenerse como líquido en el intervalo nominal del termómetro.
• No “mojar” por fuerzas adhesivas la superficie del capilar.
• Ser químicamente inerte con respecto a otros materiales en el sistema.
• Ser térmica y químicamente estable.
• Por seguridad no ser dañino durante su manufactura y uso.

Coeficiente de Expansión
Líquido Intervalo Típico
Típico Relativo en Vidrio
Termométrico °C
°C-1
Mercurio -38…+600 0.000 16
Etanol -80…+60 0.001 04
Tolueno -80…+100 0.001 03
Pentano -200…+300 0.001 45
Tabla 1.2 Líquidos Termométricos más usuales [www.metas.com]

En la mayoría de los casos encima de liquido termométrico, en el capilar, hay un relleno de gas
protector, sometido a mayor presión (con el fin de evitar evaporación del líquido termométrico)
constituido por un gas seco y libre de oxígeno (por ejemplo Nitrógeno, Argón), para dificultar
una separación o rebosamiento del líquido termométrico y para elevar su punto de ebullición,
solo en algunos termómetros especiales se evacúa el capilar por encima de la columna. Así,
todos los termómetros para altas temperaturas deben ser llenados con un gas inerte seco tal
como nitrógeno presurizado para prevenir separación del mercurio a cualquier temperatura
indicada en la escala. Los termómetros de inmersión total graduados sobre los 150 °C deben

8
ser llenados con gas para minimizar la destilación del mercurio de la parte superior de la
columna. Para termómetros graduados debajo de los 150 °C el llenado con gas es opcional
pero altamente recomendado.

Las especificaciones para el diseño y elaboración de termómetros se pueden encontrar es

normas internacionales como las de la ASTM y la ISO. Estos termómetros están diseñados
para ser usados con métodos de prueba específicos, y son identificados por tener el acrónimo
ASTM y un número inscrito. Importante la norma de especificaciones ASTM E-1 en donde se
enlistan todas las características, los puntos de calibración y las temperaturas de la columna
emergente para termómetros de inmersión parcial y un procedimiento especial que debe
seguirse para su prueba y calibración. Las normas más importantes mas a tener en cuenta son:

 ASTM E1. (2007). Standard Specification for ASTM Liquid-in-Glass Thermometers.

 ASTM E77. (2007). Standard Test Method for Inspection and Verification of
Thermometers.
 ASTM E2251. (2007). Standard Specification for Liquid-in-Glass ASTM Thermometers
with Low-Hazard Precision Liquids.
 ISO 386. (1977). Liquid-in-glass laboratory thermometers—Principles of design,
construction and use.
 ISO 1770. (1981). Solid-stem general purpose thermometers.
 ISO 4795. (1996). Glass for thermometer bulbs.

1.1.4 Clases

Los tipos básicos del termómetro de vidrio dependen fundamentalmente de la profundidad de

inmersión y de la ubicación de la escala de medición, las que se relacionan estrechamente con
la medición específica requerida o funcionalidad.

1.1.4.1 Profundidad de inmersión

La elección del tipo de inmersión depende de la medición requerida, la profundidad del fluido y
el tipo de montaje.

9
  Termómetro de inmersión parcial: indica la temperatura correctamente cuando el
bulbo y una porción específica de la columna están inmersos en el medio a la
temperatura que va a ser medida, por lo que tiene una línea que indica la profundidad
de inmersión exacta.

 Termómetro de inmersión total: indica la temperatura correctamente cuando el bulbo

y la columna entera del líquido (unos cuantos milímetros por arriba del nivel del líquido)
son inmersas en el medio para iniciar la medida.

 Termómetro de inmersión completa: Está diseñado para indicar la temperatura

correctamente cuando todo el termómetro, incluyendo la cámara de expansión, están
expuestos en el medio a la temperatura que va a ser medida.

Figura 1.3 Tipos de termómetros por profundidad [www.metas.com]

1.1.4.2 Escala de medición

 Termómetro de barra: La escala se encuentra directamente sobre la pared gruesa del

capilar.

10
  Termómetro de oclusión: La escala se encuentra sobre un soporte independiente
separado del capilar; el capilar y el soporte de la escala se encuentra encerrados en un
tubo de revestimiento.
 Existen variaciones de los tipos anteriores.

1.1.5 Calibración

Esta se realiza por comparación contra un termómetro patrón calibrado, en sistemas térmicos
con recirculación de líquido o sales, o en un lecho fluidizado, realizando la medición directa de
los puntos fijos secundarios de fusión del hielo o ebullición del agua. De acuerdo a la exactitud
que se desee obtener, se utiliza durante la calibración un termómetro de resistencia de platino
u otro tipo de termómetros trazables y con baja incertidumbre.

Cuando se requiera una precisión menor que el 1% de la escala requerida, es aconsejable una
calibración, la cual señala las correcciones para agregar o sustraer a las lecturas. Las
calibraciones pueden ser hechas por organismos de normalización.

Otra importante fuente de error en los termómetros es la mala lectura en línea del menisco, la
que puede ser provocada por la posición adoptada por el ojo cuando se hace la lectura, y que
puede ser eliminada con el uso de un teodolito. Las lecturas de termómetros muy sensibles son
también fuertemente dependientes de la presión, donde los coeficientes de presión pueden ser
tan altos como de 0.1 °C/atm.

1.1.5.1 Examen visual

Permite la detección de: separación de la columna de mercurio, pequeñas bolitas de mercurio a

lo largo de la columna capilar, oxidación y fallas o fisuras en el vidrio. Estas fallas pueden
perjudicar las lecturas del termómetro. También debe ser examinada en la escala: el grosor de
las líneas, la posible desigualdad en las divisiones y graduaciones borrosas.

1.1.5.2 Puntos de calibración

11
El termómetro debe calibrarse en todo su intervalo, incluso si lo posee, un punto de referencia.
Por lo usual, se calibra en puntos distribuidos uniformemente en el intervalo de indicación la
escala principal.

La cantidad de puntos de calibración depende del intervalo, división “mínima” de la escala y

exactitud deseada. El intervalo entre los puntos de calibración no debe ser demasiado largo o
innecesariamente pequeño tal que destruya la confiabilidad en correcciones interpoladas en
valores de temperatura. La recomendación es que sean distribuidos entre 40 a 100 divisiones,
hasta 40 para patrones y hasta 100 para ordinarios.

1.1.5.3 Medición de columna emergente

A veces se requiere utilizar termómetros de inmersión total como inmersión parcial y la

temperatura del ambiente sobre el baño, o equipo donde esté colocado el termómetro, puede
ser diferente de la temperatura del bulbo del termómetro. Debe por tanto hacerse una
corrección por la diferencia de temperaturas entre el bulbo y la columna emergente. Lo mismo
ocurre para los termómetros de inmersión parcial, las correcciones aplican solo para la
profundidad de inmersión indicada y una sola temperatura de la columna emergente. Si el
termómetro es utilizado en otras condiciones es necesario corregir la columna emergente. La
columna emergente puede ser medida utilizando un termómetro Faden o un determinado
número de termómetros de vástago corto.

1.1.5.4 Corrección de la separación de columna

Estos son varios métodos para unir la columna del líquido:

 Golpear: Se sujeta, ligeramente por debajo de la mitad y con el bulbo abajo, el

termómetro con una mano, y luego se golpea suavemente la orilla con la palma de la
otra mano.

 Centrifugar: Se sujeta igual que en el ítem anterior y con el brazo extendido en

posición horizontal se balancea rápidamente hacia abajo, y enseguida se para cuando
el termómetro alcance la posición vertical; también es válido usar una centrifugadora
con el bulbo orientado hacia afuera.

12
  Calentar: Es frecuentemente el método más rápido para unir la columna separada y
simplemente consiste en el calentamiento del bulbo del termómetro. Para termómetros
que no están graduados sobre los 250 °C, la cámara de expansión también puede servir
para unir la columna: se sujeta el termómetro en posición vertical y se calienta
lentamente, sin usar flama directa, hasta que el segmento de columna separada y una
porción de la columna principal, entra a la cámara de expansión; luego debe ser
removido antes de que el líquido termométrico suba demasiado a la cámara de
expansión ya que puede provocar un daño en el termómetro debido a un exceso de
presión.

 Enfriar: Sujetando verticalmente el termómetro, se enfría el bulbo hasta una

temperatura que permita retraer todo el mercurio hacia el bulbo; luego se regresa lenta
y cuidadosamente a temperatura ambiente.

1.1.6 Termómetro para calorímetro

Existen termómetros especiales para los calorímetros, donde se desea una medición de
diferencia de temperatura con alta precisión (de 0.01 a 0.001 °C). Para estos termómetros la
finura de la escala de graduación tiende a ser pequeña debido a la precisión que se desea
obtener en las medidas de temperaturas. La escala puede tener un error de varias decimas de
grado, pero este error se cancela cuando se hace la diferencia. Un termómetro típico de una
bomba calorimétrica tiene un rango de 19 a 35 ° , con graduaciones de 0.02 ° . Para
mediciones de los descensos del punto de congelamiento con agua o benceno como solventes,
un rango de -2 a 6 °C con graduaciones de 0.01 °C es conveniente. Tales termómetros
requieren un adecuado manejo. No son solamente frágiles, también son susceptibles de ciertos
daños (separación de la columna de mercurio, burbujas en el bulbo) derivados principalmente
de la extrema finura del hilo, cuando se mantienen verticales durante un rápido calentamiento o
enfriamiento. Si la columna de mercurio llega a separarse, se debe enfriar el bulbo en una
mezcla de hielo-sal para que el mercurio quede enteramente en el bulbo y así darle salida a las
burbujas hacia la parte superior de la columna. Luego se debe permitir que el termómetro
llegue a la temperatura ambiente en posición vertical.

13
 Figura 1.4 Termómetro en una bomba calorimétrica [SHOEMAKER]

1.2 COEFICIENTES DE TRANSFERENCIA DE CALOR Y SU IMPORTANCIA

Para la ingeniería en general y para la ingeniería química en particular y especial es

sumamente importante determinar la razón de transferencia de calor para una determinada
diferencia de temperaturas. Ello sucede con frecuencia cuando se diseñan componentes y
sistemas en los que sucede este proceso. Para estimar el costo, la factibilidad y el tamaño de
equipos requerido para transferir cierta cantidad de calor en un determinado tiempo hay que
realizar un riguroso análisis de transferencia en los cuales el ingeniero debe identificar las
resistencias entre los puntos de alta y baja temperatura, muchas de las cuales se encuentra en
fase fluida donde la convección juega un papel predominante.

Las dimensiones de refrigeradores, calentadores, evaporadores, condensadores e

intercambiadores de calor no solo dependen de la cantidad de calor sino también de la
velocidad a la que se va a transferir esta energía. Para conocer entonces uno de estos flujos
requeridos se hace necesario disponer del valor de transferencia de calor por convección.

14
De la misma manera, la operación adecuada de los componentes de un equipo tales como las
paredes o las aspas de una turbina depende de la posibilidad de enfriar rápidamente una
superficie por evacuación continua del calor en determinado tiempo.

La optimización de ciertos procesos exige conocer el valor de las resistencias que obstaculizan
uno de las fenómenos de transferencia ya sea de masa, calor o momento, porque si se quiere
acelerar o desacelerar cierta transferencia se debe actuar en el diseño sobre la resistencia
controlante. En el caso del calor, si esta resistencia se encuentra en fase fluida, dicho valor es
el reciproco del coeficiente de transferencia de calor.

El parámetro que describe todo el fenómeno de transferencia de calor por convección es el

coeficiente de transferencia de película el cual depende de gran cantidad de variables. Se
considera un fenómeno bastante complejo debido a la conducción a través de la película de
fluido adherida a la superficie del sólido por las fuerzas atractivas y también a la transferencia
de masa (advección). Este fenómeno se puede dificultar en el caso de los gases en los cuales
puede ser importante también la transferencia de calor por radiación de manera simultánea a la
convección.

Así que para su determinación y estudio, es necesario clasificar la convección de acuerdo a las
siguientes posibilidades:

 Forzada, Natural o Mixta: Si el principal generador de la turbulencia es un agente

externo o un gradiente de densidades al interior del fluido, o una combinación de
ambos.
 Interna o Externa: El fluido puede circular al interior o alrededor de un ducto.
 Con cambio o sin cambio de fase: Los cambios de fase entre líquido y gas.
 Régimen Laminar, Transición o Turbulento: Se refiere al grado de turbulencia que
puede ocurrir y que es cuantificado por el número de Reynolds o el número de Grashof
ya sea por convección forzada o natural respectivamente.
 Naturaleza del fluido: Este puede ser un líquido (incompresible) o un gas. De la misma
manera se pueden tener casos para ciertos valores de números de Prandtl.
 Tipo de Superficie: Plana o Curva.

15
En la siguiente tabla se presentan algunos valores aproximados y típicos de coeficientes
medios de transferencia de calor.

MECANISMO Y FLUIDO h, W/m2K

Condensación mixta de vapor de agua 5.000-150.000
Ebullición mixta de agua 4.000-100.000
Condensación natural de vapor de agua 5.000-35.000
Ebullición natural de agua 2.000-25.000
Condensación mixta de vapores orgánicos 900-20.000
Ebullición mixta de líquidos orgánicos 600-20.000
Convección forzada, metales líquidos 500-25.000
Convección forzada, aceites 50-3.000
Convección forzada, aire a 200 bar 250-1.000
Convección forzada, aire a presión atmosférica 15-300
Convección natural sin cambio de fase, metales líquidos 100-1.000
Convección natural sin cambio de fase, agua 80-800
Convección natural son cambio de fase, aceites 20-200
Convección natural sin cambio de fase, aire a presión atmosférica 5-20
Tabla 1.3 Rangos aproximados de valores medios del coeficiente de película [SIGALES]

1.3 ANTECEDENTES SOBRE EL COEFICIENTE DE TRANSFERENCIA DE CALOR

Isaac Newton debió su interés por la temperatura, el calor y el punto de fusión de los metales al
ejercicio de la supervisión de la calidad de la acuñación mientras trabajó en la casa de la
moneda de Inglaterra. Observó que al calentar al rojo un bloque de hierro y tras retirarlo del
fuego, el bloque se enfriaba más rápidamente cuando estaba muy caliente, y más lentamente
cuando su temperatura se acercaba a la temperatura del aire. Así, a principios del siglo XVIII,
publicó en un trabajo sobre este hecho [SIGALÉS]:

“Partes iguales de aire eran calentadas en intervalos iguales de tiempo, recibiendo un grado de
calor proporcional al calor del hierro.”

Sin embargo la expresión sencilla de esta observación, conocida como “Ley de Enfriamiento de
Newton”, la formuló 121 años después Joseph Fourier en su Theorie Analytique de la Chaleur:

16
Esta fórmula propuesta, a fin de facilitar el tratamiento matemático, exige suponer que , el
coeficiente superficial de transferencia de calor, es constante, aunque ya anteriormente se
sabía que este no lo era en el seno del fluido, como lo demuestra los experimentos con
resultados cualitativos de John Leslie en 1804 junto con el primer estudio que distinguió
cuantitativamente la radiación de la convección en el aire y la variabilidad de con la
temperatura de la superficie realizado por Louis Dulong y Alexis Thérèse Petit en 1818.

A continuación se tratan las distintas formas de determinar dicho coeficiente. En la mayoría de

los casos en que está implicada la transferencia de calor en un fluido, para diversas
configuraciones geométricas y arreglos de flujo, los coeficientes de transferencia de calor son
calculados por correlaciones empíricas y semiempíricas entre grupos adimensionales de las
variables pertinentes. Esas correlaciones de datos experimentales están reportadas en la
literatura, en las cuales se especifica la geometría, los regímenes de flujo, el tipo de
convección y la naturaleza química del fluido (fase, polaridad, etc.). Adicionalmente se indica el
grado de precisión de la correlación. Generalmente se acude a la determinación empírica del
número de Nusselt, debido en parte, a que las ecuaciones diferenciales que describen el
modelo de transporte de calor en convección, no tienen solución analítica, y sólo la tienen para
los casos sencillos.

No es práctico, pero si costoso, hacer ensayos de experimentos con muchos líquidos y

gases bajo, una gran variedad de condiciones experimentales. Ello, para tener disponibles los
valores numéricos de las magnitudes de los coeficientes de transferencia de calor, necesarios
para aplicaciones industriales.

Así la obtención de los valores mediante otros métodos se hace evidente. Una de las formas
más populares y empleadas consiste en las expresiones obtenidas mediante el análisis
dimensional que conlleva a una experimentación con menor número de pruebas. El análisis
dimensional asocia las variables en grupos adimensionales y propone una estructura de cómo
se relacionan dichos grupos, conocidos como números adimensionales. Los experimentos
permiten establecer los valores de los exponentes y coeficientes de las ecuaciones
adimensionales para un determinado intervalo de trabajo. Entonces el coeficiente de

17
transferencia de calor (que se encuentra inmerso dentro de unos de los grupos) puede ser
calculado a partir de la correlación para cualquier combinación de velocidad, tamaño de la
geometría y propiedades de fluido que esté dentro del intervalo.

Para la determinación del coeficiente de transferencia de calor para líquidos que fluyen dentro
de tuberías o tubos, un aparato típico usado consiste en un sistema formado por un
intercambiador de calor principal, -sección de tubería de prueba encerrada por un tubo
concéntrico o de doble tubo -, y un intercambiador auxiliar que sirve para efectuar la operación
opuesta de la sección de prueba, como se muestra en la figura [KERN].

Figura 1.5 Equipo para determinar coeficientes de película en convección forzada interna [KERN]

Una vez que se generan los datos experimentales se diseñan los experimentos para determinar
las correlaciones. Lo usual, es reducir el número de variables, utilizando análisis adimensional,
esto conduce a que el coeficiente de transferencia de calor, que está involucrado en el número
de Nusselt, quede principalmente en función del Número de Reynolds, el número Prandtl, etc..
y sí se presenta convección libre, interviene adicionalmente el número de Grashof en la
correlación.

Existen cinco métodos generales para determinar los coeficientes de transferencia de calor por
convección [KREITH]:

18
 MÉTODO DEFINICIÓN VENTAJAS DESVENTAJAS
Análisis Determinación de grupos Es matemáticamente simple. Resultados obtenidos

dimensional adimensionales para Es fácil de usar. incompletos e inútiles sin

combinado correlacionar datos Facilita la interpretación de datos experimentales.

experimentales, es decir datos experimentales al Contribuyen poco a la
con
combina varias variables en correlacionarlos en función de comprensión del fenómeno,
experimentos
grupos adimensionales y no los grupos adimensionales. es decir no proporciona
produce ecuaciones que información relativa a la
puedan ser resueltas. naturaleza del fenómeno.
No se pueden omitir variables
pertinentes, es decir se deben
seleccionar adecuadamente.
Se debe disponer de una
teoría preliminar o de un
conocimiento del fenómeno.
Las correlaciones no son
universales en cuanto al tipo
de flujo y geometría.
Soluciones Solución simultánea de las Las suposiciones hechas en Requiere la comprensión

matemáticas ecuaciones de cantidad de el análisis pueden suficiente de los mecanismos

movimiento y transferencia de especificarse con precisión y físicos para describirlos en
exactas de
energía que describen el comprobar su validez lenguaje matemático.
las
fenómeno. mediante experimentación. La mayoría de ecuaciones
ecuaciones
Sirve como base de son bastante complejas.
de capa
comparación para métodos
límite aproximados más simples.
Se pueden establecer
soluciones numéricas con
ordenadores de alta
velocidad.
Análisis Utiliza una ecuación simple Evita la detallada descripción Se pierde precisión y

aproximados pero que describe matemática del flujo en la exactitud en los resultados.
esencialmente las capa límite.
de las
distribuciones de velocidad y Aporta soluciones a
ecuaciones
temperaturas en la capa situaciones que no pueden
de capa
limite, para hacer un análisis tratarse mediante un análisis
límite
macroscópico. matemático exacto.
mediante Se extiende a flujo turbulento.
métodos Concuerda en lo fundamental
integrales. con los resultados obtenidos
por otros métodos.

19
La analogía Un movimiento compuesto Útil para análisis de Este modelo no

entre con dirección perpendicular al transferencia turbulenta. necesariamente debe

transferencia flujo medio explica la El mecanismo de concordar con las condiciones

transferencia de energía y transferencia puede reales existentes en la
de calor y
momentum. Este movimiento describirse como un modelo naturaleza.
transferencia
puede describirse simplificado.
de cantidad
estadísticamente por medio Los resultados
de de un método similar al de la experimentales concuerdan
movimiento. teoría cinética de los gases. con las predicciones
analíticas basadas en el
modelo hipotético.
Solución Resolución numérica de las Resuelve la necesidad de

numérica ecuaciones exactas de expresar las variables de

movimiento, a través de la campo en función del espacio
discretización de las mismas. y/o del tiempo.
Calcula con facilidad las
soluciones para distintas
condiciones de frontera.
Tabla 1.4 Métodos Generales para evaluar coeficientes de transferencia de calor [El Autor].

La evolución del uso de computadores para el análisis de la convección, puede sintetizarse así:
a finales de los cincuenta para la resolución de ecuaciones diferenciales ordinarias asociadas a
flujos autosimilares, se usan los recién inventados computadores digitales; en los sesenta se
desarrollan eficientes métodos numéricos para la resolución de ecuaciones diferenciales
parciales que rigen el comportamiento del flujo en capas límite y en un conducto, al punto que
las soluciones se vuelven rutina una vez se tengan los programas montados; en los setenta se
concentra la atención en la resolución de las ecuaciones de conservación que rigen los flujos
con vórtices, ondas de choque y características complejas; en la actualidad se ha desarrollado
software especializado en la dinámica de fluidos que permiten estudiar una amplia gama de
flujos, además de la programación interactiva compleja que permite implantar nuevas
configuraciones geométricas y de contorno [MILLS].

Sin embargo en necesario aclarar, que como el objeto de la transferencia de calor es el cálculo
del campo de temperaturas y del flujo de calor, el proceso de determinar el coeficiente de
transferencia de calor a partir de los balances fundamentales es superfluo, como lo señaló
Eugene F. Auditori, quien en 1974 argumentó que el coeficiente al ser función de la
temperatura y altamente no lineal, sobre todo en los gases y cambios de fase, no es racional

20
como parámetro, por lo que esgrimió que debían plantearse ecuaciones que arrojaran
directamente el flujo de calor, pero no tuvo la perspicacia de enfocar el problema a partir de los
balances diferenciales. De hecho los métodos modernos y avanzados de cálculo de
transferencia de calor a través de un fluido prescinden de . Actualmente, se presentaría un
inconveniente práctico si se decide abandonar el uso del coeficiente de transferencia de calor:
la inmensa cantidad de información con base experimental acumulada en expresiones que
permiten obtener un valor fiable y muy fácil de usar [SIGALÉS].

La predicción del coeficiente entre una superficie curva y un fluido, se dificulta en comparación
a la predicción en superficies planas, por lo que se hace necesario manejar métodos
numéricos, que pueden conducir a fórmulas aproximadas, o en caso de preferir una solución
más exacta, exige el uso de métodos de solución con funciones ensayo, que permita obtener
los perfiles de temperatura necesarios, para hallar el cambio local de los perfiles en la interfase
fluido-sólido.

Para resolver las ecuaciones diferenciales del transporte de calor y de cantidad de movimiento,
por lo general se suele proponer soluciones de prueba en series de potencias o con parámetros
de perturbación. Otra forma es obtener una expresión, a partir de la solución del problema
utilizando los datos generados por diferencias finitas o elementos finitos, que puede conducir a
una formula al correlacionar los resultados numéricos generados, los cuales pueden ser
presentados en tablas y gráficas.

1.4 CONVECCIÓN LIBRE O NATURAL

La transferencia de calor por convección se incrementa por el aumento de la turbulencia; esta

turbulencia puede ser generada por un agente externo como una bomba, un ventilador, un agitador,
un compresor, etc. caso que se conoce como la convección forzada, o puede generarse por un
desequilibrio de temperaturas en el seno del fluido que provoca por tanto un gradiente de
densidades que da origen a una fuerza de empuje ascensional o de flotación, en donde el fluido
caliente tiende a elevarse; claro está que en algunos casos las diferencias de densidad también
pueden ser provocadas por cambios en la composición [MILLS]. Los gradientes de densidad de
un fluido puro se relacionan con los de temperatura por medio de los coeficientes de expansión
volumétrica.

21
En el presente trabajo se propone un sistema de estudio, que corresponde a la transferencia de
calor por convección libre, alrededor del bulbo de un termómetro rodeado por aire. El
movimiento del aire es generado por la diferencia de temperatura en el aire, lo cual origina un
gradiente de densidad en el fluido, ocasionando fuerzas de flotación en el aire caliente, el cual
asciende originándose la corriente convectiva natural. Para el estudio de la transferencia de
calor alrededor del bulbo del termómetro, puede asumirse, que el sistema a estudiar puede
representarse aproximadamente como la convección libre alrededor de un cilindro.

Las magnitudes de las velocidades asociadas a la convección natural suelen ser sumamente
pequeñas, sin superar los 2 m/s [MILLS], por lo que los coeficientes de transferencia de calor
suelen ser normalmente inferiores a los de convección forzada; un valor típico en gases es de
alrededor de 5 W/m2K, y se debe analizar si la convección por radiación que sucede
simultáneamente es despreciable o no. En convección libre también hay grados de turbulencia
que determinan diferentes regímenes flujo: laminar, transición y turbulento.

Aunque estos coeficientes de transferencia de calor por convección libre son relativamente
pequeños, muchos dispositivos dependen de este modo de transferencia para su enfriamiento
como el caso de las unidades electrónicas; debido al calor interno las temperaturas de estos
cuerpos se elevan por encima de la ambiental, conforme la diferencia de temperaturas se
incrementa, también lo hace la transferencia de calor, hasta que se alcanza un estado de
equilibrio en el que la razón de generación de calor es igual a la razón de disipación de calor.

La convección natural es el mecanismo de flujo de calor predominante en los radiadores de

vapor, los muros de edificios, el cuerpo humano en aire -ambos en reposo-, el almacenamiento
de materiales, etc., por lo que para determinar la carga de calor de dispositivos calefactores,
aire acondicionado y artefactos electrónicos, y también para cuantificar las pérdidas de calor en
los tubos que transportan vapor u otros fluidos calientes, o en el enfriamiento de superficies
contenedoras de reacciones nucleares por fisión de calor, etc. se hace necesario conocer los
coeficientes de transferencia de calor por convección libre [KREITH].

En lugar del número de Reynolds empleado en convección forzada, en la convección Libre o

Natural la turbulencia es cuantificada por el número de Grashof, , o el número de Rayleigh
.

22
El número de Grashof indica la proporción o razón entre las fuerzas de empuje y las fuerzas
inerciales, y se puede calcular a partir de las propiedades del fluido y la geometría, de esta
manera:
| | | |

Mientras que el Número de Rayleigh se define como el producto entre los número de Grasohf,
, y Prandtl, .

| | | | | |
· ·

| | | |

El coeficiente de transferencia de calor por convección se encuentra a partir del número de

Nusselt, , que relaciona el calor convectivo con el calor conductivo de la siguiente manera:

Usualmente la longitud característica, , puede ser la longitud, , o el diámetro, , del cilindro

dependiendo si se orienta de manera vertical u horizontal, respectivamente.

Si el cilindro se encuentra orientado verticalmente, entonces lo más común es manejar un

número de , y , donde la longitud, , corresponde a la altura del cilindro.

| | | |
;

Ahora si se encuentra orientado horizontalmente, entonces se trabaja un número de , y

:
| | | |
;

1.5 DESCRIPCIÓN DE ALGUNOS MÉTODOS EXPERIMENTALES PARA DETERMINAR

COEFICIENTES DE TRANSFERENCIA DE CALOR EN CONVECCIÓN LIBRE

Las soluciones analíticas son complejas y difíciles, para muchas configuraciones, por lo que
mucha de la información de diseño viene de estudios experimentales. De lo consultado en la

23
literatura referente a los métodos experimentales para determinar el coeficiente de
transferencia de calor, puede observarse y clasificarse dos métodos generales:

 Métodos desarrollados únicamente con base en la ley de enfriamiento de Newton

.
 Métodos basados en las mediciones del campo de velocidad y de temperaturas del fluido y
combinados con la ley de enfriamiento de Newton .

1.5.1 Métodos desarrollados únicamente con base en la ley de enfriamiento de Newton

.

Figura 1.6 Experimento para medir el coeficiente promedio de transferencia de calor por convección.
[INCROPERA]

Para flujo externo, un arreglo experimental que suele usarse para una determinada geometría
de la superficie se muestra en la Figura 1.6. El sistema es calentado eléctricamente por una
resistencia interna, para mantener una diferencia de temperatura entre la superficie sólida y el
fluido, la cual se mide utilizando termocuplas: una situada en la superficie sólida y otra situada
en el seno del fluido, con lo cual puede determinarse la diferencia de temperatura .
Adicionalmente es necesario medir experimentalmente directa o indirectamente:

- El flujo de calor que atraviesa el sistema en estado estacionario, que corresponde a la

potencia suministrada por la red eléctrica al elemento calefactor (que es función de la
corriente y el voltaje del circuito).

24
- El área de transferencia de calor y la diferencia de temperatura entre la superficie sólida y el
fluido .

Posteriormente utilizando la Ley de Enfriamiento de Newton, se determina experimentalmente

el coeficiente como

Figura 1.7 Experimento de Convección Natural [KARLEKAR]

Para la determinación experimental del coeficiente de transferencia de calor en convección libre

(así como los factores que influyen en el mismo) entre un gas y una placa plana vertical u
horizontal, se muestra en la figura 1.7 un dispositivo que consta de un gran tanque dentro del cual
se encuentra suspendida verticalmente u horizontalmente una superficie. Se llena el tanque con
un gas transparente como nitrógeno, oxígeno o aire. Para calentar uniformemente la
superficie se inserta una resistencia eléctrica, la cual está conectada a una fuente de voltaje.
Se colocan termopares sobre la superficie y al interior del tanque, el cual está equipado con dos
válvulas para purgar y cargar un gas cualesquiera. Se varía la transmisión de energía eléctrica con
un reóstato y se miden la corriente y el voltaje para evaluar el flujo de calor suministrado
[KARLEKAR].

Se cierra el interruptor, se entrega una determinada carga de calor y luego se deja transcurrir el
tiempo suficiente para que la temperatura se estabilice. Despreciando las transferencias de
calor por los alambres y elementos de soporte por conducción, se calcula el flujo de calor por

25
radiación, el cual se descuenta de la potencia calorífica total entregada, para así hallar el flujo
de calor por convección, el cual permite hallar el coeficiente de transferencia de calor, con una
exactitud satisfactoria.

Si se incrementa el flujo de calor, aumenta la temperatura de la superficie y por consiguiente

debe también crecer el valor del coeficiente de transferencia de calor. Una variación del
coeficiente también se observa cambiando el gas, la geometría o el tamaño de la superficie, por
ejemplo para una superficie vertical doblar la altura puede llevar a una reducción del 20% del
coeficiente aproximadamente [KARLEKAR]

1.5.2 Métodos basados en las mediciones del campo de velocidades y de temperaturas

del fluido y combinados con la ley de enfriamiento de Newton .

En estos métodos mediante interferogramas y métodos de visualización de flujos se pueden

visualizar tanto las líneas isotermas como las líneas de corriente y con ello es posible inferir el
campo vectorial de velocidades asociado. Así puede medirse experimentalmente la velocidad y
la temperatura del fluido en contacto con una superficie de temperatura diferente. Conocida
experimentalmente la distribución de temperatura y la distribución de velocidad, puede
determinarse el espesor de la capa límite térmica e hidrodinámica y con toda esta información
utilizando las ecuaciones de balance de cantidad de movimiento y la ecuación de energía y
combinando con la ley de enfriamiento de Newton se logra finalmente determinar el coeficiente
de transferencia de calor.

Similar a la convección forzada, hay casos donde, por la configuración estudiada, normalmente
se dificulta hallar analíticamente los perfiles de temperaturas y velocidades, se hace entonces
necesario recurrir a las medidas experimentales para obtener dichas relaciones. La convección
natural o forzada, en régimen turbulento, requiere de datos experimentales, pero el problema
es más agudo en la convección natural que en la forzada, dado que las velocidades son
normalmente tan pequeñas que son realmente difíciles de medir [HOLMAN].

A pesar de las dificultades experimentales, se han realizado medidas de velocidad utilizando

técnicas de burbuja de hidrógeno [VLIET], anemometría de hilo caliente [CHEESEWRIGHT] y
anemómetros de fibra de cuarzo. La anemometría láser [FLACK] es particularmente útil para
las medidas de convección natural, ya que no se perturba el campo fluido.

26
Para la determinación del campo de temperaturas se usa le técnica del interferómetro de
Zehnder-Mach. Un interferómetro típico, que se utiliza para estudios de transferencia de calor
consta de una fuente de luz monocromática, una fuente para colimar el haz de luz, un par de
placas separadoras, un par de espejos y una pantalla. El haz colimado pasa a través de una
placa separadora A de tal modo que cambia su curso 90°. Este rayo desviado, rayo 1, se refleja
90° y pasa a través de la sección de prueba, que contiene un cuerpo caliente, tal como una
placa. Ya que la placa se encuentra a una temperatura mayor que la del aire de la vecindad,
existe una variación de temperatura dentro del aire que rodea la placa. Además, ya que la
densidad del aire depende de su temperatura, el haz 2 pasa a través de aire que tiene
gradientes de densidad. La longitud de onda de la luz, varia con la densidad del aire, así que la
distancia óptica para el haz 2 es diferente que para el haz 1. Estos dos haces contienen rayos
que difieren en una longitud de onda o mas, y además de otros rayos que difieren en media
longitud de onda o múltiplo. La placa separadora B combina estos rayos y el haz combinado se
proyecta en la pantalla. Además, ya que los haces 1 y 2 parten de la misma fuente, ambos son
coherentes y pueden interferir el uno con el otro. El haz combinado que parte de la segunda
placa separadora es una onda nueva de la misma frecuencia y longitud de onda pero del doble
de amplitud del haz colimado original, dando por resultado franjas oscuras y claras alternadas
sobre la pantalla [KARLEKAR].

Figura 1.8 Esquema de un interferómetro [KARLEKAR]

27
El interferómetro es básicamente un instrumento que se usa para medir el índice de refracción,
el cual está relacionado con la densidad. Por medio de una ecuación de estado de los gases,
se calcula la temperatura partiendo de la densidad y presión conocidas. Una región de líneas
oscuras muy juntas representa una región de cambios de temperatura rápidos y viceversa.

Estos campos se pueden apreciar en la Figura 1.9. Cabe anotar que en superficies calientes
verticales los interferogramas muestran que en el inferior de la superficie existe una capa límite
cuyo espesor no es cero, como suele asumirse en la solución analítica. Además muestran que
para una lamina o un cilindro vertical las isotermas presenta un perfil parabólico respecto al eje
axial, y que en cilindros horizontales calientes el campo de temperatura no es simétrico
respecto al eje, y las líneas tienden a alejarse a la parte superior del cilindro.

Figura 1.9 Interferogramas de isotermas de convección libre en una pared

vertical (izq.) [HOLMAN] y en un cilindro inclinado (der.) [MARTINENKO]

28
 2. PREDICCIÓN DEL COEFICIENTE

El Capítulo 1 trató sobre generalidades y algunos elementos conceptuales descriptivos

relacionados con la parte teórica y experimental de la transferencia de calor por convección. En
este capítulo, el trabajo continúa con los métodos teóricos para determinar coeficientes de
transferencia de calor en convección sobre sistemas cilíndricos. Se aplicaran estos métodos
para generar algunas ecuaciones teóricas del número de Nusselt para convección libre y
externa en cilindros, siendo expresiones que se encuentran disponibles en la literatura.

En el Capítulo 3 los proponentes de esta tesis generaran sus propias correlaciones teóricas
originadas al aplicar los métodos teóricos al problema específico que se quiere modelar en este
proyecto.

Se expone la revisión del estado del arte en torno a la predicción y luego la descripción de los
métodos diferencial e integral aplicados a algunos casos de interés del cilindro con el fin de
mostrar el dominio de los métodos.

2.1 ESTADO DEL ARTE

A partir de una minuciosa revisión de las ecuaciones, correlaciones y expresiones de la

literatura se extrajeron las distintas opciones ya sean de origen teórico (Análisis de la
Convección) o de origen práctico (Correlaciones Empíricas) que permiten predecir el
coeficiente de transferencia de calor para el caso de estudio: Convección libre alrededor de
cilindros verticales y horizontales.

2.1.1 Análisis de la Convección

L. Lorenz (1881) fue el primero en resolver el problema de transferencia de calor en una pared
vertical calentada, por convección libre, postulando que el fluido cercano a la superficie se
mueve en forma vertical hacia arriba, donde el componente horizontal de la velocidad son
insignificantes.

29
Schmidt y Beckman en 1930 realizaron un trabajo donde mostraron que los postulados de
Lorenz debían refinarse. Los datos de Saunder confirmaron dicho cuestionamiento.

Las ecuaciones diferenciales para convección libre en una pared plana fueron resueltas por
Ostrach (1952) con número de Prandtl desde 0.01 a 1000. Pohlhaunsen resolvió el modelo
para el aire introduciendo una función corriente y una parámetro de similitud, y encontró que el
valor máximo de velocidad se encontraba dentro de la capa límite. Este último también convirtió
las ecuaciones diferenciales parciales que gobiernan el fenómeno en una ecuación diferencial
ordinaria. Squire presentó una solución aproximada usando el método integral [KARLEKAR].

Para un cilindro vertical y basándose en la hipótesis aproximada de Langmuir de capa limite

estacionaria de una superficie curva, Elenbaas (1948) planteó:

/
⁄
.

Si ∞ se tiene que

⁄
.

solución que corresponde a una pared plana.

Los métodos de similaridad fueron luego abordados para resolver las ecuaciones diferenciales
parciales de conservación en su forma integral, en el caso de un cilindro.

Así Sparrrow y Gregg (1956) obtuvieron la solución a las ecuaciones de conservación

empleando métodos de similaridad y encontraron que cuando ⁄ no es tan grande no se
pueden ignorar los efectos de curvatura. De hecho el criterio para tratar un cilindro vertical igual
que una pared plana vertical fue dado además por Gebhart:

Las soluciones distan menos del 5% de la solución de pared plana para Pr desde 0.72 a 1.0.

LeFevre y Ede (1956) emplearon un método integral con el fin de solucionar las ecuaciones de
gobiernan el fenómeno obteniendo la siguiente expresión del número de Nusselt promedio en
régimen laminar [BEJAN]:

30
 ⁄
·

Donde tanto y estan basados en la altura, , del cilindro.

Los cilindros horizontales y las esferas son configuraciones que no dan lugar a similitud para
encontrar la solución de las ecuaciones que gobiernan el fenómeno, por lo que se hace
necesario abordar otros métodos.

Hermann (1936) dedujo una expresión para cilindros horizontales en régimen laminar
resolviendo las ecuaciones diferenciales por el método de no similaridad, estableciendo que
para aire:

⁄
.

Merk y Prins (1953-1954) obtuvieron para cilindros horizontales en régimen laminar, una
solución empleando el método integral y asumiendo espesores iguales entre las capas limite
térmica e hidrodinámica [BEJAN]:

⁄ ⁄

Para valores de de 0.7 (aire) el valor de la constante, fue calculado como 0.436.

Minkowycz y Sparrow (1974) obtuvieron resultados analíticos para cilindros verticales

isotérmicos, aplicando el método de no similaridad local y así lograron una solución numérica
⁄
para Prandtl 0.733. Demostraron que si ⁄ . el coeficiente de transferencia de
calor de un cilindro vertical isotérmico es a lo mas 4 veces el calculado para una placa vertical.

Cebeci (1974) generó resultados numéricos sobre cilindros delgados, mientras que cuerpos de
simetría axial fueron abordados en una revisión realizada por Gebhart (1988).

2.1.2 Correlaciones Empíricas

Para obtener una correlación empírica el procedimiento experimental debe ser repetido para
una variedad de condiciones de prueba, en las que se puede cambiar la velocidad del fluido
, la(s) longitud(es) característica(s) de la geometría, la naturaleza del fluido, tal que tengan
sustancialmente números distintos de números de Prandtl como aire, agua y aceite. Los

31
distintos resultados de los números de Nusselt se linealizan con respecto a Grashof y Prandtl,
que pueden ser ajustados en una correlación cuya forma más usual es:

, Pr ó Pr

Donde , y son constantes determinadas en la mayoría de casos experimentalmente y en

contadas configuraciones son deducidas. Estas constantes son independientes de la
naturaleza del fluido, ya que ésta se cuantifica en el número de Prandtl.

La limitación principal de las correlaciones es que raramente proveen valores exactos para la
coeficientes convectivos, ya que una gran variedad de condiciones que dependen del tipo del
fluido, la turbulencia del mismo y la rugosidad de la superficie, no son muy tenidas en cuenta,
por lo que se pueden incurrir en errores hasta del 25% al usar estas expresiones [BLAIR].

Se aclara que los valores de estos grupos adimensionales pueden ser locales o promedios.
Para aclarar esto se escriben los valores locales y promedios de la siguiente manera:

Axiales:

| | | |
;

| | | |
;

Radiales:

| | | |
;

| | | |
;

Mientras que para los valores promedio del número de Nusselt, se acostumbra o no a colocar
una barra:

2.1.2.1 Cilindros Verticales

32
Para cilindros verticales con un flujo de calor constante se tiene [JANNA]:

C n

0.05 0.93 0.05

4
0.05-10 1.37 0.16
4
10 0.6 0.25
Tabla 2.1 Constantes para la ecuación de convección libre en cilindros verticales

Todos los valores de las propiedades se determinan a la temperatura de película, , que

corresponde al promedio aritmético entre la temperatura de pared, , y la temperatura del
seno del fluido , excepto el coeficiente de expansión térmica, , que en el caso de gases se
determina a .

Se recuerda que para un gas ideal como el aire

Una correlación general recomendada por Yang para regiones en régimen laminar y turbulento
tiene la forma [POPIEL]:

. ⁄
.
. ⁄ ⁄
. ⁄

Los valores de los cálculos numéricos obtenidos por Cebeci para un parámetro de curvatura

√
.

pueden ser representados por la ecuación simple, para . , en la región laminar

[POPIEL]:

33
 .
. .
, .

donde , corresponde al de una pared plana vertical, dado por la correlación teórica
[KARLEKAR]:

⁄ ⁄ ⁄
, . .

o también por la correlación empírica desarrollada por Churchill y Chu para una placa vertical
de altura, , y con un borde frontal abrupto para flujo laminar [MILLS] :

⁄
, . .

donde la función del número de Prandtl, , definida por Churchill y Usagi se define como:

⁄ ⁄
.

Los valores de Cebeci también pueden ser aproximados por una ecuación validada para
números de Prandtl entre 0.01-100, es decir, . [POPIEL]:

. .
,

donde:

.
. .

. ⁄
. . . .

Popiel obtuvo una correlación a partir del enfriamiento de un cilindro en estado transitorio que
tiene una incertidumbre promedio de 3.5% con un 95% de confianza, válida para:

, , y, ⁄ :

donde:

. . . .

34
y

. . .

La dependencia generalizada de transferencia de calor por convección libre en un gas tanto

para flujo laminar y como turbulento alrededor de un cilindro vertical está dada por la expresión
de Kuehn [MARTINENKO]:

⁄
⁄ ⁄ ⁄
. ⁄ . ⁄ .

que es válida para:

. ⁄ , , /

2.1.2.2 Cilindros Horizontales

Para y . se correlacionaron los datos experimentales de Eckert y

Soehnghen [KREITH]:

⁄
.

En forma más general, para un cilindro isotérmico, Morgan sugiere una expresión de la forma
[INCROPERA]:

RaD C n
10-10-10-2 0.675 0.058
-2 2
10 -10 1.020 0.148
102-104 0.850 0.188
4 7
10 -10 0.480 0.250
7 14
10 -10 0.125 0.333
Tabla 2.2 Constantes para la ecuación de convección libre en cilindros horizontales

35
Churchill y Chu dan la siguiente correlación válida para [MILLS]:

⁄
.
. ⁄ ⁄
. ⁄

Cuando ocurre transición de una capa limite laminar a una turbulenta y el aumento
del número de Nusselt con el número de Rayleigh es mayor [MILLS]:

⁄
.
. ⁄ ⁄
. ⁄

Por el método de la capa limite, Kuehn deduce [MARTINENKO]:

⁄
⁄ ⁄ ⁄
. . ⁄ .

En la región de régimen laminar para , la aproximación que realiza Churchill se puede

aplicar con un error del 4.5% [MARTINENKO]:

⁄
⁄ ⁄
. ⁄ ⁄
. ⁄

2.2 CONSIDERACIONES GENERALES

Existe una apreciable cantidad de problemas de interés práctico acerca de la transferencia de

calor por convección natural sobre cuerpos de forma arbitraria y compleja. En la mayoría de los
casos, la superficie en la que se desarrolla el flujo es curvilínea y el ángulo , entre el vector de
gravedad, y su proyección a un plano tangente a la superficie, está cambiando como lo hace
el flujo que se mueve desde el borde de la entrada a lo largo de la superficie. En algunos
casos, la curvatura de la superficie es pequeña, y entonces la superficie se puede tratar
aproximadamente como plana y los numerosos datos disponibles en superficies planas
verticales e inclinadas se pueden utilizar.

36
 Figura 2.1 Patrón de Flujo en superficies curvas [MARTINENKO]

Las características especiales de transferencia de calor en una superficie curvada también se

manifiestan en una delgada capa límite. Los efectos causados por la curvatura del cuerpo son
especialmente significativos en pequeños y moderados números de Grashof. En la capa
adyacente a la superficie, el principal mecanismo de transferencia de calor es por conducción.

El coeficiente de transferencia de calor por convección del problema en cuestión, tiene las
siguientes características:

 Convección: Libre
 Flujo: Externo
 Régimen: Laminar
 Sin cambio de fase
 Superficie caliente,
 Fluido: Newtoniano, Aire
 Superficie: Cilindro Circular
 Estado: Estacionario

Consideraciones Generales:

 No se tiene en cuente el efecto del extremo semiesférico del bulbo, ni tampoco la unión
entre el bulbo y la columna, además se asume la transferencia de calor sólo por la
superficie lateral del cilindro.
 Si el número de Rayleigh es inferior de alrededor de 109, el flujo de capa limite será
laminar [MILLS].

37
  El aire se asume como flujo incompresible, lo que implica que el número de Mach sea
menor de 0.3, es decir: . . Si se toma la velocidad estándar del sonido en aire de
340 m/s, / , entonces la máxima velocidad que podría tener el aire para ser
flujo incompresible, seria de:

. :

. . . ·

Luego la velocidad máxima permisible para tratar el aire como flujo incompresible es
á ⁄ . Como ya se había enunciado en el Capítulo 1, “las velocidades
asociadas a la convección natural suelen ser sumamente pequeñas, sin superar los 2
m/s” [MILLS], por lo tanto si el aire estuviera a esa velocidad estamos apenas en un
1.96% del límite para dejar de ser incompresible.

El flujo de convección libre sobre un cilindro vertical se

encuentra en muchos casos de aplicación tales como haces
de tubos en los reactores nucleares, resistencias cilíndricas,
etc. por lo tanto las características especiales de estos flujos
son de interés especial para la investigación.

Si se considera una superficie caliente rodeada de un fluido

frío, el flujo es impulsado por una fuerza de empuje debida a
los gradientes de densidad en la región adyacente a la
superficie caliente. Se parte del hecho de que la velocidad en la pared tiene un valor de cero y
que empieza a incrementarse alcanzando un máximo y luego tiende a cero al alejarnos de la
pared hacia el seno del fluido; este máximo ocurre dentro de la capa límite térmica. Este
comportamiento de la velocidad se explica porque al distanciarse de la superficie la viscosidad
disminuye (por ende la velocidad tiende a aumentar), y la diferencia de densidad se va
reduciendo (por ende la velocidad también tiende a disminuir); a la larga la fuerza de flotación
(que es la fuerza impulsora de la convección libre) disminuye conforme la densidad del fluido
tiende al valor del fluido circundante, lo que provoca que la velocidad primero alcance un
máximo y luego tienda a cero al alejarse de la superficie caliente.

38
 Figura 2.2 Perfil de la capa limite sobre un cilindro vertical caso de caliente (der.) y frío (izq.) [BEJAN]

Las velocidades locales ascendentes van creciendo al aumentar la distancia desde la parte
inferior, ya que el fluido que se aloja por debajo se encuentra estancado. Cerca a la superficie
las velocidades aumentan con mayor rapidez por el marcado efecto de la viscosidad, mientras
que la diferencia de densidades afecta con menos impacto. Lejos de la superficie la viscosidad
es mucho más grande, por ende la velocidad tiene a reducirse y además los gradientes de
temperatura también.

Los campos de temperatura son similares a los de convección forzada, como también lo son
las características de ambos tipos de capas limites. La influencia de la curvatura transversal se
incrementa con la distancia axial desde el fondo del cilindro (axial) [MARTINENKO].

39
 Figura 2.3 Interferogramas de convección libre externa en cilindros verticales [MARTINENKO]

2.3 MÉTODO DIFERENCIAL

Las ecuaciones diferenciales de conservación de materia (continuidad), cantidad de

movimiento y energía en coordenadas cilíndricas, pueden ser encontradas haciendo los
respectivos balances sobre un elemento diferencial de volumen, o simplemente retomadas de
una referencia [BIRD].

Continuidad:

Cantidad de Movimiento:

En función de los gradientes de velocidad para un fluido newtoniano con propiedades

constantes.

Componente r:

40
Componente :

Componente z:

Energía:

Para fluidos newtonianos en función de las propiedades de transporte constantes.

Consideraciones adicionales a las de ítem 2.2:

 Fluido Newtoniano.
 Se desprecia la disipación viscosa y de compresibilidad en la ecuación de energía.
 Flujo laminar (laminas delgadas) en dirección z.
 Propiedades constantes, excepto la densidad y viscosidad que cambian solo en
dirección z.
 Estado estacionario.
 La pared tiene mayor temperatura que el fluido, (superficie caliente).

41
  La componente radial de la velocidad (en r) es mas pequeña que la componente axial
(en z), , es decir, el fluido asciende casi verticalmente, por tanto se desprecia la
componente radial de la velocidad. La componente angular de la velocidad es cero.
.
 La temperatura no es función del ángulo.
 Las velocidades y son funciones exclusivas de r y z.

Figura 2.4 Capas laminares alrededor del cilindro vertical [MARTINENKO]

Por todo lo anterior, las ecuaciones de conservación se reducen a:

Continuidad:

Cantidad de Movimiento:

Componente z

El gradiente de presión de la ecuación es una derivada total ya que solo depende de z y la

gravedad es negativa porque va hacia abajo.

42
Lejos de la pared, los gradientes de velocidad son cero y esta ecuación se reduce a

Lo que permite encontrar el gradiente de presión:

Sustituyendo:

Dividiendo entre :

A partir de la definición del coeficiente de expansión, se realiza la aproximación de Boussinesq:

v
v

Considerando que el volumen específico cambia muy poco con la presión, ya que los líquidos
son incompresibles, y que en el caso de los gases, los cambios de presión son despreciables
porque la altura de trabajo es muy pequeña, y como también se había señalado anteriormente
al ser el número de Mach menor de 0.3 se asume flujo incompresible, entonces es válido
aproximar la densidad solo como función de la temperatura, por lo que es posible escribir:

v v
v v

Y sabiendo que ó v
v

se tiene que v

Al reemplazar:

43
 v v
v

Aproximando el diferencial:

∆
∆

Despejando:

Así la ecuación de cantidad de movimiento queda:

Energía:

Dividiendo por :

Introduciendo la difusividad térmica:

El cambio de temperatura interviene en la ecuación de movimiento y la distribución de

velocidad en la ecuación de energía, por tanto las ecuaciones están acopladas.

44
Hay que considerar que el transporte (tanto de energía como de cantidad de movimiento en la
dirección axial, z) debido a los mecanismos moleculares, es pequeño en comparación con el
convectivo, razón por la cual se considera despreciar los términos de segundo orden:

Así el sistema de tres ecuaciones a resolver es:

Continuidad:

[2-1]

Cantidad de Movimiento:

Componente z

[2-2]

Energía:

[2-3]

Estas 3 ecuaciones de conservación se resuelven sujetas a las siguientes condiciones de

frontera:

i. En la pared:
En se tiene: , ,
ó
, , , ,
ii. En el seno del fluido:

45
 En ∞ se tiene: , , ,

lim , lim , , lim

A continuación se explicarán varios métodos para resolver estas ecuaciones de conservación;

el desarrollo en detalle de la resoluciones se encuentra en el Apéndice A.

2.3.1 Resolución Por Semejanza o Análisis Dimensional

La resolución del sistema anterior, de 3 ecuaciones diferenciales y acopladas, de manera

analítica es muy difícil por su no linealidad.

Dado un cilindro vertical de altura, , y diámetro, . Se introducen las siguientes variables

adimensionales:

46
2.3.1.1 Ecuaciones y Condiciones Límite Adimensionales

Adimensionalizando las ecuaciones de conservación [2-1], [2-2] y [2-3]:

Las condiciones limite de contorno o frontera adimensionales de estas ecuaciones son:

 En la pared:
Para

⁄ ⁄ ⁄

para todo

se tiene

 En el seno del fluido:

En ∞
para todo
se tiene

47
2.3.1.2 Solución

El propósito principal al resolver las ecuaciones es obtener el flujo de calor transferido y para
ello se requiere el coeficiente de transferencia de calor a través del Número de Nusselt medio,
.

ó ó

⁄ ⁄

se tiene

⁄ ⁄
·
⁄ ⁄

En el sistema de ecuaciones de conservación se observa que en la ecuación de cantidad de

movimiento aparece el parámetro adimensional, número de Prandtl, , y en las condiciones
limite aparece el número de Rayleigh, , además se recuerda que ; como lo que
se requiere, en número de Nusselt, , contiene el número de Rayleigh, , por ende es
correcto establecer que:

2.3.2 Resolución por Análisis de Escala

Para el análisis de escala se requiere que todas la variables adimensionales queden

parametrizadas, lo que se logra eligiendo valores característicos que garanticen que todas
estas variables tengan un orden de magnitud de la unidad o menor. Este análisis indicara que
términos de las ecuaciones diferenciales no son indispensables, además que también indica
algo sobre la estructura y forma de la solución.

48
Se definen entonces las variables adimensionales:

; ; ; ;
∆ ∆⁄

Observe que la escala de surge de la ecuación de continuidad:

es la velocidad de escala y su magnitud apropiada se determina a partir de las ecuaciones.

Se usa el espesor de la capa limite térmica, ∆, en la escala de y en todas las ecuaciones, en
lugar de usar el espesor de la capa limite hidrodinámica, , para las ecuaciones de continuidad
y cantidad de movimiento. Para , y ∆ tienen de hecho una magnitud similar, puesto que
la fuerza de empuje produce un flujo en la capa límite térmica, si bien la velocidad máxima
ocurre cerca de la pared. Para , ∆ es mucho menor que , ya que las fuerzas viscosas
ejercen una fuerza de arrastre sobre el fluido que se encuentra fuera de la capa límite térmica,
aunque es la fuerza de empuje la que impulsa al flujo, y la velocidad máxima, , debe ocurrir
para un valor de ∆. Usando ∆ en las ecuaciones de continuidad y cantidad de movimiento
se garantiza que ⁄ ,y ⁄ sean del orden de la unidad.

Las ecuaciones de cantidad de movimiento y de energía a escala se obtienen sustituyendo las

respectivas equivalencias de las varibles adimensionales en ecuaciones de conservación [2-2]
y [2-3]:

Ecuación de Cantidad de Movimiento:

∆ ∆

Ecuación de Energía:

2.3.2.1 Número Prandtl Grande,

49
Este caso es el límite cuando Prandtl es grande; para estos fluidos, ∆ ; la transferencia de
calor está controlada por el interior de la capa límite, donde las fuerzas de inercia (aceleración)
son despreciables y por tanto las fuerzas de empuje son equilibradas por las fuerzas viscosas.

Por tanto en la ecuación adimensional de cantidad de movimiento se tiene:

~
∆

de donde:

~
∆

Y por la ecuación de energía:

~
∆

de donde:

∆ ~

Resolviendo las ecuaciones

~
∆

∆ ~

con variables y ∆.

Se llega a

∆
~ ⁄

∆
~ ⁄

2.3.2.2 Número Prandtl Pequeño,

50
Este caso es el límite cuando Prandtl es pequeño; la alta difusividad térmica de estos fluidos,
hace que ∆ y por tanto la transferencia de calor está controlada por el exterior de la capa
límite, donde las fuerzas viscosas son muy pequeñas y por tanto las fuerzas de empuje son
equilibradas por las fuerzas inerciales.

Así de la ecuación de cantidad de movimiento:

de donde:

⁄
~

Y relacionando con la de energía se obtiene:

∆
~ ⁄

∆
~ ⁄

Donde es el número de Boussinesq.

2.3.2.3 Número de Prandt Unitario .

En este caso el Número de Prantl tiende a uno, o tiene orden de la unidad. Para fluidos en los
que ∆ , tanto la fuerza inercial como viscosa tienen una magnitud comparable respecto a la
fuerza de empuje, la que está impulsando el flujo en todo momento.

De la ecuación de cantidad de movimiento,

~ ~
∆

se llega a:

∆
~ ⁄

51
 ∆
~ ⁄

2.3.2.4 Resultado

La magnitud del número de Prandtl determina si la convección libre se caracteriza por el

número de Rayleigh, el de Boussinesq o el de Grashof. Teniendo en cuenta que:

para el caso en que , los tres grupos adimensionales son equivalentes.

Determinando el Número de Nusselt

De acuerdo al balance de calor en la interfase sólido-fluido, :

ó ó

se tiene

de esta manera,

~
∆

Como se define el número de Nusselt:

⁄ ⁄ ⁄ ⁄ ⁄

⁄
~
⁄ ⁄ ⁄

o de otra manera:

⁄ ⁄
~

52
Donde la constante se encuentra experimental o teóricamente, mientras que el exponente
predicho teóricamente por el análisis de escala asume los siguientes valores:

Exponente

⁄ ⁄
Tabla 2.3 Valores del exponente para el número de Prandtl

Para el caso del aire como fluido de trabajo . , el análisis lleva a:

⁄
~

Estas expresiones sirven para correlacionar datos experimentales. Tanto el análisis

dimensional por el teorema Pi de Buckingham como el principio de semejanza solo conduce al
resultado , , mientras que el análisis de escala permite obtener la forma de esta
función.

2.3.3 Resolución por Teoría de la Capa Límite

 Sistema adimensional

Se suponen perfiles semejantes de velocidad y de temperatura

Si las dos capas límites se relacionan:

Δ

De la ecuación de continuidad:

53
 Figura 2.5 Adimensionalización de los perfiles [El Autor]

Reemplazando esto junto con los perfiles propuestos en la ecuación de movimiento [2-2]:

Integrando dentro de la capa limite, respecto a ,

Definiendo:

La ecuación de movimiento queda:

Análogamente la ecuación de energía [2-3] queda:

54
Ahora si
Δ

Se tiene que:

Δ Δ

Integrando dentro de la capa limite :

Δ Δ

Definiendo:

Así la ecuación de energía queda como:

E D

Se resuelve el sistema de ecuaciones diferenciales ordinarias de cantidad de movimiento y

energía considerando un perfil de capa límite:

E D

Combinando estas ecuaciones al dividirlas entre sí:

55
Una vez que se conoce Δ, se puede encontrar a partir de las ecuaciones anteriores y así
saber el perfil de temperatura con el cual, finalmente, se determina el coeficiente de
transferencia de calor.

2.4 MÉTODO INTEGRAL

El desarrollo más detallado y paso a paso de este método se puede observar en el Apéndice A
de este trabajo.

Figura 2.6 Esquema para tomar el elemento diferencial de volumen sobre la capa limite

2.4.1 Planteamiento

2.4.1.1 Ecuación de Cantidad de Movimiento

Aplicando la Segunda Ley de Newton:

56
 Figura 2.7 Balance de cantidad de movimiento sobre el elemento diferencial de volumen [El Autor].

Desarrollando lo correspondiente al cambio o flux de cantidad e movimiento; el término de

cantidad de movimiento que sale del elemento diferencial:

∆ ∆ ∆

de la misma manera, el flux de cantidad de movimiento que entra:

Ahora identificando todas las fuerzas externas que actúan sobre el volumen de control
elemental en dirección z para encontrar la fuerza neta:

Fuerza Viscosa en la Pared:

| ∆ ∆

Fuerza de Flotación:

57
 Peso:

Calculo de la fuerza neta:

, ∆

Entonces como ya se conoce cada uno de los términos aplicamos la segunda Ley de Newton:

Dividendo entre ∆ y tomando el límite cuando ∆ y por definición de derivada:

recordando que

queda:

Anteriormente habíamos establecido la aproximación de Bousinesq (con en todos los

términos excepto en los de la fuerza de empuje), dividiendo entre :

donde:

Para finalmente obtener la ecuación de conservación de cantidad de movimiento en su forma

integral:

58
[2-4]

2.4.1.2 Ecuación de Energía

Aplicando la Primera Ley de Termodinámica al volumen de control elemental:

h , h ,

Se considera que el flujo es estacionario, se desprecia la conducción en la dirección , y la

conducción en la dirección es cero por fuera de la capa limite. La transferencia de calor se da,
de la pared al volumen de control elemental, por conducción.

Donde se tiene que:

- el sistema está en estado estacionario:

- no hay intercambio de trabajo con los alrededores:

Existen dos corrientes de salida y una de entrada en las cuales despreciamos los cambios en
su energía potencial y cinética ya que el elemento es muy pequeño.

el flujo de calor en la pared hacia el fluido se da por conducción en la dirección radial:

∆ ∆

así la primera Ley se puede escribir como:

h , h ,

59
 h , h ,

Reemplazando se tiene que:

h | h | h | ∆

y tomando el límite cuando ∆ y por definición de derivada

h h

y recordando que

entonces dividiendo entre :

h h

a partir de:

como las propiedades se asumen constantes:

h h

reemplazando este resultado y dividiendo entre :

Por definición del coeficiente de difusividad térmica se obtiene la ecuación de conservación de

energía en su forma integral

[2-5]

60
2.4.2 Solución

Resolución simultánea de ecuaciones del sistema:

 Conservación Cantidad de Movimiento:

 Conservación Energía:

 Las condiciones de frontera para este caso son:

2.4.2.1 Introducción de perfiles de velocidad y temperatura propuestos

Se establece la relación constante (independiente de z) entre dos las capas límites

(hidrodinámica y térmica):

∆ ∆

61
 Figura 2.8. Relación entre los espesores de capas térmica e hidrodinámica y sus variaciones [El Autor]

Si se plantea que

Con base en los perfiles trabajados usualmente para una pared plana vertical [MILLS] y
también de acuerdo a la Figura 2.8, se proponen las siguientes formas para los perfiles de
velocidad y temperatura:

62
 Figura 2.9 Perfil de la capa limite hidrodinámica [El Autor]

Donde es una velocidad de escala. Estos perfiles están acordes con las condiciones de
frontera y además son suaves en el borde de las respectivas capas limites, es decir:

pero no se cumple que:

sin embargo, estos perfiles en paredes planas han resultado sorprendentemente buenos al
comparar con datos experimentales [MILLS].

En función de las nuevas variables adimensionales estos perfiles son:

Ahora se sustituyen estos perfiles en las ecuaciones de conservación en su forma integral.

63
La ecuación de cantidad de movimiento [2-4]:

queda:

∆ ∆ ∆ ∆

Y la ecuación de conservación de energía [2-5]:

queda:

∆ ∆ ∆ ∆ ∆ ∆
∆

∆ ∆ ∆ ∆ ∆ ∆
∆

El sistema de ecuaciones se reduce a:

Movimiento:

∆ ∆ ∆ ∆

Energía:

∆ ∆ ∆ ∆ ∆ ∆
∆

Suponiendo que los espesores de las capas térmica, , e hidrodinámica, , son próximos, es
decir, , se tiene entonces que ∆ , lo cual es valido para números de cercanos a 1
(caso aire). Por tanto si ∆ , el sistema queda:

64
Figura 2.10 Perfiles axiales de temperatura y velocidad y capas limites térmica e hidrodinámica [El Autor]

Dividiendo cada ecuación por :

[2-6]

[2-7]

65
2.4.2.2 Solución del sistema y

El sistema de ecuaciones diferenciales ordinarias de primer orden conformado por las

ecuaciones [2-6] y [2-7] tiene a y como funciones dependientes de y se resuelve por el
método de perturbaciones, donde estas dos variables se expanden en series de potencias
alrededor de ⁄ .

Esto se sustenta en el hecho de que dicho método es aplicable a cualquier ecuación diferencial
que contenga un parámetro de perturbación, , el cual también podría aparecer en las
condiciones de frontera. Dicho parámetro debe ser una constante positiva, tal que .

En este caso, dicho parámetro es el reciproco del radio del cilindro, , el cual se encuentra en
el sistema de ecuaciones y en la condición de frontera, luego:

que en las ecuaciones [2-6] y [2-7] se observa:

si hacemos que ∞, el sistema de ecuaciones [2-6] y [2-7] queda:

que son las ecuaciones que permiten encontrar a y para una placa vertical (KARLEKAR,
MILLS).

Se observa que si el parámetro de perturbación tiende a cero, , se anula el efecto de la

curvatura, lo que implica que, ∞. El mismo resultado se puede obtener asumiendo que el
espesor de la capa limite, , es pequeño comparado con el radio, , de tal manera que la
relación / también tiende a cero y los mismos términos del sistema de ecuaciones de
conservación se cancelan quedando el modelo de la pared plana vertical.

66
Entonces se resuelve nuestro sistema de ecuaciones diferenciales con el método de
perturbaciones en / , lo cual requiere expresar la solución de la siguiente manera:

Tomando solo los dos primeros términos , queda:

Esto es válido como una aceptable aproximación, ya que al aumentar el término siguiente va
a ser mucho más pequeño. Es de esperar que el primer término, en , sea de una pared
plana ya que carece del parámetro de curvatura, / .

Estos desarrollos de y en serie de potencias con dos términos exigen averiguar los
valores de , , . Para lograr dicho objetivo se procede a sustituir estas expresiones
en el sistema de ecuaciones de conservación [2-6] y [2-7]:

donde

Reemplazando las expansiones en la ecuación de cantidad de movimiento [2-6]:

67
Desarrollando el algebra y organizando de acuerdo a las potencias de / , para luego
factorizar por el grado de las potencias y finalmente como aproximación despreciar las
potencias de / , para las cuales . Igualando las respectivas potencias, se obtiene:

Simplificando,

Reemplazando ahora de la misma manera la serie de potencias en la ecuación de energía [2-

7]:

Se obtiene:

68
 :

Simplificando,

Agrupando, de acuerdo al orden, y , los resultados arrojados tanto por la ecuación de

Cantidad de Movimiento como por la de Energía, se tiene:

2.4.2.2.1 Resolución del sistema para

Tomando las ecuaciones correspondientes al parámetro de curvatura de grado cero:

Cantidad de Movimiento

[2-8]

Energía

[2-9]

69
Es claro que este sistema corresponde a una pared plana, que no contempla el efecto de la
curvatura; como en ese caso comúnmente se hace, se supone una solución de prueba que
corresponde a variaciones de la ley de potencias para y , de esta forma:

Cuyas derivadas respectivamente son:

Donde , , y son constantes a determinar.

Sustituyendo estas soluciones de prueba en el sistema de orden cero para el parámetro de

perturbación.

En la ecuación de movimiento [2-8]:

En la ecuación de energía [2-9]:

En estas ecuaciones la potencia de debe ser idéntica en cada término lo que nos lleva a tener
un sistema de ecuaciones en el cual podemos hallar los exponentes y :

que se resuelve satisfactoriamente para

Sustituyendo estos valores en las dos ecuaciones [2-8] y [2-9] para encontrar los coeficientes
y :

Momentum:

70
Energía:

Resolviendo el sistema:

Se obtiene:

En forma aproximada:

como se presenta en Karlekar.

Ahora buscando el valor de en

o aproximadamente, según Karlekar:

71
 .

De esta manera las expresiones para y son

y de la misma forma

2.4.2.2.2 Resolución del sistema para

Tomando las ecuaciones correspondientes al parámetro de curvatura de primer grado:

Cantidad de Movimiento

[2-10]

Energía

[2-11]

Es claro que la resolución de este sistema ya será una primera aproximación al hecho de que
el cilindro es superficie curva, de la misma manera se supone, como en el anterior caso, una
solución de prueba y en modelo potencial, teniendo ya conocidas y

72
Cuyas derivadas respectivamente son:

Donde , , y son constantes a determinar.

Resolviendo el sistema de potencia cero se encuentra que:

Derivando:

Sustituyendo estas soluciones de prueba en el sistema de primer orden para el parámetro de

perturbación, conformado por las ecuaciones [2-10] y [2-11]:

En la ecuación de movimiento [2-10]:

En la ecuación de energía [2-11]:

La dependencia de se elimina en estas ecuaciones si los exponentes y cumplen que:

73
que se resuelve satisfactoriamente para

Sustituyendo estos valores en las dos ecuaciones [2-10] y [2-11] para encontrar los coeficientes
y :

Momentum:

Energía:

Resolviendo el sistema con incógnitas y :

Se tiene la solución

2.4.2.3 Número de Nusselt

Una vez encontrados todos los parámetros de los perfiles de velocidad y temperatura
propuestos entonces se procede a hallar el número de Nusselt.

74
2.4.2.3.1 Local

Un balance de calor en la interfase, , proporciona:

ó ó

. .

Despejando :

Del perfil parabólico de temperatura asumido, se tiene:

Como se asumió que ∆ y despejando para derivar parcialmente respecto a la coordenada

radial:

Evaluando en :

Entonces el coeficiente de transferencia de calor queda:

75
Se había aproximado que:

Por tanto el coeficiente local de transferencia de calor será:

con los valores de

desarrollando el algebra se llega a:

[2-12]

√
√

Por tanto el Nusselt local será

√
√

2.4.2.3.2 Promedio

Para hallar ahora el Nusselt promedio se requiere encontrar el coeficiente promedio de

transferencia de calor por el teorema integral del valor medio, de esta forma:

76
Se integra por sustitución y se obtiene:

Si

Se llega a:

[2-13]

donde

El Nusselt promedio será entonces:

[2-14]

77
una expresión más sencilla consiste en tomar en el valor promedio del coeficiente de
transferencia [2-13], una aproximación para el logaritmo natural:

ln

Para | |

Bajo este criterio se puede expresar que:

Aproximando a los cuatro primeros términos:

Reemplazando en la expresión del coeficiente promedio:

[2-15]

Hallando el Nusselt promedio:

78
Si se expresa en función del diámetro, :

[2-16]

Ya que . Esta solución fue obtenida por Le Fevre y Ede, donde el primer término
representa la solución de una pared plana vertical, mientras que del segundo en adelante, los
términos contienen el efecto de la curvatura, ya que llevan el radio o el diámetro, y por ende su
relación con la altura. Se puede observar matemáticamente que el coeficiente y el Nusselt de
un cilindro es mayor que el de una pared plana, gracias al efecto de la curvatura.

79
 3. CONSTRUCCIÓN Y SOLUCIÓN DE LOS MODELOS MATEMÁTICOS

El problema fundamental que se formuló en esta tesis fue determinar experimentalmente el

coeficiente de transferencia de calor alrededor del bulbo de un termómetro. Para cumplir con
este objetivo básico es necesario modelar el funcionamiento del termómetro de vidrio para
generar una ecuación, que pueda utilizarse para calcular el coeficiente de transferencia de
calor. Para poder confrontar los resultados del modelo del termómetro, es necesario comparar
los resultados experimentales emanados por el mismo modelo con los cálculos del coeficiente
de transferencia de calor producidos por los modelos teóricos basados en la teoría de la capa
límite. Por ello el presente capítulo se divide en dos partes: construcción del modelo de
funcionamiento del termómetro de vidrio y construcción de modelos para predecir el coeficiente
de transferencia de calor.

Antes de modelar un fenómeno físico es necesario hacer una descripción fenomenológica de la

situación que se quiere estudiar. Estos modelos deben estar disponibles a recibir la información
experimental para procesarla y producir como salida el coeficiente de transferencia de calor. La
materia prima necesaria para esto son las curvas temporales de temperatura, las propiedades
de los fluidos (tanto el termométrico como el fluido que está rodeando al bulbo), y los datos
geométricos del mismo bulbo.

3.1 FUNCIONAMIENTO DEL TERMÓMETRO DE VIDRIO: FLUIDO TERMOMÉTRICO

Una estrategia para el modelamiento es: Primero describir fenomenológicamente la situación

física que se quiere modelar. Posteriormente con base en lo que se observe en la descripción
fenomenológica, se plantean modelos teniendo en cuenta restricciones impuestas. Los
modelos más sencillos son los que tienen más restricciones y los más complejos los que tiene
menos restricciones.

A continuación se hace una descripción fenomenológica de la situación que se quiere estudiar, en

este caso la descripción del funcionamiento del termómetro de vidrio.

3.1.1 Descripción Fenomenológica

80
El funcionamiento de los termómetros de vidrio con líquido se basa en la dilatación térmica de un
líquido termométrico introducido en un recipiente de vidrio de paredes delgadas (recipiente
termométrico el cual es el bulbo) comunicado con un capilar de paredes gruesas (capilar del
termómetro). Al capilar se le ha asignado una escala. El líquido termométrico que se encuentra en el
capilar se denomina columna y su menisco (menisco de la columna) sirve como marca de lectura.
Existen varios líquidos que pueden servir como fluidos termométricos, algunos de ellos son:
Mercurio, mercurio-Talio, Toluol, Alcohol, Mezcla de Pentano.

Para el propósito del presente trabajo, el termómetro se considera como un sistema que consta
de un subsistema llamado bulbo o extremo expuesto al fluido, estando dentro de una posición
determinada dentro de un fluido o ambiente, en la cual se mide la temperatura local, y además
del otro subsistema delimitado por el capilar de vidrio que se denomina columna, dentro de la
cual se mueve el fluido termométrico bien sea dilatándose o contrayéndose.

La transferencia de calor desde el bulbo hacia el medio que lo rodea se da primero por
conducción-convección dentro del fluido termométrico seguida por la transferencia de calor
conductiva a través de su pared y finalmente convección-radiación de la superficie exterior del
bulbo hacia el medio.

Cuando el bulbo está recibiendo o perdiendo calor durante el funcionamiento del termómetro en
la columna hay una diferencia de temperatura entre el extremo de la columna cercana al bulbo
y el otro extremo de la columna; entonces en términos generales hay transferencia de calor
conductiva y convectiva a lo largo de la columna. Si la columna de fluido termométrico y la
envoltura de la columna de vidrio tienen diferentes temperaturas existirá transferencia de calor
a través de la envoltura y entre ésta y el medio externo.

Los fenómenos de transferencia que se identifican son:

 Transporte de calor conductivo (dentro del bulbo y dentro de la columna) y flujo

convectivo (dentro y fuera del bulbo), se desprecia la radiación desde la superficie del
bulbo hacia el medio.
 Transporte de cantidad de movimiento en el capilar del termómetro.

81
Una vez realizada la descripción fenomenológica se identifican todos los subsistemas que
conforman el sistema para proceder a realizar el modelamiento del funcionamiento del
termómetro. Entonces de la descripción fenomenológica anterior puede identificarse, que el
termómetro consta de al menos tres subsistemas:

 El bulbo del termómetro (que es la frontera por la que se transfiere calor).

 La columna de fluido (que se asume adiabática).
 La envoltura de vidrio del fluido termométrico.

3.1.2 Construcción

El sistema está formado por varios subsistemas. Como primer modelo, el más sencillo,
considera el sistema formado por dos subsistemas: el bulbo y la columna de fluido
termométrico.

Bulbo, b Columna, c

Figura 3.1 Subsistemas básicos del termómetro

Según lo anterior el termómetro puede idealizarse, como formado por una columna de mercurio
cilíndrica unida al cilindro de mercurio con mayor diámetro que el de la misma columna. Los
dos cilindros de mercurio están inicialmente a una temperatura ambiente, y a partir de un
instante dado, se calienta o se enfría el cilindro que sirve de bulbo exponiéndolo a un ambiente
convectivo de temperatura constante, , (temperatura del medio a medir). Para modelar
va a suponerse que el bulbo gana calor (el balance diferencial de energía da la misma ecuación
bien sea calentando o enfriando), debido al calentamiento el fluido termométrico se expande.
Se hacen las siguientes restricciones adicionales:

 En el fluido termométrico sólo se considera flujo conductivo de calor y el calor se

propaga instantáneamente (difusividad térmica infinita) y no se considera flujo
convectivo dentro del fluido termométrico. Esto equivale a que la columna tenga la
misma temperatura del bulbo que en la realidad correspondería a un termómetro de
inmersión total.

82
  Con respecto a lo anterior puede decirse que para el bulbo se aplica el método de
resistencia interna despreciable, que es válido para números de Biot menores a 0.1. Si
el mercurio a 25 °C tiene una conductividad térmica aproximada de 8.81 W/m°C, con un
coeficiente de transferencia de calor por convección natural en aire de alrededor de 5
W/m2°C, y para un bulbo del termómetro asumido con un cilindro con una longitud
característica (Relación Volumen-Área, V/A, ver en características del termómetro,
Tabla 4.1) de 1.85 x 10-3 m, se tiene que el numero de Biot estimado es:
· .
. .
.
 La transferencia de calor entre el medio externo y el mercurio sólo se da a través del
bulbo del termómetro, mientras que la columna de vidrio que contiene el resto de
mercurio se supone que no transfiere calor desde o hacia el ambiente externo que lo
rodea, es decir, el capilar de vidrio se considera adiabático, perfectamente aislado.
 Se supone que el espesor de la pared de vidrio del bulbo del termómetro es tal que la
resistencia conductiva de calor en la pared de vidrio es pequeña, ello debido a la poca
inercia que tiene el termómetro, es decir, el menisco del mercurio se mueve tan pronto
como la temperatura cambia en el bulbo. Esto se verifica calculando el número de Biot
con las especificaciones de la pared de vidrio del bulbo dadas en la Tabla 5.2 y
tomando un valor de coeficiente de transferencia de calor por convección libre en el aire
alrededor de 5 W/m2°C; se encuentra que este valor está entre 0.0017 y 0.0054,
mientras que para un valor de 7 W/m2°C el máximo valor alcanzado es de 0.012. Lo
que indica que como el Número de Biot es menor de 0.1 la resistencia interna es
despreciable y esta consideración es válida para el caso del aire.
 El fluido termométrico no moja la pared, despreciándose los efectos de tensión
superficial.
 Se ignoran los efectos debidos a las fuerzas capilares en la columna termométrica.
 No hay condensación de los vapores del fluido termométrico en las paredes de la
columna capilar del termómetro.
 Como consideración final, en este modelo va a suponerse que el fluido en el
termómetro se expande a presión constante dentro del capilar, lo cual implica suponer
que se trabaja con el calor específico a presión constante.

El desarrollo paso a paso de la construcción del modelo del fluido termométrico se puede
observar en el Apéndice B de este libro.

83
Por las relaciones termodinámicas, se puede lograr una relación entre el calor específico a
presión constante, , y el coeficiente de expansión térmica, , del fluido termométrico:

Dado que

h v

Derivando parcialmente respecto a a constante:

h v

Se define que

v v
v
v

Luego:

Despejando:

Restricciones adicionales: va a suponerse que el valor del calor específico, , y el valor del
coeficiente de expansión térmica, , son constantes en el intervalo de tiempo de trabajo del
termómetro.

El considerar como fluido termométrico el mercurio, el cual es un metal líquido, permite

aproximadamente considerar que la temperatura se distribuye instantáneamente en él, es decir
no existe distribución espacial de temperatura tanto en el bulbo como en el columna del
termómetro (ver resistencia interna despreciable).

Como anteriormente ya se mencionó el hecho de que al expandirse isobáricamente todo el

fluido termométrico y la implicación derivada por ello de asumir el calor específico a presión
constante, , y por ende también se debe asumir el coeficiente de expansión térmica, .

84
Para modelar el termómetro puede iniciarse por aplicar el principio de conservación de la
energía, el cual se puede expresar en términos de un balance de energía sobre un volumen de
control correspondiente al sistema a modelar.

Según los fundamentos de la mecánica de fluidos, el balance de energía aplicado a un volumen

macroscópico de control [FOX]

que es la primera ley de la termodinámica, donde se representan los siguientes términos:

ó : í
í í á

: ó

: ́

3.1.2.1 Termómetro

El balance de energía global o en todo el termómetro se desarrolla de la siguiente manera. Se

selecciona un volumen de control, que en este caso es el volumen limitado por el contorno del
sistema (paredes del termómetro) y la superficie libre de mercurio. En este modelo es el flujo
de calor convectivo que atraviesa la superficie del bulbo y el término que en este caso
representa:

85
Como no hay trabajo de eje, es decir trabajo debido a piezas que giran respecto a un eje, ,
se reduce a:

ó ó
́

Figura 3.2. Termómetro ASTM 56C-86 usado en la experimentación [El Autor]

También puede ignorarse el trabajo debido a fuerzas cortantes ejercidas por el movimiento del
fluido sobre las fronteras del volumen de control, por lo cual . De los términos de
la integral de volumen que en conjunto representan la energía total dentro del volumen de
control, el término de energía potencial es insignificante y el término de flujo de trabajo debido a
fuerzas de presión internas al volumen de control es pequeño, pues los cambios de presión
dentro del volumen de control no son grandes. Como no hay flujo de masa debido al
movimiento del fluido que atraviese la superficie de control, entonces la integral de superficie en
la ecuación es cero. Según las restricciones anteriores, la ecuación del balance de energía se
reduce a:

· ·

donde

86
 :

: ó

: í ó

: ó ́

: ó

:

La ecuación anterior se interpreta como: El flujo convectivo de calor suministrado al bulbo del
termómetro es igual al cambio de la energía interna en el bulbo y la columna más el cambio de
la energía cinética en el bulbo y la columna al moverse debido a la dilatación del fluido, más la
potencia del trabajo de expansión del volumen de control contra la presión sobre la superficie
libre de mercurio.

En la ecuación se considera constante. Desde el punto de vista termodinámico significa que

cuando el calor específico se determina a presión constante, una parte de la energía
suministrada al sistema se convierte en energía interna y la otra se utiliza para hacer el trabajo
de expansión de las fronteras del volumen de control [WEBER], mientras que cuando se
determina el calor específico a volumen constante toda la energía suministrada se convierte en
energía interna. Usualmente los términos caloríficos involucrados son mucho mayores que los
términos de energía cinética por lo cual el balance global conduce a:

[3-1]

El balance global expresado en términos de números adimensionales es:

87
[3-2]

Donde

: ́

: ́ ,

: é ó

: é ó

Y los números adimensionales resultantes junto con su significado:

ó
í

El volumen en el termómetro es la suma del volumen en el bulbo más el volumen en la

columna. Este balance de volúmenes expresado en términos de volúmenes fraccionales y
longitudes fraccionales da:

donde los volúmenes fraccionales son:

, y

: ó

88
 : ó

: ó

Y las longitudes fraccionales son:

: , ó

: , ó

Para expresar la velocidad fraccional que figura en el balance global en términos de velocidad
se realiza lo siguiente. Se deriva respecto al tiempo el volumen fraccional , y expresándolo en
términos del tiempo adimensional conduce a

donde y el balance global de energía en términos de la velocidad fraccional es

El fluido termométrico en el termómetro sufre dilatación y contracción. La relación entre la

temperatura y la expansión o contracción de un fluido puede determinarse mediante algunas
relaciones termodinámicas según la literatura [VAN NESS]:

Definiendo el volumen específico, v, como: v= por termodinámica se tiene:

v=v ,

la derivada total será:

v v
dv=

dividiendo entre v

89
 v v v
=
v v v

Se tiene que el coeficiente de dilatación o expansión térmica, , se define como:

v
v

Y el coeficiente de compresibilidad, , como:

v
v

La expresión de volumen específico queda:

v
=
v
v
≈
v

Si se trabaja con expansión a presión constante, , por lo cual la ecuación anterior

conduce a:

v≈v

Esta última ecuación puede expresarse en términos del volumen global del fluido
termométrico y de la masa total en el termómetro como

Lo cual conduce a

La ecuación anterior relaciona el volumen del fluido termométrico con la temperatura y puede
expresarse en términos del volumen fraccional como

Κ Κ Κ
≈

Donde Κ es el coeficiente de dilatación térmica adimensional, definido como Κ .

Desarrollando en series de Taylor el miembro derecho de la última ecuación

90
 ∞
Κ
Κ
≈
n!

Κ
Κ Κ Κ
≈ Κ
! ! i!

Según los manuales de las casas fabricantes de termómetros [HAUSER] un buen termómetro
de fluido líquido termométrico debe satisfacer que el volumen del fluido sea proporcional a los
cambios de temperatura con respecto a una tempertura inicial. Por lo que en la serie de Taylor
solo se tiene en cuenta el término lineal y se descartan los términos no lineales, ello conduce a:

Κ
≈ Κ Κ

≈ Κ

Con esta ecuación puede determinarse la velocidad con que cambia el volumen de líquido en el
termómetro y también la velocidad de desplazamiento de la superficie libre del mercurio. Para
ello es necesario derivar respecto al tiempo esta ecuación y a continuación expresarla en
términos del tiempo adimensional,

≈ Κ

Pero anteriormente se había encontrado

Combinando estas dos ecuaciones la velocidad fraccional queda expresada en términos de la

temperatura, como:

Sustituyendo esta velocidad en términos de la temperatura en la ecuación de balance [3-2],

finalmente el balance global en términos de la temperatura adimensional es

[3-3]

91
Este modelo está compuesto por dos subsistemas: la columna termométrica y el bulbo. A
continuación va a efectuarse un balance de energía en cada uno de ellos.

3.1.2.2 Columna

Para este subsistema puede seleccionarse como volumen de control, el volumen del capilar en
el cual la única frontera que es atravesada por fluido, es que une el bulbo con el capilar.
Teniendo en cuenta las restricciones iniciales junto con la consideración de expansión del fluido
a presión constante en la columna, y al aplicar el balance de energía macroscópico, este se
reduce a:

Figura 3.3 Columna (izq.) y capilar (der.) del termómetro usado [El Autor]

El balance de la energía para la columna en términos de números y variables adimensionales

es:

[3-4]

donde los números adimensionales y son y significan

í
í

92
 í

3.1.2.3 Bulbo

Figura 3.4 Bulbo del termómetro usado [El Autor]

Aplicando el balance macroscópico de energía y tomando como volumen de control el volumen

del bulbo del termómetro, la ecuación de balance de energía para el bulbo queda:

Donde

Desarrollando y simplificando la ecuación anterior se obtiene finalmente el balance de energía

en el bulbo del termómetro:

La ecuación del bulbo en términos de números adimensionales es:

[3-5]

Los números adimensionales son

donde

el significado del número adimensional es

93
 í
í

Expresando el balance de energía de cada subsistema en la forma:

Columna:
Bulbo:

Sumando estas dos ecuaciones término a término se obtiene

Donde la suma de los flujos fraccionales de acumulación es

Pero de donde por lo que la ecuación general del termómetro puede

reducirse a

Que corresponde efectivamente al balance global en el termómetro.

Comprobado el balance global en el termómetro, y expresando en términos de la temperatura

adimensional queda:

3.1.3 Solución

Para resolver el balance global de energía obtenido [3-3] y obtener la variación de la

temperatura del fluido termométrico, se asume que los parámetros Κ, y de esta ecuación
son constantes, ello conduce a:

[3-6]

94
o en forma exponencial:

θ
θ

Para hacer estimativos de algunos términos de la ecuación anterior, en la literatura [KREITH]

se consultaron los siguientes datos para el mercurio a la temperatura promedio entre
° , y ° , es decir ° : . ° ; / ; ⁄ ° .
Teniendo en cuenta que el bulbo del termómetro usado tiene un diámetro interior aproximado
de 4.73 mm y una longitud de 47.40 mm, se desprecia además la masa y volumen de la
columna al compararla con la del bulbo. Asumiendo además presión atmosférica reinante
dentro de la superficie libre de mercurio y con los datos anteriores se obtienen el valor de los
términos:

.
· ° °
°

Κ . ° ° ° .

Que permiten calcular el término

Κ . . .

De donde puede observarse que el número adimensional Κ , lo cual significa que el flujo
de trabajo fraccional (número ) es muy pequeño comparado con los términos calóricos, por
tanto la respectiva expresión se reduce a:

Y de esta manera la ecuación de balance queda finalmente como:

[3-7]

Se recuerda que esta expresión se obtuvo para el caso en el cual el fluido termométrico se
dilata por calentamiento. Si se repite todo el procedimiento aplicado pero para el caso en el que
la columna termométrica se contrae por enfriamiento en el bulbo del termómetro, se obtiene la
misma ecuación, lo cual indica que el balance de energía es invariante de la dirección de flujo
de calor al entrar o salir del termómetro.

95
Transformando el modelo logrado a la forma dimensional:

Si :

de manera que

[3-8]

Otra manera más simplificada de lograr el mismo modelo se explica a continuación.

Un balance de energía sobre el sistema cerrado y en estado permanente correspondiente al

fluido termométrico:

Anteriormente se había establecido que el modelo de dilatación del fluido termométrico se

aproximaba a:

Derivando esta expresión respecto al tiempo con constante:

96
 ≈

Reemplazado en el balance de energía:

introduciendo las variables adimensionales:

O en forma más general:

Desarrollamos anteriormente el balance hasta este punto:

despejando :

donde todos los parámetros son medibles, lo que exige desarrollar:

97
3.2 CONVECCIÓN LIBRE ALREDEDOR DEL BULBO DEL TERMÓMETRO

Parte de los objetivos del presente trabajo es generar algunas ecuaciones para predecir
teóricamente el coeficiente de transferencia de calor. En el presente numeral, con base en los
principios y fundamentos de fenómenos de transporte, se plantearan y solucionarán algunos
modelos teóricos para calcular el coeficiente de Transferencia de calor para convección libre
alrededor del bulbo de un termómetro de vidrio. La forma expedita y explicita de construcción y
solución de los tres siguientes modelos se puede encontrar en el Apéndice C del presente
trabajo.

Situación física a modelar: El problema o situación física de la cual quiere predecirse el

coeficiente de transferencia de calor corresponde básicamente a los dos experimentos que se
enuncian a continuación:

Primer experimento: Se tiene un termómetro de vidrio colgado y fijo verticalmente. A partir de

un determinado instante se calienta el bulbo del termómetro con un mechero hasta que alcance
la temperatura cercana a la máxima de la escala del termómetro. Seguidamente se suspende
el calentamiento. Entonces el bulbo del termómetro comienza a enfriarse por convección libre
en el ambiente que lo rodea.

Segundo experimento: Se tiene un termómetro de vidrio situado horizontalmente y mantenido

fijo. A partir de un determinado instante se calienta el bulbo del termómetro con un mechero
hasta que alcance la temperatura cercana a la máxima de la escala del termómetro.
Seguidamente se suspende el calentamiento. Entonces el bulbo del termómetro comienza a
enfriarse por convección libre en el ambiente que lo rodea.

3.2.1 Termómetro En Posición Vertical Con Flujo Transitorio

3.2.1.1 Descripción Fenomenológica

El fenómeno de convección natural del fluido alrededor del bulbo para las configuraciones
vertical fue ampliamente discutido en la sección 2.2.

98
Aquí se parte del hecho que despreciando la resistencia térmica del vidrio, el modelo transitorio
del fluido termométrico debe coincidir con el modelo transitorio del fluido alrededor del bulbo en
un valor radial igual al radio del bulbo.

Figura 3.5 Relación instantánea entre los perfiles de temperatura en la interfaz [El Autor]

Al observar el bulbo del termómetro, en la cual sucede el transporte de calor, puede verse que
el bulbo geométricamente está compuesto por dos partes: el cuerpo del bulbo de forma
cilíndrica y el extremo del bulbo conformado por una semiesfera. Para efectos de modelamiento
la estrategia es inicialmente simplificar el problema asumiendo algunas restricciones. Y a
medida que se plantea el modelo van asumiéndose otras restricciones.

Restricción geométrica: Asumir que el bulbo del termómetro es básicamente un cilindro.

Suposiciones:

 Flujo transitorio de transporte de calor y de cantidad de movimiento.

 Velocidad radial : No existe desplazamiento lateral de materia. A partir de la
observación de fotografías del campo de flujo de velocidad en cilindros calentados,
Figuras 1.9 y 2.3 se ve que básicamente el flujo solo tiene componente axial de
velocidad.
 No existe movimiento angular .

 Gradiente de velocidad axial : Ello en virtud de las pequeñas dimensiones del

bulbo, el recorrido que hace el aire en la longitud de bulbo es pequeño.

 Gradiente de velocidad radial : Es el más importante.

99
  Se tienen en cuenta las fuerzas de flotación.

Considerando que la velocidad del aire es pequeña alrededor del bulbo, según la literatura para
números de Mach menores a 0.3 el gas puede tratarse como fluido incompresible y para efecto
en la ecuación de energía se considera propiedades medias constantes tales como la densidad
y la viscosidad. Bajo esta restricción y añadiendo la restricción de suponer que no existe
gradiente axial de temperatura además de las restricciones dadas anteriormente, se procede a
la construcción del modelo.

3.2.1.2 Construcción

Va a utilizarse el método diferencial. Partiendo de las ecuaciones diferenciales de conservación

de movimiento y energía y bajo las consideraciones establecidas anteriormente y en 2.2,
exceptuando la de estado estacionario:

Cantidad de Movimiento:

Componente z

Energía:

Si el componente radial de la velocidad es insignificante comparado con el axial y además

considerando como aproximación adicional que , , lo que puede inducir a despreciar los
cambios de velocidad y temperatura axiales o, dicho de otra manera, estos son muy pequeños
comparados con los cambios radiales, se pueden despreciar los términos que contengan los
factores:

100
Lo que permite reescribir las ecuaciones de la siguiente forma:

Cantidad de Movimiento:

Componente z

Energía:

Que se pueden escribir como:

Cantidad de Movimiento:

Componente z

[3-9]

Energía:

[3-10]

Aquí se puede resolver la ecuación de energía desacoplada a la de cantidad de movimiento y

con esta solución bajo la condición inicial y de frontera se procede a resolver la ecuación de
movimiento o flujo.

Resolviendo entonces la ecuación [3-10]:

101
Bajo las condiciones limites:

Inicial:
,
Frontera
,
en la superficie del bulbo ó ó

∞, ∞

3.2.1.3 Solución

El desarrollo en detalle de estas ecuaciones se puede encontrar en el Apéndice C. Por el

método de separación de variables planteando una solución de la forma:

hallando las derivadas parciales respectivas y reemplazando en la ecuación de energía:

′′ ′ ′

Dividiendo entre :
′
′′ ′

Si esto es igual al parámetro, :

′
′′ ′

Se desprenden las dos ecuaciones diferenciales ordinarias:

′

102
 ′′ ′

Resolviendo cada ecuación independientemente,

empezando con la de :
′

Despejando :

ahora se prosigue con la segunda ecuación de :

′′ ′

′′ ′

A partir de la sustitución:

con lo cual se cambia de variable independiente en cada uno de los términos de la ecuación.
Sustituyendo en la ecuación y dividiendo por :

La cual es una ecuación de Bessel con parámetro, , siendo un número real no negativo.

La forma estándar de dicha ecuación es [KREYSZING]:

′′ ′

Si , cuando toma el valor del entero cero, la ecuación se reduce a:

′′ ′

Que tiene como solución:

Donde es la función de Bessel de primera clase de orden cero y es la función de

Bessel de segunda clase de orden cero.

103
Para el caso de estudio la solución es:

entonces:

por tanto el perfil espacio-temporal de temperatura:

,
,

Llamando las nuevas constantes como

y

queda:
,

aplicando las condiciones iniciales y de frontera para hallar las constantes , :

Espaciales (de frontera):

∞, ∞

∞ ,

Temporales (iniciales):
,
∞

104
Resolviendo el sistema para hallar y :

Por regla de Cramer:

Por tanto la solución queda:

, ∞

Sustituyendo las constantes en la ecuación se obtiene:

[3-11]

∞
∞

Si , , en (superficie del bulbo) la evolución en el tiempo es

[3-12]
∞
∞

o en su forma logarítmica:

∞
∞

o despejando:

∞ ∞

105
La ecuación anterior, representa la evolución en el tiempo de la temperatura en la superficie del
bulbo, obtenida a partir del perfil de temperatura del fluido que rodea al bulbo. Ecuación que
tiene la misma forma exponencial de la evolución de temperatura que se obtuvo anteriormente
en 3.1 (en Fluido Termométrico) dentro del bulbo [3-8]:

∞
0 ∞

3.2.1.4 Parámetro

3.2.1.4.1 Despreciando la resistencia de pared ,y

Por comparación entre los dos perfiles [3-8] y [3-12]:

∞ ∞
∞ 0 ∞

Se obtiene que:

o de otra manera:

′ ′

Donde es la difusividad térmica del fluido termométrico. No olvidar que el número de Fourier,
, está dado por:
′

3.2.1.4.2 Considerando la resistencia de la pared

El balance de calor establece:

ó ó

106
 ∞
∞

∞
∞

Igualmente, al inicio:
∞
∞

Dividiendo estas dos últimas expresiones, extrayendo logaritmos, aplicando propiedades y

haciendo sustituciones se llega a:

∞ ∞
∞ ∞

Recordando los perfiles [3-8] y [3-12] para sustituir en los dos primeros términos logarítmicos
de la ecuación anterior:

donde es un coeficiente global de transferencia de calor (en realidad ), ya que tiene

incluida la resistencia de pared:

Despejando:

3.2.1.5 Relación entre y

Por Fourier y Ley de Omh de Calor se tiene para el coeficiente global :

∆
∆

De donde

107
Despejando :

Como:
∆ ∞

Las resistencias presentes entre estas dos temperaturas:

Figura 3.6 Circuito de Resistencias Térmicas [El Autor]

Se puede evaluar la resistencia térmica total equivalente del circuito anterior,

no olvidar que la resistencia de pared, , es conductiva en un cilindro hueco:

luego el coeficiente global será:

también:

o también despejando , que es el coeficiente individual (que no contiene la resistencia

conductiva de la pared) es:

108
[3-13]

Con todo lo anterior se puede reescribir la expresión de obtenida inmediatamente atrás en

3.2.1.4.2:

[3-14]
′

Se observa que si , se obtiene:

Esto puede suceder si se desprecia ante y y si además .

3.2.1.6 Relación con el perfil de velocidad

La constante surge durante el proceso de solución de la ecuación de energía, y ésta

depende de la velocidad del fluido, como se observara al sustituir el perfil de temperatura en la
ecuación de cantidad de movimiento. El parámetro puede obtenerse de los datos
experimentales de temperatura en el bulbo y por mínimos cuadrados por la ecuación de la
evolución de la temperatura.

De los datos de temperatura del termómetro puede determinarse y conocido se tiene el

perfil de temperaturas con el cual se resuelve la ecuación de cantidad de movimiento [3-9]
para determinar el perfil de velocidad:

109
 ∞

Condición inicial

Condición en la superficie

Retomando el perfil de temperatura [3-11]:

∞
∞

∞ ∞

El parámetro queda vinculado a la velocidad cuando la solución anterior se reemplace en la

ecuación de cantidad de movimiento .

La solución de esta ecuación permite relacionar la velocidad con el parámetro .

3.2.1.7 Número de Nusselt

Aplicando la condición en , ∞ al perfil de temperatura [3-11] se obtiene:

Finalmente el Nusselt para todo el sistema es:

[3-15]

Obsérvese que así el perfil de temperatura varia en el tiempo, la solución produce un solo
número de Nusselt para el sistema que no depende del tiempo, si el parámetro no depende
del tiempo, lo que implica que la resistencia de pared es despreciable y que el coeficiente de
transferencia de calor por convección es constante a través del tiempo, lo que en realidad no
sucede, pero que es salvable en la medida de que se podría asumir el valor promedio como

110
constante. Como no se conoce el espesor de la capa límite, , y dado que se obtuvo
experimentalmente un valor de y de tanto instantáneo como promedio, con estos
resultados y a través de la ecuación anterior es posible encontrar dicho espesor tanto
instantáneo como promedio.

3.2.2 Termómetro En Posición Vertical Con Flujo Cuasiestacionario

3.2.2.1 Planteamiento

Durante los experimentos se observó que la temperatura en la escala del termómetro a partir
de terminado instante varía muy lentamente. Definiendo una escala de tiempo adecuada puede
considerarse en este intervalo de tiempo que la derivada de la temperatura y velocidad con
respecto al tiempo no cambian apreciablemente y el problema puede tratarse como en estado

estacionario . Manteniendo las restricciones del modelo anterior más la

consideración de Estado estacionario las ecuaciones de cantidad de movimiento y de energía

se reducen a:

3.2.2.2 Solución

Las condiciones para la ecuación de energía son

En

Considerando que , , para que esta ecuación de energía sea diferencial ordinaria y
multiplicando dicha ecuación por :

111
Lo cual es una ecuación de Cauchy-Euler [ZILL] de la forma:

Cuya solución está dada de la forma

Aplicando a la ecuación de energía, la solución será:

Derivando

y sustituyendo en la ecuación de energía:

de orden de multiplicidad 2.

Luego la solución es de la forma:

Bajo las condiciones de frontera:

En :

En

112
Despejando :

Luego para :

Reemplazando en el perfil:

La solución de la ecuación de energía bajo estas condiciones es también de la forma:

[3-16]

De la segunda condición se puede definir el número de Nusselt, definiendo como longitud

característica el diámetro del bulbo del termómetro, aunque se está en libertad de elegir la
longitud del cilindro como longitud característica.

3.2.2.3 Ecuación cantidad de movimiento

3.2.2.3.1 Relación con el perfil de temperaturas

Se procede a sustituir el perfil de temperatura obtenido en la ecuación de cantidad de

movimiento:

113
Para obtener:

Donde y son:

Condiciones límite:
En
En

3.2.2.3.2 Solución

Solucionando la ecuación de cantidad de movimiento junto con las condiciones límite:

Es también una ecuación no homogénea de Cauchy-Euler análoga a la de energía. Por lo tanto

la solución de la homogénea:

es la función complementaria:

Ahora se encuentra la particular por variación de parámetros [ZILL]:

114
Donde

Aplicando el Wronskiano y los determinantes para determinar los parámetros

. Se conoce que la función no complementaria es:

Así

De esta manera:

Encontrando

Integrando por partes:

De la misma forma para :

115
Así

entonces el perfil general de velocidad esta dado por:

3.2.2.3.3 Perfil de velocidad

Ajustando a las condiciones de frontera:

[3-17]

3.2.2.4 Número de Nusselt

Obsérvese que si se utiliza la ecuación del perfil de temperaturas conociendo (el gradiente de
temperatura en la superficie del bulbo) entonces quedaría determinado el número de Nusselt.
Esta constante puede determinarse al aplicar la condición en en la ecuación del
perfil de velocidad pues se obtiene el gradiente de temperatura en la superficie del bulbo, es
decir .

Realizando lo anterior se obtiene , al ajustar la condición de en el perfil de

velocidad:

116
El recordar que y sustituirlo en la ecuación del perfil de temperatura, permite obtener

finalmente el número de Nusselt:

[3-18]

Nota: las dos barras indican valor absoluto.

Es decir que Nusselt depende solo del espesor de la capa límite térmica.

3.2.3 Termómetro En Posición Horizontal Con Flujo Cuasiestacionario

3.2.3.1 Consideraciones y Construcción

Para plantear el modelo correspondiente al cilindro horizontal, es necesario considerar que este
problema es más complejo que el de convección en cilindro vertical, dado que en este caso, la
fuerza de flotación al trabajar el problema en coordenadas cilíndricas no tiene una sola
componente sino, tiene dos componentes (la angular y la radial). Pese a que en coordenadas
rectangulares la fuerza de flotación tiene una componente, el problema se planteará en
coordenadas cilíndricas pues al trabajarlo en coordenadas rectangulares el problema es más
difícil de plantear, y sobre todo de resolver.

En coordenadas cilíndricas se tiene el vector velocidad:

Inicialmente se asumen las siguientes restricciones:

Al observar en el interferograma (Figura 3.7) de la visualización de las líneas de corriente del

fluido en convección libre alrededor de un cilindro horizontal se considera que cerca a la
superficie del bulbo (dentro de la capa límite) el movimiento es principalmente angular pues el

117
flujo tiende a seguir el contorno del cilindro; despreciando por ello las componentes radial y
axial de velocidad ante la angular, y por ello , para que el vector velocidad se
reduzca a:

Figura 3.7 Interferograma de cilindro enfriándose por convección natural [HOLMAN]

la fuerza de flotación por metro cúbico es

3.2.3.2 Caso de gradiente angular de temperatura constante

Al aplicar la ecuación de cantidad de movimiento expuestas en el numeral 2.1.1 se obtiene:

Componente radial:

Componente angular:

118
Considerando gradiente angular de temperatura constante e ignorando la generación de calor
por fricción viscosa y propiedades constantes bajo el supuesto de flujo incompresible (que es
válida si el Número de Mach es <0.3), la ecuación de energía expuesta en el mismo numeral
2.2.1 conduce a:

Donde Gradiente de temperatura angular

Combinando las componentes angular y radial de cantidad de movimiento y utilizando la

ecuación de energía, se obtiene la ecuación:

Condiciones en y en

Al Solucionar la ecuación anterior, en serie de potencias la distribución de velocidad angular

tiene la forma:

[3-19]

Combinando la ecuación de energía con el perfil anterior de velocidad, se obtiene

119
Bajo las condiciones límite en la superficie del bulbo

En
En

Solucionando y forzando condiciones límite se obtiene el perfil de temperatura

[3-20]

donde

Al derivar el perfil de temperatura y evaluarlo en

Se definen los números de Nusselt basados en el gradiente radial y angular de temperatura:

Y sustituyendo en la ecuación anterior del gradiente radial en , se obtiene

Pero

120
 √

Como una aproximación podría considerarse , lo cual conduce a

[3-21]

3.2.3.3 Caso sin gradiente angular de temperatura: ,

El caso anterior consideró gradiente angular constante de temperatura . Si ahora se considera

la temperatura angular constante, es decir el gradiente angular de temperatura igual a cero. Al
repetir todo el procedimiento anterior, El número de Nusselt es

[3-22]

Solución idéntica a la obtenida por Langmuir para una superficie curva [MARTINENKO]

NOTA: El valor de usado en las secciones 3.2.1 hasta 3.2.3 en realidad no corresponde al espesor de la
capa límite, , sino que corresponde a la distancia radial desde el eje axial del cilindro hasta el borde
de dicha capa, es decir que lleva incluido el radio, a saber: .

En el capítulo 2 si se usó correctamente este parámetro.

121
 4. DETERMINACIÓN EXPERIMENTAL

Es este capítulo se da una descripción detallada de la manera como partiendo de unos

ensayos previos se realizó la determinación experimental para obtener los coeficientes de
transferencia de calor de las situaciones propuestas, además de la exposición de los montajes
usados y el procedimiento a seguir. El desarrollo de este capítulo, junto con apartes de los
Capítulos 2 y 5, es la base para la construcción de un modelo de guía como práctica de
laboratorio y que se expone en el Capítulo 6.

4.1 ENSAYOS Y SELECCIÓN DE EXPERIMENTOS

Para realizar la parte experimental, de la convección libre alrededor del bulbo de un termómetro
de vidrio, se pensó inicialmente en calentar el bulbo del termómetro, y registrar en la cámara de
video como aumentaba la temperatura según el ascenso del menisco. Para calentar el bulbo
del termómetro se hicieron las siguientes pruebas preliminares.

4.1.1 Ensayo-1: Calentamiento del bulbo con agua a temperatura máxima permisible

Se calentó agua, agitándola magnéticamente, en un vaso de

precipitados de 500 mL, colocando el vaso sobre una plancha eléctrica
de calentamiento, hasta alcanzar la temperatura máxima del
termómetro, que es de 35C. El ascenso de temperatura en el vaso, se
monitoreaba con un termómetro de vidrio de los que se usan
comúnmente en las prácticas de operaciones unitarias. Una vez
alcanzada la temperatura de 35C se sumergía el bulbo del termómetro
de precisión en el vaso con agua caliente y se observaba como
ascendía el menisco en el termómetro; aquí el agua caliente no se
agitaba durante el calentamiento del bulbo. En este experimento de
calentamiento, se observó, que pese a que el agua no se agitaba, la
temperatura ascendía muy rápidamente en el termómetro de vidrio, es decir el menisco del
fluido termométrico subía a alta velocidad. Luego, se hicieron varios intentos de tratar mover la
cámara para seguir el movimiento del menisco, pero éste se movía de tal forma que

122
manualmente era prácticamente imposible seguir su movimiento. De esta manera no se
lograba registrar todo el ascenso de temperatura en el termómetro de vidrio.

4.1.2 Ensayo-2: Calentamiento del bulbo con agua a menor temperatura

Como resultó imposible seguir el ascenso de la temperatura en el experimento del bulbo

sumergido y calentado con el agua, se procedió inicialmente a calentar el agua, ya no hasta
35C, si no hasta una temperatura menor, como 30C; así como la diferencia de temperatura
inicial entre la superficie del bulbo y la temperatura inicial del agua, es menor que en el ensayo-
1, se esperaba que la velocidad de ascenso del menisco fuese menor en este caso que en el
ensayo-1. Ello debido a que a menor diferencia de temperatura menor velocidad de
calentamiento. Efectivamente el menisco en el termómetro ascendía a menor velocidad, pero
aun así, resultaba difícil seguir con la cámara el ascenso del menisco.

4.1.3 Ensayo-3 calentamiento del bulbo con aceite

En los dos experimentos anteriores no se fue posible registrar en el termómetro la variación

temporal de temperaturas, por ser muy rápido el movimiento, se pensó de los dos experimentos
anteriores, que el coeficiente de Transferencia de calor era relativamente alto en las
vecindades del bulbo y el agua. Para disminuir el coeficiente de transferencia de calor, y por
ende las corrientes convectivas naturales, se pensó en cambiar de fluido en el vaso: en lugar
de agua se utilizó un fluido con mayor viscosidad, con el ánimo de disminuir el movimiento del
mismo fluido, como el caso de aceite mineral que se calentó en el vaso de precipitados hasta
30C. Efectivamente el menisco ascendió a menor velocidad que en el caso del agua, para la
misma diferencia de temperatura inicial entre el bulbo y el fluido .y si se pudo, en este caso,
registrar la variación de temperatura con la cámara, como se encuentra registrado.

4.1.4 Ensayo-4: Enfriamiento del bulbo en agua y aceite a temperatura ambiente

De los tres experimentos anteriores se concluye, mediante observación cualitativa, que el

calentamiento del bulbo con fluidos líquidos es relativamente rápido, para seguir cómodamente
el registro de temperatura con la cámara. De esta manera, ahora en vez de registrar el
calentamiento del bulbo, se registrará su enfriamiento por lo que se calentó el bulbo del
termómetro con un mechero hasta alcanzar la temperatura máxima permisible en el
termómetro, pero sin registrar la temperatura; una vez se alcanzaba esta máxima temperatura,

123
hasta la cual ascendió el menisco, se inicia el registro fílmico y también simultáneamente se
enfría el bulbo al sumergirlo en un vaso con agua a temperatura ambiente. Se observó que en
el caso de enfriamiento del bulbo con fluidos líquidos el menisco en este caso descendía muy
rápidamente y no se podía hacer el registro con la cámara para registrar el cambio de la
temperatura respecto al tiempo.

4.1.5 Ensayo-5: Enfriamiento del bulbo con aire a temperatura ambiente

Se hicieron pruebas de enfriar el bulbo, pero en lugar de utilizar líquidos a temperatura

ambiente, se utilizó un gas, en este caso aire a temperatura ambiente. Se observó que en este
caso, si podía seguirse cómodamente la variación de temperatura con la cámara, es decir, la
cámara podía moverse casi que a la misma velocidad a la que descendía el menisco. Indicando
ello que el coeficiente de Transferencia de calor es, cualitativamente, menor con gases que con
líquidos. Por tal razón, los experimentos en los cuales se hicieron registros fílmicos de la
evolución de la temperatura, fueron de enfriamiento del bulbo utilizando aire en reposo a
temperatura ambiente.

4.2 DESCRIPCIÓN DEL EXPERIMENTO

El experimento consiste básicamente en tomar las lecturas de cambio de temperatura en

función del tiempo registradas por un termómetro de vidrio. Se considera como sistema que
transfiere calor con los alrededores, el bulbo del termómetro, el cual es asumido como un
cilindro.

El bulbo se calienta con una llama hasta llevar al termómetro a su límite

superior de temperatura. Se expone al medio ambiente externo (más frío, )
y una vez comienza el enfriamiento del bulbo con el aire estancado, se
registra el descenso de temperatura con una video cámara que sigue el
menisco de mercurio de la columna hasta que prácticamente llega a equilibrio
térmico con el ambiente. Para evitar corrientes de aire, el termómetro es
encerrado en un recipiente cilíndrico de muchísimo mayor diámetro. La temperatura ambiente y
la humedad relativa del momento del ensayo también son registradas.

124
En los siguientes dos experimentos se hizo entonces el registro con la cámara de la velocidad
de descenso de temperatura, para luego hacer la determinación experimental de los
coeficientes de transferencia de calor por convección libre; sobre los valores obtenidos de los
coeficientes, se prueban tanto el modelo del termómetro como algunas expresiones teóricas
reportadas por la literatura y generadas por los autores que trabajaron en esta tesis.

A continuación se describen estos experimentos seleccionados, y sobre los cuales se centra el

presente trabajo.

4.2.1 Experimento-1: Convección libre en aire alrededor del bulbo de un termómetro

orientado verticalmente (enfriamiento)

Figura 4.1 Orientación vertical del termómetro [El Autor]

Inicialmente se cuelga, en un soporte, el termómetro de vidrio con el bulbo en posición vertical

como se muestra en la Figura 4.7, situando adecuadamente una cámara digital de video al
frente del menisco del termómetro. Seguidamente, con un mechero se calienta el bulbo del
termómetro e inmediatamente el menisco del fluido termométrico comienza a ascender, se
prosigue el calentamiento hasta que el menisco casi que alcance la lectura máxima de la
escala. Instantes antes de que el menisco alcance la temperatura máxima, se activa la cámara
para que comience a filmar. Una vez alcanzada la temperatura máxima permisible se retira el
mechero y se sitúa rápidamente un vaso bajo el termómetro de tal forma que el bulbo del
termómetro quede sumergido dentro del vaso (ver Figura 4.7). Una vez se retira el mechero, el
menisco del termómetro comienza a descender, debido a que el bulbo del termómetro se enfría
por convección, en este caso natural, y a medida que desciende el menisco, es necesario
mover hacia abajo la cámara filmadora, con la misma velocidad a la que desciende el menisco,
y con el auxilio de una guía de cremallera que sirve también de soporte de la video-cámara. De
esta forma queda registrado en video la evolución temporal de la temperatura. El vaso es

125
necesario incluirlo para minimizar perturbaciones externas que puedan inducir corrientes en el
aire que rodea las vecindades próximas al termómetro. Se hizo el experimento por duplicado.

Inicialmente se realizo el experimento de enfriamiento pero sin colocar el bulbo dentro del vaso,
se registraron las temperaturas en la cámara. El experimento se hizo varias veces debido a que
cuando se repetía sin el vaso, se encontraba que la evolución de la temperatura en el
enfriamiento del bulbo difería bastante. Al analizar la situación se encontró que las mediciones
de temperatura en el experimento sin el vaso, eran muy sensibles a condiciones como, si la
puerta del laboratorio estaba abierta o cerrada, los extractores estaban o no prendidos y a si las
personas que se encontraban en el laboratorio cerca al termómetro, se movían o no. De esta
manera, una vez estuviese caliente el bulbo del termómetro, se opta por introducir el bulbo
dentro de un vaso, para así garantizar la condición de convección libre.

4.2.2 Experimento-2: Convección libre en aire alrededor del bulbo de un termómetro

orientado horizontalmente (enfriamiento)

Figura 4.2 Orientación horizontal del termómetro [El Autor]

En este experimento el termómetro se sitúa horizontalmente sobre un soporte elaborado para

tal fin como se muestra en la Figura 4.8. A continuación para lograr el registro del descenso de
temperatura del bulbo del termómetro, el experimento se desarrolla de forma similar al
experimento-1. En esta orientación, cualitativamente se observa que se alcanza más rápido el
equilibrio térmico con respecto al experimento-1, de lo que se infiere: los coeficientes de
transferencia de calor por convección libre alrededor de un cilindro horizontal son mayores que
los correspondientes a un cilindro orientado verticalmente.

4.3 MATERIALES E IMPLEMENTOS

A parte de los elementos de soporte diseñados y/o adecuados (universal, trípode y pinzas de
aro de nuez) para lograr los montajes que se ven en las figuras, las guías para la cámara, el

126
vaso de precipitados (de 500 mL o mayor capacidad), un mechero y un cronometro, para
materializar el experimento se requirieron los siguientes implementos esenciales.

4.3.1 Termómetro de Precisión

El termómetro usado es uno de repuesto de la bomba calorimétrica ubicada en el Laboratorio

de Ingeniería Química de la Universidad Nacional. Las características fundamentales de dicho
termómetro son las siguientes:

ASTM 56C – 86

INTERVALO: 19 °C – 35 °C
PRECISIÓN: 0.02 °C o °C (división)

Figura 4.3 Termómetro de la bomba calorimétrica [El Autor]

DIMENSIONES Bulbo
Columna Capilar Termómetro
mm Total Pared Punta

Longitud 47.40 1.3273 2.10 541.10 541.10 588.50

Diámetro 7.38 ‐ ‐ 7.40 0.20 ‐

Tabla 4.1 Dimensiones del termómetro usado [El Autor]

Estas dimensiones se midieron con un calibrador. Las especificaciones de este termómetro

pueden verificarse en la norma ASTM E1. (2007): Standard Specification for ASTM Liquid-in-
Glass Thermometers, pág. 1-6, 22, 49, 50, 54, 58-64; En la tabla 4.2 se muestran las más
importantes junto con la Figura 4.4. Según la norma este termómetro debe calibrarse cada 2 °C
desde los 19 °C.

127
 ASTM 56C‐86
Nombre Bomba Calorimétrica
Rango 10 a 35°C
Para pruebas de
A Inmersion, mm total
Graduaciones:
Subdivisiones 0.02°C
Líneas largas cada 0.1°C
Números cada 0.2°C

Error máximo de escala 0.10 °C

(sobre un intervalo de 2°C el cambio en la corrección no excede 0.02°C)
ASTM
Inscripción Especial
56C‐86
Cámara de Expansión
Permite calentamiento a 66°C
B Longitud total, mm 570 a 600
C Diámetro Externo de la Columna 7.0 a 8.0
D Longitud del Bulbo 35 a 55
E Diámetro Externo del Bulbo 7.0 a 8.0
(El diámetro del bulbo no puede tener más de 0.5 mm que el de la columna)
Localización de la Escala
Desde el fondo del bulbo a la línea de 19°C
F Distancia, mm 165 a 187
G Longitud de la parte graduada 323 a 385
Cámara de Contracción
I Distancia máxima hasta arriba,mm 76
Tabla 4.2 Especificaciones ASTM del termómetro usado [ASTM E1 (2007)]

Figura 4.4 Termómetro ASTM 56C -86 [ASTM E1 (2007)]

El grosor o espesor de la pared del bulbo del termómetro se midió con la ayuda de un
Estereoscopio Advanced Optical del Instituto de Ciencias Naturales de la Universidad Nacional
de Colombia como se muestra a continuación.

128
 Figura 4.5 Medición del espesor de la pared del Bulbo con el Estereoscopio

4.3.2 Video grabadora y Software

Las grabaciones se realizaron con una handycam Sony. Se realizaron dos registros fílmicos
para cada orientación. La filmación de cada uno se llevó al programa Vegas Pro, donde se
digitalizaron 12 cuadros por segundo, cada uno de los cuales quedó en una imagen. Para la
lectura se tomaron sólo 4 imágenes por segundo, para un cambio de ¼ de segundo.

Figura 4.6 Elementos de grabación y edición [El Autor]

129
4.4 MONTAJE

4.4.1 Cilindro Vertical

Para este montaje se requiere de un soporte universal y pinzas de nuez para mantener el
termómetro aplomado verticalmente, y de otro soporte móvil para la video cámara con el fin de
guiar el seguimiento de la filmación del menisco.

Figura 4.7 Montaje para el termómetro vertical [El Autor]

4.4.2 Cilindro Horizontal

Para este montaje se diseñó una guía para la cámara, en forma de riel, véase la Figura 4.8.

130
 Figura 4.8 Montaje para termómetro horizontal [El Autor]

4.5 PROCEDIMIENTO

A partir de la experiencia tenida se diseña la siguiente secuencia de pasas para llevar a cabo el
experimento en forma satisfactoria:

 Realizar el montaje que se muestra esquemáticamente en la figuras, ya sea para cilindro

vertical u horizontal.

 Cerrar la puerta del laboratorio para eliminar corrientes de aire que podrían influir en el
ensayo. Y garantizar que no se induzca movimiento de aire durante el ensayo.

 Acercar lentamente al bulbo del termómetro la llama del mechero y observar que el menisco
del mercurio en el termómetro alcance la temperatura máxima que se fije para el ensayo.
Esta temperatura no debe sobrepasar nunca la temperatura máxima de la escala del
termómetro, pues si esto ocurre el termómetro se descalibra y podría dañarse.

131
 Una vez alcanzada la temperatura máxima fijada, la cual se toma como temperatura inicial,
retirar el mechero y colocar el vaso de tal forma que el bulbo del termómetro quede dentro
del vaso como muestra la figura.

 Al final del procedimiento anterior ya la cámara de video debe estar filmando.

 Al enfriarse el bulbo del termómetro, para la orientación vertical el menisco del mercurio
desciende por lo cual se debe mover la cámara sobre el soporte de la cámara y en la caso
de la orientación horizontal la columna de mercurio se contrae por lo cual se debe mover la
cámara sobre el riel; esto es importante tenerlo en cuenta para seguir el movimiento del
menisco y así registrar y grabar la evolución de la temperatura en el tiempo.

 Durante el ensayo si hay personas alrededor, es importante que éstas no caminen cerca
del ensayo pues al desplazarse se originan corrientes forzadas de aire y esto afecta el
experimento pues ya no se tendría convección libre.

132
 5. ANALISIS DE RESULTADOS

Según lo expuesto en el capítulo anterior, se realizaron unos ensayos preliminares que

permitieron seleccionar los dos experimentos sobre los cuales se hace el estudio de la
convección. Los experimentos seleccionados fueron los de enfriamiento del aire por convección
libre alrededor del bulbo del termómetro. Uno para el termómetro en posición vertical y el otro
en posición horizontal.

Se presentan los datos de las mediciones realizadas para cada uno de los experimentos, tanto
verticales como horizontales, realizados para la determinación de los coeficientes de
transferencia de calor. Luego se exponen los tratamientos matemáticos que se efectuaron a
dichos datos y que se compilan de tres maneras: regresión por mínimos cuadrados, derivación
e integración numérica y corrección por resistencia de pared. Seguidamente se muestran los
resultados obtenidos por aplicar tales procesos para finalmente ser discutidos y analizados,
comparando contra los valores esperados por distintos autores de la literatura.

5.1 MEDICIONES

En el capítulo anterior se mencionó como se registraron los datos de temperatura versus

tiempo. Para la lectura del tiempo se utilizó el cronometro que trae incorporada la cámara y que
se activa cuando comienza la filmación.

Una vez obtenidos con la cámara los registros temporales de temperatura, se transfirió la
información digital de la cámara al computador. Posteriormente, se siguió el proceso de extraer
del video las series de tiempo de temperatura, para ello se corrió el video con uno de los
reproductores de video que trae el software de video del computador. Este proceso de lectura,
que se esperaba que fuese fácil, resultó en la realidad muy laborioso, pues si se quería leer la
temperatura para cierto intervalo regular de tiempo, no se lograba siempre lo mismo, ya que al
parar la ejecución del video para registrar el tiempo que aparecía en la pantalla del computador,
y haciendo varias repeticiones no se obtenía duplicidad. Así que fue necesario recurrir a un
software especial, Vegas Pro, que permite producir hasta 24 cuadros por segundo de
grabación, con los cual se obtuvo una serie de fotogramas, 12 por segundo, de los cuales se
tomó la decisión de sólo hacer la lectura de 4 por segundo.

133
Los 4 experimentos registrados en convección natural en aire para flujo externo alrededor de
un cilindro son:

∞
No. EXPERIMENTO ,s % H.R. No. Fotogramas
°C °C °C (leídos)

V‐1 Cilindro Horizontal 30.025 20.5 21.405 1557.50 43 6231

V‐2 Cilindro Horizontal (replica) 30.200 20.6 21.439 1317.75 43 5272

H‐1 Cilindro Vertical 34.880 20.6 21.280 742,50 43 2971

H‐2 Cilindro Vertical (replica) 34.825 20.7 20.960 1104.25 43 4418

Tabla 5.1 Experimentos trabajados

Los datos están compilados en el Apéndice D. Aquí se exponen las gráficas de dichas curvas
temporales, donde el intervalo entre uno y otro fotograma es de ¼ de segundo, es decir 0.25 s.

5.1.1 Cilindro Vertical

Figura 5.1 Fotogramas cilindro vertical

35
TEMPERATURA, TW (°C)

30
25
20
15
10
5
0
0 200 400 600 800 1000 1200 1400 1600
TIEMPO, t (s)

Figura 5.2 Experimento V-1

134
 35,000

TEMPERATURA, TW (°C)
30,000
25,000
20,000
15,000
10,000
5,000
0,000
0 120 240 360 480 600 720 840 960 1080 1200 1320

TIEMPO, t (s)

Figura 5.3 Experimento V-2.

5.1.2 Cilindro Horizontal

Figura 5.4 Fotogramas cilindro horizontal (sentido inverso)

40,000
35,000
TEMPERATURA, TW (°C)

30,000
25,000
20,000
15,000
10,000
5,000
0,000
0 60 120 180 240 300 360 420 480 540 600 660 720 780
TIEMPO, t (s)

Figura 5.5 Experimento H-1

135
 40
35
TEMPERATURA, TW (°C)
30
25
20
15
10
5
0
0 120 240 360 480 600 720 840 960 1080 1200
TIEMPO, t (s)

Figura 5.6 Experimento H-2

5.2 TRATAMIENTO DE DATOS

5.2.1 Determinación del intervalo de temperatura a trabajar

Usando el concepto de constante de tiempo, , como el tiempo requerido para lograr

aproximadamente al · % . % del cambio total de temperatura [SMITH], que
corresponde a . De esta manera podemos plantear que la temperatura inferior del
intervalo, , está dada por:

Figura 5.7 Significado de la Constante de Tiempo [SMITH]

136
Con ello se busca rechazar la zona de la curva que obedece a una respuesta estacionaria, ya
que de acuerdo a la teoría de control cuando se tiene un · % . % del cambio
transitorio prácticamente se estabiliza en un estado estacionario. Cuando se acerca al equilibrio
termodinámico el valor de los coeficientes instantáneos de transferencia de calor tiende a
disminuir y a ser cero, debido a que la diferencia de temperatura también tiende a reducirse y a
ser cero. De esta manera se asegura el estado transitorio, al trabajar en la región exponencial
definida por , .

Es decir que de la serie de datos obtenidos, correspondientes a la evolución de la temperatura

para cada experimento, sólo se toman desde hasta , lo que implica rechazar desde
hasta .

5.2.2 Parámetros

5.2.2.1 Propiedades Físicas

Las propiedades se calcularon teniendo en cuenta la temperatura promedio entre la inicial y

final, que se evalúa por el teorema integral del valor medio, y la integral a su vez se resuelve
mediante la Regla de Simpson ⁄ :

∆
...

Y la temperatura de película promedio con la que se determinan las propiedades del aire es:

Los valores de las variaciones de las propiedades con la temperatura se extrajeron de la

literatura [KRETIH], con los cuales se generaron (por regresión) las siguientes correlaciones en
SI (la temperatura en °C, excepto para el cp,  y  del aire que además estarán en K):

Mercurio:

. .

. . . .

para el intervalo: ° °

137
Aire:
. .

. .

para el intervalo: ° °

. . . .

para el intervalo:

con [CENGEL]:

. .

5.2.2.2 Condiciones de Trabajo

Las condiciones de presión atmosférica local y aceleración gravitacional son:

5.2.2.3 Especificaciones del Bulbo

Para la envolvente cilíndrica o pared del bulbo se tienen las siguientes especificaciones:

, mm 7.38
, mm (espesor) 1.3273
, mm 4.73
, mm 45.3
Material Vidrio Pyrex
, W/m°C 1.09 [MILLS]
Tabla 5.2 Especificaciones de la pared del bulbo

138
5.2.3 Regresión por mínimos cuadrados

Para determinar los coeficientes de transferencia de calor y comprobar el modelo va a utilizarse

el método de mínimos cuadrados aplicando a la solución de la ecuación [3-7]

θ
θ

Si se arregla de la siguiente forma:

donde:

θ
θ

como θ entonces:

Aplicando el método de mínimos cuadrados puede obtenerse el coeficiente convectivo en

términos de datos experimentales:

∑
∑

esta última asume que el intercepto con el eje de la variable dependiente es cero, de lo
contrario se debe asumir que [CHAPRA]:

∑ ∑ ∑
∑ ∑

El índice de correlación, o coeficiente de determinación se encuentra mediante la fórmula

[CHAPRA]:

∑ ∑ ∑
∑ ∑ ∑ ∑

139
Aplicando lo anterior a todos los datos se obtiene un coeficiente de transferencia de calor
promedio para todo el proceso de enfriamiento de acuerdo a la definición del número de Biot
dada 3.1:

Y el Nusselt promedio será:

5.2.4 Derivación e integración numérica

Para cada pareja de datos , o punto de la curva de enfriamiento se determina

numéricamente la derivada, ⁄ , o cambio instantáneo de temperatura usando las formulas
de primera derivada por diferencias finitas con cinco puntos: hacia adelante para los dos
primeros puntos, las de hacia atrás para los dos últimos puntos y centradas para el resto de
puntos [CHAPRA]. Adaptando y aplicando dichas formulas para el caso:

Luego se calcula cada coeficiente de transferencia de calor instantáneo a través de la

definición:

Ahora por el teorema integral del valor medio se determina el valor promedio del coeficiente:

140
Dicha integral definida se resuelve por la Regla de Simpson de ⁄ , la cual requiere un número
par de particiones, . Si los datos arrojan número impar de particiones, se elimina el último par
de mediciones , .

El error de esta determinación numérica es:

Para la determinación de este error, se procede con la serie temporal de coeficientes de calor
instantáneos obtenidos anteriormente de la siguiente manera. Para cada tiempo se calcula
numéricamente la cuarta derivada usando las siguientes formulas adaptadas a la presente
situación [CHAPRA, MATHEWS]:

Para los tres primeros puntos se aplica la formula de diferencias divididas finitas hacia
adelante, para los tres últimos puntos la formula de diferencias divididas finitas hacia atrás y
para el resto de puntos la formula de diferencias divididas finitas centradas. Luego se ubica
tanto el máximo y como el mínimo de estos valores de 4ª derivada, y entre estos dos se toma la
4ª derivada de mayor valor absoluto, con la cual se calcula en la formula el error máximo
debido a la aplicación de la integración numérica.

141
5.2.5 Corrección por Resistencia de Pared

Finalmente la corrección por haber incluido la resistencia térmica conductiva de la pared del
bulbo (en realidad un coeficiente global aire-vidrio), y que está dada en la ecuación [3-13]
permite hallar el coeficiente de transferencia de calor individual del aire:

donde es la resistencia de pared dada por la ecuación de la resistencia conductiva en un

cilindro hueco:

Este valor del largo del bulbo tiene excluido el grosor de la punta, que también es una
envolvente de vidrio de espesor 2.1 mm.

. . .

De esta manera la resistencia conductiva de pared es:

. °C
.
· . ⁄ °C · . .

Como la superficie lateral exterior del cilindro es: . . Entonces la corrección

por resistencia de pared para hallar el verdadero coeficiente de transferencia de calor es:

. · . °C⁄ .

5.3 RESULTADOS

Los valores de Nusselt promedio obtenidos experimentalmente para los cilindros verticales
están referenciados con respecto a la altura, L, mientras que para cilindros horizontales están
referenciados con respecto al diámetro, D.

142
5.3.1 Intervalos de temperatura

Exp. No. V‐1 V‐2 H‐1 H‐2

, 24.004 24.068 25.79011 25.772

, 333.00 340.50 196.00 199.50
Tabla 5.3 Determinación de la temperatura correspondiente a la constante de tiempo

5.3.2 Regresión por Mínimos Cuadrados

Y Z YZ Y2 Z2

Exp.
n tm, s
No.

V‐1 1333 ‐666.5000 ‐705.2689 463.8816 444.5001 484.5755 333.00

V‐2 1363 ‐681.5000 ‐722.3564 474.6259 454.5001 496.1782 340.50
H‐1 785 ‐392.5000 ‐417.2790 274.0391 261.8335 287.1623 196.00
H‐2 799 ‐399.5000 ‐414.2404 273.3962 266.5002 280.6190 199.50

Exp. Con intercepto en el origen Sin intercepto en el origen

R2
No. Bi h*, W/m2K Nu* Bi h*, W/m2K Nu*
V‐1 0.998341 1.043603 4.2768 8.0000 0.999973 4.0980 7.6656
V‐2 0.998226 1.044281 4.1851 7.8275 0.997341 3.9970 7.4756
H‐1 0.997935 1.046616 7.2783 2.1111 0.997195 6.9347 2.0114
H‐2 0.999218 1.025876 7.0087 2.0327 0.992896 6.7834 1.9674
Tabla 5.4 Evaluación del coeficiente por regresión con mínimos cuadrados

5.3.3 Derivación e integración numérica

Exp.
hici tm, s h*, W/m2K Nu* Error Máx
No.

V‐1 16387.05 333.00 4.1009 7.6710 0.34

V‐2 16385.35 340.50 4.0101 7.5002 0.29
H‐1 16363.02 196.00 6.9571 2.0179 0.38
H‐2 16349.15 199.50 6.8292 1.9806 0.53
Tabla 5.5 Evaluación del coeficiente por derivación e integración numérica

143
5.3.4 Valores corregidos por resistencia de pared

Regresión
Derivación e integración
Exp. Con intercepto Sin intercepto
numérica
No. en el origen en el origen

h, W/m2K Nu h, W/m2K Nu h, W/m2K Nu

V‐1 4.3059 8.0544 4.1247 7.7155 4.1276 7.7210
V‐2 4.2130 7.8796 4.0224 7.5231 4.0357 7.5480
H‐1 7.3630 2.1357 7.0114 2.0337 7.0343 2.0404
H‐2 7.0872 2.0554 6.8569 1.9887 6.9037 2.0022
Tabla 5.6 Valores corregidos por resistencia de pared

Los cambios porcentuales en el coeficiente al introducir esta corrección son:

Regresión
Exp. Derivación e integración
Con intercepto Sin intercepto
No. numérica
en el origen en el origen

V‐1 0.68 0.65 0.65

V‐2 0.67 0.64 0.64
H‐1 1.16 1.11 1.11
H‐2 1.12 1.08 1.09
Tabla 5.7 Aumentos porcentuales debidos a la corrección por resistencia de pared

5.3.5 Valores esperados

5.3.5.1 Cilindro Vertical

EXPERIMENTO V‐1 V‐2

2
MODELO Nu h, W/m K Nu h, W/m2K
Predicción Teórica 10.7897 5.7681 10.8201 5.7852
Cebeci 12.5340 6.7006 12.5613 6.7161
Cebeci Pr=0.72 12.7041 6.7916 12.7613 6.8070
Yang 11.3330 6.0586 11.3670 6.0776
Janna* 2.2105 7.5898 2.2163 7.6110
Elebass* 2.7036 9.2830 2.7147 9.3224
Kuehn* 1.9118 6.5642 1.9166 6.5818
Tabla 5.8 Valores esperados de coeficientes de transferencia de calor para cilindros verticales
*Usan NuD, el resto son con NuL

144
Exp. No. T, °C TWpro, °C Tf pro, °C Pr GrL RaL GrD RaD
V‐1 20.50 26.35 23.43 0.73023 45883.09 33504.99 173.18 126.46
V‐2 20.50 26.45 23.48 0.73021 46640.81 34057.67 176.04 128.54
Tabla 5.9 Condiciones y valores de números adimensionales para el aire usados en cilindros verticales

5.3.5.2 Cilindro Horizontal

EXPERIMENTO H‐1 H‐2

MODELO NuD h, W/m2K NuD h, W/m2K
Merk and Prins 1.6107 5.5530 1.6141 5.5654
Hermann 1.4868 5.1260 1.4900 5.1374
Eckert 1.9579 6.7502 1.9621 6.7653
Churchill 1.8121 6.2474 1.8152 6.2587
Churchill Laminar 1.8904 6.5175 1.8934 6.5286
Morgan 2.2709 7.8291 2.3136 7.9873
Kuehn 2.3268 8.0220 2.3301 8.0343
Tabla 5.10 Valores esperados de coeficientes de transferencia de calor para cilindros horizontales

Exp. No. T, °C TWpro, °C Tf pro, °C Pr GrD RaD

H‐1 20.50 29.31 24.90 0.7298 255.19 186.25
H‐2 20.50 29.39 24.95 0.7298 257.36 187.83
Tabla 5.11 Condiciones y valores de números adimensionales para aire usados en cilindros horizontales

5.3.6 Comparación

Se toma como base de referencia para el cilindro vertical la regresión de la solución de Cebeci,
mientras que para el cilindro horizontal la de Churchill (Laminar con 4.5% de error), que son las
formulas más complejas y las que ofrecen mayor precisión, además de su reconocimiento.

5.3.6.1 Cilindro Vertical

Se enlistan los valores del coeficiente, en W/m2K, de transferencia de calor considerando las
correcciones por resistencia de pared.

145
 REFERENCIA EXPERIMENTAL TEÓRICO
Exp. No. Regresión con Regresión sin Derivación e
Cebeci Predicción
intercepto cero intercepto cero integración

V‐1 6.7006 4.3059 4.1247 4.1276 5.7681

V‐2 6.7161 4.2130 4.0224 4.0357 5.7852
Tabla 5.12 Coeficientes de transferencia de calor para cilindros verticales

Los errores relativos porcentuales tomando como base la referencia son:

Regresión con Regresión sin Derivación e
Exp. No. Predicción
intercepto cero intercepto cero integración

V‐1 ‐35.74 ‐38.44 ‐38.40 ‐13.92

V‐2 ‐37.27 ‐40.11 ‐39.91 ‐13.86
Promedio ‐36.50 ‐39.28 ‐39.15 ‐13.89
Tabla 5.13 Errores relativos porcentuales de los coeficientes para cilindros verticales

Tomando el promedio del error relativo porcentual de cada uno de los dos ensayos en
comparación respecto a cada uno de los modelos:
Regresión con Regresión sin Derivación e
MODELO Predicción
intercepto cero intercepto cero integración

Predicción Teórica ‐26.26 ‐29.48 ‐29.34 ‐

Cebeci ‐36.50 ‐39.28 ‐39.15 ‐13.89
Cebeci Pr=0.72 ‐37.35 ‐40.09 ‐39.97 ‐15.04
Yang ‐29.80 ‐32.87 ‐32.73 ‐4.80
Janna* ‐43.96 ‐46.40 ‐46.30 ‐24.00
Elebass* ‐54.21 ‐56.21 ‐56.12 ‐37.90
Kuehn* ‐35.20 ‐38.02 ‐37.90 ‐12.12
Tabla 5.14 Errores relativos porcentuales promedio por ensayo para cilindros verticales

5.3.6.2 Cilindro Horizontal

REFERENCIA EXPERIMENTAL
Exp. No.
Churchill
Con intercepto cero Sin intercepto cero Derivación e Integración
(Laminar)
H‐1 6.5175 7.3630 7.0114 7.0343
H‐2 6.5286 7.0872 6.8569 6.9037
Tabla 5.15 Coeficientes de transferencia de calor para cilindros horizontales

146
 Regresión con Regresión sin
Exp. No. Derivación e integración
intercepto cero intercepto cero

H‐1 12.97 7.58 7.93

H‐2 8.56 5.03 5.75
Promedio 10.76 6.30 6.84
Tabla 5.16 Errores relativos porcentuales de los coeficientes para cilindros horizontales

Regresión Regresión
MODELO Derivación e integración
intercepto cero no intercepto cero
Merk and Prins 29.97 24.73 25.36
Hermann 40.80 35.13 35.80
Eckert 6.92 2.61 3.13
Churchill 15.55 10.89 11.45
Churchill Laminar 10.76 6.30 6.84
Morgan ‐8.61 ‐12.30 ‐11.86
Kuehn ‐10.00 ‐13.63 ‐13.19
Tabla 5.17 Errores relativos porcentuales promedio por ensayo para cilindros horizontales

5.4 ANALISIS Y DISCUSIÓN

Desde el punto de vista meramente cualitativo, en los experimentos preliminares se observó

que cuando se utilizan líquidos bien sea para calentar o enfriar el bulbo del termómetro, la
temperatura variaba más rápidamente que cuando se utiliza un gas (aire). Ello es indicio que el
transporte de calor es mayor en los fluidos líquidos que en el fluido gaseoso (aire) y como el
flujo de calor convectivo es proporcional al coeficiente convectivo de transferencia de calor
entonces el coeficiente de transferencia de calor bajo las mismas condiciones es mayor en los
fluidos líquidos que transportaron el calor.

Al comparar el experimento denominado ensayo-1 con el experimento denominado ensayo-2,

el cual se hizo con una diferencia de temperatura mayor entre fluido y bulbo, se
observó que el menisco descendía más rápidamente cuando la diferencia de temperatura
era menor, por lo cual el coeficiente convectivo de transferencia de calor es menor en

147
este caso. Por lo cual es de esperarse que en el ensayo-2 los coeficientes de Transferencia de
Calor sean menores con respecto al ensayo-1.

Ahora desde un punto de vista más cuantitativo se observa que la totalidad de los valores de
los coeficientes de calor encontrados experimentalmente están dentro de los valores
esperados, por lo menos en cuanto al orden de magnitud, caso de convección natural en aire,
tanto para el cilindro vertical como para el vertical ya que todos están oscilando alrededor del
orden de magnitud esperado de 5 W/m2K de convección libre en aire. Sin embargo, los valores
obtenidos del coeficiente de transferencia de calor para los cilindros verticales presentan una
mayor desviación de los valores esperados por la literatura frente a los cilindros horizontales.

Los índices de correlación muestran un excelente ajuste de los datos experimentales al modelo
planteado para el termómetro. Como el menor valor de correlación fue de 0.997935 para el
experimento H-1, podemos decir que el modelo de ajuste del termómetro por lo menos explica
aproximadamente el 99.79% de la incertidumbre original. Este ajuste en promedio del 99.8% se
pudo lograr gracias a la exclusión y rechazo de la parte estacionaria del intervalo de
temperaturas trabajado, y por ende solo se maneja el 63.21% del cambio total teórico posible
de temperatura, desde hasta . Si esto no se hubiera efectuado el menor valor de
correlación caería hasta un 91%.

Se observa también una notoria concordancia entre los coeficientes de transferencia de calor
obtenidos por diferentes procesos numéricos, es decir hay una buena convergencia entre los
coeficientes arrojados por la regresión por mínimos cuadrados y los que se procesaron con
derivación e integración numéricas. Comparando numéricamente, se tiene que los menores
valores de coeficientes se obtienen con la regresión sin intercepto en el origen, siendo
levemente mayores (entre un 0.1 y 0.7% aproximadamente) aquellos que resultan con la
derivación e integración numéricas, mientras que los de regresión con intercepto en el origen
aumentan (entre un 3 a 5% aproximadamente), es decir son los mayores. Esto indica que el
método de regresión sin intercepto en el origen resulta equivalente al método numérico
derivativo-integrativo.

Es notorio que los errores máximos generados al realizar la integración numérica se

encuentren entre el 5.5 y el 8.3%. La confrontación de los coeficientes obtenidos experimental y
teóricamente con respecto a los valores esperados se da a continuación; para ello se verificó
que los números de , y de las Tablas 5.9 y 5.11 de los experimentos realizados

148
pertenecieran a los intervalos de las restricciones de cada una de las correlaciones empleadas
para así poder habilitar su uso; también se aclara que no todas las correlaciones extraídas de
la literatura declaran la incertidumbre y margen de error, solo la de Churchill Laminar con 4.5%.

Para cilindros verticales se tiene que los valores determinados experimentalmente están por
debajo de los esperados tanto por la predicción teórica como por los de otros autores de la
literatura. El error oscila deficitariamente entre un 26.26 (Teórica) hasta un 56.21% (Elenbass);
frente a los modelos de predicción teórica y de Yang se tienen las menores desviaciones
(alrededor del 30%), mientras que ante los modelos de Elebass y Janna se presentan los
mayores alejamientos (del 50% aproximadamente). Los menores errores frente a los modelos
siempre se logran cuando se realiza el tratamiento de datos experimentales con regresión con
intercepto en el origen. Con respecto al modelo que ofrece mayor confiabilidad y exactitud,
Cebeci, los valores obtenidos están casi que un 38% por debajo de los esperados.

Una posible explicación de esta desviación, es una reducción muy notable de la turbulencia
natural producida por la interacción del vaso alrededor del bulbo, tal vez lo mas recomendado
sea colocar un recipiente de mucho mayor diámetro o hacer el experimento dentro de una
cámara o campana. También por comparación con los correspondientes coeficientes en una
pared plana vertical, se observa que los coeficientes determinados experimentalmente son
levemente superiores en alrededor un 3.5%, lo que no está siendo notable el efecto de la
curvatura; sin embargo, este debe ser tenido en cuenta ya que se encuentra apenas en un
6.5% para cumplir con el criterio de Gebhart y así tratar el cilindro vertical como una pared
como una plana.

Analizando ahora el valor del coeficiente para el cilindro vertical obtenido por predicción con el
método integral, la máxima desviación frente a los resultados experimentales es de alrededor
del 30.47%, la condición más favorable se da cuando se compara con los arrojados por el
procesamiento de regresión-intercepto cero, del 25.35%. Comparando la predicción contra los
otros modelos de la literatura se tiene un error relativo desde el -4.8% (Yang) hasta un -37.9%
(Elebass). Frente al modelo de Cebeci, los valores predichos teóricamente están alrededor del
13.9% por debajo. El error medio frente al conjunto de modelos es del 18.8%.

El valor predicho se encuentra por debajo de los valores esperados por la literatura, lo cual
aparece como obvio, ya que durante el desarrollo del modelo teórico (empezando que se utilizó
un método integral aproximado) se hicieron varias aproximaciones, como por ejemplo el uso de
apenas dos términos para los perfiles de capa limite y velocidad de escala propuestos para la

149
resolución del sistema de ecuaciones diferenciales por el método de perturbación, además de
una siguiente aproximación con cuatro términos al usar una expansión del logaritmo natural en
serie de Maclaurin para el coeficiente promedio de transferencia de calor en las ecuaciones [2-
14] y [2-15]; es de esperarse que el valor del número de Nusselt y por ende de los coeficientes
predichos por la ecuación obtenida sean menores a los esperados, debido a que se descartan
términos que aportan en el incremento del efecto de la curvatura.

En términos generales hay una muy buena convergencia entre los valores de coeficientes
predichos para el cilindro vertical con los coeficientes que se esperan por la literatura, como ya
se mencionó, su valor es siempre inferior a los esperados, y se ajustan mucho al modelo de
Yang (solo un 4.8% de diferencia), mientras que a modelos como el de Cebeci y Kuehn que
son más exactos la variación va de un 12 a 14%. Es decir que hay convergencia del modelo de
predicción obtenido conforme a las correlaciones de la literatura para hallar el coeficiente de
transferencia de calor.

Lógicamente los valores de los coeficientes deben ser mayores que a los de una pared plana
para las mismas condiciones, el coeficiente de transferencia de calor promedio es de alrededor
de 4 W/m2K por el modelo de Pared Plana de Churchill y Chu. Lo que se debe al efecto de la
curvatura convexa que estimula el fenómeno de transferencia. Es de esperar que en superficies
de curvatura cóncava, los coeficientes de transferencia de calor sean menores a los de una
superficie plana, ya que presentan mayores zonas de estancamiento, es decir que la forma de
la superficie puede favorecer o desfavorecer el flujo. La expresión teórica que se logró para
predecir coeficientes de transferencia de calor por convección libre alrededor de cilindros
verticales aplicada al bulbo del termómetro en las condiciones de trabajo, indica un aumento
del 45% del valor del coeficiente respecto a si su hubiera asumido como una pared plana
correspondiente. El segundo término de la expresión que indica el aporte de la curvatura al
fenómeno representa para estos casos alrededor de un 31% del total del coeficiente.

Para cilindros horizontales los resultados arrojados indican que hay una mayor convergencia
con los modelos de la literatura, la desviaciones están desde -13.63% (Kuhen) hasta 40.80%
(Hermann). En general se tiene un error promedio del 9.4% ante la totalidad de los modelos.
Los valores obtenidos están en la mayoría de caso por encima de los esperados, a excepción
de Morgan y Kuehn que reportan valores más altos (Error del 32%). El modelo al que más se
ajusta la determinación es al de Eckert ya que el valor obtenido por procesamiento con
regresión por intercepto arroja un 2.61% de error y por derivación e integración numérica un

150
3.13% de error. El modelo más alejado a los resultados es el de Hermann ya que tiene
alrededor de un 37% de error; no olvidar que el modelo de Hermann es meramente teórico.
Enfrentado al modelo de Kuehn, que es el más complejo y que viene de una deducción de
teoría de capa límite, los resultados de la experimentación están alrededor de un 12.3% por
debajo de lo esperado, mientras que al de Churchill Laminar están 8% por encima.

El análisis efectuado a la comparación entre los valores obtenidos experimentalmente a través

de los tres procedimientos para cilindros verticales es similar al de horizontales: Los mayores
valores son los de regresión con intercepto en el origen (4% en promedio mayores a los de sin
intercepto en el origen) y a la vez también son los que presentan mayor desviación (alrededor
de 12.2% en promedio) frente a los valores esperados por la literatura; los valores obtenidos a
través de derivación e integración numérica son 0.5% superiores a los de regresión sin
intercepto en el origen y en conjunto estos dos métodos arrojan resultados alejados en
promedio tan solo un 7.9% de los esperados.

Al realizar la corrección de los valores de los coeficientes obtenidos por la inclusión en estos de
la resistencia de pared, lo que los hacen globales vidrio-aire, se observa lógicamente que al
lograr lo coeficientes individuales aire, estos son mayores que cuando estaban sin corregir. En
los experimentos para el cilindro vertical las correcciones permitieron un aumento relativo
porcentual en casi un 0.65%; mientras que para los efectuados con el cilindro horizontal fue de
alrededor de 1.10% aproximadamente. Estos incrementos se hacen mayores a medida que
también sean mayores los coeficientes globales a corregir (para convertirlos en individuales).

En conjunto, y tomando un promedio, los modelos de la literatura predicen que los coeficientes
para el mismo cilindro y a las mismas condiciones son ligeramente superiores para un cilindro
horizontal, que para uno vertical en apenas un 2%. Cualitativamente se observó que los valores
obtenidos horizontalmente fueron mayores –ya que el tiempo de llegada al equilibrio fue menor-
y cuantitativamente también se verifica este hecho, pero de manera exagerada ya que los
horizontales obtenidos son un 70% mayores a los verticales. Con esto se refuerza lo ya
expresado anteriormente: los resultados experimentales para cilindro vertical presentan más
error que los horizontales.

Entonces la desviación de los valores arrojados experimentalmente para los coeficientes de

transferencia de calor en cilindros verticales encierran alguna explicación física, explicación
diferente a la de los errores de medición, ya que el déficit de los valores está en alrededor de
un 35%.

151
 6. PRÁCTICA: COEFICIENTE DE TRANSFERENCIA DE CALOR
POR CONVECCIÓN LIBRE ALREDEDOR DE UN CILINDRO

6.1 OBJETIVOS

 Determinar el coeficiente de transferencia de calor en situaciones de convección libre

alrededor de un cilindro.
 Comparar el coeficiente experimental con los obtenidos a través de las correlaciones
teóricas y empíricas.
 Determinar el número de Nusselt para las configuraciones trabajadas.

6.2 FUNDAMENTO TEÓRICO

Usualmente el coeficiente de transferencia de calor viene reportado en términos del número de

Nusselt, el cual se define como:

Donde : Número de Nusselt basado en la longitud característica ; : Coeficiente de

transferencia de calor (vatios/m2K); : conductividad térmica del fluido.

El experimento que permite determinar el coeficiente de Transferencia de Calor, está basado

en la ley de enfriamiento de Newton:

Donde: : Flujo de calor [W]; : Área de Transferencia de Calor Convectiva [m2]; :

diferencia de temperatura entre la superficie del bulbo del termómetro y la temperatura en el
espesor de la capa límite térmica [°C ó K].

Para transferencia de calor alrededor de un cilindro (vertical u horizontal) en contacto con un

fluido de diferente temperatura existen varias correlaciones en la literatura. Algunas de ellas
son teóricas y otras completamente empíricas. A continuación se listan algunas de las
anteriores además de las desarrolladas en esta tesis.

152
6.2.1 Correlaciones Teóricas

Para un cilindro vertical y basándose en la hipótesis aproximada de Langmuir de capa limite

estacionaria de una superficie curva, Elenbaas (1948) planteó:

/
⁄
.

Si ∞ se tiene que

⁄
.

solución que corresponde a una pared plana.

De hecho el criterio para tratar un cilindro vertical igual que una pared plana vertical fue dado
por Gebhart:

Las soluciones distan de menos 5% de la solución de pared plana para Pr desde 0.72 a 1.0

LeFevre y Ede (1956) emplearon un método integral para solucionar las ecuaciones de
gobiernan el fenómeno obteniendo la siguiente expresión del número de Nusselt promedio para
régimen laminar [BEJAN]:

⁄
·

Donde tanto y estan basados en la altura, , del cilindro vertical.

Hermann (1936) dedujo una expresión para cilindros horizontales en régimen laminar
resolviendo las ecuaciones diferenciales por el método de no similaridad, estableciendo que
para aire:

⁄
.

Merk and Prins (1953-1954) obtuvieron para cilindros horizontales en régimen laminar una
solución empleando el método integral y asumiendo espesores iguales entre las capas limites
térmica e hidrodinámica [BEJAN]:

153
 ⁄ ⁄

Para valores de de 0.7 (aire) el valor de la constante, fue calculado como 0.436.

6.2.2 Correlaciones Empíricas

 Cilindros Verticales

Para cilindros verticales con un flujo de calor constante se tiene [JANNA]:

C n

0.05 0.93 0.05

4
0.05-10 1.37 0.16
104 0.6 0.25
Tabla 6.1 Constantes para la ecuación de convección libre en cilindros verticales

Todos los valores de las propiedades se determinan a la temperatura de película, , que

corresponde al promedio aritmético entre la temperatura de pared, , y la temperatura del
seno del fluido , excepto el coeficiente de expansión térmica, , que en el caso de gases se
determina a .

Se recuerda que para un gas ideal como el aire

Una correlación general recomendada por Yang para regiones en régimen laminar y turbulento
tiene la forma [POPIEL]:

. ⁄
.
. ⁄ ⁄
. ⁄

154
Los valores de los cálculos numéricos obtenidos por Cebeci para un parámetro de curvatura

√
.

pueden ser representados por la ecuación simple, para . , en la región laminar

[POPIEL]:

.
. .
, .

donde , corresponde al de una pared plana vertical, dado por la correlación teórica
[KARLEKAR]:

⁄ ⁄ ⁄
, . .

o también por la correlación empírica desarrollada por Churchill y Chu para una placa vertical
de altura, , y con un borde frontal abrupto para flujo laminar [MILLS] :

⁄
, . .

donde la función del número de Prandtl, , definida por Churchill y Usagi se define como:

⁄ ⁄
.

Los valores de Cebeci también pueden ser aproximados por una ecuación validada para
números de Prandtl entre 0.01-100, es decir, . [POPIEL]:

. .
,

donde:

.
. .

. ⁄
. . . .

Popiel obtuvo una correlación a partir del enfriamiento de un cilindro en estado transitorio que
tiene una incertidumbre promedio de 3.5% con un 95% de confianza, válida para:

155
 , , y, ⁄ :

donde:

. . . .

. . .

La dependencia generalizada de transferencia de calor por convección libre en un gas para

flujo laminar y turbulento sobre un cilindro vertical está dada por la expresión de Kuehn
[MARTINENKO]:

⁄
⁄ ⁄ ⁄
. ⁄ . ⁄ .

Que es válida para:

. ⁄ , , /

 Cilindros Horizontales

Para y . se correlacionaron los datos experimentales de Eckert y

Soehnghen [KREITH]:

⁄
.

En forma más general Morgan para un cilindro isotérmico sugiere una expresión de la forma
[INCROPERA]:

156
 RaD C n
-10 -2
10 -10 0.675 0.058
-2 2
10 -10 1.020 0.148
102-104 0.850 0.188
4 7
10 -10 0.480 0.250
107-1014 0.125 0.333
Tabla 6.2 Constantes para la ecuación de convección libre en cilindros horizontales

Churchill y Chu dan la siguiente correlación válida para [MILLS]:

⁄
.
. ⁄ ⁄
. ⁄

Cuando ocurre transición de una capa límite laminar a una turbulenta y el aumento
del número de Nusselt con el número de Rayleigh es mayor [MILLS]:

⁄
.
. ⁄ ⁄
. ⁄

Por el método de la capa limite Kuehn deduce [MARTINENKO]:

⁄
⁄ ⁄ ⁄
. . ⁄ .

En la región de régimen laminar para , la aproximación que realiza Churchill se puede

aplica con un error del 4.5% [MARTINENKO]:

⁄
⁄ ⁄
. ⁄ ⁄
. ⁄

6.2.3 Modelo Experimental

En la determinación experimental del coeficiente de transferencia de calor es necesario

disponer del modelo de la transferencia de calor en el bulbo de un termómetro. El modelo se
obtiene al aplicar un balance de energía al sistema (columna y bulbo del termómetro). La
evolución de la temperatura, , con el tiempo, , instantes después de someter a enfriamiento o
calentamiento el bulbo de un termómetro de vidrio es:

157
Donde es la temperatura inicial en el termómetro, la temperatura en el seno del fluido,
el área del bulbo del termómetro, el coeficiente de transferencia de calor, y el calor
específico y la masa del fluido termométrico respectivamente.

En la forma adimensional de esta ecuación interviene el Número de Biot, , el cual

contiene el coeficiente de Transferencia de Calor:

Donde es la temperatura adimensional, , es el tiempo adimensional con

como el tiempo en que tarda el termómetro en llegar al equilibrio térmico práctico.

El siguiente modelo es válido bajo las siguientes restricciones:

 Se desprecia la radiación en el aire debido a la baja diferencia de temperaturas. Es

decir que mecanismo predominante en el aire es la radiación.
 En el fluido termométrico sólo se considera flujo conductivo de calor y el calor se
propaga instantáneamente (difusividad térmica infinita) y no se considera flujo
convectivo dentro del fluido termométrico. Esto equivale a que la columna tenga la
misma temperatura del bulbo que en la realidad correspondería a un termómetro de
inmersión total.
 Con respecto a lo anterior puede decirse que para el bulbo se aplica el método de
resistencia interna despreciable, que es válido para números de Biot menores a 0,1.
 La transferencia de calor entre el medio externo y el mercurio sólo se da a través del
bulbo del termómetro, mientras que la columna de vidrio que contiene el resto de
mercurio se supone que no transfiere calor desde o hacia el ambiente externo que lo
rodea, es decir, el capilar de vidrio se considera adiabático, perfectamente aislado.
 Se supone que el espesor de la pared de vidrio del bulbo del termómetro es tal que la
resistencia conductiva de calor en la pared de vidrio es pequeña, ello debido a la poca
inercia que tiene el termómetro, es decir, el menisco del mercurio se mueve tan pronto
como la temperatura cambia en el bulbo.

158
  El fluido termométrico no moja la pared, despreciándose los efectos de tensión
superficial.
 Se ignoran los efectos debidos a las fuerzas capilares en la columna termométrica.
 No hay condensación de los vapores del fluido termométrico en las paredes de la
columna capilar del termómetro.
 Como consideración final, en este modelo el fluido en el termómetro se expande a
presión constante dentro del capilar, lo cual implica suponer que se trabaja con el calor
específico a presión constante.

6.3 EQUIPO

El equipo consiste de:

1) Un termómetro de precisión: Puede utilizarse el termómetro de vidrio de un calorímetro.

Se recomienda el termómetro de referencia ASTM 56C – 86. Ver en las siguientes dos
figuras las especificaciones.

Figura 6.1 Termómetro ASTM 56C -86

Figura 6.2 Termómetro de la bomba calorimétrica

2) Un vaso de precipitados de 500ml o mayor capacidad.

3) Un mechero.
4) Soportes: universal, trípode y pinzas de aro de nuez.
5) Un cronómetro.

159
6) Una Cámara de video.

Figura 6.3 Elementos de grabación y edición

7) Un tablero adecuado para desplazar la cámara.

Tabla 6.3 Especificaciones del termómetro ASTM 56C – 86 [ASTM E1. (2007)]

160
6.4 PROCEDIMIENTO

El experimento consiste básicamente en tomar las lecturas de cambio de temperatura en

función del tiempo registradas por un termómetro de vidrio. Se considera como sistema que
transfiere calor con los alrededores, el bulbo del termómetro, el cual es asumido como un
cilindro.

El bulbo se calienta con una llama hasta llevar al termómetro a su

límite superior de temperatura. Se expone al medio ambiente
externo (mas frío, ) y una vez comienza el enfriamiento del
bulbo con el aire estancado, se registra el descenso de
temperatura con una video cámara que sigue el menisco de
mercurio de la columna hasta que prácticamente llega a equilibrio
térmico con el ambiente. Para evitar corrientes de aire, el
termómetro es encerrado en un recipiente cilíndrico de muchísimo
mayor diámetro. La temperatura ambiente y la humedad relativa
del momento del ensayo también son registradas. Se debe conocer la mayor cantidad posible
de características del termómetro. Tome las mediciones de las longitudes del bulbo con un
calibrador y mida el espesor de la pared apoyado con un estereoscopio.

Se recomienda realizar por lo menos dos registros fílmicos por cada experimento replicando las
mismas condiciones con el ánimo de aplicar la teoría de errores. La filmación de cada registro
se puede procesar para digitalizar en fotos por medios de software como, Vegas Pro, donde se
generan hasta 24 cuadros por segundo.

6.4.1 Ensayos preliminares cualitativos

Para estudios cualitativos se proponen los siguientes ensayos

6.4.1.1 Ensayo-1: Calentamiento del bulbo con agua a temperatura máxima permisible

Caliente agua y agítela en un vaso de precipitados de 500 mL, colocando el vaso sobre una
plancha eléctrica de calentamiento, hasta alcanzar la temperatura máxima del termómetro, que
es de 35C. Monitoree el ascenso de temperatura en el agua del vaso con otro termómetro.
Una vez alcanzada la temperatura de 35C se sumerja el bulbo del termómetro de precisión en

161
el vaso con agua caliente y observe como asciende el menisco en el termómetro; no agite el
agua durante el calentamiento del bulbo.

6.4.1.2 Ensayo-2: Calentamiento del bulbo con agua a menor temperatura

Repita todo el Ensayo-1 pero caliente inicialmente el agua, ya no hasta 35C, si no hasta una
temperatura menor, 25C. Compare los dos ensayos anteriores, ¿Qué observa respecto de la
velocidad con que se mueve el menisco?.

6.4.1.3 Ensayo-3 calentamiento del bulbo con aceite

Repita el ensayo anterior pero cambie el fluido en el vaso, por uno de mayor viscosidad como
aceite mineral. Anote la diferencia cualitativa que se presenta al cambiar de fluido de
calentamiento.

6.4.1.4 Ensayo-4: Enfriamiento del bulbo en agua y aceite a temperatura ambiente

Ahora registre el enfriamiento, por lo que debe calentar el bulbo del termómetro con un
mechero hasta alcanzar la temperatura máxima permisible en el termómetro, pero sin registrar
la temperatura; una vez se alcance esta máxima temperatura, inicia el registro fílmico y también
simultáneamente enfríe el bulbo al sumergirlo en un vaso con agua a temperatura ambiente.

Nota: Para calentar el bulbo acerque y retire el mechero al bulbo para que la temperatura no
supere la temperatura máxima de la escala del termómetro, de lo contrario podría dañar el
termómetro.

6.4.1.5 Ensayo-5: Enfriamiento del bulbo con aire a temperatura ambiente

En cambio de líquidos ahora utilice un gas, como el caso aire a temperatura ambiente.

6.4.2 Experimento Enfriamiento Cilindro Vertical

Convección libre alrededor del bulbo del termómetro orientado verticalmente en aire.

Figura 6.4 Orientación vertical del termómetro

162
Para este montaje se requiere de un soporte universal y pinzas de nuez para mantener el
termómetro aplomado verticalmente, y de otro soporte móvil para la video cámara con el fin de
guiar el seguimiento de la filmación del menisco.

- Realice el montaje que se muestra esquemáticamente en la Figura 6.5

Figura 6.5 Montaje para el termómetro vertical

- Cerrar la puerta del laboratorio para eliminar corrientes de aire que podrían influir en el
ensayo. Y garantizar que no se induzca movimiento de aire durante el ensayo.
- Acercar lentamente al bulbo del termómetro la llama del mechero y observar que el menisco
del mercurio en el termómetro alcance la temperatura máxima que se fije para el ensayo.
Esta temperatura no debe sobrepasar nunca la temperatura máxima de la escala del
termómetro, pues si esto ocurre el termómetro se descalibra y podría dañarse.

163
- Una vez alcanzada la temperatura máxima fijada, la cual se toma como temperatura inicial,
se retira el mechero y se coloca el vaso de tal forma que el bulbo del termómetro quede
dentro del vaso como muestra la figura.
- Al final del procedimiento anterior ya la cámara de video debe estar filmando.
- Al enfriarse el bulbo del termómetro, el menisco del mercurio desciende por lo cual mueva
la cámara sobre el soporte de la cámara para seguir el movimiento del menisco y así
registrar y grabar la evolución de la temperatura en el tiempo.
- Durante el ensayo si hay personas alrededor, es importante que éstas no caminen cerca
del ensayo pues al desplazarse se originan corrientes forzadas de aire y esto afecta el
experimento pues ya no se tendría convección libre.

6.4.3 Experimento Enfriamiento Cilindro Horizontal

Convección libre alrededor de un cilindro orientado horizontalmente en aire

Figura 6.6 Orientación horizontal del termómetro

 Realice el montaje que se muestra esquemáticamente en la figura 6.7

 Use una guía en forma de riel para la cámara.

Figura 6.7 Montaje para termómetro horizontal

164
- Cerrar la puerta del laboratorio para eliminar corrientes de aire que podrían influir en el
ensayo. Y garantizar que no se induzca movimiento de aire durante el ensayo.
- Acercar lentamente al bulbo del termómetro la llama del mechero y observar que el menisco
del mercurio en el termómetro alcance la temperatura máxima que se fije para el ensayo.
Esta temperatura no debe sobrepasar nunca la temperatura máxima de la escala del
termómetro, pues si esto ocurre el termómetro se descalibra y podría dañarse.
- Una vez alcanzada la temperatura máxima fijada, la cual se toma como temperatura inicial,
se retira el mechero y se coloca el vaso de tal forma que el bulbo del termómetro quede
dentro del vaso.
- Al final del procedimiento anterior ya la cámara de video debe estar filmando.
- Al enfriarse el bulbo del termómetro la columna de mercurio se contrae por lo cual se debe
mover la cámara sobre el riel para seguir el movimiento del menisco y así registrar y grabar
la evolución de la temperatura en el tiempo.
- Durante el ensayo si hay personas alrededor, es importante que éstas no caminen cerca
del ensayo pues al desplazarse se originan corrientes forzadas de aire y esto afecta el
experimento pues ya no se tendría convección libre.

6.4.4 Ensayos Para Calentamiento (Opcional)

Como práctica opcional se propone: Bajo el bulbo del termómetro vertical situar a una distancia
conveniente un bombillo de 100 vatios. Al encender el bombillo el calor radiado calentará el
bulbo. Y la temperatura comenzará a aumentar. Filmar el ascenso de la columna termométrica.

6.5 GRÁFICAS NECESARIAS PARA EL INFORME

Construya las siguientes gráficas:

 Evolución de temperatura junto con su réplica para un cada experimento (Temperatura-

tiempo).
 Evolución adimensional de temperatura junto con su réplica para cada experimento.
 Linealización de la temperatura adimensional respecto al tiempo adimensional.
 Evolución del número de Biot.
 Coeficiente de transferencia de Calor instantáneo.
 Número de Nusselt instantáneo.

165
6.6 CÁLCULOS Y REQUERIMIENTOS PARA EL INFORME

6.6.1 Parámetros

6.6.1.1 Propiedades Físicas

Las propiedades se calcularon teniendo en cuenta la temperatura promedio entre la inicial y

final, que se evalúa por el teorema integral del valor medio, y la integral a su vez se resuelve
mediante la Regla de Simpson ⁄ :

∆
...

Entonces la temperatura de película promedio con la que usualmente se encuentran las

propiedades del aire es:

Los valores de las variaciones de las propiedades con la temperatura se extrajeron de la

literatura [KRETIH], con los cuales se generaron (por regresión) las siguientes correlaciones en
SI (la temperatura en °C, excepto para el cp,  y  del aire que además estarán en K):

Mercurio:

. .

. . . .

para el intervalo: ° °

Aire:
. .
. .

para el intervalo: ° °

. . . .

para el intervalo:

166
con [CENGEL]:

. .

6.6.1.2 Condiciones de Trabajo

Las condiciones de presión atmosférica local y aceleración gravitacional son:

.
.

6.6.1.3 Especificaciones del Bulbo

Con la información suministrada sobre el termómetro y/o las mediciones realizadas al bulbo
complete la siguiente tabla las especificaciones para la envolvente cilíndrica o pared del bulbo:

, mm
, mm (espesor)
, mm
, mm
Material
, W/m°C
Tabla 6.4 Especificaciones de la pared del bulbo

Con estos valores calcule para el bulbo: el área lateral del bulbo, ; el volumen, y la longitud
característica .

6.6.2 Determinación del intervalo de temperatura a trabajar

Usando el concepto de constante de tiempo, como el tiempo requerido para lograr

aproximadamente al · % . % del cambio total de temperatura, que

167
corresponde a . De esta manera podemos plantear que la temperatura inferior del
intervalo, , está dada por:

Con ello se busca rechazar la zona de la curva que obedece a una respuesta estacionaria, ya
que de acuerdo a la teoría de control cuando se tiene un · % . % del cambio
transitorio prácticamente se estabiliza en un estado estacionario. Cuando se acerca al equilibrio
termodinámico el valor de los coeficientes instantáneos de transferencia de calor tiende a
disminuir y a ser cero, debido a que la diferencia de temperaturas también tiende a reducirse y
a ser cero. De esta manera se asegura el estado tansitorio, al trabajar en la región exponencial
definida por , .

Es decir que de la serie de datos obtenidos, correspondientes a la evolución de la temperatura

para cada experimento, sólo se toman desde hasta , lo que implica rechazar desde
hasta .

6.6.3 Determinación del Coeficiente de Transferencia de Calor y el Número de Nusselt

Con los datos experimentales de temperatura Vs. tiempo puede determinarse el coeficiente de
Transferencia de Calor y con este el Nusselt promedio, que será de esta manera:

Los métodos para lograr dichos coeficientes son dos:

6.6.3.1 Regresión por Mínimos Cuadrados

Aplicando el método de mínimos cuadrados puede obtenerse el coeficiente convectivo en

términos de datos experimentales:

∑
∑

esta última asume que el intercepto con el eje de la variable dependiente es cero, de lo
contrario se debe asumir que [CHAPRA]:

∑ ∑ ∑
∑ ∑

168
El índice de correlación, o coeficiente de determinación se determina mediante la fórmula
[CHAPRA]:

∑ ∑ ∑
∑ ∑ ∑ ∑

Aplicando lo anterior a todos los datos se obtiene se obtiene un coeficiente promedio para todo
el proceso de enfriamiento.

6.6.3.2 Derivación e integración numérica

Para cada pareja de datos , o punto de la curva de enfriamiento se determina

numéricamente la derivada, ⁄ , o cambio instantáneo de temperatura usando las formulas
de primera derivada por diferencias finitas con cinco puntos: hacia adelante para los dos
primeros puntos, las de hacia atrás para los dos últimos puntos y centradas para el resto de
puntos [CHAPRA]. Adaptando y aplicando dichas formulas para el caso:

Luego se calcula cada coeficiente de transferencia de calor instantáneo a través de la

definición:

Ahora por el teorema integral del valor medio se determina el valor promedio del coeficiente:

169
Dicha integral definida se resuelve por la Regla de Simpson de 1/3, la cual requiere un número
par de particiones. Si los datos arrojan número impar de particiones, se elimina la última pareja
de medidas.

El error de esta determinación numérica es:

Para la determinación de este error, se procede con la serie temporal de coeficientes de calor
instantáneos obtenidos anteriormente de la siguiente manera. Para cada tiempo se calcula
numéricamente la cuarta derivada usando las siguientes formulas adaptadas a la presente
situación [CHAPRA, MATHEWS]:

Para los tres primeros puntos se aplica la formula de diferencias divididas finitas hacia
adelante, para los tres últimos puntos la formula de diferencias divididas finitas hacia atrás y
para el resto de puntos la formula de diferencias divididas finitas centradas. Luego se ubica
tanto el máximo y como el mínimo de estos valores de 4ª derivada, y entre estos dos se toma la
4ª derivada de mayor valor absoluto, con la cual se calcula en la formula el error máximo
debido a la aplicación de la integración numérica.

170
6.6.3.4 Corrección por Resistencia de Pared

Corrija el coeficiente global (aire-vidrio) que incluye la resistencia térmica conductiva de la

pared del bulbo para hallar el coeficiente de transferencia de calor individual del aire, usando:

donde es la resistencia de pared correspondiente a la conducción en un cilindro hueco:

6.6.3.5 Valores Esperados

De acuerdo a la información suministrada en el Fundamento Teórico y de ser necesaria otra

complentaria, calcule los coeficientes esperados y los valores de los números de Nusselt.

6.6.3.6 Errores Relativos Porcentuales

Encuentre el error relativo de cada uno de los coeficientes encontrados experimentalmente con
respecto al valor esperado del modelo seleccionado que presente mayor confiabilidad.

6.7 RESULTADOS

Presente adecuadamente los resultados obtenidos en las tablas y gráficas que considere
pertinentes.

6.8 ANÁLISIS Y DISCUSIÓN DE RESULTADOS

Para desarrollar esto tenga en cuenta los cuestionamientos presentados en el numeral 6.9.

6.9 ASPECTOS, ITEMS Y PREGUNTAS DE DISCUSIÓN

¿Es indispensable el coeficiente de transferencia de calor para los cálculos de transferencia de

calor? Señale las ventajas y desventajas de su utilización en la ingeniería de hoy ¿Cuál es la
tendencia actual?

171
Elija el modelo teórico y experimental que mejor predice el coeficiente de transferencia de calor
por convección libre y externa tanto para un cilindro vertical como para uno horizontal. De sus
razones.

Analice si se obtendrá experimentalmente el mismo coeficiente de transferencia de calor si en

vez de enfriar, se calienta el bulbo, considerando el resto de condiciones idénticas.

Explique cómo son generados los gradientes de densidad para la convección libre y su relación
con las capa limite.

De acuerdo a los ensayos cualitativos explique las diferencias entre los coeficientes que se
obtendrían en cada uno con respecto a las modificaciones como temperatura inicial, fluido, etc.

Dibuje la forma de la capa límite para una superficie cilíndrica fría y caliente, tanto vertical como
horizontal, y de hecho discuta las diferencias entre el calentamiento y el enfriamiento para una
misma orientación. También realice un análisis comparativo de la capa límite de un cilindro
horizontal Vs. vertical ¿A qué atribuye las diferencias entre los coeficientes obtenidos para cada
una de estas configuraciones?

¿Cuáles pueden ser las posibles causas del error de la determinación experimental? ¿El(los)
valor(es) obtenido(s) es (son) aceptado(s) o rechazado(s)? De sus argumentos.

¿El tomar un recipiente para evitar la la turbulencia no altera la convección libre?

¿A qué otros fluidos es aplicable esta práctica? ¿ aceite, por ejemplo?

¿Qué mejoras sugiere para hacer la práctica más eficaz y con menor grado de incertidumbre?

¿Es posible diseñar y construir un dispositivo que mida directamente el coeficiente de

transferencia de calor? Explique sus razones.

6.10 CONCLUSIONES Y SUGERENCIAS

Con base en los objetivos planteados y el análisis de resultados realice este ítem.

6.11 REFERENCIAS

172
 7. CONCLUSIONES Y RECOMENDACIONES

El modelo planteado y resuelto para el funcionamiento del termómetro de vidrio permite

determinar experimentalmente los coeficientes de transferencia de calor bajo las restricciones
impuestas. Los resultados de los coeficientes obtenidos experimentalmente fueron validados
para el funcionamiento del termómetro en posición vertical y horizontal rodeado de aire
estancado, que representa la transferencia de calor por convección natural en aire (Pr=0.72
aprox.) alrededor de un cilindro vertical (flujo externo en el bulbo) para un régimen laminar.
Estos valores se alejan de los pronosticados por la literatura en promedio un 35% para los que
se encuentran orientados verticalmente y en un 9 % para los que están horizontalmente; cabe
notar que los modelos de la literatura también tienen un margen de incertidumbre.

En estas tesis se obtuvo el modelo (adimensional) simplificado del termómetro [3-7]:

o con variables dimensionales [3-8]:

Se concluye que el termómetro puede utilizarse para medir coeficientes locales de

transferencia de calor en sitios del fluido donde se cumplan las restricciones del modelo: Flujo
reptante (en este caso corresponde a movimiento convectivo natural lento), como por ejemplo
si se tiene una placa vertical inicialmente caliente que se está enfriando por convección natural
en aire. Sería necesario situar el bulbo del termómetro a una altura con respecto al borde
inferior de la placa en sitios donde sea muy lento el movimiento del aire, esta condición puede
darse en cerca al borde inferior de la placa, pues allí comienza a acelerarse el aire.

Se determinó teóricamente el coeficiente de transferencia para convección natural alrededor de

cilindros verticales usando el método integral que permite plantear un sistema de dos
ecuaciones diferenciales ordinarias, cuya solución se desarrolla mediante la técnica de la

173
perturbación. La expresión que permite encontrar el número de Nusselt promedio [2-16] y por
ende el coeficiente de transferencia de calor promedio es

La cual está en concordancia con las correlaciones de la literatura, aun cuando genera valores
de coeficiente levemente menores a los de las correlaciones empíricas debido principalmente a
que solamente se trabaja con dos términos en la soluciones por el método de perturbación. El
primer término de esta ecuación representa una pared plana, mientras que los siguientes
describen la curvatura.

Se desarrollaron ecuaciones específicas para el problema propuesto en esta tesis de

determinar el coeficiente de transferencia de calor para convección libre con movimiento
laminar lento (flujo reptante), obteniéndose así la ecuación [3.18] para flujo cuasiestacionario y
una para flujo transitorio ecuación [3.15] para el termómetro orientado verticalmente.

Para flujo cuasiestacionario:

Para flujo transitorio:

Tanto para flujo cuasiestacionario como para flujo transitorio no fue posible efectuar la
confrontación de los resultados de estas dos ecuaciones para calcular le número de Nusselt
debido a que se requiere el conocer el valor del espesor de la capa límite, , o del radio de la
capa límite, . No olvidar que estos dos se relacionan mediante: .

174
Del modelo de flujo transitorio, se esperaba que para cada instante de tiempo correspondiera
un número de Nusselt, es decir que durante el enfriamiento del bulbo el Nusselt dependiera de
el tiempo, pero al desarrollar el modelo y calcular el Nusselt se encontró que la variable tiempo
se eliminaba por lo cual el número de Nusselt obtenido es global para todo el sistema. Por ello
con los datos de temperatura versus tiempo se determinó experimentalmente un coeficiente de
transferencia de calor para el enfriamiento en el bulbo del termómetro.

Se encontró que cuando el flujo convectivo natural es muy lento (flujo reptante), el número de
Grashof se elimina de la ecuación del número de Nusselt. Entonces puede preguntarse en que
términos del número de Nusselt se representa el movimiento del fluido. Para el caso de flujo
transitorio figura un parámetro , el cual está incluido en la ecuación de cantidad de
movimiento:

que al expresarla en forma adimensional aparecerá el número de Grashof y al solucionarla, el

parámetro quedará en función del número de Grashof y también el radio de la capa límite.

Para flujo de calor cuasiestacionario, donde Nusselt está en función únicamente del radio de la
capa límite, este espesor queda expresado implícitamente en términos del número de Grashof
en la ecuación que permite calcular la velocidad axial:

| |

donde y son:

Para el cilindro horizontal en estado cuasiestacionario se dedujeron las siguientes dos

expresiones:

Caso de gradiente angular de temperatura constante [3-21]:

175
 √

Caso sin gradiente angular de temperatura [3-22]:

Como los resultados obtenidos con el termómetro de precisión fueron relativamente aceptables
se propone la creación de un instrumento de medición directa de coeficientes locales de
transferencia de calor, utilizando según las restricciones del modelo, que fundamente su
funcionamiento por ejemplo en la localización del bulbo del termómetro sobre una gran
superficie que se encuentre en contacto con un fluido. Dicho sistema de medición debe tomar
como entrada la evolución temporal de la temperatura, la que puede ser detectada y codificada
por software que digitalice cambios en los pixeles de las imágenes; estos valores traducidos a
números son tratados también por otro software (como por ejemplo MATLAB) a través de la
derivación e integración numérica, y de acuerdo a los parámetros del mismo sistema, debe
generar como salida el coeficiente de transferencia de calor. Es un primer estudio para más
adelante utilizar un termómetro de precisión como medidor local de coeficiente de transferencia
de calor dentro de un fluido que transfiera calor.

Por lo anterior, es factible realizar a mediano plazo un diseño básico y en detalle de dicho
artefacto. De hecho se tiene proyectado construir una unidad de medición experimental de
coeficientes de transferencia de calor que contenga los anteriores elementos y que sirva como
apoyo didáctico para el aprendizaje de la Transferencia de Calor a nivel universitario.

El procedimiento experimental descrito en el cuarto capítulo de este documento, además de las

recomendaciones efectuadas, sirve como base y guía para abordar distintas prácticas
relacionadas con fenómenos de transporte (calor, cantidad de movimiento y masa). Gracias a
la experiencia recogida en este trabajo se propone y estructura una primera versión de Práctica
de Laboratorio de Fenómenos de Transporte para la determinación de los coeficientes de
transferencia de calor por convección libre en aire, la cual se encuentra en el Capítulo 6 de este
trabajo. Así una aplicación inmediata es el diseño de una práctica académica de fenómenos de
transporte, específicamente en transferencia de calor: esta práctica servirá para ilustrar los

176
rasgos típicos de la transferencia de calor por convección, tanto cualitativa como
cuantitativamente, como observarse en los ensayos experimentales 1 a 5; además permite
realizar una práctica de flujo de calor en flujo transitorio para determinar el número de Biot y el
número de Nusselt.

Para seguir validando el modelo se propone realizar otros experimentos conservando la misma
estructura básica pero adicionando ciertas variaciones: si se usan además como fluidos agua y
aceite se barre un amplio intervalo de números de Prandtl; si se usan termómetros con otras
dimensiones de bulbo y rangos distintos se amplía la gama de números de Grashof trabajados;
si se usan además otras situaciones, distintas a la convección libre, como la forzada cruzada o
tangencial al bulbo, u orientaciones inclinadas del termómetro, o sumersión de éste en un fluido
liquido con agitación o también situándolo a cierta distancia de una fuente calefactora, etc. se
podrían obtener entonces correlaciones más generales de números de Nusselt respecto a
Rayleigh en convección libre o respecto a Reynolds y Prandtl en forzada.

Buscando un refinamiento de la parte experimental se requiere elaborar montajes propios para

la medición específica (que a futuro van a servir para el diseño del instrumento de medición) del
coeficiente de transferencia de calor, que estén apoyados por otros equipos como
interferómetros así como termopares, que conecten en línea la información medida a sistemas
computarizados como proponen algunos autores [POPIEL]. De la misma manera la necesidad
de trabajar en un recinto más controlado, usando cámaras o campanas, para evitar
interferencias o la incidencia de otras fuentes de turbulencia que pueden alterar el tipo y
régimen de la convección.

Es posible que, la relativamente alta desviación de los resultados experimentales obtenidos

cuando se orienta el termómetro verticalmente, se deba a una reducción muy notable de la
turbulencia natural producida por la interferencia del vaso alrededor del bulbo, debe haber una
relación entre las dimensiones del elemento (bulbo) con respecto al confinamiento (vaso); tal
vez lo mas recomendado sea colocar un recipiente de mucho mayor diámetro o hacer el
experimento dentro de una cámara o campana.

Para el estudio de ciertos fenómenos de transporte, como la Transferencia de Masa y de

Cantidad de Movimiento asociadas a la de Transferencia de Calor se sugiere introducir ciertas
modificaciones y/o enfoques al experimento realizado.

177
En la determinación de coeficientes de transferencia de masa: se calienta el bulbo en un medio
fluido para luego extraerlo y dejarlo enfriar “húmedo” en el ambiente, situación en la que
también se miden los perfiles temporales de temperatura; luego se realiza una comparación
con los perfiles obtenidos de un enfriamiento “seco”, pero con el resto de condiciones iguales.

En la determinación de perfiles de velocidad: se asocia la solución desacoplada del modelo de

flujo transitorio para el termómetro vertical [3-11], con la ecuación de cantidad de movimiento,
que junto con el perfil de temperaturas encontrado a través de la determinación experimental
del coeficiente, permite encontrar el parámetro, , [3-14] (véase Capitulo 3) que a su vez
conecta con la solución para el perfil de velocidades; en este aspecto es interesante observar
que los modelos obtenidos para los perfiles de evolución de la temperatura, tanto del fluido
termométrico como del fluido alrededor del bulbo, se encuentran perfectamente articulados y
concatenados, claro que intermediados por la resistencia térmica conductiva de la pared, que
de acuerdo a la condición puede o no despreciarse. De la misma manera las expresiones
teóricas halladas en el Capitulo 3 desde [3-15] hasta [3-22] se pueden asociar con los
resultados experimentales para encontrar el cambio temporal del espesor de la capa límite, ya
que dicho parámetro resultó siendo la incógnita que impidió llegar por estos modelos al valor
del número de Nusselt.

Si los valores de los coeficientes obtenidos con las mediciones experimentales del estado
transitorio no se corrigen por la inclusión de la resistencia conductiva de pared, se obtendrán
coeficientes globales vidrio-aire; el proceso de corrección mediante la expresión obtenida por el
balance en el circuito térmico, permite disminuir el error en aproximadamente un 1% y obtener
coeficientes individuales de fluido (aire) que excluyen por tanto la fase sólida (vidrio); en un
inicio se consideró erróneamente que la pared de vidrio del bulbo no era tan gruesa con
respecto al diámetro del mismo bulbo, pero la medición de su espesor con el estereoscopio lo
desmintió, sin embargo, al tener dicha pared un número de Biot menor de 0.1 se puede asumir
como un sistema de resistencia interna despreciable. Este valor de la pared tiene una mayor
influencia en el cálculo de la masa del fluido termométrico dentro del bulbo que es usada para
la determinación experimental de los coeficientes.

Para la confrontación-validación del método experimental y la compresión del fenómeno, se

hace necesario disponer de modelos teóricos de mayor complejidad, que se pueden construir a

178
través del uso de otros métodos de solución del sistema de ecuaciones de conservación; para
la convección libre en cilindros verticales usar el menor número de restricciones posibles en el
método diferencial en la búsqueda de una solución exacta por métodos numéricos con
diferencias y/o elementos finitos; mientras que para los horizontales la deducción del
coeficiente a través del método aproximado integral (usando la técnica de perturbación) y por
soluciones no similares; también la posibilidad de usar la función corriente en el método
diferencial que permite transformar el sistema de ecuaciones diferenciales parciales en un
sistema de ecuaciones diferenciales ordinarias no lineal que puede resolverse por ejemplo con
Simulink de MATLAB; quedó entonces pendiente la resolución numérica del sistema de
ecuaciones de conservación para lograr una solución exacta aunque no analítica.

La pertinencia de la convección libre alrededor de un cilindro circular es indiscutible, ya que el

papel de las cantidades físicas que describen este fenómeno es relevante en muchas
aplicaciones actuales y futuras de la ingeniería, como por ejemplo: las pérdidas de calor en
tuberías verticales, la carga de calor requerida para el enfriamiento o calentamiento de
elementos cilíndricos, en los que se puede contar los conductores en dispositivos electrónicos y
la barras de fisión de reactores nucleares. No en vano, se encuentra que muchas de las
investigaciones asociadas a este tema han sido promovidas y desarrolladas por la NASA, en su
versión anterior NACA de los años 50’s.

179
 REFERENCIAS

BIBLIOGRAFIA

AHMAD, S., ARIFIN, N. M. , NAZAR, R. and POP, I. Free Convection Boundary Layer Flow
Over Cylinders of Elliptic Cross Section with Constant Surface Heat Flux. European Journal of
Scientific Research, Vol.23 No.4, 2008. Pag. 613-625,

AMBROSINI D., PAOLETTI, D., and SPAGNOLO G.S. Study of free-convective onset on a
horizontal wire using speckle pattern interferometry, Int. J. Heat Mass Transfer 46 (22) (2003)
4145–4155.

ARNAS, O. A. On laminar free convection stagnation heat transfer from an isothermal cylinder
with internal sources-sinks. Appl. Sci. Res. 27.

ASTM E1. (2007). Standard Specification for ASTM Liquid-in-Glass Thermometers. Pág. 1-6,
22, 49, 50, 54, 58-64.

ATMANE, M.A., CHAN, V.S.S. and MURRAY, D.B. Natural convection around a horizontal
heated cylinder: the effects of vertical confinement, Int. J. Heat Mass Transfer 46 (19) (2003)
3661–3672.

ASHJAEE, M., AMIRI, M. and ROSTAMI J. A correlation for free convection heat transfer from
vertical wavy surfaces. Heat Mass Transfer 44:101–111, 2007.

AZIZ, and NA, T.Y. Improved Perturbation Solutions for Laminar Natural Convection on a
Vertical Cylinder. Warme and Stoffubertragung Springer-Verlag 1982. Pag 83-87.

BARKER, Adam R. and FFOWCS WILLIAMS, John E. Transient measurements of the heat
transfer coefficient in unsteady, turbulent pipe flow. Int. J. Heat Mass Transfer 43 (2000), Pag.
3197-3207.

BHATIA, J.C. Laminar Boundary Layer on a Cylinder With Symmetric Cross-Section Bounded
by Two Equal Circular Arcs-Part I.
nd
BEJAN, A. Convection Heat Transfer. 2 edition, John Wiley and Sons Inc., New York, 1995.

BEJAN, Adrian. Heat transfer. John Wiley, New York, 1993.

BEJAN, Adrian and KRAUS, Allan D. Heat Transfer Handbook. John Wiley and Sons Inc., 2003.
Pag. 466-469, 482-483, 533-548.

BETANCOURT, Grajales R. Transferencia Molecular de Calor, Masa y de Cantidad de

Movimiento. Universidad Nacional de Colombia, Sede Manizales, 2002. Pag. 285 y 286.

BIRD, R. B; STEWART, W. E. y LIGHTFOOT, E. N. Fenómenos de Transporte. Editorial

Reverté, 1982.

180
BLAIR, M.F. Heat Transfer, 105. 1983.

BOJKOVSKI, Jovan; DRNOVSEK, Janko and PUSNIK, Igor. Analysis of equalising blocks in
calibration of thermometers. Measurement 23 (1998), 145–150.

BUCHLIN, Jean-Marie. Natural and forced convective heat transfer on slender cylinders. Rev.
Gen. Therm. Elsevier (1998) 37, 653-660.

BURMEISTER, LOUIS. Convective Heat Transfer. John Willey and Sons. Second Edition.
Canada, 1993. Pag 400-405.

CARSON, James K.; WILLIX, Jim and NORTH, Mike F. Measurements of heat transfer
coefficients within convection ovens. Journal of Food Engineering 72 (2006), 293–301.

CEBECI, Tuncer, QASIM, J. and NA, T. Y. Free convective heat transfer from slender cylinders
subject to uniform wall heat flux. Letters in heat and mass transfer Vol. I, Pag. 159- 162, 1974.

ÇENGEL, Yunus A. Heat Transfer: A Practical Approach. 2nd edition, McGraw-Hill, New York,
2003. Pag.408-414, 510-514.

CHAND, Jagdish and VIR, Dharam. Natural convection heat transfer from horizontal cylinders.
Journal of Chemical Engineerng of Japan vol. 12 No. 3, 1979, Pag. 242-247.

CHAPRA, Steven C. Métodos Numéricos para Ingenieros. Quinta Edition, Mc. Graw-Hill,
México 2006. Pag. 668-679.

CHEESEWRIGHT, R. Turbulent Natural Convection from a Vertical Plane Surface. J. Heat

Transfer, vol. 90, p. 1, febrero 1968.

CHILDS, P.; GREENWOODS, J and LONG C. Review of temperature measurements. Rev. Sc.
Instrum. 71, 2959 (2000).

CHILDS, Peter. Heat flux measurement. Practical Temperature Measurement, 2001, Pages
324-361.

CHILDS, P.R.N. Advances in temperature measurement. Advances in Heat Transfer, Volume

36, 2003, Pages 111-181.

CHILDS, P.R.N. Thermocouples. Practical Temperature Measurement, 2001, Pages 98-144.

CHOHAN, R. K.. Forced convection from liquids to industrial thermometers in cross-flow.

DULUC M.C., XIN, S., and LEQUERE P. Transient natural convection and conjugate transients
around a line heat source, Int. J. Heat Mass Transfer 46 (2) (2003) 341–354.

181
EDE, A.J. Advances in Free Convection. In IRVINE, T. F. Advances in Heat Transfer, 1967.
USA, 1967. Pag. 1-62.

ERDOGDU, Ferruh. A review on simultaneous determination of thermal diffusivity and heat

transfer coefficient. Journal of Food Engineering Volume 86, Issue 3, June 2008, Pages 453-
459.

FLACK, R. D., and WITT, C. L. Velocity Measurements in Two Natural Convection Air Flows
Using a Laser Velocimeter. J. Heat Transfer, vol. 101, p. 256, 1979.

FOX, Robert W and McDONALD, Alan T. Introduction to Fluid Mechanics. 2a Edition, John
Wyley & Sons, 1978.

FROLOVSKAYA, O. A. Self-similar solutions of unsteady boundary layers. Journal of Applied

Mechanics and Technical Physics, Vol. 43, No. 1, pag. 54-58, 2002.

GOLDSTEIN, R.J.; IBELE, W.E.; PATANKAR, S.V.; SIMON, T.W.; KUEHN, T.H.;
STRYKOWSKI, P.J.; TAMMA, K.K.; HEBERLEIN, J.V.R.; DAVIDSON, J.H.; BISCHOF, J.;
KULACKI, F.A.; KORTSHAGEN, U.; GARRICK, S. and SRINIVASAN, V. Review Heat transfer
—A review of 2003 literature. International Journal of Heat and Mass Transfer 49 (2006), 451–
534.

GOLDSTEIN, Aaron S. A Computational Model for Teaching Free Convection. Chemical

Engineering Education, 2004. Pag. 272-278.

GORI, F., SERRANO, M. G., and WANG, Y. Natural Convection along a Vertical Thin Cylinder
with Uniform and Constant Wall Heat Flux. International Journal of Thermophysics, Vol. 27, No.
5, September 2006.

GOROBETS, V. G. Free-convective heat transfer from cylindrical surfaces with transverse

sectional finning. Journal of Engineering Physics and Thermophysics, Vol. 75, No. 6, 2002.

HAGISHIMA, Aya and TANIMOTO, Jun. Field measurements for estimating the convective heat
transfer coefficient at building surfaces. Building and Environment 38 (2003), 873 – 881.

HALDAR, S.C., Conjugate analysis of heat transfer from a horizontal insulated cylinder, Int.
Commun. Heat Mass Transfer 30 (1) (2003) 139–147.

HALDAR, S.C., KOCHHAR, G.S., MANOHAR, K. and SAHOO, R.K.. Numerical study of
laminar free convection about a horizontal cylinder with longitudinal fins of finite thickness.
International Journal of Thermal Sciences 46 (2007) 692–698.

HAUSER, W. Reglamentos de Prueba del PBT. Tomo 2, Termómetros de Vidrio con Líquidos.
Deutscher Eichverlag GmbH Berlin, 1983.
HEGGSAND, P.J. and BURNS D. Single-blow experimental prediction of heat transfer
coefficients: A comparison of four commonly used techniques. Engineering Applications of
Artificial Intelligence, Volume 2, Issue 1, March 1989, Pages 40-48.

182
HERMANN R. Heat Transfer by Free Convection from Horizontal Cylinders in Diatomic Gases.
NACA Technical Memorandum 1366, 1954.

HOLMAN J. P. Transferencia de Calor. 8ª Edición, Mc. Graw-Hill, 1998, España. Pag. 149-191,
229-226.

INCROPERA, F. P., DEWITT D. P., BERGMAN, T. and LAVINE, A. S. Fundamentals of Heat

and Mass Transfer. 6th Edition, John Wiley and Sons Inc., New York, 2007. Pag. 403, 579-580.

IOAN, Pop and INGHAM Derek B.. Convective Heat Transfer: Mathematical and Computational
Modelling. Elseiver Science, 2001. Pag. 209-244.

ISLAM, Md; RAISUL, WIJEYSUNDERA, N.E. and HO, J.C. Evaluation of heat and mass
transfer coefficients for falling-films on tubular absorbers. International Journal of Refrigeration
26 (2003), 197–204.

JABBAR, Abdul and KHALIFA, N. Natural convective heat transfer coeficient - a review I
Isolated vertical and horizontal surfaces. Energy Conversion and Management 42 (2001) 491-
504.

JALURIA, Yogesh. Natural Convection. Pag. 525-570.

JANNA, William S. Engineering Heat Transfer. Second Edition. CRC Press LLC, 2000. Pag.
406-413, 426-429.

KAISER, E. Dynamic measuring error correction of encapsulated auxiliary wall heat flux
sensors made of film resistance thermometers.

KARLEKAR, B.V. Y DESMOND, R.M. T Transferencia de calor. 2ª Edición. Nueva Editorial

Interamericana, México, 1985 Pág. 590-613.

KERN, DONALD Q. Procesos de Transferencia de Calor. México, 1965. Ed. CECSA

KHALIFA, Abdul-Jabbar N.. Natural convective heat transfer coefficient -a review I. Isolated
vertical and horizontal surfaces. Energy Conversion and Management 42 (2001), 491-504.

KMONL, Vladimir. The determination of unsteady heat flow to wall from the measurements of
surface temperature made with thin film resistance thermometers.

KOBUS, C.J. Utilizing disk thermistors to indirectly measure convective heat transfer
coefficients for forced, natural and combined (mixed) convection. Experimental Thermal and
Fluid Science 29 (2005) 659–669.
KOH, J.C.Y and PRICE, J.F. Laminar Free Convection from a Non-Isothermal Horizontal
Cylinder. Mathematics Research, 1963.

183
KOUSAR, Nabeela y LIAO ShiJun. Series solution of non-similarity natural convection
boundary-layer flows over permeable vertical surface. Science China: Physics, mechanics &
astronomy. February Vol. 53 No. 2: 360–368, 2010

KREITH, Frank y BOHN, Mark S. Principios de Transferencia de Calor. Thomson Learning.

Sexta Edición, 2001. Pág. 233-286, 301-342

KREYSZIG, Erwin. Matemáticas Avanzadas para Ingeniería. Limusa-Wilei. Tercera Edición,

México, 2003. Vol. I: pág. 260-276 y Vol. II: pág. 92-145.

KUMAR, M. and NAth, g. Unsteady Free Convection MHD Boundary Layer Flow Near a Three-
Dimensikonal Stagnation Point. Indian J. Pure Appl. Math, 17 (7): Pag. 957-968, 1986.

KUZNETSOV, V. V. and FROLOVSKAYA, O. A. Boundary layers in free convection. Journal of

Applied Mechanics and Technical Physics, Vol. 41, No. 3, 2000.

LARSEN, Poul S. and ARPACI, Vedat S. On the similarity solutions to laminar natural
convection boundary layers. Int. J. Heat and Mass Transfer. Vol. 29, No. 2, pp, 342-344, 1986

LOVEDAY, D. L. and TAKI, A. H. Convective heat transfer coefficients at a plane surface on a

full-scale building facade. lnt. J. Heat Mass Transfer. Vol. 39, No. 8, pp. 1729-1742, 1996.

MACHADO COTTA, Renato. Integral Transforms in Computational Heat and Fluid Flow. CRC
Press, Inc., 1993.

MAHFOUZ, F.M., Transient free convection from a horizontal cylinder placed in a micropolar
fluid, Heat Mass Transfer 39 (5–6) (2003) 455–462.

MAHFOUZ, F.M. and KOCABIYIK, S., Transient numerical simulation of buoyancy driven flow
adjacent to an elliptic tube, Int. J. Heat Fluid Flow 24 (6) (2003) 864–873.

MARTINI, William R and CHURCHILL, Stuart W. Natural Convection Inside a Horizontal

Cylinder. A.I.Ch.E. Journal, 1960, Pag. 251-257.

MARTYNENKO, Oleg G. and KHRAMTSOV, Pavel P. Free-Convective Heat Transfer With

Many Photographs of Flows and Heat Exchange. Springer, Berlin, 2005. Chapter 1: Basic
Statements and Equation of Free Convection, pag. 1-79. Chapter 3: Free Convection on Curved
Surfaces, pag. 219-246. Chapter 6: Nonstationary Processes in Free convection, pag. 393-401.

MATHEWS, John H. y FINK, D. Kurtis. Métodos Numéricos con MATLAB. 3ª Edición, Pearson,
Prentice Hall, Madrid, 2000. Pág.356.

MERK, J. H. and PRINS, J. A. Thermal Convection in Laminar Boundary Layers II y III. Appl
Sci. Res., Section A, Vol 4, 1953, pag. 195-221.

MERRICK DUSINBERRE, George. Heat Transfer Calculations by Finite Differences.

International Textbook Company, 1961.

184
MILLS, Anthony F. Transferencia de calor. México: Irwin, 1995. Pag. 305-307, 385-434, 488-
492.

MINKINA, Waldemar. Theoretical and experimental identification of the temperature sensor unit
step response non-linearity during air temperature measurement . Sensors and Actuators 78
(1999), 81–87.

MOHAMMED, Hussein A and SALMAN, Yasin K. Free and forced convection heat transfer in
the thermal entry region for laminar flow inside a circular cylinder horizontally oriented. Energy
Conversion and Management 48 (2007) 2185–2195.

MUNOZ-COBO J.L., CORBERAN, J.M., and CHIVA,S. Explicit formulas for laminar natural
convection heat transfer along vertical cylinders with power-law wall temperature distribution,
Heat Mass Transfer 39 (3) (2003) 215–222.

NA, Tsung-Yen. Effect of Wall Conduction on Natural Convection over a Vertical Slender Hollow
Circular Cylinder. Applied Scientific Research 54: Pag. 39-50, 1995.

NA, T.Y. and CHIOU, J.P. Turbulent Natural Convection Over a Slender Circular Cylinder.
Warme- und Stofftibertragung 14, 157-164, Springer Verlag, 1980.

NIEUWVELT, C.; BESSEM, J.M. and TRINES, G.R.M.. A rapid thermometer for measurement
in turbulent flow.

OOSTHUIZEN, P.H. and PAUL, J.T. Natural convective heat transfer from an isothermal vertical
cylinder with an exposed upper surface mounted on a flat adiabatic base with a flat adiabatic
surface above the cylinder. 5th European Thermal-Sciences Conference, The Netherlands,
2008.

OSTRACH, Simon. An Analysis of Laminar Free-Convection Flow and Heat Transfer About a
Flat Plate Parallel to the Direction of the Generating Body Force. Report 1111-National Advisory
Committee for Aeronautics, NACA, 1952. Pag. 63-69.

PEREZ CRUZ, Justo Roberto. La Termodinámica de Galileo a Gibbs. Fundación Canaria

Orotava de Historia de la Ciencia, Tenerife, 1998. Pág. 17-20.

POPIEL C.O. and WOJTKOWIAK J. Experiments on free convective transfer from side walls of
a vertical square cylinder en air. ELSEIVER, ScienceDirect, 2004. Pag. 1-8.

POPIEL C.O., WOJTKOWIAK J. and BOBER K. Laminar free convective heat transfer from
isothermal vertical slender cylinder. Experiments on fre convective transfer from side walls of a
vertical square cylinder en air. ELSEIVER, ScienceDirect, 2007. Pag. 607-613.

RANI, Hari Ponnamma and NYUNG, Chang Kim. A numerical study of the Dufour and Soret
effects on unsteady natural convection flow past an isothermal vertical cylinder. Korean J.
Chem. Eng., 26(4), 946-954, 2009.

185
RANI, H.P. Transient natural convection along a vertical cylinder with variable surface
temperature and mass diffusion. Heat and Mass Transfer 40: 67–732, 2003.

RANI, Hari Ponnamma and NYUNG, Chang Kim. Transient free convection flow over an
isothemal vertical cylinder with temperature dependent viscosity. Korean J. Chem. Eng., 25(1),
34-40 (2008).

RANI, H.P and GANESAN, P. Transient natural convection along vertical cylinder
with Heat and Mass transfer. Heat and Mass Transfer 33 (1998) 449-455, Springer-Verlag
1998.

RAY, AK. Natural Convection, Pag. 1-12.

REES, S.J.; McGUIRK, J.J. and HAVES, P. Numerical investigation of transient buoyant flow in
a room with a displacement ventilation and chilled ceiling system. International Journal of Heat
and Mass Transfer 44 (2001), 3067-3080.

RODI W. and MULAS, M. Engineering Turbulence Modelling and Experiments 6. Elseiver,

2005.

SHOEMAKER, David P. Experiments in physical chemistry. 5th ed. New York, McGraw-Hill,
1989. Pag. 113, 399-400.

SIGALES, B. Transferencia de Calor Técnica. Vol. 1. Editorial Reverté, S.A., Barcelona, 2002.
Pag. 30-35, 197-206, 508-514.

SMITH, Carlos A. y CORRIPIO, Armando B. Control Automático de Procesos. Teoría y

Práctica. 1ª Edición. Noriega Editores, Limusa. México, 1991. Pág. 96, 151, 161.

TAKHAR, H.S., CHAMKHA A.J.,and NATH, G. Natural convection on a vertical cylinder

embedded in a thermally stratified high-porosity medium. Int. J. Therm. Sci. 41 (2002) 83–93.

THOMAS, L.C. Heat Transfer – Professional Version, 2nd edition. Capstone Publishing Corp,
1999.

UZEL, Ali R.; EDWARDS, Richard J. and BUTTON, Bryan L. An expert system for convective
heat transfer measurements using a transient analysis.

VAN NESS H. C. and ABBOTT M. M. Thermodynamics. McGraw-Hill, 1972.

VLIET, G. C., and LIN C. K. An Experimental Study of Turbulent Natural Convection Boundary
Layers. J. Heat Transfer, vol. 91, p.517, 1969.

WEBER, Robert L. and WHITE, Marsh W. College Physics. McGraw Hill Book Company, Inc.,
1960.

186
WEDEKIND, G. L. and KOBUS, C. J. Technical notes: Predicting the average heat transfer
coefficient for an isothermal vertical circular disk with assisting and opposing combined forced
and natural convection. Int. J. Heat Mass Transfer. Vol. 39, No. 13, pag. 2843-2845, 1996.

WEDEKIND, G. L. and KOBUS, C. J Modeling the local and average heat transfer
coefficient for an isothermal vertical flat plate with assisting and opposing combined forced
and natural convection. Int. J. Heat Mass Transfer. Vol. 39, No. 13, pp. 2723-2733, 1996.

XUA, Hang, LIAO, Shi-Jun and POP, Ioan. Series solution of unsteady boundary layer flows of
non-Newtonian fluids near a forward stagnation point. J. Non-Newtonian Fluid Mech. 139 (2006)
31–43.

ZHANG, N., ZHENG, Z.C., and ECKELS S. Study of heat-transfer on the surface of a circular
cylinder in flow using an immersed-boundary method. International Journal of Heat and Fluid
Flow 29 (2008) 1558–1566.

ZILL, Dennis G. Ecuaciones Diferenciales con Aplicaciones de Modelado. Thomson Learning,

Séptima Edición. Colombia, 2002. Pág. 292-300.

ZUBAIR, Syed M. A Laboratory Manual for the ME 316: Thermo-Fluids Laboratory (Part b: Heat
Transfer). Mechanical Engineering Department. King Fahd University of Petroleum & Minerals
Dhahran 31261, Saudi Arabia.

INFOGRAFIA

http://asi-se-hace.seriespepito.com/capitulos-doceava-temporada-12/capitulo-3-termometros-
balanzas-comerciales-pintado-de-aviones-y-bombones-de-lujo/

www.adsabs.harvad.edu

www.deepblue.lib.umich.edu/handle

www.metas.com.mx/guiametas/La-Guia-MetAs-08-09-termometros-liquido-en-vidrio.pdf

www.sinab.unal.edu.co/contenido/recursos/bases.php

www.sinab.unal.edu.co/ntc/NTC4505.pdf

www.sinab.unal.edu.co:2048/login?user=mamayorgab&pass=__ &url=http://enterprise.astm.org

www.unidata.ucar.edu.co/staff/blynds/acerca.html

187
 APENDICE A

DEDUCCIÓN DEL NÚMERO DE NUSSELT PROMEDIO

PARA CONVECCIÓN NATURAL ALREDEDOR DE UN CILINDRO VERTICAL

I. PRINCIPIO DE SEMEJANZA

La resolución de tal sistema acoplado de manera analítica es muy difícil por su no linealidad.

Dado un cilindro vertical de altura, , y diámetro, . Se introducen las siguientes variables adimensionales:

Ecuaciones y Condiciones Límite Adimensionales

Adimensionalizando ahora las ecuaciones de conservación:

188
Las condiciones limite de contorno o frontera adimensionales de estas ecuaciones son:

 En la pared:
Para

∞ ∞ ∞ ⁄

⁄ ⁄ ⁄

para todo

se tiene

 En el seno del fluido:

En ∞
para todo
se tiene

Solución

El propósito principal al resolver las ecuaciones es obtener el flujo de calor transferido y para ello se requiere el coeficiente de transferencia de
calor a través del Número de Nusselt medio, .

ó ó

en forma adimensional, dado que

∞
⁄

189
de donde:

∞
⁄
⁄

se tiene

∞
⁄ ∞
⁄

⁄
⁄

⁄ ⁄
·
⁄ ⁄

Se observa que en el sistema de ecuaciones de conservación, en la ecuación de cantidad de movimiento aparece el parámetro adimensional,
número de Prandtl, , y en las condiciones limite aparece el número de Rayleigh, , pero como se tiene que , y lo que se
requiere contiene el número de Rayleigh, , por ende es correcto establecer que:

II. ANÁLISIS DE ESCALA

Para el análisis de escala se requiere que todas la variables adimensionales queden parametrizadas, lo que se logra eligiendo valores
característicos que garanticen que todas estas variables tengan un orden de magnitud de la unidad o menor. Este análisis indicara que términos
de las ecuaciones diferenciales no son indispensables, además que también indica algo sobre la estructura y forma de la solución.

Se definen entonces las variables adimensionales:

∞
; ; ; ;
∆ ∆⁄ ∞

Observe que la escala de surge de la ecuación de continuidad:

es la velocidad de escala y su magnitud apropiada se determina a partir de las ecuaciones. Se usa el espesor de la capa limite térmica, ∆, en la
escala de y en todas las ecuaciones, en lugar de usar el espesor de la capa limite hidrodinámica, , para las ecuaciones de continuidad y
cantidad de movimiento. Para , y ∆ tienen de hecho una magnitud similar, puesto que la fuerza de empuje produce un flujo en la capa
límite térmica, si bien la velocidad máxima ocurre cerca de la pared. Para , ∆ es mucho menor que , ya que las fuerzas viscosas ejercen
una fuerza de arrastre sobre el fluido que se encuentra fuera de la capa límite térmica, aunque es la fuerza de empuje la que impulsa al flujo, y
la velocidad máxima, , debe ocurrir para un valor de ∆. Usando ∆ en las ecuaciones de continuidad y cantidad de movimiento se garantiza
que ⁄ ,y ⁄ sean del orden de la unidad.

De acuerdo a la definición de las variables, se tiene que:

190
 ∆
; ∆ ; ; ; ∞ ∞

∆
; ∆ ; ; ; ∞

Las ecuaciones de cantidad de movimiento y de energía a escala quedan:

Ecuación de Cantidad de Movimiento:

Entonces:

Se sustituyen estas equivalencias:

∆
∆ ∞
∆ ∆ ∆ ∆

∞
∆ ∆

Ecuación de Energía:

De la misma manera

∆ ∞ ∞ ∞
∆
∆ ∆ ∆ ∆

∞ ∞
∆

Número Prandtl Grande,

Este caso es el límite cuando Prandtl es grande; para estos fluidos, ∆ ; la transferencia de calor está controlada por el interior de la capa
límite, donde las fuerzas de inercia (aceleración) son despreciables y por tanto las fuerzas de empuje son equilibradas por las fuerzas viscosas.

Por tanto en la ecuación adimensional de cantidad de movimiento se tiene:

~ ∞
∆

de donde:

∞
~
∆

Y por la ecuación de energía:

∞ ∞
~
∆

191
de donde:

∆ ~

Resolviendo las ecuaciones

∞
~
∆

∆ ~

con variables y ∆.

Multiplicando las ecuaciones:

∞
· ∆ ~ ·
∆

∞
~

⁄
∞
~
⁄

Reemplazando este valor de en la segunda ecuación (energía):

∆ ~

∆ ~

∆ ~ ⁄
∞
⁄

∆ ~ ⁄
·
∞
⁄ ·

∆ ~ ⁄
∞ ·

∆
~ ⁄
∞

∆
~ ⁄
∞

donde el grupo adimensional

192
es el número de Rayleigh, , por tanto:

∆
~ ⁄

Número Prandtl Pequeño,

Este caso es el límite cuando Prandtl es pequeño; la alta difusividad térmica de estos fluidos, hace que ∆ y por tanto la transferencia de
calor está controlada por el exterior de la capa límite, donde las fuerzas viscosas son muy pequeñas y por tanto las fuerzas de empuje son
equilibradas por las fuerzas inerciales.

Así de la ecuación de cantidad de movimiento:

~ ∞

de donde:

⁄
~ ∞

Que al sustituir en la de energía:

∆ ~

se obtiene:

∆ ~

∆ ~ ⁄
·
∞

⁄
∆
~
∞

∆
~ ⁄
∞

donde el grupo adimensional

es el número de Boussinesq, , por tanto:

∆
~ ⁄

Número de Prandt Unitario .

En este caso el Número de Prantl tiende a uno, o tiene orden de la unidad. Para fluidos en los que ∆ , tanto la fuerza inercial como viscosa
tienen una magnitud comparable respecto a la fuerza de empuje, la que está impulsando el flujo en todo momento.

193
De la ecuación de cantidad de movimiento,

~ ~ ∞
∆

que se puede desglosar en las dos ecuaciones:

~
∆

~ ∞
∆

Despejando de de:

~
∆

se obtiene:

~
∆

y reemplazando en

~ ∞
∆

se llega a:

· ~ ∞
∆ ∆

~ ∞
∆

∞
· ~ ·
∆

∞
~
∆

∆
Despejando :

∆
~
∞

∆
~ ⁄
∞

donde el grupo adimensional

194
es el número de Grashof, , por tanto:

∆
~ ⁄

Resultado

La magnitud del número de Prandtl determina si la convección libre se caracteriza por el número de Rayleigh, el de Boussinesq o el de Grashof.
Teniendo en cuenta que:

para el caso en que , los tres grupos adimensionales son equivalentes.

Determinando el Número de Nusselt

De acuerdo al balance en la frontera:

ó ó

en forma adimensional, dado que

∞
∆

se tiene

de esta manera,

~
∆

Como se define el número de Nusselt:

Entonces

~ ·
∆

195
 ~
∆

Se obtuvieron los valores de ∆⁄ para los tres casos de Prandtl, y por esto:

⁄ ⁄ ⁄ ⁄ ⁄

⁄
~
⁄ ⁄ ⁄

o de otra manera:

⁄ ⁄
~

Donde la constante se encuentra experimental o teóricamente, mientras que el exponente predicho teóricamente por el análisis de escala
asume los siguientes valores:

Exponente

⁄ ⁄

Para el caso del aire como fluido de trabajo . , el análisis lleva a:

⁄
~

Estas expresiones sirven para correlacionar datos experimentales. Tanto el análisis dimensional por el teorema Pi de Buckingham como el
principio de semejanza solo lleva al resultado , , mientras que el análisis de escala nos permite obtener la forma de esta función.

III. RESOLUCIÓN POR TEORÍA DE LA CAPA LÍMITE

 Sistema adimensional

Se suponen perfiles semejantes de velocidad y de temperatura

Si las dos capas límites se relacionan:

donde

De la ecuación de continuidad:

196
Reemplazando en la ecuación de movimiento junto con los perfiles:

Integrando dentro de la capa limite, respecto a ,

Definiendo:

La ecuación de movimiento queda:

Análogamente en la ecuación de energía,

De la ecuación de continuidad se tenía que

Introduciendo esto en la ecuación de energía junto a los perfiles queda:

Ahora si

Se tiene que:

Δ Δ

Integrando dentro de la capa limite :

197
 Δ Δ

Definiendo:

Así la ecuación de energía queda como:

E D

Se resuelve el sistema de ecuaciones diferenciales ordinarias considerando un perfil de capa límite:

E D

Combinando estas ecuaciones al dividirlas entre sí:

Una vez que se conoce Δ, se puede encontrar a partir de las ecuaciones anteriores y así saber el perfil de temperatura.

IV. MÉTODO INTEGRAL

Planteamiento del modelo

198
Ecuación de cantidad de movimiento

Aplicando la Segunda Ley de Newton:

Desarrollando lo correspondiente al cambio o flux de cantidad e movimiento:

Y por continuidad ,

Así:

Se encuentra la cantidad de movimiento que entra y que sale del elemento diferencial:

∆ ∆

Si

Y derivando respecto a r:

Donde es el espesor de la capa limite, y es el radio del cilindro.

∆ ∆ ∆

De la misma manera el flux de cantidad de movimiento que entra

o en términos generales:

199
Luego el flux de cantidad de movimiento neto sobre el elemento diferencial será:

donde

Ahora identificando todas las fuerzas externas que actúan sobre el volumen de control elemental en dirección z para encontrar la fuerza neta:

Fuerza Viscosa en la Pared:

| ∆ ∆

Donde es el área lateral del cilindro.

Fuerza de Flotación:

Peso:

Calculo de la fuerza neta:

, , ,

, ∞ ∆

, ∞ ∆

dado que:

derivando respecto al área

200
 ∆

y aplicando que

queda:

Además:

sustituyendo en la ecuación de fuerza neta

, ∞ ∆ ∆

, ∆ ∞ ∆

, ∆ ∞

Ahora ya conocidos cada uno de los términos, se aplica la segunda Ley de Newton:

, ∆ ∞
∆

y queda:

∆ ∞
∆

Dividendo entre ∆ y tomando el límite cuando ∆ :

∆ ∞
∆
∆ ∆

∆
lim lim ∞
∆

Por definición de derivada:

201
 ∞

Recordando que

de esta manera se tiene:

Anteriormente se había establecido la aproximación de Bousinesq (con en todos los términos excepto en los de la fuerza de empuje),
dividiendo entre :

donde:

∞
∞

Para finalmente obtener:

Ecuación de energía

Aplicando la Primera Ley de Termodinámica al elemento de control elemental:

h , h ,

Se considera que es flujo es estacionario, se desprecia la conducción en la dirección z, y la conducción en la dirección r es cero por fuera de la
capa limite. La transferencia de calor se da de la pared al volumen de control elemental por conducción.

Donde:

el sistema está en estado estacionario:

no hay intercambio de trabajo con los alrededores:

202
Existen dos corrientes de salida y una de entrada en las cuales se desprecían los cambios en sus energía potenciales y cinéticas porque el
elemento es muy pequeño.

h , hi , h |

h , h , h | Δ h |

Y por continuidad ,

Para la dirección z se tiene que:

h h

Donde es el espesor de la capa límite térmica

h h h h

Para la dirección r se tiene que:

h | h | h∞ ∆ | ∆ h∞ |

Ahora de la ecuación de continuidad en coordenadas cilíndricas:

Podemos establecer en función de al despejar:

Separando variables:

Integrando bajo la condición de frontera:

203
Calculado en

h | ∆ h∞ | ∆ h∞

h | ∆ ∞

intercambiando el orden de integración y diferenciación:

h | ∆ ∞

El flujo de calor en la pared hacia el fluido se da por conducción en la dirección radial:

∆ ∆

Así la primera Ley se puede reescribir como:

h , h ,

h , h ,

Reemplazando se obtiene que:

h | Δ h | h | ∆

dividiendo entre ∆ :

| Δ | | ∆
∆ ∆

Además sabiendo que

h | ∆ ∞

204
 h | Δ h | ∆ h∞
∆ ∆

y tomando el límite cuando ∆ :

h | Δ h |
lim ∞ lim
∆ ∆ ∆

h | Δ h |
lim ∞
∆ ∆

Por definición de derivada

h h∞

Y recordando que

entonces:

h h∞

h h∞

dividiendo entre :

h h∞

h h∞

Se sabe que a partir de:

se llega a:

Integrando definidamente con la condición:

∞ h h∞

h h∞
∞

como las propiedades se asumen constantes:

h h∞ ∞
∞

205
reemplazando este resultado tenemos:

Dividiendo entre

Por definición del coeficiente de difusividad térmica

Solución del modelo

Resolución simultánea de ecuaciones:

Resolviendo el sistema:

Conservación Cantidad de Movimiento:

Conservación Energía

Las condiciones de frontera para este caso son:

Se establece la relación constante (independiente de z) entre las capas límites hidrodinámica y térmica:

∆ ∆

Si se plantea que

∞ ∆

Donde se tiene que

206
 ∆

Se tiene que para la ecuación de movimiento:

Derivando esta última:

Y para la ecuación de energía:

∞ ∞

Derivando esta última:

Se halla para la ecuación de movimiento:

∞ ∞

∞ ∞ ∞

Para la ecuación de Energía se tiene:

∞ ∞ ∆ ∆

∞ ∞ ∆ ∆

207
 ∞ ∞ ∆ ∆

∞
∞ ∞
∆ ∆

Se proponen las siguientes formas para los perfiles de velocidad y temperatura:

Donde es una velocidad de escala. Estos perfiles están acordes con las condiciones de frontera y además son suaves en el borde de las
respectivas capas limites, es decir:

pero no se cumple que:

aunque estos perfiles en paredes planas han resultado sorprendentemente buenos [MILLS]

Llevando a las nuevas variables adimensionales estos perfiles son:

Sabiendo que

δT ∆δ

y como

Dividiendo estas dos últimas:

Tenemos que

208
 ∆ ∆

∆ ∆ ∆ ∆

∆ ∆ ∆ ∆

∆ ∆ ∆ ∆

∆ ∆ ∆ ∆

Para la de energía:

209
 ∆ ∆

∆ ∆

∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆ ∆

∆ ∆ ∆ ∆ ∆ ∆

|
∆ ∆ ∆

Sustituyendo estos resultados entonces,

La ecuación de cantidad de movimiento

con

210
 ∆ ∆ ∆ ∆

queda:

∆ ∆ ∆ ∆

Y la ecuación de conservación de energía:

con

∆ ∆ ∆ ∆ ∆ ∆

queda:

∆ ∆ ∆ ∆ ∆ ∆
∆

∆ ∆ ∆ ∆ ∆ ∆
∆

El sistema de ecuaciones se reduce a:

Movimiento:

∆ ∆ ∆ ∆

211
Y

Energía:

∆ ∆ ∆ ∆ ∆ ∆
∆

Suponiendo que los espesores de las capas térmica, , e hidrodinámica, , son próximos, es decir que las diferencias son despreciables,
, se tiene que ∆ , lo cual es valido para números de cercanos a 1. Por tanto si ∆ , el sistema queda:

Dividiendo cada ecuación por :

Este sistema de ecuaciones diferenciales ordinarias de primer orden con y como variables en función de se resuelve por el método de
perturbaciones, donde estas variables dependientes se expanden en series de potencias alrededor de ⁄ .

Esto se puede sustentar en el hecho de que este método se puede aplicar a cualquier ecuación diferencial que contenga un parámetro de
perturbación, , el cual podría aparecer en las condiciones de frontera. Dicho parámetro debe ser una constante positiva, tal que .

En este caso dicho parámetro es el reciproco del radio del cilindro, , el cual se encuentra en el sistema de ecuaciones y en el condición de
frontera, luego:

Que en las ecuaciones se observa:

Haciendo que ∞, se tiene que el sistema de ecuaciones queda:

Que son las ecuaciones que permiten encontrar a y para una placa vertical (KARLEKAR, MILLS).

212
Se observa que si el parámetro de perturbación tiende a cero, , se anula el efecto de la curvatura, lo que implica que, ∞. El mismo
resultado se puede obtener asumiendo que el espesor de la capa limite, , es pequeño comparado comparado con el radio, , de tal manera
que la relación / también tiende a cero y los mismos términos del sistema de ecuaciones de conservación se cancelan quedando el modelo
de la pared plana vertical.

Entonces se resuelve el sistema de ecuaciones diferenciales con el método de perturbaciones en / , lo cual requiere expresar la solución de la
siguiente manera:

Tomando solo los dos primeros términos , queda:

Esto es válido como una aproximación aceptable, ya que al aumentar cada término siguiente va a ser mucho más pequeño. Es de esperar que
el primer término, en , sea de una pared plana ya que carece del parámetro de curvatura, / .

Estos desarrollos de y en serie de potencias con dos términos exigen averiguar los valores de , , .Para lograr dicho
objetivo se procede a sustituir estas expresiones en el sistema de ecuaciones de conservación:

donde

Multiplicando la ecuación de movimiento por y la de energía por . Llamando , a la cantidad , es decir,

213
Reemplazando en la ecuación de cantidad de movimiento las expansiones:

214
Organizando de acuerdo a las potencias de / :

215
Se factoriza por potencias y como aproximación se desprecian las potencias de / , para las cuales

Igualando las respectivas potencias, se obtiene:

Simplificando,

De la misma manera, ahora se reemplaza la serie de potencias en la ecuación de energía:

216
Organizando por las potencias de / :

Se factoriza por potencias y como aproximación se desprecian las potencias de / , que corresponde a

Igualando las respectivas potencias, obtenemos:

217
 :

Simplificando,

Resolución del sistema para

Cantidad de Movimiento

Energía

Es claro que este sistema corresponde a una pared plana, que no contempla el efecto de la curvatura; como en ese caso comúnmente se hace,
se supone una solución de prueba que corresponde a variaciones de la ley de potencias para y , de esta forma:

Cuyas derivadas respectivamente son:

Donde , , y son constantes a determinar.

Sustituyendo estas soluciones de prueba en el sistema de orden cero para el parámetro de perturbación.

En la ecuación de movimiento:

218
En la ecuación de energía:

En estas ecuaciones la potencia de debe ser idéntica en todos los términos lo que lleva a tener un sistema de ecuaciones en el cual se puedan
hallar los exponentes y :

Que se resuelve satisfactoriamente para

Sustituyendo estos valores en las dos ecuaciones para encontrar los coeficientes y :

Momentum:

Energía:

· ·

Resolviendo el sistema:

219
Despejando en las dos ecuaciones e igualando:

Despejando :

Amplificando por / y recordando que

/
·
/

En forma aproximada:

Como se muestra en KARLEKAR

Ahora buscando el valor de en

220
 ·

O aproximadamente según Karlekar:

De esta manera las expresiones para y son

Donde

Así:

√ √ ·

√ √

221
Resolución del sistema para

Cantidad de Movimiento

Energía

Es claro que la resolución de este sistema ya será una primera aproximación al efecto de tener una superficie curva para el cilindro, de la misma
manera se supone como en el anterior caso una solución de prueba y en modelo potencial teniendo ya conocidas y

Cuyas derivadas respectivamente son:

Donde , , y son constantes a determinar.

Resolviendo el sistema de potencia cero se encuentra que:

Derivando:

Sustituyendo estas soluciones de prueba en el sistema de primer orden para el parámetro de perturbación.

En la ecuación de movimiento:

222
En la ecuación de energía:

La dependencia de se elimina en estas ecuaciones si los exponentes y si cumplen que:

Que se resuelve satisfactoriamente para

Sustituyendo estos valores en las dos ecuaciones para encontrar los coeficientes y :

Momentum:

· ·

223
Si tenemos:

Energía:

· ·

Resolviendo el sistema con incógnitas y :

Por la Regla de Cramer (determinantes):

∆ · ·

∆ ·

Si:

224
 √

Entonces:

∆ ·

∆ · ·

∆ ·

∆ · ·

∆ ·

∆ · ·

225
De esta manera:

∆
∆

∆
∆

Número de Nusselt

Un balance de calor en la interfase, , proporciona:

ó ó

. .

Despejando :

Del perfil parabólico de temperatura asumido tenemos:

Como se asumió que ∆ y despejando para derivar parcialmente respecto a la coordenada radial

Evaluando en :

226
Entonces el coeficiente de transferencia de calor queda:

Se había aproximado que

Donde se encontró:

Así:

Por tanto el coeficiente local de transferencia de calor será

con los valores de

como

√ √

√ √

√ √

227
 √
√

Por tanto el Nusselt local será

√
√

√
√

Para hallar ahora el Nusselt promedio requerimos encontrar el coeficiente promedio de transferencia de calor por el teorema integral del valor
medio, de esta forma:

Definiendo

Resolviendo la integral indefinida

Por sustitución:

De donde se despeja

Se divide por :

228
De la sustitución despejamos :

Así el diferencial queda:

Entonces sustituyendo

en la integral se tiene:

| |

Retomando el problema inicial:

· ·

Con

229
Si

donde

El Nusselt promedio será entonces:

una expresión más sencilla consiste en tomar el valor promedio del coeficiente de transferencia una aproximación para el logaritmo natural, y
retomando:

230
Si se tiene que

ln

Para | |

Aplicando esto a:

Si

Bajo este criterio se expresar que:

Aproximando a los cuatro primeros términos:

Reemplazando en la expresión del coeficiente promedio:

231
 √

Hallando el Nusselt promedio:

232
Si se expresa en función del diámetro, , amplificando el segundo término por 2

233
 APENDICE B

CONSTRUCCIÓN DEL MODELO PARA EL TERMÓMETRO: FLUIDO TERMOMÉTRICO

Se deduce v

De la siguiente manera,

dado que

h v

derivando parcialmente respecto a a constante se tiene:

h v

se define que

v v
v
v

Luego:

Despejando:

Balance de energía en el volumen de control

Balance de Energía aplicado a un volumen macroscópico de control (Robert W. Fox y Alan T. McDonald “Introduction to Mechanics” segunda
edición, 1978 John Wiley & Sons)

Que es la primera ley de la termodinámica donde se representan los siguientes términos:

ó : í
í í á

: ó

234
 :

: í

FLUJO DE CALOR,

El flujo de calor, es principalmente por convección ya que la radiación es despreciable porque el ∆ (entre la temperatura del medio ∞
y la temperatura de pared, es decir la de la superficie de control, ) es bajo.

Se tiene entonces que por la Ley de Enfriamiento de Newton:

Donde , es el coeficiente de transferencia de calor convectivo del medio circundante; es el área transversal al flujo, área que corresponde a
la de la superficie de control.

Se supone que la resistencia conductiva de la pared es despreciable y por tanto la temperatura del seno del mercurio, , es próxima a la
temperatura de la pared, , donde es la temperatura de pared.

Definiendo como el área de transferencia de calor, , la cual corresponde al área lateral del cilindro.

La Ley de Enfriamiento de Newton queda:

Ahora, si la resistencia de la pared no es despreciable, tendríamos que hallar el circuito equivalente entre las dos temperaturas:

Ahora como se considera que (El termómetro se va a calentar), para que se escribe:

FLUJO DE TRABAJO,

Siendo como no hay esfuerzos tangenciales sobre la superficie de control

Y el trabajo debido a piezas móviles:

Donde:

ya que no hay ejes que atraviesen la frontera del sistema.

ó
y ó

235
Luego:
ó

Finalmente, el único trabajo es el debido a la expansión:

CAMBIOS DE ENERGÍA DEL SISTEMA

Los cambios de energía potencial son despreciables comparados con los cambios de energía cinética y energía interna del volumen de control.

Y también

(El sistema está en estado permanente, ∆ ∆ )

Como no hay flujo de masa que atraviese la superficie de control, se tiene:

De esta manera tenemos:

Balance global: en todo el termómetro

Dado que el sistema, , se compone de dos subsistemas y :

si es constante como y es isotérmico y por ende en cada volumen de control

· ·

Entonces, como

236
 · ·

· ·

ya que y y como es un líquido: entonces:

· ·

donde

∞:

: ó

: í ó

: ó í

: ó

:

El término de la energía cinética es despreciable respecto al de energía interna. Se tiene que:

· ·

Así de esta manera:

· ·

Por lo que la ecuación de balance se escribe como:

Desarrollando la derivación:

237
Ecuación que describe el balance global.

Llevando esta ecuación en términos adimensionales:

Variables adimensionales:

donde : ó

donde

: ó ,

Relacionemos ahora con :

Derivando cada variable adimensional:

Luego:

Definiendo otra variable adimensional

Derivando:

238
Así se tiene:

Reemplazando estos diferenciales en la ecuación de balance global:

multiplicando por

Organizando:

Y definiendo los números adimensionales:

ó
í

Así la ecuación adimensional queda:

ahora se desglosa el factor en función de la altura de la columna:

Definiendo los subíndices:

Para el volumen, , del sistema se tiene:

239
y dado que para un cilindro

siendo : á ó y :

para que así:

y definiendo la condición inicial tanto para el bulbo como para la columna:

donde tanto como son invariantes con el tiempo, se despejan:

para reemplazar en la expresión de balance de volúmenes:

dividiendo entre queda:

Retomando la definición de volumen adimensional podemos establecer para la condición inicial:

donde

: ó ó

: ó ó

Y definiendo las longitudes fraccionales en el bulbo y en la columna:

respectivamente, para que finalmente la expresión de volumen adimensional en términos de las longitudes fraccionales:

Donde:

: , ó

: , ó

Si se deriva respecto a , el volumen adimensional, , sabiendo que son constantes en el tiempo, ya que el volumen del bulbo siempre
permanece lleno, es decir permanece invariante en el tiempo; así se obtiene:

240
definiendo la velocidad lineal adimensional del fluido en la columna:

queda:

así que:

Definiendo el volumen específico, v, como: v= por termodinámica se tiene:

v=v ,

la derivada total será:

v v
dv=

dividiendo entre v

v v v
=
v v v

Se tiene que el coeficiente de dilatación o expansión térmica, , se define como:

v
v

Y el coeficiente de compresibilidad, , como:

v
v

La expresión de volumen específico queda:

v
=
v

Integrando, teniendo en cuenta que y cambian muy poco en el proceso desde , , v hasta , , v ya que no se cumple que , se
asumen constantes.

v
=
v

v≈
v

v v
v ≈

v
≈
v

Como el proceso se asume isobárico tenemos que , lo cual se reduce a

v
≈
v

Despejando v:

241
 v
≈
v

v≈v

dado que v= yv =

·
≈

·
≈

·
≈

·
≈

·
≈

Recordando que

y de la misma manera: ,

y también definiendo el coeficiente de dilatación térmica adimensional, Κ, como:

Podemos escribir la expresión del volumen termométrico adimensional del fluido en términos de la temperatura adimensional como:

Por series de Taylor la función exponencial es:

! ! ! ! !

Así queda que:

Κ
≈
!

Κ Κ Κ
≈ Κ
! ! j!

Ahora si es un buen termómetro el volumen del fluido debe ser proporcional a los cambios de temperatura, es decir solo se toman en cuenta
los términos lineales, por tanto si para y .

≈ Κ Κ

242
 ≈ Κ

Derivando respecto al tiempo adimensional, , con Κ constante:

≈ Κ

≈ Κ

Pero anteriormente se había encontrado que

Igualando se obtiene

Así:

Que se reemplaza en la ecuación de balance total

Y esta queda:

Balance en la columna del termómetro

Los términos para la Primera Ley de la Termodinámica son:

FLUJO DE CALOR,

, ya que la columna se considera adiabática

FLUJO DE TRABAJO,

Dado que el volumen de la columna está determinado por:

243
Derivando respecto al tiempo:

Donde

Así el trabajo queda:

CAMBIOS DE ENERGÍA DEL SISTEMA

Teniendo en cuenta que los cambios de energía potencial y cinética del sistema con el tiempo son nulos porque es un sistema en estado
permanente.

CORRIENTES DE ENTRADA Y/O DE SALIDA

Ahora en la columna hay una corriente de entrada de materia como se muestra en la figura:

· v

Como se define la entalpía:

h= v

· h h h

Por continuidad tenemos:

ahora un balance de masa sobre la columna nos dice:

Entonces la continuidad nos queda:

244
Y entonces el término queda:

· h h h

Donde es la entalpia de referencia.

Así la Primera Ley sobre la columna queda:

Como es un líquido

Y por tanto

Volviendo adimensional y sabiendo que:

Y se tiene además:

Con y tenemos que y derivando:

Dividiendo:

Luego:

Reemplazando:

245
Si y dividiendo por :

Definiendo

Dado que

í
í

Y retomando

O de otra manera:

Retomando tal que

con

246
Se deduce que:

Balance de energía en el bulbo del termómetro

Resolviendo los términos de la 1ª Ley:

FLUJO DE CALOR,

FLUJO DE TRABAJO,

CAMBIOS DE ENERGÍA DEL SISTEMA

Teniendo en cuenta que los cambios de energía potencial y cinética del sistema con el tiempo son nulos porque es un sistema en estado
permanente.

CORRIENTES DE ENTRADA Y/O DE SALIDA

Ahora en la columna hay una corriente de salida de materia como se muestra en la figura:

· v

· h h h

Por continuidad:

ahora un balance de masa sobre el bulbo dice:

247
Entonces la continuidad queda:

Y entonces el término queda:

· h h h

Así la Primera Ley sobre la columna queda:

Volviendo adimensional sabiendo que:

Reemplazando:

Si

y dividiendo entre :

248
Si , y tenemos:

í
í

Se tienen las tres ecuaciones de balance

TERMÓMETRO: Κ

COLUMNA:

BULBO:

Se demuestra que:

BULBO + COLUMNA = TERMOMETRO

COLUMNA:

BULBO:

+__________________________________

dado que Κ

249
 Κ

TERMÓMETRO:

Que corresponde al balance en el termómetro.

Resolviendo la ecuación

Por separación de variables:

Asumiendo , Κ y constantes en el intervalo de operación del termómetro, para integrar teniendo en cuenta las condiciones iniciales: en
se tiene θ θ

θ
θ Κ

Despejando:

250
 θ
θ

Para el mercurio se tiene que para el intervalo de trabajo que

Entonces: Κ

Así:

θ
θ

Transformar en forma dimensional

, , ,

Despejando, :

Si :

Otra manera más simplificada:

Un balance de energía con base de tiempo sobre el sistema cerrado y en estado permanente correspondiente al fluido termométrico:

Donde:

251
Anteriormente se había establecido que el modelo de dilatación del fluido termométrico se aproximaba a:

dado el comportamiento lineal del termómetro, se puede realizar aproximar:

! ! ! ! !

entonces:
≈

Derivando esta expresión respecto al tiempo con constante:

Reemplazado en el balance de energía:

introduciendo las variables adimensionales:

Se tiente:

Y dado que

Dividiendo por :

252
De acuerdo a los grupos adimensionales:

Κ de donde

Así queda:

Κ
·

·Κ

También recordando el grupo adimensional

Queda:

Que es la misma ecuación obtenida anteriormente de balance general:

En forma más general:

Se desarrolla el balance anterior hasta este punto:

despejando :

253
 APENDICE C

CONVECCIÓN LIBRE ALREDEDOR DEL BULBO DEL TERMÓMETRO

I. TERMÓMETRO EN POSICIÓN VERTICAL CON FLUJO TRANSITORIO

Al observar el bulbo del termómetro, en la cual sucede el transporte de calor. Puede verse que el bulbo geométricamente está compuesto por
dos partes el cuerpo del bulbo de forma cilíndrica y el extremo del bulbo con formado por una semiesfera. Para efectos de modelamiento la
estrategia es inicialmente simplificar el problema asumiendo algunas restricciones. Y a medida que se plantea el modelo van asumiéndose otras
restricciones

Restricción geométrica: Asumir que el bulbo del termómetro es básicamente un cilindro.

Suposiciones:

Flujo transitorio de transporte de calor y de cantidad de movimiento.

Velocidad radial : No existe desplazamiento lateral de materia. A partir de la observación de fotografías del campo de flujo de velocidad
en cilindros calentados, se ve que básicamente el flujo solo tiene componente axial de velocidad.

No existe movimiento angular .

Gradiente de velocidad axial : Ello en virtud de las pequeñas dimensiones del bulbo, el recorrido que hace el aire en la longitud de bulbo

es pequeño.

Gradiente de velocidad radial : Es el más importante.

Se tienen en cuenta las fuerzas de flotación.

Va a utilizarse el método diferencial.

Al aplicar la Ecuación de c.mov. bajo las restricciones impuestas. El componente axial de la ecuación de cantidad de movimiento se reduce a:

Considerando que la velocidad del aire es pequeña alrededor del bulbo, según la literatura para números de Mach menores a 0.3 el gas puede
tratarse como fluido incompresible y para efecto en la ecuación de energía se considera propiedades medias constantes tales como la densidad
y la viscosidad. Bajo esta restricción y añadiendo la restricción de suponer que no existe gradiente axial de temperatura más las restricciones
dadas anteriormente, Al aplicar la ecuación de energía, ésta se reduce a:

Aquí se puede resolver la ecuación de energía desacoplada de la de cantidad de movimiento y con esta solución bajo la condición inicial y de
frontera se procede a resolver la ecuación de movimiento o flujo.

254
Resolviendo entonces la ecuación

Bajo las condiciones:

 Inicial:
,
 Frontera
,
Condición en la superficie del bulbo ó ó

∞, ∞

Solución

Por el método de separación de variables planteando una solución de la forma:

hallando las derivadas parciales respectivas:

Reemplazando en la ecuación de energía:

Dividiendo entre :

Si esto es igual a la constante, :

Se desprenden las dos ecuaciones diferenciales ordinarias:

255
Resolviendo cada ecuación independientemente, empezando con la de :

Separando las variables:

Integrando,

| |
Despejando :

ahora se prosigue con la segunda ecuación de :

A partir de la sustitución:

se tiene que:

y también

con lo cual se cambia de variable independiente en cada uno de los términos de la ecuación:

· ·

Sustituyendo en la ecuación

Dividiendo por :

256
La cual es una ecuación de Bessel con parámetro, , siendo un número real no negativo. La forma estándar de dicha ecuación es
[KREYSZING]:

Si , cuando toma el valor del entero cero, la ecuación se reduce a:

Que tiene como solución:

Donde
es la función de Bessel de primera clase de orden cero y tiene un valor de:
∞

y es la función de Besse de segunda clase de orden cero, que está dada por [KREYSZING]:
∞

donde es la llamada consten de Euler la cual se define como [KREYSZING]:

lim lim . …
∞ ∞

Aplicando al caso, la solución de

es:

Si

entonces:

por tanto el perfil espacio‐temporal de temperatura:

,
,
,

Llamando las nuevas constantes como

y

queda:

aplicando las condiciones iniciales y de frontera para hallar las constantes , :

 Espaciales (de frontera):

∞, ∞

∞, ∞ ∞ ∞

257
 ∞ ∞ ∞

· · ∞

· ∞

,
, ∞

∞ ,

A partir del perfil:

, ∞

Para realizar la anterior derivada hay que tener en cuenta que [ZILL]:

 Temporales (iniciales):
,
, ∞

Resolviendo el sistema para hallar y :

Por regla de Cramer:

258
∆
∆

∆ ∞
∞

∆ ∞
∞

así el valor de las constantes es

∆ ∞
∞
∆

∆ ∞
∞
∆

Por tanto la solución queda:

, ∞

Con

Sustituyendo las constantes en la ecuación se obtiene:

259
En (superficie del bulbo) la evolución en el tiempo es

Si ,

o en su forma logarítmica:

o despejando

La ecuación anterior, representa la evolución en el tiempo de la temperatura en la superficie del bulbo, obtenida a partir del perfil de
temperatura del fluido que rodea al bulbo. Ecuación que tiene la misma forma exponencial de la evolución de temperatura que se obtiene más
adelante (en Fluido Termométrico) dentro del bulbo:

Despreciando la resistencia de pared ,y y por comparación:

260
Donde es la difusividad térmica del fluido termométrico. No olvidar que:

Retomando:

donde el número de Fourier, , está dado por:

despreciando la resistencia de pared la constante tiene el valor de:

considerando la resistencia de la pared:

ó ó

Los perfiles arrojados:

Donde es un coeficiente global de transferencia de calor (en realidad U), ya que tiene incluida la resistencia de pared:

Igualmente:

Reemplazando:

261
Ahora si se considera que :

Despejando:

Dado que se tiene:

Por definición de Fourier del coeficiente global:

∆
Este coeficiente sería el coeficiente ,
Igualando:
∆
∆

De donde

Despejando :

Ahora como:
∆

Las resistencias presentes entre estas dos temperaturas:

262
no olvidar que la resistencia de pared, , es conductiva en un cilindro hueco:

luego el coeficiente global será:

también:

o también despejando :

así:

ó ó

se puede reescribir:

Se observa que si , se obtiene:

263
Esto puede suceder si se desprecia ante y y si además .

La constante surge durante el proceso de solución de la ecuación de energía, y ésta depende de la velocidad del fluido, como se observara al
sustituir el perfil de temperatura en la ecuación de cantidad de movimiento. El parámetro puede obtenerse de los datos experimentales de
temperatura en el bulbo y por mínimos cuadrados por la ecuación de la evolución de la temperatura.

De los datos de temperatura del termómetro puede determinarse y conocido se tiene el perfil de temperaturas con el cual se resuelve la
ecuación de cantidad de movimiento para determinar el perfil de velocidad:

Condición inicial

Condición en la superficie

Retomando el perfil de temperatura:

El parámetro queda vinculado a la velocidad cuando la solución anterior se reemplace en la ecuación de cantidad de movimiento .

La solución de esta ecuación permite relacionar la velocidad con el parámetro .

Aplicando la condición en , al perfil de temperatura se obtiene:

Finalmente el Nusselt para todo el sistema es:

Obsérvese que así el perfil de temperatura varia en el tiempo, la solución produce un solo número de Nusselt para el sistema que no depende
del tiempo, si el parámetro no depende del tiempo, lo que implica que la resistencia de pared es despreciable y que el coeficiente de
transferencia de calor por convección es constante a través del tiempo, lo que en realidad no sucede, pero que es salvable en la medida de que
se podría asumir el valor promedio como constante. Como no se conoce el espesor de la capa límite, , y dado que se obtuvo
experimentalmente un valor de y de tanto instantáneo como promedio, con estos resultados y a través de la ecuación anterior es posible
encontrar dicho espesor tanto instantáneo como promedio.

264
II. TERMÓMETRO EN POSICIÓN VERTICAL CON FLUJO CUASIESTACIONARIO

Durante los experimentos se observó que la temperatura en la escala del termómetro a partir de terminado instante varía muy lentamente.
Definiendo una escala de tiempo adecuada puede considerarse en este intervalo de tiempo que la derivada de la temperatura y velocidad con

respecto al tiempo no cambian apreciablemente y el problema puede tratarse como en estado estacionario . Manteniendo las

restricciones del modelo anterior más la consideración de Estado estacionario las ecuaciones de cantidad de movimiento y de energía se
reducen a:

Las condiciones para la ecuación de energía son

En

Considerando que , , para que esta ecuación de energía sea diferencial ordinaria y multiplicando dicha ecuación por :

Lo cual es una ecuación de Cauchy‐Euler [ZILL] de la forma:

Cuya solución está dada de la forma

Aplicando a la ecuación de energía, la solución tendrá la forma:

Derivando

y sustituyendo en la ecuación de energía:

265
de orden de multiplicidad 2.

Luego la solución es de la forma:

Bajo las condiciones de frontera:

En :

En

Despejando :

Luego para :

Reemplazando en el perfil:

La solución de la ecuación de energía bajo estas condiciones es también de la forma:

De la segunda condición se puede definir el número de Nusselt, definiendo como longitud característica el diámetro del bulbo del termómetro,
aunque se está en libertad de elegir la longitud del cilindro como longitud característica.

266
Se procede a sustituir el perfil de temperatura obtenido en la ecuación de cantidad de movimiento:

Para obtener:

Donde y b son:

Condiciones límite:
En
En

Solucionando la ecuación de c.mov junto con las condiciones límite

Es también una ecuación no homogénea de Cauchy‐Euler análoga a la de energía. Por lo tanto la solución de la homogénea:

es la función complementaria:

Ahora se encuentra la particular por variación de parámetros [ZILL]:

Donde

Aplicando el Wronskiano y los determinantes para determinar los parámetros . Se conoce que la función no
complementaria es:

Así

267
De esta manera:

Encontrando

Integrando por partes:

retomando:

268
De la misma forma para :

Así

Por tanto el perfil general de velocidad esta dado por:

Ajustando a las condiciones de frontera:

Obsérvese que si se utiliza la ecuación del perfil de temperaturas conociendo (el gradiente de temperatura en la superficie del bulbo)
entonces quedaría determinado el número de Nusselt. Esta constante puede determinarse al aplicar la condición en en la
ecuación del perfil de velocidad pues se obtiene el gradiente de temperatura en la superficie del bulbo, es decir .

Realizando lo anterior se obtiene, , al ajustar la condición de en el perfil de velocidad:

El recordar que y sustituirlo en la ecuación del perfil de temperatura, permite obtener finalmente el número de Nusselt:

269
Nota: las dos barras indican valor absoluto.

Es decir que Nusselt depende solo del espesor de la capa límite térmica.

III. TERMÓMETRO EN POSICIÓN HORIZONTAL CON FLUJO CUASIESTACIONARIO

Para plantear el modelo correspondiente al cilindro horizontal, es necesario considerar que este problema es más complejo que el de
convección en cilindro vertical, dado que en este caso, la fuerza de flotación al trabajar el problema en coordenadas cilíndricas no tiene una
sola componente sino, tiene dos componentes (la angular y la radial). Pese a que en coordenadas rectangulares la fuerza de flotación tiene una
componente, el problema se planteará en coordenadas cilíndricas pues al trabajarlo en coordenadas rectangulares el problema es más difícil de
plantear, y sobre todo de resolver.

Inicialmente se asumen las siguientes restricciones:

En coordenadas cilíndricas se tiene el vector velocidad:

Al observar en el interferograma de la visualización de las líneas de corriente del fluido en convección libre alrededor de un cilindro horizontal
se considera que cerca a la superficie del bulbo (dentro de la capa límite) el movimiento es principalmente angular pues el flujo tiende a seguir
el contorno del cilindro; despreciando por ello las componentes radial y axial de velocidad ante la angular, y por ello , para que el
vector velocidad se reduzca a:

Interferograma de cilindro enfriándose por convección natural [HOLMAN]

la fuerza de flotación por metro cúbico es

270
 ∞ ∞

Caso de gradiente angular de temperatura constante

Al aplicar la ecuación de cantidad de movimiento expuestas en el numeral 2.1.1 se obtiene:

Componente radial:

Componente angular:

Considerando gradiente angular de temperatura constante e ignorando la generación de calor por fricción viscosa y propiedades constantes
bajo el supuesto de flujo incompresible (que es válida si el Número de Mach es <0.3), la ecuación de energía expuesta anteriormente conduce
a:

Donde Gradiente de temperatura angular

Combinando las componentes angular y radial de cantidad de movimiento y utilizando la ecuación de energía, se obtiene la ecuación:

Condiciones en y en

Al Solucionar la ecuación anterior, en serie de potencias la distribución de velocidad angular tiene la forma:

Combinando la ecuación de energía con el perfil anterior de velocidad, se obtiene

271
 √

Bajo las condiciones límite en la superficie del bulbo

En
En

Solucionando y forzando condiciones límite se obtiene el perfil de temperatura

donde

Al derivar el perfil de temperatura y evaluarlo en

Al utilizar este resultado en la ecuación

se obtiene

Donde : Numero de Nusselt basado en el gradiente radial de temperatura y

es el Numero de Nusselt basado en el gradiente angular de temperatura.

Como una aproximación podría considerarse , lo cual conduce a:

Pero

272
Caso sin gradiente angular de temperatura

El caso anterior consideró gradiente angular constante de temperatura . Si ahora se considera la temperatura angular constante, es decir

el gradiente angular de temperatura igual a cero. Al repetir todo el procedimiento anterior, El número de Nusselt es

273
 APÉNDICE D
DATOS Y CALCULOS EXPERIMENTALES

I. EXPERIMENTO V-1

274
275
II. EXPERIMENTO V-2

276
277
III. EXPERIMENTO H-1

278
279
IV. EXPERIMENTO H-2

280
281
V. PARÁMETROS: ESPECIFICACIONES GEOMÉTRICAS Y PROPIEDADES

282
283
284