Gasificadores

Trabajo Fin de Grado en
Ingeniera de Tecnologas Industriales

Desarrollo de modelos matemticos y diseo de los
principales equipos de una planta de gasificacin de
biomasa para produccin de electricidad.
Autor
Guillermo Figueruelo Malo
Director
Alberto Gonzalo Callejo
Escuela de Ingeniera y Arquitectura. Universidad de Zaragoza.

Septiembre 2014
TTULO: Desarrollo de modelos matemticos y diseo de los principales equipos

de una planta de gasificacin de biomasa para produccin de electricidad
AUTOR: Guillermo Figueruelo Malo
DIRECTOR: Alberto Gonzalo Callejo
DEPARTAMENTO: Ingeniera Qumica y Tecnologas del Medio Ambiente
Desarrollo de modelos matemticos y diseo de los principales equipos de una

planta de gasificacin de biomasa para produccin de electricidad.
RESUMEN
Este proyecto sigue la lnea marcada por el trabajo de la asignatura Procesos Qumicos
Industriales y que llev por ttulo Diagrama de bloques y balances para una planta de
gasificacin de biomasa en lecho fluidizado para 600kg/h de alimentacin de madera de pino
con el objetivo de producir electricidad en un motor.
El proyecto realizado es eminentemente terico. Se ha generado una herramienta que permite
predecir los productos resultantes de procesar cualquier biomasa va gasificacin.
Tras haber realizado una revisin bibliogrfica y seleccin de datos de partida utilizando
artculos cientficos y bases de datos como Phillis del Energy research Center of the Netherlands
(ECN), se ha desarrollado un modelo terico del gasificador a partir de los balances de materia y
energa y de los principales equilibrios de las reacciones de gasificacin. Para la realizacin de
esta tarea se ha utilizado el software Engineering Equation Solver (EES).
Tambin se ha realizado el dimensionamiento del resto de equipos principales como son el
separador ciclnico o el intercambiador de calor, as como un anlisis del comportamiento
mecnico del gasificador implementado en Abaqus.
Tras la fase de diseo, se ha procedido a la simulacin del sistema completo, tanto en EES como
en Aspen Hysys, comparando los parmetros y resultados obtenidos entre estas dos
simulaciones y resultados experimentales encontrados en la bibliografa revisada. El sistema
completo se trata de una planta de produccin de electricidad a partir de la combustin, en un
motor de combustin interna, del gas de sntesis generado a partir de la gasificacin de 600
kgh de biomasa.
NDICE
1.
OBJETIVO Y ALCANCE DEL PROYECTO.............................................................. 1
2.
ANTECEDENTES ........................................................................................................... 2
3.
4.
2.1
CONTEXTO GLOBAL ........................................................................................... 2
2.2
TIPOS DE BIOMASA ............................................................................................. 2
2.3
VENTAJAS DEL USO ENERGTICO DE LA BIOMASA ............................... 3
2.4
INCONVENIENTES DEL USO ENERGTICO DE LA BIOMASA................ 4
2.5
PROCESO DE GASIFICACIN DE BIOMASA ................................................ 5
EQUIPOS DE LA PLANTA ........................................................................................... 6

3.1
REACTOR DE GASIFICACIN .......................................................................... 6
3.2
SEPARADOR CICLNICO .................................................................................. 7
3.3
INTERCAMBIADOR DE CALOR........................................................................ 7
DESARROLLO DE LOS MODELOS MATEMTICOS ........................................... 9

4.1
4.1.1
REACTOR DE GASIFICACIN ........................................................................ 10

DISEO MECNICO DEL REACTOR DE GASIFICACIN .................................................. 14
4.2
SEPARADOR CICLNICO ................................................................................ 17
4.3
INTERCAMBIADOR DE CALOR...................................................................... 20
5.
RESULTADOS ............................................................................................................... 23
6.
CONCLUSIONES .......................................................................................................... 27
7.
REFERENCIAS ............................................................................................................. 28
8.
ANEXOS ......................................................................................................................... 30
8.1
TABLAS PARA EL CLCULO DE LAS CONSTANTES DE EQUILIBRIO30
8.2
ENTORNO DE ABAQUS ..................................................................................... 31
8.3
ENTORNO DE HYSYS ......................................................................................... 33
8.4
ENTORNO DE EES............................................................................................... 35
1. OBJETIVO Y ALCANCE DEL PROYECTO
1. Objetivo y alcance del proyecto
1. OBJETIVO Y ALCANCE DEL PROYECTO

El Grado en Ingeniera de Tecnologas Industriales tiene una clara singularidad que lo diferencia
del resto de ttulos de grado de Ingeniera por su carcter multidisciplinar e integrador de todas
las tecnologas industriales y su amplia base formativa generalista.
Este Trabajo Fin de Grado intenta seguir la senda marcada durante estos cuatro aos de
formacin y huye de la especializacin en un tema concreto, para abarcar en cantidad y fondo
suficiente, varias de las ramas incluidas dentro del plan de estudios como son la Ingeniera
Qumica, Energas, Ingeniera Trmica o la Ingeniera Mecnica entre otras.
En vista del considerable inters que est adquiriendo mundialmente el proceso de gasificacin
de biomasa, y ms recientemente de residuos, se ha hecho preciso el desarrollo de modelos
matemticos que permitan predecir el comportamiento de los equipos incluidos en dicho
proceso. Igualmente, se pretende realizar el diseo bsico de alguno de los equipos utilizados
para el acondicionamiento del gas producto como el separador ciclnico o el intercambiador de
calor.
2. ANTECEDENTES
2. Antecedentes
2. ANTECEDENTES
2.1
CONTEXTO GLOBAL
El incremento demogrfico y el constante aumento de la produccin de bienes y servicios,

generan un alza sostenida de la demanda de energa a escala global. Si queremos satisfacer esta
creciente necesidad y a la vez respetar el equilibrio de los ecosistemas de la Tierra, resulta
imprescindible investigar nuevas formas de produccin y gestin de la energa que nos permitan
conciliar el crecimiento con el respeto al medio ambiente.
Desde la revolucin industrial, acaecida a finales del Siglo XVIII, la capacidad humana de alterar
su entorno ha crecido exponencialmente. En la ltima dcada, ms de dos siglos despus de que
Watt desarrollase su mquina de vapor, comenzamos a ser conscientes de los problemas
ecolgicos que se han derivado de recorrer la senda industrial. Los combustibles fsiles siguen
siendo actualmente, al igual que en el inicio del proceso de industrializacin, las fuentes bsicas
de energa. La continua emisin a la atmsfera de toneladas de CO2, derivadas de la quema de
combustibles fsiles est dando lugar a una transformacin del clima de la Tierra cuyas
consecuencias ecolgicas, econmicas y sociales son imprevisibles.
El viento, el sol y el oleaje son recursos que hemos aprendido a aprovechar para producir
energa de forma limpia. Pero existe otro recurso que tambin se renueva continuamente, la
biomasa. Del total de las energas renovables, que suponen el 13% del consumo energtico
mundial, el 80% corresponde al consumo de biomasa suponiendo el 10% del consumo
energtico mundial. La siguiente energa renovable ms utilizada es la hidroelctrica con un
16% [1]. En la Comunidad Autnoma de Aragn, se disponen actualmente de aproximadamente
3
de potencia instalada de gasificacin con biomasa, existiendo algunos proyectos en
distintos estados de ejecucin, lo que hace que puedan estimarse para 2020 los 15
de
potencia instalada, lo cual supondra una produccin acumulada de energa a lo largo de todo el
periodo de 94.864
, con 30.000
en 2020 [2].
Fig 1. Prospectiva de potencia instalada y de produccin de energa elctrica en Aragn en

gasificacin con biomasa. Periodo 2013 2020.
2.2
TIPOS DE BIOMASA
La biomasa puede definirse como todo material orgnico no mineralizado que pueda ser usado
para la produccin de energa.
2. Antecedentes
Habitualmente la biomasa es clasificada en dos grandes grupos: primaria y residual [3]. La

primaria se refiere a cultivos energticos, es decir, la produccin y gestin de diversos tipos de
vegetales con el fin ltimo de obtener energa. Con biomasa residual o de proceso, denominamos
a aquella que se extrae del material de desecho resultante de diversas actividades productivas,
ya sean agrcolas, industriales o forestales.
En cuanto al origen y tipos de residuos generados, se distinguen cuatro categoras generales:
2.3
Residuos agrcolas, forestales y cultivos energticos. Los residuos agrcolas se obtienen

de la planta cultivada al separar la parte utilizada por la industria alimentaria de la que
no se destina a ese fin, que es la interesante desde el punto de vista energtico. Por su
parte los residuos forestales provienen del mantenimiento de bosques y de su
explotacin para el uso industrial de la madera. Por ltimo, los cultivos energticos,
clasificados en herbceos, de madera y acuticos; se generan con el nico fin de producir
biomasa transformable en combustible. Es por ello que su rentabilidad energtica y
competencia con el sector alimenticio son bastantes discutidas sobre todo en pases en
vas de desarrollo como Brasil.
Residuos de industrias forestales y agroalimentarias. Estos residuos pueden ser
valorados energticamente tras ser originados en industrias madereras, papeleras,
agroalimentarias
Residuos urbanos. Los denominados RSU (Residuos Slidos Urbanos) se generan en los
ncleos de poblacin como consecuencia de la actividad humana diaria y a pesar de
tener una composicin muy variable, su contenido en materia orgnica suele ser del
orden del 50%.
Residuos ganaderos. La actividad ganadera genera los comnmente denominados
purines, compuestos por deyecciones lquidas, excrementos, slidos y el agua utilizada
en las granjas. Se puede obtener una valoracin energtica del residuo mediante el
tratamiento biolgico al que debe ser sometido para ser empleado como abono en
campos de cultivo.
VENTAJAS DEL USO ENERGTICO DE LA BIOMASA
La energa renovable con biomasa tiene un enorme potencial de desarrollo en Espaa, segundo
pas de la Unin Europea en superficies forestales. De hecho, es la nica energa renovable que
presenta un balance econmico positivo por los beneficios generados [4].
La produccin de energa elctrica a partir de biomasa tiene un rendimiento neto prcticamente
cero en lo referido al ciclo del carbono (no es nulo debido a la energa invertida en transporte,
pretratamientos,). Los cultivos forestales energticos explotan determinadas especies que por
su veloz ritmo de crecimiento, mayor capacidad de rebrote, alta acumulacin de energa por
unidad de peso y gran adaptabilidad a diversas condiciones edafolgicas, son energticamente
eficientes y econmicamente rentables. Los cultivos energticos estn todava poco
desarrollados en nuestro pas, sin embargo el aprovechamiento de biomasa residual est
cobrando cada vez un papel ms importante. Espaa cuenta con millones de toneladas de
residuos forestales que se pudren en los bosques y que podran ser empleados en la produccin
de energa.
La biomasa lleva implcitas las tres patas de la sostenibilidad.
2. Antecedentes
- La sostenibilidad ambiental en todos aquellos aprovechamientos de biomasa residual forestal

como combustible, incide en:
Cumplimiento de los objetivos ambientales fijados internacionalmente.
La reduccin de riesgos de incendios, pues supone la retirada de masa forestal de los
bosques. El aprovechamiento de los residuos de origen biolgico es adems econmico y
aumenta la puesta en valor del bosque.
Captura de CO2.
- La sostenibilidad econmica es muy importante porque:
Estas nuevas actividades que se incorporan al sistema productivo primario, generan
rentas, sobre todo en el medio rural.
El cumplimiento del Plan de Energas Renovables 2020 (PER 2020) en cuanto a la
ampliacin de la potencia instalada para el aprovechamiento de la biomasa supondr la
creacin de 15.000 nuevos empleos directos e indirectos en zonas rurales, a travs de
una inversin de 3.000 millones de euros que generarn 671 millones de euros anuales
de nuevos ingresos y ahorros.
- La sostenibilidad social se fundamenta en el desarrollo rural de este tipo de energa, que fija la
poblacin en ncleos de tierras agrcolamente no rentables, que no compiten con usos
agroalimentarios. Adems, la biomasa podra participar en la sustitucin progresiva del carbn
nacional creando empleo en las cuencas mineras afectadas, dada la proximidad de dichas
cuencas con zonas forestales con potencial de produccin de biomasa.
2.4
INCONVENIENTES DEL USO ENERGTICO DE LA BIOMASA
El gas natural o el gasleo de origen fsil son sustancias fciles de transportar, provocan en su
combustin pocos problemas de cenizas e inquemados y son muy homogneos en su
composicin. Ello significa que la automatizacin de la produccin de calor es fcilmente
conseguible, la cadena comercial es fiable y los costes asequibles. Sin embargo, desde el punto de
vista tecnolgico, usar los biocombustibles slidos es lo ms difcil de llevar a cabo [3]. La
estacionalidad de los restos biomsicos, su variabilidad, el hecho de que sean slidos, y sus
caractersticas fisicoqumicas (bajo poder calorfico, alto contenido en agua y cenizas,
heterogeneidad, produccin dispersa en el territorio y en pequeas cantidades) hacen que la
competencia comercial real por el mercado calefactor sea un reto casi insuperable.
Adems, para poder emplear el gas obtenido de forma eficiente, son necesarias una serie de
operaciones posteriores a su salida del gasificador para acondicionarlo (enfriamiento,
limpieza) y es que uno de los principales problemas a resolver en el proceso de gasificacin es
la formacin de alquitranes.
Por ltimo, cabe resear que las ventajas ecolgicas de la biomasa no favorecen la economa de
escala y no son suficientes para su penetracin en el mercado. Pero su conversin en electricidad
o su gasificacin pueden ser buenas rutas de desarrollo. Sin embargo, debe garantizarse la
recuperacin de la biomasa al mismo ritmo que se consume, puesto que si se hace un mal uso de
la misma, pierde su carcter renovable y puede llegar a ser devastador para el medio ambiente
(deforestacin, prdida de hbitats, desertizacin).
2. Antecedentes
2.5
PROCESO DE GASIFICACIN DE BIOMASA
Entre todos los procesos de conversin de biomasa, la gasificacin es uno de los ms

prometedores, ya que la eficiencia energtica es mayor que en el proceso de combustin. Se
trata de una tecnologa muy compleja que tiene como principal ventaja un potencial de alto
rendimiento, pero por el contrario se trata de plantas con altos costes de inversin y que
requieren un aprovisionamiento de biomasa muy homogneo.
La gasificacin es un proceso termoqumico en el que la materia carbono es transformada en un
gas combustible (syngas) con cierto poder calorfico mediante una serie de reacciones que
ocurren a elevadas temperaturas en presencia de un agente gasificante (existen una gran
variedad, como son vapor de agua, oxgeno o aire).
Fig 2. Fases existentes en un proceso de gasificacin [3].

La operacin de un gasificador de biomasa depende de las complejas reacciones qumicas que se
llevan a cabo en su interior. Una partcula de biomasa que entra a un gasificador sufre
inicialmente un proceso de secado seguido de una descomposicin trmica (pirlisis o
desvolatilizacin). Posteriormente tienen lugar diversas reacciones slido-gas o gas-gas entre
todos los componentes presentes. Finalmente, se obtiene el syngas, que se compone de CO2, CO,
H2, H2O y otros hidrocarburos junto con otros subproductos que es necesario eliminar antes de
su uso final (slidos parcialmente quemados o char, cenizas y compuestos condensables o tars).
Puesto que la etapa de pirlisis es ms rpida que la de gasificacin, es esta ltima la que
controla el proceso global.
En comparacin con la combustin, la gasificacin es ms eficiente en la produccin de energa
elctrica y trmica [5], por lo que, considerndose neutra al balance de CO2, es la nica que
puede remplazar directamente a los combustibles fsiles. Adems se puede considerar su
instalacin tanto en aplicaciones avanzadas de pases desarrollados como para la electrificacin
rural de instalaciones aisladas o pases en vas de desarrollo.
3. EQUIPOS DE LA PLANTA
3. Equipos de la planta
3. EQUIPOS DE LA PLANTA
A continuacin se presentan los equipos que se van a considerar en este trabajo de la planta de
produccin de electricidad a partir de la combustin en un motor de combustin interna del gas
de sntesis generado a partir de la gasificacin de 600
3.1
de biomasa.
REACTOR DE GASIFICACIN
Los gasificadores se utilizan para transformar la biomasa en un gas combustible. El movimiento

relativo entre el agente gasificante y el slido marca el tipo de reactor del que se trata:
Lecho mvil: el slido se mueve lentamente.

o
Gasificador downdraft o de corrientes paralelas: el aire entra por la parte

superior y el gas se extrae por la inferior. Este es la configuracin de gasificador
ms empleada en la obtencin de electricidad y la que se toma como referencia
en este trabajo.
Gasificador updraft o de contracorriente: el agente gasificante se introduce por la
parte inferior y atraviesa el reactor de lecho mvil mientras que la biomasa se
introduce por arriba y desciende alcanzando cada vez temperaturas ms altas.
Este tipo de gasificador es sencillo y su eficiencia es alta, adems de que los gases
no tienen tiempo de craquearse por lo que se transforman en fracciones ligeras
no condensables, lo que tambin tiene su punto negativo, la alta tendencia a
formar alquitranes.
Fig 3. Gasificador downdraft y updraft [6]
Lecho fluidizado: las partculas slidas se encuentran en suspensin gracias a una

corriente ascendente de fluido, de gas. Idealmente, cuando una partcula entra en este
tipo de reactores, se supone que sufre instantneamente todas las fases del proceso de
gasificacin, por lo que se considera que la temperatura del reactor es constante y que la
transformacin se produce de forma homognea dentro del mismo. Dependiendo de su
diseo, existen dos tipos:
o Burbujeante.
Circulante: la fluidizacin es rpida, sin burbujas y se recircula toda la material

slida que no es expulsada a travs del cicln.
Fig 4. Gasificadores de lecho fluidizado: (a) Burbujeante; (b) Circulante [7].
3.2
SEPARADOR CICLNICO
Para la limpieza de la corriente de gas contaminada por partculas, la mejor opcin entre los
equipos de coleccin es un separador ciclnico [8]. Los ciclones usan la fuerza centrfuga
generada por un torbellino para separar las partculas del gas con una baja prdida de carga, un
diseo sencillo y barato y una alta eficiencia para partculas con un dimetro mayor a 10
,
como es el caso de este trabajo (en la bibliografa revisada se habla de tamaos de partculas de
250 350
[9] y hasta incluso de 500
[10]).
Es un equipo que no presenta mayores complicaciones en su diseo, que se detalla en el
siguiente captulo.
3.3
INTERCAMBIADOR DE CALOR
Tras la salida del gasificador y el paso del gas de sntesis por el separador ciclnico, se tiene una
corriente gaseosa caliente que se considera que mantiene la temperatura a su salida del reactor.
Por ello, tanto para la produccin de electricidad como para el precalentamiento del aire de
entrada al secadero, es necesario enfriar dicho flujo.
Se requiere de un intercambio de calor lquido-gas, por lo se rehye la instalacin de un
intercambiador de carcasa y tubos, los ms abundantes en el mercado, pues su uso se centra
principalmente en aquellos casos en que ambos fluidos son lquidos. La mejor opcin en este
caso se trata de un intercambiador compacto recuperador de calor de tubos con aletas
exteriores anulares. El gas de sntesis, el fluido caliente, circula por el exterior de los tubos
mientras que el fluido fro, en este caso agua de refrigeracin, circula por el interior de los tubos.
En estos equipos el orden de magnitud del coeficiente de conveccin en el lado fro puede estar
en torno a 10
mientras que en el lado caliente en torno a 100
[11]. Es esta
descompensacin la que hace conveniente el uso de las aletas en el lado del gas. Si el gas es
suficientemente limpio, se puede utilizar una densidad de aletas elevada (200-240 aletas por
metro), pero si la diferencia entre ambos coeficientes de conveccin no es tan acusada, conviene
rebajar esta relacin puesto que ya no resulta efectiva, hasta valores de 75-120 aletas por metro.
4. DESARROLLO DE LOS MODELOS

MATEMTICOS
4. Desarrollo de los modelos matemticos
4. DESARROLLO DE LOS MODELOS MATEMTICOS

El objeto de este proyecto no es estudiar con detalle el proceso que tiene lugar en el gasificador,
sino generar modelos que puedan predecir la composicin de salida y cmo vara esta ante
cambios en las condiciones de entrada.
Teniendo esto presente, la mejor opcin son los modelos 0-D. Son modelos tericos en los que se
incluyen algunas correlaciones experimentales, que permiten establecer una relacin entre las
variables de entrada y salida del volumen de control sin considerar el proceso que tiene lugar en
el interior del mismo. [12]
Este modelo considera que se alcanza el equilibrio termodinmico, lo que es una hiptesis
aceptable teniendo en cuenta que las reacciones principales son rpidas e irreversibles. Los
modelos de equilibrio predicen el mximo rendimiento de un sistema a partir de los productos
de entrada, mientras que los cinticos se centran en el progreso de la composicin en las
diferentes zonas del reactor. Los modelos de equilibrio son computacionalmente ms sencillos y
sirven para este trabajo pues son una buena aproximacin para simular gasificadores
downdraft, sin embargo se debe tener en cuenta que el modelo, por lo general, sobrestima los
rendimientos de H2 y CO, subestima los de CO2 y predice una salida de gases libre de CH4, tars y
char.
Este trabajo presenta un modelo terico matemtico para gasificacin de biomasa
implementado en el programa Engineering Equation Solver (EES) [13]. Se ha escogido este
programa porque contiene todas las funciones termodinmicas necesarias. El modelo est
basado en los modelos de equilibrio termodinmico aunque incluye algunas modificaciones
obtenidas a partir de resultados experimentales para adaptarlo a un proceso real, en el cual solo
se alcanza parcialmente el equilibrio.
Es de resear que se han introducido conceptos propios de programacin al escribir el cdigo en
EES puesto que ha sido necesario definir algunas funciones para calcular distintas propiedades
de las materias presentes en el proceso, procedimientos para simplificar clculos que se repetan
a lo largo del programa o funciones que se han denominado espera y que sirven para separar
el cdigo en distintos bloques de forma que la compilacin resulte ms rpida y eficiente. Tanto
la potencia de clculo como la imposibilidad de llamar a funciones dentro de funciones,
caractersticas de la versin del EES facilitada por la Universidad de Zaragoza con la que he
trabajado, han requerido de esta solucin. Adems, en esta fase de diseo cada equipo se ha
desarrollado en un archivo distinto, por lo que se han generado las funciones $Export y $Import
necesarias para exportar e importar respectivamente los datos de origen con los que se trabaja,
o bien las incgnitas resueltas de un equipo y que suponen los datos de partida del equipo
siguiente. Estas lneas de cdigo se muestran como comentario en el trabajo puesto que de lo
contrario no apareceran en la impresin.
A continuacin, se exponen los distintos cdigos de diseo de equipos realizados en EES,

comenzando por el equipo principal, el reactor de gasificacin, y continuando en el orden en el
que se encuentran en la planta de produccin de electricidad, el separador ciclnico y el
intercambiador de calor. Los cdigos de los mismos se encuentran en el Anexo 4.
4.1
REACTOR DE GASIFICACIN
El reactor de gasificacin se plantea tal y como se muestra en la Figura 5, de forma que se tienen
dos corrientes de entrada (F1 es la biomasa de origen cuyas composiciones se corresponden con
las obtenidas en el conversor explicado a continuacin, y F2 es el caudal de agente gasificante, en
este caso, aire) y otras dos de salida (F3 es la salida de slidos mientras que F4 es la corriente
que nos interesa puesto que se trata del gas de sntesis o syngas).
Fig 5. Esquema del reactor de gasificacin.

El primer bloque del cdigo en EES incluye las ecuaciones necesarias para expresar la biomasa
de origen sobre una base comn. Las composiciones de las maderas y residuos forestales con las
que se trabaja se obtienen de la base de datos de ECN [14]. Los datos de la composicin de la
biomasa utilizada, introducidos en el modelo de EES son los siguientes:
Tabla 1. Composicin de la biomasa de la simulacin.
Pine sawdust (#1267)

%
6,00
0,10
! " (wt%)
# ! $ (wt%)
83,18
%&! $ (wt%)
16,82
&! $ (wt%)
51,05
'! $ (wt%)
6,01
(! $ (wt%)
0,10
)! $ (wt%)
0,01
*! $ (wt%)
42,83
+'#! $ (MJ/kg)
19,02
''#! $ (MJ/kg)
20,33
En la base de datos de ECN se expresan las composiciones en masa (% en peso, expresado en la

bibliografa y en lo que sigue como % wt) de la biomasa seleccionada en varias bases:
10
As Received (ar), que expresa a la biomasa en su totalidad tal y como se encuentra en la

naturaleza contando su humedad y cenizas.
Dry, las composiciones no tienen en cuenta la humedad del material.
Dry Ash Free (daf), las composiciones no tienen en cuenta ni la humedad ni las cenizas
del material.
Fig 6. Base de datos para biomasa y residuos forestales de ECN.

En la base de datos utilizada, la mayor parte de las biomasas tienen datos de composicin en las
diferentes bases. Sin embargo ha sido necesario incluir este mdulo de conversin ya que en
algunos casos, como en el ejemplo de la tabla 1, los diferentes valores de composicin son
expresados en diferentes bases.
11
Tabla 2. Bases empleadas en la expresin de la composicin de madera.
Anlisis inmediato (% wt)

Humedad
ar
Cenizas
dry
Materia voltil
daf
Carbono fijo
daf
Anlisis elemental (% wt)
C
daf
H
daf
N
daf
S
daf
O
daf
Cenizas
dry
Tras ello, se sita un procedimiento de espera (Procedure CONVERSOR. Ver Anexo 4), de forma
que la compilacin del cdigo se realice en bloques porque de lo contrario el programa comete
errores debido a la acotacin incorrecta que realiza de las variables incgnitas. As, se pueden
comenzar los clculos siguientes con la composicin de la biomasa de origen como dato y no
como incgnita.
En segundo lugar se realiza el ajuste estequiomtrico de la reaccin de combustin de la

biomasa,. Hay que tener en cuenta que durante la combustin de una especie orgnica, el S
produce SO2 y el N NOx (mezcla de NO2 y NO), que se puede considerar como NO2 en su
totalidad, para estos clculos. De esta forma se conoce el aire estequiomtrico necesario para la
combustin completa, dato necesario, posteriormente, para calcular el defecto de aire con el que
se realiza la gasificacin. Las ecuaciones planteadas para conocer la relacin estequiomtrica,
consistentes en aplicar los balances de materia al proceso de combustin, son las siguientes:
Relacin entre O2 y N2 en la entrada de aire.

(
Balance al carbono.
Balance al hidrgeno.
Balance al oxgeno.
*-4 + 2*
Balance al azufre.
,- ./01
79
*
21
&-4 = &*
,- ./01
'-4 = 2' *-
,- ./01
= 2&*
,- ./01
)-4 = )*
,- ./01
./01
+ ' *-
./01
+ 2)*
,- ./01
,- ./01
A continuacin, se procede a expresar las ecuaciones linealmente independientes necesarias

para resolver las 9 incgnitas restantes en el sistema.
12
Se plantean los siguientes balances de materia:
Relacin entre O2 y N2 en la entrada de aire.
79
*
21 ,Prdida de materia slida desglosada en la del carbono fijo y cenizas. Se supone una
prdida del 5% del carbono fijo con respecto al existente en el caudal de entrada
(aunque este dato se puede modificar puesto que no es una incgnita) y que se extraen
todas las cenizas.
Balance al carbono.
&-4 = &*- + &* ,- + &' ,- + &-6
Balance al hidrgeno.
'-4 = 4&' ,- + 2' ,- + 2' *- + 3('6,- + 2)' ,Dos balances al nitrgeno. Hay que tener en cuenta que durante la gasificacin de una
especie orgnica el N de la biomasa va a NH3 mientras que el N2 del aire acta como una
sustancia inerte, es decir, no reacciona.
(-4 = ('6,2( ,- = 2( ,Balance al azufre. Hay que tener en cuenta que durante la gasificacin de una especie
orgnica el S va a SH2.
)-4 = )' ,Balance al oxgeno.
*-4 + 2* ,- = &*- + 2&* ,- + ' *(
,-
Por su parte, se plantea el balance de energa al reactor, de forma que el calor del mismo es
proporcional a la entalpa de formacin aparente de la biomasa, al aplicar un coeficiente
denominado El clculo de las entalpas de cada compuesto est detallado en el Anexo 4 mientras
que se consideran unas prdidas de calor del 10% [15].. La ecuacin del balance de energa
resulta:
. 809 : /;$;8 8;=
>
? 9!A809/
'? %?
9
'; = $,;
+ ;,DE ;,
>
. 80;B9/
' %
.$
La entalpa de los slidos, compuestos principalmente por carbono, que salen del gasificador se
calcula a partir de una relacin experimental [5] basada en la entalpa del grafito:
8F
I = 0,0004J + 0,8679
381,61
J
Por ltimo, se tienen dos ecuaciones ms para completar el sistema. Se trata de las ecuaciones de
equilibrio, de dos de las reacciones mas representativas del proceso de gasificacin son:
Reaccin de metanacin.
13
& + 2' &' ,
Reaccin water-gas shift.
&* + ' * &* + ' ,
' = 74.900 H
LM
N' = 41.200 H
LM
La constante de equilibrio se puede relacionar con el cambio en la energa libre de Gibbs

mediante las siguientes ecuaciones:
;
NPJL
O QJ
S
S
S
S
PR
J
'R
J
'R
J 'R
0
= %S J +
U
J
J
J
A partir de las ecuaciones anteriores y la tabla 4 (ver Anexo 1) se pueden calcular las constantes
de equilibrio de las dos reacciones para una temperatura determinada.
Los parmetros obtenidos del diseo han sido los siguientes:
En el reactor gasificador, para un flujo msico de entrada de 600
de dicha biomasa y una
temperatura de salida del reactor de 1.073 , con un porcentaje de aire respecto del
estequiomtrico para la combustin de WQ = 25,92% mientras que el ratio biomasa/aire es de
%X = 2,569.
El poder calorfico del gas de sntesis es:
+'#- = % ,; '; = 3,678 10Z H = 1.021,67
Por lo tanto, la eficiencia del gasificador es de [

15,06 106 H
600
= 9,036 10
. 809
\]^E_
\]^`a
6,Zbc4Sd
= e,S6Z4Sd = 40,7%. ('+# =
H . La eficiencia total de la planta, definida como
[090 f = [R909 ,gh [:.=. !9 [ . 809 , teniendo en cuenta que la eficiencia tipo de un motor de
combustin interna es [R909 ,gh = 30% mientras que la de un generador elctrico se puede
considerar [:.=. !9 = 95%, se tiene [090 f = 11,6%.
Teniendo todo lo anterior en cuenta, la potencia elctrica que genera la planta en este caso es
i. = [R909 ,gh [:.=. !9 '+#- = 291,175
.
4.1.1
DISEO MECNICO DEL REACTOR DE GASIFICACIN
Se ha implementado el diseo mecnico del reactor de gasificacin en el software Abaqus que

permite realizar clculos estructurales estticos mediante la aplicacin de elementos finitos.
14
El diseo se presenta como un sistema axisimtrico, es por ello que se simula 1 8 del reactor
cuyo proceso se detalla en el Anexo 2. Se ha diseado, a partir de la bibliografa consultada sobre
plantas de gasificacin experimentales situadas en laboratorios ([38, [16], [17], [18]) un reactor
gasificador downdraft de las siguientes dimensiones:
1.400
de altura total.
8 jm de dimetro en la zona del lecho, con una altura de 420
.
12 cm de dimetro en el resto del reactor, con una altura de 980 mm.
1,5
de espesor de la pared y 50
de espesor de las tapas.
El material empleado es acero inoxidable, mientras que se ha simulado la presin interior

imponiendo una carga uniforme en sentido positivo hacia el exterior de valor 1,05
, es decir,
algo ms alta que la presin atmosfrica. Se ha realizado adems un mallado cuadrtico
estructurado A 4-node doubly curved general-purpose shell, finite membrane strains, idneo
para definir las paredes del cilindro y la membrana o shell, puesto que lo se simula en este
modelo es la superficie media, es decir, el plano que est a la misma distancia del dimetro
interior que del exterior.
A continuacin se exponen los resultados obtenidos de la simulacin. En cuanto a la distribucin
de la presin a la que se ven finalmente sometidas las paredes, se observa en la figura 7 que
estas son mayores en la zona del freeboard donde tiene lugar la oxidacin, se estrecha la seccin
y se introduce el gas gasificante. Es por ello que esta zona deber estar reforzada o tener un
espesor mayor al del resto de la pared y similar al de las tapas. Por su parte, en la figura 8 se
observa que las tensiones son prcticamente constantes a lo largo del gasificador, sin embargo
estas se disparan en la zona cercana a las tapas (Figura 9). Se produce una excesiva distorsin
del mallado en las zonas cercanas a las tapas debido a las altas tensiones que all aparecen y que
no reflejan la realidad y anulan el modelo en esa zona. La solucin, que se escapa del objetivo de
este trabajo, sera realizar un afinado de la malla en esas zonas, con distintos tipos de elementos
de malla.
15
Fig 7. Distribucin de presiones.
Fig 8. Distribucin de tensiones
Fig 9. Detalle de la distribucin de tensiones en la zona cercana a la tapa.

Por ltimo, los desplazamientos (Figura 10) tambin muestran un comportamiento constante a
lo largo del gasificador, sin embargo se puede decir que se alcanzan los mximos
desplazamientos en las zonas de inflexin de la membrana.
16
Fig 10. Distribucin de desplazamientos.
4.2
SEPARADOR CICLNICO
Por lo general, las dimensiones de un separador de cicln se basan en:
Nmero de giros efectivos (. . Se demuestra que es igual a 8, para un separador tipo de

alta eficiencia, a partir de las dimensiones que se enumeran posteriormente.
1
'4 5 '
2
2m8 5 2m8
(. = k'4 5
l=
k2m8 5
l=8
'
2
m8
2
Densidad de las partculas de entrada al cicln, n? .
Se considera que es n? = 1.050
6,
tomada como valor medio de la bibliografa
revisada [7].
Caudal de entrada al cicln, 7- .
Tomados como relaciones a estas dimensiones, se tienen el resto que definen el diseo del
separador [19]:
;
m9
4 ; m. =
m9
2;' =
m9
2 ; '4 = 2m9 ; ) = 0,625m9 ; ' 2 2m9 ; m! 2
17
m9
Fig 11. Esquema de un separador de cicln [NEVERS].

Por ltimo, se definen parmetros como el dimetro de corte m?pq , el rea efectiva r8;8f= , la
velocidad del gas s en la entrada del cicln y la eficiencia [; del equipo:
r8;8f= 2
m?pq 2 t
; '; s
9 ;
2u(. s n/
; [; 2
7r8;8f=
m?x
ym z
?pq
15w
m?x
ym z
?pq
El clculo de las propiedades de la mezcla gaseosa de entrada al cicln requiere conocer la

composicin y los valores de dos propiedades intensivas independientes tales como la
temperatura y la presin [20]. Las propiedades extensivas de la mezcla, tales como la energa
interna U, la entalpa H o la entropa S, se pueden obtener sumando la contribucin de cada
componente en la condicin en la que ste existe en la mezcla. Por ejemplo, la entalpa de los
gases ideales es funcin exclusiva de la temperatura y se tiene:
18
' = '4 + ' + + '| = > ';

;}4
Sin embargo, para conocer otras propiedades de la mezcla de gases, como la viscosidad o la
conductividad, se requiere de otros clculos que se aproximen mejor al modelo. A finales de los
aos 40, el ingeniero qumico C.R Wilke describi junto con su compaero Pin Chang una
correlacin experimental que determinaba los coeficientes de difusin en mezclas de gases [21].
Se ha implementado la misma en un procedimiento de EES (Ver Anexo 4, ya en el cdigo del
separador ciclnico) para calcular con mayor exactitud la viscosidad de una mezcla de gases de
siete molculas distintas, en concreto la del gas de entrada al cicln (se ha eliminado el SH2 de
estos clculos pues generaba problemas de convergencia por sus valores de viscosidad y
conductividad tan cercanos a cero y representaba un porcentaje muy pequeo de la composicin
del gas de sntesis).
Con los datos expuestos anteriormente, se tienen las siguientes caractersticas del separador
ciclnico:
Tabla 3. Dimensiones principales del separador ciclnico.
m9
;
m.
'
'4
)
'
m!
2
0,5
1
1
1
1,25
4
0,5
El dimetro de corte resulta ser de m?~S = 36,36

mientras que se logran separar partculas
de hasta 109,1
con una eficiencia del 90%. En la figura X se puede observar la curva de
eficiencia del separador ciclnico diseado.
Fig 12. Eficiencia del separador ciclnico obtenida en el diseo en EES.

19
4.3
INTERCAMBIADOR DE CALOR
Para realizar un primer diseo, se emplea el mtodo MLDT en el que se conocen todas las
temperaturas. En este caso, se conoce la temperatura de entrada caliente, que se supone que
coincide con la de salida del reactor y la temperatura de entrada fra, que es agua a 20C. Se debe
fijar por otra parte la temperatura de salida del gas, que marcar el volumen de intercambio de
calor y el flujo de agua necesarios para alcanzar las condiciones deseadas. Por ltimo, a partir de
los balances de energa aplicados a cada fluido, se obtiene la temperatura de salida de la
corriente fra.
= % j8 J8,.
J8,/ 2
:A
j$ J$,/
J$,.
Tras ello, se estima un rea superficial exterior tal que:

r/A?,.0,./0 2
./0 Rf
Se estima un valor de coeficiente de transferencia de calor ./0 en base a las diferentes

sustancias involucradas en el intercambio de calor [11]. En este caso, al tratarse el fluido caliente
de un gas a baja presin ~1 X y el fro de agua de refrigeracin, se tiene un ./0 2
105
. Por su parte, el salto medio logartmico de temperaturas se calcula como el de flujo
a contracorriente corregido por un factor % que depende de las temperaturas de entrada y salida
de ambos fluidos, pero que en este caso se considera igual a la unidad porque la circulacin de
ambos fluidos es globalmente a contracorriente, hecho que se corrobora obteniendo el perfil de
temperaturas que se genera en el intercambiador diseado.
Fig 13. Perfil de temperaturas en el intercambiador de calor. Flujo a contracorriente.
Tras esto, se calcula el rea superficial exterior por unidad de longitud de tubo r/A?,.0 , pues se
relaciona con la anterior proporcionalmente al nmero de tubos por el ancho del
intercambiador de calor:
(0 +0 2
r/A?,.0,./0
r/A?,.0
Con este dato se calculan las dimensiones previas del intercambiador como son el nmero de
tubos (0 , nmero de circuitos (8 , nmero de tubos por fila (0,$;f y nmero de filas de tubos (] .
Hay que aadir, que desde el punto de vista termo-fluidodinmico es preferible una disposicin
de los tubos alternada, en tringulos, porque se optimiza el coeficiente de conveccin, de forma
que )D 2 ) .
20
Para el reclculo y diseo definitivo es necesario hallar los dos coeficientes de conveccin. El del
lado interior de los tubos se corresponde con un caso de conveccin forzada interior. Se
desarrolla un procedimiento (se sita en la parte superior del cdigo del intercambiador de
calor. Ver Anexo 4) para tener en cuenta todos los casos tal y como sigue [22]. En el caso de que
el flujo se encuentre en rgimen laminar y teniendo en cuenta que el flujo de calor en la pared
interior del tubo se puede considerar constante ? = j O al ser . apreciablemente menor que
; [11], se ha de determinar si la entrada es puramente trmica, en cuyo caso se aplica la
correlacin de Hausen; o si es combinada trmica-hidrulica, en cuyo caso se aplica la
correlacin de Kays. Si por el contrario se est ante un caso de rgimen turbulento, la
correlacin a aplicar ser la de Gnielinski. Por su parte, el clculo del coeficiente de conveccin
en el lado exterior de los tubos se corresponde con un caso de conveccin forzada exterior y la
correlacin a aplicar, con una disposicin del banco de tubos en tringulos, es la de Briggs y
Young. Las correlaciones empleadas son, respectivamente:
4
' s O: (s$ = 4,36 + 0,023 ln G1 + 0,012P$ I , P$ =

4
: (s$ = 4,36 + 0,036 ln G1 + 0,011P$ I

POM
X
QO$ iX$ m;
+0
m
8 QO$ 1000iX$
1 + k ; l
: (s$ =
S,~
+0
6 1I
1
+
12,7
k
l
GiX
$
Ls : (s8 = 0,134QO8 S,Zc4 iX8
6t
S,
S,44
Finalmente, se calcula el coeficiente global de transferencia de calor, mediante la ecuacin

general aplicada al lado caliente. Se desprecia la resistencia de enfriamiento del fluido fro (lado
interior de los tubos) pero no as la del fluido caliente (gas con restos de partculas y cenizas),
por lo que queda:
"
1
1
Q.,8
r8
=
+
+ r8 Q? +
8 [9,8 8 [9,8
[9,$ $ r$
La eficiencia de la superficie en el lado fro es la unidad puesto que no hay aletas, mientras que la
del lado caliente se basa en el concepto de aleta de punta adiabtica, anular rectangular [22]
cuyo clculo en EES se simplifica al tener incorporadas las funciones de Bessel. Por su parte, la
resistencia de enfriamiento del fluido caliente, considerado como gas de combustin, se puede
estimar como 4 10
[22] mientras que la de conduccin de la pared que separa ambos
fluidos, considerando la conductividad del material del tubo, resulta
Q? =
G I
x
.
]
Se deben realizar varias iteraciones, modificando parmetros del diseo del intercambiador de
calor para obtener un coeficiente global de transferencia de calor cercano al estimado. Por
21
ejemplo, se comprueba que la diferencia entre ambos coeficientes de conveccin no es tan

grande como la comentada anteriormente, por lo que la estimacin de 275 Aletas por unidad de
longitud del tubo se ha rebajado a 100. Aumentar la superficie de intercambio no implica
necesariamente una mayor transferencia de calor porque el coeficiente de conveccin exterior
puede disminuir con un aumento de la superficie o de la densidad de aletas. Se ha observado
esto tras varias iteraciones realizadas sobre las variables Dimetro exterior de las aletas m ,
Dimetro exterior de los tubos m. y Nmero de aletas por unidad de longitud del tubo ( f.0 /,]
en una tabla paramtrica en EES, de la cual se presenta en la siguiente figura. Del diseo del
intercambiador de calor se obtiene un UA=120,69 WK.
Fig 14. Iteraciones en EES para reclculo y diseo definitivo del intercambiador de calor.
22
5. RESULTADOS
5. Resultados
5. RESULTADOS
Tras el diseo realizado anteriormente, ahora se procede a simular la planta de produccin de
electricidad tomando como parmetros las caractersticas de los equipos diseados. A
continuacin se van a validar los modelos del reactor de gasificacin implementados tanto en
EES como en Hysys, con la bibliografa revisada.
Una primera comparacin entre el modelo implementado en EES y la simulacin llevada a cabo
en Hysys, simulados bajo las condiciones expuestas en el apartado anterior, valida los modelos
implementados pues las composiciones del gas de sntesis, a excepcin del metano, son muy
parecidas. Estos son los resultados obtenidos:
Tabla 4. Composiciones del gas de sntesis en ambas simulaciones.
EES
HYSYS
21,470
22,610
&* LM %
4,350
4,973
&* LM %
5,892
12,460
&' LM %
13,300
10,020
' LM %
6,767
4,607
' * LM %
0,003
0,000
)' LM %
41,140
45,330
( LM %
0,068
0,000
('6 LM %
Los resultados obtenidos son parecidos, excepto la concentracin de metano, que es la que ms
diferencias presenta entre los dos resultados. Esto se puede deber principalmente a que las
reacciones elegidas para el modelo del EES son diferentes de las que emplea Hysys.
Por otra parte, se decidi hacer una validacin del modelo implementado con datos
experimentales encontrados en la bibliografa sobre gasificacin de biomasa. La primera
referencia consultada corresponde con los trabajos de Magn Lapuerta y cols. [23] sobre
gasificacin en un reactor de lecho fluidizado circulante de distintos residuos forestales y
agrcolas (poda de pinos, olivos y vides, serrn) con aire como agente gasificante. Las
condiciones de operacin del proceso son a presin atmosfrica mientras que la temperatura del
reactor se mantiene constante a aproximadamente 1.523 Aunque los autores trabajan con un
reactor de tipo diferente, esto afectara a las condiciones de alimentacin de la biomasa, o a la
produccin de alquitranes o char, pero no al equilibrio. Si, los resultados experimentales se
acercasen a este, los resultados tendran que ser comparables.
Se analizan dos de las biomasas expuestas en la investigacin, Pinus pruning wastes y Olive
pruning wastes ya que se tienen dos biomasas con unas caractersticas muy similares en la base
de datos ECN, como son Pine sawdust (#1267) y Olive branches (#3347) respectivamente.
Las figuras 15 y 16 muestran la composicin en volumen del gas de sntesis para distintos
-fA|9 R9f ;9R / .=0 ! valores del ratio biomasa/aire o Frg definido como %X = -fA|9 R9f ; . .=0 ! ; . En la
parte superior de la imagen se muestran los resultados obtenidos en el modelo de EES, mientras
que en la parte inferior se sitan los resultados de la investigacin.
23
5. Resultados
Fig 15. Concentraciones de CO y H2 para las diferentes biomasas y valores de Frg.
Fig 16. Concentraciones de CH4 y CO2 para las diferentes biomasas y valores de Frg.
24
5. Resultados
Se observa que los resultados son buenos, ya que todas las tendencias son similares excepto en
el caso del CO2, donde las trayectorias son opuestas.
La segunda referencia consultada corresponde con los trabajos de Jun Xie y cols. [15] sobre la
gasificacin en un reactor de laboratorio de lecho fluidizado de Pine sawdust o Serrn de pino con
aire como agente gasificante. Las condiciones de operacin del proceso son a presin
atmosfrica mientras que la temperatura del reactor se mantiene constante a aproximadamente
1.073 K. Por su parte, en la investigacin se trabaja con un ER entre 0,21 0,25 mientras que en
el modelo el rango aumentaba desde 0,22 0,35. Adems, el flujo de entrada de biomasa se
mantiene en el rango entre 0,445
0,552
haban realizado las simulaciones anteriores.
, a diferencia de los 600
con los que se
En el modelo de EES se vuelve a analizar el Pine sawdust (#1267) ya que tiene unas
caractersticas muy similares a la biomasa expuesta en la investigacin. Las figuras 17 y 18
muestran la composicin en volumen del gas de sntesis para distintas temperaturas del reactor.
Las lneas continuas con smbolos cerrados corresponden a los datos experimentales mientras
que las lneas discontinuas y smbolos abiertos lo hacen a las predicciones [15].
Fig 17. Concentraciones del gas de sntesis para diferentes temperaturas del reactor.
25
5. Resultados
Fig 18. Concentraciones del gas de sntesis para diferentes valores de Frg.
Se observa que el modelo terico de EES se ajusta perfectamente a las predicciones obtenidas
del modelo desarrollado en la investigacin, sin embargo ste se aleja de sus resultados
experimentales especialmente con respecto al H2 y al CO, aunque es aceptable el
comportamiento que se observa tanto en el CH4 como en el CO2. Un aumento en la temperatura
conduce a un ligero aumento en las concentraciones de hidrgeno y monxido de carbono, as
como a una disminucin de las de dixido de carbono y metano (se observan estas tendencias en
el modelo en EES, sin embargo el de la bibliografa presenta alguna discrepancia). Esto se debe a
un desplazamiento del equilibrio hacia los productos en las reacciones endotrmicas y hacia los
reactivos en las exotrmicas [12].
La gasificacin es un proceso complejo en el que entran en juego numerosos factores por lo que
no hay una nica explicacin clara sobre la dependencia de los compuestos con las variables. Las
diferencias entre modelos tericos y experimentales se deben principalmente a la suposicin de
los modelos de equilibrio de que todas las reacciones transcurren al mismo tiempo, lo cual no es
cierto. A temperaturas altas (en torno 900), las concentraciones tericas se acercan mucho a
las tericas, lo que puede significar que en procesos reales nicamente a altas temperaturas se
llega a alcanzar el equilibrio.
26
6. CONCLUSIONES
6. Conclusiones
6. CONCLUSIONES
Las conclusiones que pueden deducirse del presente trabajo son:
Se ha simulado una planta de produccin de electricidad a partir de la combustin en un

motor de combustin interna del gas de sntesis generado a partir de la gasificacin de
600 de biomasa.
Se ha calculado la eficiencia global del sistema y la potencia generada por la planta de

produccin de electricidad.
El modelo se ajustan bien a los datos reales cuando se alcanzan las condiciones de
equilibrio. A temperaturas bajas, por debajo de 750, el modelo no es capaz de
reproducir de un modo tan exacto el proceso real de gasificacin. Las diferencias
observadas se deben principalmente a que los modelos implementados en este trabajo
no pueden predecir la produccin de alquitranes, mientras que estos se generan en las
simulaciones consultadas en la bibliografa.
En procesos reales, no suele ser habitual que el metano alcance su concentracin de
equilibrio aun cuando el resto de componentes del gas de sntesis s que llegan a valores
muy similares. Los resultados obtenidos con respecto al metano tienen una gran
variabilidad y, por tanto, una gran incertidumbre.
La temperatura de gasificacin tiene un importante efecto sobre la composicin y calidad

del gas de sntesis obtenido. Un aumento en la temperatura conduce a un aumento en las
concentraciones de hidrgeno y monxido de carbono, as como a una disminucin de
las de dixido de carbono y metano.
Se ha realizado el diseo mecnico del reactor de gasificacin y se ha comprobado que la

distribucin de presiones indica que estas son mayores en la zona de oxidacin o
freeboard, por lo que la solicitacin mecnica de esta parte singular del reactor obliga a
reforzar con un mayor espesor que en el resto de la pared.
Se han diseado otros equipos presentes en la planta de produccin de electricidad,

como son el separador ciclnico y el intercambiador de calor. Se ha obtenido el dimetro
de corte del separador ciclnico para el caso expuesto en el diseo y su eficiencia de
recoleccin, mientras que se han obtenido los principales parmetros de diseo y las
caractersticas trmicas del intercambiador de calor.
27
7. REFERENCIAS
7. Referencias
7. REFERENCIAS
[1] Key World Energy Statistics. [Internet] International Energy Agency (IEA) [consultado 12
junio 2014]. Disponible en:
http://www.iea.org/publications/freepublications/publication/kwes.pdf
[2] Plan energtico de Aragn 2.013-2.020. [Internet] Boletn Oficial Provincial de Aragn, 22
de mayo de 2014, nm. 98, pp. 16.208-16.608 [consultado 27 mayo 2014]. Disponible
en:http://www.boa.aragon.es/cgibin/EBOA/BRSCGI?CMD=VEROBJ&MLKOB=793166000404
[3] SEBASTIN NOGUS, FERNANDO. 2010a. Energas De La Biomasa (Volumen I). 1 edn.
Anon. Zaragoza: Prensas Universitarias de Zaragoza.
[4] Informe de sostenibilidad 2012. [Internet] Ence [consultado 16 agosto 2014]. Disponible
en:
http://www.ence.es/pdf/Ence_Informe%20Sostenibilidad%202012_baja-1705483.pdf
[5] PUIG-ARNAVAT, M., J. CARLOS BRUNO, A. CORONAS. ISSN 0887-0624. 2012. Modified
Thermodynamic Equilibrium Model for Biomass Gasification: A Study of the Influence of
Operating Conditions. Energy & Fuels. 26(2), pp 1385-1394.
[6] MORENO GONZLEZ, JAVIER. 2010. Simulacin de gasificacin de biomasa en lecho
fluidizado burbujeante. Proyecto Fin de Carrera. Universidad Carlos III de Madrid.
Departamento de Ingeniera Trmica y de Fluidos. Ingeniera Industrial: Tecnologas
Energticas.
[7] GOMEZ-BAREA, A., LECKNER, B. ISSN 0360-1285. 2010. Modelling of Biomass
Gasification in Fluidized Bed. Progress in Energy and Combustion Science. 36(4), pp 444509.
[8] PROFESORES DEL DEPARTAMENTO DE INGENIERA QUMICA Y TECNOLOGA DEL
MEDIO AMBIENTE. Ingeniera del Medio Ambiente. Escuela de Ingeniera y Arquitectura
de Zaragoza. Zaragoza, 2013. Tema 5, Coleccin de partculas.
[9] AZNAR, M., MANY, J. J., GARCA, G., SNCHEZ, J. L., MURILLO, M. B. ISSN 0887-0624.
2008. Influence of Freeboard Temperature, Fluidization Velocity, and Particle Size on Tar
Production and Composition during the Air Gasification of Sewage Sludge. Energy &
Fuels. 22(4), pp 2840-2850.
[10]
TINAUT, F. V., MELGAR, A., PREZ, J. F., HORRILLO, A. ISSN 0378-3820. 2008.
Effect of Biomass Particle Size and Air Superficial Velocity on the Gasification Process in
a Downdraft Fixed Bed Gasifier. An Experimental and Modelling Study. Fuel Processing
Technology. 89(11), pp 1076-1089.
[11]
MARN, J. M., GUILLN, S. Diseo y clculo de intercambiadores de calor
monofsicos. 2013. Ed. Paraninfo.
[12]
BORNIQUEL MORENO, DANIEL. Desarrollo de modelos matemticos para
gasificacin de biomasa. Proyecto Fin de Carrera de Ingeniera Qumica. Universidad de
Zaragoza. 2014.
[13]
F-Chart Software. Engineering Equation Solver (EES), Academic Professional
Version V 8.649; F-Chart Software: Madison, WI, 2010.
28
7. Referencias
[14]
Database for biomass and waste. [Internet] ECN [consultado 20 febrero 2014].
Disponible en:
https://www.ecn.nl/phyllis2/Browse/Standard/ECN-Phyllis
[15]
XIE, J., et al. ISSN 0960-8524. 2012. Simulation on Gasification of Forestry
Residues in Fluidized Beds by Eulerian-Lagrangian Approach. Bioresource technology.
121 pp 36-46.
[16]
JAIN, V., D. GROULX, P. BASU. ISSN 0306-2619. 2010. Study of Heat Transfer
between an Over-Bed Oil Burner Flame and a Fluidized Bed during Start-Up:
Determination of the Flame to Bed Convection Coefficient. Applied Energy. 87(8), pp
2605-2614.
[17]
PROFESORES DEL DEPARTAMENTO DE INGENIERA MECNICA. Mecnica del
Slido Deformable. Escuela de Ingeniera y Arquitectura de Zaragoza. Zaragoza, 2012.
[18]
Dossier del laboratorio. [Internet] CIRCE [consultado 5 junio 2014]. Disponible
en: http://www.fcirce.es/web/page.aspx?id=labgas
[19]
DE NERVERS, N. Ingeniera de control de la contaminacin del aire. 1998. Ed.
McGraw-Hill Interamericana, Mxico.
[20]
MORAN, M. J., SHAPIRO, H. N. Fundamentos de termodinmica tcnica, vol. 2, p.
695-705. 1999. Revert, S.A.
[21]
WILKE, C. ISSN 0021-9606. 1950. A Viscosity Equation for Gas Mixtures. Journal
of Chemical Physics. 18(4), pp 517-519.
[22]
MARN, J. M., MONN, C., UCHE, J. Transferencia de calor. rea de Mquinas y
Motores Trmicos. Universidad de Zaragoza. Ed. Kronos.
[23]
LAPUERTA, M., et al. ISSN 0378-3820. 2008. Gasification and Co-Gasification of
Biomass Wastes: Effect of the Biomass Origin and the Gasifier Operating Conditions. Fuel
Processing Technology. 89(9), pp 828-837.
29
8. ANEXOS
8. Anexos
8. ANEXOS
8.1
TABLAS PARA EL CLCULO DE LAS CONSTANTES DE EQUILIBRIO
La obtencin de los valores de las constantes de equilibrio de la Reaccin de metanacin y la Reaccin

water-gas shift se lleva a cabo mediante la energa libre de Gibbs.
En la tabla 4 se encuentran los datos de energa libre de Gibbs de formacin, entalpa de formacin y
entalpa de combustin estndar de los distintos compuestos presentes durante las reacciones de
gasificacin. De la misma forma, se muestra en la tabla 4 las funciones de energa libre necesarias. A
partir de estos datos se puede obtener una ecuacin para la constante de reaccin dependiente de la
temperatura.
Tabla 4. Energas libre de Gibbs, entalpas de formacin y funciones de energa libre.
H
9
9
P$9 H
'$9 H
'R
J 'R
0
LM
LM %S
LM
C
0
0
11,6
CO
-137,24
-110,6
204,1
CO_2
-394,64
-393,78
226,4
H_2
0
0
137
H_2O
-228,76
-241,98
196,7
CH_4
-50,85
-74,86
199,4
30
LM
1050
8670
9360
8500
9910
10030
8. Anexos
8.2
ENTORNO DE ABAQUS
Se ha implementado el diseo mecnico del reactor de gasificacin en el software Abaqus

siguiendo el siguiente procedimiento.
Se ha dibujado el perfil del reactor con las siguientes dimensiones:
1.400
de altura.
12 j de dimetro en la zona de freeboard, con una altura de 980
8 j de dimetro en la zona del lecho, con una altura de 420
.
1,5
de espesor de la pared y 50
de espesor de las tapas.
Como el diseo se trata como un problema axisimtrico, no hace falta revolucionar 360 la
figura, sino que sirve con hacerlo 90 tras haber escogido la opcin de membrana deformable en
3D. Tras ello, se generan las secciones Pared y Tapa, cada una con el espesor correspondiente. Se
les aplica a ambas el material creado Acero Inoxidable y se ensambla el conjunto.
Para simular la presin interior uniforme de 1,05
y las condiciones de contorno segn la
figura 19, es necesario crear un paso de simulacin posterior al inicial e indicar que dichas
propiedades se propagan linealmente.
Fig 19. Condiciones de contorno tipo para un diseo de membrana o shell.

Por ltimo, se debe mallar el conjunto para poder compilar el trabajo, ya que Abaqus resuelve
los problemas por elementos finitos. El mallado escogido se puede ver en la Figura X y se trata
de un mallado cuadrtico estructurado A 4-node doubly curved general-purpose shell, finite
membrane strains, idneo para definir las paredes del cilindro y la membrana o shell.
31
8. Anexos
Fig 20. Diseo del reactor de gasificacin en Abaqus. Detalles del perfil revolucionado
90, perfil revolucionado, presin interior y mallado.
32
8. Anexos
8.3
ENTORNO DE HYSYS
Se ha implementado la planta descrita a lo largo del trabajo en Aspen Hysys, una herramienta
informtica que permite disear y modelar distintos procesos y operaciones qumicas. El
objetivo de desarrollar el modelo en Hysys es aportar otro modelo terico, aparte del
desarrollado en EES, para validarlo conjuntamente con la bibliografa revisada.
En la figura 21 se muestra el diagrama de bloques del entorno de Hysys generado para el
modelo, cuyo proceso se detalla a continuacin.
Fig 21. Simulacin de la planta en Hysys.

Primero, se incorporan a la lista todos los componentes con los que se va a trabajar as como las
dos reacciones de equilibrio, Water-gas shift y Metanacin, que se introducirn en el gasificador.
Despus se define el paquete termodinmico que usar el programa y se asocia a las reacciones
definidas anteriormente. Se escoge el PRSV como mejora del Peng-Robinson.
Tras estos pasos se pasa ya a la pantalla de diseo. Se crean dos corrientes que se corresponden
con las de entrada al gasificador. La primera, a 25 y 101,3 i es la biomasa de entrada, cuyo
flujo msico y composicin son conocidas. La segunda, la del aire de entrada, se define con las
mismas caractersticas que la anterior, pero no se conoce su flujo. Estas dos corrientes se
incluyen en el gasificador como entradas, as como se crean otras dos corrientes, F3 y F4, de las
cuales no se conoce nada, como salidas. Es de resear el mdulo Set introducido entre estas dos
corrientes, de forma que hace las veces del Frg definido previamente, puesto que relaciona el
flujo molar de ambas corrientes. De este modo, se asegura la operacin del gasificador en
defecto de aire. De lo contrario, si la corriente F2 se dejara sin ninguna otra restriccin, se
producira realmente una combustin pues el aire de entrada sera el estequiomtrico, pues esa
es la solucin que da Hysys al aplicar los balances de materia.
El gasificador es un reactor de Gibbs, que simplemente opera minimizando la energa libre de
Gibbs del sistema a resolver aunque en el trabajo se ha intentado modificar ligeramente esas
condiciones para acercarlo ms al modelo de EES incluyndole como tipo las dos ecuaciones de
equilibrio.
La corriente F4, el gas de sntesis que sale por la parte superior del reactor, se lleva a un
intercambiador de calor, cuya corriente fra es agua a 20 y 101,3 i . El calor absorbido por la
misma se utiliza para precalentar la corriente de aire, colocando un Cooler (la corriente de agua
33
8. Anexos
caliente que sale del intercambiador de calor pasa a la temperatura inicial, 20) y un Heater
posteriormente (la energa del Cooler, el calor cedido, es el que se emplea para calentar la
corriente de aire). Adems es de resear que al existir una recirculacin, puesto que por
simplicidad se ha decidido disear un circuito cerrado de agua de refrigeracin, se debe colocar
un mdulo Recycle en la corriente de agua.
34
8. Anexos
8.4
ENTORNO DE EES
A continuacin se expone el cdigo generado en EES para el diseo del reactor de

gasificacin, separador ciclnico e intercambiador de calor respectivamente.
35
File:C:\Users\guillermo\Desktop\TFG\TFGGFM_GASIFICADOR.EES
16/09/2014 0:02:26 Page 1
EES Ver. 9.698: #3470: For use only by students and faculty in the Departamento de Ingenieria Mecanica Universidad de Za
~~~~~~~~~~~~~~~~~
************** PROCEDURES **********************
~~~~~~~~~~~~~~~~~
ESPERA CONVERSOR
Procedure CONVERSOR (alimentacion; C %ar; N %ar; H %ar; S %ar; O %ar; ash %ar; HHV ar : C F1; N F1; H F1; S F1; O F1; ash F1; F1; HHV F1)
Masas molares que entran en juego
MC := MolarMass C
[kg/kmol]
MH := MolarMass H
[kg/kmol]
MN := MolarMass N
[kg/kmol]
MS := MolarMass s
[kg/kmol]
MO := MolarMass O
[kg/kmol]
MW := MolarMass water
MA := MolarMass Air ha
[kg/kmol]
[kg/kmol]
1 Entrada reactor
C F1 := alimentacion
C %ar
100 MC
[kmol/h]
N F1 := alimentacion
N %ar
100 MN
[kmol/h]
H F1 := alimentacion
H %ar
100 MH
[kmol/h]
S F1 := alimentacion
S %ar
100 MS
[kmol/h]
O F1 := alimentacion
O %ar
100 MO
[kmol/h]
ash F1 := alimentacion
ash %ar
100
[kg/h]
F1 := C F1 + N F1 + H F1 + S F1 + O F1
[kmol/h]
HHV F1 := HHV ar
End CONVERSOR
ESPERA PORCENTAJES
Procedure PORCENTAJES F4 (CO F4; CO2 F4; CH4 F4; H2 F4; H2O F4; N2 F4; NH3 F4; SH2 F4 : CO%F4; CO2%F4; CH4%F4; H2%F4;
H2O%F4; N2%F4; NH3%F4; SH2%F4)
F4 := CO F4 + CO2 F4 + CH4 F4 + H2 F4 + H2O F4 + N2 F4 + NH3 F4 + SH2 F4
CO%F4 := 100
CO F4
F4
16/09/2014 0:02:26 Page 2
CO2%F4 := 100
CO2 F4
F4
CH4%F4 := 100
CH4 F4
F4
H2%F4 := 100
H2O%F4 := 100
N2%F4 := 100
H2 F4
F4
H2O F4
F4
N2 F4
F4
NH3%F4 := 100
NH3 F4
F4
SH2%F4 := 100
SH2 F4
F4
End PORCENTAJES F4
ESPERA MASAS
Procedure MASAS F4 (CO F4; CO2 F4; CH4 F4; H2 F4; H2O F4; N2 F4; NH3 F4; SH2 F4 : CO m4CO2 m4; CH4 m4; H2 m4; H2O m4; N2 m4;
NH3 m4; SH2 m4; m4)
Masas molares que entran en juego
MCO := MolarMass CarbonMonoxide
MCO2 := MolarMass CarbonDioxide
[kg/kmol]
[kg/kmol]
MCH4 := MolarMass Methane
[kg/kmol]
MH2 := MolarMass Hydrogen
[kg/kmol]
MH2O := MolarMass water
[kg/kmol]
MN2 := MolarMass Nitrogen
[kg/kmol]
MNH3 := MolarMass Ammonia
[kg/kmol]
MSH2 := MolarMass HydrogenSulfide

*************
CO m4 := MCO CO F4
[kg/h]
CO2 m4 := MCO2 CO2 F4
[kg/h]
CH4 m4 := MCH4 CH4 F4
[kg/h]
H2 m4 := MH2 H2 F4
[kg/h]
H2O m4 := MH2O H2O F4

N2 m4 := MN2 N2 F4
[kg/h]
[kg/h]
[kg/kmol]
16/09/2014 0:02:26 Page 3
NH3 m4 := MNH3 NH3 F4
[kg/h]
SH2 m4 := MSH2 SH2 F4
[kg/h]
m4 := CO m4 + CO2 m4 + CH4 m4 + H2 m4 + H2O m4 + N2 m4 + NH3 m4 + SH2 m4

End MASAS F4
~~~~~~~~~~~~~~~~~
************** DATOS DE PARTIDA **********************
~~~~~~~~~~~~~~~~~
Se introducen los datos de la biomasa a tratar, obtenido de la base ECN
COMPOSICIN: Porcentajes dados en el analisis inmediato, proximate analysis
W %ar
= 13,83
ash %dry
= 2,23
VM %daf
= 82,23
FC %daf
= 16,77
COMPOSICIN: Porcentajes dados en el analisis elemental, ultimate analysis

C %daf
= 48,77
H %daf
= 5,98
N %daf
= 0,59
S %daf
= 0,05
O %daf
= 44,55
Poderes calorficos [MJ/kg]

LHV daf
= 18,03 [MJ/kg]
HHV daf
= 19,35 [MJ/kg]
alimentacion = 600 [kg/h]

Datos de presiones y temperaturas
P atm
= 101,3 [kPa]
T amb
= 298 [K]
T ref
= 298 [K]
T F1
= 298 [K]
T F2
= 298 [K]
T F3
= 1523 [K]
[kg/h]
16/09/2014 0:02:26 Page 4
= 1073 [K]
T reactor
~~~~~~~~~~~~~~~~~~
************** CONVERSOR daf -> ar **********************
~~~~~~~~~~~~~~~~~~
Los clculos finales se deben realizar con la composicin de la madera tal y como se recibe (AR) y no como se mide (DAF).
Se pasan las composiciones a AR para lo que hace falta tenerlas tambin en DRY,AF&DAF
Balances ar
100 = W %ar + ash %ar + VM %ar + FC %ar

Balances dry
100 = ash %dry + VM %dry
+ FC %dry
ash %dry
= 100
ash %ar
+ VM %ar + FC %ar
VM %dry
= 100
ash %ar
VM %ar
ash %ar + VM %ar + FC %ar
Balances af
100 = W %af + VM %af + FC %af
W %af
VM %af
= 100
= 100
W %ar
W %ar
+ VM %ar + FC %ar
VM %ar
W %ar + VM %ar + FC %ar
Balances daf
VM %daf
= 100
VM %ar
VM %ar + FC %ar
******
Composiciones elementales AR
C %ar
= C %af
100 ash %ar

100
H %ar
= H %af
100 ash %ar

100
N %ar
= N %af
100 ash %ar

100
S %ar
= S %af
100 ash %ar

100
O %ar
= O %af
100 ash %ar

100
Composiciones elementales DRY
16/09/2014 0:02:26 Page 5
C %dry
= C %daf
100 ash %dry

100
H %dry
= H %daf
100 ash %dry

100
N %dry
= N %daf
100 ash %dry

100
S %dry
= S %daf
100 ash %dry

100
O %dry
= O %daf
100 ash %dry

100
Composiciones elementales AF
C %af
= C %daf
100 W %af
100
N %af
= N %daf
100 W %af
100
S %af
= S %daf
100 W %af
100
Se parte de 1kg de entrada en AF
H %af
10
H daf;inherente
H af;humedad
H daf;inherente
H %daf
100
+ H af;humedad
1000 10 W %af
W %af
9
= 10
100 = C %af + H %af + N %af + S %af + O %af

************
HHV ar
100 ash %ar

100
= HHV af
HHV dry
100 ash %dry

100
= HHV daf
HHV af
100 W %af
100
= HHV daf
***********
H2O;25C
= 2,44
[MJ/kg]
Balances ar
m ar;H2O
LHV ar
= 9
H %ar
100
W %ar
100
= HHV ar m ar;H2O H2O;25C
16/09/2014 0:02:26 Page 6
Balances dry
m dry;H2O
H %dry
100
= 9
= HHV dry m dry;H2O H2O;25C
LHV dry
Balances af
m af;H2O
= 9
H %af
100
W %af
100
= HHV af m af;H2O H2O;25C
LHV af
Call CONVERSOR alimentacion ; C %ar ; N %ar ; H %ar ; S %ar ; O %ar ; ash %ar ; HHV ar : C F1 ; N F1 ; H F1 ; S F1 ; O F1
; ash F1 ; F1 ; HHV F1
***************AJUSTE ESTEQUIOMTRICO DE LA COMBUSTIN DE LA BIOMASA***************
4 Gases salida
C F1
= CO2 F4;estq
H F1
= 2 H2O F4;estq
Balance a la especie C
O F1 + 2 O2 F2;estq
S F1
= SO2 F4;estq
N2 F2;estq
ER = 100
[kmol/h]
Balance a la especie H
[kmol/h]
= 2 CO2 F4;estq + H2O F4;estq + 2 SO2 F4;estq

Balance a la especie S
Balance a la especie O
[kmol/h]
[kmol/h]
79
O2 F2;estq
21
O2 F2
O2 F2;estq
****************REACCIN GASIFICACIN EN DEFECTO DE AIRE*******************

2 Entrada aire
N2 F2
79
O2 F2
21
[kmol/h]
Air = O2 F2 + N2 F2
3 Salida reactor a depsito de residuos
perdida%FC = 5
C F3
ash F3
Prdida de FC del 5% en masa del caudal de entrada
perdida%FC
C F1
100
= ash F1
[kmol/h]
[kg/h]
4 Syngas
C F1
= CO F4 + CO2 F4 + CH4 F4 + C F3
H F1
= 4 CH4 F4 + 2 H2 F4 + 2 H2O F4 + 3 NH3 F4 + 2 SH2 F4
Balance a la especie C
[kmol/h]
Balance a la especie H
[kmol/h]
16/09/2014 0:02:26 Page 7
N F1
= NH3 F4
2 N2 F2
S F1
[kmol/h]
Balance a la especie N (I)
= 2 N2 F4
= SH2 F4
Balance a la especie N (II)
Balance a la especie S
O F1 + 2 O2 F2
[kmol/h]
[kmol/h]
= CO F4 + 2 CO2 F4 + H2O F4
Balance a la especie O
[kmol/h]
********************
CLCULOS CONSTANTES DE EQUILIBRIO

Fo,xx: funciones de ENERGA LIBRE necesarias para los clculos
F o;C = 11,6 [kJ/kmol-K]
F o;CO
= 204,1 [kJ/kmol-K]
F o;CO2
= 226,4 [kJ/kmol-K]
F o;CH4
= 199,4 [kJ/kmol-K]
F o;H2
= 137 [kJ/kmol-K]
F o;H2O
= 196,7 [kJ/kmol-K]
difHm,xx=DELTAHm(T*)-DELTAHm(0K): diferencias de ENTALPA DE REACCIN necesarias para los clculos

difH m;C
= 1050 [kJ/kmol]
difH m;CO
= 8670 [kJ/kmol]
difH m;CO2
= 9360 [kJ/kmol]
difH m;CH4
= 10030 [kJ/kmol]
difH m;H2
difH m;H2O
= 8500 [kJ/kmol]
= 9910 [kJ/kmol]
REACCIN DE METANACIN (met): C+2H2<->CH4

H met;298K
F o;met
= 74900
[kJ/kmol]
= F o;CH4 F o;C 2 F o;H2
difH m;met
= difH m;CH4 difH m;C 2 difH m;H2
StMolarGibbs met
ln MET
K met
K met
= F o;met +
H met;298K
T reactor
difH m;met
T reactor
StMolarGibbs met
8,314 [kJ/kmol-K]
1
MET
CH4 F4
H2 F4
REACCIN DE WATER-GAS SHIFT (wgs): CO+H2O<->CO2+H2
16/09/2014 0:02:26 Page 8
H wgs;298K
F o;wgs
= 41200
[kJ/kmol]
= F o;CO2 + F o;H2 F o;CO F o;H2O
difH m;wgs
= difH m;CO2 + difH m;H2 difH m;CO difH m;H2O
StMolarGibbs wgs
ln WGS
= F o;wgs +
H wgs;298K
T reactor
difH m;wgs
T reactor
StMolarGibbs wgs
8,314 [kJ/kmol-K]
K wgs
1
WGS
K wgs
CO2 F4 H2 F4
CO F4 H2O F4
********************
Clculo de la entalpa aparente del fuel (biomasa) a partir de: -PCS=(nDELTAHf,productoscombustioncompleta)-DELTAHf,aparen
e,biomasa
h f;SO2
= 296900 [kJ/kmol]
HHV F1
1000
kJ
MJ
h f;aparente;biomasa
alimentacion = CO2 F4;estq h f;CO2 + H2O F4;estq h f;H2O;l + SO2 F4;estq h f;SO2
alimentacion
********************
BALANCE DE ENERGA
CALOR del reactor
perdidascalor% = 10
Q =
E F4
perdidascalor%
h f;aparente;biomasa
100
[kJ/kg]
= CO F4 H f;COF4
+ CO2 F4 H f;CO2F4
+ CH4 F4 H f;CH4F4
+ H2 F4 H f;H2F4
H f;H2OF4
+ N2 F4 H f;N2F4
+ NH3 F4 H f;NH3F4
+ SH2 F4 H f;SH2F4
+ H2O F4
ECUACIN BdE COMO: CALOR PERDIDO = SUMATORIO ENTALPAS DE PRODUCTOS CON SUS FLUJOS MOLARES
- SUMATORIO ENTALPAS DE REACTIVOS CON SUS FLUJOS MOLARES
Q alimentacion = E F4 + C F3 H f;C;F3
+ ash F3 h solidos
T F3 alimentacion h f;aparente;biomasa
[kJ/h]
********************
CLCULOS DE LAS ENTALPAS NECESARIAS PARA LOS BdE

ENTALPA DE LOS GASES CALCULADA COMO: ENTALPA DE FORMACIN COMPUESTO + ENTALPA A TEMPERATURA
GAS (TF4) - ENTALPA A TEMPERATURA REFERENCIA (298K)
COF4
16/09/2014 0:02:26 Page 9
H f;COF4
h f;CO
= h f;CO + h COF4;TF4
h COF4;T;ref
[kJ/kmol]
= 110400 [kJ/kmol]
h COF4;TF4
= h CarbonMonoxide ; T = T reactor ; P =P atm
h COF4;T;ref
= h CarbonMonoxide ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
CO2F4
H f;CO2F4
h f;CO2
= h f;CO2 + h CO2F4;TF4
h CO2F4;T;ref
[kJ/kmol]
= 393500 [kJ/kmol]
h CO2F4;TF4
= h CarbonDioxide ; T = T reactor ; P =P atm
h CO2F4;T;ref
= h CarbonDioxide ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
CH4F4
H f;CH4F4
h f;CH4F4
= h f;CH4F4
+ h CH4F4;TF4
h CH4F4;T;ref
[kJ/kmol]
= 74800 [kJ/kmol]
h CH4F4;TF4
= h Methane ; T = T reactor ; P =P atm
h CH4F4;T;ref
= h Methane ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
H2F4
H f;H2F4
h f;H2
= h f;H2 + h H2F4;TF4
h H2F4;T;ref
[kJ/kmol]
= 0
h H2F4;TF4
= h Hydrogen ; T = T reactor ; P =P atm
h H2F4;T;ref
= h Hydrogen ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
H2O(l)F4
H f;H2OF4
h f;H2O;l
= h f;H2O;l + h H2OF4;TF4
h H2OF4;T;ref
[kJ/kmol]
= 285500 [kJ/kmol]
h H2OF4;TF4
= h water ; T = T reactor ; P =P atm
h H2OF4;T;ref
= h water ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
N2F4
H f;N2F4
h f;N2
= h f;N2 + h N2F4;TF4
h N2F4;T;ref
[kJ/kmol]
= 0
h N2F4;TF4
= h Nitrogen ; T = T reactor ; P =P atm
h N2F4;T;ref
= h Nitrogen ; T = T ref ; P =P atm
NH3F4 (gas)
[kJ/kmol]
[kJ/kmol]
16/09/2014 0:02:26 Page 10
H f;NH3F4
h f;NH3
= h f;NH3 + h NH3F4;TF4
h NH3F4;T;ref
[kJ/kmol]
= 46200 [kJ/kmol]
h NH3F4;TF4
= h Ammonia ; T = T reactor ; P =P atm
h NH3F4;T;ref
= h Ammonia ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
SH2F4
H f;SH2F4
h f;SH2
= h f;SH2 + h SH2F4;TF4
h SH2F4;T;ref
[kJ/kmol]
= 20200 [kJ/kmol]
h SH2F4;TF4
= h HydrogenSulfide ; T = T reactor ; P =P atm
h SH2F4;T;ref
= h HydrogenSulfide ; T = T ref ; P =P atm
[kJ/kmol]
[kJ/kmol]
CF3
h solidos
= 0,84 [kJ/kg-K] ENTALPA DE LOS SLIDOS, OBTENIDA DE ARTCULO [26]
MC = MolarMass C
H f;C;F3
= h solidos
[kg/kmol]
MC T F3
[kJ/kmol]
N2F2
H f;N2F2
= h f;N2 + h N2F2;TF2
h N2F2;T;ref
[kJ/kmol]
h N2F2;TF2
= h Nitrogen ; T = T F2 ; P =P atm
[kJ/kmol]
h N2F2;T;ref
= h Nitrogen ; T = T ref ; P =P atm
[kJ/kmol]
O2F2
H f;O2F2
h f;O2
= 0
= h f;O2 + h O2F2;TF2
h O2F2;T;ref
[kJ/kmol]
[kJ/kmol]
h O2F2;TF2
= h Oxygen ; T = T F2 ; P =P atm
[kJ/kmol]
h O2F2;T;ref
= h Oxygen ; T = T ref ; P =P atm
[kJ/kmol]
Call PORCENTAJES F4 CO F4 ; CO2 F4 ; CH4 F4 ; H2 F4 ; H2O F4 ; N2 F4 ; NH3 F4 ; SH2 F4 : CO%F4 ; CO2%F4

; CH4%F4 ; H2%F4 ; H2O%F4 ; N2%F4 ; NH3%F4 ; SH2%F4
$export 'C:\Users\guillermo\Desktop\TFG\EXPORTARPORCENTAJES.txt' CO%F4;CO2%F4;CH4%F4;H2%F4;H2O%F4;N2%F4;N
3%F4;SH2%F4
Call MASAS F4 CO F4 ; CO2 F4 ; CH4 F4 ; H2 F4 ; H2O F4 ; N2 F4 ; NH3 F4 ; SH2 F4 : CO m4 ; CO2 m4 ; CH4 m4 ; H2 m4
; H2O m4 ; N2 m4 ; NH3 m4 ; SH2 m4 ; m4
$export 'C:\Users\guillermo\Desktop\TFG\EXPORTARMASAS.txt' m4
$export 'C:\Users\guillermo\Desktop\TFG\EXPORTARGASIFICADOR.txt' Patm;T,amb;T,ref;T,F1;T,F2;T,F3;T,reactor;CO,F4;CO2,F4;CH4,F4;H2,F
;H2O,F4;N2,F4;NH3,F4
CLCULOS PARA GRAFOS

RELATIVE BIOMASS/AIR RATIO
16/09/2014 0:02:26 Page 11
frg =
F1
Air
RATIO CO-H2
ratioCOH2 =
CO m4
H2 m4
File:C:\Users\guillermo\Desktop\TFG\TFGGFM_CICLON.EES
16/09/2014 20:46:35 Page 1
******* PROCEDURE WILKE PARA EL CLCULO DE LA VISCOSIDAD Y CONDUCTIVIDAD MEDIA DE UNA MEZCLA
DE GASES********
Procedure WILKE (F4; CO F4; CO2 F4; CH4 F4; H2 F4; H2O F4; N2 F4; NH3 F4; T reactor; P atm : media)
PROPIEDADES DE MATERIALES NECESARIAS EN EL PROCEDURE
MW 1 := MolarMass CarbonMonoxide
1 := Visc CarbonMonoxide ; T = T reactor ; P =P atm
x 1 :=
CO F4
F4
MW 2 := MolarMass CarbonDioxide
2 := Visc CarbonDioxide ; T = T reactor ; P =P atm
x 2 :=
CO2 F4
F4
MW 3 := MolarMass Methane
3 := Visc Methane ; T = T reactor ; P =P atm
x 3 :=
CH4 F4
F4
MW 4 := MolarMass Hydrogen
4 := Visc Hydrogen ; T = T reactor ; P =P atm
x 4 :=
H2 F4
F4
MW 5 := MolarMass water
5 := Visc water ; T = T reactor ; P =P atm
x 5 :=
H2O F4
F4
MW 6 := MolarMass Nitrogen
6 := Visc Nitrogen ; T = T reactor ; P =P atm
x 6 :=
N2 F4
F4
MW 7 := MolarMass Ammonia
7 := Visc Ammonia ; T = T reactor ; P =P atm
x 7 :=
NH3 F4
F4
Duplicate i=1;7
CLCULO MATRIZ DEL NMERO ADIMENSIONAL
Duplicate j=1;7
16/09/2014 20:46:35 Page 2
If
i = j
Then
PHI i ; j := 1
Else
16/09/2014 20:46:35 Page 3
PHI i ; j :=
1 +
mu i
mu j
1 +
MW i
MW j
1 /
MW j
MW i
1 /
EndIf
End
End
Duplicate i=1;7
CALCULO VECTOR DENOMINADOR
Duplicate j=1;7
den i :=
d j := x j PHI i ; j
W @} +
End
% En
951 ; 1 ; 7
End
RESULTADO FINAL
Duplicate i=1;7
media :=
mu media
W r 1
:= x i
mu i
den i
End
En
953 ; 1 ; 7
End WILKE
*************************************
******************FIN PROCEDURE WILKE*******************
$import 'C:\Users\guillermo\Desktop\TFG\exportargasificador.txt' /skiplines=15 Patm;T,amb;T,ref;T,F1;T,F2;T,F3;T,reactor;CO,F4;CO2,F4;CH4
F4;H2,F4;H2O,F4;N2,F4;NH3,F4
*************Clculo volmenes especificos en F4 necesarios para calcular QF4*************

vCOF4 = v CarbonMonoxide ; T = T reactor ; P =P atm
vCO2F4 = v CarbonDioxide ; T = T reactor ; P =P atm
vCH4F4 = v Methane ; T = T reactor ; P =P atm
vH2F4 = v Hydrogen ; T = T reactor ; P =P atm
vH2OF4 = v water ; T = T reactor ; P =P atm
vN2F4 = v Nitrogen ; T = T reactor ; P =P atm
vNH3F4 = v Ammonia ; T = T reactor ; P =P atm
~~~~~~~~
************** CICLN **********************
~~~~~~~
Dimensiones fijadas
16/09/2014 20:46:35 Page 4
Ne
= 8
F4;p
Nmero de giros efectivos. Demostrado a partir de su definicin

[kg/m3] Densidad de las partculas de entrada al cicln [14]
= 1050
Q F4
0,000277778
h
s
CO F4 vCOF4 + CO2 F4 vCO2F4 + CH4 F4 vCH4F4 + H2 F4 vH2F4 + H2O F4
vH2OF4 + N2 F4 vN2F4 + NH3 F4 vNH3F4

Caudal gas de entrada al cicln: [m3/s]=kmol/h*1h/3600s*m3/kmol
*******************************
Viscosidad del gas de entrada al cicln PROCEDURE WILKE
F4 = CO F4 + CO2 F4 + CH4 F4 + H2 F4 + H2O F4 + N2 F4 + NH3 F4
Call WILKE F4 ; CO F4 ; CO2 F4 ; CH4 F4 ; H2 F4 ; H2O F4 ; N2 F4 ; NH3 F4 ; T reactor ; P atm : media
*******************************
Dimensionamiento
Do Dimetro total
Wi =
Do
4
De =
Do
2
H =
Do
2
H1 = 2 Do
S = 0,625 D o
H2
= 2 Do
Dd
Do
4
Parmetros
Ne
A ciclon
H1 + H2
2
= Wi H
u =
Q F4
A ciclon
D p50
H1 +
= 90
Velocidad del gas a la entrada del cicln

9 media W i
2 3,142 N e u F4;p
1 /
Dimetro de corte
16/09/2014 20:46:35 Page 5
i
100
=
1 +
D p\D;p50
D pi
D p50
D pi
D p50
D pi
D p50
Defino la variable Dp\D,p50 para poder representar en un grafo la eficacia del cicln en funcin del tamao
de partcula
File:C:\Users\guillermo\Desktop\TFG\TFGGFM_INTERCAMBIADORDECALOR.EES
18/09/2014 18:51:36 Page 1
******* PROCEDURE NUSSELTF PARA EL CLCULO DEL NUSSELT EN EL INTERIOR DE LOS TUBOS********
Procedure NUSSELTF (Re f; D i; Pr f; L t : Nusselt f)
If
Re f < 2300
Then
Rgimen LAMINAR
x cdh := 0,056 Re f D i
Longitud entrada hidrulica
x cdt := 0,043 Re f Pr f D i
Gz f :=
If
Longitud entrada trmica
Re f Pr f D i
Lt
x cdt > 10 x cdh
Then
Entrada puramente trmica
Nusselt f := 4,36 + 0,023 ln 1 +

Else
0,0012 Gz f
Correlacin de Hausen
Correlacin de Kays
Entrada combinada trmica/hidrulica
Nusselt f := 4,36 + 0,036 ln 1 +
0,0011 Gz f
EndIf
Else
Rgimen TURBULENTO
:= 0,00004
f :=
[m] Rugosidad de la pared del tubo interior de acero comercial
0,8685 ln
1,964 ln Re f
Re f
f
Nusselt f :=
3,8215
1 + 12,7
Re f 1000
f
8
3,71 D i
Pr f
0,5
Pr f
2 /
1 +
Factor de friccin
Di
Lt
2 /
Correlacin de Gnielinski
para rgimen
turbulento
EndIf
End NUSSELTF
******************FIN PROCEDURE NUSSELTF*******************

*************************************
******* PROCEDURE WILKE PARA EL CLCULO DE LA VISCOSIDAD Y CONDUCTIVIDAD MEDIA DE UNA MEZCLA
DE GASES********
Procedure WILKE (F4; CO F4; CO2 F4; CH4 F4; H2 F4; H2O F4; N2 F4; NH3 F4; T c;m; P atm : c)
PROPIEDADES DE MATERIALES NECESARIAS EN EL PROCEDURE
MW 1 := MolarMass CarbonMonoxide
1 := Visc CarbonMonoxide ; T = T c;m ; P =P atm
x 1 :=
CO F4
F4
18/09/2014 18:51:36 Page 2
MW 2 := MolarMass CarbonDioxide
2 := Visc CarbonDioxide ; T = T c;m ; P =P atm
x 2 :=
CO2 F4
F4
MW 3 := MolarMass Methane
3 := Visc Methane ; T = T c;m ; P =P atm
x 3 :=
CH4 F4
F4
MW 4 := MolarMass Hydrogen
4 := Visc Hydrogen ; T = T c;m ; P =P atm
x 4 :=
H2 F4
F4
MW 5 := MolarMass water
5 := Visc water ; T = T c;m ; P =P atm
x 5 :=
H2O F4
F4
MW 6 := MolarMass Nitrogen
6 := Visc Nitrogen ; T = T c;m ; P =P atm
x 6 :=
N2 F4
F4
MW 7 := MolarMass Ammonia
7 := Visc Ammonia ; T = T c;m ; P =P atm
x 7 :=
NH3 F4
F4
Duplicate
Duplicate
i=1;7 j=1;7
If i =DEL
j Then
PHI i ; j := 1
CLCULO MATRIZ
NMERO ADIMENSIONAL
PHI i ; j :=
1
1 +
MW i
MW j
End
Duplicate
End Duplicate
i=1;7 j=1;7
d j := xDENOMINADOR
CALCULO VECTOR
PHI i ; j
j
End
1 +
mu i
mu j
Else
1 /
MW j
MW i
1 /
EndIf
18/09/2014 18:51:36 Page 3
den i :=
W ! \
En951 ; 1 ; 7
End
Duplicate
i=1;7
RESULTADO FINAL mu i
mu c i := x i
den i
End
c :=
W ! \
En953 ; 1 ; 7
End WILKE
******************FIN PROCEDURE WILKE*******************

*************************************
******* PROCEDURE PROPIEDADES MEDIAS DEL GAS. FLUIDO EXTERIOR********
Procedure PROPIEDADESC (CO%F4; CO2%F4; CH4%F4; H2%F4; H2O%F4; N2%F4; NH3%F4; SH2%F4; T c;m; P atm : c;
Pr c)
CONDUCTIVIDAD
CO := k CarbonMonoxide ; T = T c;m ; P =P atm
CO2 := k CarbonDioxide ; T = T c;m ; P =P atm
CH4 := k Methane ; T = T c;m ; P =P atm
H2 := k Hydrogen ; T = T c;m ; P =P atm
H2O := k water ; T = T c;m ; P =P atm
N2 := k Nitrogen ; T = T c;m ; P =P atm
NH3 := k Ammonia ; T = T c;m ; P =P atm
SH2 := k HydrogenSulfide ; T = T c;m ; P =P atm
c :=
CO CO%F4 + CO2 CO2%F4 + CH4 CH4%F4 + H2 H2%F4 + H2O H2O%F4 + N2 N2%F4 + NH3
100
PRANDTL
Pr CO := Pr CarbonMonoxide ; T = T c;m ; P =P atm
Pr CO2
:= Pr CarbonDioxide ; T = T c;m ; P =P atm
Pr CH4
:= Pr Methane ; T = T c;m ; P =P atm
Pr H2 := Pr Hydrogen ; T = T c;m ; P =P atm

Pr H2O
:= Pr water ; T = T c;m ; P =P atm
Pr N2 := Pr Nitrogen ; T = T c;m ; P =P atm
18/09/2014 18:51:36 Page 4
Pr NH3
:= Pr Ammonia ; T = T c;m ; P =P atm
Pr SH2
:= Pr HydrogenSulfide ; T = T c;m ; P =P atm
Pr c :=
Pr CO CO%F4 + Pr CO2 CO2%F4 + Pr CH4 CH4%F4 + Pr H2 H2%F4 + Pr H2O H2O%F4 + Pr N2 N2%F4 + Pr
100
End PROPIEDADESC
******************FIN PROCEDURE PROPIEDADESC*******************

*************************************
$import 'C:\Users\guillermo\Desktop\TFG\exportargasificador.txt' /skiplines=15 Patm;T,amb;T,ref;T,F1;T,F2;T,F3;T,reactor;CO,F4;CO2,F4;CH4,F4
H2,F4;H2O,F4;N2,F4;NH3,F4
$import 'C:\Users\guillermo\Desktop\TFG\exportarporcentajes.txt' /skiplines=9 CO%F4;CO2%F4;CH4%F4;H2%F4;H2O%F4;N2%F4;N

3%F4;SH2%F4
DATOS DE PARTIDA
FLUIDO CALIENTE. GASES DE SALIDA DEL REACTOR. EXTERIOR DE LOS TUBOS
T c;e = T reactor
T c;s = 473 [K]
T c;m =
T c;e + T c;s
2
Viscosidad del gas de entrada al Intercambiador de calor PROCEDURE WILKE

F4 = CO F4 + CO2 F4 + CH4 F4 + H2 F4 + H2O F4 + N2 F4 + NH3 F4
Call WILKE F4 ; CO F4 ; CO2 F4 ; CH4 F4 ; H2 F4 ; H2O F4 ; N2 F4 ; NH3 F4 ; T c;m ; P atm : c
cp c = Cp Air ha ; T = T c;m ; P =P atm
FLUIDO FRO. AGUA A EVAPORAR
T f;e = 293 [K]
mf = 1
[kg/s]
u f = 1,3 [m/s] MARN

f
= k water ; T = T f;m ; P =P atm
= Visc water ; T = T f;m ; P =P atm
v f = v water ; T = T f;m ; P =P atm

cp f = Cp water ; T = T f;e ; P =P atm
MW
f
= MolarMass water
= water ; T = T f;m ; P =P atm
18/09/2014 18:51:36 Page 5
Pr f = Pr water ; T = T f;m ; P =P atm

Balances de Energa
q BdE
q BdE
0,000277778
mf
cp f
MW
h
s
F4 cp c
T c;e T c;s
T f;s T f;e
T f;e + T f;s
2
T f;m =
DISEO PREVIO
REA SUPERFICIAL EXTERIOR ESTIMADA
A sup;ext;est
1000
U est
= 150
ml
= F ml;cc
ml;cc
q BdE
U est ml
Se estima un valor de coeficiente global de transferencia, para hallar un rea estimada de intercambio.
MARN p.10
T c;s T f;e
W
kW
ln
T c;e T f;s
T c;s T f;e
T c;e T f;s
F = 1
REA SUPERFICIAL EXTERIOR POR UNIDAD DE LONGITUD DE TUBO
A sup;ext
1 N aletas;L t
3,142 D e + N aletas;L
3,142
2
Da
De
+ 3,142 D a t
Dimensiones
D i = 0,0138 [m] Dimetro interior de los tubos
De
= 0,0164 [m] Dimetro exterior de los tubos
Da
= 0,022 [m] Dimetro exterior de las aletas. El interior es igual a De
La
Da De
2
t = 0,000254 [m] Espesor de las aletas

= 100 [m-1] Nmero de aletas por unidad de longitud del tubo
N aletas;L
s =
1 N aletas;L t
N aletas;L
= 0,1
=
Distancia entre aletas
sigma=Amin,transv/Asup,ext; Amin,transv es rea mnima transversal & Asup,ext es el rea total de intercambio de calor
A min;transv
A sup;ext
18/09/2014 18:51:36 Page 6
A c;min
= A min;transv
Lt
k = k Aluminum ; T = T c;m
Clculo del banco de tubos y de la carcasa
ST
= 1,5 D e
Distancia entre centro de los tubos en disposicin tringulo
Nmero de tubos
N t = Ceil
A sup;ext;est
A sup;ext L t
Nmero de tubos por ancho IdQ
Nmero de circuitos
N c = Ceil 4
mf
f MW v f 3,142 D i
Carcasa
N t;fila
= Nc
MARN. Se puede empezar considerando que cada circuito se encuentra en posicin horizontal, de forma que el nmero
de tubos por fila coincida con el nmero de circuitos
N t;fila S T = L t
NL
= Ceil
Nt
N t;fila
Nmero de filas de tubos
RECLCULO Y DISEO DEFINITIVO

COEFICIENTE DE CONVECCIN DEL CALOR EN EL INTERIOR DE LOS TUBOS
hf =
Nusselt f f
Di
Re f = u f D i
MW
f vf
Call NUSSELTF Re f ; D i ; Pr f ; L t : Nusselt f

COEFICIENTE DE CONVECCIN DEL CALOR EN EL EXTERIOR DEL BANCO DE TUBOS
hc =
Nusselt c c
De
Call PROPIEDADESC CO%F4 ; CO2%F4 ; CH4%F4 ; H2%F4 ; H2O%F4 ; N2%F4 ; NH3%F4 ; SH2%F4 ; T c;m ; P atm
: c ; Pr c
Nusselt c = 0,134 Re c
Re c =
0,000277778
0,681
h
s
Pr c
1 /
s
La
0,2
s
t
0,11
Correlacin de Briggs y Young para disposicin

del banco de tubos en tringulo
m4 D e
A c;min c
$import 'C:\Users\guillermo\Desktop\TFG\exportarmasas.txt' /skiplines=2 m4

CLCULO COEFICIENTE GLOBAL DE TRANSFERENCIA DE CALOR, U. Ecuacin general: 1/Ui*Ai=Rtotal=1/(etao,c*h,c
18/09/2014 18:51:36 Page 7
A,c)+R''e,c/(etao,c*A,c)+Rp+R''e,f/(etao,f*A,f)+1/(etao,f*h,f*A,f)
Se desprecia la resistencia de enfriamiento del fluido fro (lado interior de los tubos) pero no as la del fluido caliente (gas
con restos de partculas y cenizas), quedando
1
Uc
o;c h c
R'' e;c
o;c
+ Ac Rp +
Ac
o;f h f A f
A c = A sup;ext L t
A f = 3,142 D i L t
Aa
3,142
2
Da
De
+ 3,142 D a t
N aletas;L L t
Eficiencia de una superficie con aletas. En el lado fro es etao,f=1 porque no hay aletas
o;f
= 1
o;c = 1
Aa
Ac
1 a
Eficiencia de una aleta. Se calcula con las funciones de Bessel
m =
1 /
hc
c t
mr1 =
m De
2
mr2 =
m Da
2
De
=
Da
2
De
2
K1
mr1
I1
mr2
I1
mr1
K1
mr2
K0
mr1
I1
mr2
+ I0
mr1
K1
mr2
Resistencia de enfriamiento del fluido caliente

R'' e;c = 0,0004 [m2*K/W]
Resistencia de conduccin de la pared que separa ambos fluidos
ln
Rp
De
Di
2 3,142 L t k

Gasificadores

Cargado por

Copyright:

Formatos disponibles

Gasificadores

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Gasificadores

Cargado por

Copyright:

Formatos disponibles

Trabajo Fin de Grado en

Ingeniera de Tecnologas Industriales

Guillermo Figueruelo Malo

Alberto Gonzalo Callejo

Escuela de Ingeniera y Arquitectura. Universidad de Zaragoza.

TTULO: Desarrollo de modelos matemticos y diseo de los principales equipos

AUTOR: Guillermo Figueruelo Malo

DIRECTOR: Alberto Gonzalo Callejo

DEPARTAMENTO: Ingeniera Qumica y Tecnologas del Medio Ambiente

Desarrollo de modelos matemticos y diseo de los principales equipos de una

OBJETIVO Y ALCANCE DEL PROYECTO.............................................................. 1

CONTEXTO GLOBAL ........................................................................................... 2

TIPOS DE BIOMASA ............................................................................................. 2

VENTAJAS DEL USO ENERGTICO DE LA BIOMASA ............................... 3

INCONVENIENTES DEL USO ENERGTICO DE LA BIOMASA................ 4

PROCESO DE GASIFICACIN DE BIOMASA ................................................ 5

EQUIPOS DE LA PLANTA ........................................................................................... 6

REACTOR DE GASIFICACIN .......................................................................... 6

SEPARADOR CICLNICO .................................................................................. 7

DESARROLLO DE LOS MODELOS MATEMTICOS ........................................... 9

REACTOR DE GASIFICACIN ........................................................................ 10

SEPARADOR CICLNICO ................................................................................ 17

TABLAS PARA EL CLCULO DE LAS CONSTANTES DE EQUILIBRIO30

ENTORNO DE ABAQUS ..................................................................................... 31

ENTORNO DE HYSYS ......................................................................................... 33

1. OBJETIVO Y ALCANCE DEL PROYECTO

1. Objetivo y alcance del proyecto

1. OBJETIVO Y ALCANCE DEL PROYECTO

El incremento demogrfico y el constante aumento de la produccin de bienes y servicios,

Fig 1. Prospectiva de potencia instalada y de produccin de energa elctrica en Aragn en

Habitualmente la biomasa es clasificada en dos grandes grupos: primaria y residual [3]. La

Residuos agrcolas, forestales y cultivos energticos. Los residuos agrcolas se obtienen

VENTAJAS DEL USO ENERGTICO DE LA BIOMASA

- La sostenibilidad ambiental en todos aquellos aprovechamientos de biomasa residual forestal

INCONVENIENTES DEL USO ENERGTICO DE LA BIOMASA

PROCESO DE GASIFICACIN DE BIOMASA

Entre todos los procesos de conversin de biomasa, la gasificacin es uno de los ms

Fig 2. Fases existentes en un proceso de gasificacin [3].

Los gasificadores se utilizan para transformar la biomasa en un gas combustible. El movimiento

Lecho mvil: el slido se mueve lentamente.

Gasificador downdraft o de corrientes paralelas: el aire entra por la parte

Fig 3. Gasificador downdraft y updraft [6]

Lecho fluidizado: las partculas slidas se encuentran en suspensin gracias a una

Circulante: la fluidizacin es rpida, sin burbujas y se recircula toda la material

Fig 4. Gasificadores de lecho fluidizado: (a) Burbujeante; (b) Circulante [7].

4. DESARROLLO DE LOS MODELOS

4. Desarrollo de los modelos matemticos

4. DESARROLLO DE LOS MODELOS MATEMTICOS

A continuacin, se exponen los distintos cdigos de diseo de equipos realizados en EES,

4. Desarrollo de los modelos matemticos

Fig 5. Esquema del reactor de gasificacin.

En la base de datos de ECN se expresan las composiciones en masa (% en peso, expresado en la

4. Desarrollo de los modelos matemticos

As Received (ar), que expresa a la biomasa en su totalidad tal y como se encuentra en la

Fig 6. Base de datos para biomasa y residuos forestales de ECN.

4. Desarrollo de los modelos matemticos

Tabla 2. Bases empleadas en la expresin de la composicin de madera.

Anlisis inmediato (% wt)