Conversión
y conservación
energética
Programa:
Maestría en Electromecánica
Cohorte: 2019
Nombre del docente que imparte el módulo:
PhD. Enrique Torres Tamayo
Alumno:
Milton Xavier Remache Sasig
Latacunga - Ecuador
Fecha en que se comienza el módulo: 26-05-2020
ix
1. Título del proyecto
Dimencionamiento de una central termosolar de colectores parabólicos para
producir 20-MW.
2. Introducción.
El Ecuador tiene la ventaja al tener una geografía variada y de radiación solar alta,
es decir es conveniente la utilización de la energía solar para el desarrollo del
cálculo de una central termoeléctrica, los colectores parabólicos nos ayudaran a
concentrar la radiación solar, la cual mediante la trasferencia de calor permitirá el
calentamiento del líquido refrigerante o fluido caloportador (aceite térmico)
permitiendo así tener un fluido a alta temperatura, este fluido será trasportado por
medio de una bomba o un hidroneumática al intercambiador de color de flujo
contrario permitirá por un lado ingresara el aceite térmico y por el otro agua
precalentada, el agua se evaporara en el intercambiador de calor dirigiéndose a la
turbina esta a su vez permite el paso al generador, el agua que sale del generador
ingresa al condensador donde se da el paso de vapor a alta presión y alta temperatura
a líquido a baja temperatura y baja presión, el presente trabajo está orientado a
desarrollar el cálculo de una planta de generación de energía termosolar y el mismo
puede servir de base para proponer esta alternativa de generación eléctrica en el
país.
3. Problema de la investigación.
En la actualidad a nivel mundial se conoce sobre la utilización de energías
renovables y la aplicación en los distintos campos de la ingeniería, en nuestro país
el desarrollo de las energías renovables tienen mucha importación ya que en la
mayoría de instituciones de educación superior implementan Bosques solares,
parques eólicos, entre otros sistemas de generación, sin embargo el desarrollo de
proyectos para el aprovechamiento de la energía solar por medio de sistemas de
concentradores parabolicos para generar energía termosolar es muy limitada.
2
El cálculo de los sistemas de concentración de la energía solar térmica en nuestro
país es muy bajo, el presente calculo es muy importante en Ecuador, siendo
necesario para la investigación y posterior análisis de los datos obtenidos.
3.1. Formulación del problema
¿El cálculo para el dimensionamiento de una central termosolar de colectores
parabólicos para producir 20-MW ayudara para la generación de un modelo
aproximado de central de generación eléctrica termo solar?
4. Objetivos.
4.1. General.
Determinar el dimensionamiento de una central termosolar de colectores
parabólicos.
4.2. Específicos.
• Desarrollar el cálculo y las características que permitan el ensamblaje del
sistema de los colectores solares parabolicos.
4.3. Sistemas de tareas en relación a los objetivos específicos.
Para realizar la investigación se desarrollará el cálculo del colector solar, mediante
el conocimiento de estos se procederá al cálculo mediante un software Mathcad
prime 6.0.
3
5. Descripción del objeto de estudio.
Él estudia que realizaremos se basa en la generación de energía eléctrica utilizando
el ciclo termodinámico de Rankine, los colectores solares parabólicos recolectan o
concentran la energía térmica la cual es pasada mediante la trasferencia de calor al
aceite térmico con las propiedades de lograr grandes temperaturas.
El primer circuito (Bloque solar), será el del aceite térmico el cual está compuesto
de tubería metálica que compone el captador, esta tubería a demás servirá como
fuente de almacenaje del aceite, débenos de tener presente que los fluidos
refrigerantes tienen puntos de trabajo tanto de congelación como de evaporación,
(este último nos da un rango de trabajo la alta temperatura de 400°C al pasar esta
temperatura perderá sus propiedades), el punto de congelación también es
importante ya que podríamos taponar las tuberías siendo un daño serio, los puntos
de congelación dependen de los fluidos en los aceites varían desde los 10 grados
hasta los 200°C aproximadamente, este circuito contara:
• Bomba de fluido térmico: Es la que nos permite tener una presión constante
en el circuito del aceite refrigerante.
• Campo de colectores: Es el área donde se dispondrá o se instalaran los
sistemas de captación solar.
• Calentador de respaldo: Este puede ser un caldero nos servirá para mantener
la temperatura por arriba del punto de congelación del fluido refrigerante o
aceite térmico.
• Intercambiador de sales, (tanque caliente y tanque frio): El intercambiador
de sales realiza una mezcla en proporciones adecuadas de los minerales para
evitar la corrosión de la tubería además en algunos sistemas sirve como
anticongelante evitando la congelación del fluido.
• Tanque de expansión (hidroneumático): Este nos permite tener la presión
constante en el sistema de captación.
El segundo circuito principal (Bloque de potencia): Es el circuito del agua que
también trabaja como fluido intercambiador de calor, consta de la tubería de
4
circulación, tenemos presenta que el punto de trabajo del agua es, punto de
ebullición y el punto de congelamiento, (aproximadamente 100°C para el punto de
ebullición y 0°C para el punto de congelación), este circuito está integrado por:
• Condensador: El flujo de agua sale del generador a presión baja y
temperatura baja, pero teniendo que disminuirlas aún más.
• Bomba de condensado: Esta recibe el agua a baja temperatura y presión que
mediante el desairado sube las presiones a una más alta de 120 psi y la
temperatura permanece constante.
• Desgasificador: Recibe el aire particulado imposible de comprimir y lo
desecha en forma de gas y partículas mediante esta la presión y temperatura
suben.
• Turbina: Recibe el vapor de agua a alta temperatura y presión para así mover
las aletas.
• Generador: luego de que pase el vapor de agua por la turbina de dirige al
recalentador volviendo de nuevo a una presión y temperatura alta.
• Torre de enfriamiento, (bomba de agua de enfriamiento).
Los elementos que realizan la transferencia de calor en los dos circuitos cerrados
son:
• Precalentador.
• Evaporador.
• Sobrecalentador.
• Recalentador.
Los cuatro elementos mencionados son intercambiadores de calor dependiendo de
la configuración se logrará las temperaturas en los rangos deseados.
5
6. Esquema y parámetros de trabajo.
Gráfico 1: Esquema de generación de energía con concentradores cilíndricos parabólicos.
429 °C =703 °K
19 °C =293 °K
39 °C =313 °K
Fuente: Fundamentals of Heat and mass transfer. Pag 732.
6
Gráfico 2: Esquema de una central de generación termoeléctrica ideal.
429 °C =703 °K
19 °C =293 °K
39 °C =313 °K
Fuente: [Link]/bitstream/20.500.12423/1083/3/TL_GomezChavarryBilly.pdf.
7
Mediante los esquemas de la parte superiore describiremos nuestros parámetros de
trabajo la planta de generación consta de las siguientes áreas:
• Campo solar
• Sistema de transferencia de calor
• Bloque de potencia
4.1. Parámetros de partida para el dimensionamiento del campo solar.
Como parámetro de inicio de diseño tenemos que la potencia será de 20 MW, las
temperaturas de ingreso del aceite térmico serán de 400°C y la de salida de 450°C,
tenemos una diferencia de temperatura de 50°C, el flujo de radiación solar incidente
del área de instalación será de 20. 000 W/m², el diámetro de la tubería del
concentrador será de 70 mm, el flujo masico de flujo térmico es de 2.5 kg/s.
Conociendo que nuestro medio de transferencia de calor será el aceite térmico
describiremos sus propiedades donde la densidad es 700 kg/m³, la conductividad
térmica es de 0.078 W/m*K, la capacidad calorífica es de 2590 J/kg*k, la viscosidad
dinámica es de 0.15*10ˉ³ N*s/ m², mediante los datos presentados calcularemos:
• Calcular la longitud de la tubería en cada concentrador.
• Cálculo del flujo de calor transferido.
• Cálculo de la temperatura de la superficie del concentrador en la salida
Número de Reynolds.
• Cálculo del número de prandtl.
• Cálculo del número de Nusselt.
• Cálculo de la temperatura de la superficie del concentrador en la salida.
• Cantidad mínima de concentradores para general una potencia de 20MW de
electricidad, eficiencia máxima en la producción de vapor sobrecalentado
(Eficiencia de Carnot).
• La cantidad de energía mínima requerida.
• Cantidad mínima de concentradores solares parabólicos.
8
5. Desarrollo del cálculo.
Datos iniciales del sistema
Flujo de radiación solar incidente
Diámetro de la tubería en el concentrador
Flujo másico de aceite térmico
Temperatura media de entrada del aceite térmico en el colector
Temperatura media de salida del aceite térmico en el colector
Potencia generada
Propiedades del aceite térmico
Densidad de flujo Viscosidad dinámica
Conductividad Térmica
Capacidad calorífica
Gráfico 3: Datos de la tubería.
Fuente: Fundamentals of Heat and mass transfer. Pag 733.
9
1. Cálculo de la longitud de la tubería en cada concentrador
Para calcular la longitud de la tubería de cada concentrador multiplicamos el flujo
masico de aceite térmico por la capacidad calorífica del aceite esta lo dividimos
para el perímetro de la tubería que es la superficie del contacto con el fluido, el flujo
de radiación solar incidente esta asu vez le multiplicamos por la variación de la
temperatura de salida menos la temperatura de ingreso dandomos como resultado
la longitud de la tubería siendo 73.6 metros.
2. Cálculo del flujo de calor transferido
Para el calculo de flujo de calor trasferido multiplicamos el flujo de radiación solar
incidente por el perímetro de la tubería por la longitud de la tubería del concentrador
dándonos como resultado 0.324 MW.
3. Cálculo de la temperatura de la superficie del concentrador en la salida
Número de Reynolds
El número de Reynols es adimensional este es igual la densidad de flujo por la
velocidad del fluido por el diámetro de la tubería del líquido sobre la viscosidad
dinámica.
Para el flujo en un tubo circular (𝐴 = 𝜋𝐷2 /4), el número de Reynolds se reduce a:
De acuerdo a la teoría revisada la diferencia entre el régimen laminar y turbulento
se obtiene por el número de Reynolds, siendo para el régimen R<2000, el régimen
10
laminar y R>4000 es turbulento, la zona entre 2000 y 4000 se le denomina Régimen
de Transición, donde el flujo es inestable pasando de régimen turbulento a laminar
o viceversa.
4 Cálculo del número de prandtl
El número de Prandtl es adimensional y describe el espesor relativo de las capas
límite de la velocidad del flujo y térmica, el número de Prandtl es igual a la
difusividad molecular de la cantidad del movimiento sobre la difusividad molecular
del calor, que es igual también a la viscosidad dinámica por la capacidad calorífica
específico del fluido sobre la conductividad térmica.
Pr va desde 0.01 para metales líquidos hasta 100 000 para aceites pesados. Es decir
que en los primeros el calor se difunde con mucha rapidez (Pr<<1) y en los
segundos con mucha más lentitud (Pr>>1).
Como referencia se menciona que el agua tiene un numero de Pr 10 y alrededor de
1 en gases, es decir que tenemos un aceite térmico adecuado para el sistema ya que
tenemos un numero de Prandtl de 4.98.
5. Cálculo del número de Nusselt
El coeficiente del número de Nusselt sale de comparar la trasferencia de calor por
conducción con la de convección.
Tasa de trasferencia de calor por conducción: Cuando no hay movimiento en las
moléculas y pasa el calor de una a otra molécula.
Tasa de trasferencia de calor por convección: Cuando ay movimiento de moléculas
y pasado una a una por el movimiento.
11
Por lo tanto, el NuD representa el mejoramiento de la trasferencia de calor a través
de una capa de fluido como resultado de la convección con la conducción, si el
número de NuD =1 significa que la conducción es igual a la convección es decir ay
una conducción pura.
Calculamos el coeficiente de trasferencia de calor por convección.
6. Cálculo de la temperatura de la superficie del concentrador en la salida
La temperatura de la superficie a la salida del concentrador se encuentra dividiendo
el flujo de radiación solar incidente sobre el coeficiente de trasferencia de calor mas
la temperatura media a la salida del aceite termico. Teniendo como resultado que a
la salida del concentrador una temperatura de 466.9 °C.
7. Cantidad mínima de concentradores para general una potencia de 20MW de
electricidad, eficiencia máxima en la producción de vapor sobrecalentado
(Eficiencia de Carnot)
Temperatura de entrada del agua de enfriamiento en el condensador
Diferencia de temperatura en el condensador
Temperatura mínima del agua de enfriamiento a la salida del
condensador
12
Rendimiento del ciclo de vapor: Este lo obtenemos dividiendo la temperatura
mínima del agua sobre la temperatura máxima obteniento un valor de 0.555 es decir
tendremos que aumentar la capacidad en un 45%.
8. La cantidad de energía mínima requerida
La capacidad inicial es de 20 MW para cumplir con la generación y cubrir el índice
de seguridad tendremos que instalar 16 MW más adicionales.
9. Cantidad mínima de concentradores solares
parabólicos
En total tendremos que instalar 111 concentradores.
6. Conclusiones.
1. Mediante los datos analizados debemos mencionar que trabajaremos en un
sistema con flujo turbulento que tiene buena capacidad de de transferencia
de calor entre el aceite térmico y los captadores solares.
2. Para cumplir con el rendimiento del sistema de generación eléctrica tanto
en el día despejado, los días nublados y la noche tendemos que instalar 16
MW adicionales.
3. Debemos tener cuidado con el aceite térmico ya que a partir de 450°C pierde
sus propiedades térmicas, esto debido que a la salida de los concentradores
de calor la temperatura llegaría a 466.9°C.
4. En los procesos de calculo en la parte superior se describe para que sirve
cada uno de los paramatros anallizados.
13
Datos ejercicio 1.
Se desea determinar el comportamiento de las pérdidas de energía calorífica desde la
superficie exterior de un horno conociendo los siguientes datos.
Cálculo del flujo de calor por convección forzada por el exterior de un conducto
en superficies planas.
Temperatura del fluido
Temperatura de la superficie
Velocidad del aire atmosférico
Longitud de la superficie, Lc
Cálculo de la temperatura pelicular.
Propiedades del fluido.
Área de la superficie Número de Prandtl
Coeficiente cinemático de Viscosidad
Coeficiente de conductividad térmica del fluido
14
Cálculo del número de Reynolds
Cálculo del número de Nusselt
Cálculo del coeficiente de transferencia de calor por convección
15
Cálculo del flujo de calor transferido
16
Comportamiento del número de Reynolds en función de la longitud característica.
Comportamiento del coeficiente de transferencia de calor en función del número
de Reynolds.
17
Comportamiento del flujo de calor transferido en función del número de
Reynolds.
18
Ejercicio 2.
Cálculo del flujo de calor por convección forzada por el exterior de un conducto
en superficies planas
Se desea determinar el comportamiento de las pérdidas de energía calorífica desde la
superficie exterior de un horno conociendo los siguientes datos.
Temperatura de la superficie
Velocidad del aire atmosférico
Temperatura del fluido
Longitud de la superficie
Área de la superficie
Cálculo de la temperatura pelicular
Propiedades del fluido
Coeficiente cinemático de Viscosidad
Número de Prandtl
Coeficiente de conductividad térmica del fluido
Cálculo del número de Reynolds
El valor del Reynolds crítico en superficies planas es 500000
19
Cálculo del número de Nusselt
Cálculo del coeficiente de transferencia de calor por convección
Cálculo del flujo de calor transferido
20
Ejercicio 3.
Disipadores de calor (Aletas)
Datos del problema
Profundidad de las aletas Eficiencia de las aletas
Altura de las aletas
Número de aletas Velocidad del aire atmosférico
Espesor de las aletas
Longitud de las aletas Temperatura del fluido
Temperatura de la
superficie
Cálculo de la temperatura pelicular
Propiedades del fluido
Coeficiente cinemático de
Viscosidad
Coeficiente de conductividad térmica del fluido
Número de Prandtl
21
Cálculo del número de Reynolds
El valor del Reynolds crítico en superficies planas es 500000
Cálculo del coeficiente de transferencia de calor por convección
Cálculo de la longitud característica de las aletas
Cálculo del Área superficial de las aletas
Cálculo del Área total de las aletas
22
Cálculo del flujo de calor total disipado por las aletas
Cálculo del flujo de calor total disipado por la superficie sin aletas
Cálculo del incremento del flujo de calor disipado
23
Ejercicio 4.
Un intercambiador de calor compacto tiene hileras de placas metálicas separadas por N
aperturas de ancho "a". El aceite y el agua fluyen por aperturas separadas N/2 en el
intercambiador de calor, en direcciones opuestas alternando las aperturas. El aceite y el
agua entran a temperaturas de 100 y 30°C, respectivamente. La temperatura de salida del
aceite es de 60°C.
Datos Iniciales
Las propiedades del aceite se escogen de la Tabla A.5 (Incropera, 2013)
24
Interpolación Lineal para las propiedades del aceite y el agua
y=
Propiedades del agua
25
El número de Nusselt se calcula usando la tabla 8.1, conociendo que en los intercambiadores de
calor de flujos cruzados el régimen es laminar y "a" es mucho menor que "L".
Para determinar el área del intercambiador de calor es necesario establecer la relación entre el
número de aperturas "N" y la longitud "L". El ancho de la apertura a=L/N.
Área del intercambiador de calor Diámetro Hidráulico
A=L²(N-1)
Cálculo del coeficiente de transferencia de calor por convección en función del número de apertura y la
longitud
26
Cálculo del flujo de calor transferido
Cálculo de la temperatura de salida del agua
Cálculo de la temperatura media logarítmica
Cálculo de la longitud del intercambiador de calor
Cálculo de la longitud del intercambiador de calor
Cálculo del ancho entre las placas
Cálculo del diámetro hidráulico
Cálculo de la velocidad media del agua
Cálculo del número de Reynolds para el agua
Régimen Laminar
Cálculo de la velocidad media del aceite
27
Cálculo del número de Reynolds para el aceite
Régimen Laminar
Cálculo de las pérdidas de presión para el agua
Cálculo de las pérdidas de presión para el aceite
Cálculo del coeficiente de convección para el agua
= 7.54kN/2L =7.54 X 0.625 W/m·K X N/2L=(2.36W/m·K) N/L
Cálculo del coeficiente de convección para el Cálculo del coeficiente
global de transferencia de calor
=7.54kN/2L=7.54X0.138W/m·K X N/2L=(0.520W/m·K) N/L
Cálculo de la cantidad de calor transferida por unidad de tiempo
28
Cálculo de la longitud del intercambiador de calor
Gráfico del comportamiento de la longitud en función del número de Aletas
29
Ejercicio 5.
Se utiliza un intercambiador de calor de tubo concéntrico a contracorriente para enfriar el
aceite lubricante de una turbina industrial a gas. El flujo másico del agua de enfriamiento a
través del tubo interior es de 0,2 kg/s y el diámetro es de 25 mm, mientras que el caudal
de aceite a través del anillo exterior es 0.1 kg/s y el diámetro exterior 45 mm. El aceite y el
agua entran a temperaturas de 100 y 30, respectivamente. Que longitud deben tener la
tubería para que la temperatura de salida del aceite sea menor de 60 °C.
Datos Iniciales
Las propiedades del aceite se escogen de la Tabla A.5 (Incropera, 2013)
30
Interpolación Lineal para las propiedades del aceite y el agua.
y=
31
Cálculo del flujo de calor transferido
Cálculo de la temperatura de salida del agua
Cálculo de la temperatura media logarítmica
Cálculo del número de Reynolds para el agua
El régimen es turbulento porque el número de Reynolds es mayor de
2300
Cálculo del número de Nusselt para el agua
Cálculo del coeficiente de transferencia de calor por convección para el agua
32
Procedimiento para el cálculo del coeficiente de transferencia de calor por
convección para el aceite
Diámetro hidráulico
Con el valor de la relación Se determina el valor del número de Nusselt en la Tabla
8.2
Coeficiente de transferencia de calor por convección para el aceite
33
Cálculo del coeficiente global de transferencia de calor
Cálculo de la longitud del intercambiador de calor
34
