Microturbina de generación de electricidad y calor

Contenido
1. Introducción............................................................................................................... 4
2.2. Justificación. ...................................................................................................... 4
2.2. Alcance. ............................................................................................................. 4
2. Proceso de diseño de una turbo máquina. .......................................................... 6
2.1. Teoría de las turbomaquinas. ......................................................................... 9
2.2. Turbomaquinas centrifugas. ......................................................................... 14
2.2.1. Flujo a través de un TC centrifugo. ..................................................... 15
2.2.2. Parámetros principales de caracterización de una TMT.................. 16
2.3. Turbomaquinas axiales. ................................................................................ 19
2.3.1. Funcionamiento. ..................................................................................... 19
2.3.2. Diseño de las turbomaquinas motoras. .............................................. 20
2.3.3. Teoría adimensional aplicada a turbomaquinas axiales. ................. 21
2.3.4. Teoría del equilibrio radial..................................................................... 21
2.4. Cálculos previos. ............................................................................................ 22
2.4.1. Condiciones del compresor. ................................................................. 23
2.4.2. Condiciones del recuperador. .............................................................. 24
2.4.3. Calculo de mezcla de cámara de combustión. .................................. 24
2.4.4. Condiciones de la turbina. .................................................................... 27
2.4.5. Condiciones del intercambiador de agua-gas. .................................. 27
3. Compresor. .............................................................................................................. 29
3.1. Cálculos preliminares. ............................................................................... 29
3.2. Triángulos de velocidad. ........................................................................... 32
3.3. Trabajo interno............................................................................................ 34
3.4. Rendimiento volumétrico........................................................................... 38
3.5. Cálculos finales. ......................................................................................... 41
3.6. Calculo de la zona de difusión. ................................................................ 45
4. Turbina. .................................................................................................................... 48
4.1. Salto energético, velocidad periférica y diámetro medio. .................... 48
4.2. Diseño de la corona fija. ............................................................................ 49
4.3. Diseño de la corona móvil......................................................................... 53
4.4. Pérdidas, salto y rendimientos periféricos.............................................. 54

1
David Espada García

4.5. Perdidas por rozamiento de disco, salto y rendimientos internos. ..... 55

4.6. Tabla de resultados. .................................................................................. 56
4.7. Parámetros de diseño de los alabes. ...................................................... 57
4.8. Diseño del alabe en la corona fija. .......................................................... 57
4.9. Diseño del alabe en la corona móvil. ...................................................... 61
5. Cámara de combustión. ........................................................................................ 63
5.1. Introducción. ................................................................................................ 63
5.2. Calculo de la temperatura de combustión. ............................................. 65
5.3. Cálculos dimensionales para la cámara de combustión. ..................... 67
5.4. Tipos de cámaras diseñadas. .................................................................. 75
5.5. Conclusión y resultados. ........................................................................... 78
6. Intercambiadores de calor. ................................................................................... 81
6.1. El uso de los tubos aleteados. ..................................................................... 81
6.2. Intercambiador agua/gas. ............................................................................. 82
6.2.1. Datos de entrada. ................................................................................... 83
6.2.2. Balance térmico. ..................................................................................... 83
6.2.3. Características constructivas del intercambiador.............................. 84
6.2.4. Características geométricas de los tubos aleteados. ....................... 86
6.2.5. Dimensiones y paso en el conducto de aire. ..................................... 87
6.2.6. Calculo de la superficie de intercambio de calor. .............................. 90
6.2.7. Resultados. ............................................................................................. 95
6.2.8. Calculo de caída de presión externa................................................... 96
6.3. Recuperador turbo máquina. ........................................................................ 97
6.3.1. Datos de entrada. ................................................................................... 97
6.3.2. Balance térmico. ..................................................................................... 97
6.3.3. Características constructivas del intercambiador.............................. 98
6.3.4. Calculo del área de la superficie de intercambio. ............................. 98
6.3.5. Cálculos dimensionales....................................................................... 100
6.3.6. Caída de presión externa. ................................................................... 100
7. Estructura de la maquina. ................................................................................... 101
7.1. Zona de generación eléctrica y admisión. ............................................ 101
7.2. Zona de intercambio de calor. ................................................................ 101
7.3. Zonas de componentes periféricos. ...................................................... 102

2
 Microturbina de generación de electricidad y calor

7.4. Zonas eléctricas (Potencia y control). ........ ¡Error! Marcador no definido.

8. Presupuesto y planificación. ............................................................................... 103
9. Consideraciones ambientales. ........................................................................... 107
9.1. Etapa de fabricación ....................................................................................... 107
9.2. Etapa de servicio ............................................................................................. 107
9.3. Fin de vida de la maquina ............................................................................... 107
10. Bibliografía. ........................................................................................................... 108

3
David Espada García

1. Introducción.
2.2. Justificación.

Hoy en día el cálculo y diseño de este tipo de maquinas, está muy automatizado
con procesos de alta tecnología como CFD’s que incorporan los propios
programas de elementos finitos que vienen en cualquier software de diseño
industrial. Pese a eso, siempre se tiene que hacer un diseño previo, ir desde lo
más general a lo más preciso y cambiar cosas.

En este proyecto, se cubren las etapas de pre diseño, diseño unidimensional de

cada componente principal de la turbo máquina y finalmente en aquellos
componentes más complejos o los cuales no existe mucho dato empírico ni
bibliografía de cálculo al respecto, se intenta simular mediante CFD, para
comprobar si el funcionamiento es correcto.

Tal y como se puede ver, se va de un nivel más genérico a un nivel mas de

detalle y no se hace un uso constante de las herramientas CFD, sino que
aplicamos conocimientos teóricos del funcionamiento de estas turbo máquinas,
tanto ecuaciones como graficas y apoyados en ocasiones con este tipo de
herramientas.

2.2. Alcance.

En el presente proyecto, se lleva a cabo el diseño completo de la parte

termodinámica y mecánica, de un grupo de generación de electricidad mediante
una turbo máquina.

El proceso a seguir es un cálculo previo de cada uno de los componentes

mecánicos que forman parte del corazón de la maquina generadora de
electricidad y calor. Por lo tanto, el pre cálculo de la turbina, compresor, cámara
de combustión, recuperador e intercambiador de calor. Estos componentes
mecánicos son los que quedan dentro del diseño además de la estructura y
distribución de esos componentes dentro de la maquina general.

Una vez tenemos los cálculos previos, se calcula componente a componente con
mayor detalle mientras se va empezando con los prediseños. Una vez tenemos
los valores finales, se modifica el diseño hasta encontrar el ideal según las
necesidades. Si es necesario se hace un cálculo de CFD para obtener
resultados más precisos.

4
 Microturbina de generación de electricidad y calor

Cuando ya están calculados todos las partes citadas anteriormente, se

ensamblan cálculos y diseños para encajar el conjunto que estará formado por la
turbo máquina, recuperadores, CC y estructura general de la maquina.

La parte electrónica de control y eléctrica de potencia, así como sistemas de

refrigeración externos y lubricación, quedan fuera del alcance del proyecto.

5
David Espada García

2. Proceso de diseño de una turbo máquina.

Para el cálculo de una turbo maquina, el primer paso es un pre cálculo con todos
los componentes. De ahí aplicando las condiciones de contorno, puedes saber
que potencias son las necesarias de obtener con la turbina y cuanto consumo
tendrá el compresor.

El proyecto termodinámico en sí, constituye el modelo matemático con el que se

aproxima el conjunto de procesos reales del gas que tienen lugar en el motor.

Una vez elegido el tipo de ciclo en que va a operar la maquina, se debe hacer un
estudio del ciclo termodinámico que realizara el fluido de trabajo. Con la ayuda
de tablas y datos de otros componentes de motores similares, se hace una
primera estimación del rendimiento de los principales componentes (Compresor,
Cámara de combustión y Turbina). Con estos datos, se puede obtener en
primera aproximación las temperaturas en los puntos de control más importantes
y aproximar el rendimiento total del ciclo.

A partir de estos datos, empieza la segunda fase, perteneciente al diseño

aerodinámico de la turbina y el compresor.

Esta sección, se divide a su vez en:

1. Diseño preliminar: Teniendo fijado algún parámetro crítico (rpm, diámetro

de la turbina, etc.), se calculan los triángulos de velocidad por el método
unidimensional. Se utilizan tablas y formulas para corregir los ángulos de
entrada y salida del fluido de trabajo, calculados mediante la ecuación de
Euler. Posteriormente, se hace un primer dimensionado de la maquina.
Es también común utilizar gráficos experimentales para deducir en
primera aproximación las dimensiones principales y el rendimiento de la
maquina. Una vez se han obtenido las dimensiones principales de la
maquina y los rendimientos de la turbina y el compresor, ya más acordes
a la realidad debido a que se conoce en primera aproximación la
velocidad del fluido y las dimensiones aproximadas de los conductos, se
deben corregir estos en el ciclo termodinámico. En el caso de las
antiguas turbinas de vapor de acción de alabes cilíndricos y en
compresores centrífugos que no exigen gran rendimiento, este es
frecuentemente el único método empleado para su cálculo.
2. Estudio tridimensional del fluido. En el caso de maquinas axiales, se debe
calcular la variación de las propiedades del fluido y de su velocidad de la

6
 Microturbina de generación de electricidad y calor

base a la punta del alabe. Debido a esta variación de velocidad, si se

desea diseñar un alabe sin perdidas por choque, se deberán variar sus
ángulos materiales. Se suelen calcular tres triángulos de velocidades; en
la base, mitad y punta del alabe, pero se pueden calcular tantos como se
quiera, dependiendo del control que se desee tener sobre la superficie
del alabe. Como primera aproximación al análisis tridimensional de las
propiedades del fluido, es habitual partir de la ecuación de equilibrio
radial. En el caso de maquinas radiales con buen rendimiento, se suelen
utilizar métodos teóricos para determinar la distribución de velocidades a
lo largo de la superficie del alabe, pero son bastante más complicados
que en maquinas axiales.

3. Diseño de la sección: La tercera parte, es el diseño de la sección de la

paleta. Este se suele hacer mediante generadores semiautomáticos de
forma, en los que introduciendo una serie de parámetros característicos
calculados, se genera dicha sección. También se puede hacer utilizando
perfiles de ala de avión en el caso de maquinas axiales.

4. Verificación mediante programas de CFD. Actualmente, es habitual

predecir el flujo y mejorar las formas de la maquina mediante programas
de CFD, tanto en radiales como en axiales. Es la última fase del diseño
aerodinámico.

Toda turbomaquina, hidráulica o térmica, se basa en la ecuación de Euler, es

decir, imparte o absorbe energía al fluido mediante los cambios de dirección y de
velocidad.

En el caso de las turbomaquinas hidráulicas, se pueden despreciar los efectos

de la compresibilidad del fluido, pero no ocurre así en las turbomaquinas
térmicas.

Las turbinas de gas pueden funcionar según varios ciclos, pero el aquí estudiado
es el ciclo de Brayton abierto no regenerativo de simple compresión y simple
expansión.

La turbina como maquina, consta de tres órganos fundamentales donde se

producen las transferencias energéticas. A continuación, se describen por orden
de paso del fluido a través de la maquina:

1. Compresor: Absorbe aire atmosférico y eleva su presión a costa de un

cierto trabajo proporcionado por el eje. Existen diversos tipos de

7
David Espada García

compresores. En su clasificación más general, atendiendo al recorrido del

fluido relativo al eje de la maquina, se clasifican en compresores axiales y
centrífugos.

El compresor puede constar de uno o varios escalonamientos. Los

compresores centrífugos dan una relación de compresión mayor por
escalonamiento que los axiales, pero tienen un rendimiento algo peor. El
número de compresiones en un motor turbina de gas, no suele exceder
de dos en el caso de TC centrífugos, y puede ser incluso de dieciocho o
más si monta compresor axial. Es también frecuente la combinación de
ambos tipos, montando primero un número variable pero relativamente
reducido de escalonamientos de compresor axial, seguido por el
escalonamiento final centrifugo. Los rodetes de los compresores
centrífugos, se clasifican a su vez en abiertos, semiabiertos y cerrados.
Los rodetes abiertos casi no se utilizan debido a su mal rendimiento y
baja relación de compresión. Los que dan mayor rendimiento son los
cerrados, pero la relación de compresión por escalonamiento no puede
ser muy elevada debido a la limitación de velocidad periférica del rodete.
Los semiabiertos, dan un buen rendimiento y alta relación de compresión
debido a que pueden alcanzar elevadas velocidades periféricas. La
relación de compresión por escalonamiento en los TC radiales suele ser
de 1,5:1 a 3:1, pudiendo llegar a valores de hasta 12:1 en modelos
experimentales. En cambio, en los TC axiales los valores normales son
de 1,05:1 a 1,2:1, y como máximo realizable 2:1.

Por el contrario, los compresores axiales son en muchos casos, la opción

favorita a pesar del alto número de escalonamientos que requieren, el
menor rango de variación de revoluciones con buen rendimiento y la
mayor dificultad constructiva. Sin embargo, estos compresores dan un
mayor caudal, un mejor rendimiento, y en el caso de turborreactores,
disminuyen el área frontal del motor disminuyendo así la resistencia de
forma o de arrastre, que a altas velocidades es mucho más importante
que la de fricción incrementada debido a la mayor longitud del motor.
Además, el mal funcionamiento a diferentes vueltas del punto de diseño
se soluciona parcialmente con motores de doble eje.

2. Cámara de combustión: Eleva la temperatura del fluido a presión

constante, y por lo tanto, eleva también su volumen específico. Es por lo
tanto, el órgano de la maquina donde se produce la adición de calor al
ciclo. Toda turbina de gas, debe proporcionar un rendimiento de la

8
 Microturbina de generación de electricidad y calor

combustión elevado, una pérdida de presión mínima, y funcionamiento

estable. También deben satisfacer una fácil inspección.

Existen fundamentalmente cuatro tipos de cámaras de combustión. En

aviación o en propulsión en general, se utilizan las cámaras anulares
(turborreactores), tubulares (turbohélices) y anulo tubulares.

Las cámaras anulares, constan de un único conducto anular en el cual se

produce la combustión de los gases que se dirigen a la turbina. Las
cámaras tubulares, constan de un número variable de cámaras de
combustión situadas en la periferia.

En cambio, en las unidades pesadas de TG industriales, se suelen utilizar

cámaras donde el espacio y peso no son tan críticos, resultando cámaras
mucho más voluminosas.

3. Turbina: Es el órgano motor de la maquina. Se encarga de acelerar el

gas y transformar parte de este en trabajo útil en el eje. Las turbinas de
aviación, entregan toda la potencia útil al compresor, y el salto entalpico
restante se produce en la tobera de expansión, incrementando así la
velocidad de los gases y produciendo empuje. En cambio, las TG
utilizadas para propulsión directa o generación de energía, entregan toda
la potencia útil al eje, utilizada para accionar el compresor y el generador
o eje de transmisión.

Es común el montaje de turbinas en dos ejes; uno para la turbina de alta

presión, que acciona el compresor, y otro eje para la turbina de baja
presión que acciona el turboalternador.

Atendiendo al camino recorrido por el fluido, también se clasifican en

turbinas radiales y turbinas axiales, parecidas en forma a los TC.

Para un mismo salto de presión, las turbinas axiales requieren de un

numero de escalonamientos significativamente menor que los
compresores axiales.

2.1. Teoría de las turbomaquinas.

Ecuación de Euler.

Para el estudio energético del flujo a través de una turbomáquina (figura 10.14)
se aplican las ecuaciones de conservación en forma integral al dominio de
estudio entre una sección de entrada al rodete y otra de salida del rodete:

9
David Espada García

Figura 1. Esquema de una turbomaquina.

Si se considera la hipótesis de fluido ideal (coeficientes de transporte nulos:

viscosidad, µ = 0 y conductividad térmica, k = 0) y que la máquina está
perfectamente equilibrada., el momento de giro del eje tiene sólo una
componente principal en la dirección de rotación (dirección axial), la ecuación
vectorial del momento cinético para flujo estacionario, se puede reducir a una
ecuación escalar en dicha dirección, cuya expresión es:

Figura 2. Momentos sobre la turbomáquina.

10
 Microturbina de generación de electricidad y calor

Figura 3. Flujo elemental de entrada a través de un elemento de área.

Dado que el momento en el eje sólo tiene componente en la dirección del propio
eje (Mz en dirección axial), será esa la única dirección considerada en el
intercambio energético del rodete. Desde luego, existirán esfuerzos y momentos
en las otras dos direcciones (radial y tangencial), pero serán mucho menores
que los correspondiente a la dirección z (figura 10.15). Aunque el efecto de
dichos esfuerzos en las otras direcciones pudiera constituir fuente de problemas
en los apoyos, no deberían tener valores comparables a los del par en la
dirección axial.

En la entrada, la sección de paso, es la corona circular entre el cubo (RC) y la

punta (RP) (figura 10.16). En coordenadas cilíndricas, la integral de variación
temporal del momento cinético en la sección de entrada es:

En flujo estacionario, el caudal másico, debe ser igual en la entrada y en la

salida:

El radio medio de entrada por la velocidad tangencial media de entrada es:

11
David Espada García

En cilíndricas, el producto vectorial del vector de posición radial y la velocidad de

entrada es:

Con lo que la componente axial es:

El producto escalar del vector velocidad por el vector área elemental, es menos
el caudal elemental:

De forma totalmente análoga, se obtienen expresiones para la sección de

salida:

El radio medio de salida por la velocidad tangencial media de salida es:

Con todo el momento provocado por el paso del flujo por el rodete es:

12
 Microturbina de generación de electricidad y calor

Por simplificar la expresión anterior, habitualmente se encuentra sin los símbolos

del promedio, y el caudal másico se expresa por la densidad media y el caudal
volumétrico:

La potencia en el eje será el producto de la velocidad de giro por el momento

sobre el eje:

Quedando, la denominada ecuación de Euler de las turbomáquinas:

En donde:

- u1 es la velocidad tangencial del radio medio del borde de ataque o

de entrada del álabe.
- u2 es la velocidad tangencial del radio medio del borde de estela o
de salida del álabe.
- vt1 es la velocidad tangencial del flujo en el radio medio del borde de
ataque.
- vt2 es la velocidad tangencial del flujo en el radio medio del borde de
estela.

La componente tangencial, que hemos denotado por vt; se suele denotar

también, en componente de la velocidad tangencial de un punto del álabe
(u=ωr); es decir por vu. Esta notación es muy útil, en la representación de los
triángulos de velocidades (tanto en la entrada como en la salida):

Velocidad absoluta = velocidad tangencial + velocidad relativa

𝑉 =𝑈+𝑊

Las componentes tangenciales de la velocidad absoluta, de entrada y de salida:

vu1 y vu2 son las que aparecen en la Ec. de Euler. Las componentes normales,
vn1 y vn2, son las que determinan el caudal volumétrico que atraviesa un

13
David Espada García

determinado elemento de área.

Figura 4. Triángulo de velocidades.

2.2. Turbomaquinas centrifugas.

El turbocompresor centrífugo, absorbe energía del eje y la imparte al fluido. Por

ello, está clasificada como Turbomaquina Generadora (TMG).

El incremento de energía en el rodete, esta dado por el incremento de presión

estática y el de presión dinámica. El primero se origina debido al retraso de la
corriente por difusión, y el segundo debido al incremento de energía cinética.

Se define el grado de reacción R de un TC como:

El trabajo de cualquier TMG, según el método unidimensional, viene dado por la

siguiente expresión:

Donde el incremento de presión dinámica viene expresado por el primer término

y el de presión estática por los dos últimos. A veces, pero, es más cómoda la
siguiente expresión, equivalente a la anterior para la energía comunicada al
fluido en un TC.

Una vez el flujo abandona el rodete, pasa al sistema difusor, que se compone de
una corona de alabes fijos que disminuyen la presión dinámica, transformándola
en presión estática a cambio de una disminución de la velocidad. En algunos TC,
existe una voluta para disminuir aun más la velocidad del fluido, pero no tiene
aplicación en el campo de propulsión.

14
 Microturbina de generación de electricidad y calor

2.2.1. Flujo a través de un TC centrifugo.

La corriente que atraviesa un TC centrifugo, no está perfectamente guiada por

los alabes, es decir, las líneas de corriente no siguen perfectamente el contorno
del alabe. Esto es debido a que el flujo que lo atraviesa, con una velocidad
relativa W, se ve afectado por la aceleración de Coriolis. Entonces se genera una
cara de presión y una de succión, y se forman torbellinos en sentido contrario al
de rotación, que en la punta del alabe se traducirá en un decremento de la
velocidad teórica C2u. Este fenómeno, se conoce como deslizamiento o “slip”.

De esta manera, el Angulo de salida de la corriente no es igual al material β2

trazado por el alabe, y por lo tanto variara el triangulo de salida. Se le llama
factor de deslizamiento o “slip” a la relación entre la velocidad real tangencial del
fluido y la velocidad tangencial ideal, y se puede calcular con la ayuda de
ecuaciones. Las más utilizadas son el factor de deslizamiento definido por
Socola y el definido por Stanitz.

También se utiliza el denominado factor de disminución del trabajo, ez, definido

por:

Hay varias formulas para estimar este coeficiente. Una de las que más
concuerda con la experiencia, es la siguiente, desarrollada por Eckert:

Donde z, es el numero de alabes del TC.

15
David Espada García

Figura 5. Esquema de las velocidades de salida.

Podemos observar como hay una disminución de trabajo debido al efecto “slip.”

2.2.2. Parámetros principales de caracterización de una TMT.

Es frecuente en el diseño de los TC, el uso de tablas y diagramas para estimar

las formas preliminares del compresor rápidamente. Estos diagramas, se
obtienen a partir de deducciones teóricas y experimentales, y utilizan
coeficientes adimensionales para la caracterización de la maquina a proyectar.

Los parámetros principales para la caracterización y dimensionado del

compresor son los siguientes:

- ψ Coeficiente de presión: Es la relación entre el salto adiabático

isentropico comunicado al fluido y la energía cinética del rotor en el radio
máximo (aunque puede variar el radio considerado dependiendo de la
formula desarrollada).

- υ Coeficiente de caudal: Muy utilizado en el diseño de TC. Es la relación

entre el caudal volumétrico en la arista de baja presión del
escalonamiento y un caudal hipotético de valor Área periférica del rodete
por su velocidad tangencial.

16
 Microturbina de generación de electricidad y calor

- Número específico de revoluciones. Se definen varios números

específicos de revoluciones. Los mas empleados son:

a) Nq: Numero específico de revoluciones en función del caudal

volumétrico

b) Ns: Numero específico de revoluciones en función de la potencia

c) No: Numero especifico de revoluciones adimensional

d) σ: Numero específico de revoluciones de los TC

El más utilizado en el diseño de TC, son n0 y σ definidos por

Prefiriéndose en muchos casos σ por depender de otros coeficientes

adimensionales. Utilizando las ecuaciones (5) y (6), se obtiene que

Y si sin cometer un gran error, suponemos que la raíz de la división entre

densidades es prácticamente la unidad, se tiene que:

Los gráficos que dependen de coeficientes adimensionales, son muy utilizados

en el proyecto preliminar de los turbocompresores. En el cálculo de este
proyecto, se utiliza el diagrama de Eckert, válido para ángulos de entrada de
aproximadamente 30o.

17
David Espada García

Figura 6. Diagrama de Eckert.

Para estimar el numero optimo de alabes, también es frecuente la utilización de

gráficos y formulas. El siguiente grafico es de gran ayuda para estimar este
parámetro.

Figura 7. Diagrama de correlación dimensional.

18
 Microturbina de generación de electricidad y calor

Si se conoce la relación de diámetros, se obtiene la curva característica del

grafico, mediante la fórmula:

2.3. Turbomaquinas axiales.

La turbina, está compuesta por una corona fija de toberas en el estator y una
corona móvil de toberopaletas en el rotor. Está clasificada como turbomaquina
motora (TMM), debido a que entrega potencia al eje.

2.3.1. Funcionamiento.

El fluido de trabajo, que ha salido del compresor a cierta presión, y se le ha

aumentado la temperatura a presión constante en la cámara de combustión,
entra en la primera

Corona de paletas del estator. La función de estas, es acelerar y cambiar la

dirección

Del fluido. Lo primero se hace a costa de una reducción de la sección de paso, y

el cambio de dirección, importante para obtener la componente tangencial de la
velocidad que es la componente que entrega trabajo útil a la turbina, se hace
gracias al cambio de Angulo material de la tobera.

Al abandonar la corona de toberas, el fluido entra al rotor con una cierta

velocidad relativa igual a la velocidad de salida C1 menos la velocidad periférica
de la toberopaleta U1.

Finalmente, la toberopaleta, vuelve a acelerar el fluido y a aplicarle un cambio de

dirección, produciendo así una fuerza sobre la toberopaleta, que multiplicado por
su velocidad periférica, da el trabajo útil entregado por el fluido.

La componente axial de la velocidad, a pesar de no dar trabajo útil a la turbina,

es el que permite el flujo de masa, y por lo tanto es esencial para el
funcionamiento de la turbina, ya que sin él, no existiría flujo másico y la potencia
seria nula. Además, en los motores de aviación, es la velocidad axial del chorro
la componente responsable del empuje del motor.

Igual que en los TC, definiremos el grado de reacción como

19
David Espada García

En las turbinas de gas, los valores más usuales de R, oscilan entre 0,2 y 0,5.

Las TMM, en general gozan de un mayor rendimiento que las TMG, debido a
que en ellas se acelera el gas a costa de una reducción de presión, y por lo tanto
existe un gradiente de presiones negativo que disminuye los efectos de
separación de la capa limite.

2.3.2. Diseño de las turbomaquinas motoras.

El proyecto preliminar de los alabes, se hace siguiendo los siguientes pasos:

- Un primer estudio mediante el método unidimensional, o ecuación de

Euler, que da los parámetros principales en la sección media del alabe.

- Posteriormente los alabes se diseñan en 3D. Esto suele hacerse

combinando la teoría del equilibrio radial con alguna condición impuesta.
La más utilizada es la condición de diseño con Caxial cte, conocido como
diseño de torbellino libre, aunque también es muy utilizada la condición
de α cte de la base a la punta ya que los alabes son cilíndricos y por lo
tanto económicos. Es habitual calcular las propiedades del fluido en las
base, sección media y punta del alabe, aunque cuando el proyectista
desea un gran control de las formas, se pueden calcular para muchas
más aun cuando la longitud del alabe sea corta.

- Una vez definidos los principales parámetros en las diferentes secciones

del alabe, se introducen en un generador de formas en 2D. Normalmente,
suelen estar dotados de aplicación para predecir la velocidad del fluido.
Estos programas deben gozar de una fácil modificación de las formas,
debido a que es normal que en la cara de succión del alabe de un primer
diseño se formen zonas de flujo supersónico, y se debe poder corregir
fácilmente.

- Finalmente, se crea una estructura tridimensional del alabe con las

secciones salidas del generador 2D de formas, y se hace un análisis
computacional para verificar que se cumplen las predicciones, ya que el
comportamiento del fluido interactuando con todo el alabe, en ocasiones
puede ser muy diferente de un análisis en 2D.

20
 Microturbina de generación de electricidad y calor

2.3.3. Teoría adimensional aplicada a turbomaquinas axiales.

Como se ha dicho, en el proyecto preliminar de toda TMT, se utiliza la ecuación

de Euler. En las TM axiales, al considerar el radio del rodete constante, la fuerza
ejercida sobre la paleta, vendrá dada por:

Pero

Y como u1=u2

Y el trabajo, vendrá dado por

Con esta ecuación, se conocen los triángulos de velocidades de la turbina, y el

estado termodinámico del gas ya que puede ser calculado conociendo el
coeficiente de velocidad de las toberas. Para conocerlo, simplemente se deberá
suponer Cp. constante, y la variación de entalpia vendrá expresada por la
siguiente ecuación:

Posteriormente, se establecen las dimensiones principales mediante la ecuación

de continuidad y conociendo el estado termodinámico del fluido en cada sección.

2.3.4. Teoría del equilibrio radial.

En una fase ya más avanzada del proyecto de la turbina, se considera la

variación de las propiedades del fluido en el sentido radial del alabe. La fuerza
centrifuga que actúa en el conducto debido a la componente tangencial de la
velocidad, debe anularse con un gradiente positivo de presiones del fluido en

21
David Espada García

dirección radial. Para llegar a una solución analítica del problema, se hacen las
siguientes suposiciones:

- Flujo estacionario
- Flujo no viscoso y adiabático
- Las propiedades del fluido en un corte radial son iguales
- No existe velocidad radial

Si finalmente, se supone que la entalpia total no varía respecto al radio

Se deberá cumplir la ecuación de equilibrio radial:

Figura 8. Teoría del equilibrio radial.

Imagen 3. Elemento de fluido. Teoría del equilibrio radial

2.4. Cálculos previos.

Antes de empezar a calcular los componentes con más detalle, se hace unos
cálculos previos sencillos y que tienen en cuenta todos los componentes, para
situar la maquina dentro de unas condiciones de contorno.

Condiciones de entrada:

- Temperatura ambiente: 15ºC

- Gasto másico de una maquina de 100 kW: 0.78 kg/s
- Presión atmosférica: 1 bar.

22
 Microturbina de generación de electricidad y calor

- Presión salida del compresor: 4.5 bar.

- Calor especifico de gases genérico: 1.01 kJ/kg
- Constate de gases: 0.287 kJ/kg
- Gamma: 1.4

Se diseña una maquina de este tipo:

Figura 9. Esquema general de la máquina.

Con las condiciones de contorno expuestas anteriormente, pasamos al cálculo

general de cada componente para situar los parámetros necesarios para
empezar a posterior, a calcular detenidamente todos los componentes y así,
poder diseñarlos.

2.4.1. Condiciones del compresor.

Suponemos una temperatura de entrada de 15ºC, según lo establecido.

El compresor trabajara a una presión de 4,5 bar y con un caudal de aire

aproximado de 0.78 kg/s. El calor específico será de 1 kJ/kg,

Suponemos un rendimiento isotrópico de 0.85

4.5
0.287·ln
1
1
𝑇2𝑠 = 𝑒 · 288.15 = 441.8 𝐾

Por lo tanto:

441.8 − 288.15
𝑇2 = + 288.15 = 468.9 𝐾 (195.8º𝐶)
0.85

Aplicando la ecuación de balance de energía:

23
David Espada García

𝑊𝑐 = 1 · 0.78 · 468.9 − 288.15 = 142.4 𝑘𝑊

2.4.2. Condiciones del recuperador.

Fijamos una temperatura de salida del recuperador de unos 270/290ºC, en los
gases de escape. Esta temperatura debe ser así de elevada, ya que luego se
necesita suficiente energía calorífica para generar calor y calentar el agua.

Por lo tanto tenemos:

Temperatura de entrada de aire a 195.8ºC y una temperatura de gases que

provienen de la turbina a unos 590ºC. Fijamos la temperatura de salida de los
gases a unos 290ºC.

Asignamos un rendimiento del 80% al recuperador, para ser conservadores.

Hacemos un balance energético, y obtenemos:

Calor cedida de los gases de escape:

𝑄𝑐 = 1.01 · 078 · 590 − 290 = 189 𝑘𝑊

Como tiene un rendimiento del 80%, la temperatura que vera justo antes de la
cámara de combustión, será:

189 · 0.8
𝑇3 = + 195.8 = 435.7 º𝐶
0.78 · 1.01

2.4.3. Calculo de mezcla de cámara de combustión.

Este apartado es muy importante, para calcular exactamente la mezcla del

sistema y por lo tanto la cantidad de fuel y aire necesario para la combustión a
las condiciones exigidas.

La composición de los gases de escape, es la siguiente:

Siguiendo la formulación de la combustión para exceso de aire, tenemos que:

79 79
𝐶𝐻4 + 2𝜆 𝑂2 + 𝑁 𝐶𝑂2 + 2𝐻2 𝑂 + 2𝜆 · 𝑁 + 2 𝜆 − 1 𝑂2
21 2 21 2

Efectuando el balance en la cámara de combustión, sabemos que:

𝑚𝑘 · 𝑕𝑘 − 𝑚𝑘 · 𝑕𝑘 = 𝑊 + 𝑄𝑙𝑜𝑠𝑡
𝑃𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑅𝑒𝑎𝑐𝑡𝑖𝑣𝑜𝑠

24
 Microturbina de generación de electricidad y calor

Y sabemos que para cada producto y reactivo, tenemos que:

Reactivos:
𝑇𝑓𝑖𝑛
𝑓
𝑕𝑓𝑢𝑒𝑙 = 𝑕𝑓𝑢𝑒𝑙 + 𝐶𝑝𝑓𝑢𝑒𝑙 𝑇 𝑑𝑇
298𝐾

𝑇𝑎𝑖𝑛
𝑓
𝑕𝑂2 = 𝑕𝑂2 + 𝐶𝑝𝑂2 𝑇 𝑑𝑇
298𝐾

𝑇𝑎𝑖𝑛
𝑓
𝑕𝑁2 = 𝑕𝑁2 + 𝐶𝑝𝑁2 𝑇 𝑑𝑇
298𝐾

Productos:
𝑇𝑐𝑐
𝑓
𝑕𝐶𝑂2 = 𝑕𝐶𝑂2 + 𝐶𝑝𝐶𝑂2 𝑇 𝑑𝑇
298𝐾

𝑇𝑐𝑐
𝑓
𝑕𝐻2 𝑂 = 𝑕𝐻2 𝑂 + 𝐶𝑝𝐻2 𝑂 𝑇 𝑑𝑇
298𝐾

𝑇𝑐𝑐
𝑓
𝑕𝑂2 𝑝 = 𝑕𝑂2 + 𝐶𝑝𝑂2 𝑇 𝑑𝑇
298𝐾

𝑇𝑐𝑐
𝑓
𝑕𝑁2 𝑝 = 𝑕𝑁2 + 𝐶𝑝𝑁2 𝑇 𝑑𝑇
298𝐾

Si cada calor específico, no depende de la integral, porque son valores fijos, y

aplicamos la ecuación de la combustión para productos y reactivos, obtenemos
que:

288.15 𝑇4
Δ𝐻 = − Δ𝐻r + Δ𝐻r288 + Δ𝐻r = 0 (𝐴𝑑𝑖𝑎𝑏𝑎𝑡𝑖𝑐𝑜)
𝑇3 288.15

Suponiendo una lambda λ=5.369, tenemos que:

ΔHr -144463,401 cal/mol

ΔHr298 -191760 cal/mol
ΔHp 362,59484
cal/mol
-108107,652
T4 1225,42023 k

Con estas condiciones, podemos obtener la temperatura de escape deseada

que son los 1225 K, es decir aproximadamente 940 ºC.

Aplicando la lambda obtenida, en la admisión nos encontramos con lo siguiente:

25
David Espada García

Mezcla admisión
Masa molar
aire 28,84 g/mol
metano 16,04 g/mol
Mezcla adm. 28,5944757 g/mol

Siendo, la mezcla a la entrada de:

(%) Mezcla según lambda

0,01918159 1 Mol metano
0,98081841 51,1333333 Mol aire
1

Para la lambda encontrada, tenemos la siguiente distribución de productos:

Mezcla escape
masa molar mezcla (%)
CO2 44 44 0,02951663
H2O 18 36 0,02414997
N2 28 1131,06933 0,75875794
O2 32 279,616 0,18757547
TOTAL 1490,68533 1
Mezcla esc. 28,9810682 g/mol

Podemos ver en la tabla, los productos más representativos de la combustión

con metano, que es la que aproximamos por su similitud al gas natural.

Llegados a este punto y con los porcentajes obtenidos, podemos calcular cuánto
caudal de fuel es necesario:

Caudal total y dividido

ṁ 27,2779962 mol/s 0,78 Kg/s
ṁ fuel 0,52323522 mol/s --> 0,00839269 Kg/s
ṁ air 26,754761 mol/s --> 0,77160731 Kg/s

Haciendo el siguiente calculo, sabemos el porcentaje global dentro del sistema,

que existe entre el fuel y el aire:

26
 Microturbina de generación de electricidad y calor

0.7716 · 1
𝐶𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑎𝑖𝑟𝑒 𝑡𝑜𝑡𝑎𝑙 = = 91.96
0.00839

Es decir, tenemos un exceso de aire en el sistema de 92:1. Es un valor

perfectamente válido dentro del mundo de las turbinas de gas de estas
dimensiones.

2.4.4. Condiciones de la turbina.

Para la turbina, según otras marcas que fabrican este tipo de maquinas, se
supondrán los siguientes valores:

Temperatura de entrada de los gases de la combustión de 940ºC. Este valor lo

obtenemos del cálculo del proceso de combustión de la cámara.

El rendimiento isentrópico, lo marcamos al igual que en el compresor, al 85%.

Por lo tanto:

1
0.287·ln
4.5
1.01
𝑇5𝑠 = 𝑒 · 1225.4 = 799 𝐾

Aplicando el rendimiento insentropico, tenemos:

799 − 1225.4
𝑇5 = 0.85 + 1225.4 = 863.2 𝐾 (590º𝐶)
288.15

Así por lo tanto, el trabajo de la turbina, será:

𝑊𝑡 = 1.01 · 0.78 · 1225.4 − 799 = 285 𝑘𝑊

2.4.5. Condiciones del intercambiador de agua-gas.

Fijando una temperatura de entrada y salida de los gases de escape, podemos

obtener cual es la potencia calorífica que generamos:

Fijamos unas temperaturas de entrada y salida de 290 y 70ºC respectivamente,

por lo tanto la potencia de calor, es:

𝑄𝑐 = 0.78 · 1.01 · 290 − 70 = 168.4 𝑘𝑊

27
David Espada García

Es el calor que generamos, ahora bien, si usamos el agua de entrada a una

temperatura de 50ºC, precalentada y un caudal de agua de 1.6 kg/s, para
mantener una velocidad alrededor de los 2 m/s en una tubería convencional,
obtenemos una temperatura de salida de:

168.4 · 0.8
𝑇2 𝑎𝑔𝑢𝑎 = + 50 = 70º𝐶
1.6 · 4.18

Por lo tanto podemos calentar como mínimo el agua a 70ºC, para un caudal de
1.6 kg/s. Si bajamos el caudal y por lo tanto la velocidad del agua, podemos
calentar a mayor temperatura y siempre estaremos por debajo de la velocidad de
confort del flujo dentro de una cañería.

28
 Microturbina de generación de electricidad y calor

3. Compresor.
A continuación se lleva a cabo el cálculo del compresor en todos sus aspectos,
es decir, calculo de presiones absolutas, velocidades de salida y entrada así
como formas geométricas de los alabes y la voluta.

3.1. Cálculos preliminares.

Datos de entrada
Qe 0,63246501 m3/s 2276,87402 m3/h
Pe 1 bar
Te 15 ºC
he 15,07 kJ/kg Tablas
Єoe 1,205 Tablas
n 65000 rpm
RC 4,5

La tabla anterior, nos muestra los datos de entrada del problema, los cuales son
la relación de compresión, la velocidad de giro, la presión, la temperatura de
entrada y caudal.

Los valores remarcados en azul son obtenidos mediante tablas a la temperatura

de 15º del ambiente.

Para el cálculo de la relación de compresión total, usamos la siguiente fórmula:

𝑃2 4.5
𝜖𝑐 = = = 4.5
𝑃1 1

Multiplicando la relación de compresión total por la presión inicial obtenemos la

final:

𝜖03 = 𝜖𝑐 · 𝜖0𝐸 = 4.5 · 1.205 = 5.423

Con este valor, entramos en las tablas y buscamos h3s total:

𝑕3𝑠 𝑡𝑜𝑡 = 169.6 𝑘𝐽/𝑘𝑔

Por lo tanto, la energía aportada del proceso de compresión es:

𝑘𝐽
𝑌𝑠 = 𝑕3𝑠 𝑡𝑜𝑡 − 𝑕𝑒 = 169.6 − 15.07 = 154.3
𝑘𝑔

29
David Espada García

Pasamos a calcular el número específico de revoluciones. Es un coeficiente muy

usado como ya se ha explicado anteriormente, para la comparación y tablas de
los compresores.

65000 −0.75
𝜍 = 2.108 · · 0.632 · 154.3 · 1000 = 0.233
60

Con el coeficiente de número especifico de revoluciones, entramos en el

diagrama de Eckert y obtenemos lo siguiente

Figura 10. Diagrama de Eckert en nuestro diseño.

Tal y como salen marcadas las líneas rojas, obtenemos los valores:

β1 30 º
β2 50 º
ψ 1,09 (coef. Pres)
ϕ 0,0625 (coef.Caud)
d1/d2 0,513 (relacion de diámetros)

Por el mero hecho de usar el grafico anterior, usamos un Angulo β1=30º, a la

entrada de los alabes.

30
 Microturbina de generación de electricidad y calor

Aplicando la ecuación de la velocidad periférica, obtenemos:

2 · 𝑌𝑠 · 1000 2 · 154.3 · 1000 𝑚

𝑈2 = = = 532.5
𝜓 0.0625 𝑠

Y por lo tanto, el diámetro exterior del rodete será de:

𝑈2 532.5
𝑑2 = = = 0.1565𝑚 = 156.5 𝑚𝑚
𝜋 · 𝑛 𝜋 · 65000

Teniendo la relación de diámetros obtenida con el diagrama de Eckert, podemos

encontrar el diámetro interno:

𝑑1 = 0.513 · 0.1565 = 0.0803 𝑚 = 80.3 𝑚𝑚

Ahora de nuevo usamos otro diagrama, pero en primer lugar calculamos el

Angulo medio y la inversa de proporcionalidad de diámetros:

30 + 50
𝛽𝑚 = = 40º
2

𝑑2 1
= = 1.949
𝑑1 0.513

Por lo tanto, entrando en el siguiente grafico, obtenemos el numero de alabes.

Figura 11. Número de álabes para nuestro diseño.

31
David Espada García

De esta grafica, podemos leer que para la división t/l, según valores de 0.45 o
0.35, obtenemos un resultado u otro. Es decir algo más de 13 alabes y algo más
de 17 alavés respectivamente.

Así que para desempatar, usaremos la siguiente ecuación basada en la

experiencia en donde:

𝛽2 50
𝑧= = = 16.66
3 3

Por lo tanto seleccionamos Z=17 alabes para nuestro rotor.

3.2. Triángulos de velocidad.

El triangulo de entrada será rectángulo, tal que así:

𝑑1 𝑚
𝑈1 = 𝑈2 = 273.2
𝑑2 𝑠
𝑚
𝐶1 = 𝐶1𝑚 = 𝑈1 · tan 𝛽1 = 273.2 · 0.5774 = 157.7
𝑠
𝑚
𝑊1 = 𝑈12 + 𝐶12 = 315.4
𝑠

Para obtener el triangulo de salida, usaremos la siguiente simplificación:

𝑚
𝐶2𝑚 = 𝐶1𝑚 = 157.7
𝑠

Calculo de C2u, para número infinito de alabes.

𝑚
𝐶2𝑢∞ = 𝑈2 − 𝐶2𝑚 · 𝑐𝑡𝑔 𝛽2 = 532.5 − 157.7 · 𝑐𝑡𝑔 50 = 400.15
𝑠

Con numero finito de alabes, mediante el siguiente diagrama podemos encontrar

C2

Calculamos el siguiente coeficiente:

𝑑1
𝑧 1− = 17 · 1 − 0.513 = 8.279
𝑑2

Accedemos a la tabla por el coeficiente calculado y obtenemos ez.

32
 Microturbina de generación de electricidad y calor

Figura 12. Gráfica del valor ez.

Como se explica anteriormente, existe una ecuación que calcula “ez” en función
de diferentes parámetros como el numero de alabes, grados, etc.

Con “ez” calculado, procedemos a hacer el cálculo de la velocidad real:

𝑚
𝐶2𝑢 = 𝑒𝑧 · 𝐶2𝑢∞ = 349.4
𝑠
𝑚
𝐶2 = 𝐶2𝑚2 + 𝐶2𝑢2 = 383.3
𝑠

Por lo tanto, representamos los triángulos de velocidades:

Triangulo de velocidades de entrada:

33
David Espada García

Triangulo de velocidades de salida:

3.3. Trabajo interno.

Despreciando las pérdidas por rozamiento de disco, el trabajo total será:

𝑘𝐽
𝑌𝑖 = 𝑌𝑢 = 𝑈2 · 𝐶2𝑢 = 186.04
𝑘𝑔

Las perdidas internas, o hidráulicas, serán:

𝑘𝐽
𝑦𝑖 = 𝑌𝑖 − 𝑌𝑠 = 186.04 − 157.5 = 31.5
𝑘𝑔

Con estos datos podemos ya conocer el estado termodinámico del aire en las
siguientes secciones de control:

- Entrada del compresor, subíndice E.

- Dentro del compresor, en la boca de aspiración A antes del rodete,
subíndice A.
- Entrada de los alabes, subíndice 1.
- Salida de los alabes, subíndice 2.

A las magnitudes totales se les añadirá el subíndice “tot”

Por lo tanto:

𝑃𝐸 = 𝑃𝐸 𝑡𝑜𝑡 = 1 𝑏𝑎𝑟 ; 𝑇𝐸 = 𝑇𝐸 𝑡𝑜𝑡 = 288.15 𝐾

286.9 · 288.15 𝑚3
𝑣𝐸 = 𝑣𝐸 𝑡𝑜𝑡 = = 0.8267
1 · 105 𝑘𝑔

Si se desprecian las pérdidas entre E y A (entrada del compresor y boca de

aspiración), tendremos que:

𝑃𝐴 𝑡𝑜 𝑡 = 𝑃𝐸 𝑡𝑜𝑡 = 1 𝑏𝑎𝑟 ; 𝑇𝐴 𝑡𝑜𝑡 = 𝑇𝐸 𝑡𝑜𝑡 = 288.15 𝐾

34
 Microturbina de generación de electricidad y calor

𝑚3
𝑣𝐴 𝑡𝑜𝑡 = 𝑣𝐸 𝑡𝑜𝑡 = 0.8267
𝑘𝑔

Las magnitudes estáticas en la sección A se obtienen suponiendo una velocidad

conveniente en dicha sección.

Para evitar desprendimientos de la corriente a la entrada de los alabes se

procura que la corriente se acelere suavemente desde la entrada del rodete. La
práctica aconseja hacer Ca=0.9·C1, es decir:

𝑚
𝐶𝑎 = 0.9 · 157.7 = 141.9
𝑠

Teniendo en cuenta que ɣ=1.4 y que Ra=286.9 N·m/kg·K, obtenemos:

𝛾 1.4 𝐽
𝐶𝑝 = 𝑅𝑎 = 286.9 = 1004
𝛾−1 0.4 𝑘𝑔 · 𝐾

Además se toma el proceso como adiabático e irreversible, lo cual:

𝑃𝑎 · 𝑣𝑎 𝛾 = 𝑃𝑎 𝑡𝑜𝑡 · 𝑣𝑎 𝑡𝑜𝑡 𝛾

Podemos deducir que:

𝑣𝑎 𝑡𝑜𝑡 0.826 𝑚3
𝑣𝑎 = 1 = 1 = 0.9036
2 𝛾−1 2 1.4−1
𝑘𝑔
𝐶𝑎 𝛾−1 141.9 1.4 − 1
1− 1−
𝐴𝑎 2 20 · 288.15 2

𝛾 1.4
𝑣𝑎 𝑡𝑜𝑡 0.826
𝑃𝑎 = 𝑃𝑎𝑡𝑜𝑡 =1· = 0.883 𝑏𝑎𝑟
𝑣𝑎 0.9036

𝑃𝑎 · 𝑣𝑎 0.883 · 105 · 0.9036

𝑇𝑎 = = = 278.08 𝐾
𝑅𝑎 286.9

Para calcular el estado del aire a la entrada y salida de los alabes, es

conveniente estimar las perdidas internas en cada zona del compresor, tal y
como indicamos a continuación.

El sistema difusor constara de una sola caja espiral (sin corona directriz previa)
que se diseñara para reducir un 25% la velocidad de salida del rodete.

- Perdidas hasta la entrada de los alabes:

𝐶12
𝑌(𝑎 − 1) = 𝜁1
2

- Perdidas en el rodete:

35
David Espada García

𝑊12
𝑌(1 − 2) = 𝜁2
2

- Perdidas en el difusor:

𝐶22 − 𝐶32
𝑌(2 − 3) = 𝜁3
2

Para los coeficientes de perdidas, la experiencia recomienda los siguientes

valores:

𝜁1 = 0.1 − 0.15

𝜁2 = 0.2 − 0.25

𝜁3 = 0.25

En nuestro caso, adoptaremos unos valores de:

𝜁1 = 0.135

𝜁2 = 0.25

𝜁3 = 0.25

Por lo tanto calculando las ecuaciones anteriores, obtenemos el siguiente

resultado:

157.72 𝐽
𝑌 𝑎 − 1 = 0.135 = 1678.9
2 𝑘𝑔

315.42 𝐽
𝑌(1 − 2) = 0.25 = 12436.5
2 𝑘𝑔

2
1
383.32 − 4 · 383.3 𝐽
𝑌(2 − 3) = 0.25 = 17218.7
2 𝑘𝑔

Por lo tanto, se tiene que cumplir la ecuación siguiente, es decir, que todo el
trabajo útil, sea equivalente al consumo energético de cada una de las partes:

𝑘𝐽
𝑌𝑖 = 𝑌𝑆 + 𝑌𝑎−1 + 𝑌1−2 + 𝑌2−3 = 185.86
𝑘𝑔

Se cumple el valor, ya que Yi calculado anteriormente nos dio un resultado de

186.03 kJ/kg. Asi, podemos dar por buenos los valores obtenidos.

36
 Microturbina de generación de electricidad y calor

A partir de aquí, se calcularan las presiones, temperaturas y velocidades,

teniendo en cuenta las pérdidas, para obtener valores más precisos por cada
etapa.

En el punto 1, obtenemos:

𝛾 1.4
𝐶𝑎2 + 2 · 𝑌𝑎−1 𝛾−1 141.92 + 2 · 1678.9 0.4
𝑃1 𝑡𝑜𝑡 = 𝑃𝑎 1+ = 0.883 · 1 + = 1.0199
2 · 𝐶𝑝 · 𝑇𝑎 2 · 1004 · 278.07

𝑇1 𝑡𝑜𝑡 = 𝑇𝑎 𝑡𝑜𝑡 = 288.15 𝐾

𝑅𝑎 · 𝑇1 𝑡𝑜𝑡 286.9 · 288.15 𝑚3

𝑣1 𝑡𝑜𝑡 = = = 0.8226
𝑃1 𝑡𝑜𝑡 1.0199 · 105 𝑘𝑔

Y el estado estático en el punto 1, será:

𝛾 1.4
𝐶12 𝛾−1 157.72 0.4
𝑃1 = 𝑃1 𝑡𝑜𝑡 1− = 1.0199 1 − = 0.8746 𝑏𝑎𝑟
2 · 𝐶𝑝 · 𝑇1 𝑡𝑜𝑡 2 · 1004 · 288.15

𝛾−1 0.4
𝑃1 𝛾 0.8746 1.4
𝑇1 = 𝑇1 𝑡𝑜𝑡 = 288.15 = 284.8 𝐾
𝑃1 𝑡𝑜𝑡 1.0199

𝑅𝑎 · 𝑇1 286.9 · 284.8 𝑚3
𝑣1 = = = 0.9046
𝑃1 0.8746 · 105 𝑘𝑔

Para hallar el estado del aire a la salida del rodete se deberá tener en cuenta el
trabajo interno de este. Por tanto:

𝑕2 𝑡𝑜𝑡 = 𝑕𝐸 𝑡𝑜𝑡 + 𝑌𝑖

𝑌𝑖
𝑇2 𝑡𝑜𝑡 = 𝑇𝐸 𝑡𝑜𝑡 +
𝐶𝑝

Substituyendo valores en la ecuación de la temperatura total, obtenemos:

186.03 · 103
𝑇2 𝑡𝑜𝑡 = 288.15 + = 473.4 𝐾
1004

Para encontrar la presión total en el punto 2, aplicamos la siguiente ecuación:

𝑃𝐹 𝑡𝑜𝑡
𝑃2 𝑡𝑜𝑡 = 𝛾
𝑌 𝛾−1
1 − 2 · 𝐶𝑝2−3
· 𝑇2 𝑡𝑜𝑡

37
David Espada García

 Donde Pf tot, es la presión total a la salida del compresor. Sabiendo que

las perdidas Y(2-3), tienen lugar entre la salida del rodete y la salida del
compresor. Aplicando valores:

4.5
𝑃2 𝑡𝑜𝑡 = 1.4 = 5.12 𝑏𝑎𝑟
4622 0.4
1 − 2 · 1004 · 473.4

𝑅𝑎 · 𝑇2 𝑡𝑜𝑡 286.9 · 473.4 𝑚3

𝑣2 𝑡𝑜𝑡 = = = 0.2653
𝑃2 𝑡𝑜𝑡 5.12 · 105 𝑘𝑔

Y para el estado estático en el punto 2, obtenemos:

𝐶22 383.32
𝑇2 = 𝑇2 𝑡𝑜𝑡 − = 473.4 − = 400.3 𝐾
2 · 𝐶𝑝 2 · 1004

𝛾 1.4
𝑇2 𝛾−1 400.3 0.4
𝑃2 = 𝑃2𝑡𝑜𝑡 = 5.12 · = 2.845 𝑏𝑎𝑟
𝑇2 𝑡𝑜𝑡 473.4

𝑅𝑎 · 𝑇2 286.9 · 400.3 𝑚3
𝑣2 = = = 0.6705
𝑃2 2.845 · 105 𝑘𝑔

3.4. Rendimiento volumétrico.

En este caso el rendimiento volumétrico, siempre dependerá de la forma

laberíntica adoptada. Podemos prescindir de la fuga exterior que suele ser muy
pequeña. (ge)

38
 Microturbina de generación de electricidad y calor

Figura 13. Gráfica de relación de tipo de cierre.

En esta figura, escogemos el coeficiente de caudal ζ=1.29, y adoptaremos para

el cierre, Z=4 elementos.

En nuestro caso:

𝑃2
= 3.85 𝑏𝑎𝑟
𝑃1

De la figura, podemos obtener lo siguiente:

𝑔𝑖 𝑇2 𝑡𝑜𝑡
· = 0.246
𝜁𝐴 𝑃2

Para el cálculo del área transversal A, del intersticio, será:

𝐴 = 𝜋 · 𝑑𝑜 · 𝛿𝑖

39
David Espada García

 Donde:

do – diámetro medio de la corona intersticial.

δi – longitud radial del intersticio.

Con el fin de obtener una buena seguridad en marcha del compresor, haremos
δi = 0.2 mm.

Ahora bien, siendo el caudal en la sección A mencionada al principio:

𝜋
𝑄𝑎 = 𝑑𝑎2 − 𝑑𝑒 2 𝐶𝑎
4

Siendo:

4 · 𝑄𝑎
𝑑𝑎 = + 𝑑𝑒 2
𝜋 · 𝐶𝑎

Por otra parte, siendo el caudal másico:

𝑄𝑒 0.6325
𝐺= = = 0.765 𝑘𝑔/𝑠
𝑣𝑒 0.8267

Juntando todos los componentes, obtendremos:

𝑚3
𝑄𝑎 = 𝐺 · 𝑣𝑎 = 0.765 · 0.9036 = 0.6913
𝑠

Para el cálculo del eje, adoptaremos un esfuerzo de torsión máximo admisible de

σ = 19.6 MN/m2.

La ecuación que expresa el par transmitido por el eje del compresor:

𝑃𝑎 · 60 · 103
𝑀=
2𝜋𝑛

 Donde:

Pa – potencia de accionamiento en kW.

n – revoluciones por minuto.

Para hallar la potencia de accionamiento estimamos un rendimiento total del

compresor de un 80%, para ser conservadores. Por lo tanto:

𝐺 · 𝑌𝑠 0.765 · 154.53
𝑃𝑎 = = = 147.7 𝑘𝑊
𝜂𝑡𝑜𝑡 0.8

40
 Microturbina de generación de electricidad y calor

Por lo tanto, el par será:

147.7 · 60 · 103
𝑀= = 21.71 𝑁𝑚
2 · 𝜋 · 65000

Aplicando la formula conocida de la torsión, buscamos el diámetro mínimo del

eje del compresor:

3 16 · 𝑀 3 16 · 21.71
𝑑𝑒 = = = 14.8 𝑚𝑚
𝜋𝜍 𝜋 · 19.6 · 106

Ahora, substituyendo los valores hallados para Qa y de, en las ecuaciones

anteriores, obtenemos:

4 · 0.6913
𝑑𝑎 = + 14.82 = 0.08 𝑚 = 80.7 𝑚𝑚
𝜋 · 141.9

𝐴 = 𝜋 · 𝑑𝑎 · 𝛿𝑖 = 𝜋 · 0.087 · 0.2 = 50.7 𝑚𝑚2

Y el caudal intersticial es:

0.246 · 1.29 · 0.507 · 10−4 · 2.85 · 105

𝑔𝑖 = = 0.023 𝑘𝑔/𝑠
9.81 · 400.27

Por lo tanto, una vez calculados todos estos parámetros, nos vemos con la
capacidad de calcular el rendimiento volumétrico, que será el siguiente:

𝐺 0.765
𝜂𝑣 = = = 0.9704 = 97.04 %
𝐺 + 𝑔𝑖 0.765 + 0.023

3.5. Cálculos finales.

En la recta final del diseño del compresor, se calcularan todas las dimensiones
gracias a los parámetros calculados anteriormente, los cuales tienen importancia
para el propio dimensionado del compresor.

Para el ancho de la entrada de los alabes, tendremos:

𝐺 · 𝑣1 0.765 · 0.9046
𝜂𝑣 0.9704
𝑏1 = = = 0.0179 𝑚 = 17.93 𝑚𝑚
𝜋 · 𝑑1 · 𝐶1𝑚 𝜋 · 0.0803 · 157.7

41
David Espada García

El ancho del rodete a la salida será:

𝐺 · 𝑣2 0.765 · 0.4036
𝜂𝑣 0.9704
𝑏2 = = = 0.0041 𝑚 = 4.11 𝑚𝑚
𝜋 · 𝑑2 · 𝐶2𝑚 𝜋 · 0.156 · 383.3

El rendimiento interno del compresor, sin tener en cuenta las perdidas

intersticiales, llamado también, en ocasiones, por analogía con las bombas
hidráulicas, rendimiento hidráulico ηh, será:

𝑌𝑠
𝜂𝑕 = = 0.8314 = 83.14 %
𝑌𝑖

Por lo tanto, conviene comprobar ahora que el rendimiento total estimado al

principio, se ajusta a los cálculos realizados. Si suponemos un rendimiento
mecánico del 98% se tendrá:

𝜂𝑡𝑜𝑡 = 𝜂𝑕 · 𝜂𝑣 · 𝜂𝑚 = 0.8314 · 0.9704 · 0.98 = 0.7966 = 79.7 %

Habíamos estimado un rendimiento total del 80%, por lo tanto, después de esta
comprobación, podemos continuar los cálculos por que el error cometido no llega
ni al 0.4 %. Lo damos por correcto.

Los alabes se diseñaran mediante arco de circunferencia y el radio del arco se

calculara mediante la siguiente ecuación:

0.156 2 0.0803 2
2
𝑟2 − 𝑟1 2 −
𝜌= = 2 2 = 0.15
2 · (𝑟2 · cos 𝛽2 − 𝑟1 · cos 𝛽1) 0.156 0.0803
2· · cos 50 − · cos 30
2 2

= 145.2 𝑚𝑚

Ahora, el lugar geométrico de los centros de todos los alabes, se encontrara en

la circunferencia de centro del rodete y radio (ρ’)

𝜌′ = 𝑟1 2 + 𝜌2 − 2 · 𝑟1 · 𝜌 · cos 𝛽1

2
0.0803 0.0803
= + 0.14522 − 2 · · 0.1452 · cos 50 = 0.112 𝑚 = 112.2 𝑚𝑚
2 2

42
 Microturbina de generación de electricidad y calor

Representación grafica de las medidas del rotor:

Figura 14. Esquema de medidas y dimensiones del rotor.

Aquí tenemos unas imágenes, donde podemos ver el diseño final del rotor,
carcasa y voluta completa.

Figura 15. Rotor diseñado del compresor.

Este sería el rotor con las medidas obtenidas en los cálculos anteriores. Se
puede ver la inclinación de los alabes para conseguir un arco de circunferencia
adecuado y unos ángulos de entrada y salida establecidos inicialmente.

Los alabes están peinados hacia delante porque de este modo aumenta la
eficiencia del compresor.

43
David Espada García

Figura 16. Voluta calculada y seccionada.

En esta imagen se ve una vista seccionada de la voluta y carcasa del compresor.

En la izquierda, podemos ver la boca de admisión que llega desde el generador

y el filtro de aire.

En ambas imágenes se puede observar, la entrada a la voluta que es una

pequeña ranura y la variación de tamaño de la misma.

En la imagen derecha, se puede ver la pista del rotor y la zona de anclaje del
compresor al cuerpo de la turbomaquina.

Figura 17. Diseño del encaje con el rotor y generador.

La parte del extremo que coge al rotor del compresor, va atornillada con una
tuerca y fijada mediante presión. La zona estriada, es la encargada de encajar y
mover el alternador.

44
 Microturbina de generación de electricidad y calor

Figura 18. Conjunto de compresor, admisión y generador.

En esta última imagen, se puede ver todo el conjunto del compresor

ensamblado, junto con el generador y la admisión.

El detalle interesante es que la admisión, guía el flujo de aire por el interior de

una cámara donde está situado el generador. Este sistema hace que el aire que
absorbe el compresor, sirva a su vez para refrigerar el generador eléctrico.

El componente de color lila, es la tapa del compresor y cuerpo de la

turbomaquina. Se encarga de conectar este conjunto con el resto de la maquina.

3.6. Calculo de la zona de difusión.

Para las dimensiones de la zona de difusión, tendremos en cuenta valores

basados en la experiencia y ecuaciones explicitas para ello.

Figura 19. Esquema de diseño de la voluta.

45
David Espada García

La experiencia marca que la zona del conducto, tiene que tener un diámetro
máximo (d3) comprendido entre 1.05 y 1.1 respecto a d2.

En este caso tenemos un valor de d2=156.5 mm y nosotros adoptamos un valor

d3=172 mm, lo cual equivale a un 1,09d2 y está dentro de los parámetros
fijados.

Para proceder al cálculo dimensional de la voluta, primero se calcula el

coeficiente de circulación de la velocidad a través del perímetro cerrado que
abarca el rodete (T). Se trata de una magnitud constante, para una voluta dada y
para un régimen determinado de funcionamiento del compresor.

𝑇 = 2𝜋 · 𝑟2 · 𝐶2

Se dividirá la voluta en ángulos de 45º, por lo tanto tendremos 9 dimensiones

diferentes de voluta desde 0º hasta 360º ambos incluidos.

Para ello tenemos que:

𝜃º
𝑞𝜃 = 𝑞
360

Donde q, es el caudal del compresor y qΘ, equivale al caudal que pasa según el
ángulo al cual se encuentre.

Figura 20. Medidas para el diseño de la voluta.

Deduciendo la ecuación que define una voluta circular, en función de las

dimensiones y el caudal, obtenemos la siguiente ecuación para obtener el radio
de la voluta según el ángulo en el cual nos encontramos:

46
 Microturbina de generación de electricidad y calor

𝜃º 𝑞 2 · 𝜃 · 𝑞 · 𝑟3
𝜌= · ±
360 𝑇 360 · 𝑇

Como ejemplo, para un ángulo de 180º, substituyendo la ecuación, obtenemos el

siguiente resultado:

𝑇 = 2𝜋 · 156.5 · 349.4 = 171.7

180º 0.6913 2 · 180 · 0.6913 · 0.172

𝜌= · ±
360 171.7 360 · 171.7

= + → 0.02833𝑚 = 28.33 𝑚𝑚 (∨)

= − → −0.024𝑚 = −24 𝑚𝑚 (𝑋)

Escogemos el valor del “+”, que es físicamente valido y significa que el radio de
la voluta en ese punto es de 28,33 mm.

Extrapolando para todos los ángulos definidos, obtenemos la siguiente tabla:

rho (+) rho(-) (mm) Área

q0 0 0 m3/s 0 0 0 0
q45 45 0,086 m3/s 0,0136 -0,0126 13,660 586,21
q90 90 0,172 m3/s 0,0196 -0,0176 19,61 1208,4
q135 135 0,259 m3/s 0,0242 -0,0212 24,298 1854,7
q180 180 0,345 m3/s 0,0283 -0,0243 28,326 2520,8
q225 225 0,432 m3/s 0,0319 -0,0269 31,935 3204,0
q270 270 0,518 m3/s 0,0352 -0,0292 35,246 3902,9
q315 315 0,604 m3/s 0,0383 -0,0312 38,332 4616,1
q360 360 0,691 m3/s 0,0412 -0,0331 41,239 5342,7

Por lo tanto ya tenemos definido todo el compresor completo, desde el rotor a la

carcasa y sus conductos y acabando por el difusor de salida para reducir la
velocidad e incrementar la presión para conseguir mantener los 4,5 bares
específicos.

47
David Espada García

4. Turbina.
Para el diseño de la turbina, como se ha explicado anteriormente, podíamos
escoger dos tipo de conceptos. Por un lado la turbina centrifuga de similar
funcionamiento con el compresor pero a la inversa o por otro, una turbina axial.
Para poder desarrollar ambas tecnologías, se ha decidido por una turbina axial.

Para empezar, expondremos los datos de entrada que son los equivalentes a la
salida de la cámara de combustión:

Cálculos turbina axial de 1 etapa

Presión total 4,41 bar --> Po tot

Presión salida 1,1 bar
To in turbina 1223,15 k --> To tot
Velocidad giro 70000 rpm
Caudal másico 0,78 kg/s

Se diseñara la turbina con un grado de reacción en la mitad del alabe σm=0.3 y

se diseñaran los alabes con torsión, según la ley del remolino libre.

Supondremos ɣ = cte, para simplificar los cálculos y no tener que usar tablas. El
resultado es más directo y menos lioso sin dejar de lado unos valores bastante
precisos.

Asi por lo tanto, en los cálculos de la turbina se usaran las siguientes

simplificaciones:

- ɣ=cte=1.34
- Coef. Obstrucción a la salida de los alabes fijos y móviles, τ=1.
- La constante de los gases en la turbina será Ri=287.4 Nm/kg·K

4.1. Salto energético, velocidad periférica y diámetro medio.

Salto en la turbina:

𝛾−1
𝛾 𝑃2 𝛾
𝑌𝑠 = 𝑇𝑜𝑡𝑜𝑡 𝑅𝑖 · 1 − =
𝛾−1 𝑃𝑜𝑡𝑜𝑡

Substituyendo valores:

48
 Microturbina de generación de electricidad y calor

1.34−1
1.34 1.1 1.34 𝑘𝐽
𝑌𝑠 = 1223.15 287.4 · 1 − = 411.41
1.34 − 1 4.41 𝑘𝑔

Por otro lado, obtenemos la velocidad según la siguiente ecuación:

𝑚
𝐶𝑠 = 2 · 𝑌𝑠 = 2 · 411.41 · 103 = 907.09
𝑠

La velocidad periférica a la mitad del alabe Um, se elegirá teniendo en cuenta la

siguiente ecuación. Y adoptando luego, un valor previo de Voa, en el valor medio
de 0.4.

Este valor, suele designarse según el tamaño de la turbina, siendo 0.33, para
turbinas de muy baja potencia y 0.5 para turbinas de uso industrial de medio y
gran tamaño.

Por lo tanto, el coeficiente óptimo de velocidad periférica será aproximadamente:

𝑈𝑚
= 0.4 · 1 + 0.8 · 𝜍𝑚 = 0.4 · 1 + 0.8 · 0.3 = 0.496
𝐶𝑠

Por lo tanto:

𝑈𝑚 𝑚
𝑈𝑚 = 𝐶𝑠 · = 907.09 · 0.496 = 449.92
𝐶𝑠 𝑠

Ahora podemos calcular el diámetro en la mitad del alabe. Sera:

𝑈𝑚 449.92
𝑑𝑚 = = = 0.1227 𝑚 = 122.7 𝑚𝑚
𝜋 · 𝑛 𝜋 · 65000

4.2. Diseño de la corona fija.

Este diseño es importante, ya que la corona fija, hace se guía del flujo para los
alabes móviles. Encaran el flujo para obtener el mayor resultado posible.

En este caso se trata de una turbina de reacción, tipo así:

49
David Espada García

Figura 21. Esquema de variación de P-V en la turbina.

En este tipo de turbinas, por cada etapa disminuye la presión linealmente,

mientras que la velocidad oscila.

Como se trata de una turbina de 1 sola etapa, solo habrá una corona fija y otra
móvil. La móvil será como esta:

Figura 22. Tipo de diseño del rotor y alabes.

Después de la pequeña introducción, se procede al cálculo del salto energético

en la corona fija:

kJ
Δ𝑕sf = 1 − σm Ys = 1 − 0.3 · 411.41 = 287.98
kg

Velocidad real de salida del gas después de la corona fija:

𝑚
𝐶1 = 𝑘𝑓 · 2 · Δ𝑕sf = 0.95 · 2 · 287.989 = 790.98
𝑠

50
 Microturbina de generación de electricidad y calor

El calor especifico a presión constante, será:

𝛾 1.34 𝐽
𝐶𝑝 = 𝑅𝑖 = 287.4 = 1132.69
𝛾−1 0.34 𝑘𝑔

La temperatura ideal total (después de una expansión a.i.) a la salida de la

corona fija, será:

Δ𝑕sf 287.98
𝑇1𝑡𝑜𝑡 = 𝑇𝑜𝑡𝑜𝑡 − = 1223.15 − = 968.9 𝐾
𝐶𝑝 1132.69

Con lo cual, podemos calcular la relación de presiones:

𝛾−1 1.34−1
𝑃1 𝑇1𝑡𝑜𝑡 𝛾 968.9 1.34
= = = 0.3992
𝑃𝑜𝑡𝑜𝑡 𝑇𝑜𝑡𝑜𝑡 1223.15

Lo cual, según la tabla siguiente, es menor que la crítica. Por lo tanto los alabes
fijos serán convergentes.

Figura 23. Tabla de valores máximos.

Continuando con los cálculos la presión 1, la obtendremos según el caculo:

𝑃1 = 𝑃𝑜𝑡𝑜𝑡 · 0.3992 = 1.76 𝑏𝑎𝑟

Que sería la presión justa a la salida de la corona fija.

La temperatura real a la salida de la tobera, seria:

51
David Espada García

𝐶12 790.982
𝑇1 = 𝑇𝑜𝑡𝑜𝑡 − = 1223.15 − = 993.69 𝐾
2 · 𝐶𝑝 2 · 1132.69

Y por analogía, el volumen especifico en este mismo punto, viene a ser:

𝑅𝑖 · 𝑇1 287.4 · 993.69 𝑚3
𝑣1 = = = 1.622
𝑃1 1.76 · 105 𝑘𝑔

Para empezar a calcular los triángulos de velocidades y diseñar los alabes,

escogeremos el Angulo frontal de α=17º, que está dentro de lo recomendado
para este tipo de turbinas (15º-25º).

Por lo tanto para un Angulo α=17º, tenemos:

𝑚
𝐶1𝑎 = 𝐶1 · sin 𝛼1 = 790.98 · sin 17 = 210.79
𝑠
𝑚
𝐶1𝑢 = 𝐶1 · cos 𝛼1 = 790.98 · cos 17 = 689.48
𝑠

La longitud radial de los alabes directrices a la salida de la corona fija, será de:

𝐺 · 𝑣1 0.78 · 1.622
𝑙′ 1 = = = 0.01557 𝑚 = 15.57 𝑚𝑚
𝜋 · 𝑑𝑚 · 𝐶1𝑎 𝜋 · 0.1227 · 210.79

Esa misma longitud radial, se adoptara para la entrada, es decir, construiremos

los alabes con longitud radial constante.

En este punto, se comprobara el grado de reacción en la base del alabe, y para

ello aplicaremos la siguiente ecuación, siendo el radio de la base del alabe:

𝑑𝑚 − 𝑙 ′ 1 0.1227 − 0.01557
𝑟𝑏 = = = 0.05359 𝑚 = 53.6 𝑚𝑚
𝑠 2

Y en la zona media:

𝑑𝑚 0.1227
𝑟𝑚 = = = 0.06138 𝑚 = 61.4 𝑚𝑚
2 2

Ahora, se comprueba que el grado de reacción no sea inferior entre 0.05 y 0.1,
para proceder al diseño de los alabes con torsión.

2
1 − 0.3 0.05359
𝜍𝑏 = 1 − 2 · 1 − sin 17 · 1 − = 0.109
0.05359 0.06138
0.06138

Por lo tanto podemos continuar con el diseño en base a alabes con torsión.

52
 Microturbina de generación de electricidad y calor

4.3. Diseño de la corona móvil.

Calculamos la velocidad relativa con los valores obtenidos anteriormente:

𝑊1 = 𝑈𝑚2 + 𝐶12 − 2 · 𝑈𝑚 · 𝐶1 · cos 𝛼1 =

𝑚
449.922 + 720.982 − 2 · 449.92 · 720.98 · cos 17 = 319.1
𝑠

Y el Angulo de salida, equivaldrá a:

𝐶1𝑎 210.79
𝛽1 = tan−1 = tan−1 = 138.6º
𝑈𝑚 − 𝐶1𝑢 449.92 − 720.98

Llegados a este punto, tomamos los coeficientes Km=0.95 y Θw=0.916, que son
parámetros que dependen de W1 y el Angulo de salida, con lo cual:

Figura 24. Gráfica para el cálculo de Kf.

Figura 25. Valor de Kc para buscar Kf.

Con ambas graficas, podemos obtener la km según la siguiente ecuación:

53
David Espada García

𝐾𝑚 = 𝐾𝑎 · 𝐾𝑏 · 𝐾𝑐

Pero para el nivel de proyecto que se está llevando a cabo, basta con estimar
Km=0.95.

𝑚
𝑊2 = 𝐾𝑚 𝜍𝑚 · 𝐶𝑠 2 + Θ𝑤 · 𝑊12 = 0.95 0.3 · 907.092 + 0.916 · 319.12 = 554
𝑠

La temperatura real del gas a la salida, será la siguiente:

𝑊22 − 𝑊12 5542 − 319.12

𝑇2 = 𝑇1 − = 993.69 − = 903.14 𝐾
2 · 𝐶𝑝 2 · 1132.69

Y el volumen especifico en el mismo lugar, será de:

𝑅𝑖 · 𝑇2 287.4 · 903.14 𝑚3
𝑣2 = = = 2.36
𝑃2 1.1 · 105 𝑘𝑔

Para este proyecto, se admite un pequeño solape, entre los alabes de la corona
fija y los de la móvil, por lo tanto usaremos la siguiente ecuación, para encontrar
la longitud radial:

𝑑𝑚 122.8
𝑙1′ = 𝑙 ′ 1 + 2 · = 15.57 + 2 · = 15.8 𝑚𝑚
1000 1000

Como simplificación, al igual que en la corona fija, partiremos de la base de una

longitud radial constante tanto a la salida como a la entrada de la corona así por
lo tanto, l2’=l1’=15,8 mm.

Por lo tanto:

𝐺 · 𝑣2 0.78 · 2.36
𝛽2 = sin−1 = sin−1 = 33.01º
𝜋 · 𝑑𝑚 · 𝑙2 · 𝑊2 𝜋 · 0.1228 · 15.8 · 554
𝑚
𝐶2𝑎 = 𝑊2𝑎 = 𝑊2 · sin 𝛽2 = 554 · sin 33.01 = 301.83
𝑠
𝑚
𝐶2𝑢 = 𝑈𝑚 − 𝑊2𝑢 = 𝑈𝑚 − W2 · cos 𝛽2 = 449.2 − 554 · cos 33.01 = −14.7
𝑠

Uniendo los dos vectores, tenemos una velocidad de salida:

𝑚
𝐶2 = 𝐶2𝑢2 + 𝐶2𝑎2 = 302.19
𝑠

4.4. Pérdidas, salto y rendimientos periféricos.

Perdidas en los alabes fijos:

54
 Microturbina de generación de electricidad y calor

kJ
𝑌1 = 1 − 𝐾𝑓 2 Δ𝑕𝑠𝑓 = 1 − 0.952 287.98 = 28.08
kg

La velocidad relativa teórica a la salida de los alabes móviles, será:

𝑚
𝑊2′ = 𝜍𝑚 · 𝐶𝑠 2 + 𝑊12 = 0.3 · 907.092 + 319.12 = 590.5
𝑠

Perdidas en los alabes móviles:

2
𝑊2′ − 𝑊22 590.52 − 5542 𝑘𝐽
𝑌2 = = = 20.86
2 2 𝑘𝑔

Perdidas por velocidad de salida:

𝐶22 302.92 𝑘𝐽
𝑌3 = = = 45.66
2 2 𝑘𝑔

Para calcular las perdidas intersticiales, necesitamos calcular parámetros

dimensionales de las coronas y aéreas de paso, además de otros factores
relacionados con la tipología constructiva de los alabes.

Este punto lo simplificamos y no lo calculamos por que las perdidas intersticiales

suponen apenas un 0.5% de las pérdidas totales de la turbina.

El porcentaje dado es un valor calculado, ya que en cálculos ha quedado

demostrada la baja importancia de ese factor y por lo tanto, se toma la decisión
de simplificarlo.

Por lo tanto, según la siguiente ecuación, el salto periférico, será:

𝑘𝐽
𝑌𝑢 = 𝑌𝑠 − 𝑌1 + 𝑌2 + 𝑌3 = 411.41 − 28.07 + 20.86 + 45.66 = 316.4
𝑘𝑔

Y el rendimiento periférico será:

𝑌𝑢 316.4
𝜂𝑢 = = = 76.9%
𝑌𝑠 411.41

4.5. Perdidas por rozamiento de disco, salto y rendimientos internos.

Perdidas por rozamiento de disco:

1
𝑁𝑟 = 0.0095 · 𝜌 · 𝑛3 · 𝑑𝑚5 = 0.0095 · · 650003 · 0.12285 = 178.2 𝑊
2.36

Y esto vale a una entalpia perdida de:

55
David Espada García

𝑁𝑟 178.2 𝑘𝐽
𝑌5 = = = 0.228
𝐺 0.78 𝑘𝑔

Por lo tanto, el salto interno es:

𝑘𝐽
𝑌𝑖 = 𝑌𝑠 − 𝑌1 + 𝑌2 + 𝑌3 + 𝑌5 = 316.4 − 0.228 = 316.16
𝑘𝑔

Y su rendimiento interno:

𝑌𝑖 316.16
𝜂𝑖 = = = 76.84%
𝑌𝑠 411.41

Suponiendo un rendimiento mecánico del 98%, la potencia útil de la turbina,

será:

𝑃𝑢 = 𝐺 · 𝑌𝑖 · 𝜂𝑚 = 0.78 · 316.16 · 0.98 = 241.7 𝑘𝑊

El cálculo de torsión de los alabes, será:

𝑙1 15.81
= = 0.13
𝑑𝑚 122.8

Este valor de relación entre el diámetro medio y la longitud radial del alabe, debe
de dar un resultado mayor entre 0.08 y 0.1, por lo tanto estamos en la zona
correcta y significa que los alabes torsionados son necesarios.

4.6. Tabla de resultados.

Calculamos los alabes de la base a la punta según la ley del torbellino libre,
siguiendo el procedimiento correcto.

Magnitud Base Medio Punta

r (mm) 53,59 61,38 69,16
r/rm 0,87 1,00 1,13
u (m/s) 392,86 449,92 506,97
C1u (m/s) 789,61 689,48 611,88
Tg α1 0,267 0,306 0,345
α1 14,95 17,00 19,01
C1a (m/s) 210,79 210,79 210,79
Tg β1 -0,531 -0,880 -2,009
β1 152,02 138,65 116,46
C2u (m/s) -16,81 -14,68 -13,03
C2a (m/s) 301,83 301,83 301,83
Tg β2 0,74 0,65 0,58
β2 36,38 33,01 30,13

56
 Microturbina de generación de electricidad y calor

Nos hemos limitado a hacer los triángulos de velocidades en la entrada y la

salida sobre la base y la punta, que junto con los valores ya calculados sobre la
sección media, tal y como se ha podido observar en la tabla anterior.

4.7. Parámetros de diseño de los alabes.

Para el diseño de los alabes, nos hemos orientado, gracias a los datos teóricos y
empíricos del libro de diseño y eficiencia de turbinas de gas de David Gordon
Wilson.

En el trata las siguientes variables y proceso, para diseñar el perfil del alabe
según los parámetros de entrada:

4.8. Diseño del alabe en la corona fija.

En primer lugar, definimos los parámetros necesarios para poder construir los
alabes. En la siguiente imagen, se pueden ver cuáles son los que cobran mayor
importancia para definirlos:

Figura 26. Medidas para el diseño del álabe.

Se puede observar, que hay varios parámetros y que como veremos a

continuación, bastantes dependen de uno básicamente, que es el ángulo λ.

Para encontrar ese valor, podemos acceder a la grafica siguiente mediante los
ángulos de entrada y salida del alabe.

57
David Espada García

Figura 27. Angulo delta según ángulos de entrada y salida del álabe.

Mediante el ángulo de entrada y el de salida del alabe encontramos el valor de

ángulo λ. Según si trabajamos en la base, mitad o parte superior, tendremos
unos valores u otros, lo cual nos definirá un perfil de alabe diferente con
diferentes resultados.

Una vez tenemos en Angulo lambda, calculamos la división b/c, que es:

𝑏
= cos 𝜆
𝑐

Por otro lado, buscamos el cociente b/s óptimo:

𝑏
= 2 · 𝑐𝑜𝑠(𝛼2)2 · tan 𝛼1 − tan 𝛼2
𝑠 𝑜𝑝𝑡

Llegados a este punto, nos toca suponer un valor de “b” o “c”, en este caso,
hemos supuesto el valor de b, que recomiendan que sea 2/3 de la altura de la
aleta como máximo, así por lo tanto:

𝑏 = 10 𝑚𝑚

Calculamos c:

𝑏
𝑐=
cos 𝜆

Calculamos S, mediante la siguiente ecuación:

𝑏
𝑆=
𝑏
𝑠 𝑜𝑝𝑡

Calculamos el incremento de desviación a la entrada del alabe con la siguiente

fórmula:

58
 Microturbina de generación de electricidad y calor

𝛼1 𝑐
Δ𝜃 𝑖𝑛𝑑 = 14 · 1 − + 9 · 1.8 −
70 𝑠

Para calcular el valor del radio exterior “e” del alabe, debemos usar la siguiente
ecuación, la cual trabaja en unos ciertos valores:

0.25 < 𝑠 𝑒 < 0.625

Nosotros escogemos un valor medio de 0.33 y calculamos e, mediante:

𝑠
𝑒=
0.33

Calculamos la relación entre la distancia “o” de los extremos de alabe y la

distancia radial entre ellos con la siguiente fórmula:

𝑜
= cos 𝛼1
𝑠

Como ya sabíamos el valor de la distancia radial entre alabes “S”, encontramos

el valor de “o”.

Finalmente para acabar, buscamos el espesor del final de la aleta mediante la

siguiente ecuación:

0.05
Δ𝑠𝑡𝑒 = ·𝑠
cos 𝜆

Si substituimos los valores, podemos rellenar la siguiente tabla. En ella hay que
tener en cuenta que el valor de α1, se toma respecto a otro ángulo, por lo tanto,
el nuevo valor del ángulo será:

𝛼1 = 𝛼1′ − 90º

Para el caso de la corona fija, α2 = 0º y la referencia se toma de la misma

manera, por lo tanto no hace falta sumar ni restar ángulos.

Los resultados obtenidos son:

λ 60 58 55
b/c 0,5 0,529 0,573
b/s (opt) 0,498 0,559 0,615
b 10 10 10
c 20,00 18,87 17,43
s 20,06 17,88 16,24

ΔΘind 21,23 20,70 20,54

59
David Espada García

Match salida 1,16 1,16 1,16

Ang. Rectif.2

s/e 0,33 0,33 0,33

e 60,25 53,70 48,76
o/s 0,26 0,29 0,33
o 5,18 5,23 5,29

Δste 2,01 1,69 1,42

A continuación podemos ver imágenes del diseño de los alabes de la corona fija.

Figura 28. Álabe del estator.

Podemos ver cómo queda dentro de la cavidad de la carcasa de la turbina,

ensamblados los álabes. Existe un aro exterior de guía y sujeción para cada
alabe, tal y como se puede apreciar en la imagen. Además hace de función de
“camisa” para el paso del flujo de gas.

Figura 29. Turbina con los álabes del estator montados.

60
 Microturbina de generación de electricidad y calor

4.9. Diseño del alabe en la corona móvil.

Para el diseño de los alabes de la corona móvil, se sigue exactamente el mismo

procedimiento, la única diferencia son los valores de entrada que darán otros
resultados del ángulo λ.

Así por lo tanto, los resultados son los siguientes:

λ 5 20 35
b/c 0,996 0,939 0,819
b/s (opt) 2,279 1,588 1,119
b 10 10 10
c 10,04 10,64 12,21
s 4,39 6,30 8,94

ΔΘind -2,80 5,26 12,61

Match salida 0,51 0,51 0,51

Ang. Rectif.2

s/e 0,33 0,33 0,33

e 13,17 18,91 26,83
o/s 0,59 0,54 0,50
o 2,60 3,43 4,49

Δste 0,22 0,34 0,55

A continuación se exponen las imágenes definitivas de cómo es el diseño de la

turbina en la corona móvil.

Figura 30. Álabes del rotor montados.

61
David Espada García

Se pueden ver montadas en el rotor de la turbina de gas. Se puede observar la

forma variable en torsión radialmente.

Aquí el detalle de los alabes.

Figura 31. Forma según diseño del alabe torsionado del rotor.

Para acabar, se muestra la turbina ensamblada con la zona de difusión que llega
a los intercambiadores de calor.

Figura 32. Montaje de la turbina completa con la zona de escape.

62
 Microturbina de generación de electricidad y calor

5. Cámara de combustión.
Para el diseño de la cámara de combustión, se echo mano de la literatura
especifica de turbo maquinas térmicas y turbinas de gas.

5.1. Introducción.

Las cámaras de combustión son dispositivos que se encargan de introducir una

cierta cantidad de energía calorífica sobre un fluido. Normalmente, en las
cámaras de combustión de una turbina de gas, este proceso suele realizarse a
presión constante y la adición de calor se ve reflejada en el incremento de
temperatura del fluido.

A día de hoy existen tres tipos de cámaras de combustión para turbinas de gas:

Figura 33. Tipos de cámara de combustión.

- Cámaras tubulares o “can”. Este tipo de cámaras se componen por un

revestimiento cilíndrico montado concéntricamente en una carcasa
cilíndrica tal y como muestra la figura 2.1. La principal ventaja de este tipo
de sistemas es su reducido coste de fabricación. Sin embargo, su
longitud y peso excesivos ha hecho que este tipo de cámaras no se
utilice en aeronaves, sino más bien en unidades industriales.

- Cámaras tubo-anulares. El diseño de este tipo de cámaras se basa en un

grupo de forros tubulares, generalmente de 6 a 10, dispuesto en el
interior de una carcasa anular tal como ilustra la figura 2.1. Este concepto
de cámaras trata de combinar la compactación de la cámara anular con
la resistencia mecánica de la cámara tubular. Su principal inconveniente
deriva de la necesidad de tener interconexiones de llama (tubos de fuego
cruzado).

63
David Espada García

- Cámaras anulares. Este tipo de cámaras poseen un revestimiento anular

montado concéntricamente en una carcasa anular tal y como muestra la
figura 2.1. Esta es sin duda, una de las cámaras ideales en cuanto a
limpieza aerodinámica. El resultado de este diseño deriva en una unidad
compacta y unas pérdidas de presión menores que en los otros tipos de
cámaras. Su principal inconveniente viene ocasionado por su excesiva
carga de pandeo en el revestimiento exterior.

Aun así, existen un gran número de variaciones dentro de estos tres grupos.
Tratan de jugar con el diseño de la geometría, los tipos de inyectores de
combustible o incluso la orientación de la propia cámara. Todo ello con tal de
adaptarse a las necesidades del tipo de unidad de motor que se necesite. No
obstante todas las cámaras tienen en común las zonas en las que se divide el
proceso de combustión tal como se ilustra en la figura 2.2. Como observaremos
en grados de modelizaciones dos y tres, la diferenciación de estas zonas resulta
fundamental para el cálculo del valor simulado de emisiones o la temperatura de
transferencia de la pared de cámara.

Figura 34. Zonas características de una CC.

Para un correcto diseño de las cámaras de combustión se debe tener en cuenta

que entran en juego aspectos de transferencia de calor, transferencia de masa,
termodinámica, dinámica de gases y dinámica de fluidos.

Para nuestro caso, se diseñara una cámara de combustión derivada de la tubo

anular, pero únicamente tubular, tal y como se muestra en la figura XXXXX.

Las cámaras de combustión tubulares tienen forma cilíndrica y están montadas

de manera concéntrica en el interior de otro cilindro. Las principales ventajas que
presentan son su simplicidad, su fácil diseño y su fácil acceso.

64
 Microturbina de generación de electricidad y calor

Como problema presentan que son grandes y pesadas en comparación a otros

tipos de cámara de combustión y por ello su aplicación está relegada únicamente
a la industria. A diferencia de otros casos, va montada de manera externa, de ahí
el fácil mantenimiento.

El flujo de comburente puede darse de manera directa o de manera inversa. En

el caso del flujo inverso el aire comprimido entra a la cámara de combustión a
través de un anillo que es concéntrico y exterior al cuerpo del quemador. Se
conoce como flujo inverso debido a que el aire accede en dirección opuesta a la
llama. El aire que circula entre la camisa y el quemador tiene varios puntos de
entrada formando las distintas zonas de llama. Las cámaras de combustión de
flujo inverso presentan la ventaja de ser más cortas que las de flujo directo y se
utilizan en aplicaciones de gran tamaño, aun que más complejas de diseñar
desde el punto de vista fluido térmico.

5.2. Calculo de la temperatura de combustión.

Mediante el algoritmo de un programa en C++, se calculara la temperatura de la

llama aproximada, según el siguiente procedimiento.

El algoritmo se adjuntara en los anexos del proyecto.

Formulación

Igual que en el anterior ejercicio, solo que hay que buscar la Cp media de los
productos:

𝑇𝑐𝑐
1
𝐶𝑝𝑋 = · 𝐶𝑝𝑋 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298𝐾

Algoritmo de cálculo

Es un proceso inverso al anterior.

En primer lugar, mediante los caudales de aire y de fuel, obtenemos la 

necesaria, es decir como es la mezcla.

En segundo lugar, accedemos a la ecuación global donde se relaciona la lambda

con los productos y los reactivos.

Mediante el siguiente proceso iterativo, calcularemos la CP media de cada

producto con una temperatura de combustión propuesta y esa CP calculada la
relacionaremos con las entalpias de cada producto y reactivo, donde allí
aparecerá la Tcc.

65
David Espada García

Esa Tcc, se aislara de la ecuación, la calcularemos y la compararemos con la

propuesta dentro de los márgenes propuestos del error. Si no se cumple, esta
nueva Tcc será la siguiente para iterar y así hasta que converja y este dentro del
error que hayamos fijado.

Proceso de cálculo:

1- Buscamos :

(𝑚𝑎 𝑚𝑓 ) (𝑚𝑎 𝑚𝑓 )
𝑡𝑕 𝑡𝑕
𝜆= = (𝑚𝑎 𝑚𝑓 )
(𝑚𝑎 𝑚𝑓 ) 𝑏 𝑊𝑎 𝑡𝑕
𝑠𝑡 𝑎+ 3.76 + 1 ·
4 𝑊𝑓

2- Proceso iterativo global:

Suponemos Tcc*

Calculamos la Cp media de cada producto (x) con la Tcc*

𝑇𝑐𝑐
1
𝐶𝑝𝑋 = · 𝐶𝑝 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298𝐾 𝑋

Volvemos a calcular Tcc, aislándola de la ecuación de la

energía:
𝑚𝑘 · 𝑕𝑘 − 𝑚𝑘 · 𝑕𝑘 = 𝑊 + 𝑄𝑙𝑜𝑠𝑡
𝑃𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑅𝑒𝑎𝑐𝑡𝑖𝑣𝑜𝑠

| Tcc -Tcc*| < δ NO

SI Cálculos finalizados

66
 Microturbina de generación de electricidad y calor

Cp media de los productos:

𝑇𝑐𝑐
1
𝐶𝑝𝐶𝑂2 = · 𝐶𝑝𝐶𝑂2 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298

𝑇𝑐𝑐
1
𝐶𝑝𝐻2 = · 𝐶𝑝𝐻2 𝑂 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298

𝑇𝑐𝑐
1
𝐶𝑝𝑂2 = · 𝐶𝑝𝑂2 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298

𝑇𝑐𝑐
1
𝐶𝑝𝑁2 = · 𝐶𝑝𝑁2 𝑇 𝑑𝑇
𝑇𝑐𝑐 − 𝑇𝑟𝑒𝑓 298

Los datos de entrada son:

- Caudal de aire: 0.7716 kg/s

- Caudal de fuel: 0.0814 kg/s
- Caudal total: 0.78 kg/s
- Datos de Cp(t) de cada uno de los componentes de los productos y
reactivos que forman parte de la combustión.

Resultado de la temperatura obtenida, es:

- Temperatura de llama = 1734.24 K

-

5.3. Cálculos dimensionales para la cámara de combustión.

No existe ningún tipo de bibliografía que explique el cálculo exacto o

metodología de cálculo para saber el comportamiento de la cámara, pero si
existen muchos datos empíricos derivados de la experiencia de grandes
industrias, las cuales nos basamos para el diseño dimensional de la cámara.

Los principales factores a tener en cuenta en el diseño de la cámara de

combustión de una turbina de gas, dependen de sus condiciones operativas, de
entre las que podemos destacar las siguientes:

67
David Espada García

- La combustión tiene que ser estable, para permitir las fuertes variaciones
de la relación aire-combustible que para los ciclos regenerativos está
entre 60/1 y 120/1, y para los no regenerativos entre 100/1 y 200/1.

- La velocidad del fluido oscila, en la mayor parte de los casos, entre 30 y

60 m/s.

- La masa de aire teóricamente necesaria para la combustión de 1 kg de

combustible se determina a partir de las reacciones estequiometrias de la
combustión; para los combustibles líquidos utilizados normalmente en las
turbinas de gas, la relación aire-combustible estequiométrica está entre
14,7 y 15, pudiéndose tomar 14,9 como valor medio.

- La cantidad de aire real suministrada al combustible es mayor que la

teórica, definiéndose un coeficiente de exceso de aire a, como la relación
entre la cantidad real de aire y la teórica, por kg de combustible.

Para que la combustión tenga lugar totalmente dentro de la cámara de

combustión, es necesario que todos los procesos se realicen con suficiente
rapidez, ya que se tienen que efectuar en una corriente de aire a una velocidad
determinada. Por esta razón, la combustión en estas cámaras sólo es posible
cuando la velocidad de propagación de la llama sea del mismo orden que la
velocidad de la corriente.

La temperatura de la combustión y, en consecuencia, la velocidad de

propagación de la llama, dependen del coeficiente de exceso de aire
obteniéndose su valor máximo para un coeficiente de exceso de aire a = 1,
aproximadamente.

Para mezclas ricas, a < 1, la temperatura de combustión disminuye debido a la

combustión incompleta. En las mezclas pobres, a > 1, también disminuye debido
a la dilución por el aire de los productos de la combustión. Como la relación aire-
combustible de funcionamiento global es del orden de 60/1 para las condiciones
de diseño, mientras que la estequiométrica es de aproximadamente, 15/1, es
necesario que el aire suministrado por el compresor se introduzca
progresivamente en la cámara de combustión.

En el proceso de inyección de aire se pueden distinguir tres fases que dan lugar
a lo que se conoce como aire primario, aire secundario y aire terciario.

68
 Microturbina de generación de electricidad y calor

Aire primario.- Se corresponde aproximadamente con un 15% a 20% del aire

total y se introduce alrededor del chorro de combustible creando una mezcla de
aire-combustible relativamente rica con el objeto de obtener una temperatura
elevada, necesaria para una combustión rápida, r = 1.

Aire secundario.- Se corresponde aproximadamente con un 30% del aire total;

se introduce a través de orificios practicados en el tubo de llama para completar
la combustión; para que el rendimiento sea elevado, se tiene que procurar que
este aire se inyecte en los puntos adecuados del proceso a fin de evitar que la
llama se enfríe localmente dando lugar a una drástica disminución de la
velocidad de combustión en esa zona. El aire secundario proporciona un
coeficiente de exceso de aire del orden de, a »1,5, riqueza 1/1,5, Þ 0,66 a 0,7.

Aire terciario.- El aire restante, 50% a 55%, se mezcla con los productos de la
combustión en la zona de dilución, con el objeto de reducir su temperatura hasta
la requerida a la entrada de la turbina. Hay que procurar una turbulencia
suficiente para que las corrientes caliente y fría se mezclen a fondo y así
conseguir una distribución de temperaturas a la salida prefijada de antemano. En
las turbinas de gas industriales, el tamaño de la cámara de combustión no es un
factor preponderante en el diseño. Las dimensiones principales de una cámara
de combustión se eligen, fundamentalmente, tomando como base datos
experimentales de cámaras similares ya existentes. Una vez construido un
primer prototipo se ensaya y se introducen las modificaciones necesarias para
conseguir un comportamiento óptimo. Este comportamiento viene definido por su
rendimiento, pérdida de presión de remanso, límites de estabilidad de la
combustión, perfil de temperaturas del flujo a la salida de la cámara, etc.

Aunque la forma y disposición de las cámaras de combustión son bastante

flexibles, se puede establecer una normativa basada en datos experimentales,
que sirva de ayuda para el diseño del volumen, sección transversal, longitud, y
tamaño y disposición de los orificios del tubo de llama. El volumen de la cámara
de combustión depende de la intensidad de la combustión, de forma que
definidos los valores de este parámetro, se puede hacer una primera
aproximación del volumen de la misma.

Para seleccionar el área de la sección transversal de la cámara d combustión se

define el concepto de velocidad de referencia cref en la forma:

69
David Espada García

𝐺𝑎𝑖𝑟𝑒
𝐶𝑟𝑒𝑓 =
𝜌2 · Ω𝑚

 Donde:
- G: es el caudal másico de aire.
- ρ2: es la densidad del aire a la entrada a la cámara de combustión.
- Ωm: el area de la sección transversal máxima de la cámara que para el
diseño tubular, vale:
𝜋 · (𝐷𝑒 2 − 𝐷𝑖 2 )
Ω𝑚 =
4
 Donde:
- De: diámetro de la carcasa de la cámara de combustión.
- Di: diámetro interior, de la celda de la combustión.

La velocidad de referencia Cref, es un parámetro teórico y puede que no llegue a

alcanzarse nunca, en ninguna sección de la cámara. Sus valores máximos serán
de unos 35 m/s, siendo para este, las turbinas de menor tamaño y mayor
velocidad de giro. Si se trata de una turbina de grandes dimensiones, industrial,
aun pueden ser más bajos.

Desde el punto de vista del diseño es conveniente expresar la pérdida de

presión de remanso en la cámara de combustión mediante un parámetro
adimensional denominado coeficiente de pérdida de carga CL de la forma:

2 · Δ𝑃2
𝐶𝑙 =
𝜌2 · 𝐶𝑟𝑒𝑓 2

Para la perdida de carga en la cámara de combustión, se supondrá un valor de

un 2,5% de caída de presión respecto a la presión de entrada en ella.

El coeficiente de pérdida de carga, puede ser muy variable, aun que

normalmente se debe de ajustar por debajo de 15, para valores de velocidad
relativa de unos 35 m/s.

En cuanto al área transversal del tubo interno de combustión, suele situarse

entre el 55/65% respecto de la máxima área de combustión de la propia cámara.
La longitud optima, supone que para una velocidad de referencia dada, el
coeficiente de carga, será mínimo.

70
 Microturbina de generación de electricidad y calor

Figura 35. Secciones transversales de una CC.

Por lo que respecta la longitud, la experiencia demuestra que las cámaras

conservan siempre unos mismos valores entre longitud/diámetro. Habitualmente,
estos parámetros, oscilan entre 3 y 6 veces la longitud respecto al tubo de llama,
o entre 2 y 4 veces si se trata del diámetro máximo de la cámara. Los valores
más bajos suelen ser para casos donde existe torbellinador.

Figura 36. Dimensiones generales de una CC.

Los orificios de la cámara de combustión, deben de estar situados tal y como se

comenta anteriormente. En cuanto a la zona primaria es muy difícil de calcular si
no imposible. Normalmente se establece que los agujeros suelen ser 1/10 del
diámetro de la cámara de combustión.

Cálculos:

Datos de entrada:

T3 435,78 ºC
P3 4,5 bar
ρ3 2,202943664 kg/m3
Gaire 0,7716 kg/s

71
David Espada García

Calculamos Ωm, y para ello fijamos primero unos valores de diámetros

adecuados (han sido tanteados con una hoja de cálculo hasta dar con los mas
correctos).

- De = 0.178 m
- Di = 0.138 m

Comprobamos que la relación de áreas se mantenga entre 0.55 y 0.65.

Aplicando los diámetros descritos anteriormente, la relación de áreas es 0.6, por
lo tanto estamos dentro de los márgenes.

𝜋 · (0.1782 − 0.1382 )
Ω𝑚 = = 0.00992743 𝑚2
4

Ahora podemos calcular la velocidad de referencia, que recordemos que no

puede superar los 35 m/s.

0.7716 𝑚
𝐶𝑟𝑒𝑓 = = 34.98
2,2029 · 0.00992743 𝑠

Nos encontramos dentro del valor adecuado.

Si aplicamos una caída de presión de un 2,5%, obtenemos:

∆𝑃3 = 4.5 · 0.025 = 0.1125 𝑏𝑎𝑟

Por lo tanto el coeficiente de pérdida de carga, será:

2 · 0.1125
𝐶𝑙 = = 8.205
2.2029 · 34.982

Es un valor que está por debajo de 15, por lo tanto, es un valor aceptable y
coherente.

Ahora calcularemos las zonas de combustión y dimensiones generales.

Min Max
Relacion con Dint 3 6
Relacion con Dext 2 4
Longitud según Dint 0,414 0,828 m
Longitud según Dext 0,356 0,712 m
* Debe de estar dentro de los valores max/min
Adoptamos una Long. De 0,4 350 mm

72
 Microturbina de generación de electricidad y calor

Como se trata de una cámara de combustión con torbellinador, podemos ajustar

las dimensiones. Por lo tanto nos centraremos en el valor más alto de los
mínimos.

Longitud de la cámara = 420 mm.

Ahora queda dividir la cámara en diferentes secciones de zona primaria y

secundaria + zona de dilución. Para ello usamos la siguiente característica:

Figura 37. Distribución de zonas en la CC.

Como se puede observar, se recomienda que menos del 30% de la longitud,

pertenezca a la zona primaria y secundaria y el otro más de 70% pertenezca a la
zona de dilución.

Por lo tanto, las medidas de diseño serán:

Zona primaria + secund. 117,6 mm 28%

Zona de dilución 302,4 mm 72%
Longitud total se cumple 420 mm

Para acabar con el proceso de diseño dimensional, se hizo una grafica

suponiendo diámetros variables. Como existe una relación de diámetros que va
desde 0,55 a 0,65, se juega con los valores de diámetro para encontrar la
cámara de combustión más equilibrada:

73
David Espada García

Cref 35 35 35 35 35 35 35 35 35 35 35
Ωm 0,01 0,01 0,01 0,01 0,01 0,01 0,01 0,01 0,01 0,01 0,01
RA 0,55 0,56 0,57 0,58 0,59 0,6 0,61 0,62 0,63 0,64 0,65
De 0,168 0,170 0,172 0,174 0,176 0,178 0,181 0,183 0,186 0,188 0,191
Di 0,125 0,127 0,130 0,133 0,135 0,138 0,141 0,144 0,147 0,151 0,154
Δd 21,74 21,41 21,09 20,76 20,44 20,11 19,79 19,47 19,14 18,81 18,49
CL 8,34 8,34 8,34 8,34 8,34 8,34 8,34 8,34 8,34 8,34 8,34
Lmin 0,374 0,382 0,390 0,398 0,406 0,415 0,424 0,433 0,442 0,452 0,461
Lmax 0,673 0,681 0,689 0,697 0,705 0,714 0,723 0,732 0,742 0,753 0,763
Zona 1+2 126
Zon. Diluc. 294

Podemos ver el comportamiento de diferentes parámetros en la grafica:

0,200 22
0,190 21,5
0,180 21
0,170 20,5
Ø int.
0,160 20
Ø ext.
0,150 19,5
ΔØ
0,140 19
0,130 18,5
0,120 18
0,55 0,57 0,59 0,61 0,63 0,65

Figura 38. Gráfica diametros respecto dimensiones.

Escogemos el punto óptimo personal, en el 0.6, ya que se considera el punto en

el que la diferencia entre diámetros es media, y las dimensiones son contenidas.

Si nos fijamos en la izquierda, la diferencia de diámetros es muy grande, lo que

hace una cámara de combustión muy descompensada y la longitud de la cámara
es demasiado corta para desarrollar el flujo.

En cuanto a la zona derecha, ocurre al revés, ambos diámetros están muy

juntos, creando poca zona para el flujo de aire comprimido. Por otro lado, la
proporcionalidad con la longitud, hace que la cámara tenga dimensiones
demasiado amplias.

74
 Microturbina de generación de electricidad y calor

5.4. Tipos de cámaras diseñadas.

Se han diseñado dos cámaras de combustión distintas que se basan en la figura

XXXXX. Así por lo tanto, ambas tienen sus zonas delimitadas y sus
características parecidas.

Se ha querido hacer 2 diseños diferentes, haciendo hincapié en la economía del

diseño frente al funcionamiento del mismo.

Cámara con torbellinador:

Se trata de una cámara de combustión más elaborada, la cual consta con un

elemento en la entrada. Un torbellinador. No es más que un simple orificio con
forma de aspas de ventilador que redirige el flujo de aire hacia adentro del tubo
de llama en forma de espiral.

Figura 39. Corrientes de flujo en la CC.

Como se puede ver, se sitúa alrededor del inyector y además lo protege de

recibir el aire directo, para no bajar demasiado la temperatura de la llama.

75
David Espada García

Figura 40. Diseño interno de la CC.

La cámara tiene unos orificios distribuidos en el perímetro del tubo de llama, para
así generar la mezcla adecuada entre combustible y comburente.

Tal y como se puede ver anteriormente, existen dos zonas, la de flujo primario y
la de dilución. En cuanto a la zona primaria, se divide en dos entradas. Los
orificios para flujo primario y para los orificios de flujo secundario.

Figura 41. Zonas diseñadas dentro de la CC.

En la zona primara, encontramos dos subzonas, que son exactamente la zona

de flujo primario y la de flujo secundario.

76
 Microturbina de generación de electricidad y calor

Como se ve en la siguiente imagen, los agujeros en forma espiral frontales que

están situados alrededor del torbellinador, filtran el flujo para hacerlos pasar por
una cámara en donde sale el flujo por los pequeños agujeros perimetrales
frontales y la ranura que apunta hacia las paredes del tubo de llama. Todo este
sistema, es el sistema de aire primario, junto con el flujo que pasaría por el
torbellinador.

Seguidamente tenemos los agujeros en disposición radial, en grupos de 2 en 2.

Están en una colocación cercana a la llama, y son de menos diámetro que los de
la zona de dilución. Estas perforaciones, ayudan a bajar la temperatura de
combustión y aumentar la mezcla y transferencia de calor al igual que el
torbellinador.

Figura 42. Torbellinador y agujeros de dilución.

Como se ve en la imagen XXXX, al final del tubo de llama, casi en contacto con
la carcasa exterior, existe un cerramiento, tipo laberintico, el cual deja pasar un
pequeño flujo de aire. Esto se hace para que prácticamente todo el flujo, se vaya
por los agujeros anteriores y poca masa de aire se escape por ahí.

La razón por la cual no se cierra totalmente, es para conseguir que el poco flujo
que pase, que lo hace a una velocidad elevada, refrigere esa zona critica la cual
puede encontrarse a altas temperaturas (del orden de 800º).

Aun que compleja, es una cámara eficiente con un funcionamiento muy correcto.

77
David Espada García

Cámara simple de tubo:

Se opto por probar un nuevo diseño sencillo, para hacer una comparación entre
ambas cámaras, pero los resultados no eran satisfactorios, así por lo tanto, se
desecho esta opción. Todo y eso, se hizo el diseño que se explica a
continuación:

Figura 43. Diseño de cámara simple.

Tal y como se ve en la figura anterior, podemos ver que la cámara de

combustión es más corta, lo cual beneficia en el espacio dentro de la maquina.

La situación de los agujeros del flujo secundario y de dilución, es muy similares

al diseño anterior, y en cuanto a los del flujo primario, son recolocados para tener
una posición igual a los citados anteriormente.

Pese a que la entrada del aire comprimido en ambas cámaras es en el mismo

lugar, en el primer caso el flujo es redirigido de manera axial hacia el tubo de
llama, mientras que en este caso entra de manera transversal.

5.5. Conclusión y resultados.

Podemos ver el comportamiento de la temperatura, que coincide con el flujo de

la llama. La temperatura máxima es bastante superior a los 1200ºC, pero se ha
escalado para poder observar resultados.

78
 Microturbina de generación de electricidad y calor

Figura 44. Resultado de temperaturas.

Podemos observar que la temperatura de entrada coincide con la fijada de más

de 400ºC, la cual proviene del compresor y el recuperador.

En la zona de la salida, nos encontramos a una temperatura de 984ºC.

Se puede ver la zona en donde, debido a la entrada de flujo por los agujeros,
hay disminución de temperatura. Un incremento de mezcla
combustible/comburente.

En la zona de los cerramientos, la temperatura esta algo por encima de los 600
en el fluido. Sin ello podría sufrir la carcasa de metal.

Figura 45. Resultado de número de Mach.

El número de Mach, es proporcional a la velocidad.

79
David Espada García

Se puede ver las zonas donde la velocidad es mayor, como la entrada o salida.
Además podemos observar la zona del torbellinador, como se acelera el flujo, al
igual que en las paredes debido a los agujeros y en la zona central.

Figura 46. Resultados de presión.

La presión es prácticamente constante a la entrada y salida de la cámara de

combustión. Por lo tanto reafirmamos la simplificación de combustión a presión
constante.

La zona central, debido a la temperatura del flujo incrementa la presión

ligeramente.

Figura 47. Resultados de velocidad.

La velocidad es proporcional a Mach, por lo tanto vemos unos valores similares

en otra escala.

La velocidad de entrada y salida, son muy similares. Las únicas variaciones son
debidas a la densidad, por lo tanto la velocidad de salida es ligeramente
superior.

80
 Microturbina de generación de electricidad y calor

6. Intercambiadores de calor.
Para el diseño del intercambiador de calor agua/gas, al igual que en el caso del
recuperador, se basara en un diseño de tubos aleteados, ya que
constructivamente son bastante simples, económicos y relativamente sencillos
de construir.

La metodología propuesta ha sido validada para aplicarse en superficies

extendidas o no, que operan con gases que pueden contener partículas así
como gases limpios o líquidos.

En el proyecto se tendrá en cuenta el cálculo de la resistencia hidráulica interna

del fluido de trabajo en los bancos de tubos, el cálculo del régimen térmico de los
tubos aleteados y el cálculo mecánico. Todos ellos basados en
recomendaciones y normas establecidas para realizar los cálculos térmicos e
hidráulicos de cualquier intercambiador.

La metodología utilizada permite calcular valores del coeficiente de convección

en el intervalo de numero de Reynolds= 3000 – 200000, con tuberías en
triangulo y en cuadro y con un coeficiente de aleteado que va de 1,7 a 6,5. Por lo
tanto podemos abarcar el problema siguiendo un procedimiento relativamente
simple y sencillo con resultados óptimos.

6.1. El uso de los tubos aleteados.

Los tubos aleteados son el elemento principal de cualquier intercambiador tanto

gas/gas como gas/liquido. Dependiendo de los requisitos de cada equipo, se
escoge el tipo de tubo aleteado. Actualmente los más utilizados son los tubos
circulares con aletas en espiral.

Figura 48. Distribución de los tubos.

81
David Espada García

Figura 49. Tipos de tubos.

6.2. Intercambiador agua/gas.

Al igual que en el recuperador, en esta ocasión se aprovechan el calor residual

de los gases de escape (por debajo de 300ºC) que ya han calentado el aire que
sale del compresor y que ahora tienen la función de generar el calor de proceso
necesario para obtener el mayor rendimiento a la turbo máquina.

En nuestro caso la temperatura de los gases en esta zona se encuentra

alrededor de los 270ºC y a una presión casi atmosférica, ya que después de la
turbina, se descarga al exterior.

82
 Microturbina de generación de electricidad y calor

Los pasos para calcular el intercambiador paso a paso, teniendo en cuenta la

perdida de carga, la transferencia de calor y las dimensiones de cada uno de los
componentes del intercambiador, son las siguientes:

6.2.1. Datos de entrada.

Tg1 95 ºC 368,15 k
Tg2 40 ºC 313,15 k
Ta1 25 ºC 298,15 k
Ta2 X ºC X k

mg 2,7 kg/s CPg 1010

ma 3 kg/s

Pa1 0,2 Mpa

Perdidas aire 0,001 Mpa
Pg1 0,5 Mpa
Perdidas agua 0,01 Mpa
Dimensiones 0,6 (a) m
del conducto 0,5 (b) m
flujo gas 1 (c) m

En los datos de entrada, imponemos unas condiciones de entrada y salida de

temperaturas del gas y la temperatura de entrada del agua. A partir de los datos
obtenemos la temperatura de salida del agua y la construcción necesaria para
que se cumpla.

6.2.2. Balance térmico.

En primer lugar realizamos un balance térmico, fijando una temperatura de

entrada y salida de los gases de escape. Este será el modo de calcular, es decir,
fijamos unas condiciones y diseñamos el intercambiador para que se cumplan
dichas condiciones.

Por lo tanto, aplicando la ecuación de balance de calor:

𝑄 = 𝑚𝑔 · 𝐶𝑝𝑔 · 𝑇𝑔1 − 𝑇𝑔2 = 173.3 𝑘𝑊

Igualando el calor obtenido al calor absorbido por el agua, obtenemos la entalpia

de salida del agua mediante la ecuación:

𝑄
𝑕𝑎2 = + 𝑕𝑎1 = 2941758 𝐽/𝑘𝑔
𝑚𝑎

83
David Espada García

Con la entalpia ha2, buscamos en la tabla de las propiedades del agua y

encontramos una temperatura Ta2 = 51ºC.

Llegados a este punto, a modo de simplificación para el cálculo del

intercambiador, usaremos las temperaturas medias tanto del gas como para el
agua:

𝑇𝑚𝑔 = 0.5 · 𝑇𝑔1 + 𝑇𝑔2 = 180º𝐶

𝑇𝑚𝑎 = 0.5 · 𝑇𝑎1 + 𝑇𝑎2 = 40.5º𝐶

Con estos valores medios, buscamos sus propiedades que se exponen en la

siguiente tabla:

densidad m 994,59 kg/m3 0,7833 kg/m3

vol. M 0,0010076 m3/kg 1,27665007 m3/kg
viscosidad m 0,000000643 m2/s 0,00002886 m2/s
conduct. M 0,628 W/mk 0,03707 W/mk
Pr m 4,34 0,683
* En función de P y T CP m 1010 J/KgK
difusión 0,00004222 m2/s

A partir de este momento disponemos de los datos necesarios para el cálculo y

dimensionado del intercambiador.

En un primer momento se desprecia el rendimiento al hacer el balance, ya que

normalmente tienen rendimientos elevados y en los cálculos más precisos
siguientes se tiene en cuenta.

6.2.3. Características constructivas del intercambiador.

El diseño del intercambiador está compuesto por un banco de tubos aleteados y

en una disposición triangular, ya que se favorece a la turbulencia y por lo tanto a
un mayor intercambio de calor.

84
 Microturbina de generación de electricidad y calor

Figura 50. Diseño final de la caja porta tubos.

Los tubos se agrupan en serpentines longitudinales en la dirección del flujo del

aire (no gas). Los codos del serpentín se ubican fuera de los límites para el
conducto del aire. En nuestro caso no existen codos.

Se ha escogido un sistema de tubo aleteado circular de construcción bimetálica.

En la parte interior, los tubos están fabricados de acero al carbono ideal para
generadores de vapor. En la zona exterior, aluminio con aletas en espiral.

Figura 51. Tubo aleteado seleccionado.

Las dimensiones del tubo son las siguientes:

Tubo interno Tubo exterior con aletas

material acero al carbono material aluminio
d. interior 0,02 m diámetro 0,028 m
espesor 0,0025 m altura aleta 0,0135 m

85
David Espada García

paso aletas 0,003 m

espesor aleta 0,0008 m
resistencia 0,000189 m2 k/W

6.2.4. Características geométricas de los tubos aleteados.

Para el cálculo de la superficie del intercambio de calor, se deja de lado la zona

de los codos. Por lo tanto

𝐴𝑡𝑜𝑡 = 𝐴

Según la fórmula que aplicamos, en un metro de de longitud de tubo, la

superficie total será:

𝜋 2 𝐿𝑎
𝐴𝐴1 = 𝐷 − 𝑑2 + 2𝐷𝛿𝐴 𝑧 = 1.219 𝑚2
2 𝑆𝑎

En el caso contrario, la parte de tubo que no tiene aleta por metro es:

𝛿𝑎
𝐴 𝑇1 = 𝜋 · 𝑑 𝐿𝑎 1 − 𝑧 + 𝐿𝑡 = 0.0645 𝑚2
𝑆𝑎

Por lo tanto sumando ambas áreas, obtenemos el área total de tubo por 1 metro
de longitud:

𝐴1 = 𝐴𝐴1 + 𝐴 𝑇1 = 1.284 𝑚2

Para hacer comparaciones entra intercambiadores o tener una idea de relación

de superficies que existen en uno, es ideal calcular las áreas relativas.

Relación de Área aleteada:

𝐴𝐴1
𝑅𝐴𝑎 = 𝐴1 = 0.9497

Relación de Área de tubo sin aleta:

𝐴 𝑇1
𝑅𝐴𝑡 = 𝐴1 = 0.0502

Otros cálculos de área necesarios para el cálculo del dimensionado e

intercambio de calor del aparato, son las siguientes:

Área total de la superficie del tubo en contacto con la aleta, es decir, superficie
exterior de tubo, para 1 metro (La=1):

𝐴𝑠𝑜𝑠𝑡 = 𝜋 · 𝑑 · 𝐿𝑎 = 0.0879 𝑚2

86
 Microturbina de generación de electricidad y calor

El área interna del turbo, por donde circula el agua, para 1 metro:

𝐴𝑖𝑛𝑡 = 𝜋 · 𝑑𝑖𝑛𝑡 · 𝐿𝑎 = 0.0628 𝑚2

Una vez tenemos los valores de áreas calculados, se procede al cálculo del
coeficiente de aleteado:

𝐴1
𝜑𝐴 = 𝐴𝑠𝑜𝑠𝑡 = 14.596

Y para hacer los cálculos de balances y dimensionado, calculamos la relación de

área total externa, teniendo en cuenta el área de aleta y el área interior:

𝐴1
𝑅𝑎𝑟𝑒𝑎𝑠 = 𝐴𝑖𝑛𝑡 = 20.435

Llegados a este punto, podemos pasar al siguiente punto y empezar a

dimensionar el interior del intercambiador.

6.2.5. Dimensiones y paso en el conducto de aire.

Considerando las limitaciones en las dimensiones del intercambiador se elige la

simplificación de “lados iguales” con la densidad máxima, en el cual los pasos
transversal y longitudinal entre tubos tienen la siguiente relación:

3
𝑆2 = 𝑆
2 1

Siendo S2 y S1, los lados longitudinal y transversal respectivamente tal y como

marca la siguiente imagen:

Figura 52. Sistema seleccionado.

Se elige esta distribución, porque de este modo se incita al flujo a que sea
turbulento y por lo tanto a aumentar la transferencia de calor entre los materiales.

87
David Espada García

Las dimensiones que seleccionamos para nuestro intercambiador son las

siguientes:

Dimensiones 0,6 (a) m

del conducto 0,5 (b) m
flujo gas 1 (c) m

La magnitud del paso transversal entre tubos S1, se elige considerando la

dimensión transversal máxima del tubo aleteado, es decir, el diámetro extremo
de la aleta:

𝐷 = 𝑑 + 2 · 𝑙𝑎 = 0.055 𝑚

Otra de las premisas del diseño, consiste en respetarla cantidad de tubos (z1) en
todas las filas transversales del banco:

𝑎
𝑆1 = = 0.0667 𝑚
𝑧1 + 0.5

En nuestro diseño, fijamos una z1=10, que será el número máximo de tubos
aleteados que se acomodan con una separación tecnológica suficiente en un
ancho predefinido “a”.

El paso longitudinal entre los tubos que considera la relación anterior, es:

3
𝑆2 = 𝑆 = 0.0577 𝑚
2 1

Por lo tanto el paso diagonal, usando el método de “lados iguales”, equivale a:

𝑆′2 = 𝑆1 = 0.0667 𝑚

Los pasos relativos entre tubos y la relación entre paso transversal y longitudinal:

𝑆1
𝜍1 = = 2.38
𝑑

𝑆2
𝜍2 = = 2.062
𝑑

𝑆′2
𝜍′2 = = 2.38
𝑑

𝜍1 𝑆1
𝑅𝑒𝑙. 𝑑𝑒 𝑝𝑎𝑠𝑜𝑠 = = = 1.155
𝜍2 𝑆2

88
 Microturbina de generación de electricidad y calor

Una vez llegados a este punto, se procederá a calcular la velocidad y el área

libre del paso del aire, según las ecuaciones específicas para tuberías en
posición triangular.

En primer lugar se calcula el diámetro relativo del tubo aleteado:

2 · 𝑙𝑎 · 𝛿𝑎
𝑑𝑟𝑒𝑙 = 𝑑 + = 0.0352 𝑚
𝑆𝑎

Y seguido, el parámetro del banco de tubos con la siguiente fórmula:

𝑆1 − 𝑑𝑟𝑒𝑙
𝜑𝑟𝑒𝑙 = =1
𝑆′2 − 𝑑𝑟𝑒𝑙

Cuando 𝜑𝑟𝑒𝑙 es ≤ 2el área mínima libre se ubica en el plano del paso transversal
y su valor es:

𝐹 = 𝑎 · 𝑏 − 𝑧1 · 𝐿𝑠𝑐 · 𝑑𝑟𝑒𝑙 = 0.2436 𝑚2

Aquí la longitud de los tubos en los límites del conducto del aire se toma igual a
su altura:

Lsc=b=0.7 m

Con el valor promedio del volumen especifico del aire con las propiedades a
Temp. gas medio, podemos calcular la velocidad del aire con la siguiente
expresión:

𝑚𝑔 · 𝑉𝑜𝑙. 𝐸𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐𝑜 𝑚
𝑢𝑔 = = 4.09
𝐹 𝑠

Ahora se procede al cálculo idéntico pero para el interior del tubo, es decir, para
la velocidad del agua que circula en el interior.

Para ello, necesitaremos calcular Zp, que equivale al número de tubos

conectados en paralelo, que se determinara por el número de tubos en la fila
transversal del banco de tubos z1 y por el número de vueltas de los serpentines
nx. Para la simplificación en tubos en posición triangular, se adopta que:

𝑛𝑥 ≥ 2

Por lo tanto, tomamos el valor de nx=2

𝑍𝑝 = 𝑛𝑥 · 𝑧1 = 20

2
𝜋 · 𝑑𝑖𝑛𝑡
𝑓 = 𝑍𝑝 · = 0.00628 𝑚2
4

89
David Espada García

La velocidad media a la que fluye el agua en el interior del tubo será:

𝑚𝑎 · 𝑉𝑜𝑙. 𝐸𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐𝑜 𝑚
𝑢𝑎 = = 0.32
𝑓 𝑠

6.2.6. Calculo de la superficie de intercambio de calor.

El área de la superficie de intercambio, se puede encontrar mediante la siguiente

ecuación:

𝑄 · 1000
𝐴=
𝑈 · ∆𝑇

- U  Es el coeficiente global de transmisión de calor y e determina

mediante la fórmula:

𝜑
𝑈=
𝐴 1 𝐴 1
· + · 𝑅𝑡 +
𝐴𝑖𝑛𝑡 𝑕2 𝐴𝑖𝑛𝑡 𝑕1𝑟𝑒𝑙

Donde ψ es el coeficiente de efectividad térmica y se suele tomar un valor de

0,95, para un flujo guiado y limpio. Seria 0,85 si las características del flujo se
viesen afectadas de alguna manera.

Para determinar la resistencia térmica de la pared de los tubos bimetálicos las

magnitudes δ y k, se pueden despreciar debido a que son muy pequeñas y por
lo tanto podemos tomar Rt=Rc, es decir, la resistencia total es prácticamente la
resistencia térmica por contacto.

Para nuestro caso Rc = 0.000189 m2K/W

Como el tipo de parámetro escogido es respecto a flujo de gas limpio, las

ecuaciones deben de ir referenciadas a lo mismo, así, por lo tanto calcularemos
h1rel y h2, mediante las siguientes ecuaciones.

Para calcular h1rel:

𝐴𝐴 𝐴𝑇
𝑕1𝑟𝑒𝑙 = · 𝐸 · 𝜇𝑎 · 𝜑𝑒 + 𝑕𝑐
𝐴 𝐴

Donde, hc se calcula mediante las siguientes ecuaciones:

- Parámetro de forma del banco de tubos, que para tubos en triangulo es:

𝜍1 1.26
𝑋= − − 2 = −0.9316
𝜍2 𝜑𝑎

90
 Microturbina de generación de electricidad y calor

- Exponentes n y Cq, se calculan como:

𝑛 = 0.7 + 0.08 · tanh 𝑋 + 0.005 · 𝜑𝑎 = 0.7144

1.1
𝐶𝑞 = 1.36 − tanh 𝑋 · − 0.014 = 0.0725
𝜑𝑎 + 8

- Consideramos el numero de filas del banco z2 ≥ 10, entonces el

Coeficiente Cz = 1 (Según tabla).

Figura 53. Gráfica para el cálculo de Cz.

Con estos tres parámetros calculados, podemos obtener hc, de la siguiente

ecuación:
𝑛
𝐾𝑔 𝑢𝑔 · 𝑑 𝑊
𝑕𝑐 = 1.13 · 𝐶𝑧 · 𝐶𝑞 · · · 𝑃𝑟𝑔 0.33 = 35.62
𝑑 𝑣𝑔 𝑚2 𝐾

Obtenido ya hc, pasamos a calcular el coeficiente de eficiencia teórica de la

aleta.

El valor de tal coeficiente, se determina multiplicando E·μa·ψE. Para determinar

E, que es el coeficiente de eficiencia teórica de aleta, se necesita calcular el
valor del parámetro de la aleta m. Para ello, a su vez, requiere conocer la
conductividad térmica de la aleta correspondiente a su temperatura media Ta,
evaluando Ta de acuerdo con la siguiente ecuación:

𝑇𝑎 = 𝑇𝑚𝑔 − 𝑇𝑚𝑔 − 𝑇𝑚𝑎 · 𝐸 ′ = 39.29º𝐶

91
David Espada García

Como depende de E, usamos un valor inicial E’=0,9. El resultado debe de ser

muy parecido al supuesto, sino se debería de corregir hasta hacer converger.

Figura 54. Gráfica para el cálculo de E.

Una vez calculada Ta, calculamos el factor m, que depende de las propiedades
del material (conductividad térmica), según a la temperatura a la que se
encuentra:

2 · 𝑕𝑐
𝑚= = 19.85
𝛿𝑎 · 𝐾𝑎

En este caso tenemos una conductividad térmica del aluminio de Ka=225 W/mK,
para la Ta calculada.

Se calcula la altura relativa de la aleta, según la siguiente expresión:

𝐷 𝐷
𝑙 ′ 𝑎 = 𝑙𝑎 · 1 + 0.191 + 0.054 · · ln = 0.0162 𝑚
𝑑 𝑑

92
 Microturbina de generación de electricidad y calor

Una vez llegados aquí, podemos calcular E, el coeficiente de aleta teórica de la

aleta:

(𝑚 · 𝑙 ′ 𝑎)
tanh⁡
𝐸= = 0.967
𝑚 · 𝑙′ 𝑎

Como el valor calculado de E≈E’, no hace falta re calcular nada ni variar Ka.

Podemos pasar al siguiente punto y calcular el coeficiente de corrección ψE:

𝐷
𝜑𝐸 = 1 − 0.016 · − 1 · 1 + 𝑡𝑎𝑛𝑕 2𝑚 · 𝑙𝑎 − 1 = 0.9717
𝑑

Cuando tenemos aletas de espesor constante, el valor de μa es 1.

Substituyendo los valores calculados anteriormente así como los de relaciones

de áreas con y sin aletas, de antes, podemos obtener el valor de coeficiente de
convección relativo h1rel:

𝐴𝐴 𝐴𝑇 𝑊
𝑕1𝑟𝑒𝑙 = · 𝐸 · 𝜇𝑎 · 𝜑𝐸 + · 𝑕𝑐 = 33.57 2
𝐴 𝐴 𝑚 𝐾

Calculo de h2

Para la transferencia de calor hacia el fluido interior, se calculara el coeficiente

de transferencia h2. Para ello necesitaremos calcular Reynolds interior, Prandtl,
viscosidad del fluido interior (ya buscada anteriormente) y la velocidad del fluido.

𝑈𝑓 · 𝑑𝑖𝑛𝑡
𝑅𝑒𝑓 = = 9976
𝑣𝑓

Siendo uf, la velocidad media del agua dentro de los tubos, que para nuestros
valores y cálculos, obtenemos un resultado de 0.32 m/s.

Por lo tanto, aplicando los valores obtenidos de Ref y Prf, para calcular el
coeficiente de convección h2, es necesario buscar la magnitud de K y ζ, que
intervienen en la expresión.

900
𝐾 =1+ = 1.09
𝑅𝑒𝑓

𝜁 = (1.82 · 𝑙𝑜𝑔 𝑅𝑒𝑓 − 1.64)−2 = 0.0315

Además, necesitaremos buscar el factor de corrección Ctem que se evalúa

mediante la temperatura media de la superficie interna del tubo Tw, dando
valores que posteriormente se corrigen, del área de la superficie externa de
intercambio de calor A’=95 m2 y del coeficiente de convección h’2=3500 W/m2K.

93
David Espada García

Por lo tanto considerando estos valores y aplicándolos a las siguientes

ecuaciones:

𝐴′
𝐴′𝑖𝑛𝑡 = = 4.649 𝑚2
𝐴
𝐴𝑖𝑛𝑡

𝑄 103
𝑇′ 𝑤 = 𝑇 + = 51.15 º𝐶
𝐴′𝑖𝑛𝑡 𝑕′2

Para el caso de estudio, cuando el fluido de trabajo es agua y tiene lugar su

calentamiento (el flujo de calor se dirige de la pared al agua) el factor de
corrección Ctem se determina con la siguiente fórmula:

0.11
𝜇𝑓
𝐶𝑡𝑒𝑚 =
𝜇𝑤

Para la presión del agua ya impuesta al principio del apartado y la temperatura

media del fluido, aplicamos su viscosidad ya calculada y para la T’w calculada
buscamos su viscosidad, tal que:

𝜇𝑤 = 1.8 · 10−5 𝑃𝑎 · 𝑠

Por lo tanto, aplicando la ecuación de Ctem:

𝐶𝑡𝑒𝑚 = 0.693

Substituyendo todos los valores obtenidos anteriormente en la para el cálculo de

h2:

𝑘𝑓 0.125 · 𝜁 · 𝑅𝑒𝑓 · 𝑃𝑟𝑓 · 𝐶𝑡𝑒𝑚 𝑊

𝑕2 = = 1535.25
𝑑𝑖𝑛𝑡 𝜆 · +4.5 · 𝜁 0.5 · 𝑃𝑟𝑓 0.666 − 1 𝑚2 𝐾

Por lo tanto regresando al principio del apartado, con todos los valores en la
mano, podremos calcular el coeficiente global de transferencia de calor:

𝜑 𝑊
𝑈= = 20.23 2
𝐴 1 𝐴 1 𝑚 𝐾
· + · 𝑅𝑡 +
𝐴𝑖𝑛𝑡 𝑕2 𝐴𝑖𝑛𝑡 𝑕1𝑟𝑒𝑙

Diferencia de temperaturas medias

∆𝑇𝑚𝑎𝑦𝑜𝑟 = 239º𝐶

∆𝑇𝑚𝑒𝑛𝑜𝑟 = 21.15º𝐶

Por lo tanto la temperatura media de ellas:

94
 Microturbina de generación de electricidad y calor

∆𝑇𝑚𝑎𝑦𝑜𝑟 − ∆𝑇𝑚𝑒𝑛𝑜𝑟
∆𝑇 = = 89.84º𝐶
∆𝑇𝑚𝑎𝑦𝑜𝑟
ln ∆𝑇
𝑚𝑒𝑛𝑜𝑟

6.2.7. Resultados.

Llegados a este punto, calcularemos los resultados para obtener las

dimensiones de nuestro intercambiador, así como número de tubos, etc.

En primer lugar calculamos el área de superficie de intercambio:

𝑄 · 103
𝐴= = 95.35 𝑚2
𝑈 · Δ𝑇

Se comprueba que no hay variación en Ctem, y por lo tanto, podemos usar el

mismo valor y seguir calculando.

𝐴
𝐴𝑖𝑛𝑡 = = 4.67𝑚2
𝐴
𝐴𝑖𝑛𝑡

𝑄 103
𝑇𝑤 = 𝑇 + = 64.7º𝐶
𝐴𝑖𝑛𝑡 𝑕2

Como se puede ver, Tw≈T’w, lo cual quiere decir que Ctem no requiere
corrección ya que depende de ese valor.

Por lo tanto seguimos con el cálculo, y buscamos el número total de tubos del
intercambiador:

𝐿𝑎𝑙
𝑧= = 106.09
𝐿𝑠𝑐

Numero de filas de tubos transversales:

𝑧
𝑧2 = 10.609
𝑧1

Se toma el valor inmediatamente superior de z2, es decir z2=11

Profundidad del conducto de gases:

𝑐 = 𝑧2 − 1 𝑆2 = 0.5547

Número real de tubos del intercambiador:

𝑍 ′ = 𝑧1 · 𝑧2 = 106.9

Longitud real de los tubos aleteados del intercambiador:

95
David Espada García

𝐿′ 𝑎𝑙 = 𝐿𝑠𝑐 · 𝑍 ′ = 74.26 𝑚

6.2.8. Calculo de caída de presión externa.

El cálculo de la caída de presión, lo realizaremos según la longitud relativa de la

𝐴
superficie extendida 𝑡𝑜𝑡 𝐹 y el diámetro equivalente de la sección más estrecha
del banco de tubos deq .

𝐴𝑡𝑜𝑡
La magnitud 𝐹 se calcula según la fórmula:

𝐴𝑡𝑜𝑡 𝜋 · 𝑑 · 𝑆𝑎 + 2 · 𝑙𝑎 · 𝛿𝑎 + 2 · 𝑙𝑎 · (𝑙𝑎 + 𝑑)
= = 40.8
𝐹 𝑆1 · 𝑆𝑎 − (𝑑 · 𝑆𝑎 + 2 · 𝑙𝑎 · 𝛿𝑎)

El diámetro equivalente para los bancos de tubos con una distribución en

triangulo con υ ≤ 2, se determina según:

2 · 𝑆𝑎 · 𝑆1 − 𝑑 · 2 · 𝑙𝑎 · 𝛿𝑎
𝑑𝑒𝑞 = = 0.00629 𝑚
2 · 𝑙𝑎 + 𝑆𝑎

Para el cálculo del coeficiente de perdidas ζo, se utiliza la siguiente ecuación,

que es válida en el caso de un banco de tubos con posición en triangulo. Se
determinan las magnitudes del exponente n, del coeficiente Cr y el C’z, que
intervienen en la formula.

0.25 0.57
𝐴 𝑆1 𝑆1
𝑛 = 0.17 · · · 𝑒 −0.36· 𝑆2 = 0.308
𝐹 𝑆2
0.53 1.3
𝐴 𝑆1 𝑆1
𝐶𝑟 = 2.8 · · · 𝑒 −0.9· 𝑆2 = 8.525
𝐹 𝑆2

Como z2 ≥ 8, tomamos C’z = 1.

De acuerdo con la ecuación, se determina el coeficiente de pérdidas del banco

de tubos:
−𝑛
′
𝑢𝑔 · 𝑑𝑒𝑞
𝜁𝑜 = 𝐶 𝑧 · 𝐶𝑟 · = 0.906
𝑣𝑔

Llegados a este punto, con la siguiente formula, podemos calcular la caída de

presión externa en el intercambiador:

𝜌𝑔 · 𝑢𝑔2
Δ𝑃 = 𝐶𝑜𝑝 · 𝜁𝑜 · 𝑧2 · = 0.07 𝑃𝑎
2

96
 Microturbina de generación de electricidad y calor

6.3. Recuperador turbo máquina.

Los recuperadores de calor aprovechan la energía térmica en los gases de

escape de la combustión, en este caso de nuestra turbina de gas.

El procedimiento de cálculo, es el mismo que en el intercambiador de calor agua-

gas. El motivo, es para unificar los sistemas de intercambio y hacer un mismo
diseño.

Las ventajas de aprovechar un mismo tipo de intercambiador son básicamente

constructivas y económicas. Se pueden pedir más metros de un mismo tubo en
el cual pueda funcionar con agua o aire por igual y la caja del intercambiador
puede ser similar, con pocas variaciones.

Por lo tanto, en el cálculo se simplificara todo el proceso por haber sido ya

explicado previamente.

6.3.1. Datos de entrada.

Datos de entrada
Tg1 590 ºC 863,15 k
Tg2 290 ºC 563,15 k
Ta1 196 ºC 469,15 k
Ta2 X ºC X k

mg 0,78 kg/s CPg 1150

ma 0,76 kg/s

Pg1 0,15 Mpa

Perdidas gas 0,001 Mpa
Pa1 0,45 Mpa
Perdidas aire 0,002 Mpa
Dimensiones 0,6 (a) m
del conducto 0,5 (b) m
flujo gas 1 (c) m

6.3.2. Balance térmico.

Potencia calorífica aportada por el salto térmico del gas de escape:

97
David Espada García

Q=269.1 kW

Obtenemos una entalpia de salida del aire de h2=3209.2 kJ/kg

Por lo tanto una temperatura de salida del aire de 435ºC

Las temperaturas del aire y gas medias son 316 y 440ºC respectivamente y sus
propiedades a esas temperaturas son:

Propiedades físicas del aire a MEDIA: Prop. Físicas del gas m

densidad m 1,83 kg/m3 1,506 kg/m3
vol. M 0,54638 m3/kg 0,66405 m3/kg
viscosidad m 0,00001756 m2/s 0,00002886 m2/s
conduct. M 0,02932 W/mk 0,03707 W/mk
Pr m 0,701 0,683
* En función de P y T CP m 1100 J/KgK
difusión 0,00004222 m2/s

6.3.3. Características constructivas del intercambiador.

Utilizaremos el mismo tubo exactamente que en el caso anterior, para unificar

productos. La ventaja es reducir el coste.

La única variación tiene que ver con el funcionamiento de los flujos que se
mueven dentro, que están en otras temperaturas y presión, por lo tanto, cambia
la velocidad de cada uno de ellos

Por otro lado, el valor de Z cambia, ya que hay más tubos para intercambiar
calor, porque el aire es peor conductor del calor que el agua.

Asi por lo tanto Z=15 y obtenemos una velocidad media del gas de 4.3 m/s y una
velocidad media del aire de 30,43 m/s.

6.3.4. Calculo del área de la superficie de intercambio.

Las variables como X, n, Cq, que dependen de valores constructivos, siguen

igual, pero varía lo siguiente:

hc=36.95 W/m2K

98
 Microturbina de generación de electricidad y calor

En cuanto al apartado del coeficiente de eficiencia teórica de la aleta, tenemos

que, la temperatura media de intercambio es:

Ta=306.04ºC

En este caso la variación de la conductividad del aluminio debido a la

temperatura, es insignificante, por lo tanto, cogemos el mismo valor: Ka=225
WmK.

El cálculo del coeficiente de aleta, nos da un resultado muy similar al supuesto

anteriormente, por lo tanto, continuamos con los cálculos:

𝐸 ≈ 𝐸′ 0.95 ≈ 0.9

Asi por lo tanto, el coeficiente de corrección calculado será de:

υe= 0.9715

Con todo ello el valor de coeficiente de convección relativo h1rel, que será:

h1rel=34.8 W/mK

Para el flujo interior, el número de Reynolds se ve incrementado debido al paso

del aire a alta velocidad por el interior de los tubos:

Ref=34667

Con ello, calculamos varios parámetros como K,ζ y suponemos un area A’, de 95
m2, como en el anterior caso, ya que nos interesa mantener un tamaño
determinado de intercambiador.

Con todo esto, encontramos el area interior A’int, que será 4.65 m2 y una
temperatura de pared T’w de 332ºC. El factor de corrección Ctem, será 0.982.

Aplicando los datos y parámetros calculados, podemos acceder al cálculo del

coeficiente de convección interno h2.

h2= 112.95 W/m2K

Por lo tanto aplicando las ecuaciones anteriores, encontramos un coeficiente

global de transferencia de calor U, que tendrá un resultado de: 4.45 W/m2K.

La diferencia de temperaturas medias será:

ΔT mayor=155ºC.

ΔT menor=136ºC.

99
 David Espada García

El cálculo de temperatura media logarítmica será

ΔT=145ºC.

Por lo tanto, el cálculo de area total de intercambio, será:

Area=416.25 m2.

6.3.5. Cálculos dimensionales.

Se adjunta la tabla del Excel con los resultados dimensionales finales:

Aint 20,369333 m2
Tw 432,4600855 ºC * Si Tw ~ T'w --> Ctem no requiere corrección
Lal 324,1880035 m * Longitud de tubos con aletas
z 463,1257193 num. Tubos
Z2 30,87504795 numero filas * cogemos superior
C 1,168437763 m * Profundidad conducto gases
Zr 463,1257193 * Numero de tubos en intercambiador
Lr alet. 324,1880035 m * metros de tubo aleteado dentro

Habrá un total de 432 tubos que equivalen a 15 columnas y 31 filas

aproximadamente, siempre cogiendo el número superior.

6.3.6. Caída de presión externa.

El resultado es prácticamente insignificante para los cálculos. Es de 1.15 Pa, de

caída de presión exterior.

100
 Microturbina de generación de electricidad y calor

7. Estructura de la maquina.

7.1. Zona de generación eléctrica y admisión.

Esta sería la zona principal de la maquina. El corazón de la misma. En este
lugar, el compresor, se encarga de absorber aire y comprimirlo, para después
calentarlo haciéndolo pasar por el recuperador y finalmente para por la cámara
de combustión que genera los productos a alta temperatura que servirán para
que la turbina genere energía.

Figura 55. Zona principal.

De izquierda a derecha, tenemos la caja de filtro y la admisión con el generador.

Seguidamente compresor y cámara de combustión y finalmente la turbina
conectada al recuperador.

7.2. Zona de intercambio de calor.

El Intercambiador de tubos consta de 2 zonas independientes pero unidas en un
mismo cuerpo.

Figura 56. Zona térmica de la máquina.

En la imagen está abierto para poder diferenciar ambas zonas. La grande de la

izquierda, es el recuperador. La destapada, seria la zona de intercambiador agua
gas.

101
David Espada García

7.3. Zonas de componentes periféricos.

Debajo de la maquina, quedan varios huecos libres. En ellos se han situado
algunos elementos.

Figura 57. Zona general.

En la zona central, irían situados todos los elementos auxiliares, es decir,

elementos que trabajan para el resto, como por ejemplo el sistema de lubricación
de la maquina. Podemos ver la bomba de aceite, el radiador y el ventilador. Todo
ello para mantener unas condiciones adecuadas del aceite.

Tambien se ve un ventilador. En este caso es de refrigeración de todos los

componentes eléctricos y/o electrónicos.

El espacio también se puede usar para almacenar recambios o herramientas

necesarias para el mantenimiento de la maquina.

Figura 58. Zona global de control y potencia.

En ambos laterales se han situado los elementos electrónicos (izquierda) de

control y los eléctricos de potencia a la derecha. El objetivo de ambos, es
controlar el estado y funcionamiento de la maquina mientras que el otro se
encarga de dar energía a los componentes y transformar la energía obtenida del
generador.

102
 Microturbina de generación de electricidad y calor

8. Presupuesto y planificación.
Planificación:

103
 David Espada García

Presupuesto del proyecto de ingeniería.

El presupuesto, va en función de las horas invertidas en cada uno de los

apartados del proyecto, sin tener en cuenta la fabricación.

Pre-proyecto FASE 1 Calculos/Diseño previo FASE 2

Tarea Sub-tarea Horas Tarea Sub-tarea Horas

Informacion 6 Calculos generales bypass
Marco legal 4 Compresor 1
TOTAL 1 10 Turbina 1,5
Intercamb. 0,5
Recuper. 1,5
C.Comb. 2
TOTAL 6,5
Pre-diseño
Componentes 0,5
Distribucion 0,5
medidas 0,5
TOTAL 1,5

TOTAL 2 8

Diseño y calculo FASE 3 Proyecto FASE 4

Sub-
Tarea Sub-tarea Horas Tarea tarea Horas
Compresor Memoria 100
Calculo
velocidades 3 Condiciones tecnicas/Anexo 20
Geometria
alabes 3 Planos 60
Diseño rotor 8 Presentacion Maquina 20
Dimensionado
eje 1 TOTAL 4 200
Dimensionado
tobera 2
Diseño voluta 10
Diseño cuerpo 4
TOTAL 31
Generador

104
 Microturbina de generación de electricidad y calor

S/ Espec. ABB 1
TOTAL 1
Cam. Comb.
Calculo proceso 8
Seleccion
Inyector 2
Diseño de
camara 10
TOTAL 20
Turbina
Calculo
velocidades 4
Geometria
alabes 4
Diseño rotor 6
Dimensionado
eje 0,5
Dimensionado
tobera 5
Diseño voluta 8
Diseño cuerpo 0,5
TOTAL 28
Recuperador
Caudales 2
Calculo forma 6
Diseño 6
TOTAL 14
Intercambiador
Caudales 1,5
Calculo forma 5
Diseño 5
TOTAL 11,5
Estructura
Calculo
estructura 3
Diseño partes 8
Componentes 10
TOTAL 21
Gases/Liquidos
Instalacion 5
Diseño
recipientes 4
Ubicación 2
Tuberias 5
TOTAL 16
Sistema elect.

105
 David Espada García

Controles
varios 6
Alimentacion 6
TOTAL 12

TOTAL 3 154,5

Horas totales proyecto 372,5 h

Coste ingenieria 60 €/h
Coste total 22350 €

106
 Microturbina de generación de electricidad y calor

9. Consideraciones ambientales.
9.1. Etapa de fabricación
La maquina está compuesta en el 90% de componentes metálicos, ya sean
aluminio o aceros. En estos casos, la mayoría de los procesos de fabricación de
estas piezas son mediante arranque de viruta. También existen piezas plasticas
de los rodillos que son fabricados por el mismo método que los metales o no.

En estos casos, la viruta es separada para cada material, recogida y enviada al

punto de reciclado (“deixalleria”). Tanto el aluminio, como el acero, como algunos
tipos de plásticos, son productos reciclables.

Por último, comentar la necesidad de reciclar los envoltorios de los componentes

de serie que se adquieren (tornilleria, etc.), así como elementos de cartón,
mediante los contenedores adecuados para los diferentes residuos.

9.2. Etapa de servicio

Durante el servicio de la maquina en el, es indispensable hacer los cambios de
aceite estipulados. El sistema es igual que un cambio de aceite de un vehículo
convencional o maquina industrial y el aceite residual debe de llevarse al punto
de reciclaje idóneo para este tipo de residuos. El aceite es muy contaminante y
no debe de desecharse en cualquier lugar.

9.3. Fin de vida de la maquina

Una vez llegados a este punto, la maquina debe ser desmontada y clasificada
según materiales o componentes para desecharlos donde sea conveniente.Los
materiales citados anteriormente, son perfectamente fundibles con nuevos y por
lo tanto son reutilizables.

Proceso de despiece y clasificación de materiales:

- Turbina y alternador, son llevados a puntos de reciclaje, en donde los reparan y

reutilizan o desechan debidamente.

- Lubricante llevado al punto de reciclaje adecuado

- Despiece de todas las partes del intercambiador y planchas y demás soportes.

En el desmontaje, separar toda la tornillería, elementos de acero y planchas, de
los elementos de aluminio del reductor. Llevar ambos materiales por separado al
punto de reciclado idóneo.

107
David Espada García

10. Bibliografía.

Libros.

Society of Automotive Engineers. "Supercharger Testing Standard". Code: J1723,

1995.

Richard G. Budynas, J. Keith Nisbett. “Diseño en ingeniería mecánica de Shigley”

(octava edición). Ed. Mc Graw Hill, 2008.

Virgil Moring Faires. “Diseño de elementos de maquinas (edición 1)”. Ed. Limusa
Noriega editores, 1999.

Dr. M. Muñoz, Dr. F. Payri. “Motores de combustión interna alternativos” Fundación

general – U.P.M. 1989.

Carles Riva i Romeva. “Disseny de maquines IV. Selección de materials 1” Edicions

UPC. 1999.

Emilio Bianco. “Manuale practico lavorazioni meccaniche”. Editorial Rosenberg &

Sellier, 1956.

F. M. Golden, L. Batres, G. Terrones M. “Termofluidos, turbomaquinas y maquinas

térmicas”. Ed. CECSA, 1989.

Olga Alcaraz, Jose Lopez, Vicente Lopez. “FISICA, Problemas y ejercicios resueltos”.
Ed. Pearson Prentice Hall, 2006.

Webs.

http://www.turbobygarrett.com/turbobygarrett/

http://www.turbodriven.com/es/Default.aspx

http://www.skf.com/

http://www.corneplas.com/nylon.php

http://www.transmisionesbayod.es/indice.htm

http://www.rbracing-rsr.com/turbosupermaps.html#kkk

http://www.dotmar.com.au

http://www.fag.de

http://www.rotrex.com

http://www.opac.net

http://www.klueber.com/

http://www.castrol.com

http://www.bosch.com.mx

108