Transferencia de Calor Con Colectores Planos

Anlisis del mecanismo de
transferencia de calor en
colectores solares trmicos y de la
influencia de las condiciones de
contorno

Titulacin:
Ingeniero Industrial
Alumno/a:
Juan Francisco Cern Muoz
Director/a/s:
Jos Prez Garca
Ruth Herrero Martn

Cartagena, 23 de Septiembre de 2012

2

Contenido

ndice de figuras ........................................................................................................................................... 4
ndice de tablas ............................................................................................................................................ 8
CAPTULO 1: Introduccin ............................................................................................................................ 9
1. Inters y justificacin del proyecto ................................................................................................. 9
2. Antecedentes ................................................................................................................................ 10
3. Objetivos ....................................................................................................................................... 10
CAPTULO 2: Mecanismos de transferencia de calor en captadores solares trmicos. ............................. 11
1. Caractersticas geomtricas y constructivas de los captadores solares planos ............................ 11
1.1. Introduccin a los colectores solares planos ........................................................................ 11
1.2. Componentes de un captador solar plano ........................................................................... 11
1.3. reas caractersticas de un colector plano. .......................................................................... 14
1.4. Tipos de absorbedor ............................................................................................................. 14
2. Modelo analtico de captadores solares planos tipo parrilla ........................................................ 16
2.1. Ecuacin de balance de potencia ......................................................................................... 16
2.2. Distribucin de la temperatura en el colector ...................................................................... 16
2.3. Producto transmitancia-absortancia .................................................................................... 17
2.4. Coeficiente global de prdidas de calor del colector ........................................................... 18
2.5. Factores de eficiencia de las aletas y del colector ................................................................ 19
2.6. Factor de remocin de calor ................................................................................................. 22
2.7. Correlaciones para la transferencia de calor en tubo liso .................................................... 23
2.8. Ecuacin para las prdidas de presin ................................................................................. 23
3. Mecanismos de transferencia de calor en captadores solares ..................................................... 25
3.1. Transferencia de calor con el fluido caloportador en rgimen laminar ............................... 25
3.2. Conveccin Natural en Cavidades ........................................................................................ 37
3.3. Transmisin de calor por radiacin ...................................................................................... 44
CAPTULO 3: Modelos numricos implementados. Resultados obtenidos. ............................................... 60
1. Tubo horizontal con flujo de calor uniforme ................................................................................. 61
1.1. Definicin de la geometra y mallado en GAMBIT ................................................................ 61
1.2. Hiptesis de simulacin en FLUENT. ..................................................................................... 63
1.3. Metodologa de procesado de resultados ............................................................................ 69
1.4. Validacin del modelo y metodologa del procesado........................................................... 75
2. Modelo de tubo con absorbedor .................................................................................................. 82
2.2. Hiptesis de simulacin y condiciones de contorno ............................................................ 85
2.3. Resultados obtenidos ........................................................................................................... 88
3

3. Cavidad natural en la cavidad de aire entre absorbedor y cubierta de vidrio .............................. 93
3.1. Simulaciones en 2-dimensiones ........................................................................................... 93
4. Modelo completo sin radiacin para un slo conducto ................................................................ 99
4.2. Hiptesis de simulacin ...................................................................................................... 103
4.3. Anlisis de resultados ......................................................................................................... 109
5. Modelo completo con carga solar para un solo conducto .......................................................... 121
5.1. Introduccin a los modelos de radiacin solar ................................................................... 121
5.2. Modelo de carga solar ........................................................................................................ 134
5.4. Anlisis de resultados ............................................................................................................. 141
5.4.1. Anlisis de la transmisin de calor ..................................................................................... 142
5.4.2. Eficiencia del captador........................................................................................................ 145
6. Modelo completo 9 tubos con carga solar .................................................................................. 148
6.1. Definicin de la geometra y mallado en GAMBIT .............................................................. 148
6.3. Anlisis de resultados ......................................................................................................... 155
CAPTULO 4: Conclusiones y trabajos futuros .......................................................................................... 163
1. Conclusiones del proyecto .......................................................................................................... 163
2. Trabajos futuros .......................................................................................................................... 165
Anexo I: Propiedades de los materiales ................................................................................................... 166
1. Propiedades del agua .................................................................................................................. 169
2. Propiedades del aire .................................................................................................................... 169
3. Propiedades de materiales slidos .............................................................................................. 170

ndice de figuras
4

ndice de figuras
Captulo 2
Figura 2.1: Componentes de un captador solar plano ............................................................................... 11
Figura 2.2 reas caractersticas de un colector .......................................................................................... 14
Figura 2.3: Absorbedor de tipo serpentn .................................................................................................. 15
Figura 2.4: Absorbedor de tipo parrilla ...................................................................................................... 15
Figura 2.5: Esquema de un colector solar plano. Duffie and Beckman (1974). .......................................... 17
Figura 2.6: Distribucin de las temperaturas en un absorbedor plano. Duffie and Beckman (1974). ....... 17
Figura 2.7: Balance de energa en un colector solar plano. De Agbo y Okorolgwe (2007) ........................ 18
Figura 2.8: Geometra del tubo y de la aleta de un absorbedor. De Eisnman (2004). ............................... 19
Figura 2.9: Red trmica de una aleta .......................................................................................................... 20
Figura 2.10: Nodos del modelo .................................................................................................................. 21
Figura 2.11: Direccin de flujo en un colector solar de parrilla. ................................................................. 24
Figura 2.12: Desarrollo de la capa lmite hidrodinmica en un tubo ......................................................... 26
Figura 2. 13: Balance de fuerzas sobre un elemento diferencial para el flujo laminar completamente
desarrollado en un tubo circular ................................................................................................................ 27
Figura 2.14: Evolucin del perfil de temperatura del fluido en el flujo desarrollado con condicin de
contorno q
p
'' uniforme ............................................................................................................................... 29
Figura 2.15: Evolucin del perfil de temperatura del fluido en el flujo plenamente desarrollado con
condicin de contorno T
p
constante .......................................................................................................... 29
Figura 2.16: Conveccin forzada: solucin numrica del flujo para q
p
'' uniforme y T
p
constante. Nmeros
de Nusselt local y medio. De A. Garca (2006). ........................................................................................... 31
Figura 2.17: Conveccin forzada con flujo de calor constante. Nmero de Petukhow y Polyakov [1988]
frente a la solucin numrica. D A. Garca (2006). ..................................................................................... 32
Figura 2.18: Perfiles de flujo en conveccin mixta laminar en un tubo horizontal. Arriba: recirculaciones
simtricas en una seccin transversal; abajo: trayectoria en espiral del fluido a lo largo del tubo. .......... 33
Figura 2.19: Nmero de Nusselt de Petukhov y Polyakov [1988] en convenccin mixta .......................... 35
Figura 2.20: Desarrollo de la capa lmite sobre una placa vertical caliente ............................................... 38
Figura 2.21: Geometra del recinto ............................................................................................................ 40
Figura 2.22: Nmero de Rayleig Crticos para una cavidad horizontal ....................................................... 42
Figura 2.23: Comparacin de la ecuacin I con los resultados de Garon y Goldstein , Chu y Goldstein y
Rossby para agua (Pr=6) ............................................................................................................................. 43
Figura 2.27: Espectro electromagntico ..................................................................................................... 44
Figura 2.28: Distribucin de energa emitida por un cuerpo negro ........................................................... 46
Figura 2.29: Naturaleza direccional de la radiacin. (a) Emisin de radiacin desde un rea diferencial
dA1 en un ngulo slido dw subtendido por dAn en un punto sobre dA1. (b) Sistema de coordenadas
esfricas ...................................................................................................................................................... 47
Figura 2.30: Definicin de ngulos (a) plano y (b) slido ........................................................................... 48
Figura 2.31: ngulo slido subtendido por dAn en un punto sobre dA1 en el sistema de coordenadas
esfricas ...................................................................................................................................................... 49
Figura 2.32: Proyeccin de dA1 normal a la direccin de la radiacin ....................................................... 49
Figura 2.33: Emisin de un elemento diferencial de rea dA1 en un hemisferio hipottico centrado en un
punto sobre dA1 ......................................................................................................................................... 50
Figura 2.34: Comparacin de la emisin de un cuerpo negro y de una superficie real. Distribucin
espectral (a). Distribucin direccional (b) .................................................................................................. 52
Figura 2.35: Intercambio radiativo en un recinto isotrmico ..................................................................... 55
Figura 2.36: Reflexin difusa y especular ................................................................................................... 56
Figura 2.37: Sistemas de coordenadas para la funcin de reflexin .......................................................... 56
5

Figura 2.38: Representacin esquemtica de un experimento para la medida de la funcin de reflexin
.................................................................................................................................................................... 57
Figura 2.39: Superficie selectiva ideal ........................................................................................................ 59
Figura 2.40: Reflectividad vs longitud de onda para diversas superficies reales ....................................... 59
Figura 2.42: Ecuacin de transferencia de radiacin. Trminos. ............................................................. 122
Captulo 3
Figura 3.1: Geometra para el modelo del tubo horizontal ........................................................................ 61
Figura 3.2: Malla transversal en el tubo horizontal .................................................................................... 62
Figura 3.3: Algoritmo segregado. User's Manual, Chapter 25, FLUENT. .................................................... 63
Figura 3.4: Condicin de contorno de flujo de calor uniforme. FLUENT .................................................... 65
Figura 3.5: Propiedades en la entrada de flujo al dominio......................................................................... 67
Figura 3.6: Parmetros a definir en la salida del fluido .............................................................................. 68
Figura 3.7: Secuencia de superficies creadas a lo largo del conducto ........................................................ 71
Figura 3.8: Temperatura de la pared interior del conducto ....................................................................... 72
Figura 3.9: Nusselt en conveccin forzada. Datos obtenidos con mass-weighted average ....................... 73
Figura 3.10: Nusselt en conveccin forzada. Datos obtenidos con area-weighted average ...................... 73
Figura 3.11: Coeficiente de friccin frente a la distancia. Tubo liso horizontal con flujo de calor 1000
W/m2 y temperatura de entrada de 300 K. ............................................................................................... 74
Figura 3.12: Reynolds para los casos simulados para tubo liso horizontal con flujo de calor constante ... 75
Figura 3. 13: Resultados para la conveccin mixta en tubo liso con flujo de calor constante. Nmero de
Nusselt local, en funcin del Rayleigh y de la distancia adimensional x*. Comparacin con la correlacin
de Petukhov y Polyakov [1988]. ................................................................................................................. 76
Figura 3.14: Comparacin Nusselt en tubo horizontal con flujo de calor uniforme para alto nmero de
Reynolds. .................................................................................................................................................... 77
Figura 3.15: Error relativo para el nmero de Nusselt en tubo simple. ..................................................... 78
Figura 3.16: Nmero de Nusselt plenamente desarrollado, en funcin del nmero de Rayleigh Ra.
Comparacin con las correlaciones de Petukhov y Polyakov [1988] y A.Garca [2006] ............................. 78
Figura 3.17: Conducto horizontal liso, con flujo de calor uniforme, Ra=2.48E+05, Seccin z=0.9 m. Campo
de temperaturas (K) y vectores de velocidad. ............................................................................................ 79
Figura 3.19: Conducto horizontal liso, con flujo de calor uniforme, Ra=2.78E+05, Seccin z=1.80 m.
Campo de temperaturas (K) y vectores de velocidad. ............................................................................... 79
de temperaturas (K) y vectores de velocidad ............................................................................................. 80
Figura 3.25: Geometra modelo para 1 tubo simplificado .......................................................................... 82
Figura 3.26: Caractersticas geomtricas del tubo y aleta del absorbedor ................................................ 83
Figura 3.27: Seccin del plano XY, malla para fluido y conducto. .............................................................. 83
Figura 3.28: Detalle mallado de la soldadura. ............................................................................................ 84
Figura 3.29: Mallado del absorbedor. ........................................................................................................ 85
Figura 3.30: Condiciones de operacin para el modelo completo para 1 solo tubo con absorbedor. ...... 86
ndice de figuras
6

Figura 3.31: Condiciones trmicas sobre condiciones de contorno desacopladas (Original y Shadow) .... 87
Figura 3.32: Comparacin del nmero de Nusselt de tubo simple y tubo con absorbedor. ...................... 88
Figura 3. 33: Componente de la velocidad en Y. Tubo con absorbedor. Seccin 1.7m. ............................. 89
Figura 3.34: Componente de la velocidad en Y. Tubo simple. Seccin 1.7m. ............................................ 89
Figura 3. 35: Componente de la velocidad en X. Tubo con absorbedor. Seccin 1.7m. ............................ 89
Figura 3. 36: Componente de la velocidad en X. Tubo simple. Seccin 1.7m. ........................................... 89
Figura 3.37: Nmero de Nusselt desarrollado para el tubo con absorbedor. Comparacin con la
correlacin de Petukhov y Polyakov. .......................................................................................................... 90
Figura 3.38: Comparativa para el nmero de Nusselt con una condicin de contorno de flujo de calor
uniforme y una condicin de generacin interna de calor en el absorbedor, en conveccin forzada.
Reynolds 470-900 ....................................................................................................................................... 91
Figura 3. 39: Comparativa para el nmero de Nusselt con una condicin de contorno de flujo de calor
uniforme y una condicin de generacin interna de calor en el absorbedor, en conveccin mixta.
Reynolds 470-900. ...................................................................................................................................... 92
Figura 3.40: Malla 2D espacio de aire superior. Pruebas para la conveccin natural. Malla 1. ................. 94
Figura 3.41: Detalle de la malla junto a las paredes. Pruebas para la conveccin natural. ....................... 94
Figura 3.42: Seleccin de la temperatura de operacin y de la gravedad. ................................................ 95
Figura 3.43: Campo de temperaturas y velocidades. Malla 1.Conveccin natural. Ra=11000 .................. 97
Figura 3.44: Campo de temperaturas y velocidades. Malla 2.Conveccin natural. Ra=11000. ................. 97
Figura 3.45: Campo de temperaturas y velocidades. Malla 3. Conveccin natural. Ra=11000. ................ 98
Figura 3.46: Geometra del modelo completo para un tubo. ..................................................................... 99
Figura 3.47: Malla para el aire entre absorbedor y cubierta de vidrio. Modelo completo 1 solo tubo (sin
radiacin) .................................................................................................................................................. 101
Figura 3.48: Malla modelo completo 1 solo tubo (sin radiacin). Seccin transversal. ........................... 102
Figura 3.49: Vector y magnitud para definir el volumen completo.......................................................... 102
Figura 3.50: Condiciones de operacin para el modelo completo para 1 solo tubo (sin radiacin). ....... 103
Figura 3.51: Frontera superior del modelo. Condicin de contorno de tipo convectiva. ........................ 105
Figura 3.52: Fina capa de material creada sobre la misma condicin de contorno. ................................ 106
Figura 3.53: Fuente de generacin interna de calor para un flujo de 1000 W/m
2
. .................................. 107
Figura 3.54: Flujo de calor sobre la superficie interior del conducto. Radiacin absorbida 600 W/m2. . 110
Figura 3.55: Campo de temperaturas en fluido y tubera. Seccin 1.1 metros. ....................................... 110
Figura 3.56: Variacin del flujo de calor en la coordenada Y. Caso 600_Te_50 ....................................... 111
Figura 3.57: Variacin del flujo de calor en funcin de la coordenada Z. Caso 600_Te_50. .................... 112
Figura 3.58: Nmero de Nusselt frente a la distancia adimensional. Simulaciones para modelo completo
con 1 tubo................................................................................................................................................. 112
Figura 3.59: Resultados experimentales de transmisin de calor en conveccin forzada y mixta en tubo
liso de la Tesis de A.Garca [2006]. ........................................................................................................... 113
Figura 3. 60: Resultados numricos de transmisin de calor en conveccin mixta. Nmero de Nusselt
desarrollado en conveccin mixta para modelo completo. Comparacin con los resultados numricos en
tubo+absobedor ....................................................................................................................................... 113
Figura 3.61: Streamlines en la cavidad de aire. ........................................................................................ 115
Figura 3.62: Campo de temperaturas para las seccin 0.2 ,0.9, 1.2 y 1.6 m. Caso 600_Te_20................ 115
Figura 3.64: Nmero de Nusselt frente al nmero de Rayleigh para la conveccin natural en la cavidad
superior de aire. Comparacin con las predicciones de Hollands [1984]. ............................................... 116
Figura 3. 63: Campo de temperaturas para las seccin 0.2 ,0.9, 1.2 y 1.6 m. Caso 1000_Te_50. ............ 116
Figura 3.65: Distribucin de temperaturas en la seccin de salida del colector para un flujo de 1.36 g/s.
Validation of CFD simulation for flat plate solar energy collector, Selmi et al [2007].............................. 117
Figura 3.66: Campo de temperaturas para una generacin interna de 600 W/m
2
. Casos para temperatura
de entrada a 20 , 30, 40, 50, 60 y 70C .................................................................................................... 118
7

Figura 3.67: Eficiencia del captador para modelo completo de1 tubo sin radiacin. Comparacin con
banco de ensayos y modelo analtico adaptado. ..................................................................................... 120
Figura 3. 68: Flujo de calor por radiacin ................................................................................................. 122
Figura 3. 69: Condiciones de contorno en Paredes Opacas en el Modelo DO ......................................... 123
Figura 3. 70: Condiciones de contorno en paredes semi-transparentes del Modelo DO ........................ 125
Figura 3. 71: Condiciones de contorno en pared-semitransparente exterior en el Modelo DO .............. 125
Figura 3. 72: ngulos solares para una superficie vertical y horizontal. .................................................. 131
Figura 3.73: Introduccin de los parmetros de radiacin....................................................................... 135
Figura 3. 74: Lugar, fecha, hora y orientacin de la malla para las simulaciones. Solar Calculator. ........ 136
Figura 3.75: Superficies no irradiadas ...................................................................................................... 137
Figura 3.76: Carga Solar para paredes opacas. Definicin de la absortividad .......................................... 138
Figura 3. 77: Absortividad y transmisividad para una superficie semitransparente ................................ 138
Figura 3.78: Nmero de Nusselt para las simulaciones con colector completo 1 tubo con radiacin.
Comparacin con las correlaciones de Petukhov y Polyakov [1988] de conveccin mixta y forzada. ..... 142
Figura 3. 79: Nmero de Nusselt desarrollado para modelo de 1 tubo con radiacin. Comparacin con
los resultados numricos para el tubo con absorbedor. .......................................................................... 143
Figura 3. 80: Campos de temperatura y vectores de velocidad para la cavidad de aire entre absorbedor y
cubierta de vidrio. Caso completo 1 tubo con radiacin. ......................................................................... 144
Figura 3.81: Eficiencia del captador para modelo completo de1 tubo con radiacin para velocidades de
viento de 4 y 8 m/s. Comparacin con banco de ensayos. ...................................................................... 146
Figura 3.82: Opcin "Split Volume". GAMBIT ........................................................................................... 149
Figura 3. 83: Opcin "Bidireccional" en la herramienta "Split". GAMBIT ................................................. 149
Figura 3. 84: Geometra resultante despus de utilizar "Split Volume" ................................................... 150
Figura 3. 85: Geometra completa ............................................................................................................ 150
Figura 3. 86: Seccin transversal de la malla para el captador completo. ............................................... 151
Figura 3.87: Comparacin del nmero de Nusselt entre el conducto central y el conducto extremo para
un caso con temperatura de entrada de 20C, y gasto igua a 144 kg/h. ................................................. 156
Figura 3.88: Nmero de Nusselt desarrollado frente al nmero de Rayleigh en modelo de 9 tubos. ..... 157
Figura 3. 89: Nmero de Nusselt local frente a la distancia adimensional. Comparacin con las
correlaciones de Petukhov y Polyakov [1988]. ......................................................................................... 157
Figura 3.90: Distribucin de sondas de temperatura en el captador. ...................................................... 158
Figura 3.91: Comparacin del campo de temperaturas del tubo central para Q=144 kg/h. Datos
experimentales (derecha) y datos numricos (izquierda). ....................................................................... 160
Figura 3.92: Curvas de eficiencia para 216 kg/h, 144 kg/h y 72 kg/h ....................................................... 161
Figura 3.93: Comparacin de los resultados numricos para Q=144 kg/h, con los resultados
experimentales. ........................................................................................................................................ 162

ndice de tablas
8

ndice de tablas

Tabla 3.1: Caractersticas malla modelo tubo simple ................................................................................. 62
Tabla 3.2: Condiciones de contorno en tubo simple .................................................................................. 65
Tabla 3. 3: Simulaciones realizadas para tubo simple ................................................................................ 69
Tabla 3.4: Caractersticas malla de tubo con absorbedor .......................................................................... 85
Tabla 3.5: Condiciones de contorno para el modelo de tubo con absorbedor .......................................... 88
Tabla 3.6: N de nodos para mallas 2D en cavidad de aire. ....................................................................... 95
Tabla 3.7: Error cometido en cada malla para simulaciones 2D en cavidad de aire .................................. 96
Tabla 3.8: Volumenes incluidos en el modelo completo .......................................................................... 100
Tabla 3.9: Materiales en el modelo completo .......................................................................................... 104
Tabla 3.10: Condiciones de contorno para el modelo completo con un solo tubo .................................. 108
Tabla 3.11: Simulaciones realizadas para el modelo completo con un solo tubo .................................... 108
Tabla 3.12: Irradiacin solar extraterrestre .............................................................................................. 130
Tabla 3.13: Condiciones de contorno para el modelo completo con un solo tubo con carga solar......... 140
Tabla 3. 14: Materiales en el modelo completo de un solo tubo con carga solar .................................... 140
Tabla 3.15: Simulaciones realizadas en modelo completo con un solo tubo ........................................... 141
Tabla 3.16: Caractersticas de la malla en el modelo completo de nueve tubos ..................................... 151
Tabla 3. 17: Condiciones de contorno para el modelo completo con nueve tubos ................................. 154
Tabla 3.18: Simulaciones realizadas con el modelo completo de nueve tubos ....................................... 155

9

CAPTULO 1: Introduccin

1. Inters y justificacin del proyecto
El aprovechamiento de la energa solar como una energa renovable ha estado entre los intereses
principales de muchos ingenieros e investigadores a lo largo de estos dos ltimos siglos debido a sus
amplias aplicaciones, entre ellas la generacin de electricidad y la produccin de agua caliente sanitaria.
Hoy en da, los sistemas para produccin agua caliente sanitaria (ACS) se usan el sector residencial,
como el terciario y en establecimientos industriales. La produccin de ACS mediante la utilizacin de la
energa solar puede reducir los costes hasta en un 70%.
Desde la dcada de los 70 se le prest una atencin especial a este tipo de sistemas, sobre todo, a partir
de la crisis global de 1973. Desde entonces la eficiencia de los colectores solares se ha mejorado
progresivamente. Estas mejoras se pueden atribuir al uso de nuevos materiales y recubrimientos
selectivos, la mejora de los aislamientos trmicos y el desarrollo de tratamientos adecuados a las
cubierta de vidrio para mejorar el efecto invernadero en el interior del captador.
Debido a todas estas mejoras, en la actualidad, una de las pocas alternativas para incrementar la
eficiencia de los colectores solares planos pasa por mejorar la transferencia de calor entre el absorbedor
y el fluido calo-portador. Est puede mejorarse bien reduciendo la resistencia trmica de la soldadura
absorbedor-tubo o bien mejorando el coeficiente interno de transmisin de calor por conveccin. En el
primer caso existen serias limitaciones debido a que la lmina del absorbedor es de un espesor muy
pequeo y la tcnica de soldadura es compleja. En el segundo caso se puede conseguir una mejora de la
transferencia de calor mediante la utilizacin de dispositivos insertados tipo muelle o lminas en espiral.
La caracterizacin experimental de las mejoras implementadas en un captador solar requieren: fabricar
el prototipo, disponer de un banco de ensayos correctamente instrumentado, y de que se produzcan
durante un periodo de tiempo suficiente las condiciones climatolgicas adecuadas para que se pueda
obtener de forma fiable la curva de eficiencia del captador. Todo esto supone unos costes muy elevados
para los fabricantes.
La simulacin numrica es una herramienta que podra reducir costes de forma significativa en estas
fases de implementacin y prueba de las mejoras introducidas. Por todo ello, disponer de un modelo
validado resultara de gran utilidad y reducira los costes de implementacin y ensayo experimental de
las mejoras producidas.

CAPTULO 1: Introduccin
10

2. Antecedentes
Este proyecto final de carrera es continuacin de otro PFC anterior del mismo autor titulado Estudio de
la transmisin de calor en colectores trmicos mediante simulacin numrica, donde se realiz un
modelo simplificado de un solo tubo en el que se estudiaba la evolucin del nmero de Nusselt y el
coeficiente de friccin para una condicin de contorno aplicada en el absorbedor de temperatura
constante y un modelo ms complejo el que se el absorbedor se encontraba sometido directamente a
los efectos de la radiacin solar, pero no se incluan el resto de elementos del captador.
Todo el trabajo preliminar desarrollado en este primer PFC, en cuanto a metodologa de mallado y
conocimiento de los modelos de radiacin implementados en Fluent, as como las conclusiones
obtenidas han sido la base del trabajo realizado en este segundo PFC.
El PFC se enmarca dentro de la lnea de trabajo de Mejora de la eficiencia en captadores solares
trmicos del Grupo de Investigacin de Mecnica de Fluidos e Ingeniera Trmica. Los trabajos
realizados en materia de simulacin numrica, en la que se centra este PFC, se complementa con los
trabajos experimentales que se vienen desarrollando durante los ltimos cuatro aos en los dos bancos
de ensayo de los que se dispone, uno exterior, en el que se pueden ensayar dos captadores a escala real
simultneamente montados sobre un seguidor solar a dos ejes, y otro interior, en construccin, en el
que se podrn ensayar en condiciones mucho mas controladas tambin dos captadores solares, en este
caso a escala, formados por slo tres tubos. En esta instalacin se tiene previsto realizar ensayos de
caracterizacin del flujo para diferentes ngulos de inclinacin, diferentes fluidos, incluyendo
nanofluidos, etc.

3. Objetivos
El objetivo general del Proyecto Final de Carrera es desarrollar una metodologa de simulacin numrica
que permita cuantificar las mejoras introducidas en un captador solar trmico, sin necesidad de
desarrollar un prototipo y probarlo en banco de ensayos. La metodologa de simulacin consistir en:
- Desarrollo de la geometra a simular con Gambit. Se desarrollarn diversos modelos,
dependiendo de la complejidad y las variables que se quieran tratar, partiendo de un modelo
simplificado (absorbedor y conducto) hasta llegar un modelo completo (varios tubos,
absorbedor, aislamiento, aire intersticial y lmina de vidrio)
- Estudio de los modelos de transferencia de calor aplicables al problema en estudio
implementados en Fluent y definicin de las condiciones de contorno y propiedades de los
slidos y fluidos
- Mallado de la geometra.
- Ajuste y validacin de los modelos numricos implementados en Fluent y la metodologa de
procesado de resultados.
- Anlisis de resultados. Obtencin de correlaciones del nmero d Nusselt.
- Obtencin del campo de temperaturas del absorbedor y comparacin con datos
experimentales.
- Obtencin de la curva de eficiencia trmica del captador y comparacin con datos
experimentales.
11

CAPTULO 2: Mecanismos de transferencia de calor en captadores
solares trmicos.

1. Caractersticas geomtricas y constructivas de los captadores
solares planos

En este apartado se darn a conocer algunas nociones bsicas sobre colectores solares planos, as como
las ecuaciones de comportamiento de un colector solar plano de tipo parrilla.
1.1. Introduccin a los colectores solares planos
Estos colectores actan como un receptor que recogen la energa procedente del Sol, calentando una
placa, que calienta a su vez al fluido. Se utilizan en aplicaciones en las que se requiere que la energa se
libere a temperaturas bajas, como en la edificacin, el sector terciario, etc. Su temperatura de operacin
no suele pasar de los 100 C.
Se caracterizan adems, por usar tanto la radiacin solar directa como la difusa, por lo que no
necesitan seguir continuamente el recorrido del Sol. Adems, su mantenimiento no es nada
costoso y son de construccin simple.
Existen muchos diseos de colectores solares planos, profundizndose en algunos en apartados
posteriores. A continuacin se exponen las caractersticas y elementos propios de un colector
solar plano, esto es, sus componentes y sus reas caractersticas.
1.2. Componentes de un captador solar plano
1. Cubierta exterior
2. Carcasa
3. Absorbedor
4. Aislamiento
5. Aislamiento
6. Parrilla de tubos
7. Aislamiento lateral

1.2.1. Cubierta exterior
Su misin es producir efecto invernadero dentro del captador, reducir las prdidas por conveccin y
asegurar la estanqueidad del colector al agua y al aire, en unin con la carcasa y las juntas. El efecto
invernadero producido impide que toda la radiacin que la ha cruzado hasta alcanzar al absorbedor
salga al exterior, reteniendo la radiacin en su interior.
Una cubierta exterior debe cumplir las siguientes caractersticas:
Figura 2.1: Componentes de un captador solar plano
CAPTULO 2: Mecanismos de transferencia de calor en captadores solares trmicos.
12

- Alto coeficiente de transmisin de de la radiacin solar.
- Bajo coeficiente de transmisin para las ondas largas.
- Bajo coeficiente de conductividad trmica.
- Alto coeficiente de reflexin para la longitud de onda larga de la radiacin emitida por la
placa captadora.
- Coeficiente de dilatacin pequeo.
- Materiales plsticos.

Adems, la cubierta suele ser tratada, para mejorar sus prestaciones, mediante tratamientos anti-
reflectantes sobre la superficie exterior o interior de la cubierta.

1.2.2. Absorbedor
El elemento fundamental del captador. Su funcin es absorber la radiacin solar y transformarla
en energa trmica, que se transmite al fluido circulante por tubos.
Existen distintos modelos, entre ellos:
- Dos placas metlicas de cobre separadas por unos milmetros, entre las cuales circula el fluido
caloportador.
- Placa metlica de cobre sobre la cual estn soldados o embutidos los tubos por los que circula
el fluido caloportador. En lugar de una placa metlica se puede dotar de unas aletas de cobre a
los tubos de cobre.
- Dos laminas de metal de cobre unidas a gran presin excepto en los lugares que forman el
circuito del fluido calo portador, los cuales han sido abombados mediante insuflacin de aire.
- Lminas de plstico.
Las caractersticas e indicaciones de la placa absorbedora son las siguientes:
- Los tratamientos de la superficie mediante pinturas son ms econmicos que los tratamientos
selectivos pero tienen menos duracin.
- Si la instalacin va a funcionar mediante termosifn las prdidas de carga no deben ser
superiores a 3 mm.c.a por m
2
de superficie de captacin, para que la circulacin sea la
adecuada y no se produzcan grandes saltos trmicos.
- No se debe mezclar el cobre y el acero, para evitar la corrosin de este ltimo.
- El fluido caloportador que circula por la placa debe tener un reparto de temperaturas
equilibrado, es decir, el campo de temperaturas debe ser lo ms homogneo posible.
- En los paneles con doble placa, la transmisin de calor es directa, no ocurriendo lo
mismo para los que poseen los tubos soldados o embutidos. En este ltimo caso la
transferencia de calor va a depender de: la conductividad de la placa, la separacin, dimetro
y espesor de los tubos, el rendimiento y rgimen del flujo (laminar, transicin o turbulento)
y de la buena ejecucin de las soldaduras o de los acoplamientos a presin.
- Las entradas y salidas del fluido deben ser suaves de manera que las prdidas locales sean
bajas y que las soldaduras no estn forzadas para impedir posibles fugas.
- Se deben calorifugar bien las entradas y salidas para evitar prdidas importantes
debido a la creacin de puentes trmicos entre la placa y los elementos no aislados.
- Debe ser capaz de soportar la presin interna del fluido incluso en condiciones extremas.
- Se deben instalar elementos de filtrado adecuados para evitar obstrucciones.
Los materiales ms utilizados en su fabricacin son:
13

- Acero.
- Acero inoxidable.
- Aluminio.
- Cobre.
1.2.3. Aislamiento trmico
Es esencial su presencia, ya que reduce las perdidas trmicas. El aislamiento recubre la pared trasera y
los laterales del colector.
Los aislantes trmicos se caracterizan por:
- Resistir altas temperaturas sin deteriorarse, lo que muchas veces se consigue colocando entre
la placa y el aislante una capa reflectante, que impida que el aislante reciba directamente la
radiacin.
- No desprender vapores al descomponerse por el calor y en caso de ocurrir que no se adhieran a
la cubierta.
- No degradarse por el envejecimiento u otro fenmeno a la temperatura habitual de
trabajo.
- Soportar la humedad que se pueda producir en el interior de los paneles sin perder sus
cualidades.
- Poseer un espesor de entre 40 y 70 mm.
Los materiales ms usados son:
- Lana de vidrio y mineral.
- Espuma rgida de poliuretano.
- Poliestireno expandido.
1.2.4. Carcasa
Protege y soporta los elementos constituyentes del colector solar.
Cumple las siguientes caractersticas:
- Rigidez y resistencia estructural que asegure la estabilidad. Es de suma importancia ya
que debe resistir las cargas de viento.
- Resistencia de los elementos de fijacin mecnica para los esfuerzos a transmitir.
- Resistencia qumica a la corrosin y a los efectos corrosivos de la atmsfera y a la inestabilidad
qumica debido a las inclemencias del tiempo, ya que deben estar a la intemperie.
- Aireacin del interior del colector para evitar la condensacin del agua.
- Evitar toda geometra que permita la acumulacin de agua, hielo o nieve en el exterior del
colector.
- Permitir un fcil desmontaje de la cubierta para poder tener acceso a la placa captadora.

14

1.3. reas caractersticas de un colector plano.

Figura 2.2 reas caractersticas de un colector
- El rea total, A
G
, es el rea mxima proyectada por el captador completo. Se trata del rea
entre los lmites exteriores del captador, generalmente los bordes externos de la carcasa del
mismo.
- El rea de apertura, A
a
, se trata del rea proyectada mxima a travs de la cual penetra la
radiacin solar sin concentrar en el captador. Es la de la superficie visible o abierta del captador
para la radiacin solar. Por lo general coincide con el rea de la cubierta transparente visible
(sin contar la junta).
- rea del absorbedor, A
A
, se trata del rea mxima de proyeccin del absorbedor. En los
captadores solares planos es la suma de las reas de las aletas y de las tuberas internas de
distribucin expuestas a la radiacin.

1.4. Tipos de absorbedor
Se distinguen varios modelos de absorbedores en funcin de la configuracin de las tuberas que los
constituyen, destacando los absorbedores de parrilla y tubos y los de serpentn.
En cualquiera de las configuraciones se ha de asegurar:
- Buen equilibrio hidrulico para asegurar que el rendimiento del captador no disminuya. Para
ello el caudal ha de ser el mismo a travs del absorbedor.
- Circulacin del fluido en rgimen turbulento para favorecer la transferencia de calor. Aunque
no siempre se consigue este punto, lo normal es que el flujo del conducto sea laminar o de
transicin.

1.4.1. Absorbedor tipo serpentn
En los absorbedores de tipo serpentn todo el fluido de trabajo pasa por un solo tubo en forma de
serpentn. Debido a su configuracin presentan una mayor prdida de carga que el de tipo parrilla.
Las ventajas de este tipo de absorbedor son:
- La sencillez en fabricacin.
- La posibilidad de conexionado en paralelo de mucho captadores.
- La factibilidad de operacin de captadores individuales bajo condiciones de bajo flujo.

15

Figura 2.3: Absorbedor de tipo serpentn
1.4.2. Absorbedor de tipo parrilla
Se trata del tipo de colector a emplear en este proyecto.
A pesar de la mayor dificultad que requiere su fabricacin, el absorbedor tipo parrilla ha logrado una
cuota mayor en el mercado que el de serpentn, debido al uso extendido de las aletas absorbedoras.
Adems de producir menos prdidas de presin el incremento de temperaturas es menor y por tanto el
rendimiento de operacin es mayor.
En este diseo varios tubos finos en paralelo, unidos a tubos de distribucin en cada extremo del
captador, se combinan para formar una disposicin en forma de parrilla.

Figura 2.4: Absorbedor de tipo parrilla
Las ventajas que presenta este tipo de absorbedor son:
- La posibilidad de bajas prdidas de carga en captadores individuales. La prdida de carga
depende de los dimetros de los tubos empleados y del conexionado en paralelo o en serie.
- Posibilidad de uso en sistemas de termosifn.
- Versatilidad para el conexionado del campo de captadores.
- Mejor rendimiento debido a un incremento de temperaturas menor entre la entrada y la salida.

16

2. Modelo analtico de captadores solares planos tipo parrilla
Las ecuaciones que definen el modelo analtico de los captadores trmicos son:
2.1. Ecuacin de balance de potencia
El comportamiento de un colector solar puede determinarse mediante la ecuacin de balance de
potencia que describe la relacin entre la potencia til,
, la potencia absorbida por el colector, S, y las

prdidas trmicas y pticas.
La potencia absorbida, S, es igual a la radiacin solar incidente en la superficie del colector menos las
prdidas pticas.
Las prdidas trmicas de los alrededores se deben a las radiaciones infrarrojas, conveccin y
conduccin. Las prdidas trmicas dependen del coeficiente global de transmisin de calor, U
L
, a los
alrededores, de la temperatura ambiente, T
a
, y de la temperatura media del absorbedor,
abs
T .

La radiacin incidente global, G, est compuesta por las radiaciones directa y difusa, y las radiaciones
reflejadas de los alrededores y la tierra. La forma en que esta radiacin incidente es absorbida por el
colector depende de la radiacin incidente, la fraccin de la radiacin incidente,
g
t , que es transmitida
a travs de la cubierta transparente y la fraccin de radiacin incidente absorbida por el absorbedor,
abs
. La potencia absorbida puede ser definida como:

Sin embargo, en poder conocer el valor de la temperatura en el absorbedor es bastante
dificultoso, ya que depende del diseo del colector, la radiacin solar incidente, y las condiciones
del fluido que circula por dentro. Por ello, en el estudio de los colectores solares se introduce un
trmino, denominado factor de remocin del calor F
R
, para obtener una expresin con datos ms
asequibles; por lo que la ecuacin se reescribe de la siguiente manera:

A continuacin se detalla la obtencin de los distintos elementos de esta expresin. Sin embargo,
antes hay que analizar el comportamiento de la temperatura en un colector, con el fin de poder
establecer consideraciones que permitan obtener expresiones vlidas para el estudio del colector.

2.2. Distribucin de la temperatura en el colector
Para poder apreciar la distribucin de la temperatura en un colector, nos basaremos en la Figura 2.5 y la
Figura 2.6.

17

Figura 2.5: Esquema de un colector solar plano. Duffie and Beckman (1974).

Figura 2.6: Distribucin de las temperaturas en un absorbedor plano. Duffie and Beckman (1974).
El funcionamiento de un colector consiste en que la radiacin solar absorbida se dirija hacia los tubos
con el fluido caloportador. De esta manera, el paso de energa al tubo genera un gradiente de
temperaturas en el absorbedor en el sentido del flujo del fluido. Esto se observa en la Figura 2.5 y Figura
2.6. En la 2.4 se observa la radiacin incidente sobre el colector entre dos tubos; y en la 2.5 se observa
cmo vara la distribucin de temperaturas segn las coordenadas descritas en la figura a), en b)
se muestra la distribucin de temperaturas tridimensionalmente, en c) cmo evoluciona a lo largo de x
para una y fija; y en d), para una x fija, cmo evoluciona a lo largo de y.
2.3. Producto transmitancia-absortancia
Se define el producto global de transmitancia-absortancia como:
( ) ( ) ( )
G
G G G
g
abs g h g
s
abs g d
b
abs g b
abs g
o t
|
o t
|
o t
o t
|
.
|
\
|
+
|
.
|
\
| +
+
=
2
cos 1
2
cos 1
) (

Dicho producto relaciona la radiacin total absorbida por el captador con la radiacin total incidente.

18

Los productos (
abs g
o t )
b,
(
abs g
o t )
s
y (
abs g
o t )
g
se definen en funcin del ngulo de incidencia de la
radiacin y de las caractersticas constructivas del colector. En esta funcin el colector viene descrito por
un nmero de cristales, N
G
, cada uno con un ndice de refraccin n
g
y una longitud de extincin kL, y la
absortancia de la placa del colector
abs
.
2.4. Coeficiente global de prdidas de calor del colector

Las prdidas trmicas de cualquier sistema trmico solar presentan las tres formas de transferencia de
calor: radiacin, conveccin y conduccin.
- Las prdidas de calor por conduccin tienen lugar en los laterales y en la cara trasera del
absorbedor.
- Las prdidas de calor por conveccin tienen lugar desde el absorbedor hasta la cubierta de
cristal y pueden verse reducidas por el espacio de evacuacin entre el absorbedor y la cubierta
de cristal y por la optimizacin del hueco entre ellos.
- Las prdidas de radiacin ocurren desde el absorbedor, debido a la temperatura de ste.

Figura 2.7: Balance de energa en un colector solar plano. De Agbo y Okorolgwe (2007)
En la Figura 2.7 puede observarse el patrn de comportamiento que siguen las prdidas de calor en un
colector solar plano tpico.
Las prdidas de calor desde la cubierta transparente hasta el aire ambiente se deben a los intercambios
de radiacin y conveccin afectados por la velocidad del viento, la tierra, condiciones de los alrededores
y por la onda de radiacin desde el cielo.
El coeficiente global de prdidas de calor, U
L
, en W/m
2
K, es la suma de los coeficientes de prdidas por
la parte de arriba, laterales y parte trasera del colector.

Las prdidas en los laterales, si el diseo del colector es ptimo, son despreciables, pues en relacin
con la parte superior y la inferior es bastante pequea. Se suele agrupar con el coeficiente de la
parte trasera. Este, viene dado por la conductividad trmica y del espesor del aislante utilizado:

19

El coeficiente de prdidas para la parte superior del colecto es mucho ms compleja, pues tiene en
cuenta los fenmenos anteriormente descritos acerca de las prdidas de calor. Klein, en 1957, obtuvo
esta expresin para dicho coeficiente:

h
w
es el coeficiente de transferencia de calor debido al viento en la cubierta del colector (W/m
2
K) cuya
expresin viene dada por:

Siendo el parmetro u la velocidad del viento (m/s);
El resto de parmetros se han calculado mediante las siguientes expresiones:

2.5. Factores de eficiencia de las aletas y del colector

La geometra tpica del caso tratado sigue el esquema lmina-y-tubo de la Figura 2.8. El absorbedor
acta como una aleta, de anchura w, y se encuentra soldada al tubo por una zona de ancho g. La
transferencia de calor a travs de la aleta, puede ser descrita mediante el modelo de Duffie and
Beckman (1991) y de Lund (1989).

Figura 2.8: Geometra del tubo y de la aleta de un absorbedor. De Eisnman (2004).
Para poder conocer la distribucin de temperaturas entre dos tubos es necesario considerar
despreciable el gradiente de temperatura en la direccin del flujo. Tras realizar todo el clculo
diferencial para el estudio del gradiente de temperaturas, obtenemos una expresin para conocer el
calor recogido por un tubo a ambos lados por unidad de longitud es:

Siendo F la eficiencia de la aleta, expresndose:

20

El parmetro M puede expresarse mediante la ecuacin:

La potencia recogida en la superficie de la zona de contacto entre la aleta y el tubo por unidad de
longitud viene dada por:

La potencia til por unidad de longitud ser por tanto suma de la energa de la aleta y la energa de la
zona de contacto:

Esta potencia til es transferida al fluido, como sigue:

Donde F es el factor de eficiencia local del colector para un elemento dl y es constante para cualquier
geometra de colector y flujo msico. Su expresin puede obtenerse aplicando las leyes de Kirchhoff a la
red trmica de la Figura 2.9, teniendo en cuenta la absorcin de la radiacin solar en la base de la aleta:

El factor de eficiencia local del colector F representa la relacin entre la resistencia trmica entre el
colector y el ambiente y la resistencia trmica entre el fluido de trabajo y el ambiente.

Figura 2.9: Red trmica de una aleta
El parmetro k
gF
de la expresin (2.17) es del coeficiente de transmisin de calor entre el fluido y la zona
de contacto, en W/m
2
K, y viene dado por (obsrvese la Figura 2.10):

21

Figura 2.10: Nodos del modelo

El nmero de Nusselt, Nu, vendr dado por una correlacin adecuada al tipo de flujo (laminar, de
transicin o turbulento). El parmetro c se conoce como parmetro de conductancia en el contacto,
cuya expresin es:

y el parmetro
es el factor de forma, que viene dado por:

El parmetro k
gd
de la ecuacin (2.19) se trata del coeficiente de transmisin del calor entre la zona de
contacto y el tubo, en W/m
2
-K. Para un contacto metlico k
gd
es excesivamente alto, por lo que tendr
poca influencia en el factor F. Puede por tanto despreciarse el trmino 1/c de la ecuacin (2.18).
En la expresin (2.20),
viene dado por:

Donde
d
se conoce como parmetro tipo Biot, cuya expresin es:

Siendo Bi el nmero de Biot, que viene dado por:

h
f
es el coeficiente de transferencia de calor del fluido (W/m
2
K), que puede relacionarse con el nmero
de Nusselt mediante la siguiente relacin:

22

Sustituyendo la ecuacin (2.24) en la (2.23), se llega a:

2.6. Factor de remocin de calor

La relacin entre el calor til del colector y el calor til que poseera si la superficie completa
del absorbedor estuviese a la temperatura del fluido de entrada se denomina factor de remocin de
calor. Su expresin se obtiene a partir de la ecuacin (2.17). Esta ecuacin ofrece el calor por
longitud de colector, as que para definir el calor en todo su conjunto, es necesario realizar una integral
a lo largo del colector. Para ello, adems de la ecuacin (2.17), contamos con que el balance de potencia
para un fluido de trabajo desde la entrada a la salida del tubo es:

Para un elemento diferencial de longitud quedara:

Y se aplica un balance de energa para ste volumen de control:

Separando variables e integrando desde y = 0 (T
f
=T
in
) hasta y = L (T
f
=T
e
), se llega a:

La diferencia de temperatura del fluido puede ser aislada aadiendo y restando T
in
en el numerador de
la ecuacin (2.29):

Igualando las ecuaciones (2.27) y (2.3):

Sustituyendo con (2.31) y despejando F
R
, obtenemos:

Podemos obtener las temperaturas de salida y del absorbedor, combinando las ecuaciones (2.32) y
(2.1), respectivamente:

23

2.7. Correlaciones para la transferencia de calor en tubo liso

La transferencia de calor viene dada por el nmero de Nusselt y depende de Reynolds y Prandtl. Las
expresiones que definen estos nmeros adimensionales se presentan a continuacin:

El nmero de Nusselt expresa la relacin entre el calor transferido por conveccin a travs del fluido y el
que se transferira si slo existiese conduccin. Es funcin del coeficiente de transferencia de calor del
fluido h
f
, en W/m
2
K; el dimetro del tubo D
i
, en m; y de la conductividad trmica del fluido k
f
, en W/mK.
El nmero de Prandtl depende exclusivamente de las propiedades fsicas del fluido, siendo la
viscosidad dinmica del fluido, en kg/sm; c
f
el calor especfico del fluido, en J/kgK; y k
f
la conductividad
trmica del fluido en las paredes, en W/mK. Relaciona la velocidad de difusin viscosa y la velocidad de
difusin trmica.
El nmero de Reynolds representa la relacin entre las fuerzas de inercia y las fuerzas viscosas. Es
funcin de la viscosidad dinmica , en kg/sm; el flujo msico del fluido m , en kg/s; y el dimetro
interior del tubo D
i
, en m.
En rgimen laminar la transmisin de calor puede estar determinada por los efecto de flotacin debidos
a gradientes trmicos en el fluido. Esto tiene como efecto la creacin de corrientes secundarias en el
fluido que benefician la transmisin de calor en un proceso de calentamiento (que ms tarde se
analizarn). Para evaluar en qu medida la transmisin de calor est influenciada por la
conveccin natura, se cuenta con el nmero de Richardson Ri (Gr/Re
2
), que se trata del ratio entre
las fuerzas de flotacin (evaluadas por el nmero de Grashof Gr), y las fuerzas convectivas evaluadas por
el nmero de Reynolds. Las correlaciones para Nusselt se discutirn ms adelante.
2.8. Ecuacin para las prdidas de presin

Para estimar las prdidas de presin de un colector solar de tipo parrilla, ha de tenerse en cuenta el
diseo de partida.

24

Figura 2.11: Direccin de flujo en un colector solar de parrilla.
Como se muestra en la Figura 2.11, el fluido de trabaja entra al colector a travs de un tubo secundario y
all es dirigido hacia los tubos principales, tubos en los cuales tiene lugar al transferencia de calor. Una
vez recorrido el tubo principal, el fluido es recogido por otro tubo secundario, que se encarga de dirigirlo
hacia el resto del circuito.
A la vista del esquema, puede deducirse que van a existir prdidas de carga debidas a la friccin tanto en
los tubos principales como en los secundarios. Adems se producirn unas prdidas de presin
adicionales en los accesorios del circuito hidrulico (conexiones en T), con diferente distribucin de
canales y tipo de flujo (unin y bifurcacin).
La magnitud de las prdidas de presin que experimenta un flujo incompresible por un tramo recto de
tubera se calcula mediante la ecuacin de Darcy-Weisbach:

Donde v es la velocidad del fluido, en m/s, que puede expresarse como

y f
D
es el factor de friccin de Darcy, que depender del tipo de flujo: laminar o turbulento.
- Flujo laminar (0<Re<2300)

- Flujo turbulento (Re>2300)

Los coeficientes de friccin de Fanning y de Darcy se relacionan como sigue:

Por cada uno de los tubos principales tendremos la misma prdida de carga, pues los tubos estn
dispuestos en paralelo y se considera igual flujo msico circulando por cada uno de ellos.
25

Para hallar las prdidas de presin debidas a la friccin de los tubos secundarios se hara necesario
obtener el flujo msico en cada seccin pues sta va variando a lo largo del tubo. Para simplificar los
clculos, se estimarn dichas prdidas hallando un valor medio para ste.
Por otro lado, la magnitud de las prdidas que se originan en un accesorio se calcula mediante la
siguiente ecuacin:

Siendo K el coeficiente de prdidas.
Se elegir una lnea de corriente y se estimar la prdida de carga debida a los accesorios:

Donde
son las constantes de prdidas.

3. Mecanismos de transferencia de calor en captadores solares
3.1. Transferencia de calor con el fluido caloportador en rgimen laminar

La transmisin de calor en rgimen laminar es un fenmeno complejo en el que intervienen las
siguientes variables:
- La condicin de contorno trmica en la pared.
- Los efectos de entrada.
- La variacin de las propiedades fsicas del fluido con la temperatura.
- La conveccin natural por efecto de la flotacin.
La condicin de contorno trmica en la pared (flujo de calor q
p
= uniforme o temperatura de calor T
p
=
constante) influye en gran medida en la transferencia de calor en rgimen laminar, al contrario de lo que
ocurre en rgimen turbulento. Las correlaciones experimentales o soluciones numricas no se van a
poder extrapolar a condiciones de contorno distintas de aquellas para las que se obtuvieron.
La dificultad en el estudio del flujo laminar tiene que ver en gran medida con el hecho de que los fluidos
que transfieren calor en este rgimen suelen tener propiedades fsicas que dependen fuertemente de la
temperatura. Los efectos de estas variaciones en las propiedades fsicas, fundamentalmente viscosidad
y densidad, van a influir notablemente en la transmisin de calor. En flujo laminar, los gradientes de
temperatura existentes en el fluido provocan variaciones espaciales de densidad que propiciarn la
aparicin de fuerzas de flotacin. Cuando estas fuerzas de flotacin son dominantes, crean un
movimiento secundario y la transferencia de calor asociada se denomina conveccin natural o libre.
Si la transferencia de calor se debe por completo a un movimiento impuesto por gradientes externos de
presiones motrices, se tiene conveccin forzada; a la accin combinada de la conveccin libre y forzada
se le denomina conveccin mixta y est presente normalmente en las aplicaciones de transmisin de
calor en rgimen laminar.

26

3.1.1. Consideraciones hidrodinmicas
3.1.1.1. Condiciones de flujo
Considere un flujo laminar en un tubo circular de radio r
0
Figura 2.12 donde el fluido entra al tubo con
velocidad uniforme. Sabemos que cuando el fluido hace contacto con la superficie, los efectos viscosos
se vuelven importantes y su produce una capa lmite al aumentar . Este desarrollo ocurre a expensas
de una regin de flujo no viscosos que se contrae y concluye con la unin de la capa lmite en la lnea
central. Despus de esa unin, los efectos viscosos se extienden sobre toda la seccin transversal y el
perfil de velocidad ya no cambia al aumentar . Se dice entonces que el flujo est completamente
desarrollado. La longitud necesaria para que se de esta situacin, es la longitud de entrada
hidrodinmica
.

Figura 2.12: Desarrollo de la capa lmite hidrodinmica en un tubo
Para flujo laminar (
, la longitud hidrodinmica de entrada se puede obtener a partir de una

expresin de la forma (Langhaar, 1942):

Esta expresin se basa en la suposicin de que el fluido entra al tubo desde una boquilla redonda
convergente y por ello se caracteriza mediante un perfil de velocidad casa uniforme a la entrada.
3.1.1.2. Perfil de velocidad en la regin completamente desarrollada
La forma del perfil de velocidad se puede determinar sin dificultad tratndose del flujo laminar de un
fluido incompresible de propiedades constantes en la regin completamente desarrollada de un tubo
circular. Una caracterstica importante de las condiciones hidrodinmicas en la regin completamente
desarrollada es que la componente radial de la velocidad y el gradiente de la componente axial de la
velocidad
son cero en todas partes.

En consecuencia, la componente axial de la velocidad depende solo de .
Para obtener el perfil de velocidades podemos realizar un simple balance de cantidad de movimiento,
que finalmente se reduce a un balance entre las fuerzas cortantes y de presin en el flujo. Haciendo
esto y aplicando la ley de la viscosidad de Newton, se obtiene:

1
Nmero de Reynolds: es un nmero adimensional que relaciona las propiedades fsicas del fluido, su
velocidad y la geometra del conducto.
27

Figura 2. 13: Balance de fuerzas sobre un elemento diferencial para el flujo laminar completamente desarrollado
en un tubo circular
Las constantes de integracin se pueden determinar al recurrir a las condiciones de frontera:

que, respectivamente, imponen los requerimientos de deslizamiento cero en la superficie del tubo y de
simetra radial alrededor de la lnea central. Es sencillo evaluar las constantes, y se sigue que:

De aqu se observa que el perfil de velocidad completamente desarrollado es parablico. Ntese que el
gradiente de presin siempre deber ser negativo.
Si ahora se calcula la velocidad media en el conducto, como:

Si se introduce la ecuacin (2.50) en la (2.51), se obtiene:

La sustitucin de este resultado en la ecuacin (5), da el perfil de velocidad:

3.1.1.3. Gradiente de presin y factor de friccin
Para determinar la cada de presin y poder compararla se trabaja con un parmetro adimensional,
factor de friccin de Darcy, que se define como:

Este valor se comparar con la solucin analtica, para flujo laminar completamente desarrollado sobre
tubo liso.
28

3.1.2. Consideraciones Trmicas
Despus de revisar la mecnica de fluidos con respecto al flujo interno, consideramos ahora los efectos
trmicos. Si entra fluido al tubo de la figura 3 a una temperatura uniforme T(r,0) que es menor que la
temperatura de la superficie, ocurre la transferencia de calor por conveccin y se comienza a producir
una capa lmite trmica. Adems, si la condicin de la superficie del tubo se fija mediante la imposicin
de una temperatura uniforme (
o un flujo de calor uniforme (
, finalmente se
alcanza una condicin trmica completamente desarrollada. La forma del perfil de temperatura
completamente desarrollada difiere segn se mantenga una temperatura superficial uniforme o
un flujo de calor constante. Para ambas condiciones la temperatura de entrada aumenta al aumentar .
Para el flujo laminar la longitud de entrada trmica se puede expresar como (Kays, 1980):

Al comparar las ecuaciones (1) y (11) es evidente que, si , la capa lmite hidrodinmica se
desarrolla ms rpido que la capa lmite trmica
, mientras que lo inverso es cierto para

. Por ello, como en las simulaciones se ha utilizado como fluido de trabajo el agua, el nmero de
Prandtl para el rango de temperaturas de trabajo siempre es mayor que la unidad, luego la capa lmite
hidrodinmica es ms corta que la trmica como efectivamente ms adelante se comprobar.
3.1.2.1. Regin completamente desarrollada

- Flujo de calor uniforme en la pared
Cuando se aplica un flujo de calor uniforme a un tubo por el que fluye un fluido la temperatura media
del fluido,
, se incrementa linealmente con la distancia axial a lo largo del tubo. El valor de
puede calcularse haciendo un balance energtico desde la entrada a la seccin de intercambio

trmico hasta una posicin genrica :

donde
es la temperatura media del fluido a la entrada de la seccin de intercambio trmico (x=0). La

temperatura media del fluido en una posicin queda:

De la ecuacin anterior, para un caudal y un aporte de calor dados y asumiendo propiedades fsicas
constantes, la temperatura media del fluido
vara linealmente con la distancia axial . Para

mantener la condicin de contorno de flujo de calor uniforme, la temperatura de la pared
debe
tambin variar linealmente con . Se tiene as los perfiles de temperaturas que se esquematizan en la
Figura 2.14:

29

Figura 2.14: Evolucin del perfil de temperatura del fluido en el flujo desarrollado con condicin de contorno q
p
''
uniforme
Para un flujo totalmente desarrollado con
, el nmero de Nusselt en conveccin forzada

puede determinarse analticamente resolviendo la ecuacin de la energa en forma diferencial. Este
nmero de Nusselt asinttico en conveccin forzada toma un valor constante, quedando:

- Temperatura de pared constante en la pared
Si se tiene una condicin de contorno de temperatura de pared constante entonces el flujo de calor ya
no es uniforme a lo largo de . La temperatura media del fluido se aproximar asintticamente a la de
la pared y el flujo de calor en la pared ir decreciendo en consecuencia. La figura siguiente esquematiza
el proceso. Con un tubo infinitamente largo, la temperatura del fluido alcanzara la de la pared y el flujo
de calor cesara.

Figura 2.15: Evolucin del perfil de temperatura del fluido en el flujo plenamente desarrollado con condicin de
contorno T
p
constante
Para un flujo totalmente desarrollado con
, el nmero de Nusselt en conveccin forzada puede

determinarse analticamente resolviendo la ecuacin de la energa en forma diferencial. Este nmero de
Nusselt asinttico en conveccin forzada toma un valor constante, quedando:

30

3.1.2.2. Regin de entrada trmica
La solucin a la ecuacin de energa, para la regin de entrada en ms difcil de obtener, pues la
velocidad y temperatura dependen ahora de as como tambin de . Aun si se ignora el trmino de
adveccin radial, el gradiente de temperatura axial tal vez ya no se puede simplificar mediante la
ecuacin . La solucin ms sencilla es suponer que las condiciones trmicas se generan en
presencia de un perfil de velocidad completamente desarrollado, es lo que se llama el problema de
longitud de entrada trmica. Se trata de una suposicin razonable para fluidos con Prandtl grande.
Desde el momento en que comienza el proceso de calentamiento, el efecto de transferencia de calor en
la pared penetra gradualmente en el flujo. Primero se aprecia una delgada capa lmite trmica que
crece en direccin aguas abajo y finalmente se vuelve tan gruesa como el propio conducto. El fluido
adopta un perfil de temperaturas que vara axial y transversalmente
con lo que la
temperatura media del fluido va a depender de la posicin axial

En un proceso de calentamiento, el salto de temperaturas pared-fluido va a ser mximo en el punto de
comienzo de la seccin de calentamiento e ir disminuyendo a medida que el fluido avanza axialmente.
De este modo se tiene que en la regin de entrada trmica el coeficiente de transferencia de calor es
local y depende de la posicin axial, tomando un valor mximo al comienzo en el punto de partida y
disminuyendo hasta alcanzar un valor constante en donde el flujo se ha desarrollado plenamente. La
longitud de entrada trmica,
, es la distancia axial necesaria para conseguir un valor de nmero de

Nusselt 1,05 veces el valor correspondiente a flujo plenamente desarrollado. De este modo en una
posicin axial x se definen dos valores de nmero de Nusselt: el nmero de Nusselt local,
, y el
nmero de Nusselt global,
. Para calcular el nmero de Nusselt global en un tubo de longitud l se

debe de promediar el nmero de Nusselt local a lo largo de toda su longitud:

Dentro de la regin de entrada los nmeros de Nusselt local y global deben darse referidos a una
longitud adimensional. La distancia adimensional del punto de medida al inicio de la seccin de ensayo
trmica es
2
:

Tambin es habitual emplear el nmero de Graetz, Gz:

3.1.2.3. Correlaciones de nmero de Nusselt en conveccin forzada

2
La distancia adimensional, x
*
, se puede considerar como la relacin entre el tiempo que un elemento
fluido ha permanecido en el tubo y el tiempo necesario para que el calor se difunda en l.
31

Para un flujo desarrollado hidrodinmicamente la solucin numrica, tanto con condicin de contorno
de flujo de calor uniforme como de temperatura de calor constante, la di por primera vez (Shah, 1978).
En la Figura 2.16 se puede ver los resultados para el nmero de Nusselt local y medio:

Figura 2.16: Conveccin forzada: solucin numrica del flujo para q
p
'' uniforme y T
p
constante. Nmeros de
Nusselt local y medio. De A. Garca (2006).
En este caso se van a utilizar correlaciones empricas que aproximan la solucin numrica. Destacan la
correlacin de (Petukhov, 1988) y las de (Churchill, 1973).
Correlacin de Churchill y Ozoe [1973].
Propusieron la siguiente correlacin para el nmero de Nusselt local para flujos laminares en desarrollo
simultneo con flujo de calor uniforme:

Un flujo en desarrollo simultaneo de las regiones trmica e hidrulica el nmero de Nusselt depende de
la distancia adimensional y de Prandtl. En los casos con altos nmero de Prandtl, se puede despreciar la
longitud de entrada hidrulica por lo que la correlacin anterior se aproxima mucho a la solucin
numrica para flujo hidrodinmicamente desarrollado.
Correlacin de Petukhov y Polyakov [1988]
Estos autores propusieron la siguiente correlacin para el nmero de Nusselt local en conveccin
forzada para flujos laminares desarrollados hidrodinmicamente con flujo de calor uniforme.

32

Esta correlacin es vlida para un flujo desarrollado hidrodinmicamente o para fluidos de alto nmero
de Prandtl. En la siguiente figura se muestra la correlacin comparada con la solucin numrica del
problema:

Figura 2.17: Conveccin forzada con flujo de calor constante. Nmero de Petukhow y Polyakov [1988] frente a la
solucin numrica. D A. Garca (2006).

3.1.2.4. Variacin de las propiedades fsicas del fluido

Cuando la diferencia de temperaturas entre la pared del tubo y el fluido es sustancial se harn patentes
las variaciones radiales en las propiedades fsicas del fluido, fundamentalmente en la capa lmite. En
lquidos la variacin que hay que tener en cuenta es la de la viscosidad dinmica, puesto que el resto de
propiedades fsicas (conductividad trmica, calor especfico y densidad) varan poco con la temperatura.
Por ello en el modelo utilizado las propiedades fsicas del fluido se han puesto en funcin de la
temperatura a travs de polinomios como se puede ver en el Anexo I.
3.1.3. Los efectos de flotacin: Conveccin Mixta

En tubos horizontales, la direccin de la corriente de conveccin forzada es transversal al vector de la
aceleracin gravitatoria. En procesos de calentamiento, el fluido en contacto con la pared est ms
caliente (ms ligero) y asciende a lo largo de sta. Se forman dos corrientes que convergen en la
generatriz superior del tubo, donde por continuidad desalojan al fluido que all se encuentra y lo obligan
a descender hacia la parte central de la seccin de paso del fluido. Se producen as dos recirculaciones
simtricas en direccin transversal al flujo principal; la superposicin de ste flujo secundario con el flujo
principal se traduce en una trayectoria en espiral del fluido.
Las corrientes que se esquematizan en la figura siguiente producen un aumento significativo del
coeficiente de transmisin de calor en la parte inferior del tubo y una disminucin en la parte superior.
El coeficiente de transmisin de calor local
es ahora variable en la seccin,
por lo
que se debe hablar de un coeficiente de transmisin de calor local promediado en la seccin
:
33

Figura 2.18: Perfiles de flujo en conveccin mixta laminar en un tubo horizontal. Arriba: recirculaciones simtricas
en una seccin transversal; abajo: trayectoria en espiral del fluido a lo largo del tubo.

Con una condicin de contorno de flujo de calor uniforme, la temperatura en la pared ser variable en la
seccin y se deber realizar un buen promediado para calcular el nmero de Nusselt local:

siendo
la temperatura media de la pared en la seccin transversal.

Al igual que en conveccin forzada, en conveccin mixta la condicin de contorno trmica influye en
gran medida en la transmisin de calor. Con una condicin de flujo de calor
, una vez
desarrollado trmicamente el flujo se mantiene constante la diferencia entre la temperatura de la pared
y la temperatura del fluido. Esta caracterstica va a hacer que los efectos de conveccin permanezcan
uniformes a lo largo de todo el tubo. Por el contrario, con una condicin de contorno de temperatura
de pared constante, el flujo secundario crece primero en intensidad en la regin de entrada trmica
para ir desapareciendo despus a medida que la temperatura del fluido se aproxima a la de la pared.
En procesos de calentamiento, la presencia de las recirculaciones tpicas de conveccin natural beneficia
notablemente la transmisin de calor, con coeficientes de transferencia que superan en varias veces los
valores en conveccin forzada.
3.1.3.1. Parmetros adimensionales

Los parmetros que se deben de tener en cuenta para la conveccin mixta son:
- En las ecuaciones de cantidad de movimiento, Grashof y Reynolds:
El nmero de Grashof, Gr, representa un balance entre las fuerzas de flotacin debidas a una diferencia
de temperaturas caracterstica entre la pared y el fluido y las fuerzas viscosas presentes.
Alternativamente al nmero de Grashof es habitual encontrar el nmero de Rayleigh, Ra = GrPr. La
diferencia de temperatura DT est relacionada con la condicin de contorno trmica y da lugar a
diferentes definiciones de los nmeros de Grashof y Rayleigh:
34

El cociente
que aparece en las ecuaciones de cantidad de movimiento indica la importancia

relativa de los efectos de flotacin frente a los efectos de conveccin forzada. Si el valor del parmetro

es mucho menor que la unidad, entonces los ltimos trminos de las ecuaciones de cantidad de
movimiento son despreciables y se tiene un problema de conveccin forzada pura. Por el contrario, para
valores de
mucho mayores que la unidad, los trminos de flotacin son dominantes y se tiene
un problema de conveccin natural pura o conveccin libre. Para valores de
entre los dos

lmites anteriores, el problema es de conveccin mixta, donde los efectos de conveccin natural y
forzada tienen una magnitud semejante.
- En la ecuacin de la energa, Reynolds, Prandtl y Eckert
El nmero de Eckert slo es importante en fluidos con alta velocidad y es por ello despreciable en flujos
laminares, donde las velocidades son normalmente bajas. Por tanto el trmino viscoso de disipacin de
energa se puede despreciar normalmente.
Si se tiene en cuenta que el coeficiente de transferencia de calor en conveccin est directamente
relacionado con el campo de temperaturas:

y que el nmero de Nusselt no es ms que el coeficiente de adimensionalizado,

entonces se puede concluir que el nmero de Nusselt local promediado en la pared depende de los
parmetros que resultan de la adimensionalizacin de las ecuaciones del movimiento:

3.1.3.2. Correlaciones de conveccin mixta para flujo de calor uniforme
Correlacin de Petukhov y Polyakov [1988]
Puede decirse que los estudios de estos autores son los ms completos hasta la fecha sobre los efectos
de la conveccin natural en la transmisin de calor en conductos. En el libro de Petukhov y Polyakov
[1988] sobre conveccin mixta en rgimen turbulento incluyeron un captulo sobre conveccin mixta
laminar, basndose en sus propios datos experimentales obtenidos con agua aos atrs. Propusieron un
nmero de Grashof lmite a partir del cual se hacen apreciables los efectos de la conveccin natural:

Para el nmero de Nusselt local en conveccin mixta laminar promediado en la seccin,
,
obtuvieron la siguiente correlacin:

35

donde
, es el nmero de Nusselt en conveccin forzada pura, ecuacin (20),
, es el nmero de
Rayleigh para una condicin de contorno de flujo de calor uniforme y , es una funcin de x*:

En la siguiente figura se muestra la variacin del nmero de Nusselt con la distancia adimensional para
distinto nmeros de Rayleigh
:

Figura 2.19: Nmero de Nusselt de Petukhov y Polyakov [1988] en convenccin mixta
Segn la figura anterior, En conveccin forzada pura la regin trmica plenamente desarrollada empieza
en
; aqu el nmero de Nusselt se hace constante e igual al conocido valor de
. Sin
embargo, en conveccin mixta la regin en desarrollo trmico es mucho menor (
).
Estos investigadores tambin dieron una correlacin para flujo desarrollado, donde el Nusselt se hace
independiente de x*:

De la ecuacin (2.66), se puede deducir que en la regin de entrad el nmero de Nusselt local es
proporcional a
. Por otro lado, la ecuacin (2.77) indica que en la regn plenamente

desarrollada se tiene un nmero de Nusselt asinttico que nicamente depende del nmero de
Rayleigh,
.

36

Correlacin de Hong et al. [1974]
Los autores realizaron un estudio experimental y analtico en tubos horizontales, desarrollando
nicamente correlaciones para flujo plenamente desarrollado. Propusieron la siguiente correlacin para
sus resultados numricos.

Las experiencias las llevaron a cabo en un tubo de vidrio y otro de acero inoxidable, con agua y
etilenglicol como fluidos de ensayo. Los nmeros de Nusselt obtenidos con el tubo metlico fueron
superiores a los obtenidos con el tubo de vidrio. Para tener en cuenta el efecto de la conductividad
trmica de la pared, incluyeron en sus correlaciones un parmetro adimensional Pw:

que representa el cociente entre las resistencias trmicas circunferencial y radial en lapared del tubo.
Propusieron la siguiente correlacin5 de sus datos experimentales denmero de Nusselt plenamente
desarrollado:

donde todas las propiedades fsicas estn basadas en la temperatura de la capa lmite,

Correlacin de Morcos y Bergles [1975]
El trabajo de estos autores es una ampliacin experimental del trabajo desarrollado por Hong et al
[1974]. Realizaron nuevos experimentos en dos tubos de distinto material: uno de acero inoxidable y
otro de vidrio, obteniendo mayores nmeros de Nusselt en el tubo metlico. nicamente realizaron
experiencias en flujo plenamente desarrollado, proponiendo la siguiente correlacin experimental:

donde las propiedades fsicas del fluido s=e evalan a la temperatura de la capa lmite.La Ecuacin (2.81)
es vlida en los rangos:

37

3.2. Conveccin Natural en Cavidades

3.2.1. Introduccin
La tasa de transferencia de calor entre dos placas inclinadas dentro de cierto ngulo tiene una gran
importancia a la hora de mejorar la eficiencia de un colector solar. El efecto de la gravedad junto con la
diferencia de temperaturas, hacen posible la aparicin de flujos convectivos debido al gradiente de
densidades del aire encerrado entre el absorbedor y la cubierta o cubiertas de vidrio. Por ello se va a
analizar las ecuaciones bsicas que gobiernan este modo de transferencia de calor, se har una revisin
de los estudios que abarcan este mbito prestando una especial atencin a la conveccin natural en el
interior de recintos cerrados rectangulares y a aquellos que se han hecho sobre colectores solares en
parrilla como el que aqu se est modelando. Finalmente se presentar el modelo que se ha utilizado y
cules han sido los resultados obtenidos.
3.2.2. Ecuaciones gobernantes
Las ecuaciones que describen la transferencia de momento y energa en la conveccin libre se originan
de los principios de conservacin relacionados. Las fuerzas inerciales y viscosas siguen siendo
importantes, pero a diferencia de la conveccin forzada las fuerzas de empuje juegan un papel principal.
Considera un flujo laminar de capa lmite gobernado por fuerzas de flotacin. Suponga condiciones de
propiedades constantes bidimensionales estables en las que la fuerza de gravedad acta en la direccin
x negativa. Tambin, con una excepcin, suponga que el fluido es incompresible. La excepcin implica
explicar el efecto de la densidad variable en la fuerza de flotacin (aproximacin de Boussinesq), pues
esta variacin es la que induce el movimiento del fluido. Finalmente, suponga que son vlidas las
aproximaciones de capa lmite.
Con estas simplificaciones, la ecuacin de momento en x es:

La ecuacin anterior se puede expresar en una forma ms conveniente advirtiendo primero que, si en la
direccin y no hay fuerza de gravedad, entonces
. De aqu el gradiente de presin en la direccin

x, en cualquier punto en la capa lmite debe ser igual al gradiente de presin en la regin en reposo
fuera de la capa lmite. Sin embargo, en esta regin u=0, y la ecuacin anterior se reduce a:

la cual se debe aplicar a cualquier punto en la capa lmite de conveccin libre.
El primer trmino en el lado derecho de la ecuacin es la fuerza de flotacin, y el flujo se origina debido
a que la densidad es variable. El origen de esta variacin se puede hacer ms explicito al introducir el
coeficiente volumtrico de expansin trmica:

38

Figura 2.20: Desarrollo de la capa lmite sobre una placa vertical caliente
Esta propiedad termodinmica de fluido proporciona una medida de la cantidad por la cual cambia la
densidad en respuesta a un cambio en la temperatura a presin constante. Si se expresa en la siguiente
forma aproximada:

se sigue que:

Al sustituir en la ecuacin de momento en la direccin x, se tiene:

donde ahora es aparente cmo la fuerza de flotacin, que impulsa al flujo, se relaciona con la diferencia
de temperaturas.
Como los efectos de empuje se confinan a la ecuacin de momento, las ecuaciones de conservacin de
masa y energa permanecen sin cambio de la conveccin forzada. El conjunto de ecuaciones que
describen el problema son:

39

Para el caso del aire en el interior del absorbedor se puede considerar como un gas ideal, y en ese caso
el coeficiente de expansin es:

donde T es la temperatura absoluta.

3.2.3. Parmetros adimensionales
Consideremos ahora los parmetros adimensionales que gobiernan el flujo de conveccin libre y la
transferencia de calor. Introduciendo a las ecuaciones que gobiernan el flujo, los siguientes parmetros:

donde L es una longitud caracterstica y
es una velocidad de referencia arbitraria, las ecuaciones de

momento en x y de energa se reducen a:

El parmetro adimensional en el primer trmino del lado derecho de la ecuacin de cantidad de
movimiento es una consecuencia directa de la fuerza de empuje. Sin embargo, como se expresa en
trminos de la velocidad de referencia
, no es conveniente en su forma acutal. Se acostumbra por

tanto a trabajar con una forma alternativa que se obtiene al multiplicar por
. El resultado es el
nmero de Grashof Gr
L
.

El nmero de Grashof indica la razn de las fuerzas de empuje a las fuerzas viscosas que actan sobre el
fluido.

40

3.2.4. Conveccin natural en el interior de recintos
Los problemas en recintos surgen cuando una superficie solida envuelve completamente un fluido,
como es el caso del aire contenido entre el absorbedor, la cubierta de vidrio, y el aislamiento lateral. La
geometra del problema se puede representar esquemticamente segn la figura siguiente. La cavidad
est limitada por dos placas paralelas e isotermas (en el caso del captador solar esta condicin no se
cumple, pero podemos tomar estas investigaciones como punto de partida) inclinadas un ngulo
desde la horizontal, espaciadas una distancia L. La temperatura T
h
es mayor que la temperatura T
c
, luego
aquella cavidad con =0, es una cavidad horizontal calentada desde abajo, cuando =90, entonces es
una cavidad vertical calentada desde un lado, y cuando =180, se trata tambin de una cavidad
horizontal calentada desde arriba.

Figura 2.21: Geometra del recinto

El rea de cada placa que est en contacto con el fluido es
y la pared que limita al fluido por

los laterales tiene un espesor uniforme b y una conductividad trmica
. La condicin de contorno
trmica en el exterior de las superficies es adiabtica. El fluido, a su vez, se pude considerar con
completamente transparente o completamente opaco a la radiacin trmica. Normalmente los lquidos
se consideran opacos y los gases transparentes. Si el fluido es transparente (como ocurre en el caso de
los colectores) el intercambio de radiacin puede afectar a la conveccin natural alterando la
distribucin de temperaturas.
Se pueden asignar dos nmeros de Nusselt a los problemas con cavidades. El primero de ellos:

y el segundo:

donde
, es el flujo de calor que no proviene de la radiacin de la superficie interior que se encuentra a

una temperatura T
h
,
, es el flujo de calor radiante que proviene tambin de esa cara y
es el flujo de
calor de conduccin desde el rea A
w
de la cara caliente . La transferencia de calor
, son los
respectivos valores de
cuando el fluido est en estado estacionario y su comportamiento

trmico es como el de un slido. Para fluidos que son opacos:
41

Si el fluido es transparente,
se deben determinar por un anlisis combinado de

radiacin y conduccin (Hollands, 1984)
El anlisis dimensional revela que, en el caso ms general:

El nmero de Rayleigh basado en la dimensin L:

3.2.5. Nmero de Rayleigh crticos para cavidades horizontales
El nmero de Rayleigh crtico gobierna el inicio del movimiento por efectos convectivos del fluido en el
interior de la cavidad. En cavidades que se encuentran restringidas en direccin horizontal (es decir,
W>>L y H>>L) el Ra
c
depende de los parmetros geomtricos que describen la cavidad y de las
propiedades trmicas de la pared, pero no del nmero de Prandtl. Para =0, Ra
c
, se encuentra limitado
por los valores de Ra
cp
y Ra
ci
. Los trabajos de (Catton, 1972) evalan los valores anteriores y los
tabularon para un paraleleppedo rectangular. El mayor de estos valores es Ra
cp,
que se aplica a paredes
perfectamente conductoras (no es nuestro caso) y el menor de ellos Ra
ci
que se aplica a paredes a
paredes adiabticas (es el que ms se acerca al modelo de captador). Segn estos estudios para una
cavidad rectangular:

Estos estudios tabularon los valores crticos del Ra (Figura 2.22). Si nos fijamos en las dimensiones para
el modelo de 1 solo tubo, 1.83x0.123x0.02, vemos que nos encontramos prcticamente en la zona
asinttica (W/L=6.15 , H/L=91) y Ra
c
=1708.
3.2.6. Conveccin en cavidades rectangulares
Cavidades extensas en direccin horizontal (H>>L y W>>L). Cuando =0, el fluido caliente, ms ligero,
se encuentra debajo del fluido frio ms pesado, lo que da lugar a una inestabilidad inherente. A pesar de
esta inestabilidad el fluido permanece estacionario si el nmero de Rayleigh (Ra) es menor que un valor
crtico, para este tipo de geometra dicho valor es 1708. Para los casos donde Ra>Ra
c
la inestabilidad
conduce a un movimiento convectivo que depende de Ra y Pr. Para un Ra ligeramente superior a Ra
c
, el
movimiento consiste en cilindros estables del orden de L, pero conforme aumenta Ra, se dan patrones
de flujo ms complejos, y eventualmente se convierte en inestable. Para valores muy altos de Ra, se
vuelve turbulento. Las caractersticas de la transferencia de calor reflejan la existencia de diferentes
regmenes de flujo: para Ra< Ra
c
el fluido esta estacionario, y por lo tanto Nu=1, el movimiento en
clulas se inicia en Ra
c
produciendo un fuerte aumento de Nussel, que finalmente se vuelve asinttico
alrededor del valor de
para valores altos de Ra.

42

Figura 2.22: Nmero de Rayleig Crticos para una cavidad horizontal
Para =0, la ecuacin recomendada para Nu (Holland, 1984):

donde los asteriscos en los corchetes indican que solamente se deben tomar los valores positivos, es
decir:

Los valores de
son funcin de Pr, y se encuentran tabulados en la tabla siguiente:

Las siguientes ecuaciones fijan la dependencia de
con Pr:

43

En la siguiente figura se puede ver la correlacin 2.106, dibujada para un rango de Prandlt propio del
agua a temperaturas moderadas comparada con otros trabajos experimentales.

Figura 2.23: Comparacin de la ecuacin I con los resultados de Garon y Goldstein , Chu y Goldstein y Rossby para
agua (Pr=6)
3.2.7. Transferencia de calor en cavidades rectangulares inclinadas
Dependiendo de la inclinacin , el flujo en el interior de una cavidad con W/H>8 puede parecerse al
correspondiente en una cavidad horizontal o vertical, es raro que sea una combinacin de las
caractersticas de ambos. Consecuentemente, salvo contadas excepciones, el nmero de Nusselt en una
cavidad inclinada se puede determinar, como una aproximacin razonable, ya sea de la solucin del
Nusselt horizontal o vertical, por medio de leyes escala angular.
Para el caso de cavidades con
, que es el que ms se acerca a las condiciones que

estamos simulando. La correlacin ms recomendada para ngulos comprendidos entre 0 y 60 grados y
para nmeros de Rayleigh por encima de 10
5
es la de Hollands et al:

Sin embargo para valores de Ra=10
6
, la correlacin anterior subestima los resultados obtenidos por
(S.M. Elsherbiny, 1982) en un 10%.

44

3.3. Transmisin de calor por radiacin

3.3.1. Transferencia de calor por radiacin
La transferencia de calor por radiacin tiene el objetivo de llevar hasta el captador la energa necesaria
para calentar el fluido calo-portador. Aunque en la mayora de los problemas de la ingeniera la
transferencia de calor por radiacin es despreciable, para el caso de colectores solares donde la flujo de
energa es normalmente dos rdenes de magnitud ms pequeo que en el resto de equipamiento
trmico convencional , la radiacin trmica es un modo de transferencia de calor importante.
3.3.2. El espectro electromagntico

La radiacin trmica es energa electromagntica que se propaga por el espacio a la velocidad de la luz.
Para la mayora de las aplicaciones de energa solar, nicamente la radiacin trmica es importante. La
radiacin trmica es emitida por cuerpos en virtud de su temperatura, los tomos, molculas o
electrones por los que estn formados, se elevan a estados excitados volviendo espontneamente hacia
estados de menor energa gracias a la emisin de energa en forma de radiacin electromagntica.
El espectro de radiacin electromagntica se divida en bandas de diferente longitud de onda. Estas
bandas y la longitud de onda asociada a cada una se representan en la siguiente figura:

Figura 2.24: Espectro electromagntico
El ancho de longitudes de onda de importancia para la radiacin solar y sus aplicaciones van desde el
ultravioleta hasta el infrarrojo, esto es desde 0.3 hasta aproximadamente 25 m. La radiacin solar
fuera de la atmsfera tiene la mayora de su energa en el rango de 0.25 a 3 m.
3.3.3. El Cuerpo Negro

Un cuerpo negro es un perfecto absorbedor de radiacin. No importa que longitud de onda o qu
direccin describa la radiacin incidente sobre l , toda la que incide sobre el mismo ser absorbida. Por
lo que un cuerpo negro es un concepto ideal ya que todas las sustancias van a reflejar una parte de la
radiacin incidente.
Adems, un cuerpo negro es un perfecto emisor de radiacin. Se puede utilizar un experimento simple
para demostrar que si un cuerpo es un perfecto emisor de radiacin, entonces tambin debe ser un
perfecto absorbedor. Supongamos que tenemos un pequeo cuerpo negro y otro no-negro dentro de
un gran recinto al que se ha hecho el vacio. Si este espacio est aislado entonces tanto el recinto, como
45

el cuerpo negro como el cuerpo real alcanzarn al mismo tiempo un temperatura de equilibrio. El
cuerpo negro por definicin absorber toda la radiacin incidente y para mantener una temperatura
constante tambin deber emitir una cantidad igual de energa. El cuerpo real absorber menos
radiacin que el cuerpo negro y consecuentemente emitir menos radiacin que el cuerpo negro.

3.3.4. La ley de Planck y la ley de desplazamiento de Wien

La radiacin en la regin del espectro electromagntico, es decir, desde 0.2 m hasta aproximadamente
1000 m se llama radiacin trmica, y es emitida por sustancias en virtud de su temperatura. La
distribucin de longitud de onda emitida por un cuerpo negro est dada por la ley de Planck:

Donde h es la constante de Planck y k es la constante de Boltzmann. El conjunto

normalmente se llaman primera y segunda constantes de radiacin y se representan por los smbolos
.
Derivando e igualando a cero la ley de Planck podemos conocer la longitud de onda para la cual
tenemos la mxima intensidad de radiacin de cuerpo negro. Esto nos lleva a la ley de desplazamiento
de Wien, que se puede escribir como:

En la figura siguiente se representa la ley de Planck que muestra la distribucin de radiacin para un
cuerpo negro a temperaturas de 6000, 1000 y 400 K. La temperatura de 6000 K, es la de la superficie
solar, por lo que representa la distribucin de radiacin solar fuera de la atmsfera terrestre.

46

Figura 2.25: Distribucin de energa emitida por un cuerpo negro
3.3.5. Ecuacin de Stefan-Boltzmann

La ley de Planck representa la distribucin del espectro radiacin en un cuerpo negro, pero
normalmente en aplicaciones prcticas se requiere el clculo de la energa total. Por lo tanto si
integramos la ley de Planck en todo el rango de longitudes de onda encontramos que la energa emitida
por un cuerpo negro debe ser:

donde es la constante de Stefan-Boltzmann que es igual a 5.667 x10
-8

.
3.3.6. Intensidad de radiacin

La radiacin que emite una superficie se propaga en todas las direcciones posibles, y a menudo estamos
interesados en conocer su distribucin direccional. Tambin, la radiacin incidente sobre una superficie
puede venir de diferentes direcciones, y la forma en que la superficie responde a esta radiacin
depende de la direccin. Tales efectos direccionales se pueden tratar mediante la introduccin del
concepto de intensidad de radiacin:

47

3.3.6.1. Definiciones

Considere la emisin en una direccin particular desde un elemento de rea
, como se muestra en la
Figura 2.26. Esta direccin se puede especificar en trminos de los ngulos cenital y acimutal, y ,
respectivamente, de un sistema de coordenadas esfricas (Figura 2.26). Una superficie diferencialmente
pequea en el espacio , a travs de la cual pasa esta radiacin, subtiende un ngulo slido cuando
se ve desde un punto sobre
. En la Figura 2.27 vemos que el ngulo plano diferencial se define

por una regin entre los rayos de un crculo y se mide como la razn del elemento de longitud de arco
sobre el crculo al radio r del crculo. De manera similar, en la Figura 2.27b el ngulo slido diferencial
est definido por una regin entre los rayos de una esfera y se mide como la razn del elemento de
rea
sobre la esfera al cuadrado del radio de la esfera. En consecuencia,

El rea
es normal a la direccin (, ), y segn se muestra en la Figura 2.28, se puede representar

como
para una superficie esfrica. Por consiguiente,

Figura 2.26: Naturaleza direccional de la radiacin. (a) Emisin de radiacin desde un rea diferencial dA1 en un
ngulo slido dw subtendido por dAn en un punto sobre dA1. (b) Sistema de coordenadas esfricas

48

Figura 2.27: Definicin de ngulos (a) plano y (b) slido
En tanto que el ngulo plano tiene unidades de radianes, la unidad del ngulo slido es el
esterorradin.
Si consideramos ahora la proporcin en que la emisin de
pasa a travs de
. Esta cantidad se
puede expresar en trminos de la intensidad espectral
de la radiacin emitida. Definimos

formalmente
, como la razn a la que se emite energa radiante a la longitud de onda en la

direccin (, ), por unidad de rea de la superficie emisora normal a esta direccin, por unidad de
ngulo slido alrededor de esta direccin, y por intervalo de longitud de onda unitaria alrededor de .
Ntese que el rea que se utiliza para definir la intensidad es la componente de
perpendicular a la
direccin de la radiacin. La Figura 2.29 muestra que esta rea proyectada es igual a
. La
intensidad espectral, que tiene unidades de
, es entonces,

donde
, es la razn a la que la radiacin de longitud de onda sale de
y pasa a
travs de
. Al reacomodar la ecuacin anterior, se sique que:

49

Figura 2.28: ngulo slido subtendido por dAn en un punto sobre dA1 en el sistema de coordenadas esfricas

Figura 2.29: Proyeccin de dA1 normal a la direccin de la radiacin
donde
tiene unidades de . Esta importante expresin nos permite calcualr la rapidez a al que
la radiacin emitida por una superficie se propgaga en la regin de espacio definida por el ngulo slido
alrededor de la direccin (, ). Al expresar la ecuacin anterior por unidad de rea de la supefcie
emisora y sustituir la ecuacin: , el flujo de radiacin espectral asociado con
es:

Si se conocen las distribuciones espectral y direccional de
, el flujo de calor asociado con la emisin

en cualquier ngulo slido finito o sobre cualquier intervalo finitio de longitudes de onda se puede
determinar al integrar la ecuacin anterior. Por ejemplo, el flujo de calor espectral asociado con la
emisin en un hemisferio hipottico arriba de
, como se muestra en la Figura 2.30, es:

50

Figura 2.30: Emisin de un elemento diferencial de rea dA1 en un hemisferio hipottico centrado en un punto
sobre dA1
3.3.6.2. Relacin con la emisin

La intensidad de radiacin est relacionada con varios flujos de radiacin importantes. La emisin ocurre
desde cualquier superficie que est a una temperatura finita. Se introduce el concepto de potencia
emisiva para cuantificar la cantidad de radiacin emitida por rea superficial unitaria. La potencia
espectral emisiva hemisfrica,
se define como la intensidad a la que se emite radiacin

de longitud de onda en todas direcciones desde una superficie por unidad de longitud de onda
alrededor de y por rea superficial unitaria. En consecuencia, se relaciona con la intensidad espectral
de la radiacin emitida mediante una expresin de la forma:

en cuyo caso es equivalente al flujo de calor dado por la ecuacin anterior. Ntese que
es un flujo
que se basa en el rea superficial real, mientras que
se basa en el rea proyectada. El trmino

que aparece en el integrando es una consecuencia de esa diferencia.
La potencia emisiva total hemisfrica,
, es la rapidez a la que se emite radiacin por unidad de

rea en todas las longitudes de onda y en todas las direcciones posibles. En consecuencia,

51

Aunque la distribucin direccional de la emsin superficial vara de acuerdo con la naturaleza de la
superficie, hay un caso especial que proporciona una aproximacin razonable para muchas superficies.
Hablamos de un emisor difuso como una superficie para la que la intensidad de la radiacin emitida es
independiente de la direccin en cuyo caso
. Entonces se tiene que:

3.3.6.3. Relacin con la irradiacin

Aunque nuestra atencin est puesta en la radiacin emitida por una superficie, los conceptos
anteriores se pueden extender a la radiacin incidente. Tal radiacin se puede originar de la emisin y
reflexin que ocurre en otras superficies y tendr distribucin espectrales y direccionales determinadas
por la intensidad espectral
. Esta cantidad se define como la proporcin en que la energa

radiante de longitud de onda incide de la direccin (, ) , por unidad de rea de la superficie
interceptora normal a esta direccin, por unidad de ngulo slido alrededor de esta direccin, y por
intervalo de longitud de onda unitaria alrededor de .
La intensidad de la radiacin incidente se puede relacionar con un flujo radiativo importante,
denominado irradiacin, que abarca la radiacin incidente desde todas direcciones. La irradiacin
espectral
se define como la velocidad a la que la radiaicn de longitud de onda incide

sobre una superficie, por unidad de rea de la superficie y por intervalo de longitud de onda unitaria
alrededor de . En consecuencia:

donde es el ngulo slido unitario. El factor se origina debido a que
es un flujo que
se basa en el rea superficial real, miestras que
se define en trminos del rea proyectada. Si la

irradiacin total
representa la razn a la que incide la radiacin por unidad de rea desde

todas las direcciones y a todas las longitudes de onda, se sigue que:

3.3.6.4. Relacin con la radiosidad

El ltimo flujo radiativo de inters, denominado radiosidad, explica toda la energa radiante que sale de
una superficie. Como esta radiacin incluye la parte reflejada de la irradiacin, as como la emisin
directa, la radiosidad es por lo general diferente de la potencia emisiva. La radiosidad espectral
representa la rapidez a la que la radiacin de longitud de onda unitaria alrededor .

Como explica la radiacin que sale en todas direcciones, se relaciona con la intensidad asociada con la
emisin y la reflexin,
, por una expresin de la forma:

De aqu la radiosidad total
asociada con todo el espectro es:

52

3.3.7. Emisin superficial

Despus de desarrollar la nocin de un cuerpo negro para describir el comportamiento de una superficie
ideal, podemos ahora considerar el comportamiento de superficies reales. Recuerde que el cuerpo
negro es un emisor ideal en el sentido de que ninguna superficie puede emitir ms radiacin que un
cuerpo negro a la misma temperatura. Es por tanto conveniente elegir el cuerpo negro como una
referencia al describir la emisin desde una superficie real. Una propiedad radiativa superficial conocida
como, emisividad, se puede definir como la razn de la radiacin emitida por la superficie a la radiacin
emitida por un cuerpo negro a la misma temperatura.
Es importante reconocer que, en general, la radiacin espectral emitida por una superficie real difiere
de la distribucin de Planck (Figura 2.31a). Adems, la distribucin direccional (Figura 2.31b) puede ser
diferente de la difusa. Por tanto, la emisividad puede tomar valores diferentes segn se est interesado
en la emisin a una longitud de onda dada o en una direccin dada, o bien en promedios integrados
sobre longitud de onda y direccin.

Figura 2.31: Comparacin de la emisin de un cuerpo negro y de una superficie real. Distribucin espectral (a).
Distribucin direccional (b)

Definimos la emisividad direccional espectral
de una superficie a una temperatura T

como la razn de la intensidad de la radiacin emitida en la longitud de onda y en la direccin de y
a la intensidad de radiacin emitida por un cuerpo negro a los mismos valores de T y . De aqu:

Observe cmo los subndices y designan el inters en una longitud de onda y direccin espefcficas
para la emisividad. Por el contrario, los trminos que aparecen dentro de los parntesis designan la
dependencia funcional respecto a la longitud de onda, direccin, y/o temperatura. La ausencia de
variables direccionales en el parntesis del denominador en la ecuacin anterior implica que la
intensidad es independiente de la direccin, que es, por supuesto, una caracterstica de la emisin del
cuerpo negro. De manera similar una emisividad direccional total
, que representa un promedio

espectral de
, se puede definir como:

53

En la mayora de los clculos de ingeniera, se desea trabajar con propiedades superficiales que
representen promedios direccionales. Una emisividad espectral hemisfrica por tanto se define como:

La emisividad total hemisfrica, que representa un promedio sobre todas las direcciones y longitudes de
onda posibles, se define como:

Por lo tanto tambin se puede definir como:

3.3.8. Absortividad

La absortividad es una propiedad que determina la fraccin de la irradiacin absorbida por una
superficie. La determinacin de la propiedad es complicada por el hecho de que, como la emisin, se
puede caracterizar por una dependencia direccional y espectral. La absortividad direccional espectral,
, de una superficie se define como la fraccin de la intensidad espectral incidente en la

direccin de y que la superficie absorbe. De aqu:

En esta expresin, ignoramos cualquier dependencia de la absortividad sobre la temperatura de la
superficie. Tal dependencia es pequea para la mayora de las propiedades radiativas espectrales.
Est implcito en los resultados anteriores que las superficies pueden exhibir una absorcin selectiva
con respecto a la longitud de onda y a la direccin de la radiacin incidente. Para la mayora de los
clculos de ingeniera, sin embargo, es deseable trabajar con propiedades superficiales que representen
promedios direccionales. Definimos por tanto una absortividad hemisfrica espectral
como:

Que tambin se puede expresar como:

Por tanto,
depende de la distribucin direccional de la radiacin incidente, as como tambin de la

longitud de onda de la radiacin y de la naturaleza de la superficie absorbente.
La absortividad hemisfrica total, , representa un promedio integrado sobre la direccin y la longitud
de onda. Se define como la fraccin de la irradiacin total absorbida por la superficie:
54

En consecuencia, , depende de la distribucin espectral de la radiacin incidente, as como de su
distribucin direccional y de la naturaleza de la superficie de absorcin. Advierta que, aunque , es
aproximadamente independiente de la temperatura superficial, no se puede decir lo mismo en cuanto a
la emisividad hemisfrica total, .
3.3.9. Ley de Kirchhoff

Considere un recinto isotrmico grande de temperatura superficial T
s
, dentro del cual estn confinados
varios cuerpos pequeos (Figura 2.32). Como estos cuerpos son pequeos en relacin con el recinto,
tienen una influencia insignificante sobre el campo de radiacin, que se debe al efecto acumulado de
emisin y reflexin por la superficie del recinto. Recuerde que, sin importar sus propiedades radiativas,
tal superficie forma una cavidad de cuerpo negro. En consecuencia, sin importar su orientacin, la
irradiacin experimentada por cualquier cuerpo en la cavidad es difusa e igual a la emisin de un cuerpo
negro a T
s

Bajo condiciones de estado estable, debe existir equilibrio trmico entre los cuerpos y el recinto. De
aqu,
y la transferencia netra de energa a cada superficie debe ser cero. Al aplicar un

balance de energa a una superficie de control alrededor del cuerpo 1, se sigue que:

o de la ecuacin:

Como este resultado se debe aplicar a cada uno de los cuerpo confinados, obetenemos:

Esta relacin se conoce como la ley de Kirchhoff. Una implicacin principal es que como
. Por tanto, ninguna superficie real puede tener una potencia emisiva que escapa la de una
superficie negra a la misma temperatura, y se confirma la nocin del cuerpo negro como un emisor
ideal.
55

.
Figura 2.32: Intercambio radiativo en un recinto isotrmico
De la definicin de la emisividad hemisfrica total, una forma alternativa de la ley de Kirchhoff es:

De aqu, para cualquier superficie en el recinto:

Es decir, la emisividad hemisfrica total de la superficie es igual a su absortividad hemisfrica total.
Sin embargo, se deben recordar las condiciones restrictivas inherentes en su derivacin. En particular, se
supone que la irradiacin superficial corresponde a la emisin desde un cuerpo negro a la misma
temperatura que la superficie.
La derivacin anterior se puede repetir para condiciones espectrales, y para cualquier superficie en el
recinto se sigue que:

Ecuacin que se aplica si la irradiacin es difusa o si la superficie es difusa. Una forma de la ley de
Kirchhoff para la cual no hay restricciones incluye las propiedades espectrales.

Esta igualdad siempre es aplicable pues
son propiedades superficiales inherentes. Es decir,

respectivamente, son independientes de las distribuciones espectral y direccional de la radiacin
emitida e incidente.
3.3.10. Reflectividad de una superficie

Considere la distribucin espacial de radiacin reflejada por una superficie. Cuando la radiacin
incidente, llega a la superficie en forma de un lpiz estrecho (esto es, contenido en un ngulo slido
pequeo), se pueden dar dos distribuciones posibles para la radiacin reflejada. Estos dos casos son el
56

difuso y el especular. La reflexin difusa ocurre si, sin importar la direccin de la radiacin incidente, la
intensidad de la radiacin reflejada es independiente del ngulo de reflexin. Por el contrario, si toda la
reflexin es en la direccin de
2
(Figura 2.33), que es igual al ngulo incidente
1
, se dice que ocurre la
reflexin especular. Aunque ninguna superficie es perfectamente difusa o especular, la ltima condicin
se aproxima ms de cerca con superficies de espejo pulidas y la primera condicin mediante superficies
speras.

Figura 2.33: Reflexin difusa y especular
En general, la magnitud de la intensidad reflejada en una direccin en particular para una superficie
dada es una funcin de la longitud de onda y de la distribucin espacial de la radiacin incidente. La
intensidad de radiacin reflejada en una direccin en particular se puede calcular con la siguiente
ecuacin:

La situacin fsica es la que se muestra en la Figura 2.34.

Figura 2.34: Sistemas de coordenadas para la funcin de reflexin

Desde un punto de vista prctico, no se suele usar el esquema de la figura anterior. Un experimento
equivalente, es irradiar una superficie con un flujo monodireccional (con un ngulo slido pequeo).
Midindose la energa reflejada en cada direccin. Esta medida divida por la medida del ngulo slido
ser aproximadamente igual a la intensidad reflejada.

57

Figura 2.35: Representacin esquemtica de un experimento para la medida de la funcin de reflexin
Existen dos tipos de reflectividad hemisfrica. La reflectividad hemisfrica angular, que se obtiene de
recoger la energa reflejada en una direccin particular cuando la superficie es irradiada desde todas
direcciones. Y la reflectividad hemisfrica monocromtica, que se designa por
donde el
subndice i indica que la radiacin incidente tiene una direccin especfica. Esta ltima reflectividad se
define como la razn entre la energa radiante monocromtica reflejada en todas direcciones y el flujo
incidente de energa radiante dentro de un ngulo slido pequeo
i
. Se puede expresar como:

Examinando la ecuacin anterior podemos observar que
es una propiedad de la superficie.

La energa incidente se puede escribir en trminos de una intensidad incidente integrada a lo largo de
todo el hemisferio:

Y la reflectividad hemisfrica monocromtica es:

donde el subndice r en
se usa para especificar la radiacin reflejada en un una direccin

especfica. Tambin se puede expresar como:

Podemos ver que
, es dependiente de la distribucin angular de la radiacin incidente por lo

que no es una propiedad de la superficie. Para el caso especial de que la radiacin incidente sea difusa,
entonces se puede demostrar que:

58

3.3.11. Relacin entre absortividad, reflectividad y emisividad

Ahora se va a demostrar como nicamente es necesario conocer una propiedad, la reflectividad angular-
hemisfrica monocromtica, y podremos obtener tanto la absortividad como la emisividad para un
cuerpo opaco.
Considere una superficie localizada en un recinto isotermo, que se mantiene a una temperatura T. La
intensidad monocromtica en la direccin (,) desde un rea infinitesimal se divide en dos partes, una
reflejada y otra emitida, y debe ser igual a

Las intensidades reflejadas y emitidas son:

pero, si consideramos la intensidad recibida difusa, entonces podemos cambiar
por
. Si cancelamos
de cada trmino entonces :

Y por la ley de Kirchhoff:

La relacin anterior se puede considerar como un principio de conservacin de la energa. La energa
monocromtica incidente en cualquier direccin es reflejada y absorbida.
3.3.12. Superficies selectivas

Los colectores solares deben tener una gran absortividad para la radiacin sobre todo en el espectro de
la energa solar. Al mismo tiempo, los colectores pierden energa por una combinacin de mecanismos
incluyendo emisin de radiacin trmica desde la superficie absorbedora, por ello es deseable que la
emisividad a grandes longitudes de onda sea lo ms baja posible para reducir las prdidas. Ya que,
normalmente, la temperatura del absorbedor suele ser menor de los 200C lo que indica que
prcticamente el espectro de emisin del sol (superficie a 6000 K) casi no se solapa con el espectro de la
superficie absorbedora, el 98% de la radiacin solar extraterrestre tiene una longitud de onda menor a
3m, mientras que solo el 1% de la radiacin emitida por la superficie absorbedora est por debajo de
ese valor.
Por lo tanto para estas circunstancias lo deseable es tener una superficie cuya absortividad sea elevada
a longitudes de onda cortas, y que tenga poca emisividad a longitudes de onda ms larga. Este concepto
se puede ver en la Figura 2.36
59

Figura 2.36: Superficie selectiva ideal
En la prctica, la dependencia de la longitud de onda de
no se acerca a la curva ideal de la Figura

2.37. En la Error! No se encuentra el origen de la referencia. podemos ver algn ejemplo de la
selectividad en superficies reales. Como vemos, estas superficies no tienen un marcado valor crtico para
la longitud de onda (
c
) o unas propiedades uniformes para longitudes de onda corta o larga.

Figura 2.37: Reflectividad vs longitud de onda para diversas superficies reales

CAPTULO 3: Modelos numricos implementados. Resultados obtenidos.
60

CAPTULO 3: Modelos numricos implementados. Resultados
obtenidos.

En este captulo de resumen los modelos numricos simulados, as como los resultados de inters y
conclusiones obtenidas en cada caso. Debido a la complejidad se ha seguido un procedimiento basado
en la validacin progresiva de modelos simples y metodologa de procesado de resultados utilizando
modelos simples, para llegar a un modelo final completo del captador solar y simulado en condiciones
muy similares a las reales.

En primer lugar se ha simulado un tubo simple sin pared al que se le aplica la condicin de contorno de
flujo de calor uniforme a lo largo de toda su superficie con objeto de validar la malla utilizada para la
regin del fluido caloportador y el procesado de resultados. Tambin ha servido para identificar los
mecanismos de transferencia de calor presentes. En segundo lugar se ha simulado un modelo formado
por un tubo simple unido al absorbedor mediante un cordn de soldadura, el objetivo en este caso era
comprobar la diferencia que supone en cuanto a la transferencia de calor, aplicar el flujo de calor en la
pared de tubo o en el absorbedor. Este modelo tambin ha servido para comparar los resultados
obtenidos con dos condiciones de contorno diferentes; flujo de calor uniforme y generacin interna de
calor.

El tercer modelo simulado incorpora los volmenes de aire superior e inferior y la cubierta de vidrio, es
decir, es un modelo de captador completo, pero para un solo tubo. Previamente se haba estudiado el
mallado necesario en la cavidad de aire superior entre el absorbedor y el vidrio para obtener resultados
suficientemente exactos al simular la conveccin natural. El objetivo de estas simulaciones era
comprobar si con un modelo simple, sin radiacin exterior y con un solo tubo, y en consecuencia de bajo
coste computacional, se poda predecir la curva de eficiencia trmica del captador.

Finalmente, sobre este tercer modelo se ha aplicado un modelo de carga solar y se han comparado los
resultados con el modelo anterior. Una vez validado el modelo de carga solar se ha aplicado a un
captador completo formado por 9 tubos de geometra estndar y se ha obtenido y comparado la curva
de eficiencia trmica con la de un captador estndar. El modelo completo validado permitir obtener
informacin acerca del campo de temperaturas en el absorbedor y de la evolucin de las propiedades
fluidas y parmetros adimensionales en el captador que resultara muy complejo obtener
experimentalmente.

61

1. Tubo horizontal con flujo de calor uniforme

El modelo que se presenta a continuacin, es un modelo muy simplificado que nicamente consta del
interior de un conducto idntico al que se encuentra en un colector solar. Sobre la pared del mismo se
ha impuesto una condicin de contorno de flujo de calor uniforme para observar cmo se produce la
transferencia de calor en el mismo.
Hay que resaltar que este modelo dista mucho del caso real, ya que por una parte el flujo de calor
proviene de una condicin impuesta artificialmente, mientras que realmente son los fenmenos de
radiacin y conveccin los que determinan el flujo calor a travs del absorbedor. Por otra parte, el flujo
de calor no se distribuye uniformemente sino que se conduce a travs de la soldadura lo que hace que
la temperatura del conducto se mucho mayor en esta zona que en la parte inferior, como
posteriormente se ver.
Por tanto el, nico objetivo de estas simulaciones es el de validar la malla elegida para el conducto y el
procesado de resultados para el clculo del Nmero de Nusselt local (
). Para ello se ha comparado

con las correlaciones presentadas en el captulo anterior tanto para el
con la distancia adimensional
, y para el Nusselt con flujo desarrollado (
) frente al nmero de Rayleigh (Ra).

1.1. Definicin de la geometra y mallado en GAMBIT

Como hemos dicho la geometra es muy simple, y consiste simplemente en un conducto de 7 mm de
dimetro por 1.83 metros de largo:

Figura 3.1: Geometra para el modelo del tubo horizontal
La malla transversal utilizada a lo largo del conducto de la Figura 3.1, es la siguiente:
62

Figura 3.2: Malla transversal en el tubo horizontal
Como vemos la malla se ha realizado procurando que sea lo ms regular posible a lo largo de toda la
seccin, hay que tener en cuenta que esta malla ha sido la misma que se ha utilizado para los casos
donde se modela el colector completo, es decir, con la soldadura incluida por lo que la concentracin de
nodos a lo largo del eje Y es mayor.
La malla anterior se ha espaciado cada 10 mm a lo largo del conducto, lo que hace que tengamos nodos
suficientes para caracterizar el comportamiento del fluido a lo largo de la direccin longitudinal (Z). Las
caractersticas de la misma son:
Grid Size
Level Cells Faces Nodes Partitions
0 98820 302856 105432 12

Memory Usage
Cells Faces nodes objps edges
Number Used 98820 302856 105432 5 0
Mbytes Used 36 32 6 0 0
Number Allocated 98820 302856 105432 1024 0
Mbytes Allocated 36 32 6 0 0
Tabla 3.1: Caractersticas malla modelo tubo simple

63

1.2. Hiptesis de simulacin en FLUENT.
1.2.1. Hiptesis de simulacin
- Solver: Pressure Based Segregated (SIMPLE)
Es el indicado para flujos con bajas velocidades e incompresibles. Con este algoritmo las ecuaciones que
gobiernan el problema se resuelven de manera secuencial, y dado que dichas ecuaciones son no
lineales y adems estn acopladas, el bucle se debe realizar de forma iterativa hasta alcanzar la
convergencia.
En el algoritmo segregado, las ecuaciones individuales que rigen las variables que son solucin (u ,v, w,
p, T, k, etc.) se resuelve una tras otra, es decir, mientras se resuelve se desacopla del resto de
ecuaciones. Esto hace que este mtodo no necesite demasiados recursos de memoria ya que las
ecuaciones discretizadas se van guardando una cada vez. Sin embargo, para alcanzar la convergencia se
necesita mayor nmero de iteraciones, debido al desacoplamiento durante la resolucin las ecuaciones.
Los pasos en el algoritmo son los siguientes:

Figura 3.3: Algoritmo segregado. User's Manual, Chapter 25, FLUENT.

El algoritmo segregado que se ha escogido en concreto ha sido el SIMPLE, por ser un problema no
demasiado complejo, en el que el flujo es laminar.

64

- Esquema de discretizacin: Segundo Orden
- Simulaciones estacionarias
- Modelo Viscoso: Laminar
- Ecuacin de la energa: resuelta
- Gravedad: (0,-9.81,0) (m/s
2
)
- Fluido: Agua
- Propiedades del fluido: (Ver 1.2.4), funcin de la temperatura
- Velocidad de entrada: la velocidad se establece fijando el gasto msico por unidad de superficie
de absorbedor. Este valor se ha tomado como
, ver 1.2.3.
- Flujo de calor a travs de la superficie del conducto, ver 1.2.3
1.2.2. Criterios de convergencia
Despus de discretizar, la ecuacin de conservacin para una variable general en una celda P se puede
escribir como:

Donde
es el coeficiente central,
son los coeficientes de influencia de las celdas vecinas, y b es la

contribucin de la parte constante del trmino fuente
donde
, y de las condiciones de
contorno. En la ecuacin (3.1):

Por lo tanto el residuo para
se define como el desequilibrio en la ecuacin (3.1) sumado para todas

las celdas. A este se le denomina residuo no-escalado:

Pero, normalmente es difcil juzgar la convergencia examinando este tipo de residuos, por lo que no se
suele utilizar los residuos sin escalar. FLUENT los escala usando un factor de escala representativo del
flujo de la variable , a travs del dominio, es decir:

Se ha considerado que una simulacin ha convergido cuando los residuos se estabilizan con un valor
inferior a 10
-5
para la ecuacin de continuidad y 10
-7
para la ecuacin de la energa.
Por otra parte adems de tener en cuenta los residuos para evaluar la convergencia de las simulaciones
tambin se han vigilado los flujos de masa y energa a travs de las fronteras que forman el modelo, de
manera que es posible establecer un balance energtico y msico en el dominio. En estas simulaciones
al no existir ningn parmetro de generacin interna ni tampoco sumideros o fuentes, cuanto menor
sea la diferencia entre el gasto entrante y saliente mayor ser la convergencia del caso.
El lmite puede establecerse como 10
-10
kg/s y 10
-4
W, para las ecuaciones de continuidad y conservacin
de la energa respectivamente.

65

1.2.3. Condiciones de contorno

A continuacin se muestra un resumen de las condiciones de contorno:
Frontera Condicin de contorno
Cara_interior_tubo Wall (flujo de calor uniforme)
Entrada_de_fluido Velocity-inlet
Salida_de_fluido Pressure-outlet
Tabla 3.2: Condiciones de contorno en tubo simple
- Cara_interior_tubo
Esta frontera es de tipo Wall, en ella especificamos un flujo de calor constante. El flujo de calor al que se
somete el conducto depende de la radiacin que incide sobre el captador solar. Por ello se ha hecho un
barrido de temperaturas de entrada para dos flujos de calor diferentes.
Podemos decir que en un colector la mayor parte de la energa que absorbe el absorbedor se transmite
por conduccin al tubo. Por lo tanto lo que se ha hecho es suponer un flujo de calor, dentro de uno
valores razonables sobre la placa absorbedora (W/m
2
) y trasponer dicho flujo a la superficie del tubo.
Se han simulado para unos valores de 500 W/m
2
y 1000 W/m
2
, esto quiere decir que si la superficie de
absorbedor que corresponde a cada conducto es de 0.123x1.83 m, los flujos aplicados a la pared son:

Esto hace que, para los flujos anteriores, en el tubo se tengan:

Figura 3.4: Condicin de contorno de flujo de calor uniforme. FLUENT

66

Al definir el flujo de calor como condicin de contorno en FLUENT, estamos especificando el flujo de
calor en dicha superficie. FLUENT usa la ecuacin
y el dato que se
proporciona para determinar la temperatura de la pared del la celda fluida adyacente como:

donde, el coeficiente convectivo del lado del fluido se calcula en base a las condiciones locales del flujo.
- Entrada_de_fluido
En esta condicin es de tipo velocity-inlet, se emplea para definir la velocidad y otras propiedades
escalares en la zona de entrada de fluido. Hay que tener en cuenta que se va modelar nicamente un
flujo de tipo laminar, con valores de Reynolds hasta un mximo de 2300.
FLUENT al definir un flujo entrando en el dominio fsico, las componentes de la velocidad y los escalares
definidos (Temperatura) se utilizan para calcular el gasto msico y los flujos de cantidad de
movimiento y energa.
El gasto msico entrando en una celda adyacente a una condicin de velocity intlet se calcula como:

Tngase en cuenta que solo la componente normal de la velocidad sobre la superficie de entrada es la
que contribuye a la entrada de gasto msico.
En cuanto a los parmetros y caractersticas que es necesario especificar se encuentran:
- Se elige el mtodo para especificar la direccin del flujo, en nuestro caso: magnitud y direccin
normal a la superficie de entrada.
- En nuestro caso, el modelo se encuentra fijo, de tal modo que la velocidad se especificara en
trminos absolutos.
- Se especifica el valor absoluto del vector velocidad
- Se tiene que especificar lo que Fluent denomina Outflow Gauge Pressure, dicho valor es
necesario para que en el caso de que durante proceso de resolucin del caso exista reflujo, la
cara por la cual retroceda el fluido se tratara como una condicin de pressure- outlet, con
este valor.
- Por ltimo tambin se deber de tomar muy en cuenta la temperatura a la que fluido entra en
el dominio computacional.
67

Figura 3.5: Propiedades en la entrada de flujo al dominio
En las simulaciones realizadas con este modelo se ha hecho un barrido de temperaturas de entrada para
cada uno de los flujos de calor. Las temperaturas analizadas han sido de 20, 50 y 70 C. Para este
modelo no se han realizado ms simulaciones ya solo vamos a comparar los parmetros adimensionales,
para modelos ms complejos, se ha hecho un estudio ms exhaustivo.
La velocidad de entrada se define fijando el flujo msico por rea de absorbedor
a 0.02
. Por tanto la velocidad se calcula en funcin de la superficie que le corresponde a cada

conducto, es decir,
, y segn la densidad de entrada, que a su vez es

funcin de la temperatura. Es decir:

Las simulaciones realizadas se han hecho para el flujo msico nominal, y para el doble y la mitad del
mismo. En el punto se detallan, los casos resueltos.
Salida_de_fluido
La salida de fluido se modela con otro tipo de condicin de contorno, pressure outlet. Requiere la
especificacin de la presin esttica en la superficie de salida del tubo. FLUENT nicamente utiliza
este valor mientras que el flujo sea subsnico, tal y como es el caso.
- La presin introducida es de 300.000 Pascales, y se trata de una presin relativa. Aparte
activamos la opcin de Radial Equilibrium Distribution, de manera que la presin anterior
solo se aplica a la posicin de radio mnimo (relativo al eje de rotacin) de la frontera. La
presin esttica en el resto de la zona se calcula a partir del supuesto de que la velocidad radial
es despreciable, de modo que el gradiente de presin viene dado por:
68

- Donde r es la distancia desde el eje de rotacin y
es la velocidad tangencial. Por supuesto

esta opcin solo se puede utilizar en el caso de que la velocidad rotacional sea 0, como es el
caso del colector solar.

Figura 3.6: Parmetros a definir en la salida del fluido

1.2.4. Materiales
En este caso el dominio solo est formado por el volumen fluido. El fluido de trabajo para las
simulaciones que se han hecho en este proyecto ha sido el agua, cuyas propiedades se han introducido
en funcin de la temperatura.
Esto ltimo es muy importante ya que para los casos simulados hay que tener en cuenta los efectos de
la flotacin, que son una consecuencia de la accin combinada del campo gravitatorio con un gradiente
de densidades provocado a su vez por un gradiente de temperaturas.
Las funciones utilizadas fueron las siguientes:

69

1.2.5. Simulaciones realizadas
La distribucin de casos se ha hecho en base a 3 parmetros, la temperatura de entrada, el flujo de calor
y el caudal.
Caso Temperatura de
entrada (K)
Flujo de calor (W/m
2
) Velocidad de entrada
(m/s)
300_1Q_600 300 3355.89 0.117394411
300_1.5Q_600 300 3355.89 0.176091617
300_0.5Q_600 300 3355.89 0.058697206
300_1Q_1000 300 5593.16 0.117394411
300_1.5Q_1000 300 5593.16 0.176091617
300_0.5Q_1000 300 5593.16 0.058697206
323_1Q_600 323 3355.89 0.118380874
323_1.5Q_600 323 3355.89 0.177571312
323_0.5Q_600 323 3355.89 0.059190437
323_1Q_1000 323 5593.16 0.118380874
323_1.5Q_1000 323 5593.16 0.177571312
323_0.5Q_1000 323 5593.16 0.059190437
343_1Q_600 343 3355.89 0.119626691
343_1.5Q_600 343 3355.89 0.179440036
343_0.5Q_600 343 3355.89 0.059813345
343_1Q_1000 343 5593.16 0.119626691
343_1.5Q_1000 343 5593.16 0.179440036
343_0.5Q_1000 343 5593.16 0.059813345
Tabla 3. 3: Simulaciones realizadas para tubo simple
Como vemos se ha impuesto dos valores extremos en el flujo de calor a travs de la pared,
normalmente en un captador solar, el flujo de calor se encontrar entre ellos. En cuanto a la
temperatura de entrada del fluido en el conducto se han fijado tres valores nicamente ya que el
objetivo final es el de validar le procesado de resultados, en los casos donde se trabaje con modelos
ms completos se realizar un barrido de temperaturas ms fino, con el propsito de conocer
finalmente la eficiencia del captador.

1.3. Metodologa de procesado de resultados
1.3.1. Evaluacin del Nmero de Nusselt
El nmero de Nusselt trata de comparar entre la transmisin de calor por conduccin y por conveccin,
para el dimetro del conducto se define como:
. Vamos a evaluarlo por secciones

transversales, es decir, promediando la temperatura de pared, el flujo de calor, la temperatura de fluido
y la conductividad radialmente en la seccin del conducto.
PROCESO ANTERIOR
Antes de iniciar el proceso seguido se va a realizar un breve recordatorio de la forma de procesar los
resultados en el PFC anterior. En aqul caso el Nusselt no se promediaba en una seccin concreta si no
en un volumen de control, definido por dos secciones cada centmetro. Es decir obtenamos la
temperatura del fluido a la entrada y la salida del volumen de control, calculbamos una temperatura
promedio de pared en ese intervalo, asimismo promedibamos el calor especfico y obtenamos el
coeficiente convectivo en la pared. El proceso es el siguiente:

70

donde
, es la temperatura promedio en el volumen de control y
es cualquier propiedad del fluido

evaluada (Calor especfico, conductividad, densidad). A su vez
, son la temperatura de entrada

promediada en rea y la temperatura de salida. El promediado en rea es el siguiente:

donde puede ser una propiedad cualquiera, A es el rea total de la superficie y
, son las reas de las

celdas que componen dicha superficie.
El coeficiente convectivo se obtendr como:

donde S es la superficie que rodea el volumen de control
, siendo
, la longitud entre la
entrada y salida, como se ha dicho, 1 centmetro. Y es el gasto msico que al igual que las
propiedades del fluido se evala en cada seccin haciendo uso de la herramienta Surface-Integrals, se
trata de una cantidad calculada sumando el producto de la densidad con el producto escalar del vector
normal a la superficie de la celda y el vector velocidad:

Por tanto ya somos capaces de obtener el coeficiente convectivo en cada volumen de control, y por lo
tanto en nmero de Nusselt.
La gran deficiencia del proceso era la de obtener el calor que radialmente recibe el fluido
indirectamente. Es decir utilizar el calor especfico suponiendo un modelo de flujo incompresible trae
consigo un error. Por ello para el procesado usado en esta nueva revisin se ha evaluado el calor en
cada seccin. La descripcin de esta nueva forma de trabajar con las simulaciones se hace a
continuacin.
PROCESADO ACTUAL
Al igual que en el procesado anterior, se van a crear una serie de superficies en la zona del tubo
espaciadas cada centmetro que nos permitan observar la temperatura del fluido, as como sus
propiedades, tambin podremos comparar perfiles de velocidad y temperatura a lo largo del tubo, as
como, recirculaciones (importante en el caso de conveccin mixta).
71

Figura 3.7: Secuencia de superficies creadas a lo largo del conducto
A continuacin pasamos a evaluar las siguientes propiedades y campos del fluido: temperatura, presin
esttica, velocidad, nmero de Reynolds, densidad, viscosidad, nmero de Prandtl, conductividad
trmica. Estos valores son promediados en cada seccin, dicho promedio se puede hacer a travs de dos
herramientas: Mass-weighted-average o Area-weighted-average, que se enuncian:
- Promediado en rea: es el promedio de una cantidad que se calcula dividiendo el sumatorio
del producto de la variable seleccionada por el rea de cada celda que forma la superficie;
entre el rea total de la superficie:

- Promediado en masa: se calcula dividiendo el sumatorio del producto de la variable
seleccionada por el valor absoluto del producto escalar del vector normal al rea de la celda y
el vector velocidad entre el valor absoluto del producto escalar de los vectores de la cara y
velocidad del flujo que la atraviesa:

Por ltimo falta por conocer el calor recibido en cada seccin y la temperatura de pared (
para
poder conocer el nmero de Nusselt. Estos dos parmetros se obtienen de forma distinta a los
anteriores. A travs de la opcin XY plot, podemos visualizar el valor en nodos de una determinada
variable. Los pasos para definirlo son:
- Para graficar una variable en el eje y como una funcin de la posicin segn en el eje x, se
marca Position on X Axis, y elegimos el parmetro a representar.
- Se elige la direccin del vector en funcin del cual vamos a graficar la variable. En este caso
como se desea que sea a lo largo del conducto, la direccin coincide con el eje Z. La direccin es
(0,0,1).
72

- Se elige la superficie donde se encuentra los datos que queremos representar. Para la
simulacin es la cara interior del conducto.

De esta manera, para la temperatura de pared y el flujo de calor, se obtiene archivos donde se muestran
cada centmetro(es cada 1 cm porque la malla est definida longitudinalmente de esta forma) una serie
de valores correspondientes a los valores de esas dos variables en los nodos que forman la superficie
interior del conductor. Por lo que si promediamos de modo acimutal obtenemos un valor medio para
esa seccin.

Figura 3.8: Temperatura de la pared interior del conducto
Por lo que si despreciamos el flujo de calor axial, se puede establecer como nica componente
importante del flujo de calor, la radial. Segn la ley de enfriamiento de Newton:

Se puede despejar el coeficiente convectivo, y evaluar en Nusselt en cada seccin del conducto como:
.
Por ltimo hay que tener en cuenta que adems de la variacin de las propiedades del fluido en sentido
longitudinal, si la diferencia de temperaturas entre la pared del tubo y el fluido es sustancial se harn
patentes las variaciones radiales en las propiedades fsicas del fluido, fundamentalmente en la capa
lmite. En lquidos la variacin que hay que tener en cuenta es la de la viscosidad dinmica, puesto que
el resto de propiedades fsicas (conductividad trmica, calor especfico y densidad) varan poco con la
temperatura.
Para tener en cuenta esto, el Nusselt calculado se corrige para considerar la influencia de la variacin de
la viscosidad del fluido en direccin radial. Se emplea un factor de correccin en la forma
donde
es la viscosidad del fluido a la temperatura media del fluido en la seccin y
es la viscosidad del
fluido a la temperatura de la pared del tubo, el exponente n para rgimen laminar y calentamiento
n=0.14.

73

Solo queda por fijar el mtodo de promediado para las variables y propiedades del fluido. Se probaron
ambas posibilidades obtenindose los siguientes resultados:

Figura 3.9: Nusselt en conveccin forzada. Datos obtenidos con mass-weighted average

Figura 3.10: Nusselt en conveccin forzada. Datos obtenidos con area-weighted average

A la vista de los resultados se opta por hacer un procesado Mass-Weighted-Average ya que coincide
mejor con los datos de referencia.
1
10
100
0.0001 0.001 0.01 0.1 1
N
u
x

x*
Petukhov y Polyakov [1988];
conveccin forzada
Datos Numricos Re500 Qcte
Conveccin forzada
1
10
100
0.0001 0.001 0.01 0.1 1
N
u
x

x*
Petukhov y Polyakov
[1988]; conveccin forzada
Datos Numricos Re500
Qcte conveccin forzada
74

1.3.2. Evaluacin del coeficiente de friccin

Hay que tener en cuenta que en el conducto con flujo de calor constante aplicado en la pared, las
variaciones de las propiedades del fluido a lo largo de eje Z son importantes, lo que hace que el
coeficiente de friccin vare longitudinalmente
Para obtener un valor adecuado del factor de friccin, se deben de tener en cuentas estas variaciones,
para ello al igual que se ha hecho en el caso del nmero de Nusselt, se ha dividido el conducto en
tramos de un centmetro de espesor, calculndose el factor de friccin en cada una de estas secciones.
El clculo del factor de friccin de cada tramo discreto se calcula tomando como cada de presin la
diferencia de presiones entre las dos caras que conforman los tramos. Como longitud se toma la
correspondiente al tramo discretizado, que en todos los casos ser de 1 centmetro y se evala la
velocidad y la densidad del fluido utilizando el promediado de Area-weighted-average.

Figura 3.11: Coeficiente de friccin frente a la distancia. Tubo liso horizontal con flujo de calor 1000 W/m2 y
temperatura de entrada de 300 K.

Como se ve en la Figura 3.11, existe una zona de crecimiento de la capa lmite hidrodinmica hasta
llegar a un valor ya coincidente con la solucin analtica. A partir de ah se llega a la zona de flujo
plenamente desarrollada donde se observa una disminucin progresiva del coeficiente de friccin
consecuencia de la variacin de la temperatura a lo largo del conducto, ya que a medida que el flujo se
calienta el coeficiente de friccin disminuye.
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0 0.5 1 1.5 2
C
o
e
f
i
c
i
e
n
t
e

d
e

f
r
i
c
c
i
n
,

f

Distancia (m)
Simulacin Caso_26_1.5Q_1000
f=64/Re
75

Tal y como se ha adelantado, el promediado de la velocidad a diferencia de lo que ocurra en el nmero
de Nusselt, da mejores resultados si usamos el Area-weighted-avg para calcular la velocidad
promedio.
Esto es as, ya que al evaluar el nmero de Nusselt estamos realizando un balance de energa trmica,
por lo que al ser una magnitud extensiva si influye la masa de fluido considerada en el clculo de la
temperatura promedio del fluido. En cambio al calcular el factor de friccin, estamos realizando un
balance de energa mecnica, en concreto de energa de presin y energa cintica, que al ser
magnitudes intensivas no se ve influenciado por la masa total de fluido considerada y la velocidad
promedio de la seccin haciendo el area-weighted-avg es independiente, aunque deberan dar
resultados ms parecidos.
1.4. Validacin del modelo y metodologa del procesado
1.4.1. Resultados en la transmisin de calor
En este apartado se van a presentar los resultados alcanzados en las simulaciones realizadas con tubo
horizontal liso con flujo de calor constante. Para ello se va a comparar con las correlaciones presentadas
en el captulo anterior, atendiendo principalmente a los efectos de la conveccin mixta, y su relacin
con el nmero de Nusselt.
Como se ha comentado en el punto anterior, los casos se han distribuido de manera que se cubran el
rango de caudales y temperaturas donde trabajan normalmente los captadores solares. Hay que tener
claro que el objetivo final del proyecto es obtener un modelo que permita obtener una curva de
eficiencia lo ms ajustada a la realidad, es por esto por lo que para algunos casos en los que la
temperatura de entrada es alta y el caudal tambin los valores de Reynolds resultan algo superiores a
los de rgimen laminar.

0 500 1000 1500 2000 2500 3000 3500
N de Reynolds
Caso_70_1Q_1000
Caso_70_1Q_600
Caso_70_0.5Q_1000
Caso_70_0.5Q_600
Caso_26_1Q_1000
Caso_26_1Q_600
Caso_26_1.5Q_1000
Caso_26_1.5Q_600
Caso_26_0.5Q_600
Caso_50_1Q_1000
Caso_50_1Q_600
Caso_50_1.5Q_1000
Caso_50_1.5Q_600
Caso_50_0.5Q_1000
Caso_50_0.5Q_600
Figura 3.12: Reynolds para los casos simulados para tubo liso horizontal con flujo de calor constante
76

En cada simulacin la temperatura va variando a lo largo del conducto, esto hace que las propiedades
del fluido cambien longitudinalmente y por tanto en cada superficie donde calculamos el nmero de
Nusselt tenemos diferentes nmeros de Reynolds. En la figura anterior se muestran la variacin de este
parmetro adimensional para cada caso simulado.
En las simulaciones se ha buscado determinar el nmero de Nusselt local promediado
circunferencialmente para cada punto separado 1 cm a lo largo del tubo,
, que adems de depender

de la posicin en el conducto depende tambin de la conveccin natural,
. Las
propiedades fsicas como ya se ha explicado, van variando adems de longitudinalmente, radialmente,
por lo que todos los nmeros de Nusselt calculados se han corregido por el factor
para obtener
una correlacin libre de la variacin de las propiedades fsicas en una seccin transversal.
Los resultados se han comparado con las correlaciones de Petukhov y Polyakov [1988], que han
estudiado la transmisin de calor en tubos lisos horizontales en conveccin mixta, dando como
resultado:

Figura 3. 13: Resultados para la conveccin mixta en tubo liso con flujo de calor constante. Nmero de Nusselt local, en funcin del Rayleigh y de
la distancia adimensional x*. Comparacin con la correlacin de Petukhov y Polyakov [1988].
1
10
100
0.0001 0.001 0.01 0.1 1
Rayleigh=2.78E+05
Rayleigh=1.12E+06
Petukhov y Polyakov [1988],
Conveccin Mixta, Ra=2.78E+05
Conveccin forzada
Conveccin Mixta, Ra=1.12E+06
77

Se puede ver en la Figura 3. 13 que los resultados obtenidos no se ajustan del todo bien en la regin de
entrada donde como sabemos se da un crecimiento simultaneo de la capa lmite trmica e
hidrodinmica. La explicacin de este efecto es que el trabajo de Petukhov y Polyakov fue realizado para
flujo con la capa lmite hidrodinmica ya desarrollada. Adems se puede observar como a medida que
crece la distancia adimensional y , antes de alcanzar un flujo plenamente desarrollado los resultados
numricos se encuentran ms cercanos a la correlacin, esto se justifica en que la longitud de entrada
hidrodinmica es menor que la trmica. Por tanto, en la regin de entrada trmica, una vez
desarrollada la capa lmite hidrodinmica los resultados siguen la tendencia de la correlacin, que se
puede describir como
.
Como era de esperar a medida que el nmero de Rayleigh se incrementa, los efectos de flotacin se
hacen mayores, y la longitud de entrada trmica se reduce. En ninguno de los casos analizados se llega a
tener una conveccin forzada pura durante todo el tubo, pero a medida que el flujo de calor aportado
aumenta ,el nmero de Reynolds se incrementa y la temperatura de entrada se reduce, nos acercamos
ms a la situacin de conveccin forzada, por lo que el nmero de Nusselt disminuye y como se ver, la
componente de la velocidad en direccin radial (recirculaciones) se va atenuando conforme la diferencia
de temperaturas en la pared del tubo se reduzca.
En la zona de flujo plenamente desarrollado los resultados se ajustan con un error del entorno al 1.5%
(Figura 3.15), hay que decir que en aquellos casos donde el nmero de Reynolds est por encima de los
2700, como era de esperar, los resultados no coinciden con la correlacin, porque se est superando el
lmite de rgimen laminar. En la siguiente figura podemos apreciar esto ltimo:

Figura 3.14: Comparacin Nusselt en tubo horizontal con flujo de calor uniforme para alto nmero de Reynolds.

1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Dixtancia adimensional, x*
Simulacin, Re=3600
Conveccin forzada
Conveccin mixta
78

Figura 3.15: Error relativo para el nmero de Nusselt en tubo simple.
Adems de comparar el Nusselt con la distancia adimensional, es importante resaltar que en conveccin
mixta, como el flujo plenamente desarrollado ocupa prcticamente la totalidad del conducto, no est
dems realizar una comparacin de el Nusselt plenamente desarrollado
, frente al nmero de
Rayleigh. Para calcular
, se ha hecho el promedio de los valores Nusselt cuyo valor asociado de x*

se encuentra por encima de 0.01. Para esta comparacin se han desechado aquellas simulaciones con
Re>2500.

Figura 3.16: Nmero de Nusselt plenamente desarrollado, en funcin del nmero de Rayleigh Ra. Comparacin
con las correlaciones de Petukhov y Polyakov [1988] y A.Garca [2006]
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
E
r
r
o
r

r
e
l
a
t
i
v
o

(
%
)

Distancia adimensional, x*
Error Relativo Ra=2.78E+05
Error Relativo Ra=1.12E+05
1
10
1.00E+05 1.00E+06 1.00E+07
N
m
e
r
o

d
e

N
u
s
s
e
l
t

d
e
s
a
r
r
o
l
l
a
d
o
,

N
u

Nmero de Rayleigh, Ra
Simulaciones Tubo
Simple
Petukhov y Polyakov
[1988]
A.Garca [2006]
79

Como vemos, los datos coinciden bastante bien, sobre todo con la correlacin de Petukhov y Polyakov
[1988]. , sin embargo hay mayores discrepancias con el trabajo de A.Garca [2006], la causa es que el
trabajo de de los primeros fue realizado con simulaciones para Prandtl del orden 2-7, muy similares a los
de las simulaciones sin embargo los casos experimentados con A.Garca, fueron realizados para nmeros
de Prandtl mayores.
Por ltimo, antes de presentar el siguiente modelo, vamos visualizar el flujo en el interior del conducto
donde podremos comprobar los efectos de flotacin. Para ellos vamos a tomar diferentes secciones a lo
largo del tubo para diferentes casos.
Se va a mostrar el campo de temperaturas junto con los vectores de velocidad para el caso simulado con
menor nmero de Rayleigh, Ra=2.78x10
5
., y por lo tanto el que tiene mayor longitud de entrada
trmica:

Figura 3.18: Conducto horizontal liso, con flujo de calor
uniforme, Ra=2.78E+05, Seccin z=0.3 m. Campo de
temperaturas (K) y vectores de velocidad.
80

Se puede apreciar como a lo largo del tubo el agua se va calentando, y debido a los efectos de la
flotacin se dan dos recirculaciones en direccin del eje haciendo que el flujo caliente ascienda por la
parte exterior y descienda un flujo ms frio por la parte central. Como se ve en la seccin cercana a la
entrada, el gradiente de temperaturas en la pared es muy uniforme lo que hace tengamos un nmero
de Nusselt propio de la conveccin forzada como se ve en la Figura 3. 13. A medida que avanzamos el
gradiente se incrementa y vemos como a pesar de que el flujo de calor es uniforme las temperaturas
ms altas se dan en la parte superior, por la conveccin mixta.

temperaturas (K) y vectores de velocidad
81

Si comparamos las cuatro figuras anteriores con las del caso de Ra=2.78E+05, podemos comprobar que
en este caso ya en la seccin de entrada z=0.3 metros tenemos un perfil de temperaturas muy parecido
al que tenemos en el caso anterior para un seccin ms avanzada. Lo anterior justifica que para un
Rayleigh mayor la zona de flujo plenamente desarrollado se alcance antes como veamos en la Figura 3.
13.
Tambin se puede ver como los componentes de velocidad en el caso de Ra=1.27E+06 son menores, si
los comparamos con la simulacin a Ra=2.78E+05. La causa radica en el nmero de Reynolds ya que en
el primero se encontraba en un rango de 1400-2000, mientras que en el segundo, 700-1200. Esta
diferencia hace que la distribucin de recirculaciones en el conducto difiera entre los casos, aunque en
la parte inicial el patrn es el mismo en la Figura 3.23, se puede ver como se dan otras dos
recirculaciones en la parte inferior del tubo haciendo que en la zona central el fluido tienda a ascender
algo que no ocurra en el caso anterior.
Se puede concluir esta apartado diciendo que, por una lado queda validado el procesado de nmero
Nusselt por la buena afinidad con los trabajos experimentales de Petukhov y Polyakov. Como era de
esperar para aquellos casos en los que Reynolds est muy por encima de 2500 el nmero de Nusselt ya
no se ajusta bien. Y adems, se ha comprobado el patrn de flujo dentro del conducto y su dependencia
con el nmero de Rayleight.

82

2. Modelo de tubo con absorbedor

Lo que se pretende hacer al crear este modelo es poder comparar los resultados con los de tubo simple.
Este modelo a diferencia del que se haca referencia en el apartado anterior adems del volumen
fluido, consta del tubo, la soldadura y el absorbedor.
Estas simulaciones todava distan mucho del modelo completo, pero nos sirven para poder distinguir las
diferencias en cuanto a la transmisin de calor hacia el fluido calo-portador, ya que en el caso anterior el
flujo de calor se mantena uniforme a lo largo de la pared del conducto mientras que en este caso el
flujo se impone a travs de una condicin de contorno sobre el absorbedor, y por lo tanto llega hasta el
conducto a travs de la soldadura. Por lo tanto el parmetro ms importante a estudiar en este caso es
el Nusselt en la pared del tubo.
Adems de esto, otro de los objetivos ha sido el de comparar la utilizacin de la condicin de flujo de
calor constante o imponer en el absorbedor un trmino de generacin interna de calor. El objetivo de lo
anterior es comprobar que efectivamente es indiferente la utilizacin de una u otra condicin, ya que
posteriormente en el modelo que implementa la cavidad de aire superior es necesario, para no
desacoplar el problema utilizar la generacin interna.
Antes de describir la malla, hay que resaltar, que a excepcin de la pared superior del absorbedor el
resto de superficies que delimitan el dominio computacional se encuentran aisladas por lo que las
prdidas son nulas, al igual que ocurra con el modelo de tubo simple.


La geometra es la siguiente:

Figura 3.25: Geometra modelo para 1 tubo simplificado

83

Las dimensiones que definen la geometra son:

Figura 3.26: Caractersticas geomtricas del tubo y aleta del absorbedor
Las dimensiones que tenemos que tener en cuenta son:
- Dimetro interior del conducto (D
i
): 7mm
- Dimetro exterior del conducto (D
e
): 8mm
- Espesor de la soldadura: (punto ms estrecho) 0.5 mm
- Ancho de la soldadura (g): 1.5 mm
- Espesor del absorbedor (
abs
): 0.5 mm
- Ancho del absorbedor (w): 123mm (se estudia un solo tubo)
- Longitud de los tubos del colector: 1830 mm

- Construccin de la malla
Se parti de la malla que se tena del fluido junto con el conducto:

Figura 3.27: Seccin del plano XY, malla para fluido y conducto.

84

Al igual que ocurra para el caso del tubo simple, se opto por generar una malla en el plano XY, para
despus crear el volumen extruyendo geometra y malla a la vez. Por ello a la imagen de la Figura 3.27,
se le aadi la parte de la soldadura junto con el absorbedor:

Figura 3.28: Detalle mallado de la soldadura.
Se puede ver en la figura anterior como la densidad de nodos en la parte de la soldadura es mayor
comparado con la parte central del conducto, la idea es que en esa zona los gradientes de temperatura
van a ser muy importantes si tenemos en cuenta que todo el flujo de calor hasta llegar al conducto
terminar discurriendo por un paso muy estrecho.
Adems del flujo de calor que pasa a travs de la soldadura, hay que tener en cuenta que toda ella
estar rodeada de aire, con lo que el salto trmico en la frontera lateral de la misma va a ser
considerable.
Como se aprecia en la Figura 3.28, el que la soldadura tenga mayor nmero de celdas hace que, si
queremos obtener una malla lo ms estructurada posible los nodos del absorbedor y del conducto que
unen ambas zonas deben estar ms juntos, por lo que incrementar la densidad de malla en la soldadura
trae consigo aumentar el coste computacional del modelo ya que es el punto central de toda la
geometra.
Una vez mallada la parte de la soldadura realizamos el mallado del absorbedor, teniendo presente en
todo momento intentar crear unas celdas lo ms regulares posibles, es decir, evitando que sean muy
oblicuas, o que tengan ngulos muy estrechos. Para ello, y para que el error numrico fuese menor, en
el absorbedor los nodos que se encuentran cercanos a la soldadura estn ms prximos unos de otros
dado que el gradiente de temperaturas va a ser mayor que hacia el extremo del modelo (Figura 3.29)
85

Figura 3.29: Mallado del absorbedor.

A modo de resumen se exponen el tamao de malla utilizado, as como el tipo de condiciones de
contorno utilizadas:
Grid Size
0 185196 599788 229300 12
Memory Usage
Number Used 185196 599788 229300 29 0
Mbytes Used 49 60 12 0 0
Tabla 3.4: Caractersticas malla de tubo con absorbedor

2.2. Hiptesis de simulacin y condiciones de contorno
2.2.1. Hiptesis de simulacin
- Esquema de discretizacin: Primer Orden
- Gravedad: (0,-9.81,0) (m/s
2
)
- Condiciones de operacin:
86

Figura 3.30: Condiciones de operacin para el modelo completo para 1 solo tubo con absorbedor.

Como se ha indicado en la introduccin los objetivos de este modelo son dos: por un lado analizar la
transmisin de calor en condiciones de flujo no uniforme en la pared del tubo, y por otro lado
comprobar se puede poner una condicin de generacin de energa interna o flujo de calor constante y
los resultados son similares.
Ya que cuando se model el captador completo sin carga solar en un primer lugar se pens en utilizar
una condicin de flujo de calor constante sobre el absorbedor. El hecho de probar con generacin
interna de calor o con flujo de calor es debido a las dificultades de acoplamiento de los problemas de
conveccin natural y radiacin que se dan entre la superficie del absorbedor y la cubierta de vidrio del
colector solar (correspondiente a un volumen de aire); y el problema de conduccin y conveccin mixta
que se dan desde el absorbedor, pasando por la soldadura hasta llegar al fluido caloportador que recoge
ese flujo calorfico.
Las dificultades surgieron por la forma de modelar la radiacin, ya que FLUENT no es capaz de asociar a
zonas slidas la absorcin de calor debida a la propia radiacin. Por ello se hace necesario utilizar una
condicin de flujo de calor sobre el absorbedor, cuyo mdulo se podr calcular como el producto escalar
del vector normal a la superficie del absorbedor por el vector de radiacin solar incidente, a su vez este
valor se multiplicar por el valor intensidad de radiacin. En captulos posteriores se explicar con ms
detalle el modelo de radiacin.
La conclusin al prrafo anterior es que el absorbedor estar sometido en su parte superior a fenmeno
de conveccin libre y a su vez absorber radiacin, como se ha dicho, a travs de una condicin de
contorno de flujo de calor constante.
Esto no es posible por la forma en que Fluent trata a la interfase entre zonas de distinto material, o lo
que llama two-sided-walls. Es decir, cuando lee una malla y se da esta caracterstica automticamente
crea una shadow, o sea, asigna a cada zona, para nuestro caso, absorbedor y aire, una condicin de
contorno diferente. Por que cabe la posibilidad de especificar diferentes condiciones trmicas en cada
zona o mantener acopladas dichas zonas.
- Si las dos zonas son acopladas (Coupled), no se requerir ninguna condicin de frontera trmica
adicional, ya que el solver calcular el flujo de calor directamente de la solucin en las celdas
87

adyacentes. Esta opcin por tanto no es viable, dado que no es posible introducir el flujo de
calor debido a la radiacin y tambin tener acoplado la cara superior del absorbedor al proceso
de conveccin natural que se da en el aire que la rodea.
- Por otra parte se pueden desacoplar las dos caras de la condicin de contorno y elegir
diferentes condiciones trmicas en cada una. Como se ha dicho antes, cada una de ellas estar
asignada al aire y el absorbedor, respectivamente. Si se elige esta opcin, la condicin de
contorno y su shadow asociada, estn totalmente aisladas. Como se puede ver en la figura,
las condiciones de contorno especificadas se asignan sobre la cara de dentro de las paredes:

Figura 3.31: Condiciones trmicas sobre condiciones de contorno desacopladas (Original y Shadow)
Por lo que se llego a la siguiente conclusin: utilizar la opcin Coupled sobre las condiciones de contorno
(Original y Shadow) introduciendo una Generacin interna de calor en la placa que acta de absorbedor
cuyo valor fuera el mismo que la condicin de flujo de calor, es decir, transformar los
, a
,
sabiendo el espesor del absorbedor. Esta es la condicin que usaremos para el modelo completo sin
modelo de carga solar.
Resumiendo, las condiciones de contorno del modelo para un tubo con absorbedor son:
Absorbedor_cara_inferior Wall (adiabtica)
Absorbedor_cara_inferior:001 Wall (coupled)
Absorbedor_cara_inferior:001_shadow Wall (coupled)
Absorbedor_cara_superior Wall (adiabtica/flujo de calor cte)
Absorbedor_frontal Wall (adiabtica)
Absorbedor_lado_I Wall (adiabtica)
Absorbedor_lado_II Wall (adiabtica)
Absorbedor_trasera Wall (adiabtica)
Cara_exterior_tubo Wall (adiabtica)
Cara_exterior_tubo:022 Wall (coupled)
Cara_exterior_tubo:022-shadow Wall (coupled)
Cara_interior_tubo Wall (coupled)
Cara_interior_tubo-shadow Wall (coupled)
Entrada_de_tubera Wall (adiabtica)
Salida_tubera Wall (adiabtica)
Soldadura_frontal Wall (adiabtica)
88

Soldadura_lado_I Wall (adiabtica)
Soldadura_lado_II Wall (adiabtica)
Soldadura_trasera Wall (adiabtica)
Tabla 3.5: Condiciones de contorno para el modelo de tubo con absorbedor
2.3. Resultados obtenidos
2.3.1. Comparacin de Nusselt para tubo simple y tubo con absorbedor
A continuacin se muestra una comparacin entre casos de tubo simple y tubo con absorbedor con la
misma temperatura de entrada, mismo flujo de calor y mismo caudal.

Figura 3.32: Comparacin del nmero de Nusselt de tubo simple y tubo con absorbedor.
La principal diferencia entre el modelo de tubo simple y el que tiene absorbedor es la distribucin del
flujo de calor, que a su vez influye en el perfil de temperaturas de la seccin del conducto. Sin embargo
y a la vista de los resultados que se ven en la Figura 3.32, hace que aunque de manera circunferencial en
la pared del tubo las diferencias en el flujo de calor sean importantes, las recirculaciones provocadas por
los efectos de flotacin en conveccin mixta hacen que prcticamente el nmero de Nusselt sea
ligeramente superior al que se tena con tubo liso, principalmente para flujo desarrollado.
El incremento es debido a la aparicin de zonas en la parte inferior del conducto que ven favorecido el
intercambio de calor gracias a elevados gradientes de temperatura debido a mayores fuerzas de
flotacin. Este aumento en la intensidad de las recirculaciones viene dado tambin por el efecto
provocado por la soldadura que hace que las temperaturas alcanzadas en la parte superior del conducto
sean mayores.
1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Petukhov y Polyakov [1988], Conveccin forzada
Petukhov y Polyakov [1988], Conveccin Mixta, 50_1Q_1000
Petukhov y Polyakov [1988], Conveccin mixta 26_1Q_600
Simulacin 50_1Q_1000, Tubo Simple
Simulacin 50_1Q_1000, Tubo con absorbedor
Simulacin 26_1Q_600, Tubo con absorbedor
Simulacin 26_1Q_600, Tubo Simple
89

Para demostrar lo anterior se muestran a continuacin una representacin grfica de las componentes
de velocidad en direccin X y las componentes en direccin Y, para una seccin a 1.7 metros de la
entrada (donde el nmero de Nusselt ya es diferente para cada caso) para el caso 50_1Q_1000 que se
mostrado en la Figura 3.32.

Figura 3. 36: Componente de la velocidad en X. Tubo simple.
Seccin 1.7m.
Figura 3. 35: Componente de la velocidad en X. Tubo con
absorbedor. Seccin 1.7m.
Figura 3.34: Componente de la velocidad en Y. Tubo simple. Seccin
1.7m.
Figura 3. 33: Componente de la velocidad en Y. Tubo con
absorbedor. Seccin 1.7m.
90

En el tubo simple es cierto que se alcanza una velocidad tanto en direccin X como en direccin Y,
mayor a la que se alcanza en el tubo con absorbedor, aunque, es muy localizada, dndose en la zona
superior del conducto. Sin embargo, en el tubo con absorbedor a pesar de alcanzarse una velocidad
mxima inferior a la que se tiene en la zona superior en el tubo simple, (tanto en direccin X como en
Y), se logran velocidades mayores y ms uniformes en el resto de la seccin.
Esto justifica lo que ya se haba adelantado, y es que en el tubo simple se tiene nmeros de Nusselt muy
uniforme a lo largo de toda la seccin, pero, cuando aadimos la soldadura y el absorbedor, el nmero
de Nusselt alcanzada en las cercanas de la soldadura es menor que la se tena en el tubo simple, aunque
en el resto de pared se tienen velocidades mayores y por lo tanto mejor transmisin de calor, esto hace
que de manera promediada, el nmero de Nusselt se incremente.
2.3.2. Nmero de Nusselt desarrollado

Figura 3.37: Nmero de Nusselt desarrollado para el tubo con absorbedor. Comparacin con la correlacin de
Petukhov y Polyakov.
En la figura anterior se ha representado el nmero de Nusselt para flujo plenamente desarrollado. Para
ello se ha hecho el promedio de este parmetro adimensional en cada seccin para aquellos valores que
tuviesen relacionado un valor de distancia adimensional por encima de 0.01. De esta manera nos
aseguramos de que todo los valores que promediamos estn en la parte de flujo desarrollado.
Como era de esperar observando el grfico de nmero de Nusselt local frente a x*, los valores quedan
por encima de lo que predecan los resultados experimentales. As mismo se ha obtenido una
correlacin para los resultados numricos que nos permita comparar con las simulaciones del captador
completo. Esta es:

y = 0.9026x
0.1759

R = 0.9684
1
10
100
1.00E+05 1.00E+06 1.00E+07
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Nusselt desarrollado
Petukhov y Polyakov , Conveccin
Mixta, Nusselt desarrollado [1988]
Potencial (Nusselt desarrollado)
91

2.3.3. Comparacin de la condicin de contorno flujo de calor frente a generacin
interna

Ya se ha comentado que para utilizar el modelo con la cavidad de aire entre absorbedor y cubierta de
vidrio sin implementar el modelo de carga solar no es posible utilizar la condicin de flujo de calor
constante para no desacoplar el problema, por ello se define la otra opcin de generacin interna de
energa. Los resultados que se obtienen al comparar una y otra condicin son:

Figura 3.38: Comparativa para el nmero de Nusselt con una condicin de contorno de flujo de calor uniforme y
una condicin de generacin interna de calor en el absorbedor, en conveccin forzada. Reynolds 470-900

1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

x*
Nusselt Qcte Conveccin forzada
Nusselt GIcte conveccin forzada
92

Figura 3. 39: Comparativa para el nmero de Nusselt con una condicin de contorno de flujo de calor uniforme y
una condicin de generacin interna de calor en el absorbedor, en conveccin mixta. Reynolds 470-900.
Los resultados para con ambas condiciones de contorno prcticamente son coincidentes, por lo que
para las sucesivas simulaciones donde se modele el colector completo se podr utilizar la Generacin
interna, y as evitar desacoplar el problema, de radiacin y conveccin. La grfica anterior est realizada
para Re 400-900.

2.3.4. Conclusiones
Queda probado que si el flujo de calor se aade sobre la pared del absorbedor en vez de hacerlo en vez
de imponer una condicin de flujo de calor uniforme en la superficie del conducto se incrementa
ligeramente la capacidad de transmisin de calor al flujo, se est hablando del orden de 0.2C de
incremento en la temperatura de salida del fluido. Y por otra parte se ha justificado que se puede
utilizar generacin interna de calor y flujo de calor constante obtenindose resultados muy similares, y
de esta forma poder crear un modelo completo sin problemas de desacoplamiento entre las diferentes
zonas que lo forman.

1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

x*
Nusselt Qcte conveccin mixta
Nusselt GIcte conveccin mixta
93

3. Cavidad natural en la cavidad de aire entre absorbedor y
cubierta de vidrio

Una vez que hemos elegido como simular una carga solar de manera simplificada, vamos a centrarnos
sobre el proceso de conveccin natural en el aire que rodea al absorbedor y al conducto. En el capitulo
anterior se ha introducido este fenmeno de manera terica, exponiendo las ecuaciones que gobiernan
el flujo y los parmetros adimensionales relacionados.
Ahora se va a explicar cmo configurar FLUENT para poder simular este efecto de manera satisfactoria.
FLUENT permite analizar casos con conveccin natural dentro de un dominio cerrado, donde la solucin
depender de la masa dentro del dominio. Pero esta masa no ser conocida hasta que no se conozca la
densidad, por lo que se ofrecen dos posibilidades para modelar este fenmeno:
- A travs de un clculo transitorio. De esta forma, la densidad inicial ser computada teniendo
en cuenta la presin inicial y la temperatura, por lo que la masa inicial es conocida. Conforme la
solucin progresa en el tiempo, la masa se conservar. No se ha elegido esta forma de modelar,
ya que por un lado en nuestro proyecto nicamente nos hemos dedicado a calcular en estado
estacionario, y por otra parte este enfoque funciona bien cuando existen grandes diferencias
de temperaturas, y en el modelo del colector no se dan tales circunstancias.

- La forma elegida es la aproximacin de Boussinesq. Desde este punto de vista, lo que se
especifica una densidad constante, por lo que la masa est claramente especificada. Este
enfoque es nicamente vlido si las diferencias de temperaturas son pequeas.
Esta eleccin est fundamentada en que para muchos problemas de conveccin natural se puede
conseguir una convergencia ms rpida con esta aproximacin que introduciendo la densidad como
funcin de la temperatura. Este modelo trata la densidad como una constante en todas las ecuaciones
excepto en la ecuacin de cantidad de movimiento:

donde
es una densidad constante,
es la temperatura de operacin, y es el coeficiente de

expansin trmica. Esta aproximacin se dice que tiene validez cuando los cambios de densidad son
pequeos y se cumple que:
.
3.1. Simulaciones en 2-dimensiones
Para conocer si la densidad de malla era la adecuada a la simulacin, y justificar la decisin de usar la
aproximacin de Boussinesq, se hicieron una serie de pruebas con una malla 2D, con las dimensiones
del recinto superior de aire (el que queda entre absorbedor y cubierta de vidrio). En este caso se fijaron
unas temperaturas superior e inferior de modo que se obtuviese un nmero de Rayleigh similar al que
se obtiene en los casos simulados en el modelo completo para un tubo.

94

3.1.1. Geometra y mallas utilizadas

Figura 3.40: Malla 2D espacio de aire superior. Pruebas para la conveccin natural. Malla 1.

Las dimensiones como vemos en la geometra anterior para todas las mallas probadas fue de 123 x 20
mm, es decir, el mismo espacio que hay entre absorbedor y cubierta.
En cuanto a la construccin de la malla fue muy simple, se fijaron los nodos en cada uno de los
segmentos que forman el espacio, dejando un espacio menor en las cercanas de las paredes donde se
fijo una condicin de temperatura constante ya que el gradiente de temperaturas sera mayor en esos
puntos (Figura 3.41).

Figura 3.41: Detalle de la malla junto a las paredes. Pruebas para la conveccin natural.
95

Finalmente se construyeron tres mallas con la misma distribucin de nodos que hemos visto en la Figura
3.40, pero decrecientes en densidad, de forma que:
Malla N de Nodos
Malla 1 1701
Malla 2 451
Malla 3 396
Tabla 3.6: N de nodos para mallas 2D en cavidad de aire.
3.1.2. Hiptesis de simulacin y condiciones de contorno

o Solver: Pressure Based Segregated (SIMPLE)
o Esquema de discretizacin: Primer Orden
o Simulaciones estacionarias
o Modelo Viscoso: Laminar
o Ecuacin de la energa: resuelta
o Gravedad: (0,-9.81,0) (m/s
2
)
o Fluido de trabajo: aire

Una vez establecido la direccin y mdulo de la gravedad, para imponer la aproximacin de Boussinesq
tenemos que realizar las siguientes acciones:

o Seleccionar dentro del campo de caractersticas del aire, en el apartado de densidad,
el modelo de Boussineq, e introducir un valor constante. Para este caso como se
trabaja entre 300 y 315 K, se ha introducido la densidad para un valor intermedio,
1.1614 kg/m
3
.
o Especificar la temperatura de operacin

o Y, tambin en el panel de materiales, para el aire, definir el coeficiente de expansin
trmica, que si tratamos al aire como gas ideal este es aproximadamente, .
En este caso no se da, pero para el modelo completo, se puede especificar que materiales y cules no
son los que participan de la aproximacin.

Figura 3.42: Seleccin de la temperatura de operacin y de la gravedad.

96

Para las paredes laterales se ha impuesto la condicin de adiabticas, tal y como ocurre en el modelo
completo. Para las paredes superior e inferior se ha impuesto una condicin de temperatura constante,
315 K para la pared inferior y 300 K para la superior. Se impone esta diferencia de temperaturas para
obtener un nmero de Rayleigh similar al obtenido en el modelo completo. Por supuesto, al poner la
temperatura mayor en el lado inferior supone la aparicin de efecto de conveccin natural, por el
gradiente de densidad junto con la fuerza gravitatoria.
El nmero de Rayleigh definido para este tipo de problema, est basado en la altura de la cavidad, L:

donde es la difusividad trmica del aire (m
2
/s) y es la viscosidad (Pa.s).
Para la diferencia de temperaturas anterior y evaluando las propiedades del fluido a la temperatura
media, el nmero de Rayleigh para las simulaciones realizadas es:

3.1.4. Resultados obtenidos
Para evaluar los resultados obtenidos se ha calculado el nmero de Nusselt, y se ha comparado con la
correlacin presentada en el captulo anterior sobre conveccin libre en cavidades inclinadas (K.G.T.
Hollands, 1976).

Para el nmero de Rayleigh simulado, .
El procedimiento de clculo del nmero de Nusselt es:

donde
, es el flujo de calor por unidad de rea evaluado en la pared ms caliente, L es la longitud

caracterstica, en este caso, la altura de la cavidad (L=20 mm), , es la diferencia de temperatura entre
las paredes superior e inferior y , es la conductividad trmica a la temperatura media del aire.
Los resultados para cada malla son:
Malla Nusselt simulacin Nusselt Correlacin Error (%)
Malla 1 2.594 2.6 0.23
Malla 2 2.53 2.6 2.76
Malla 3 2.463 2.6 5.56
Tabla 3.7: Error cometido en cada malla para simulaciones 2D en cavidad de aire
97

Como vemos en la tabla anterior las simulaciones parecen acercarse de forma muy notable a los
trabajos realizados experimentalmente, siendo para la malla ms fina para la cual se tiene un menor
error. Adems se ha conseguido una convergencia muy rpida, haciendo falta unas 1000 iteraciones
para llegar a valores residuales de 1x10
-10
en la ecuacin de energa y 1x10
-5
para la ecuacin de
continuidad.
Adems de evaluar el Nusselt, tambin se ha visualizado el campo de temperaturas y las componentes
del campo de velocidades, observndose la generacin de recirculaciones en forma de celda, con una
altura entorno a la dimensin caracterstica, L como predecan los trabajos de Hollands (celdas de
Bernard).

Figura 3.43: Campo de temperaturas y velocidades. Malla 1.Conveccin natural. Ra=11000

Figura 3.44: Campo de temperaturas y velocidades. Malla 2.Conveccin natural. Ra=11000.
98

Figura 3.45: Campo de temperaturas y velocidades. Malla 3. Conveccin natural. Ra=11000.

Como podemos observar, sobre todo en la Figura 3.43, las celdas formadas son 6, esto coincide con los
resultados experimentales que afirman, que se desarrollan recirculaciones con una dimensin
caractersticas de L (L=20 mm), por lo que si se hace un clculo rpido dividiendo el ancho de la cavidad
123 mm por L, efectivamente se obtiene aproximadamente un valor de 6.
3.1.5. Conclusiones
En primer lugar se ha comprobado que se puede utilizar la aproximacin de Boussinesq en nuestro
problema para resolver este fenmeno obteniendo una convergencia rpida. Utilizar una malla de
entorno a los 1000 nodos de forma transversal al modelo puede ser suficiente. Y por ltimo decir que
como en el modelo completo se tendrn unas condiciones de temperatura uniformes en el absorbedor,
y muy probablemente las celdas de recirculacin que se obtengan no sigan una distribucin tan
uniforme como en este caso.

99

4. Modelo completo sin radiacin para un slo conducto

Este modelo consta de todas las partes por las que est formado el colector solar, es decir, tanto de la
cubierta de cristal, el aire sobre el absorbedor y debajo del mismo, y por supuesto del absorbedor, la
soldadura y el conducto. Pero se ha simplificado y nicamente se modela la parte correspondiente a
slo tubo.
En este caso no se ha implementado un modelo de radiacin, sin embargo, para tener en cuenta la
energa que procede del Sol se ha introducido una condicin de generacin interna de calor en el
absorbedor, para simular la carga solar. Esto nos permite analizar la transmisin de calor a travs del
absorbedor por la soldadura, as como, la conveccin natural que se da en el recinto limitado por el
absorbedor y la cubierta de vidrio. Adems, para tener en cuenta las prdidas por efecto convectivo
sobre la cubierta se ha impuesto una condicin de contorno convectiva.
El objetivo del modelo fue poder obtener resultados sobre la transmisin de calor a travs de la
soldadura sin necesidad de implementar un modelo de carga solar, elegir la mejor opcin para simular el
comportamiento del aire en el interior del colector y observar el campo de temperaturas en el
absorbedor. Se trata por tanto de un paso intermedio en el proceso de modelado hasta llegar a la
solucin definitiva.

La geometra que forma el modelo viene descrita en la figura siguiente:

Figura 3.46: Geometra del modelo completo para un tubo.

100

Adems, rodeando la geometra anterior, como se ve en la Figura 3.46, se ha definido una capa de aire
tal y como existe en un colector real. En cuanto a sus dimensiones son, en la parte superior 20 mm por
encima del absorbedor, y en la parte inferior la capa tiene un espesor de 15.5 mm contando desde el
centro del conducto (origen de coordenadas).
- Zonas que forman el modelo:
Zona Slido/Fluido
Absorbedor Aluminio
Aire_superior Aire
Aire_inferior Aire
Cristal Vidrio
Fluido Agua
Soldadura Cobre
Tubo Cobre
Tabla 3.8: Volumenes incluidos en el modelo completo
Partiendo de la malla que se construyen en apartado 2, seguimos con la construccin de la misma de la
siguiente forma:
Por ltimo falta por mallar la zona del aire, tanto del que se encuentra en la parte inferior del tubo
como el que se encuentra en la parte superior del absorbedor. Hay que destacar que la zona inferior del
modelo y las paredes laterales son adiabticas, es decir, se ha supuesto que el aislamiento es perfecto.
Esto hace que verdaderamente el aire de la parte inferior no tenga demasiada importancia desde el
punto de vista energtico.
Sin embargo, la zona superior del aire, va a estar sometida a una diferencia de temperaturas entre la
frontera superior perteneciente al vidrio y la frontera inferior donde se encuentra el absorbedor,
estando este ltimo ms caliente que el vidrio. Esto, junto con la gravedad hace que el aire ms caliente
junto al absorbedor tienda a ascender por tener una menor densidad mientras que el aire en la parte
superior ,ms frio, descienda, es decir, se crea un flujo de conveccin natural que va a tener bastante
importancia en la eficiencia del captador. Por esto, la malla que se genere en esta zona debe ser
adecuada para describir los fenmenos convectivos con suficiente precisin.

Para poder conocer cul podra ser el nmero adecuado de elementos que conformaran la malla se opt
por realizar unas simulaciones con un modelo ms simple en 2 dimensiones cuyos resultados se
exponen en el punto siguiente.
101

Finalmente, con los resultados que se obtuvieron se realiz una malla de esa zona que consta de 1806
nodos con la siguiente distribucin:
Figura 3.47: Malla para el aire entre absorbedor y cubierta de vidrio. Modelo completo 1 solo tubo (sin radiacin)
Se puede apreciar en la figura anterior como tanto en las cercanas del absorbedor como en el vidrio se
ha intentado que los nodos estuviesen ms cercanos unos de otros, ya que el gradiente de temperatura
en esa zona ser considerablemente superior que el de la zona central.
El punto de partida eran los nodos que ya estaban sobre el segmento superior del absorbedor lo que
supona crear una malla demasiado fina si se optaba porque en el vidrio estuviesen igualmente
distribuidos, es decir, habra una zona central con muchos elementos. Sin embargo, si en la cubierta se
colocan pocos nodos entonces las celdas sera muy oblicuas y los ngulos muy cerrados, o sea, una malla
muy irregular. Por lo tanto se tom una solucin de compromiso de manera que se tuviera una malla lo
suficientemente regular pero que al mismo tiempo no fuese demasiado fina, ya que hay que recordar
que cada una las secciones como la que se ha mostrado se repiten a largo del conducto 183 veces, es
decir, cada centmetro.
Por ltimo falta por mallar la zona inferior del aire, teniendo en cuenta que debe conectar con todos los
nodos situados en la pared exterior del tubo y la cara inferior el absorbedor, lo que supone una cantidad
enorme de elementos si se malla de forma regular. Por ello, y dado que esta parte no afectar, como ya
se ha comentado, en exceso a la configuracin del problema, se ha optado por una malla irregular
intentando que el coste computacional no se disparase.

102

Figura 3.48: Malla modelo completo 1 solo tubo (sin radiacin). Seccin transversal.

Una vez creada la seccin transversal de la geometra y realizada la malla, solo queda extruirla para
conseguir crear el modelo completo. Para ello se especific la direccin del vector normal hacia donde
se quera expandir, la magnitud de dicho vector normal y se confirm la opcin de arrastrar la malla al
crear el volumen final a travs de la opcin Sweep faces que da GAMBIT:

Figura 3.49: Vector y magnitud para definir el volumen completo.

103

Finalmente se obtiene el modelo completo, repitiendo la malla transversal cada centmetro, es decir,
cada elemento tiene un centmetro segn el eje Z.
Grid Size
0 990030 3074869 1093296 8

Memory Usage
Number Used 990030 3074869 1093296 48 0
Mbytes Used 323 325 58 0 0

4.2. Hiptesis de simulacin

- Gravedad: (0,-9.81,0) (m/s
2
)

- Condiciones de operacin:

Figura 3.50: Condiciones de operacin para el modelo completo para 1 solo tubo (sin radiacin).
Tal y como se ha justificado en el punto anterior se ha utilizado la aproximacin de Boussinesq para
simular la conveccin natural en el aire, sin especificar la densidad de este ltimo en funcin de la
temperatura, por ello es necesario definir al temperatura de operacin
en este panel.
Como ocurra en el caso del tubo simple, las propiedades del agua van en funcin de la temperatura, es
decir, para simular la conveccin mixta en el interior del tubo no se utiliza la aproximacin de
104

Boussinesq debido a que tenemos diferencias de temperatura entre tubo y fluido elevadas adems de
nmeros de Rayleigh del orden de 1x10
6
, que hacen que esta aproximacin no converja de manera
satisfactoria. De lo anterior se deduce que efectivamente FLUENT permite utilizar la aproximacin para
unos materiales y para otros no.
Para el agua, adems de definir sus propiedades en funcin de la temperatura, es posible especificar
otro parmetro llamado, densidad de operacin
. Por defecto, FLUENT calcula este parmetro como

promediado en todas las celdas del fluido en cuestin, pero al especificarlo se puede mejorar la
convergencia.
4.2.1. Materiales
Material Tipo Propiedades
Agua Fluido Funcin de la temperatura (Ver 1.2.4) (Anexo I)
Aire Fluido Constantes excepto densidad (Boussinesq)
Cobre Slido Constantes (Anexo I)
Aluminio Slido Constantes (Anexo I)
Vidrio Slido Constantes (Anexo I)
Tabla 3.9: Materiales en el modelo completo
Solo remarcar que para el caso del aire es necesario especificar el coeficiente de expansin trmica .
Todas las propiedades se encuentran definidas en el Anexo I.

- Entrada_de_fluido
Al igual que en el caso del tubo simple, esta condicin es de tipo, velocity-intlet, y en ella especificamos
la velocidad de entrada y la temperatura del agua en el interior del dominio.
Las simulaciones realizadas para este modelo realizan un barrido de temperaturas de entrada desde los
20 hasta los 70C, para tres irradiancias diferentes (600,800 y 1000 W/m
2
). Todos estos casos , al igual,
que en el tubo simple se han ensayado para un gasto msico por unidad de rea de absobedor, esto es,
0.02 kg/s m
2
.
Por lo tanto segn la temperatura de entrada calculamos la densidad, y en funcin de la misma el
caudal para obtener finalmente la velocidad conociendo el dimetro del conducto.
- Salida_de_fluido
Aqu se da la misma situacin que en el caso con tubo simple. Es una condicin del tipo pressure-
outlet. Ver 1.2.3
- Fronteras externas del modelo
Todas aquellas superficies del modelo que son la frontera del dominio son superficies de tipo wall,
todas ellas se han considerado adiabticas, es decir, tanto la superficie inferior como los laterales (que
en la realidad formaran parte del tubo siguiente) excepto la superficie superior, la que compondra la
cubierta de vidrio.
105

Figura 3.51: Frontera superior del modelo. Condicin de contorno de tipo convectiva.

La frontera superior del captador se definido como de tipo convectiva, imponindose una
temperatura en el exterior del dominio y un coeficiente convectivo. Por supuesto que hubiese sido ms
realista crear un modelo mayor donde el colector permaneciera rodeado de aire pudiendo variar la
velocidad del viento, pero sin duda el coste computacional se habra disparado.
Se opt por utilizar una correlacin ms ampliamente utilizada para el coeficiente convectivo en
captadores solares en funcin de la velocidad del viento. La que se eligi fue:

En todos los casos se supuso una temperatura exterior de 25 C y una velocidad del viento de 4 m/s.
Una vez que introducimos los datos, FLUENT calcula el flujo de calor a travs de la pared como:

Donde
son los valores que nosotros especificamos y
es la temperatura de pared de la
cubierta.
Adems de estos parmetros trmicos en la Figura 3.51, tambin se definen el grosor y material que
forman el vidrio, ya que aunque por defecto una frontera no tiene espesor, es posible crear una fina
capa de material, de manera que FLUENT resuelve una ecuacin de conduccin en una dimensin para
establecer cul es la resistencia trmica ofrecida por la pared.
106

Figura 3.52: Fina capa de material creada sobre la misma condicin de contorno.
La resistencia trmica en la pared es x/k, donde k es la conductividad de la pared y x es el espesor de
la misma. La condicin de contorno que se ha definido se especifica sobre la superficie fuera del
dominio fluido, es decir, sobre la inner surface.
- Fronteras internas en el modelo
Todas las superficies son de tipo wall, de manera que como son fronteras que se encuentran entre dos
regiones bien slido/slido o fluido/solido, FLUENT de forma automtica a leer la el archivo .msh,
crear una superficie shadow, es decir, habr una condicin de contorno para una regin y su
correspondiente shadow para la otra regin con quien comparte la interfase.
La relacin entre una frontera y su correspondiente shadow, se puede hacer de dos formas, como ya
se ha comentado en el apartado Error! No se encuentra el origen de la referencia., finalmente se
ecidi que las superficies permaneciesen acopladas (coupled), de manera que se calcule el flujo de
calor de forma directa a partir de la solucin de las celdas adyacentes.
Por lo tanto para tener en cuenta la irradiancia que finalmente absorbe el absorbedor, no se impone
una condicin de flujo de calor constante, sino que se define un trmino de generacin interna de calor
en la zona absorbedor de la siguiente forma:
Si por ejemplo el flujo de calor que finalmente absorbe la lmina de aluminio es de
W/m
2
,
entonces el trmino de generacin interna de calor (W/m
3
) se calcula como:

donde
es la superficie de captacin correspondiente a un tubo, es decir, 0.123mx1.83m, y
es
el espesor del absorbedor (0.5 mm).
107

Por lo tanto de manera aproximada es posible calcular el flujo de calor absorbido por el captador
realizando el producto escalar del vector normal de la superficie absorbedora por el vector que indica la
orientacin del haz de radiacin directa (sun direction vector), este vector se puede obtener utilizando
la herramienta Solar Calculator (que explicaremos en el modelo con radiacin) . El resultado de dicho
producto lo multiplicamos por la radiacin directa normal y le sumamos el valor de la radiacin difusa.
A este valor lo multiplicamos por el producto absortancia-transmitancia y tendremos un valor
aproximado de la potencia por unidad de superficie que finalmente absorbe el colector solar.

Figura 3.53: Fuente de generacin interna de calor para un flujo de 1000 W/m
2
.

A continuacin se presenta una tabla resumen con las condiciones de contorno establecidas para las
fronteras del modelo:

Absorbedor_cara_inferior Wall (acoplada)
Absorbedor_cara_inferior-shadow Wall (acoplada)
Absorbedor_cara_inferior:001 Wall (acoplada)
Absorbedor_cara_inferior:001_shadow Wall (acoplada)
Absorbedor_cara_superior Wall (acoplada)
Absorbedor_cara_superior-shadow Wall (acoplada)
Absorbedor_frontal Wall (adiabtica)
Absorbedor_lado_I Wall (adiabtica)
Absorbedor_lado_II Wall (adiabtica)
Absorbedor_trasera Wall (adiabtica)
Air_inf_front Wall (adiabtica)
Air_inf_ladoI Wall (adiabtica)
Air_inf_ladoII Wall (adiabtica)
Air_inf_tras Wall (adiabtica)
Air_sup_front Wall (adiabtica)
Air_sup_ladoI Wall (adiabtica)
Air_sup_ladoII Wall (adiabtica)
Air_sup_tras Wall (adiabtica)
108

Cara_exterior_tubo Wall (acoplada)
Cara_exterior_tubo-shadow Wall (acoplada)
Cara_exterior_tubo:34 Wall (acoplada)
Cara_exterior_tubo:034-shadow Wall (acoplada)
Cara_interior_tubo Wall (acoplada)
Cara_interior_tubo-shadow Wall (acoplada)
Cristal Wall(conveccin)
Entrada_de_tubera Wall (adiabtica)
Salida_tubera Wall (adiabtica)
Soldadura_frontal Wall (adiabtica)
Soldadura_lado_I Wall (acoplada)
Soldadura_lado_I-shadow Wall (acoplada)
Soldadura_lado_II Wall (acoplada
Soldadura_lado_II-shadow Wall (acoplada)
Soldadura_trasera Wall (adiabtica)
Suelo Wall (adiabtica)
Tabla 3.10: Condiciones de contorno para el modelo completo con un solo tubo

Las simulaciones se han realizado para diferentes temperaturas de entrada, tres flujos de calor
diferentes que absorbe el captador (600,800 y 1000 W/m
2
), y manteniendo constante el gasto msico
por unidad de rea de captacin en 0.02 kg/s, as como la velocidad del viento de 4m/s.
Caso Temperatura de
entrada (C)
Generacin interna
(W/m
3
)
Velocidad de entrada
(m/s)
600_Te_20 20 1,200,000.0 0.117204696
600_Te_30 30 1,200,000.0 0.117497862
600_Te_40 40 1,200,000.0 0.117891623
600_Te_50 50 1,200,000.0 0.118380874
600_Te_60 60 1,200,000.0 0.118960782
600_Te_70 70 1,200,000.0 0.119626691
800_Te_20 20 1,600,000.0 0.117204696
800_Te_30 30 1,600,000.0 0.117497862
800_Te_40 40 1,600,000.0 0.117891623
800_Te_50 50 1,600,000.0 0.118380874
800_Te_60 60 1,600,000.0 0.118960782
800_Te_70 70 1,600,000.0 0.119626691
1000_Te_20 20 2,000,000.0 0.117204696
1000_Te_30 30 2,000,000.0 0.117497862
1000_Te_40 40 2,000,000.0 0.117891623
1000_Te_50 50 2,000,000.0 0.118380874
1000_Te_60 60 2,000,000.0 0.118960782
1000_Te_70 70 2,000,000.0 0.119626691
Tabla 3.11: Simulaciones realizadas para el modelo completo con un solo tubo

109

4.3. Anlisis de resultados
4.3.1. Anlisis de la transmisin de calor
Para el modelo que estamos estudiando se va analizar el proceso de conveccin mixta que se da en el
interior del tubo, teniendo en cuenta que las condiciones bajo las cuales se estn realizando las
simulaciones distan mucho de las que se tenan cuando se estudiaba la transmisin de calor en tubo
simple horizontal con flujo de calor constante. Hay que resaltar que ahora, el flujo de calor se aplica en
el absorbedor y se transmite por conduccin por las paredes del conducto a travs de la soldadura, con
lo que se pierde la uniformidad, concentrndose, como se ver a continuacin, sobre la parte ms
cercana a este punto, acentundose adems por los efectos de flotacin que hacen que el fluido
caliente tienda a ascender. Adems del estudio del nmero de Nusselt, tambin se va a analizar otro
parmetro adimensional muy influyente dentro de los fenmenos de conveccin natural, el nmero de
Rayleigh, y su relacin con la longitud de desarrollo de la capa lmite trmica.
El procesado de resultados ser el que se utiliz en el caso del tubo simple y el tubo con absorbedor,
dados los buenos resultados que se obtuvieron al comparar los resultados numricos obtenidos con el
trabajo experimental de Petukhov y Polyakov [1988]. Tambin compararemos con la correlacin que se
obtuvo para el Nusselt desarrollado en el tubo con absorbedor.
Adems de comprobar la transmisin de calor en el flujo por el interior del conducto, tambin se
visualizar la distribucin del campo de temperaturas sobre el absorbedor y se analizar la conveccin
natural en la cavidad superior del captador.

4.3.1.1. Nmero de Nusselt en el flujo interno
Como se ha comentado, en este modelo, el flujo de calor ya no es constante en la pared del conducto,
segn la configuracin del colector, el flujo debe pasar a travs de la soldadura, y posteriormente se
distribuye por la pared del tubo. Si tenemos en cuenta que el interior del conducto el fluido ms caliente
por los efectos de la flotacin tiende a situarse en la parte ms cercana a la soldadura, la mxima
diferencia de temperatura se dar en la parte inferior del conducto, y por lo tanto el flujo de calor ser
mayor en esa zona (Figura 3.54: Flujo de calor sobre la superficie interior del conducto. Radiacin
absorbida 600 W/m2.Figura 3.54).
Si representamos la temperatura de fluido junto con la temperatura de la pared podemos justificar lo
que se ha descrito en el prrafo anterior (Figura 3.55). Vemos como efectivamente es en la parte inferior
donde el gradiente de temperaturas es mayor. Por lo tanto los efectos de la conveccin mixta actan
como barrera a la transmisin de calor cerca de la soldadura, pero al mismo tiempo permite que en el
resto del conducto se mantenga un elevado gradiente.
110

Figura 3.54: Flujo de calor sobre la superficie interior del conducto. Radiacin absorbida 600 W/m2.

Figura 3.55: Campo de temperaturas en fluido y tubera. Seccin 1.1 metros.
Adems de la variacin circunferencial del flujo de calor, tambin cambia a lo largo del eje Z. Si
representamos esta variable para algunas secciones a lo largo del conducto podemos observar que a
medida que nos vamos alejando de la seccin de entrada, la variacin es significativa para pequeas
variaciones de la coordenada Z, sin embargo, conforme se va desarrollando la capa lmite trmica esta
variacin se atena.

111

Figura 3.56: Variacin del flujo de calor en la coordenada Y. Caso 600_Te_50

En la Figura 3.56, se comprueba que a 10 cm de la entrada el flujo de calor es prcticamente uniforme a
lo largo de la seccin, todava no se ha desarrollado un perfil definido de temperaturas y velocidades, y
la capa lmite es muy incipiente, por lo que el gradiente de temperaturas en uniforme a lo largo de toda
la seccin. En la seccin a 20 cm la variacin es importante, se puede ver como la diferencia entre la
seccin anterior y sta es significativa, sin embargo las siguientes que se han representado mantienen
valores muy prximos aunque la variacin vaya creciendo de manera muy progresiva.
Por lo tanto si representamos el promedio de flujo de calor en direccin longitudinal se observa como al
principio se tiene un pico, (no hay capa lmite y la transmisin de calor es muy buena) pero cuando
tenemos flujo plenamente desarrollado, prcticamente se mantiene constante en cada seccin XY. Tal y
como hacemos el procesado de Nusselt, donde se realiza un promediado del valor de flujo de calor que
se obtiene en cada nodo de la pared interna del conducto en cada seccin, aunque realmente el flujo no
sea uniforme a efectos de comparar con correlaciones experimentales con flujo de calor uniforme, s
que lo es. Ver Figura 3.57.
Por lo tanto al igual, que se haca con el tubo simple y el tubo con absorbedor, vamos a obtener el
nmero de Nusselt en funcin de la distancia adimensional para cada caso simulado (Figura 3.58). Los
resultados obtenidos son muy similares a los que se obtenan en las simulaciones con tubo simple,
existiendo una regin de entrada, donde se superan los valores que predicen los resultados
experimentales (zona de desarrollo simultaneo de la capa lmite trmica e hidrodinmica), una segunda
parte con transmisin de calor por conveccin forzada donde se obtienen valores muy prximos a la
correlacin (en esta zona ya se ha desarrollado la capa lmite hidrodinmica) y por ltima en la zona de
flujo plenamente de desarrollado (conveccin mixta) donde para todo Rayleigh siempre los resultados
numricos superan a los de la correlacin.

112

Esto es debido a que el flujo de calor aplicado ahora al fluido es asimtrico acimutalmente y esto se
traduce en mayores recirculaciones y por tanto un incremento en el nmero de Nusselt frente a cuando
se aplica un flujo de calor constante en la pared del tubo directamente.

Figura 3.57: Variacin del flujo de calor en funcin de la coordenada Z. Caso 600_Te_50.

2000
2200
2400
2600
2800
3000
3200
3400
3600
3800
0 0.5 1 1.5 2
F
l
u
j
o

d
e

C
a
l
o
r
,

W
/
m
2

Distancia sobre el eje Z (m)
Flujo de calor (W/m2)
1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Petukhov y Polyakov [1988], Conveccin
forzada
Mixta, Rayleigh = 1.87E+05
Mixta, Rayleigh=1.06E+06
Caso 600_Te_20, Rayleigh=1.87E+05
Caso 1000_Te_60, Rayleigh=1.06E+06
Figura 3.58: Nmero de Nusselt frente a la distancia adimensional. Simulaciones para modelo completo con 1 tubo.
113

Figura 3.59: Resultados experimentales de transmisin de calor en conveccin forzada y mixta en tubo liso de la
Tesis de A.Garca [2006].
Los resultados experimentales sobre tubo liso que obtiene A. Garca [2006] para conveccin mixta
tambin se encuentran por encima, como se ve en la figura anterior, de lo que dice la correlacin de
Petukhov y Polyakov [1988]. Se ha obtenido numricamente unos resultados muy parecidos a los
obtenidos experimentalmente, especialmente en cuanto a tendencias (aunque en estos casos
experimentales el flujo de calor tambin se aplica directamente en la pared del tubo)
Para concluir el estudio de la trasmisin de calor al flujo interno, vamos a comparar con los resultados
que han obtenido para el caso de tubo con absorbedor y los de estas simulaciones:

Figura 3. 60: Resultados numricos de transmisin de calor en conveccin mixta. Nmero de Nusselt desarrollado
en conveccin mixta para modelo completo. Comparacin con los resultados numricos en tubo+absobedor
1
10
100
100000 1000000 10000000
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Resultados Numricos Modelo completo 1 tubo
Petukhov y Polyakov [1988]
Resultados tubo+absorbedor
114

Se comprueba como resultados para el modelo completo y el tubo con absorbedor son prcticamente
coincidentes, lo que significa que los procesos de conveccin natural y las prdidas que se producen a
travs de la cubierta de vidrio no tienen influencia en los efectos de conveccin mixta en el interior del
conducto. Adems, se sigue observando como al modelar el absorbedor, que bsicamente acta como
una aleta, la transmisin de calor mejora.
4.3.1.2. Conveccin natural en el aire entre absorbedor y cubierta.
Al estar el absorbedor a ms temperatura que la cubierta de vidrio, se crea un gradiente de densidades
haciendo que el flujo ms caliente tienda a ascender por conveccin natural. Anteriormente se expuso
un pequeo modelo en 2D para explicar este fenmeno, en l las condiciones de la temperatura
superior e inferior eran constantes, sin embargo, en un colector real, la temperatura del absorbedor
deja de ser uniforme por lo que los resultados que se van obtener van a ser muy diferentes. En cambio,
aunque la distribucin de temperaturas y velocidades sea distinta sigue siendo vigente la aproximacin
de Boussinesq para tener en cuenta este fenmeno.
En los casos que analizbamos en el punto 3.3.5, se observaba como se formaban pequeas celdas de
recirculaciones de aire, con una dimensin del orden de la distancia entre las paredes isotrmicas.
Dichas celdas son comunes cuando el nmero de Ra est por encima del crtico, que para este tipo de
geometra que se sita en torno a
, aunque si el nmero de Rayleigh aumenta mucho el

patrn de flujo es ms complejo e incluso puede convertirse en transitorio. En los casos que se analizan
con el presente modelo, el nmero de Rayleigh se encuentra entre Ra=5000-19000, por lo que podemos
considerar que el flujo es estacionario. El nmero de Raylegih se ha evaluado para la temperatura media
de la cubierta y el absorbedor, es decir:

donde
, son la temperatura media obtenida mediante el area weighted average para la

superficie superior del absorbedor y la cubierta de vidrio, respectivamente. L, es la distancia entre las
paredes anteriores (L=20 mm) y y son la viscosidad y difusividad trmica del aire evaluadas a la
temperatura media
.
Aunque nos encontramos dentro de un intervalo de nmero de Rayleigh moderado, como la
temperatura del absorbedor no es uniforme, en cada seccin el campo de temperaturas de la cavidad
entre absorbedor y vidrio va variando al igual que lo hace la distribucin de recirculaciones que ya no
sigue el patrn establecido del modelo 2D (Figura 3.61).
En la Figura 3.62, se representa el campo de temperaturas para una simulacin con una irradiancia
absorbida de 600 W/m
2
, se distingue claramente como a medida que nos desplazamos hacia la salida
del conducto cada vez al temperatura del aire aumenta, aunque el patrn de flujo es muy similar. Sin
embargo, cuando se aumenta la temperatura de entrada del agua (Figura 3. 63), el patrn vara de
forma muy acusada de seccin a seccin, lo que indica que para nmero de Rayleigh ms elevados el
flujo es ms complejo, y posiblemente se necesite realizar simulaciones transitorias.
Para comprobar, al menos, si los criterios de convergencia han sido los adecuados, se ha comparado con
la correlacin presentada en el captulo anterior. La correlacin es vlida para cavidades donde la
dimensin L (altura) es mucho menor que la otras dos dimensiones y la temperatura inferior superior e
inferior son uniformes. La primera condicin s que se cumple en nuestro modelo, sin embargo, la
segunda no, por lo es previsible que las diferencias sean evidentes.

115

Figura 3.61: Streamlines en la cavidad de aire.

Figura 3.62: Campo de temperaturas para las seccin 0.2 ,0.9, 1.2 y 1.6 m. Caso 600_Te_20.
116

El nmero de Nusselt lo obtenemos como:

donde
es el flujo por unidad de rea desde la pared a alta temperatura.

Figura 3.64: Nmero de Nusselt frente al nmero de Rayleigh para la conveccin natural en la cavidad superior de
aire. Comparacin con las predicciones de Hollands [1984].

1
10
1000 10000 100000
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Nusselt Numrico
Hollands [1984]
Figura 3. 63: Campo de temperaturas para las seccin 0.2 ,0.9, 1.2 y 1.6 m. Caso 1000_Te_50.
117

A pesar de que las condiciones con las que estamos comparando la correlacin no son las mismas con
las que se efectuaron los experimentos, podemos ver en la figura anterior, como el error no es
demasiado elevado, incluso, para aquellas simulaciones donde el gradiente de temperaturas en el
absorbedor es menor, los resultados se acercan mucho ms a las predicciones experimentales. Adems
comparando con el modelo validado de colector solar presentado por Selmi et al [2007], se pueden
establecer algunas similitudes, aunque al geometra del modelo sea diferente (en nuestro caso el tubo
est debajo del absorbedor).

Figura 3.65: Distribucin de temperaturas en la seccin de salida del colector para un flujo de 1.36 g/s. Validation
of CFD simulation for flat plate solar energy collector, Selmi et al [2007].
4.3.1.3. Campo de temperaturas en el absorbedor
A continuacin se va a representar el campo de temperaturas en la pared superior del absorbedor.
Compararemos los resultados obtenidos para el caso de 600 W/m
2
, para las diferentes temperaturas de
entrada, con el objeto de que posteriormente comparar con el modelo de captador completo:

118

Como se aprecia en la figura anterior a medida que aumentamos la temperatura de entrada,
lgicamente se incrementa la temperatura del absorbedor. En cuanto al gradiente transversal, como era
de esperar, en la parte central del absorbedor, debido a la transmisin de calor hacia el fluido, se tiene
una menor temperatura, a medida que nos desplazamos hacia los extremos la temperatura va
aumentando. En el gradiente longitudinal se aprecia un aumento de la temperatura, es decir, al
principio el flujo de calor es mayor, sin embargo, conforme el fluido se va calentado, la diferencia de
temperaturas es cada vez menor. Si comparamos los diferentes casos mostrados, vemos que la
diferencia est en el rango pero no en la distribucin de temperaturas, obtenindose un gradiente muy
parecido en todos ellos.

Figura 3.66: Campo de temperaturas para una generacin interna de 600 W/m
2
. Casos para temperatura de entrada a 20 , 30, 40, 50, 60
y 70C
119

4.3.2. Eficiencia del captador
A continuacin se presenta la curva de eficiencia del colector, que aunque para un solo tubo, permitir
saber al comparar con los resultados para el captador completo si se puede reproducir la curva de
eficiencia un coste computacional mucho menor. La eficiencia trmica de informacin global acerca de
la energa que es capaz de aprovechar del total de energa que recibe del Sol. Se define como:

donde
viene dado por la expresin:
, y son la densidad
y el calor especfico respectivamente y se evalan a la temperatura media del fluido
. G
es la irradiancia globar interceptada y
es el rea del absorbedor.

Calcular el
, no supone ningn problema pues conocemos perfectamente cul es el gasto msico y

la temperatura de entrada y salida. Sin embargo, en este modelo no conocemos G, ya que directamente
se est introduciendo un valor de generacin interna sobre el absorbedor, es decir, lo que estamos
imponiendo es la energa que absorbe directamente el absorbedor, y no la que recibe del Sol.
Si recordamos el modelo analtico presentado en el captulo anterior, la potencia absorbida se defina
como:

Por lo que para poder comparar la solucin, se han tomado la transmisividad y la absortividad del banco
de ensayos de Alfonso XII, 0.93 y 0.95 respectivamente. Por lo que dividiendo el flujo de calor al que
equivale la generacin de interna de calor por este producto absortancia-transmitancia podemos un
obtener un valor ms o menos comparable.
Aunque est claro que en este modelo se cuentan con numerosas simplificaciones (se define la potencia
absorbida, no hay prdidas por conduccin, solo existe 1 tubo ) vamos a comparar con los resultados
que se obtuvieron en el PFC de Mara Jos Martnez sobre curvas de eficiencia. La curva para un colector
para un gasto de 144 kg/h (el mismo que las simulaciones) en tubo liso es:

donde
es la temperatura reducida:

donde
es la temperatura ambiente, que se ha fijado en todas las simulaciones en 25C,
es la
temperatura media del fluido
.
Adems de comparar con los resultados en un colector real, para asegurarnos de que los valores de
eficiencia pueden ser aceptables si tenemos en cuenta las simplificaciones realizadas, se ha comparado
con un modelo analtico algo rudimentario, (aunque con la ecuaciones que se definieron en el captulo
anterior) que se ha modificado para adaptarlo a las condiciones de la simulacin.
Estas modificaciones han consistido en eliminar las prdidas por conduccin a travs de los laterales y la
parte inferior, y nicamente calcular el coeficiente de transmisin para la parte superior. Adems, se ha
introducido directamente la potencia absorbida S, en vez de dejar que calcule G. Por otra parte, en vez
utilizar la correlacin de Churchill y Ozoe para el clculo del Nusselt, se ha introducido en el algoritmo de
clculo la correlacin que se ha obtenido para las condiciones experimentales.
120

Figura 3.67: Eficiencia del captador para modelo completo de1 tubo sin radiacin. Comparacin con banco de
ensayos y modelo analtico adaptado.
Las diferencias con los resultado experimental son muy grandes como ya se prevea, debido a que no se
han tenido en cuenta por ejemplo las prdidas por conveccin por lo laterales y por la parte inferior, sin
embargo, se obtiene una buena correspondencia con el modelo analtico adaptado. Podemos concluir
que aunque nos encontremos todava en una situacin muy ideal, se puede decir que seguimos la
tendencia de los estudios analticos con unas simulaciones muy simplificadas.

0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
-0.01 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
E
f
i
c
i
e
n
c
i
a

d
e
l

c
a
p
t
a
d
o
r

T*
Simulacin
Modelo analtico
Captador experimental
121

5. Modelo completo con carga solar para un solo conducto

En las simulaciones que se describen en este apartado se utiliza la misma malla que se usaba con el
modelo completo sin radiacin pero en este caso s que se ha implementado un modelo de carga solar.
Al simular con un modelo de carga solar, ya no hace falta aplicar la condicin de generacin interna
sobre el absorbedor, sino que ser la propia radiacin exterior la que actuando sobre el modelo haga
que el fluido calo portador llegue a calentarse.
En este apartado nos centraremos en describir sobre todo, como se ha implementado el modelo de
carga solar Solar Ray Tracing, as como de la descripcin terica de los modelos de radiacin
disponible en FLUENT. Tambin se expondr como se ha utilizado la herramienta Solar Calculator
para obtener los parmetros de radiacin necesarios en la simulacin. El objetivo final es comprender
como funciona el modelo de carga solar y validarla para un en un modelo simplificado de menor coste
computalcional.
Las caractersticas geomtricas y el mallado son iguales a los modelos descritos en el apartado anterior,
por lo que nos centraremos directamente en el problema de radiacin
5.1. Introduccin a los modelos de radiacin solar
5.1.1. Modelos de Radiacin Implementados
El software de CFD Fluent tiene implementados 5 modelos diferentes de radiacin, que permiten incluir
est fenmeno de transmisin de calor con o sin un medio partcipante (en el que se incluya un trmino
en la ecuacin de la energa). Dichos modelos son:
- Discrete Trasnfer Radiation Model (DTRM)
- P-1 Radiation Model
- Rosseland Radiation Model
- Surface-to Surface (S2S) Radiation Model
- Discrete Ordinates (DO) Radiation Model
Adems de estos modelos de radiacin Fluent incluye un modelo de carga solar, que va a ser uno de los
puntos importantes dentro de nuestras simulaciones.
5.1.2. Modelo DO
Despus de la revisin de todos ellos, bsicamente quedaron 2 opciones modelo DO o modelo S2S
complementados con un modelo de carga solar. En primer lugar se empez por el primero de ellos, la
diescusin sobre la su utilizacin y los problemas que se tuvieron con l se avanzaron en el trabajo
anterior a este PFC y se resumen a continuacin.
Como se ha comentado el primero que se comenz a aplicar fue el modelo DO. Este modelo se extiende
a problemas de cualquier espesor ptico, (aL) y permite resolver tanto problemas de radiacin de
superficie contra superficie como problemas de combustin donde participa la radiacin.
Este modelo hace posible contemplar la radiacin que es absorbida a travs de un medio que participa
de la misma. Por ejemplo, el aire encerrado en la parte superior del absorbedor. Es decir, no solo
tiene en cuenta la radiacin absorbida o transmitida a travs de una pared, sino, que adems, al
definir un coeficiente de absorcin para cada banda del espectro del medio que rodea dichas
paredes, se puede modelar la energa que se dispersa.

122

Esto se puede ver ms claramente en la siguiente figura:

Figura 3. 68: Flujo de calor por radiacin

Y la ecuacin de transmisin de calor por radiacin para un medio que absorbe, emite y dispersa:

Figura 2.38: Ecuacin de transferencia de radiacin. Trminos.
Donde se tiene que (a+
s
) son el espesor ptico o la opacidad del medio.
5.1.2.1. Descripcin del modelo DO
Resuelve la ecuacin de transferencia radiante, para un nmero finito de ngulos discretos slidos, cada
uno asociado a un vector en la direccin s, para un sistema de coordenadas cartesianas. Si se modela
radiacin no gris, la ecuacin de transferencia de radiacin se puede escribir como:

123

Donde es la longitud de onda,
es el coeficiente de absorcin, y
es la intensidad de cuerpo negro

que viene dada por la funcin de Planck. Tanto los coeficientes de dispersin, como el ndice de
refraccin se consideran independientes de la longitud de onda.
Se puede implementar la radiacin no-gris en el modelo DO, dividiendo el espectro en N longitudes de
banda. Los intervalos para longitudes de banda pueden ser configurados por el usuario, y se
corresponden con valores en el vaco. El comportamiento de cada banda se asume gris.
Acoplamiento entre la ecuacin de la energa y el Modelo DO
El acoplamiento entre intensidades de energa y radiacin en una celda acelera la convergencia para un
esquema de volmenes finitos para transferencia de calor por radiacin. Mejora la convergencia para
aplicaciones con espesor ptico (logaritmo neperiano entre flujo de luminoso antes de atravesar la capa
y despus de atravesarla), mayor que 10. Esto es beneficioso cuando la dispersin es significante. En
nuestro caso esta opcin podra desactivarse.
Condiciones de contorno en paredes opacas.
El modelo DO, permite la especificacin de paredes opacas en el interior del dominio (con flujo
adyacente o zonas slidas en ambos lados de la pared), o externo al dominio (con un fluido adyacente o
una zona slida, slo). Las paredes opacas se tratan como cuerpo gris, si se computa radiacin gris, o no-
gris si el modelo no-gris est siendo usado.

Figura 3. 69: Condiciones de contorno en Paredes Opacas en el Modelo DO
Como se observa en la Figura 3. 69 hay una radiacin incidente
en el lado a de la pared opaca.

Parte de esa energa radiante es reflejada difusamente y directamente, dependiendo del factor de
difusin
para el lado a de la pared que se especifica como una condicin de contorno. Parte de la
energa radiante incidente es absorbida en la superficie de la pared y otra parte de la energa incidente
es emitida desde la pared. La cantidad de energa incidente absorbida por la superficie de la pared y la
cantidad emitida posteriormente dependen de la emisividad de la superficie y de la fraccin difusa. Para
un modelo NON-GRAY, se deber especificar una emisividad interna para cada longitud de onda. Es
obvio que la radiacin no es trasmitida a travs de la pared opaca.
124

La energa radiante incidente que impacta sobre la pared opaca su puede reflejar de nuevo al medio que
la rodea y absorbida por la pared. La energa que es reflejada puede ser: reflejada de manera difusa o/y
directa, dependiendo de la fraccin difusa
.Si la
es la cantidad de energa radiante incidente en

la pared opaca, entonces las siguientes cantidades son computadas por FLUENT como:
-Emisin desde la superficie de la pared =

-Energa difusa reflejada=
(1-

-Energa especular reflejada= (1-

-Absorcin en la superficie de la pared=

Donde
es la fraccin difusin, n es el ndice de refraccin del medio adyacente,
es la emisividad,
es la constante de Boltzman y
es la temperatura de la superficie de la pared

La absorcin en la pared asume que tanto la absortividad como la emisividad son iguales. Para una
pared puramente difusa,
es igual a 1 y no hay energa radiante directa. De forma similar, para una

pared puramente especular,
es igual a 0 y no hay energa difusa reflejada.

Condiciones de contorno en paredes semi-transparentes
FLUENT permite la especificacin de paredes semi-transparentes interiores y exteriores para el
modelo DO. En el caso de paredes interiores semi-transparentes, la radiacin incidente puede pasar a
travs de la pared y ser transmitida hacia el medio adyacente (y posiblemente refractada), esta puede
ser reflejada de nuevo en el medio circundante, y absorbida a travs de espesor de la pared. La
transmisin y la reflexin pueden ser difusas o directa. Pudindose especificar la
para toda la
radiacin transmitida y reflejada; el resto se trata como directa. Para paredes semi-transparentes
exteriores, hay dos posibles fuentes de radiacin en el contorno de la pared: un haz de irradiacin desde
fuera del dominio computacional y una radiacin incidente desde celdas adyacente a fluidos o zona
slidas.
Para radiacin no gris, las condiciones de contorno en paredes semi-transparentes son aplicadas en una
banda base. La energa radiante dentro de una banda es transmitida, reflejada y refractada como gris;
no hay transmisin, reflexin o refraccin de energa radiante de una banda a otra.
Por defecto la ecuaciones DO son resueltas en todas las zonas fluidas, pero en ninguna zona solida. Por
lo que en nuestro modelo se tiene que especificar las zonas solidas como participantes de la radiacin
en el panel SOLID como parte de las condiciones de contorno.

125

Paredes interiores semi-transparentes

Figura 3. 70: Condiciones de contorno en paredes semi-transparentes del Modelo DO
La pared se trata como si no tuviera espesor. La energa incidente
puede pasar a travs de la pared

semi-transparente si y solo si el fluido contiguo o la celdas solidas participan en la radiacin, as
permitiendo a la radiacin estar acoplada. El acoplamiento de la radiacin se establece cuando la pared
se especifica como semi-transparente. Habr que explcitamente especificar el acoplamiento de la
radiacin en la pared interior cambiando el tipo condicin de contorno a semi-transparente en el panel
WALL.
La energa incidente que es transmitida a travs de la pared puede ser transmitida especularmente o
difusivamente. La radiacin tambin puede ser reflejada hacia el interior de nuevo al medio que lo
rodea si el ndice de refraccin
para la zona fluida que representa el medio a es diferente del ndice

de refraccin
para el medio b. La energa reflejada al mismo tiempo puede ser reflejada

directamente o difusamente.
Si la pared semi-transparente posee espesor, entonces el espesor y el coeficiente de absorcin
determinan la absortividad de la pared delgada. Si tanto el espesor como el coeficiente de absortividad
se establecen como 0, entonces la pared no tiene absortividad.

Figura 3. 71: Condiciones de contorno en pared-semitransparente exterior en el Modelo DO
126

La figura muestra un caso general en el que una irradiacin
se aplica al exterior de una pared

semi-transparente con espesor nulo y un coeficiente de absorcin no nulo como propiedad del material.
Por esto la energa que absorbe es 0.
El flujo irradiado pasa a travs de la pared semi-transparente desde el exterior del dominio
computacional hacia el fluido adyacente o el medio solido a. La energa transmitida se puede refractar y
dispersar especularmente o difusamente, dependiendo del ndice de refraccin y de la fraccin difusa
que proporcionas como condicin de contorno.
Hay un flujo adicional ms all de
que se aplica cuando se seleccionan las condiciones de

contorno en la barra Thermal. Este flujo externo a la pared semi-transparente se computa con Fluent
como:

El haz de irradiacin se define por su magnitud, su direccin y el ancho que se especifique. La magnitud
se especifica en trminos de flujo de calor radiante. El ancho del haz como el ngulo slido donde se
distribuye la radiacin.
A este modelo se le aadi la carga solar a travs del algoritmo DO Irradation, calculando los parmetros
de iluminacin con la herramienta Solar- Calculator. Todo esto se hizo igualmente en el modelo
siguiente donde se explicar ms detenidamente.
El modelo anterior no dio buenos resultados por los siguientes motivos:
- En primer lugar el modelo DO va ms encaminado hacia problemas donde existen un volumen
de fluido encerrado que participa (absorbe, emite y atena la radiacin). En el colector este
fenmeno no tiene apenas importancia, ya que principalmente el aire entre el vidrio y el
absorbedor se calentar por que se calientan las paredes que lo contienen y por el efecto
invernadero que se produce en el mismo al no dejar escapar la radiacin reflejada en las
paredes del absorbedor.
- Por otra parte, en este modelo no se permite la caracterizacin de las superficies mediante se
transmisividad o absortividad, valores ampliamente conocidos.
- En este modelo se aplicaba la radiacin solar en la superficie semitransparente, con la intencin
de que dicha radiacin transmitida llegara al absorbedor, pero dicha radiacin no llegaba al
absorbedor, aunque paradjicamente si era transmitida.
- Otra particularidad del modelo, es que nicamente permite designar los volmenes que
participan en la radiacin, y no las superficies. De tal manera que por ejemplo al designar que el
volumen absorbedor participa de la radiacin, se incluye tanto su parte superior como su parte
inferior, cuando nosotros nicamente queremos observar los efectos de la radiacin solar
incidente en la cara de arriba del mismo.

5.1.3. Nuevo modelo de carga solar

A la vista de los resultados que se obtuvieron utilizando el modelo DO, se buscaron otras posibles
soluciones al problema de radiacin. Se analizaron trabajos que estuviesen utilizando los modelos de
radiacin con FLUENT en condiciones similares a las que tenemos en el captador solar y se encontr una
Tesis que utilizaba un modelo de carga solar con el mismo software CFD aplicado a la cabina de un coche
(Jonsson, 2007).
127

En este trabajo se analizaba la distribucin de temperaturas en el interior de un vehculo cuando estaba
sometido a la exposicin solar. El problema es bastante similar al que se trata en este PFC: la radiacin
solar llega el dominio computacional a travs de una superficie semi-transparante (cubierta de vidrio) y
posteriormente se va distribuyendo a lo largo de las superficies interiores, absorbindose y
reflejndose, haciendo que estas paredes se calienten y a su vez calentando el aire en el interior del
habitculo.
El aire que se encuentra por encima de la superficie del absorbedor en las simulaciones realizadas
vendra a ser el aire en la cabina, y las paredes internas de la misma seran la superficie superior del
absorbedor.
El modelo que se implementa en esta tesis trata al aire como un fluido que no participa de la radiacin,
es decir, no se produce dispersin de radiacin en su interior, y toda la energa radiante que se
transmite a travs de la cubierta de vidrio llega a las superficies internas. Est suposicin es bastante
acertada si consideramos al aire como un medio con escasa opacidad, de ah que se utilice el modelo de
radiacin S2S (Surface to surface) que describiremos ms adelante junto con un modelo de carga solar
(Solar Ray Tracing).
5.1.3.1. Descripcin del nuevo modelo
Para ir ms all del modelo simplificado, donde se introduca el efecto de la radiacin a travs de una
condicin de generacin interna de calor, es necesario implementar un modelo con el que se consiga
tener en cuenta: las propiedades pticas del captador, (que a su vez son funcin de la direccin y
longitud de onda de la radiacin incidente), los efectos de la reflexin, transmisin y absorcin de la
cubierta de vidrio, la disposicin geomtrica (efecto de sombra) y adems permita calcular la radiacin
incidente en funcin del lugar y la fecha en donde se encuentre el captador. Para esto se ha hecho uso
de una herramienta de FLUENT, el Solar Ray Tracing Algorithm, que se describir posteriormente.
Hay que resaltar que la herramienta anterior no se ocupa de la emisin desde la superficies del dominio,
y la componente reflejada de la carga solar incidente original se distribuye uniformemente a travs de
las superficies. Por esto para completar el modelo de carga solar, sera necesario emplear un modelo de
radiacin. El que ms se ajusta a esta aplicacin, es el modelo S2S o surface to surface, cuyo principal
inconveniente es el clculo de los llamados factores de forma, ya que suponen una gran demanda de
recursos. Por ello en el presente PFC no se ha implementado, pero sera interesante en trabajos
posteriores tenerlo presente, y verificar si las diferencias con este estudio simplificado son importantes.
Solar Ray Tracing Algorithm
Esta herramienta se puede utilizar para predecir la fuente de energa que resulta de la radiacin solar
incidente. Para ello utiliza el haz, que ha sido modelado usando la posicin del Sol y los parmetros de
iluminacin, y lo aplica a toda aquella pared o condicin de contorno de entrada o salida que se
especifique, llevando a cabo un anlisis de las posibles sombras y calculando finalmente el flujo de calor
que resulta de la radiacin incidente.
Este flujo de calor se acopla como un nuevo trmino a la ecuacin de la energa. La fuente de calor se
aade directamente a las celdas que rodean la superficie, asignndose en el siguiente orden: shell
conduction, celdas slidas y celdas fluidas. En principio el vector de posicin del Sol y la intensidad solar
se pueden introducir directamente, o bien a travs del Solar Calculator, que explicaremos despus.
Esta herramienta permite incluir los efectos tanto de la radiacin directa como de la difusa, definiendo
dos bandas espectrales para la radiacin directa (el infrarrojo y el visible) para especificar propiedades
del material diferentes en cada una de ellas, y una nica para la radiacin difusa. En el caso de
materiales opacos es necesario definir la absortividad en las dos bandas, y para semi-transparentes se
requiere tanto la absortividad como la transmisividad. Estos valores se especifican para el caso en el que
128

los rayos inciden de forma normal a la superficie, y el FLUENT, que interpola para obtener su valor segn
el ngulo de incidencia.
El algoritmo tiene en cuenta la carga difusa transmitida a travs de la cubierta de vidrio y tambin la
radiacin que se dispersa en el interior del dominio. Una fraccin de ese flujo de calor radiante, que el
software llama internal scattered energy se aplica a todas las superficies participantes en funcin de
su rea. Este valor de este parmetro depende de scattering fraction, que es la relacin de radiacin
directa que ha sido reflejada por las superficies opacas (despus de entrar a travs de una superficie
semi-transparente) y que se puede considerar que permanece en el interior y se puede distribuir
uniformemente entre las superficies que participan dentro del modelo de radiacin. A partir del valor
de internally scattered energy , se puede obtener el parmetro de ambient flux, que es el primero
entre la superficie total que participa de la radiacin.

Como se puede ver las principales desventajas de este modelo son:
- No tiene en cuenta la emisin desde las superficies participantes.
- La radiacin reflejada es distribuida uniformemente entre las superficies que forman el modelo.

Materiales semi-transparentes
Para modelar las propiedades pticas, como ya se ha dicho, es necesario definir la transmisividad y la
reflectividad del material. Hay que tener en cuenta, que dichas propiedades son dependientes, adems
de la longitud de onda de la radiacin y del ngulo de incidencia con variaciones significativas para
ngulos mayores de 40. Fundamentalmente al incrementarse el ngulo de incidencia desde cero, la
transmisividad se reduce, la reflectividad se incrementa y la absortividad inicialmente se incrementa
debido a la trayectoria ptica alargada y despus se reduce conforme ms radiacin incidente es
reflejada. Aunque no se puede asumir que todos las superficies transparentes tengan las misma
dependencia con respecto al ngulo.
Para vidrio con recubrimiento selectivo, la transmisividad y reflectividad espectral para cualquier ngulo
incidente se pueden aproximar en el modelo con respecto al ngulo normal de incidencia como
[Finlayson and Arasteh, 1993]:
Para la transmisividad:

donde,

Y la reflectividad:

129

Las constantes usadas en las ecuaciones anteriores son para vidrio con recubrimiento selectivo y la
fuente utilizada ha sido Finlayson, E.U., Arasteh, D.K., Huizenga, C., Rubin, M.D., Reilly, M.S. 1993.
Datos de entrada
Los parmetros que hay que definir son:
- Vector director solar (Solar Calculator)
- Irradiacin solar directa (Solar Calculator)
- Irradiacin solar difusa (Solar Calculator)
- Fraccin espectral: es la fraccin de radiacin solar incidente

- Absortividad (IR y visible), en paredes opacas
- Absrotividad (IR y visible) y transmisividad (IR y visible), en paredes semi-transparentes
- Absortividad y transmisividad (difusa), en paredes semi-transparentes
- Scattering fraction, se ha dejado el valor por defecto: 1
- Reflectividad en el suelo, se ha dejado el valor por defeco:0.2
Solar Calculator
Esta utilidad de FLUENT es capaz de calcular la direccin del haz de radiacin solar y la irradiacin para
una determinada hora, fecha y posicin. Estos valores se pueden utilizar como entradas inputs, para
introducirlo en el Solar Ray Tracing Algorithm. Los valores que debemos definir son:
- Posicin Global (Latitud, Longitud y Zona Horaria). En nuestro caso se ha introducido la posicin
par el captador que se encuentra colocado en la azotea del Aulario General I en el campus del
Paseo Alfonso XII.
- Fecha y hora de comienzo de la simulacin. En este caso solo no se han realizado simulaciones
transitorias por lo que, la fecha y hora inidicada se corresponden con las condiciones
estacionarias de la simulacin.
- Orientacin de la malla, para esto es necesario especificar la direccin Norte y Este del vector
de la malla. En el caso de que se hagan simulaciones donde el captador tenga cierto ngulo con
respecto a la horizontal, ser sumamente importante especificar correctamente la direccin del
vector de la malla.
- Mtodo de clculo de la irradiacin solar. El programa de CFD da dos posibilidades Fair
Weather Conditions y Theoretical Maximun method, que se describirn a continuacin.
- Y por ltimo, un factor solar, que simplemente es una reduccin lineal para la carga incidente
calculada.
Una vez que se introducen esos datos, se devuelven los siguientes parmetros:
- Vector director solar
- Irradiacin solar directa en la superficie terrestre
- Irradiacin solar difusa, para una superficie vertical y horizontal
- Irradiacin difusa reflejada en el suelo, para una superficie vertical
Mtodos para el clculo de la irradiacin
Como se ha indicado existen dos mtodos para calcular la irradiacin solar: Fair Weather Conditions y
Theoretical Maximun method. Aunque son similares hay una diferencia clave, el primero de ellos
impone una atenuacin en la carga solar representativa de condiciones atmosfricas claras pero no
130

totalmente despejadas, mientras que el segundo simplemente aplica a la radiacin extraterrestre un
factor de atenuacin.
Fair Weather Conditions
Segn el ASHRAE Handbook, la componente irradiancia E
t
es la suma de la radiacin directa E
D
, la difusa
E
d
, y la difusa reflejada en el suelo E
r
. La componente directa de la radiacin es el producto de la
irradiacin directa normal por el coseno del ngulo de incidencia entre los rayos solares incidentes y el
vector normal de la superficie:

La componente directa normal a la superficie terrestre
, se calcula como:

Donde, A es la irradiacin solar aparente con una masa de air m=0 y B es un coeficiente de extincin
atmosfrica, que se muestran en la tabla siguiente:

Tabla 3.12: Irradiacin solar extraterrestre
Para calcular la carga difusa sobre una superficie vertical se puede hacer mediante la siguiente frmula:

Donde C es una constante cuyo valor se muestra en la tabla anterior, Y es el ratio de radiacin difusa en
una superficie vertical con respecto a una superficie horizontal.
La ecuacin para la irradiacin difusa para superficies que no sean horizontales se calcula con:

Donde en el ngulo de la superficie con respecto al plano horizontal.
La ecuacin para la irradiacin solar reflejada en el suelo en una superficie viene dada por:

131

Donde
es la reflectividad del suelo. La carga difusa total en una superficie vendr dada por la suma de
en el caso de que la radiacin difusa se introduzca a travs del Solar Calculator.

Determinacin del ngulo de incidencia
La posicin del Sol en el cielo viene expresada en funcin de la altitud solar por encima de la
horizontal y por el azimut solar medido desde el Sur (Figura 3. 72). Estos ngulo a su vez dependen de
la latitud L; de la declinacin solar , y del tiempo solar aparente, expresado como el ngulo
horario H, donde H=0.25 (nmero de minutos desde el medio da local) en grados.
Las siguientes ecuaciones relacionan y con los parmetros que se han mencionado:

Figura 3. 72: ngulos solares para una superficie vertical y horizontal.
La figura anterior muestra los ngulos de posicin del Sol y los ngulos de incidencia para una superficie
vertical y horizontal. La lnea OQ conduce hasta el Sol, la NOS es la que va de Norte a Sur, y la EOW cruza
de Este a Oeste. La recta OV es perpendicular al plano horizontal en donde se encuentraN el azimut
solar HOS y el azimut de la superficie, POS, (). El ngulo HOP se definira como surface solar azimut ,
es decir:

El ngulo de incidencia para cualquier superficie se define como el ngulo entre los rayos solares
incidentes y un vector normal a la superficie. Para una superficie horizontal como se ve en la Figura 3.
72Error! No se encuentra el origen de la referencia. , el ngulo de incidencia
es QOV, y para una

superficie vertical,
, es QOP.
Para cualquier superficie, el ngulo de incidencia est relacionado con , y con el ngulo de
inclinacin de la superficie :
132

Cuando se tenga una superficie horizontal, , entonces

Y para una superficie vertical, :

Theoretical Maximum Method
Este otro mtodo, para calcular la irradiacin directa normal, utiliza una ecuacin obtenida de un cdigo
desarrollado por Laboratorio Nacional de Energa Renovable (NREL) de EEUU, llamado Solar Position
and Intensity Code (Solpos).

Donde
es la irradiacin directa normal extraterrestre y
es un factor de correccin para

tener en cuenta el efecto de absorcin y dispersin en la atmsfera.
Modelo de radiacin S2S
Hay que comenzar indicando que finalmente no se ha implementado este modelo junto con el Solar Ray
Tracing, es decir, en los resultados obtenidos en el PFC no se ha tenido en cuenta la radiacin emitida
desde las superficies que forman el propio colector. No se ha utilizado debido al gran consumo de
recursos y a la gran cantidad de tiempo necesaria para obtener soluciones que converjan, sin emargo
se va a realizar una breve descripcin del mismo. Tambin es cierto que los materiales que forman el
absorbedor suelen tener emisividades muy bajas, del orden de 0.05, por lo que la radiacin total emitida
por la superficie tampoco ser excesiva. nicamente los laterales del captador en el espacio entre e
absorbedor y la cubierta de vidrio van recubiertos de un material de alta emisividad.
El modelo de radiacin de superficie contra superficie (S2S) en combinacin con el algoritmo de trazado
de rayos solares son los que completan el modelo de carga solar. El modelo de radiacin se usa para
incluir la energa reflejada internamente procedente del Sol.
El modelo tiene en cuenta el intercambio de radiacin en un recinto de superficies grises. El parmetro
clave que afecta a este intercambio es la distancia de separacin entre superficies, el tamao de cada
una de ellas y la orientacin de las mismas. Estas variables son las que se tienen en cuenta para el
clculo de los factores de forma. Lo que ocurre es que el clculo de los view-factors, demanda una
gran cantidad de recursos, por lo tanto las superficies que se incluyen en la radiacin deben ser las
mnimas para que el tiempo de computacin no se dispare.
El modelo S2S usa la suposicin de que las caras son superficies grises y difusas, es decir, que las
propiedades de las superficies son independientes de la longitud de onda, y que segn la ley de
Kirchhoff, la emisividad y la absortividad para cada banda son iguales. Otra de las suposiciones del
modelo de radiacin es que el medio que separa las superficies no tiene influencia en la energa
irradiada G; la energa incidente dentro de cada superficie ser una parte reflejada (G), o absorbida
(G) y otra fraccin puede ser transmitida (G). Como la mayora de las superficies (excepto el vidrio que
133

cubre el colector) son tratadas como opacos la transmisividad puede ser eliminada y entonces se
cumple:

Y, como hemos dicho, es igual a la emisividad, por lo que se tiene:

La ecuacin de la energa usada en el modelo S2S est compuesta de la energa directamente emitida y
reflejada. La energa reflectada por cada superficie es independiente de la energa incidente en la
superficie por lo que puede ser expresada como una funcin del flujo de energa que deja todas las
superficies radiantes. El flujo de energa que deja la superficie K se describe como:

donde
representa el flujo de energa que deja la superficie,
el flujo de energa incidente

sobre la superficie,
es la reflectividad de la superficie,
es la emisividad y es la constante de
Boltzman y T es la temperatura.
La cantidad de energa incidente desde una superficie a otra superficie ser directamente proporcional a
los factores de forma. Los factores de forma
es la fraccin de energa incidente que deja la superfice

k e incide sobre la superficie j. Por lo que el flujo incidente de energa
puede ser expresado en

trminos de la energa que deja todas las otras superficies:

Donde
es el area de la superficie k y
es el factor de forma entre la superficie k y la superfice j.

Para N superficie usando la denominada relacin de reciprocidad, dada por :

Por lo que:

Entonces la ecuacin del flujo de energa saliente de una superficie se puede expresar como:

134

5.1.3.2. Conclusiones y modelo utilizado

El modelo ms completo como ya se ha anticipado es el que tiene en cuenta por un lado la carga
incidente, estamos hablando de radiacin de longitud de onda corta junto con la radiacin de onda
corta que se emite en el absorbedor y se va reflejando por las superficies que forman el modelo.
Finalmente, debido al gran consumo de recursos computacionales que supone el clculo de los factores
de forma para mallas grandes, solo se ha implementado la herramienta el modelo implementado ha
sido el del Solar Ray Tracing sobre el modelo, de manera que no tenemos en cuenta la radiacin emitida
por el absorbedor y la radiacin que es reflejada junto con la radiacin difusa que entra a travs de la
cubierta de vidrio, se distribuye en funcin del rea de cada pared.

5.2. Modelo de carga solar
El modelo Solar Ray Tracing Algorithm, no tiene elevados costes computacionales y bsicamente
aplica trminos de generacin de calor dentro de la ecuacin de la energa en los elementos son
adyacentes a la condicin de contorno donde se implementa este modelo.
A la hora de aplicar el trmino de generacin de energa sobre las celdas existe un orden preestablecido,
primero se aplica sobre celdas pertenecientes a zonas slidas, y despus sobre aquellas que forman
parte de zonas fluidas, pero nicamente a una de las dos. Por ejemplo, para el caso de la superficie
superior del absorbedor, el trmino de generacin interna se aplica sobre las celdas que son parte del
slido absorbedor, y no en las celdas que forman parte del volumen superior de aire.
El trmino que introduce en la ecuacin de la energa depende de los parmetros de radiacin directa,
difusa, el vector de posicin del Sol, y tambin la fraccin del espectro de radiacin incidente que
pertenece al visible y cual al infrarrojo. Estos parmetros se pueden introducir de diversas formas,
como a continuacin veremos.
Para calcular finalmente que trmino de calor se aplica al elemento adyacente a una condicin de
contorno dada, adems de tener en cuenta los parmetros que hemos definido en el prrafo anterior,
FLUENT cuenta con un algoritmo para tener en cuenta las sombras que se generan entre las superficies
que forman el modelo. Este algoritmo funciona de la siguiente forma: para la superficie que queramos
estudiar (llammosla superficie 1), se traza desde su centroide un rayo en direccin a la posicin del Sol,
si este rayo corta a otra superficie (superficie 2) y adems est superficie 2 est enfrente de la superficie
1, entonces se dice que si la superficie 2 es opaca tapa completamente la superficie 1. Si fuera
semitransparente, entonces atena la radiacin incidente.

135

5.2.1. Implementacin del modelo de carga solar
Parmetros de la irradiacin Solar

En primer lugar dentro del men de modelos, en el apartado de radiacin, existe un subapartado
llamado Solar Load, donde seleccionamos Solar Ray Tracing.
Este modelo de carga solar se puede utilizar simultneamente con un modelo de radiacin (Rosseland,
P1, S2S) o bien, como en el caso que nos ocupa, de manera independiente. Ya se comento en el
capitulo anterior que se podra mejorar el modelo utilizando el modelo S2S para tener en cuenta la
energa emitida por las propias superficie que forman la simulacin.
Una vez seleccionado el modelo de carga solar, se deben introducir los parmetros de radiacin:

Figura 3.73: Introduccin de los parmetros de radiacin
- En primer lugar tenemos que elegir la posicin del Sol, este parmetro se puede introducir de
forma manual o usar la herramienta Solar Calculator que ya introducimos en el captulo
anterior.
- En segundo lugar se especifican los parmetros de iluminacin:
o Irradiacin solar directa, que tiene en cuenta la cantidad energa por unidad de rea
en W/m
2
debida a la radiacin solar. Este parmetro se puede introducir utilizando el
Solar Calculator, o bien especificar un valor constante.
o Irradiacin solar difusa, tiene en cuenta la cantidad de energa en W/m
2
debida a la
radiacin solar difusa. Al igual que el caso anterior se puede definir un valor constante
o bien utilizar el Solar Calculator.
o Y por ltimo, se define Spectral Fraction, que es la proporcin de radiacin incidente
que pertenece al espectro de radiacin visible. Se toma como valor por defecto 0.5.
Para las simulaciones realizadas, no se especificaron valores constantes para los parmetros anteriores
sino que se us la herramienta Solar Calculator que hemos descrito en el apartado anterior. En la
figura siguiente se muestran los datos de entrada.
136

Figura 3. 74: Lugar, fecha, hora y orientacin de la malla para las simulaciones. Solar Calculator.
Como el objetivo de estas simulaciones era el de comprobar el funcionamiento del modelo de carga
solar nicamente se han desarrollado casos estacionarios para una hora y fecha concretas, buscndose
que la orientacin del Sol fuese lo ms perpendicular posible, y en una poca del ao con alta
irradiancia.
En cuanto a la posicin, se refiere a la Latitud y Longitud del banco de ensayos de captadores solares
situado en el Campus Alfonso XIII de la universidad.
La orientacin de la malla es un punto crtico para el clculo del vector de posicin del Sol y los
parmetros de iluminacin. Para definirlo tenemos que tener en cuenta donde se sita el origen de
coordenadas de nuestra malla y como hemos definido los ejes en el espacio, para entonces establecer la
posicin relativa de la misma con respecto a los puntos cardinales. Se puede comprobar en la Figura 3.
74, que en las simulaciones la malla est en posicin horizontal. Si en sucesivos trabajos se quiere
cambiar al ngulo del captador con respecto a la horizontal ser necesario, adems de girar la malla en
s, determinar cul es su orientacin.
Para terminar definir parmetros hemos de elegir el mtodo para el clculo de la irradiacin. Se ha
elegido el Fair Weather Method, por ello puede ser que en las simulaciones a pesar de hacerlas para el
mes de Julio, la radiacin puede parecer baja.

Condiciones de contorno
Una vez que se han definido los parmetros del modelo de carga solar, se tienen que establecer cules
son las condiciones de contorno que participan de l. Teniendo en cuenta la geometra del captador y
las caractersticas de los materiales por los que est compuesto, se pueden extraer diversas
conclusiones.
En primer lugar, la zona inferior del tubo se encuentra aislada, los materiales aislantes adems de baja
conductividad, son opacos y normalmente con absortividades bastante bajas por lo que podemos
considerar que dichas zonas no van a participar del intercambio de energa radiante. Adems en las
simulaciones que estamos realizando la malla est en posicin horizontal, por lo que se supone que
nunca le llegara radiacin solar.
Si la frontera inferior del modelo no se tiene en cuenta la radiacin incidente, todas aquellas superficies
que se encuentran por debajo del absorbedor, lgicamente tampoco recibirn ningn tipo de radiacin.
137

A s mismo, la zona de entrada y salida de fluido, tampoco se ha incluido dentro de la carga solar puesto
que, en condiciones normales, existir un colector de entrada y salida que cubrirn dichas secciones. En
la Figura 3.75, se puede ver en color azul las fronteras que no forman parte de la carga solar.
Por tanto las nicas superficies sobre las que incide la carga solar, son las que forman la cavidad superior
en el colector, es decir, la cubierta de vidrio, que se trata como superficie semi-transparente, las
paredes verticales y la superficie superior del absorbedor, que son opacas. Posteriormente cuando se
calcule el trmino de generacin sobre cada una de ellas se tendr en cuenta como ya se ha dicho su
posicin relativa, y tambin generacin de sombras entre ellas.

Figura 3.75: Superficies no irradiadas

En cada una de las superficies que forman el modelo de carga solar tenemos que definir la absortividad
para el espectro de radiacin infrarroja y para el espectro visible. Esto se tiene que hacer para fronteras
opacas ( Figura 3.76). En el caso de fronteras semitransparentes, tambin es necesario definir la
transmisividad para la radiacin directa ,en el visible y en el infrarrojo, y adems se debe proporcionar
otro parmetro, lo que llama FLUENT diffuse hemipherical, que es la proporcin de radiacin difusa
que se transmite y se absorbe por parte de la cubierta de vidrio (Figura 3. 77).
138

Figura 3.76: Carga Solar para paredes opacas. Definicin de la absortividad

Figura 3. 77: Absortividad y transmisividad para una superficie semitransparente
Para calcular la reflectividad en cada superficie simplemente:

Si observamos las diferencias en la definicin de absortividad entre una superficie opaca y una
transparente podemos ver cmo a pesar de que la superficie opaca absorbe radiacin difusa no
podemos especificar el parmetro de diffuse hemispherical. Esto es por lo que se explicaba en el
apartado anterior, y es que la radiacin que se refleja en las superficies junto con la radiacin difusa
que se transmite a travs del vidrio forman la llamada internal scattered energy que se aplica a las
superficies que participan de este modelo segn el rea de cada una de ellas.

- Modelo de radiacin:
o Sin modelo de radiacin
o Carga Solar: Solar Ray Tracing
139

- Gravedad: (0,-9.81,0) (m/s
2
)
Las condiciones de contorno para el modelo con radiacin son las mismas que las que se impusieron
para el modelo sin radiacin, con la diferencia de que en este caso no se incluye el trmino de
generacin interna de calor en el absorbedor. A cambio de esto se han de especificar los parmetros de
absortividad y transmisividad para las condiciones de contorno que participan del modelo de carga
solar.
En las dos figuras anteriores se pueden comprobar estos parmetros tanto para las superficies que son
opacas como para la cubierta de vidrio (semitransparente). Estos valores se han tomado ms o menos
en concordancia con los valores normales, para posteriores simulaciones, se podrn ajustar los
parmetros para coincidir plenamente con los datos del banco de ensayos experimental.

Frontera Condicin de
contorno
trmica
Solar
Ray
Tracing
Absortividad Transmisividad
IR Visible IR Vis HD
Absorbedor_cara_inferior Wall (acoplada) NO
Absorbedor_cara_inferior-shadow Wall (acoplada) NO
Absorbedor_cara_inferior:001 Wall (acoplada) NO
Absorbedor_cara_inferior:001_shadow Wall (acoplada) NO
Absorbedor_cara_superior Wall (acoplada) SI 0.8 0.8 - - -
Absorbedor_cara_superior-shadow Wall (acoplada) SI 0.8 0.8 - - -
Absorbedor_frontal Wall (adiabtica) NO
Absorbedor_lado_I Wall (adiabtica) NO
Absorbedor_lado_II Wall (adiabtica) NO
Absorbedor_trasera Wall (adiabtica) NO
Air_inf_front Wall (adiabtica) NO
Air_inf_ladoI Wall (adiabtica) NO
Air_inf_ladoII Wall (adiabtica) NO
Air_inf_tras Wall (adiabtica) NO
Air_sup_front Wall (adiabtica) SI 0.8 0.8 - - -
Air_sup_ladoI Wall (adiabtica) SI 0.8 0.8 - - -
Air_sup_ladoII Wall (adiabtica) SI 0.8 0.8 - - -
Air_sup_tras Wall (adiabtica) SI 0.8 0.8 - - -
Cara_exterior_tubo Wall (acoplada) NO
Cara_exterior_tubo-shadow Wall (acoplada) NO
Cara_exterior_tubo:34 Wall (acoplada) NO
Cara_exterior_tubo:034-shadow Wall (acoplada) NO
Cara_interior_tubo Wall (acoplada) NO
Cara_interior_tubo-shadow Wall (acoplada) NO
Cristal Wall(conveccin) SI 0.1 0.1 0.8 0.8 0.8
Entrada_de_fluido Velocity-inlet NO
Entrada_de_tubera Wall (adiabtica) NO
Salida_de_fluido Pressure-outlet NO
Salida_tubera Wall (adiabtica) NO
Soldadura_frontal Wall (adiabtica) NO
Soldadura_lado_I Wall (acoplada) NO
Soldadura_lado_I-shadow Wall (acoplada) NO
Soldadura_lado_II Wall (acoplada NO
140

Soldadura_lado_II-shadow Wall (acoplada) NO
Soldadura_trasera Wall (adiabtica) NO
Suelo Wall (adiabtica) NO
Tabla 3.13: Condiciones de contorno para el modelo completo con un solo tubo con carga solar

5.3.2. Materiales
Material Tipo Propiedades
Agua Fluido Funcin de la temperatura (Ver 1.2.4) (Anexo I)
Aire Fluido Constantes excepto densidad (Boussinesq)
Cobre Slido Constantes (Anexo I)
Aluminio Slido Constantes (Anexo I)
Vidrio Slido Constantes (Anexo I)
Tabla 3. 14: Materiales en el modelo completo de un solo tubo con carga solar
5.3.3. Parmetros de iluminacin y vector de posicin del Sol
En las simulaciones con este modelo se han hecho para el da 1 de Julio, a las 11:00 horas, para una
Latitud de: 37.606 y Longitud 0.97854. Introduciendo estos valores en el Solar Calculator e indicando
el mtodo para el clculo de los parmetros, se devuelve la siguiente informacin:
Fair Weather Conditions:
- Sun Direction Vector: X: -0.106432, Y: 0.920328, Z: -0.376389
- Sunshine Fraction: 1
- Direct Normal Solar Irradiation (at Earth's surface) [W/m^2]: 867.887
- Diffuse Solar Irradiation - vertical surface: [W/m^2]: 142.976
- Diffuse Solar Irradiation - horizontal surface [W/m^2]: 117.454
- Ground Reflected Solar Irradiation - vertical surface [W/m^2]: 91.6194

Estos datos son los inputs para el Solar Ray Tracing, en funcin de stos, de las condiciones de contorno
que se han impuesto y de la posicin relativa de las superficies con respecto a vector solar, calcula los
trminos de generacin de calor que se aplican sobre cada una de las fronteras. Al inicializar el caso se
muestra la siguiente informacin:
Internally Scattered Energy [W]: 49.4996, Ambient Flux [W/m^2]: 93.6961
Boundary ID: 10, Integral Energy Source [W]: 3.429276e+000
Boundary ID: 15, Integral Energy Source [W]: 2.304923e-001
Total Integral Energy Source [W]: 1.844054e+002
Compute Time: 4 sec

El primer parmetro que calcula es la Internally Scattered Energy, que es la suma de la radiacin
difusa que se introduce en el interior del modelo a travs de la cubierta de vidrio y la radiacin reflejada
por todas las superficies. Y el otro parmetro Ambient Flux se calcula simplemente dividiendo el
parmetro anterior por la suma de las superficies que participan de la carga solar, es decir, las
superficies laterales (0.123*0.02*2+1.83*0.02*2), la superficie del vidrio (1.83*0.123) y la superficie del
absorbedor (1.83*0.123), la suma total es de: 0.5283 m
2
. Dividiendo est 49.996 W, entre el valor
anterior se obtiene los 93.6961 W /m
2
.
141

A continuacin se establece el trmino de fuente de generacin de energa en cada una de las
superficies que se han seleccionado para formar parte del Solar Ray Tracing. La que mayor valor tiene
es la ID:47, porque es la cara superior del absorbedor, y como era de esperar es la que ms energa va a
absorber. La cubierta de vidrio, ID:13, tiene un trmino bastante inferior al del absorbedor, a pesar de
que es la superficie con mayor radiacin incidente, esto se debe a que es semitransparente y solo una
pequea parte es absorbida, con lo que se comprueba que efectivamente la herramienta para
determinar las sombras funciona. El resto de superficies son las paredes laterales del colector, y por la
posicin del Sol (prcticamente perpendicular a la superficie el absorbedor) el trmino de energa es
muy pequeo.
Si tenemos en cuenta cmo funciona el modelo de carga solar, en todas las superficies
independientemente de que no les llegue radiacin solar directa siempre estarn afectadas por el
Ambient Flux, por ejemplo la ID:15, que es la superficie vertical sobre la entrada de fluido, su rea es
0.02*0.123=2.46x10
-3
m
2
, y multiplicada por el Ambient Flux, se obtiene precisamente el trmino de
generacin de energa: 0.2305 W. Por lo tanto es una superficie donde no llega radiacin directa.

El objetivo del modelo era la verficacin del modelo de carga solar, por lo que nicamente se han hecho
simulaciones estacionarias para la fecha indicada antes, con el gasto nominal de 0.02 kg/s.m
2
, variando
la temperatura de entrada del fluido y tambin la velocidad del viento. La distribucin de casos es la
siguiente:
Caso Temperatura de
entrada (C)
Velocidad del Viento
(m/s)
(m/s)
Julio_Te_20 20 4 0.117204696
Julio_Te_30 30 4 0.117497862
Julio_Te_40 40 4 0.117891623
Julio_Te_50 50 4 0.118380874
Julio_Te_60 60 4 0.118960782
Julio_Te_70 70 4 0.119626691
Julio_Te_20_v8 20 8 0.117204696
Julio_Te_30_v8 30 8 0.117497862
Julio_Te_40_v8 40 8 0.117891623
Julio_Te_50_v8 50 8 0.118380874
Julio_Te_60_v8 60 8 0.118960782
Julio_Te_70_v8 70 8 0.119626691
Tabla 3.15: Simulaciones realizadas en modelo completo con un solo tubo

El objetivo del modelo, como ya se ha comentado, es el de comprender como se implementaba el
modelo de carga solar sobre el mismo, verificar como se distribuye la radiacin solar en el dominio
computacional, y tratar de comparar los resultados tanto de la conveccin mixta para el flujo interno en
el conducto como los que se tienen sobre la conveccin natural en el aire.
Por ltimo, tambin se va a calcular la eficiencia del captador, teniendo en cuenta que ahora si tenemos
los parmetros de iluminacin, es decir, contamos con la energa incidente G, y no como ocurra en el
modelo anterior donde trabajbamos con la potencia absorbida S.
142

5.4.1. Anlisis de la transmisin de calor
5.4.1.1. Nmero de Nusselt en el flujo interno
Al igual que se ha hecho para los casos anteriores vamos a representar el nmero de Nusselt frente a la
distancia adimensional (Figura 3.78: Nmero de Nusselt para las simulaciones con colector completo 1
tubo con radiacin. Comparacin con las correlaciones de Petukhov y Polyakov [1988] de conveccin
mixta y forzada.Figura 3.78). La tendencia es la misma de la se obtena en los casos donde no
implementbamos el modelo de radiacin.
Si representamos ahora el Nusselt para flujo desarrollado junto con la correlacin que se obtuvo en las
simulaciones de tubo con absorbedor, se obtienen resultados que se ajustan muy bien a los resultados
de modelos anteriores por lo que se puede concluir que, las simulaciones aunque tenga implementado
un modelo de carga solar, siguen obteniendo buenas soluciones para el flujo interno.

Figura 3.78: Nmero de Nusselt para las simulaciones con colector completo 1 tubo con radiacin. Comparacin
con las correlaciones de Petukhov y Polyakov [1988] de conveccin mixta y forzada.

1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Simulacin Julio_Te_60
Simulacin Julio_Te_20
Petukhov y Polyakov [1988], Conveccin Mixta, Julio_Te_60
Petukhov y Polyakov [1988], Conveccin Mixta, Julio_Te_20
143

Figura 3. 79: Nmero de Nusselt desarrollado para modelo de 1 tubo con radiacin. Comparacin con los
resultados numricos para el tubo con absorbedor.
Los resultados son coherentes con lo obtenido anteriormente con otros modelos sin radiacin por lo
que se pude concluir que el modelo de carga solar ha sido adecuadamente implementado y proporciona
resultados coherentes.

5.4.1.2. Conveccin natural en el aire
Las condiciones de contorno que dominan este fenmeno ahora han cambiado. Las paredes verticales
que delimitan la cavidad entre la cubierta de vidrio y absorbedor siguen estando aisladas, sin embargo,
al estar incluidas dentro del modelo de carga solar, se les aade un trmino de generacin de energa a
los elementos adyacentes que forman parte de la cavidad de aire. Esto hace que la temperatura en los
extremos se incremente de manera notable ya que el nico camino de fuga del flujo de calor es
precisamente a travs del aire.
En cuanto a la pared que forma la cubierta de vidrio, la condicin de contorno sigue siendo convectiva,
pero como en el caso de las paredes verticales, dado que se ha definido un valor de absortividad de 0.1,
parte de la radiacin incidente es absorbida por la misma. Por lo que tambin hay un trmino de
generacin en las celdas de aire adyacentes a esta pared.
Como consecuencia de ello, el patrn que se obtiene cambia totalmente, ya que las temperaturas
elevadas de los extremos del absobedor favorecen el flujo ascendente, por esto en todos los casos
analizados, siempre se dan dos celdas de recirculacin.

4
1.00E+05 1.00E+06
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Modelo completo con radiacin
Poliakov y Petukhov [1988]
Resultados tubo+absorbedor
144

Figura 3. 80: Campos de temperatura y vectores de velocidad para la cavidad de aire entre absorbedor y cubierta
de vidrio. Caso completo 1 tubo con radiacin.

La figura anterior muestra diferentes secciones para un caso con una temperatura de entrada de 20C,
se encuentras a 0.6, 1.2y 1.7 metros desde la entrada. Se ve como en este caso solamente existen dos
reciculaciones alcanzndose las temperaturas ms elevadas en los extremos. A pesar de que, lo que
podra parecer es que si el absorbedor es donde el trmino de energa va a ser mayor, su superficie
tendra que estar ms caliente. Sin embargo, la accin el flujo en el conducto va removiendo el calor del
mismo, mientras que en las paredes verticales adiabticas, todo el flujo de calor va a parar al aire, que
por consiguiente termina calentndose.
En cuanto a la simetra de las circulaciones se debe a que, prcticamente la potencia que aplica el
modelo de carga solar es prcticamente el mismo en las paredes laterales. En el caso de que la posicin
del Sol cambiase, la potencia asignada entonces tambin podra cambiar, y por lo tanto variar el patrn
de flujo.

145

Al igual que se hiciera con el modelo completo con generacin interna de calor, vamos a analizar la
eficiencia del captador utilizando las mismas expresiones que ya se han expuesto. La principal diferencia
entre las simulaciones presentes y las anteriores es el modelo de carga solar.
Antes definamos directamente cual era la potencia absorbida por el absorbedor, por lo que no
podamos tener en cuenta el efecto de los parmetros pticos, la distribucin de radiacin sobre el
absobedor, el efecto de las sombras, o los parmetros de iluminacin, tanto la radiacin directa como la
difusa. El modelo de carga solar adems de lo anterior nos permite modelar el comportamiento de una
superficie semitransparente, la atenuacin que produce sobre la radiacin incidente, as como la
absorcin de radiacin de la propia superficie.
Para el clculo de la eficiencia es necesario, por un lado obtener el calor til, que se puede calcular con
facilidad, pues se disponen tanto de las temperaturas como de los caudales. Pero tambin necesitamos
conocer cul es la radiacin incidente en el plano de captacin. Este dato no lo tenemos directamente,
ya que como se explico en el punto 5.3.3, al inicializar el caso FLUENT calcula cual es la potencia (W) que
asigna a cada superficie, pero esta cantidad de energa es la que realmente llega a absorber la pared en
cuestin. Es decir, no nos sirve para obtener G, para ello tenemos que utilizar los outputs que
devuelve el Solar Calculator cuando le especificamos cul es la localizacin, la fecha y la orientacin de
la malla.
Como sabemos la radiacin incidente ser la suma de la radiacin directa y de la difusa. Para obtener la
radiacin directa incidente debemos tener en cuenta que FLUENT proporciona el dato de: irradiacin
normal directa en la superficie terrestre (W/m
2
) y el vector de posicin del Sol. La radiacin incidente
ser la proyeccin del vector normal a la superficie, en la direccin del vector solar. Por tanto hay que
realizar el producto escalar de el vector de posicin del Sol por el vector normal al plano de incidencia y
a su vez multiplicar por la irradiacin directa. Para el caso que nos ocupa, si el captador se mantiene
horizontal, tenemos:

A la radiacin directa, hay que sumarle la difusa, que al ser independiente de la direccin, simplemente
tendremos que sumar al valor anterior, la irradiacin difusa para una superficie horizontal, que es otro
de los parmetros de iluminacin que devuelve FLUENT.
De esta forma, la energa incidente para las simulaciones que se han realizado es:

146

Figura 3.81: Eficiencia del captador para modelo completo de1 tubo con radiacin para velocidades de viento de
4 y 8 m/s. Comparacin con banco de ensayos.
Como se aprecia las diferencias entre el captador experimental y las simulaciones realizadas son
elevadas, pero hay que tener en cuenta que los parmetros pticos son muy distintos. El producto
transmistanca-absortancia para el captador experimental es igual a 0.93*0.95=0.883, mientras que en
los resultados numricos se defini una absortividad en el absorbedor de 0.8 y una transmisividad
tambin de 0.8 por lo que =0.8*0.8=0.64.
Si nos fijamos el valor de las curvas de eficiencia en el eje de ordenadas, se comprueba como
efectivamente se encuentran con valores muy prximos al producto transmitancia-abosortancia, lo que
coincide con la frmula de eficiencia que se da en el libro de Duffie y Beckman:

Teniendo en cuenta que normalmente el factor de remocin de calor se mueve dentro del intervalo 0.8-
0.9.

Adems se ha comparado con el modelo analtico, que ya se utiliz en el apartado 4, se observa como
prcticamente la diferencia entras las curvas es muy parecida a la que se tena para el modelo sin
radiacin. Por otra parte tambin se comprueba como la velocidad del viento tiene una muy baja
influencia en los resultados que se obtenan, aunque lgicamente, se obtienen valores ligeramente
superiores cuando la velocidad es menor.

0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
-0.01 0.00 0.01 0.02 0.03 0.04 0.05
E
f
i
c
i
e
n
c
i
a

d
e
l

c
a
p
t
a
d
o
r

Temperatura reducida, T*
Resultados nmericos v=4m/s
Captador experimental
Resultados experimentales v=8m/s
Modelo analtico
147

Se puede concluir diciendo que la diferencia entre pendientes es igual que en el apartado anterior, es
decir, en el experimental hay prdidas por conveccin por los laterales y el fondo del captador que no se
tienen en cuenta en el modelo numrico. Y por otra parte se tiene una buena referencia al comparar el
valor de corte de la curva en eje de ordenadas con el producto , son prcticamente iguales lo que
coincide con las curvas en colectores comerciales. Y por otra pare se ha demostrado la poca incidencia
de la velocidad del viento en el modelo.

148

6. Modelo completo 9 tubos con carga solar

Una vez que hemos analizado el funcionamiento del modelo de carga solar en una malla simplificada
para un solo conducto, vamos a proceder a crear la geometra para el captador completo que consta de
9 tubos. Bsicamente lo que se busca es disponer de un modelo lo ms parecido posible al captador
solar estndar estudiado experimentalmente y por otro lado crear una malla que pueda servir en
posteriores trabajos para analizar el reparto de caudales entre los distintos tubos, y estudiar el captador
desde el punto de vista de las prdidas de presin que se producen en cada elemento, as como estudiar
la posibilidad de incluir en los tubos algn tipo de dispositivo insertado para estudiar la mejora en la
transferencia de calor que se produce.
El mayor inconveniente de crear un modelo completo es que el coste en los requisitos de memoria
necesarios y el tiempo de simulacin requerido hasta obtener una solucin convergida. Para solventar
el problema se optado por crear la mitad de la geometra e imponer un condicin de simetra, y adems,
ajustar el tamao del mallado utilizado en el modelo de un solo conducto. En los siguientes apartados se
explicar el proceso seguido.
En cuanto a las condiciones de contorno y los parmetros de radiacin son los mismos que los que se
aplicaron para el modelo de un solo conducto con carga solar. Se han realizado tambin simulaciones
siguiendo la misma planificacin y se ha ampliado el estudio considerando el que caudal total circulante
es del 50% y del 150% respecto del gasto msico nominal.

6.1.1. Definicin de la geometra
En un primer momento se pens en utilizar el mallado que se haba utilizado para el modelo completo
de 1 solo tubo (Figura 3.48), que tiene aproximadamente un milln de nodos, por lo que la geometra
completa contara con ms de nueve, algo inasumible desde el punto de vista del coste computacional.
Por lo que como solucin alternativa, pensando en poder aprovechar la malla que ya se haba utilizado y
estaba suficientemente testada, y observado que los resultados prcticamente eran simtricos (con la
excepcin de la radiacin incidente) en los casos anteriores se propuso eliminar la mitad de la geometra
y aplicar una condicin de simetra con respecto a un plano YZ que cortase al captador en dos partes
iguales.
La construccin de la geometra parti de la que ya se tena y los pasos fueron:
1) En primer lugar se deba cortar por la mitad malla anterior. Para esto en GAMBIT se cre una
superficie que divida en dos la geometra y utilizando la herramienta split volume , se
creaban dos volmenes diferenciados.
149

Figura 3.82: Opcin "Split Volume". GAMBIT
Se puede ver en la figura anterior como para dividir un volumen es necesario especificar, el
volumen que se quiere cortar y cul es el tipo de topologa con la que se usa para ello. En
nuestro caso se va a seleccionar una cara real (la que se acababa de crear), de manera que
como la superficie intercepta todo el volumen el resultado es otros dos volmenes separados.
Cuando usamos este comando se pueden seleccionar otras subopciones , una es Retain, que
hace que GAMBIT mantenga la superficie que se ha utilizado para cortar el volumen, para
nuestro modelo no es necesario. Otra opcin, Connected permite que los volmenes creados
estn conectados, pero como posteriormente se va a eliminar uno de los volmenes creados
tampoco se selecciona. Y por ltimo, Bidireccional, que hace que se pueda cortar la propia
superficie que se utiliza como herramienta de corte (Figura 3. 83), para nuestra malla no hace
falta tenerla en cuenta, ya que directamente no se mantiene la propia superficie.
La geometra resultante se puede ver en la

Figura 3. 83: Opcin "Bidireccional" en la herramienta "Split". GAMBIT
150

Figura 3. 84: Geometra resultante despus de utilizar "Split Volume"

2) Una vez que tenemos la zona central de la malla simtrica, que consta de 4,5 tubos, solo
tenemos que ir copiando la malla que ya tenamos tantas veces ms segn el sentido negativo
del eje X.
Utilizando Move/Copy Volumes, podemos copiar los cuatro mdulos que faltan. Para ello es
necesario seleccionar que operacin queremos realizar, que es la de trasladar el volumen una
distancia dada, en este caso, 0.123 metros, que es el ancho de cada mdulo. A la hora de copiar
la geometra, seleccionaremos la opcin copy mesh linked, de forma que GAMBIT copiar la
cara y su malla asociada y enlazar malla copiada con la malla en la cara original.

Si procedemos de esta forma, finalmente obtenemos la siguiente geometra:

6.1.2. Proceso de mallado

Con la geometra anterior, sin modificar el mallado que ya exista en cada mdulo obtenemos un malla
con cerca de cinco millones de nodos. Por esto, a pesar de haber reducido el tamao de la misma a la
mitad imponiendo una condicin de simetra, todava sigue siendo una malla muy pesada. A la hora de
Figura 3. 85: Geometra completa
151

aligerar la misma se propusieron diversas soluciones, en funcin de los objetivos de las simulaciones y la
precisin que se quera alcanzar.
Una de las soluciones que se implement fue la de sustituir la malla estructurada de la cavidad de aire
entre cubierta de vidrio y absorbedor por una no estructurada con elementos de mayor tamao. El
origen de esta propuesta es que, ya vimos en el apartado 3 de este captulo que reducir la malla de cada
seccin transversal de esta cavidad a 300 nodos tiene un error del orden del 5%, que es bastante
asumible si tenemos en cuenta que el principal objetivo de estas simulaciones es estudiar la eficiencia y
no tanto los fenmenos de conveccin natural que se dan en esta zona. La malla final qued como:

Figura 3. 86: Seccin transversal de la malla para el captador completo.
Sin embargo, esta solucin no supona un ahorro suficiente ya que no se poda disminuir demasiado el
nmero de nodos en el aire porque la malla del absorbedor que es ms densa, hace que en las cercanas
se tenga una elevada densidad lo impeda, pues stas se encuentran conectadas y por el contrario en las
paredes verticales los nodos se encuentren muy espaciados.
Finalmente si queramos reducir la malla an ms, hasta llegar a valores razonables del orden de tres
millones de nodos se podran coger varios caminos. Uno de ellos era el de volver a mallar la seccin
transversal de la zona del fluido, el tubo, la soldadura y absorbedor, pero se desech porque se tendra
que volver a validar utilizando por ejemplo simulaciones en tubo simple.
Una solucin mucho ms cmoda era la de reducir es espaciado longitudinal de nodos en la malla.
Actualmente se tiene un nodo cada centmetro, esto es, 183 nodos longitudinales, y se ha pasado a
tener 150 nodos. De esta manera logramos reducir el nmero de nodos hasta los 3 millones. Se deja
para estudios posteriores en anlisis de sensibilidad de la malla.
A continuacin se muestra un resumen de las caractersticas de la malla:
Grid Size
0 2797050 8716997 3115524 32
Memory Usage
Number Used 2797050 8716997 3115524 127 0
Mbytes Used 900 963 166 0 0
Tabla 3.16: Caractersticas de la malla en el modelo completo de nueve tubos

152


3

- Modera de radiacin:
o Sin modelo de radiacin
o Carga Solar: Solar Ray Tracing
- Gravedad: (0,-9.81,0) (m/s
2
)

Frontera Condicin de
contorno
trmica
Solar
Ray
Tracing
Abortividad Transmisividad
IR Visible IR Vis HD
Absorbedor_cara_inferior Wall (acoplada) NO
Absorbedor_cara_inferior-shadow Wall (acoplada) NO
Absorbedor_cara_inferior:086 Wall (acoplada) NO
Absorbedor_cara_inferior:086_shadow Wall (acoplada) NO
Absorbedor_cara_inferior:087
Absorbedor_cara_inferior:087_shadow
Absorbedor_cara_superior Wall (acoplada) SI 0.9 0.9 - - -
Absorbedor_cara_superior-shadow Wall (acoplada) SI 0.9 0.9 - - -
Absorbedor_frontal Wall (adiabtica) NO
Absorbedor_lateral Wall (acoplada) NO
Absorbedor_lateral-shadow Wall (acoplada) NO
Absorbedor_lateral:085 Wall (adiabtica)
Absorbedor_trasera Wall (adiabtica) NO
Air_inf_front Wall (adiabtica) NO
Air_inf_lateral Wall (adiabtica) NO
Air_inf_tras Wall (adiabtica) NO
Air_sup_front Wall (adiabtica) SI 0.9 0.9 - - -
Air_sup_lateral Wall (adiabtica) SI 0.9 0.9 - - -
Air_sup_tras Wall (adiabtica) SI 0.9 0.9 - - -
Cara_exterior_tubo_i Wall (acoplada) NO
Cara_exterior_tubo_i-shadow Wall (acoplada) NO
Cara_exterior_tubo_i:088 Wall (acoplada) NO
Cara_exterior_tubo_i:088-shadow Wall (acoplada) NO
Cara_interior_tubo_ii Wall (acoplada) NO
Cara_exterior_tubo_ii-shadow Wall (acoplada) NO
Cara_exterior_tubo_ii:082 Wall (acoplada) NO
Cara_exterior_tubo_ii:082-shadow Wall (acoplada) NO
Cara_interior_tubo_iii Wall (acoplada) NO
Cara_exterior_tubo_iii-shadow Wall (acoplada) NO
Cara_exterior_tubo_iii:076 Wall (acoplada) NO
Cara_exterior_tubo_iii:076-shadow Wall (acoplada) NO

3
El esquema de discretizacin de primer orden se impuso para que el tiempo de computacin no se
disparase y poder mejorar la convergencia aunque se pierda precisin.
153

Cara_interior_tubo_iv Wall (acoplada) NO
Cara_exterior_tubo_iv-shadow Wall (acoplada) NO
Cara_exterior_tubo_iv:074 Wall (acoplada) NO
Cara_exterior_tubo_iv:074-shadow Wall (acoplada) NO
Cara_interior_tubo_v Wall (acoplada) NO
Cara_exterior_tubo_v-shadow Wall (acoplada) NO
Cara_exterior_tubo_v:001 Wall (acoplada) NO
Cara_exterior_tubo_v:001-shadow Wall (acoplada) NO
Cara_interior_tubo_i Wall (acoplada) NO
Cara_interior_tubo_i-shadow Wall (acoplada) NO
Cara_interior_tubo_ii Wall (acoplada) NO
Cara_interior_tubo_ii-shadow Wall (acoplada) NO
Cara_interior_tubo_iii Wall (acoplada) NO
Cara_interior_tubo_iii-shadow Wall (acoplada) NO
Cara_interior_tubo_iv Wall (acoplada) NO
Cara_interior_tubo_iv-shadow Wall (acoplada) NO
Cara_interior_tubo_v Wall (acoplada) NO
Cara_interior_tubo_v-shadow Wall (acoplada) NO
Cristal Wall(conveccin) SI 0.1 0.1 0.8 0.8 0.8
Entrada_de_fluido_i Velocity-inlet NO
Entrada_de_fluido_ii Velocity-inlet NO
Entrada_de_fluido_iii Velocity-inlet NO
Entrada_de_fluido_iv Velocity-inlet NO
Entrada_de_fluido_v Velocity-inlet NO
Entrada_de_tubera_i Wall (adiabtica) NO
Entrada_de_tubera_ii Wall (adiabtica) NO
Entrada_de_tubera_iii Wall (adiabtica) NO
Entrada_de_tubera_iv Wall (adiabtica) NO
Entrada_de_tubera_v Wall (adiabtica) NO
plano_simtrico Simetra
plano_simtrico:077 Simetra
Salida_de_fluido_i Pressure-outlet NO
Salida_de_fluido_ii Pressure-outlet NO
Salida_de_fluido_iii Pressure-outlet NO
Salida_de_fluido_iv Pressure-outlet NO
Salida_de_fluido_v Pressure-outlet NO
Salida_tubera_i Wall (adiabtica) NO
Salida_tubera_ii Wall (adiabtica) NO
Salida_tubera_iii Wall (adiabtica) NO
Salida_tubera_iv Wall (adiabtica) NO
Salida_tubera_v Wall (adiabtica) NO
Soldadura_frontal_tubo_i Wall (adiabtica) NO
Soldadura_frontal_tubo_ii Wall (adiabtica) NO
Soldadura_frontal_tubo_iii Wall (adiabtica) NO
Soldadura_frontal_tubo_iv Wall (adiabtica) NO
Soldadura_frontal_tubo_v Wall (adiabtica) NO
Soldadura_lado_i_tubo_i Wall (acoplada) NO
Soldadura_lado_i_tubo_i-shadow Wall (acoplada) NO
Soldadura_lado_i_tubo_ii Wall (acoplada) NO
Soldadura_lado_i_tubo_ii-shadow Wall (acoplada) NO
Soldadura_lado_i_tubo_iii Wall (acoplada) NO
154

Soldadura_lado_i_tubo_iii-shadow Wall (acoplada) NO
Soldadura_lado_i_tubo_iv Wall (acoplada) NO
Soldadura_lado_i_tubo_iv-shadow Wall (acoplada) NO
Soldadura_lado_i_tubo_v Wall (acoplada) NO
Soldadura_lado_i_tubo_v-shadow Wall (acoplada) NO
Soldadura_lado_ii_tubo_i Wall (acoplada) NO
Soldadura_lado_ii_tubo_i-shadow Wall (acoplada) NO
Soldadura_lado_ii_tubo_ii Wall (acoplada) NO
Soldadura_lado_ii_tubo_ii-shadow Wall (acoplada) NO
Soldadura_lado_ii_tubo_iii Wall (acoplada) NO
Soldadura_lado_ii_tubo_iii-shadow Wall (acoplada) NO
Soldadura_lado_ii_tubo_iv Wall (acoplada) NO
Soldadura_lado_ii_tubo_iv-shadow Wall (acoplada) NO
Soldadura_lado_ii_tubo_v Wall (acoplada) NO
Soldadura_lado_ii_tubo_v-shadow Wall (acoplada) NO
Soldadura_trasera_tubo_i Wall (adiabtica) NO
Soldadura_trasera_tubo_ii Wall (adiabtica) NO
Soldadura_trasera_tubo_iii Wall (adiabtica) NO
Soldadura_trasera_tubo_iv Wall (adiabtica) NO
Soldadura_trasera_tubo_v Wall (adiabtica) NO
Suelo Wall (adiabtica) NO
Tabla 3. 17: Condiciones de contorno para el modelo completo con nueve tubos
6.2.2. Parmetros de iluminacin y vector de posicin del Sol
En las simulaciones con este modelo se han hecho para el da 1 de Julio, a las 11:00 horas, para una
Latitud de: 37.606 y Longitud 0.97854. Introduciendo estos valores en el Solar Calculator e indicando
el mtodo para el clculo de los parmetros, se devuelve la siguiente informacin:
Fair Weather Conditions:
- Sun Direction Vector: X: -0.106432, Y: 0.920328, Z: -0.376389
- Sunshine Fraction: 1
- Direct Normal Solar Irradiation (at Earth's surface) [W/m^2]: 867.887
- Diffuse Solar Irradiation - vertical surface: [W/m^2]: 142.976
- Diffuse Solar Irradiation - horizontal surface [W/m^2]: 117.454
- Ground Reflected Solar Irradiation - vertical surface [W/m^2]: 91.6194
Al igual que se hizo con el modelo de un solo tubo, a continuacin se presenta el trmino de generacin
de calor que se implementa sobre aquellos elementos que pertenecen a superficies que forman parte
del modelo de carga solar. Los valores mayores son la ID: 67, que es la superficie superior del
absorbedor y la ID:54 es la cubierta de vidrio. Las otras paredes son las superficies verticales que forman
la cavidad de aire entre vidrio y absorbedor.
Internally Scattered Energy [W]: 159.017, Ambient Flux [W/m^2]: 69.0134
Total Integral Energy Source [W]: 8.231533e+002
Compute Time: 15 sec

155

A diferencia del modelo anterior aqu se pretenda obtener la caracterizacin completa del captador y
por ello se han hecho simulaciones tanto para el gasto msico nominal como para otros dos ms,
variando la temperatura de entrada desde los 20 a los 70C. Para esta malla solo se ha probado con una
velocidad de viento de 4 m/s en vistas de la poca variacin que se obtuvo en la eficiencia en el anlisis
de los casos anteriores.
Caso Temperatura de
entrada (C)
Velocidad del Viento
(m/s)
(m/s)
Julio_Te_20_1Q 20 4 0.117204696
Julio_Te_30_1Q 30 4 0.117497862
Julio_Te_40_1Q 40 4 0.117891623
Julio_Te_50_1Q 50 4 0.118380874
Julio_Te_60_1Q 60 4 0.118960782
Julio_Te_20_0.5Q 20 4 0.058602348
Julio_Te_30_0.5Q 30 4 0.058748931
Julio_Te_40_0.5Q 40 4 0.058945811
Julio_Te_50_0.5Q 50 4 0.059190437
Julio_Te_60_0.5Q 60 4 0.059480391
Julio_Te_20_1.5Q 20 4 0.175807044
Julio_Te_30_1.5Q 30 4 0.176246793
Julio_Te_40_1.5Q 40 4 0.176837434
Julio_Te_50_1.5Q 50 4 0.177571312
Julio_Te_60_1.5Q 60 4 0.178441173
Tabla 3.18: Simulaciones realizadas con el modelo completo de nueve tubos


6.3.1. Nmero de Nusselt en el flujo interno

Al igual que se ha hecho para los casos anteriores a continuacin se va a presentar el valor del nmero
de Nusselt local frente a la distancia adimensional y el nmero de Nusselt para flujo totalmente
desarrollado. Sin embargo, a diferencia de las simulaciones anteriores, ahora tenemos diferentes
conductos, por lo que en primer lugar antes de evaluar este parmetro adimensional tenemos que
conocer si existen diferencias importantes entre los tubos.
Se debe tener en cuenta que las condiciones de entrada son las mismas para todos los conductos, y
adems, el modelo de radiacin aporta un trmino de generacin de energa que es comn para todos
los elementos que forman las paredes que son parte del modelo de carga solar. Si representamos la
evolucin del nmero de Nusselt para un caso dado, en el conducto central y en el conducto extremo se
tiene:

156

Figura 3.87: Comparacin del nmero de Nusselt entre el conducto central y el conducto extremo para un caso
con temperatura de entrada de 20C, y gasto igua a 144 kg/h.
Efectivamente, al tener unas condiciones muy parecidas, las diferencias son mnimas. De aqu en
adelante para obtener el nmero de Nusselt en cada caso, simplemente calcularemos el central, a la
vista de que el error cometido ser mnimo. Aun as, hay que remarcar que si la distribucin de caudales
fuera diferente en cada tubo, las diferencias claro est seran muy notables.
Tambin se puede comparar la transmisin de calor entre diferentes casos, por ejemplo para el mismo
gasto msico, (144 kg/h) se obtiene los resultados presentados en la Figura 3. 89. Al igual que ocurra en
las simulaciones con modelos ms simples cuando se introduce el absorbedor, en condiciones de flujo
desarrollado se tiene un nmero de Nusselt ms elevado.
Tambin se observa como a temperaturas mayores, la transmisin de calor mejora, esto se puede
justificar teniendo en cuenta que
, por lo que cuanto mayor es el nmero de

Rayleigh mayor nmero de Nussetl. Si recordamos la expresin para el Rayleigh, este es directamente
proporcional al flujo de calor pero es inversamente proporcional al cuadrado de la viscosidad, por lo que
aunque el flujo de calor disminuya con la tempeatura, la viscosiad tambin lo hace, con lo que el
nmero de Raylegih aumenta.
Adems tambin se ha calculado el nmero de Nusselt para flujo desarrollado, comparndose con la
correlacin que se obtuvo en las simulaciones de tubo con absorbedor (Figura 3.88). Se puede concluir
lo que ya se anticipaba en las simulaciones con radiacin con 1 tubo, y es que el modelo de radiacin
est bien implementado pues se obtienen resultados muy similares a los casos sin carga solar.
1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Tubo central
Tubo extremo
157

1
10
100
0.0001 0.001 0.01 0.1 1
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Resultados Te=20C, 144 kg/h
Resultados Te=50C, 144 kg/h
Petukhov y Polyakov [1988], Conveccin Mixta, Te=20C,144 kg/h
Petukhov y Polyakov [1988], Conveccin Mixta, T=50C,144 kg/h
Figura 3. 89: Nmero de Nusselt local frente a la distancia adimensional. Comparacin con las correlaciones de Petukhov y Polyakov [1988].
1
10
1.00E+05 1.00E+06
N
m
e
r
o

d
e

N
u
s
s
e
l
t
,

N
u

Nussel dessarrollado (9 tubos)
Resultados tubo con absorbedor
Petukhov y Polyakov [1988]
Figura 3.88: Nmero de Nusselt desarrollado frente al nmero de Rayleigh en modelo de 9 tubos.
158

6.3.2. Campo de temperaturas en el absorbedor

A continuacin se va a representar el campo de temperaturas en la pared superior del absorbedor, y se
va a comparar con los resultados que se obtuvieron en el PFC de Maria Jos Martnez Pujante. En este
trabajo se diseo un programa que calcula y representa mediante grficos el campo de temperaturas de
los captadores solares para un gasto msico y una temperatura reducida dados, a partir del filtrado de
datos experimentales del captador solar instalado en el campus Alfonso XII.
En este banco de ensayos se utilizan sondas RTD para medir la temperatura de pared sobre el
absorbedor, dichas sondas se colocan estratgicamente sobre los tubos y sobre las aletas del
absorbedor, suponiendo simetra en el campo de temperaturas. La distribucin de sondas se puede ver
en la siguiente figura:

Figura 3.90: Distribucin de sondas de temperatura en el captador.
Hay que destacar que en el PFC con el que vamos a comparar los resultados se tuvo que hacer frente a
algunas deficiencias en el banco de ensayos (varias sondas se estropearon durante la ejecucin de los
experimentos y no se disponen de estos valores experimentales y existen valores de temperatura en
algunas posiciones necesarios para calcular y representar el campo de temperaturas que no fueron
tenidos en cuenta a la hora de elaborar el diseo del banco experimental), por lo que fue necesario
redefinir las variables asignadas a estas posiciones mediante interpolaciones y extrapolaciones a
partir de aquellos valores de temperatura correspondientes a sondas que midieron
correctamente. Realizadas dichas interpolaciones el programa que desarrollo la autora del PFC realiza la
representacin del campo de temperaturas del tubo central con sus aletas.
Finalmente se ha comparado con las representaciones realizadas para un gasto de 144 kg/h (este
coincide con el caudal que discurre por nuestra simulacin) y para temperaturas reducidas desde 0 a
0.04.
Los campos de temperatura que se representan a continuacin estn delimitados por la zona de
influencia del tubo central, no se representa el captador completo, solo esa parte.

159

160

Por una parte se observa que el rango de temperaturas mxima y mnima para cada temperatura
reducida no difiere excesivamente entre los resultados numricos y los experimentales. La mxima
diferencia se encuentra en T*=0.0 donde el rango de temperaturas numrico es de 23C a 46 C y el
experimental es de 36 C a 52C, sin embargo para T*=0.03, el rango numrico es 52C a 70C y el
experimental 54 C a 72 C.
En cuanto a la distribucin de temperaturas hay que resaltar que al no encontrar en la misma escalara
de colores pude que a primera vista se observen diferencias mayores de las que realmente existen. Por
ejemplo si nos fijamos en el caso de T*=0.03, donde la tonalidad amarilla prcticamente coincide en el
caso experimental y en el numrico, se puede comprobar, como la distribucin es realmente parecida, a
excepcin de la zona de los extremos en las cercanas de la entrada de fluido, donde se puede
comprobar cmo en el caso numrico las temperaturas son menores que las que se alcanzan en el
captador experimental. Lo anterior se repite en todos los casos.
Se puede concluir que existe una elevada concordancia entre los datos de temperatura entre el
captador experimental y los casos numricos teniendo en cuenta que para el caso experimental se ha
realizado interpolaciones para obtener todo el campo y que para el caso numrico las condiciones de
contorno no son exactamente las mismas, sobre todo en lo que se refiere a la soldadura.

Figura 3.91: Comparacin del campo de temperaturas del tubo central para Q=144 kg/h. Datos experimentales (derecha) y datos
numricos (izquierda).
161


Las simulaciones que se han llevado a cabo han sido realizadas para 144 kg/h, 72 kg/h y 216 kg/h. A
continuacin se muestran las curvas de eficiencia obtenidas para cada uno de estos casos.

Figura 3.92: Curvas de eficiencia para 216 kg/h, 144 kg/h y 72 kg/h
Todas las simulaciones tienen los mismos parmetros pticos, lo que se demuestra observando como
todas confluyen en unos valores muy similares en el eje de ordenadas, es decir, queda demostrado que
tienen el mismo factor ptico. La pendiente en todas es la misma, ya que las condiciones de contorno
son iguales. Estos resultados coinciden bastante bien con los datos experimentales de los que se
dispone, en los que se observa que para caudales por encima del nominal la curva de eficiencia se
mantiene constante. Sin embargo a bajos caudales, como el simulado, experimentalmente si se observa
un incremento de la pendiente de la curva, debido a que al reducirse el caudal, el incremento de
temperaturas en el fluido es mayor, el absorbedor opera a mayor temperatura y por tanto las prdidas
trmicas al ambiente son mayores (hay que recordar que en el modelo simulado no se han considerado
las prdidas por conveccin de los laterales y la parte posterior del captador).
Tambin se ha comparado con los resultados experimentales del banco de ensayos del que dispone el
Grupo de investigacin de Mecnica de Fluidos e Ingeniera Trmica en la Universiad, para el caudal de
144 kg/h (Figura 3.93). La diferencia se debe a que el factor ptico no es el mismo ya que en el modelo
experimental =0.883, mientras que en las simulaciones =0.9*0.8=0.72. Hay que tener en cuenta que
el valor en ordenadas se puede relacionar con
, por lo que es algo ms bajo que el producto

absortancia-transmitancia.
Por otra parte la pendiente en las simulaciones es menos pronunciada, ya que si recordamos que la
pendiente se relaciona con a=
, queda claro que si en nuestro caso no tenemos en cuenta las

prdidas que se producen por conveccin en los laterales y la parte inferior la pendiente lgicamente
tiene que ser menor.
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
-0.01 0 0.01 0.02 0.03 0.04 0.05
E
f
i
c
i
e
n
c
i
a

Q=144 kg/h
Q=216 kg/h
Q=72 kg/h
162

Figura 3.93: Comparacin de los resultados numricos para Q=144 kg/h, con los resultados experimentales.

Podemos concluir que la simulacin alcanza un buen ajuste lineal, coherente con los parmetros pticos
definidos y las condiciones de contorno dadas, y que el modelo de carga solar realmente introduce una
energa externa al modelo ajustada a las condiciones radiativas reales.
Por tanto ajustando los parmetros pticos a los que tiene el banco de ensayos experimental y tratando
de implementar alguna correlacin para incluir las prdidas que no se tienen en cuenta, se pude llegar a
obtener unos resultados numricos muy parecidos a los experimentales.

0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
-0.01 0 0.01 0.02 0.03 0.04 0.05
E
f
i
c
i
e
n
c
i
a

Resultados numricos Q=144 kg/h
Eficiencia captador experimental
163

CAPTULO 4: Conclusiones y trabajos futuros
1. Conclusiones del proyecto

- Se ha realizado una sntesis de los procesos de transmisin de calor que intervienen en el
problema de captadores solares para produccin de agua calienta sanitaria, como son la
conveccin forzada y mixta en el interior del conducto, la conveccin natural en la cavidad
interior del colector y incidencia de la radiacin externa sobre el mismo.
- Se ha realizado una bsqueda bibliogrfica para obtener correlaciones experimentales y
analticas acerca de los procesos de transferencia de calor que se dan en el problema estudiado
y poder comparar y validar tanto la metodologa del procesado de resultados como la
exactitud de los resultados numricos obtenidos. Al mismo tiempo tambin se ha investigado
acerca de otros trabajos que tambin utilizan software CFD para la simulacin de la radiacin
solar.
- Se ha desarrollado una metodologa de simulacin numrica con un software de propsito
general que permite estudiar en profundidad diferentes aspectos del funcionamiento de los
captados solares trmicos.
- Se ha seguido un proceso consistente en generar diferentes modelos simples y evolucionarlos
validando paso a paso los resultados hasta obtener un modelo completo similar a los
captadores solares estndar.
- El modelo ms simple creado ha sido un tubo horizontal sin pared al que se aplica un flujo de
calor uniforme a lo largo de su superficie. Con este modelo se ha validado la malla utilizada en
el volumen del fluido calo-portador y la metodologa de procesado de resultados para obtener
el nmero de Nusselt local y promedio para la zona de flujo desarrollado. Los resultados se
correlacionan con la distancia adimensional y con el nmero de Reynolds, respectivamente,
comparando los resultados numricos con las correlaciones experimentales para conveccin
mixta y forzada disponibles en la bibliografa consultada.
- Se ha corroborado por comparacin con el patrn de flujo en el interior del conducto la
aparicin de las recirculaciones tpicas de la conveccin mixta para el rango de nmeros de
Rayleigh para el que se han realizado las simulaciones.
- Se han comparado los resultados del modelo con tubo simple con los obtenidos para un
modelo al que se le ha aadido la pared del conducto, la soldadura y el abosrbedor para un
tubo. Con ello se han podido verificar las diferencias existentes entre la aplicacin de un flujo
uniforme de calor en la pared del tubo y uno aplicado sobre la pared superior del absorbedor.
- Se ha obtenido una correlacin para nmero de Nusselt correspondiente a la regin de flujo
desarrollado para el modelo de un solo tubo con absorbedor donde se comprueba un ligero
aumento de ste parmetro con respecto al caso de tubo simple. Este hecho pues ser debido
al aumento en la intensidad de las recirculaciones el interior del conducto debido a la mayor
asimetra acimutal que se obtiene para la temperatura de pared.
- Se ha estudiado la utilidad de un modelo completo para un solo tubo sin carga solar en el que
se ha aplicado una condicin de generacin interna de calor sobre el absorbedor para predecir
la curva de eficiencia completa del captador con un coste computacional muy reducido. En ste
se han podido estudiar los procesos de conveccin natural en el aire entre el absorbedor y la
cubierta de vidrio, la distribucin de temperaturas sobre el absorbedor y las prdidas por
conveccin al ambiente a travs de la cubierta de vidrio, sin necesidad de implementar un
modelo de radiacin.
- Se han comparado los resultados de la conveccin natural del modelo anterior con una
correlacin para el nmero de Nusselt en cavidades con paredes isotermas, obtenindose
CAPTULO 4: Conclusiones y trabajos futuros
164

buenos resultados , al tiempo que tambin se ha comparado el nmero de Rayleigh crtico con
el nmero de Rayleigh de las simulaciones verificndose que para valores en torno al crtico se
obtiene un campo de temperaturas y unas recirculaciones estacionarias, sin embargo, para
temperaturas de entrada de fluido elevadas (50-60C) el comportamiento es transitorio.
- Se han estudiado en profundidad los diferentes modelos de radiacin implementados en
FLUENT y se han definidos las ventajas e inconvenientes de cada uno de estos modelos con
objeto de aplicar el ms adecuado al problema en estudio. Finalmente se ha escogido como
modelo de carga solar el modelo Solar Ray Tracing por su facilidad de implementacin y su
bajo coste computacional.
- Se ha creado un modelo completo para un solo tubo con carga solar, obtenindose valores para
el nmero de Nusselt de flujo desarrollado muy similares a los del modelo sin radiacin, lo que
indica que el modelo de carga solar ha sido adecuadamente implementado y proporciona
resultados coherentes.
- Para el modelo anterior se han calculado las curvas de eficiencia trmica del captador,
concluyndose que el factor ptico obtenido mediante simulacin numrica es muy parecido al
obtenido experimentalmente en banco de ensayos para un captador estndar. Sin embargo, el
coeficiente de prdidas, es decir, la pendiente de las curva es menor. Este hecho puede ser
debido a que no se han tenido en cuenta las prdidas de conveccin a travs de los laterales y
la parte posterior del captador.
- Por ltimo se ha creado una malla simtrica para modelar el captador completo de nueve tubos
estndar. Se han obtenido buenos resultados para el nmero de Nusselt de flujo desarrollado,
con lo que se remarca lo que ya se haba concluido en el modelo con carga solar para un solo
tubo.
- Tambin se han obtenido las curvas de eficiencia trmica con el modelo de 9 tubos para tres
gastos msicos diferentes, el nominal de 144kg/h, 72 kg/h y 216 kg/h que representan el 50% y
el 150% respecto del gasto msico nominal. Las curvas de eficiencia trmica obtenidas
numricamente para estos diferentes gastos msicos son muy parecidas. Sin embargo,
experimentalmente se observa que si bien el factor ptico es el mismo, el coeficiente de
prdidas es mayor para el caso de menor gasto msico. Esto es debido a que a menos gasto
msico, mayor salto trmico del fluido calo-portador y mayor temperatura del absorbedor, por
lo que las prdidas por conveccin aumentan. Como en el estudio numrico no se han
modelado las prdidas por conveccin por los laterales y fondo del captador, este efecto no se
aprecia.
- El modelo desarrollado permite estudiar en detalle el campo de temperaturas del absorbedor y
analizar la influencia sobre el funcionamiento y la eficiencia trmica del captador la
modificacin de las caractersticas de operacin y geomtricas, pudiendo estudiarse la
influencia de otros parmetros como el ngulo de inclinacin, la influencia de fluidos calo-
portadores formados por mezclas anticongelantes, la influencia de la geometra de la
soldadura, la seccin transversal de los tubos, etc.

165

2. Trabajos futuros

- Dados los resultados que se han obtenido sera muy conveniente introducir para el modelo
completo de 9 tubos la absortividad y la transmisividad del captador ensayado
experimentalmente, as como, introducir correlaciones de prdidas por conveccin para los
laterales y la parte inferior.
- Realizar un estudio de sensibilidad de la malla para el modelo completo de 9 tubos para reducir
el coste computacional sin perder exactitud en los resultados con objeto de que se puede
completar la malla con nuevos elementos y modelos de radiacin ms completos.
- Implementar el modelo S2S con el objeto de verificar si existen diferencias notables al tener en
cuenta la emisin desde al superficies internas del modelo.
- Estudiar la evolucin del nmero de Nusselt y la curva de eficiencia trmica para diferentes
ngulos de inclinacin del captador.
- Insertar al modelo de nueve tubos los colectores de entrada y salida para estudiar la
problemtica de la distribucin de caudales y obtener un modelo de prdidas de presin en el
captador.
- Estudiar la posibilidad de insertar elementos en el interior de los conductos como muelles o
lminas en espiral para mejorar la transferencia de calor.
- Estudiar la influencia de diferentes parmetros de inters como:
o Diseo y tamao de la soldadura absorbedor-tubo
o Posicin del tubo sobre el absorbedor
o Influencia de insertar elementos entre el absorbedor y la cubierta de vidrio para
suprimir la conveccin natural.
o Etc.

Bibliografa
166

Bibliografa
Bibliografa bsica

[1] G.L. Morrision, M. Behnia, M. Cook, N.K. Groenhout, D.R Mills. Optimal Design of advanced
solar water heaters. Proceedings of the Sixth Australasian Natural Convection workshop,
Murdoch University, 1999.
[2] M. Selmi, A. Marafi, M. Al-Khawaja. Performance prediction of a flat plate solar collector
using CFD software. Proceedings of IMEC2004 International Mechanical Engineering
Conference, Kuwait, December 5-8, 2004
[3] M. Selmi, A. Marafi, M. Al-Khawaja. CFD Simulation of fluid flor and heat transfer in a flan
plate solar collector. Proceedings of the 2004 World Renewable Energy Congress. Boulder, CO,
August 28- September 3, 2004
[4] D. Dovi. S. vai, A. Galovi. Estimating heat losses in solar collectors by IR thermography
and numerical simulations. Proceedings of 9
th
International Conference on Quantitative
InfraRed Thermography. July 2-5, Krakow-Poland 2008
[5] M. Selmi, M. Al-Khawaja, A. Marafia. Validation of CFD simulation for flan plate solar energy
collector. Renewable Energy 33 (2008) 383-387
Otras referencias:

[1] Hobbi, A., Siddiqui K., Experimental study on the effect of heat transfer enhancement
devices in flat-plate solar collectors, International Journal of Heat and Mass Transfer, Vol. 52,
Issues 19-20, 2009, pp 4650-4658.
[2] Prasad, B.N., and Saini J.S., Effect of artificial roughness on heat transfer and friction factor
in a solar air heater, Solar Energy, Vol. 41, Issue 6, 1988, pp 555-560.
[3] Gupta, D., Solanki, S. C. and Saini J. S., Thermohydraulic performance of solar air heaters
with roughened absorber plates, Solar Energy, Vol. 61, Issue 1, 1997, pp 33-42.
[4] Bhagoria, J. L., Saini, J. S., and Solanki S. C., Heat transfer coefficient and friction factor
correlations for rectangular solar air heater duct having transverse wedge shaped rib
roughness on the absorber plate, Renewable Energy, Vol. 25, Issue 3, 2002, pp 341-369.
[5] Momin, A-M. E., Saini, J. S., and Solanki, S. C., Heat transfer and friction in solar air heater
duct with V-shaped rib roughness on absorber plate, International Journal of Heat and Mass
Transfer, Vol. 45, Issue 16, 2002, Pages 3383-3396
[6] Sahu, M.M., and Bhagoria, J.L., Augmentation of heat transfer coefficient by using 90
broken transverse ribs on absorber plate of solar air heater, Renewable Energy, Vol. 30, Issue
13, 2005, pp 2057-2073.
167

[7] Jaurker, A.R., Saini, J.S., and Gandhi, B.K., Heat transfer and friction characteristics of
rectangular solar air heater duct using rib-grooved artificial roughness, Solar Energy, Vol. 80,
Issue 8, 2006, pp 895-907.
[8] Mittal, M.K., Varun, Saini, R.P., Singal, S.K., Effective efficiency of solar air heaters having
different types of roughness elements on the absorber plate, Energy, Vol. 32, Issue 5, May
2007, pp 739-745.
[9] Varun, Saini, R.P., Singal, S.K., A review on roughness geometry used in solar air heaters,
Solar Energy, Vol. 81, Issue 11, 2007, pp 1340-1350
[10] Kumar, A., Bhagoria, J.L., R.M. Sarviya, Heat transfer and friction correlations for artificially
roughened solar air heater duct with discrete W-shaped ribs Energy Conversion and
Management, Vol. 50, Issue 8, 2009, pp 2106-2117.
[11] Corts, A., Piacentini, R., Improvement of the efficiency of a bare solar collector by means
of turbulence promoters, Applied Energy, Vol. 36, Issue 4, 1990, pp 253-261
[12] Sparrow, E. M. & Tao, W. Q. Enhanced heat transfer in a flat rectangular duct with
streamwise-periodic disturbances at one principal wall. J. Heat Transfer, 105,1983, p. 851.
[13] Saini, R. P., Saini, J. S., Heat transfer and friction factor correlations for artificially
roughened ducts with expanded metal mesh as roughness element, International Journal of
Heat and Mass Transfer, Vol. 40, Issue 4, 1997, pp 973-986.
[14] Kumar, A., Prasad, B. N., Investigation of twisted tape inserted solar water heaters-heat
transfer, friction factor and thermal performance results, Renewable Energy, Vol. 19, Issue 3,
2000, pp 379-398.
[15] Jaisankar, S., Radhakrishnan, T.K., Sheeba K.N., Experimental studies on heat transfer and
friction factor characteristics of forced circulation solar water heater system fitted with helical
twisted tapes, Solar Energy, Vol. 83, Issue 11, 2009, pp 1943-1952.
[16] Varun, S., Thermal performance optimization of a flat plate solar air heater using genetic
algorithm, Applied Energy, In Press, Corrected Proof, Available online 20 November 2009.
[17] Klein S.A., Calculation of flat plate loss coefficients, Solar Energy 17, 1975, pp. 17:79.
[18] Hottel HC, Woetz BB. Transaction of the American Society of Mechanical Engineering
1942, Vol. 64, pp 91.
[19] Duffie, J.A., Beckman, W.A., Solar Engineering of thermal processes, second ed. Wiley
Interscience, 1991, New York, pp 268-272.
[20] Lund, K.O., General Thermal analysis of serpentine-flow flat-plate solar collector
absorbers. Solar Energy 42, 1989, 133-142.
[21] Eisenmann W., Vajen K., Ackermann H.. On the correlations between collector efficiency
factor and material content of parallel flow flat-plate solar collectors. Solar Energy 76, 2004,
pp. 381-387.

168

[22] Garca, A. Tesis Doctoral: Mejora de la transmisin de calor en tubos circulares lisos
mediante la insercin de muelles en espiral, 2007, UPCT.
[23] Tam L.M. and Ghejar A.J. Transitional heat transfer in plain horizontal tubes. Heat Tranfer
Engineering, 27, 2006, pp. 23-28.

169

Anexo I: Propiedades de los materiales
1. Propiedades del agua
Las propiedades del agua en funcin de la temperatura son:
T (C) T (K)
Densidad
(kg/m
3
)
Viscosidad
Molecular
(Pa.s)
Viscosidad
Cinemtica
(m
2
/s)
Calor
esp.
(J/kgK)
Prandtl
Cond.
Trmica
(W/mK)
0 273.15 1000 0.001788 0.000001788 4217 13.67 0.551572495
10 283.15 1000 0.001307 0.000001307 4192 9.47 0.578557973
20 293.15 998 0.001003 0.000001005 4182 7.01 0.598366049
30 303.15 996 0.000799 0.000000802 4178 5.43 0.614773849
40 313.15 992 0.000657 0.000000662 4178 4.34 0.632476037
50 323.15 988 0.000548 0.000000555 4180 3.56 0.643438202
60 333.15 983 0.000467 0.000000475 4184 2.99 0.653487625
70 343.15 978 0.000405 0.000000414 4189 2.56 0.662712891
80 353.15 972 0.000355 0.000000365 4196 2.23 0.667973094
90 363.15 965 0.000316 0.000000327 4205 1.96 0.67794898
100 373.15 958 0.000283 0.000000295 4216 1.75 0.681787429
120 393.15 943 0.00023 0.000000249 4248 1.45 0.67382069
2. Propiedades del aire
Las propiedades del agua en funcin de la temperatura son:
T (C) T (K)
Densidad
(kg/m
3
)
Viscosidad
Molecular
(Pa.s)
Difusividad
trmica
(m
2
/s)
Calor esp.
(J/kgK)
Prandtl
Cond.
Trmica
(W/mK)
0 273.15 1.292 1.72E-05 1.86E-05 1006 0.72 0.0242
20 293.15 1.204 1.81E-05 2.12E-05 1006 0.71 0.0257
40 313.15 1.127 1.90E-05 2.40E-05 1007 0.7 0.0272
60 333.15 1.059 1.99E-05 2.69E-05 1008 0.7 0.0287
80 353.15 0.999 2.09E-05 3.00E-05 1010 0.7 0.0302
100 373.15 0.946 2.18E-05 3.32E-05 1012 0.69 0.0318
120 393.15 0.898 2.27E-05 3.66E-05 1014 0.69 0.0333
140 413.15 0.854 2.34E-05 3.98E-05 1016 0.69 0.0345
160 433.15 0.815 2.42E-05 4.32E-05 1019 0.69 0.0359
180 453.15 0.779 2.50E-05 4.67E-05 1022 0.69 0.0372
200 473.15 0.746 2.57E-05 5.05E-05 1025 0.68 0.0386
220 493.15 0.715 2.64E-05 5.43E-05 1028 0.68 0.0399
240 513.15 0.688 2.72E-05 5.80E-05 1032 0.68 0.0412
260 533.15 0.662 2.79E-05 6.20E-05 1036 0.68 0.0425
280 553.15 0.638 2.86E-05 6.59E-05 1040 0.68 0.0437
300 573.15 0.616 2.93E-05 6.99E-05 1045 0.68 0.045

Anexo I: Propiedades de los materiales
170

3. Propiedades de materiales slidos

En todas las simulaciones las propiedades de los materiales slidos son constantes e iguales a:

Densidad (kg/m
3
) Calor esp. (J/kgK) Conductividad Trmica (W/mK)
Aluminio 2719 871 202.4
Cobre 8978 381 387.6
Vidrio 2700 840 0.78

Transferencia de Calor Con Colectores Planos

Cargado por

Copyright:

Formatos disponibles

Transferencia de Calor Con Colectores Planos

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Transferencia de Calor Con Colectores Planos

Cargado por

Copyright:

Formatos disponibles

Anlisis del mecanismo de

, la potencia absorbida por el colector, S, y las

es el factor de forma, que viene dado por:

viene dado por:

son las constantes de prdidas.

, la longitud hidrodinmica de entrada se puede obtener a partir de una

son cero en todas partes.

o un flujo de calor uniforme (

, mientras que lo inverso es cierto para

, se incrementa linealmente con la distancia axial a lo largo del tubo. El valor de

puede calcularse haciendo un balance energtico desde la entrada a la seccin de intercambio

es la temperatura media del fluido a la entrada de la seccin de intercambio trmico (x=0). La

vara linealmente con la distancia axial . Para

, el nmero de Nusselt en conveccin forzada

, el nmero de Nusselt en conveccin forzada puede

, es la distancia axial necesaria para conseguir un valor de nmero de

. Para calcular el nmero de Nusselt global en un tubo de longitud l se

es ahora variable en la seccin,

la temperatura media de la pared en la seccin transversal.

que aparece en las ecuaciones de cantidad de movimiento indica la importancia

entre los dos

, es el nmero de Nusselt en conveccin forzada pura, ecuacin (20),

; aqu el nmero de Nusselt se hace constante e igual al conocido valor de

. Por otro lado, la ecuacin (2.77) indica que en la regn plenamente

. De aqu el gradiente de presin en la direccin

es una velocidad de referencia arbitraria, las ecuaciones de

, no es conveniente en su forma acutal. Se acostumbra por

y la pared que limita al fluido por

, es el flujo de calor que no proviene de la radiacin de la superficie interior que se encuentra a

, es el flujo de calor radiante que proviene tambin de esa cara y

cuando el fluido est en estado estacionario y su comportamiento

se deben determinar por un anlisis combinado de

para valores altos de Ra.

son funcin de Pr, y se encuentran tabulados en la tabla siguiente:

, que es el que ms se acerca a las condiciones que

. En la Figura 2.27 vemos que el ngulo plano diferencial se define

sobre la esfera al cuadrado del radio de la esfera. En consecuencia,

es normal a la direccin (, ), y segn se muestra en la Figura 2.28, se puede representar

para una superficie esfrica. Por consiguiente,

de la radiacin emitida. Definimos

, como la razn a la que se emite energa radiante a la longitud de onda en la

, es la razn a la que la radiacin de longitud de onda sale de

. Al reacomodar la ecuacin anterior, se sique que:

, el flujo de calor asociado con la emisin

, como se muestra en la Figura 2.30, es:

se define como la intensidad a la que se emite radiacin

se basa en el rea proyectada. El trmino

, es la rapidez a la que se emite radiacin por unidad de

. Entonces se tiene que:

. Esta cantidad se define como la proporcin en que la energa

se define como la velocidad a la que la radiaicn de longitud de onda incide

se define en trminos del rea proyectada. Si la

representa la razn a la que incide la radiacin por unidad de rea desde

representa la rapidez a la que la radiacin de longitud de onda unitaria alrededor .

, por una expresin de la forma:

asociada con todo el espectro es:

de una superficie a una temperatura T

, que representa un promedio

, se puede definir como:

, de una superficie se define como la fraccin de la intensidad espectral incidente en la