CAP 7 REACTORES

Los reactores químicos, particularmente para procesos continuos, son a menudo diseñados a
medida para involucrar múltiples fases (por ejemplo, vapor, líquido, sólido reactivo y catalizador
sólido), diferentes geometrías (p. ej., tanques agitados, fl ujos tubulares, convergentes y divergentes
boquillas, flujos en espiral y transporte de membrana), y varios regímenes de impulso, calor y
transferencia de masa (p. ej., viscoso flujo, flujo turbulento, conducción, radiación, difusión y
dispersión).

Hay tantas configuraciones, que implican diferentes combinaciones de estos atributos, que intenta
desarrollar los modelos de reactores generalizados han tenido un éxito limitado. La mayoría de los
simuladores de procesos proporcionan cuatro tipos de reactores. modelos, que incluyen: (1) un
modelo estequiométrico que permite especificación de conversiones de reactivos y extensiones de
reacción para una o más reacciones especificadas; (2) un modelo para múltiples fases (vapor, líquido
y sólido) en equilibrio químico, donde el enfoque del equilibrio para las reacciones individuales
puede ser especificado; (3) un modelo cinético para un tanque de agitación continua reactor (CSTR)
que supone una mezcla perfecta de homogéneo fases (líquido o vapor); y (4) un modelo cinético
para un flujo de enchufe reactor tubular (PFTR o PFR), para fases homogéneas (líquido o vapor) y
suponiendo que no haya mezcla inversa (dispersión).

Estos modelos ideales se utilizan en las primeras etapas del proceso. síntesis, cuando los detalles de
los diseños del reactor son menores importante pero los efluentes del reactor y las tareas de calor
son necesarios. Los modelos de reactores ideales se reemplazan por modelos personalizados
modelos a medida que los detalles adquieren significación. Para este propósito, todos los
simuladores de hoja de cálculo proporcionan facilidades para la inserción de modelos generados por
el usuario. Estos se refinan a menudo como el diseño procede y como datos del reactor del
laboratorio o planta piloto regresan, con algunos de los simuladores que proporcionan instalaciones
para estimar los parámetros de modelos cinéticos por no lineal regresión.

Cuando se trabaja con los modelos de reactores ideales, el lector debe consultar los libros de texto
disponibles sobre análisis de reactores y diseño, por ejemplo,
ElementsofChemicalReactionEngineering por H.S. Fogler (2005), The Engineering of Chemical
Reacciones de L.D. Schmidt (1998), Ingeniería de reacción química por O. Levenspiel (1999), y para
los manuales de usuario y segmentos tutoriales de los simuladores de hoja de cálculo. El
seguimiento la discusión de los modelos de reactores ideales utilizados en simuladores es precedido
por una breve revisión de la estequiometría de reacción y grado de reacción, que en conjunto
proporcionan la base para cálculo del modelo de reacción de conversión en los simuladores.

Modelos avanzados para reactores tubulares que representan no enchufes el flujo se introduce en la
Sección 7.5

Estequiometria de reacción:

Para la mayoría de los modelos de reactores en los simuladores de hoja de cálculo, es necesario para
proporcionar reacciones químicas R que involucren C especies quimicas:
Donde Aj es la fórmula química para las especies j y nij es el coeficiente estequiométrico para la
especie j en la reacción i (negativo para reactivos, positivo para productos). Como ejemplo, para el
fabricación de metanol, que se ordenen los productos químicos de acuerdo con la volatilidad
decreciente, es decir, (1) H2, (2) CO, y (3) CH3OH. La reacción se puede escribir como -2H2-
CO+CH3OH= 0, y el coeficiente estequiométrico matriz es:

Alcance de la reacción:

Considera una sola reacción. En el reactor modelos estequiométricos, uno especifica la conversión
fraccional, Xk, del reactivo clave k,

donde nk; en y nk; fuera están los lunares de las especies k entrando y saliendo del reactor y 0 Xk 1,
o la extensión (número de extensión de moles) de la reacción i:

está especificado. Las tasas de flujo molar de los componentes en el el efluente del reactor, nj, out,
se calcula a partir del componente caudales molares en la alimentación del reactor, nj, in, por un
componente ecuación de balance de materiales que es consistente con la reacción estequiometría Si
una especificación de la conversión fraccional del componente clave, k, está hecho con la ecuación:

Si se especifica el alcance de la reacción:

Por ejemplo, para la conversión de CO y H2 a CH3OH, suponiendo una alimentación inicial de 100
kmol / h CO y 600 kmol / h Conversión de H2 y 70% de CO (el componente clave), utilizando Eq.
(7.4), las tasas de flujo molar de los tres componentes en el el efluente del reactor es:

La fracción molar de metanol en el efluente del reactor es Ych3OH=70/(460+30+70)=0.125. Si, en

cambio, el alcance de la reacción se especifica como 70 kmol / h, la ecuación. (7.5) da los mismos
resultados.
Para reacciones múltiples, las reacciones deben especificarse como serie o paralelo. El primero es
equivalente a tener reactores en serie con la alimentación a cada reactor excepto el primero, siendo
el producto del reactor anterior. Cada reacción puede tener un reactivo clave diferente. Para
reacciones paralelas, es preferible para especificar el alcance de la reacción para cada reacción, que
resultados en:

Equilibrio:

Una reacción química se puede escribir como un estequiométrico general ecuación, en términos de
reactivos A, B, etc., y productos R, S, etc.,

Al escribir esta ecuación, es muy importante que, a menos que de lo contrario, cada reactivo y
producto se entiende que ser el componente puro en una fase separada y designada: gas, líquido o
sólido. Una reacción se caracteriza por dos cantidades termodinámicas importantes, a saber, el calor
de reacción y la energía de reacción de Gibbs (libre). Además, estas dos cantidades son funciones de
temperatura y presión. Los datos termodinámicos están ampliamente disponibles en simuladores y
en otros lugares por más de mil componentes, para el cálculo de estas dos cantidades bajo estándar
condiciones de estado, por ejemplo, a una temperatura de referencia de 25 C y 1 barra con todos los
componentes en una fase designada, generalmente como un gas ideal. El efecto de la temperatura
sobre el calor de reacción depende de la capacidad calorífica de los reactivos y productos y el efecto
de la temperatura en esas capacidades de calor.

Para muchas reacciones, el efecto de la temperatura en el calor de reacción es relativamente

pequeño. Por ejemplo, la reacción de CO y H2 para formar metanol gaseoso,

tiene un calor de reacción estándar, ΔHrxn a 25 C de -90; 400 kJ / kmol de metanol, mientras que a
800 C, el calor de reacción, ΔHrxn; es -103; 800 kJ / kmol, un relativamente pequeño cambiar por un
cambio tan grande en la temperatura. Por el contrario, el efecto de la temperatura sobre la energía
de reacción de Gibbs puede ser muy grande Por ejemplo, para la misma formación de metanol
reacción, la energía de reacción estándar de Gibbs, ΔGrxn, es -25; 200 kJ / kmol a 25 C. Ya, a 500 C, el
Gibbs energía de reacción, ΔGrxn, ha sufrido un cambio drástico a 88; 000 kJ / kmol.

Muchas reacciones de importancia industrial están limitadas por equilibrio químico, con conversión
parcial de la limitación reactivo y, con la velocidad de la reacción inversa igual a la velocidad de la
reacción hacia adelante. Para una composición de alimentación especificada y temperatura y presión
finales, la composición del producto en equilibrio químico puede calcularse por cualquiera de dos
métodos: (1) constantes de equilibrio químico (valores K) calculado a partir de la energía de reacción
de Gibbs combinada con ecuaciones de balance de materiales para un conjunto de independientes
reacciones, o (2) la minimización de la energía de Gibbs de la sistema de reacción El primer método
es aplicable cuando el la estequiometría se puede especificar para todas las reacciones que se
consideran. El segundo método requiere solo una lista de posibles productos para el primer método,
una constante de equilibrio químico, K, se calcula para cada reacción estequiométrica independiente
en el conjunto, usando la ecuación:
donde ai es la actividad componente. Para una solución de gas, la actividad está dada por:

donde Φi, es el coeficiente de fugacidad del componente i en el gas mezcla, igual a 1.0 para una
solución de gas ideal, y Pi es el presión parcial. En general, fi es una función de T, P y composición. A
presiones bajas a moderadas, Φi=

1: 0, para que la actividad es igual a la presión parcial en bar. Es común para reemplazar la actividad
en la ecuación para K con lo anterior ecuación para dar:

donde KΦ= 1: 0 para presiones bajas a moderadas. Para una solución líquida, la actividad está dada
por:

donde yi es el coeficiente de actividad del componente i en el mezcla líquida y es igual a 1.0 para una
solución líquida ideal, Vi es el volumen molar parcial del componente i, y Ps yo es el presión de vapor
del componente i. Las presiones están en barra. En general, gi es una función de T, P y composición.
Ideal soluciones líquidas a presiones bajas a moderadas, yi= 1: 0, entonces que la actividad es igual a
la fracción molar. Es común a reemplazar las actividades en la ecuación. (7.9) con la ecuación. (7.12)
para dar:

donde Ky= 1: 0 para soluciones líquidas ideales a bajo a moderado presiones La mayoría de los libros
de texto sobre termodinámica química cuadros actuales de log10K en función de la temperatura
para

Muchas reacciones químicas. La ecuación de vant Hoff relaciona K a temperatura por

Si se supone que el calor de reacción es independiente de la temperatura sobre un rango particular

de temperatura, integración y conversión al formulario log10 da la correlación aproximada ecuación:
donde T es la temperatura absoluta. Muchos equilibrio químico las curvas se representan con una
precisión razonable por esta ecuación. Por ejemplo, para la reacción en fase gaseosa de CO y H2
para formar metanol, en un rango de temperatura de 273 a 773 K,

Típicamente, la reacción de síntesis de metanol es catalizada por óxido de cobre-zinc, a una presión
de 100 bar y una temperatura de 300 C. Se utiliza un gran exceso de hidrógeno para ayudar a
absorber El relativamente alto calor de reacción. En estas condiciones, K= 0.0002202 y KΦ=0.61. Por
lo tanto,

Si X es la conversión fraccional de equilibrio de la limitación reactivo, CO, luego usando la misma

composición de alimentación inicial que antes y la estequiometría para la reacción, el equilibrio las
fracciones molares son:

Combinando las cuatro ecuaciones anteriores para dar un no lineal ecuación en X y resolver da X=
0.7087. El segundo método para calcular el equilibrio químico es aplicar el criterio de que la energía
total de Gibbs, G, es mínimo a temperatura y presión constantes. Alternativamente, uno podría usar
la entropía, S, como máximo o el Energía de Helmholtz, A, como mínimo, pero la energía de Gibbs es
más ampliamente aplicado. Dos ventajas de este segundo método son: (1) la evitación de tener que
formular estequiométricos ecuaciones (solo los productos posibles deben especificarse), y (2) la
facilidad de formulación para múltiples fases y Equilibrio de fase simultánea. Para una sola fase, el
total la energía de Gibbs en una T y P especificada viene dada por:

Donde Ni es el número molar del componente i Gi, es el Energía de Gibbs molar parcial del
componente i en el equilibrio mezcla. Los componentes son los del feed plus aquellos que pueden
ser producidos por reacciones químicas. La energía de Gibbs se minimiza con respecto al Ni, que
comparte limitado por los equilibrios atómicos. El método es fácilmente extendido a múltiples fases.
Parecería que por simplicidad o utilidad, el segundo se preferiría el método porque cuando se utiliza
este modelo, no es necesario especificar un conjunto independiente de reacciones químicas. Por lo
tanto, el diseñador no está obligado a identificar las reacciones que tienen lugar Sin embargo, dado
que la mayoría de los reactores están diseñados para enfatizar las reacciones deseadas y reducir o
excluir indeseados reacciones, las reacciones químicas que tienen lugar en el reactor generalmente
se conocen cuando se diseña el reactor.
Por lo tanto, se puede preferir el primer modelo, con el segundo modelo siendo útil para la
exploración preliminar de la termodinámica posibilidades Para que el segundo modelo sea más útil
se deben aplicar restricciones a ciertas reacciones improbables. Si esto no se hace, el segundo
método puede producir resultados que son incorrectos porque implícitamente requieren reacciones
que son no es cinéticamente factible.

Para instrucciones sobre el uso del equilibrio modelos de reactor constante y Gibbs en el proceso
simuladores, ver los módulos multimedia, que se puede descargar desde el sitio web de Wiley
asociado con este libro (ASPEN Químico Reactores)

Cinética:

La conversión fraccional y los modelos de reactor de equilibrio son útiles en las primeras etapas del
diseño del proceso al conducir material y estudios de balance energético. Sin embargo,
eventualmente los sistemas de reactores debe estar configurado y dimensionado. Esto requiere
conocimiento de cinética de reacción, que se obtiene mediante la realización de laboratorio
experimentos para reacciones no catalíticas homogéneas, powerlaw las expresiones se usan
comúnmente para la regresión de laboratorio datos cinéticos Estas expresiones no siempre se basan
en el ecuación estequiométrica porque varias reacciones elementales los pasos pueden estar
involucrados, la suma de los cuales es el estequiométrico ecuación, pero una de las cuales puede
controlar la velocidad de reacción global. Los pasos de reacción elementales rara vez involucran más
de dos moléculas la ecuación cinética general de la ley de potencia es:

Dónde -r j es la tasa de desaparición del componente j (en mol / tiempo-volumen), Ci es la

concentración del componente i (en mol / volumen), t es tiempo, k es el coeficiente de velocidad de
reacción, αi es el orden de reacción con respecto al componente i, y C es el número de
componentes.

Para reacciones en fase gaseosa, la presión parcial, Pi, es algo veces utilizado en lugar de la
concentración, Ci, en el ecuación cinética El coeficiente de velocidad de reacción es una función de
temperatura dada por el empírico Arrhenius ecuación:

donde ko es el factor pre-exponencial, y E es la activación energía. Para reacciones catalizadas por

catalizador poroso sólido partículas, la secuencia de pasos elementales puede incluir adsorción en la
superficie del catalizador de uno o más reactivos y / o desorción de uno o más productos. En ese
caso, un La ecuación cinética de Langmuir-Hinshelwood (LH) se encuentra a menudo para ajustar los
datos cinéticos experimentales con mayor precisión que el expresión de la ley de poder de la
ecuación. (7.19). La formulación de LH es caracterizado por un término denominador que incluye
concen- Traciones de ciertos reactivos y / o productos que son fuertemente adsorbido en el
catalizador. La ecuación LH también puede incluir un prefijo, h, llamado factor de efectividad
general, que explica para resistencias de transferencia de masa y calor, tanto externas como interno,
a las partículas de catalizador. Como ejemplo, laboratorio Los datos cinéticos para la oxidación al aire
de SO2 a SO3 se ajustan bien por la siguiente ecuación de LH:
CSRT:

El modelo de reactor cinético más simple es el CSTR (tanque de agitación continua reactor), en el
que se supone que los contenidos estar perfectamente mezclado Por lo tanto, la composición y la
temperatura. se supone que son uniformes en todo el reactor volumen e igual a la composición y
temperatura de el efluente del reactor. Sin embargo, los elementos fluidos no todos tener el mismo
tiempo de residencia en el reactor. Más bien, hay Una distribución del tiempo de residencia. No es
difícil proporcionar mezcla perfecta del contenido de fluido de un recipiente para aproximarse Un
modelo CSTR en un reactor comercial. Una mezcla perfecta El reactor se utiliza a menudo para
reacciones homogéneas en fase líquida. El modelo CSTR es adecuado para este caso, siempre que
que la reacción tiene lugar bajo adiabático o isotérmico condiciones Aunque los cálculos solo
involucran algebraicos ecuaciones, pueden ser no lineales. En consecuencia, un posible complicación
que debe considerarse es la existencia de múltiples soluciones, dos o más de las cuales pueden ser
estables

PFR

Más complejo es el reactor tubular de flujo de enchufe (PFR o PFTR), en el que la composición del fl
uido fluye como un enchufe, cambia gradualmente la longitud del reactor, con sin gradientes de
composición o temperatura en la dirección radial. Además, las tasas de transferencia de masa y calor
son insignificantes en la dirección axial. Por lo tanto, el PFR es completamente sin mezclar, con todos
los elementos fluidos que tienen la misma residencia tiempo en el reactor. Si el reactor funciona
bajo adiabático o condiciones no isotérmicas, la temperatura del flujo el fluido cambia gradualmente
a lo largo del reactor. Todos los simuladores proporcionan modelos unidimensionales de flujo de
enchufe que descuidan la dispersión axial. Por lo tanto, no hay radiales gradientes de temperatura,
composición o presión; y masa la difusión y la conducción de calor no ocurren en el eje dirección. El
funcionamiento del reactor puede ser adiabático, isotérmico, o no adiabático, no isotérmico. Para
este último, calor la transferencia hacia o desde la mezcla de reacción ocurre a lo largo del longitud
del reactor. Considere el caso de la operación adiabática con un químico reacción. Un balance molar
para el reactivo limitante, A, puede ser escrito como:

donde FA0 es la velocidad de flujo molar de A que ingresa al reactor, X es la conversión fraccional de
A, V es el volumen del reactor y rA es la velocidad de reacción de A escrita en función de fraccional
conversión y temperatura. Porque los simuladores de procesos calcular entalpías referidas a los
elementos, según lo descrito por Felder y Rousseau (2000), con valores para el estándar entalpía de
formación integrada en las propiedades del componente banco de datos, el calor de reacción se
maneja automáticamente, y El balance energético para la operación adiabática se convierte
simplemente en:

donde H es la tasa de flujo de entalpía de la mezcla de reacción en energía / unidad de tiempo, y T0

es la temperatura de entrada. Ecuaciones (7.24) y (7.25) se resuelven mediante integración
numérica. El siguiente ejemplo ilustra el uso del simulador. modelos para un PFR para dimensionar
un reactor adiabático de flujo de enchufe para el hidrodealquilación no catalítica de tolueno.7

A altas velocidades de flujo (números altos de Reynolds) en un tubular largo reactor, generalmente
se supone que el modelo PFR es válido porque el flujo turbulento puede aproximarse al flujo del
enchufe sin apreciable Masa axial y transferencia de calor. Sin embargo, cuando el transporte radial
los efectos son significativos, el uso de la suposición de flujo de enchufe puede conducir a errores
significativos, basados en el trabajo de Churchill y compañeros de trabajo, como se discute a
continuación en la Sección 7.5. En Reynolds números por debajo de 2,100, el flujo laminar persiste y
el modelo PFR no es válido debido al perfil de velocidad parabólica (sin enchufe). Existe una
condición parcialmente mixta con una distribución de tiempo de residencia para elementos fluidos.
Modelos más rigurosos que explican la mecánica de fluidos y el transporte se analizan en la Sección
7.5. Dinámica de fluidos computacional (CFD), como se discutió en la Sección 7S.6, puede usarse
para resolver tales modelos.

Para reacciones en fase líquida, un solo tipo PFR o CSTR El reactor es de uso frecuente. Para una
reacción única en condiciones isotérmicas, El volumen de un PFR es menor que el de un CSTR para
misma conversión y temperatura. Sin embargo, para (1) autocatalítico reacciones, donde la
velocidad de reacción depende de la concentración de un producto, o (2) reacciones autotérmicas,
donde el alimento está frío, pero la reacción es altamente exotérmica, el El volumen de un CSTR
puede ser menor que el de un PFR, de modo que La dispersión axial en un reactor tubular puede ser
realmente beneficiosa. En general, un CSTR no se usa para una reacción en fase gaseosa porque de
la dificultad para obtener una mezcla perfecta en el gas. Para reacciones homogéneas no catalíticas,
un reactor tubular es ampliamente utilizado porque puede manejar líquidos o vapores, con o sin
cambio de fase ten el reactor. El modelo PFR generalmente es adecuado para el reactor tubular si el
flujo es turbulento y si se puede suponer que cuando un cambio de fase ocurre en el reactor, la
reacción tiene lugar predominantemente en una de las dos fases Los modos térmicos más simples
son isotérmica y adiabática. El no adiabático, no isotérmico el modo generalmente es manejado por
un perfil de temperatura especificado o por transferencia de calor hacia o desde alguna fuente de
calor especificada o sumidero y un área correspondiente de transferencia de calor y transferencia de
calor general coeficiente. Ya sea una conversión fraccional de un límite reactivo o se especifica un
volumen del reactor. Los cálculos requieren la solución de ecuaciones diferenciales ordinarias. Para
reactores catalíticos de lecho fijo, un modelo PFR con un la ecuación cinética pseudo-homogénea
suele ser adecuada y se conoce como un modelo 1-D (unidimensional). Sin embargo, si el reactor no
es adiabático con transferencia de calor ao desde la pared, el modelo PFR no suele ser adecuado y
un 2- Modelo D, que involucra la solución de ecuaciones diferenciales parciales para variaciones de
temperatura y composición en ambos las direcciones axiales y radiales son necesarias. Los
simuladores no incluyen modelos 2-D, pero pueden ser generados por el usuario y insertado en el
simulador. Los modelos para reactores catalíticos de lecho fluidizado son los más complejo y no se
puede modelar adecuadamente con el Modelos CSTR o PFR. Porque parte del gas que pasa a través
del lecho fluidizado puede pasar por alto el catalizador suspendido, la conversión en un lecho
fluidizado puede ser menor que eso predicho por el modelo CSTR.

Diseño de reactores de configuración compleja:

Como se discutió en la Sección 6.5, el control de temperatura es un importante consideración en el

diseño del reactor. La operación adiabática es siempre considerado en primer lugar porque
proporciona el más simple y reactor menos costoso. Sin embargo, cuando las reacciones son
altamente exotérmica o endotérmica, a menudo es deseable hacer ejercicio algo de control sobre la
temperatura. Métodos para lograr esto, como se muestra en la Figura 7.1, incluye: (a) transferencia
de calor ao desde el fluido que reacciona, a través de una pared, hacia o desde un exterior agente de
enfriamiento o calentamiento; (b) un portador de calor inerte o reactivo o diluyente en el fluido de
reacción; (c) una serie de lechos de reactores con un intercambiador de calor para enfriar o calentar
entre cada par de camas y (d) enfriamiento por inyección en frío (también llamado contacto directo
enfriamiento rápido), donde la alimentación combinada se divide en dos o más partes, una de las
cuales ingresa al reactor entrada mientras las partes restantes entran al reactor en otra ubicaciones.
Los siguientes son ejemplos industriales de la aplicación. de estos cuatro métodos. Al considerar
estos ejemplos, un medida útil del grado de exotermicidad o endotermicidad de una reacción es el
aumento de temperatura adiabático (ATR) para reacción completa con reactivos en la relación
estequiométrica. Un ejemplo del uso industrial de un reactor intercambiador de calor. en la Figura
7.1a está en la fabricación de anhídrido ftálico, producido por la oxidación de ortoxileno con aire en
el presencia de partículas de catalizador de pentóxido de vanadio, como se discutió por Rase (1977).
La reacción, que se lleva a cabo aproximadamente a las 375 C y 1.2 atm, es altamente exotérmico
con un ATR de aproximadamente 1; 170 C, incluso con nitrógeno en el aire proporcionando cierta
dilución. La operación adiabática no es factible. El reactor se asemeja a unintercambiador de calor
vertical de carcasa y tubo. Cientos de tubos largos de pequeño diámetro, dentro de la carcasa, están
llenos de catalizadorpartículas y a través de las cuales pasa el gas de reacción hacia abajo.Un medio
de transferencia de calor que consiste en un nitrito de sodio-potasio la sal fusionada con nitrato
circula fuera de los tubos a través del cáscara para eliminar el calor de reacción. El agua se descarta
como calor. transferir medio en este caso porque la presión de agua requerida Sería muy alto. La
distribución de la tasa de transferencia de calor no es adecuado para mantener condiciones
isotérmicas para la reacción fluido, pero su temperatura cambia en menos de 40 C. En algunos
aplicaciones, el arreglo involucra camas llenas de catalizador intercalados con tubos que transportan
un refrigerante (por ejemplo, el TVA diseño en la Figura 6.11). El estireno es producido por la
deshidrogenación catalítica de etilbenceno a 1,2 atm y una temperatura de aproximadamente 575 C,
como descrito por Smith (1981). La reacción es suficiente. endotérmica, con un ATR de
aproximadamente 460 C, de modo que si el reactor se hizo funcionar adiabáticamente con una
alimentación de puro etilbenceno, la temperatura del fluido de reacción sería disminuir hasta tal
punto que la velocidad de reacción sería indebidamente comprometido, lo que resulta en un
volumen de reactor muy grande. Para mantener una temperatura razonable, una gran cantidad de
vapor se agrega a la alimentación (típicamente con una relación molar de vapor a etilbenceno igual a
20: 1), que se precalienta a 625 C antes de entrar al reactor (Figura 7.1b). El vapor es inerte y se
recupera fácilmente del efluente del reactor por condensación. La presencia del vapor reduce la
reacción. tasa porque la concentración de estireno se reduce, pero el El reactor puede funcionar
adiabáticamente de una manera simple. El trióxido de azufre, que se usa para producir ácido
sulfúrico, es producido por oxidación catalítica de dióxido de azufre en el aire con partículas de
catalizador de pentóxido de vanadio a 1,2 atm y un temperatura de aproximadamente 450 ° C,
según lo discutido por Rase (1977). La operación adiabática no es factible ya que la reacción es
altamente exotérmico con un ATR de aproximadamente 710 C, incluso con nitrógeno en el aire
proporcionando algo de dilución. Por lo tanto, el sistema del reactor consta de cuatro lechos de
reactores adiabáticos, del mismo diámetro pero diferente altura, en serie, con un intercambiador de
calor entrecada par de camas, como se muestra en la Figura 7.1c. La temperatura se eleva
adiabáticamente en cada lecho del reactor, y el reactor caliente el efluente se enfría en el
intercambiador de calor colocado antes del cama de al lado. Cuando el ATR es mayor, como en la
fabricación de amoníaco del gas de síntesis, según lo descrito por Rase (1977), Se recomienda el
diseño de inyección en frío en la Figura 7.1d. Para reactores catalíticos de lecho fijo 1-D, es deseable
reduzca el volumen del recipiente al mínimo. Según lo presentado por Aris (1960), esto se logra por
diseño. Si z es la dirección a lo largo del reactor, la trayectoria de la masa y ecuaciones de balance de
energía para una reacción única en X (z) yT (z) el espacio se ajusta para que coincida con la
trayectoria correspondiente a la velocidad de reacción máxima, (X *, T *) [línea curva en la Figura
7.2], lo más cerca posible. Por lo tanto, los reactores refrigerados por tubo (o calentados),
convertidores de inyección en frío (o inyección en caliente) y múltiples camas adiabáticas con
intercoolers (o intercalentadores) deben diseñarse cuidadosamente de una manera que

Como ejemplo, considere una reacción reversible exotérmica en un PFR. Para este caso, la velocidad
de la reacción inversa aumenta más rápidamente al aumentar la temperatura que la tasa de La
reacción hacia adelante. Además, la reacción inversa es lenta y el Reacciona rápidamente a bajas
temperaturas. Por lo tanto, para un máximo velocidad de reacción, la temperatura debe ser alta a
bajas conversiones y bajo a altas conversiones. Esto se muestra en la Figura 7.2, tomado de Smith
(1981), donde la velocidad de reacción para una secuencia de conversiones fraccionarias, X,
comenzando con X1 0, se traza contra temperatura, T. Para cada valor de X, la velocidad de reacción
La curva de la figura 7.2 muestra un valor máximo. Un lugar de se muestran las tasas máximas,
correspondientes al sólido y discontinuo línea que pasa por los puntos C y B, con el máximo la
velocidad de reacción disminuye al aumentar la conversión fraccional. En cada nivel de conversión,
la temperatura deseada corresponde a la velocidad de reacción máxima. En la Figura 7.2, la
alimentación entra en temperatura TA, con una velocidad de reacción en el punto A. Aunque el no
se muestra la velocidad máxima de reacción para X1, está claro que TA es no la temperatura
correspondiente a la tasa máxima. Si el la temperatura de entrada no se puede aumentar, es mejor
operar isotérmicamente en TA hasta que se alcanza la conversión en el punto C, y luego siga el perfil
óptimo CB hasta la conversión deseada. Supongamos que la conversión de salida del reactor es X4.
Entonces lo deseado la trayectoria de temperatura del reactor es la línea continua ACB, con reactor
temperatura de salida TB. Correspondiente a esta trayectoria, pero no que se muestra en la Figura
7.2, es un perfil de servicio térmico, que debe ser combinado con intercambio de calor para lograr la
velocidad de reacción óptima trayectoria. Alternativamente, la operación isotérmica del reactor a TA
corresponde a la trayectoria ACD en la Figura 7.2. En este caso, las velocidades de reacción no están
en sus valores máximos, excepto en punto C, que requiere un mayor volumen de reactor. Si en lugar
de un PFR, un Se utilizaron CSTR, la temperatura óptima de operación para lograr la conversión X4
sería TB, que corresponde a la velocidad máxima de reacción para esa conversión. Al especificar un
perfil de temperatura para un PFR o una temperatura de salida para un CSTR, El volumen óptimo del
reactor se puede determinar junto con el perfil de servicio térmico correspondiente requerido.
Según lo discutido por Van Heerden (1953), la alimentación del reactor la temperatura tiene un
efecto importante en la estabilidad de un reactor autotérmico, es decir, un reactor cuya
alimentación está precalentada por su efluente. Para una reacción exotérmica reversible, como en
síntesis de amoníaco, la tasa de generación de calor varía de forma no lineal con la temperatura de
reacción, como se muestra en la curva (a) en la Figura 7.3. A bajas temperaturas, la tasa de calor la
generación está limitada por la baja tasa de reacción directa al amoniaco A temperaturas muy altas,
la velocidad de reacción es limitado por el equilibrio, de modo que, nuevamente, baja generación de
calor Las tasas son de esperar. Sin embargo, en algún intermedio temperatura, la velocidad de
reacción exhibe un valor máximo. En contraste, debido a que la transferencia de calor por
convección es dominante, el La velocidad de eliminación del calor es casi lineal con la reacción
temperatura, con una pendiente que depende del grado de calor intercambio entre la salida y la
entrada. Por lo tanto, la intersección de la línea de eliminación de calor (b) y la generación de calor la
línea (a) a veces conduce a tres posibles condiciones de operación: (O) el estado de no reacción, (I) el
punto de ignición y (S) El punto de operación deseado. Tanto el que no reacciona como el los puntos
de operación deseados son estables, ya que un pequeño positivo la perturbación en la temperatura
del reactor provoca el calor tasa de eliminación para superar la tasa de generación de calor,
disminuyendo La temperatura del reactor. Del mismo modo, una pequeña perturbación negativa en
la temperatura del reactor tiene el efecto contrario, conduciendo a un aumento de la temperatura.
Usando los mismos argumentos, el el punto de ignición es inestable debido a una pequeña
perturbación positiva en temperatura conduce a un salto al deseado, estable punto de operación,
mientras que una perturbación negativa conduce a un llamado "reventón" al estado estable y sin
reacción. Nota que un análisis similar para un CSTR adiabático en un estudio de caso 12S.1 también
detecta la posibilidad de tres estados estacionarios. Para convertidores refrigerados por tubo, Van
Heerden (1953), y para convertidores de disparo en frío, Stephens y Richards (1973), se refieren a La
diferencia de temperatura entre los puntos de operación I y S como el "margen de estabilidad".
Claramente, la operación en S con mayor los márgenes de estabilidad serían más robustos a las
perturbaciones. Así, un diseño con una mayor tasa de transferencia de calor, indicado por la línea
(b0) en la Figura 7.3, claramente tendría un margen de estabilidad más bajo. En tales casos, la
proximidad del punto de operación estable al punto de ignición inestable conduce a una mayor
probabilidad de pérdida de control frente a los trastornos del proceso. Por ejemplo, el catalizador de
síntesis de amoníaco se desactiva, principalmente por envenenamiento debido a impurezas de
alimentación, o por sinterización a alta temperatura, lo que reduce el área superficial del catalizador.
Lewin y Lavie (1984) estudiaron el efecto de la desactivación del catalizador en el funcionamiento
óptimo de un convertidor de amoníaco enfriado por tubo, lo que puede conducir a una pérdida de
estabilidad, ya que la disminución de la actividad del catalizador conduce a tasas de generación de
calor más bajas, como se muestra en la línea (a0) en la figura 7.3. El siguiente ejemplo ilustra cómo
se diseña un reactor de inyección en frío para maximizar la conversión mientras se satisfacen los
márgenes de estabilidad.

DISEñO DE REDES DE REACTORES

Esta sección describe el uso de la región alcanzable (AR), que define las composiciones alcanzables
que pueden ser obtenido de una red de reactores químicos. Esto es análogo al tema de las
composiciones de productos factibles en destilación, presentada en la Sección 8.5. La región
alcanzable en el espacio de composición fue introducido por Horn (1964), con desarrollos y
extensiones más recientes de Glasser y compañeros de trabajo (Glasser et al. 1987; Hildebrandt et
al., 1990). La figura 7.7 ilustra la región alcanzable para van de Vusse cinética (van de Vusse, 1964),
basada en las reacciones:

Las reacciones 1, 2 y 3 son de primer orden en A, B y B, respectivamente, mientras que la reacción 4

es de segundo orden en A. La velocidad las constantes a una temperatura particular son: k1=0:01 s-
1,k2 =5s-1, k3=10 s-1 y k4=100 m3/ kmol s.

El límite de la región alcanzable, que se muestra en la Figura 7.7, es compuesto de arcos, cada uno
de los cuales resulta de la aplicación de un tipo de reactor distinto, como se describe a continuación.
Para el caso de la cinética de van de Vusse con una alimentación de 1 kmol / metro3 de A, la Figura
7.7 indica que el límite AR es com planteado de un arco que representa un CSTR con derivación
(curva C), un CSTR (punto O), y un CSTR seguido de un PFR (curva D).

Dentro de la región delimitada por los tres arcos y la horizontal línea base (CB=0) , se pueden lograr
composiciones de productos con alguna combinación de estas configuraciones del reactor. Los
configuración adecuada del reactor a lo largo del límite de la región alcanzable depende de la
concentración deseada de fluyente de A. Cuando 1> CA> 0:38 kmol / m3, un CSTR con bypass (curva
C) proporciona la concentración máxima de B, mientras que cuando CA <0:38 kmol / m3 , esto se
logra utilizando un CSTR (punto O), seguido de un PFR (Curva D). Tenga en cuenta que el
concentración máxima alcanzable, CB=1:25x10-4kmol / m3 , se obtiene usando un CSTR seguido de
un PFR (en el punto M a lo largo de la curva D). Evidentemente, la región alcanzable proporciona
asistencia útil en el diseño del reactor óptimo redes. Un procedimiento para la construcción de
alcanzables regiones se discuten a continuación.

Construcción de la región alcanzable:

Un método sistemático para la construcción de la región alcanzable usando CSTR y PFR, con o sin
mezcla y de derivación, para un sistema de reacciones químicas, tal como se presenta por
Hildebrandt y Biegler (1995), se demuestra por van de cinética Vusse

1. Comience por la construcción de una trayectoria para un PFR desde el punto de

alimentación, continuando a la conver- sión completa de A o equilibrio químico. En este
caso, el PFR trayectoria se calcula mediante la resolución de simultánea- mente las
ecuaciones cinéticas para A y B:

donde t es el tiempo de residencia PFR. Tenga en cuenta que no se requieren ecuaciones

cinéticas para C y D para la cons- trucción de la región alcanzable en el espacio de dos
dimensiones debido a que sus composiciones no aparecen en las ecuaciones. (7.40) y (7.41).
La trayectoria en CUN Csi espacio se representa en la figura 7.8a como curva ABC. En este
ejemplo, un componente se convierte completamente.
2. Cuando la trayectoria de PFR limita una región convexa, esto constituye una región
alcanzable candidata. Un la región convexa es aquella en la que todas las líneas rectas
dibujado desde un punto en el límite a cualquier otro punto en el límite se encuentran
totalmente dentro de la región o en el límite Si no, la región no es convexa. cuando los
vectores de velocidad, en concentraciones fuera del candidato AR no apuntan de nuevo en
él, los límites actuales son el límite de n AR y el procedimiento termina. En este ejemplo mas
alla de la curva
3. La trayectoria PFR se expande por arcos lineales, lo que representa la mezcla entre el PFR
efluente y la corriente de alimentación, que se extiende la región alcanzable candidato.
Tenga en cuenta que un arco lineal que conecta dos puntos en una trayectoria composición
se expresa por la ecuación:

donde C 1 y C 2 son vectores de dos corrientes en el espacio de la composición, c * es la

composición de la corriente mixta, y a es la fracción de la corriente con la composición c 1 en
la corriente mixta. Los arcos lineales se ensayaron a continuación para asegurar que no hay
vectores de tasa posicionado en ellos señalan de la AR. Si hay
 tales vectores, proceder al siguiente paso, o volver al paso 2. Como se muestra en la figura
7.8a, un lineal de arco, ADB, se añade, que se extiende la región alcanzable a ADBC. Puesto
que para este ejemplo, los vectores de tasas calculadas a lo largo de este arco se encuentran
señalar de la AR extendida, proceder al siguiente paso.

4. Puesto que hay vectores que apuntan fuera de la envolvente convexa, formados por la unión
entre la trayectoria PFR y lineal arcos de mezcla, es posible que una trayectoria CSTR
agranda la región alcanzable. Después de colocar la trayectoria CSTR que se extiende la AR
se colocan los mayoría, arcos lineales adicionales que repre- envió la mezcla de las
corrientes para asegurar que los restos AR convexo La trayectoria CSTR es calculado
resolviendo la forma CSTR de la cinética ecuaciones para A y B, dadas por las ecuaciones.
(7.40) y (7.41) en función del tiempo de residencia, t

Para este ejemplo, la trayectoria CSTR que se extiende la AR más es la calculada desde el
punto de alimentación, en CA0, la mayor concentración de A. Esto se indic como curva AEF
en la Figura 7.8b, que pasa a través de punto B. Desde la unión de la AR previa y la La
trayectoria CSTR no es convexa, un arco lineal, AGO, es aumentado como se muestra en la
figura 7.8c. Este arco representa un CSTR con una secuencia de derivación.
5. Se dibuja una trayectoria PFR desde la posición donde la línea de mezcla se encuentra con la
trayectoria CSTR. Si esto La trayectoria del PFR es convexa, extiende la AR previa para formar
un candidato ampliado AR. Luego regrese a Paso 2. De lo contrario, repita el procedimiento
desde el Paso 3. Como se muestra en la Figura 7.8d, la trayectoria del PFR, OHI, conduce a
una región convexa alcanzable. Las fronteras de la región son: (a) el arco lineal, AGO, que
representa un CSTR con flujo de derivación; (b) el punto O, que representa un CSTR; y (c) el
arco OHI, que representa un CSTR seguido de un PFR en serie. Se observa que la
composición máxima de B es obtenido en el punto H, usando un CSTR seguido de un PFR.

Arly, el diseño óptimo del reactor minimiza el anualizado costo, calculado para contabilizar el
capital y operativo costos, y no simplemente el diseño que maximiza el rendimiento
selectividad. No obstante, la región máxima alcanzable identifica todo el espacio de
concentraciones ofensibles. El seguimiento El ejemplo muestra cómo se utiliza la región
alcanzable para seleccionar red de reactores más apropiada para maximizar el rendimiento de u
producto deseado donde se producen varias reacciones competitivas

El principio de las invariantes de reacción

Debido a que la región alcanzable depende de construcciones geométricas, se limita efectivamente

al análisis de sistemas involucrando dos especies independientes. Sin embargo, como lo muestra
Omtveit y col. (1994), sistemas que implican dimensiones superiores puede ser analizado usando el
enfoque AR bidimensional por aplicando el principio de reacción invariantes de Fjeld et al. (1974) La
idea básica consiste en imponer equilibrios atómicos. sobre las especies que reaccionan. Estas
restricciones lineales adicionales imponer una relación entre las especies que reaccionan, permitir-
Hacer que el sistema completo se proyecte en un espacio reducido de especies independientes. El
análisis AR anterior puede usars cuando este espacio reducido es en dos dimensiones. Deje que el
sistema de reacción consista en ni moles de cada especie. i, cada uno contiene aij átomos del
elemento j. Los cambios molares en cada una de las especies debido a la reacción se combinan en el
vector Dn, y los coeficientes aij forman la matriz atómica Anotando eso desde el número de átomos
de gramo para cada elemento permanecer constante, AH=0 Particionando Dn y A en dependiente,d,
e independientes, i, componentes:

Suponiendo que A re es cuadrado y no singular, una expresión para los cambios en el número de
moles de cada dependient la especie se obtiene por manipulación algebraica:

La dimensión de i es igual al número de especies menos El número de elementos (átomos) en la

especie. Cuando esto dimensión es dos o menos, el principio de invariantes de reacción permite
la aplicación de la región alcanzable a complejos sistemas de reacción Esto se ilustra en el
siguiente ejemplo, introducido por Omtveit et al. (1994)

MODELOS RIGOROSOS PARA REACTORES QUIMICOS TUBULARES

En las secciones anteriores de este capítulo, modelos de reactores ideales (PFR y CSTR) fueron
desarrollados y aplicados a simples configuraciones complejas y a redes. Estos modelos ideales
que representan los dos extremos de la no mezcla axial y mezcla completa, están ampliamente
disponibles en computadora programas de simulación, y solo requieren termodinámica y datos
de cinética de reacción. Reactores continuos de tanque agitado, que se usan a menudo para
reacciones en fase líquida, pueden ser con fi simulado con el modelo CSTR ideal si es adecuado
Se proporciona agitación del contenido del recipiente. Sin embargo, para la aplicación a
reactores de flujo tubular, el PFR ideal modelo, que ignora el transporte axial de calor y masa y
no asume gradientes de velocidad, temperatura o composición en la dirección radial, a veces
puede, como se muestra en Esta sección, esta en grave error. Con el advenimiento de fácilmente
métodos numéricos aplicados y / o el uso de CFD (computacional programas de dinámica de
fluidos), ahora es posible, como discutido en esta sección, para dar cuenta rutinariamente de
radial gradientes en reactores de flujo tubular. Esta sección trata sobre una introducción a
soluciones rigurosas para reactores de flujo tubular, particularmente aquellos que operan bajo
condiciones de flujo laminar. En la sección final de este capítulo, Sección 7S.6, un El diseño
detallado del reactor se presenta utilizando un modelo CFD con El programa COMSOL.

Condiciones isotermales

El caso más simple de un reactor de flujo tubular isotérmico que no se ajusta al modelo PFR es
un reactor de flujo laminar presentado por Cleland y Wilhelm (1956). Asumieron una reacción
química irreversible de primer orden que ocurre bajo isotérmica, continua, en estado estable,
densidad constante (no cambio de volumen en la reacción), con flujo laminar completamente
desarrollado condiciones en un tubo recto y circular de diámetro constante. Descuidaron los
efectos de convección (naturales) libres y difusión axial de las especies, pero tuvo en cuenta la
parabólica distribución de velocidad y difusión radial de especies. La distribución de velocidad
radial causa una distribución radial de tiempos de residencia, lo que a su vez provoca una
concentración radial degradado. Sin embargo, el gradiente de concentración disminuye por
difusión molecular radial. Por lo tanto, si la difusión radial es rápido, se aborda el supuesto de fl
ujo del enchufe; sin embargo, si es lento, el efecto sobre la longitud del reactor requerido,
comparado para eso para el fl ujo del enchufe, puede ser muy significativo, como será mostrado.
Bajo los supuestos de Cleland y Wilhelm, el ecuación de gobierno, una ecuación diferencial
parcial parabólica, es un único balance de masa diferencial, en un punto del reactor para las
especies que controlan la velocidad de reacción. Así,

donde vavg es la velocidad de flujo promedio, r es la distancia radial desde el eje del tubo, R es el
radio interior del tubo, ci es el concentración molar de la especie que controla la velocidad, z es
el distancia axial desde la entrada del reactor, Di es el molecular difusividad de las especies que
controlan la velocidad, yk es el primero velocidad de reacción constante. Si la concentración de
entrada es ci0, las condiciones de contorno para la ecuación. Eq. (7.62) no puede resolverse
analíticamente, excepto en los dos límites de la relación de difusividad molecular a reaccionar
constante, expresada como Di / kR. Cuando un se acerca infinito (difusión rápida en
comparación con la velocidad de reacción y / o tubo de diámetro pequeño), se aborda el
supuesto de flujo de enchufe porque la concentración es uniforme sobre la cruz sección del tubo
en cada posición axial. Entonces, la ecuación. (7.62) con sus condiciones límite se pueden
integrar para dar,

Cuando a se aproxima a cero (difusión lenta en comparación con la velocidad de reacción y / o

un tubo de gran diámetro), se forma un gradiente de concentración radial debido a la
distribución de los tiempos de residencia de las líneas de flujo de flujo laminar sobre la sección
transversal del tubo. Integración de la ecuación. (7.62) da un resultado más complicado, donde
la concentración depende de r y z:

donde la integral dentro de los corchetes es la exponencial integral del orden 3, E, donde Eq.
(7.64) no es aplicable a condiciones de flujo turbulento porque el perfil de velocidad no es
parabólico. Para ese caso un modelo matemático dado por Churchill y Yu (2006), que utiliza una
expresión muy precisa para la velocidad radial distribución en flujo turbulento (Churchill, 2001;
Churchill y Zajic, 2002), puede aplicarse o puede emplearse CFD. Usando las ecuaciones. (7.63) y
(7.64), los dos casos limitantes para el flujo laminar se compara en la figura 7.13, donde el i el
subíndice ha sido eliminado. A primera vista, los dos límites Los resultados no parecen ser muy
diferentes. Sin embargo, de hecho, como la conversión fraccional ci / ci0 se acerca a 1 (es decir,
como c / c0, la fracción no convertida, se acerca a 0 en la Figura 7.13), la diferencia entre los dos
resultados se vuelve muy significativo cuando el tiempo de residencia promedio requerido o su

Se considera la equivalencia, longitud del reactor. Esto es se muestra en la figura 7.14, donde la
relación de la longitud del reactor para longitud de difusión radial a reactor cero para radial
infinito la difusión (flujo de tapón) se representa contra la fracción sin reaccionar en una escala
logarítmica. Se incluyen los resultados de Cleland y Wilhelm para una reacción de primer orden ,
Juntos con resultados de Bosworth (1948) para orden cero, Denbigh (1951) para reacciones de
segundo orden Considerar Una conversión fraccional de 0.99 (es decir, 0.01 o 1% sin reaccionar).
La figura 7.14 muestra las relaciones limitantes del reactor. longitudes de 1.82, 1.53 y 1.34 para
cero, primero y segundo orden reacciones, respectivamente. Por lo tanto, el efecto de la falta de
La difusión radial en el flujo laminar puede ser sustancial. Los reactores de flujo laminar reales
operarán entre los dos extremos de la figura 7.13. Para condiciones particulares entre Los
extremos, soluciones numéricas de la ecuación. (7.62), que son obtenido fácilmente por
diferencia finita estándar, elemento finito, Colocación, o procedimientos CFD para resolver
diferencial parcial ecuaciones, pueden ser empleados. Cleland y Wilhelm usaron un método
implícito de diferencia finita, que se discute en detalle de Lapidus (1962).

Condiciones No isotermicas:
Cuando no se puede suponer el flujo isotérmico, la velocidad de reacción constante en la
ecuación. (7.62) no será constante, pero dependerá de temperatura, de acuerdo con la ecuación
de Arrhenius

donde k1 es el factor pre-exponencial, E es la energía de activación por mol para la reacción, R es

la constante de gas universal y T es la temperatura absoluta. Además, un balance de energía
térmica debe estar acoplado al balance diferencial de masa de especies, ec. (7.62). Si asumimos
constante propiedades, el balance de energía diferencial es

donde r es la densidad, Cp es el calor específico, kT es el calor conductividad, y DHrx es el calor

de reacción por mol de especie reaccioné, (positivo) por una reacción endotérmica y (negativo)
para una reacción exotérmica Las condiciones de contorno para la ecuación. (7.66) incluye T=T0
(la temperatura de entrada del fluido) az=0andqT / qr= 0 en la línea central del tubo donde r=0

La condición límite requerida adicional generalmente se impone en la pared interior del tubo, r,
Qw, con una condición adiabática de Qw= 0, donde Qw=kTre qT / qr es el flujo de calor en la
pared interior del tubo.nConsidere el caso más simple de un reactor tubular adiabáti con en
condiciones extremas de difusión radial rápida en comparación con la velocidad de reacción y / o
un diámetro de tubo muy pequeño. Esto es equivalente a la suposición de flujo de enchufe. Los
gradientes radiales en las ecuaciones (7.62) y (7.66) desaparecen para que estas dos ecuaciones
convertirse en ecuaciones diferenciales ordinarias que pueden ser combinado eliminando el
diferencial dz, resultando en La siguiente relación entre la concentración de especies y
temperatura, que es aplicable tanto a la lámina tubular como a la condiciones de flujo
turbulento.

Ahora, suponiendo una reacción irreversible de primer orden,

Eq. (7.67) se puede combinar con (7.68) para eliminar cualquier ci o T, de modo que la ecuación
resultante se pueda integrar. Como mostrado por Churchill (1974, 2005), es más conveniente
eliminar T e introducir nuevas variables:
y temperatura de salida del reactor adiabático, T, para completar conversión re

donde se supone que el término entre paréntesis es independiente de temperatura La ecuación

diferencial ordinaria única resultante es luego integrado por partes para dar una solución de
forma cerrada en términos de funciones integrales exponenciales de los siguientes formar,

donde y es una variable ficticia y Ei es diferente de función integral exponencial en la ecuación.

(7.64), pero es otro en un serie de funciones integrales exponenciales. Tanto exponencial las
funciones integrales utilizadas en esta sección están tabuladas por Abramowitz y Stegun (1964) y
se pueden calcular con programas de software como Maple, Mathematica, Matlab y Excel con
un complemento apropiado. La solución de forma cerrada para el reactor de flujo de enchufe
adiabático, tubular es

donde Y Para una conversión fraccionaria de salida especificada, se puede utilizar la ecuación
(7.72) para calcular la longitud del reactor y la ecuación. (7.67) puede dar la salida temperatura
absoluta. Soluciones de forma cerrada para irreversible Reacciones de otros órdenes
inadiabáticos, reactores de flujo enchufables por Douglas y Eagleton (1962). Estas soluciones
también implican funciones integrales exponenciales. Cuando la suposición de rápido la difusión
radial no es válida, las ecuaciones de transporte (7.62) y (7.66) debe resolverse simultáneamente
por un número apropiado métodos o con un programa CFD, como se ilustra en la Sección 7S.6.