Scribd
Cargar
35 vistas2 páginas

Polymath RX

El documento describe un archivo Polymath para modelar un reactor de flujo pistón (PFR) considerando balances de masa, energía y presión utilizando un catalizador de cromia-alúmina. Se presentan datos iniciales como flujos molares, constantes cinéticas, y condiciones iniciales de temperatura y presión. Además, se incluyen ecuaciones diferenciales para el balance de masa, energía y la ecuación de Ergun para la presión.

Cargado por

tommasielloexequiel
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
35 vistas2 páginas

Polymath RX

El documento describe un archivo Polymath para modelar un reactor de flujo pistón (PFR) considerando balances de masa, energía y presión utilizando un catalizador de cromia-alúmina. Se presentan datos iniciales como flujos molares, constantes cinéticas, y condiciones iniciales de temperatura y presión. Además, se incluyen ecuaciones diferenciales para el balance de masa, energía y la ecuación de Ergun para la presión.

Cargado por

tommasielloexequiel
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

; Archivo Polymath para modelar un reactor PFR con balance de masa, energia y presion

; Catalizador Cromia-Alúmina

; ========================

; Datos iniciales

; ========================

F0 = 107.71 1000 3600 ; Flujo molar inicial (mols)

k0 = 5.9e-3 ; Constante preexponencial (s^-1)

Ea = 96.2 1000 ; Energía de activación (Jmol)

R = 8.314 ; Constante de gas (Jmol·K)

Cp = 0.03 1000 ; Capacidad calorífica (Jmol·K)

deltaH = 129 1000 ; Calor de reacción (Jmol)

T0 = 550 + 273.15 ; Temperatura inicial (K)

P0 = 3.5 1e5 ; Presión inicial (Pa)

C0 = 1.0 ; Concentración inicial de n-butano (molL)

epsilon = 0.4 ; Porosidad del lecho

dp = 0.003 ; Diámetro de partícula (m)

mu = 1.8e-5 ; Viscosidad del gas (Pa·s)

A = 0.01 ; Área de la sección transversal del reactor (m²)

; ========================

; Definiciones auxiliares

; ========================

k = k0 EXP(-Ea (R T)) ; Constante cinética (s^-1)

r=k C ; Velocidad de reacción (molL·s)

u = F0 (rho A) ; Velocidad superficial del gas (ms)

rho = P (R T) ; Densidad del gas (kgm³)

; ========================

; Ecuaciones diferenciales

; ========================
D(C)D(tau) = -r ; Balance de masa

D(T)D(tau) = (-r deltaH) (F0 Cp) ; Balance de energía

D(P)D(tau) = -(150 (1 - epsilon)^2 mu u (dp^2 epsilon^3))

- (1.75 (1 - epsilon) rho u^2 (dp epsilon^3)) ; Ecuación de Ergun

; ========================

; Condiciones iniciales

; ========================

C = C0 ; Concentración inicial (molL)

T = T0 ; Temperatura inicial (K)

P = P0 ; Presión inicial (Pa)

; ========================

; Integración

; ========================

tau = 0 ; Tiempo espacial inicial (s)

END tau = 100 ; Tiempo espacial final (s)

DT = 0.1 ; Paso de integración

También podría gustarte