Versión 2.0 · reversionada
El reactor batch, en tiempo real.
Gemelo digital de un Tanque Agitado Discontinuo: balances de materia y energía, equilibrio químico y dimensionamiento, con un motor nuevo que corre en tu navegador. Construido sobre las bases de la tesis «Modelado y Simulación de un Reactor Batch en Python».
Reacción exotérmica adiabática calculada en vivo por el motor 2.0
Qué hay de nuevo
Resultados al instante
El motor corre en tu navegador. Sin servidor que despertar: cada cambio se recalcula en milisegundos.
Gemelo digital interactivo
Mové la temperatura, el caudal de refrigerante o la cinética y mirá cómo responde el reactor en vivo, con reproducción del batch.
Cinética generalizada
N especies, coeficientes y órdenes arbitrarios (ley de potencias), reactivos en exceso, inertes y elección de la especie clave.
Nuevas herramientas de análisis
Diagrama conversión–temperatura con equilibrio, recta adiabática y curva de T óptima; gráfico de Levenspiel; calor generado y removido.
Calibrar con datos
Cargá una corrida medida y ajustá k, E, ΔH y UA con Levenberg–Marquardt, con su incertidumbre y la comparación contra el modelo nominal.
Mapa de riesgo
Una red neuronal entrenada con 24 000 simulaciones dibuja al instante la frontera de runaway según T₀, refrigerante y UA.
Pronóstico en vivo
Una corrida llega dato a dato y dos métodos, ML y física + datos, pronostican la temperatura y avisan antes del runaway.
Validado contra la tesis
Pruebas automáticas contra soluciones analíticas y contra el motor Python original, con las diferencias documentadas.
Integradores propios
Dormand–Prince 5(4) con salida densa y eventos, y Rosenbrock para casos rígidos como el descontrol térmico.
De la tesis a la 2.0
La versión de la tesis se conserva sin cambios como trabajo académico. La 2.0 mantiene su modelo (balances, Arrhenius, Van't Hoff, Gibbs, intercambiador ICQ y dimensionamiento) y corrige por diseño los puntos detectados al revisarla.
| Aspecto | Tesis (v1) | Versión 2.0 |
|---|---|---|
| Velocidad inversa | Usaba C_C·C_D con el signo estequiométrico cambiado; para A ⇌ B valía cero. | −r = k·(Π Cᵢ^αᵢ − Π Cⱼ^βⱼ / K_c), coherente con K_eq. |
| Equilibrio con intercambiador (ICQ) | Al llegar a X_eq se congelaban X y también T. | El equilibrio surge de la cinética; el intercambiador sigue actuando. |
| Volumen en el balance de energía | Implícito en 1 L. | Volumen de mezcla V explícito. |
| Corte de la integración | Adiabático hasta t = 50 000 con paso 1 (segundos de cálculo). | Eventos: conversión objetivo o reacción completa (milisegundos). |
| Unidades de tiempo | La producción en kg/min se mezclaba con tiempos en segundos. | Unidad de tiempo explícita; producción en kg/h convertida. |
| Precisión numérica | Tolerancia por defecto de SciPy (rtol = 1e-3). | rtol = 1e-8, verificada contra soluciones analíticas. |
| Equilibrio con van't Hoff | ΔH por mol de A con K_c por reacción: con 2A ⇌ C, K_c(T) variaba a la mitad de la pendiente. | ΔH de la reacción escrita (a·ΔH); con a = 1 coincide con la tesis. |
Próximamente: asistente con inteligencia artificial
Explicación de resultados en lenguaje natural y configuración de casos describiéndolos en palabras. Más adelante: un modelo híbrido física + ML entrenado con datos medidos y múltiples reacciones simultáneas.