COMPLEXA · Ciencias de la Complejidad y Teoría del Caos
Módulo 18. Organización fuera del equilibrio y redes de reacción

171.1. Sistema, frontera y entorno

Un sistema termodinámico queda definido al especificar qué parte del mundo se estudia, qué frontera la delimita y qué intercambios la atraviesan. Esta elección determina los balances, las variables observables y el significado de las condiciones externas. Para analizar una organización sostenida fuera del equilibrio conviene identificar primero sus reservas, sus corrientes y sus mecanismos de relajación. La sola apariencia estable de una estructura informa poco sobre los procesos que permiten conservarla.

La frontera puede ser una pared material, una membrana selectiva o una superficie imaginaria de control. El entorno comprende los cuerpos y reservorios capaces de intercambiar materia, energía o trabajo con el sistema. Su descripción depende del problema; una temperatura impuesta puede bastar para representar un baño térmico, mientras que una alimentación química requiere especificar composición, caudal y, cuando corresponda, presión y temperatura. El entorno tampoco tiene que permanecer inalterado durante toda la observación.

Un sistema abierto intercambia materia y energía con su entorno. Un reactor alimentado continuamente pertenece a esta categoría, aunque su volumen y sus concentraciones se mantengan constantes. La energía puede cruzar la frontera con la materia, como calor o como trabajo. Hay que distinguir estos canales para no contar dos veces una misma contribución y para asignar correctamente el transporte de entalpía de las corrientes.

Un sistema cerrado intercambia energía, pero no materia. Un recipiente sellado puede recibir calor y efectuar trabajo sobre el exterior sin dejar de ser cerrado. Un sistema aislado no intercambia ni materia ni energía. Estas clasificaciones describen solo las condiciones de intercambio y son independientes de la complejidad interna, la homogeneidad espacial y la rapidez de evolución. Una aproximación de aislamiento puede ser adecuada durante un intervalo breve y dejar de serlo a tiempos largos.

El no equilibrio puede existir como transitorio sin apertura. Dos regiones de un recipiente aislado, inicialmente a temperaturas distintas, intercambian energía dentro del sistema y se aproximan al equilibrio. Durante ese proceso hay gradientes y producción interna de entropía. La apertura interviene cuando se trata de sostener diferencias durante tiempos prolongados, siempre que el entorno pueda compensar su relajación.

También importa la resolución de la descripción. Una concentración media constante puede ocultar perfiles espaciales, fluctuaciones o reacciones compensadas. Al elegir variables macroscópicas se acepta una reducción en la que las propiedades microscópicas restantes quedan representadas mediante relaciones constitutivas y parámetros efectivos. La validez de esa reducción exige separar escalas y comprobar que los mecanismos omitidos no controlan la respuesta observada en el régimen estudiado.

En las secciones siguientes, un reactor agitado permitirá relacionar balances, estados estacionarios y escalas temporales. Su sencillez deja ver que sostener una corriente no garantiza una inestabilidad ni una organización compleja. Estas últimas requieren mecanismos dinámicos adicionales, restricciones y condiciones de contorno específicas, cuya presencia debe establecerse con un modelo y con observaciones.

171.2. Balances extensivos de masa y energía

Para una cantidad extensiva XX dentro de una frontera, el cambio temporal es la suma de transportes y transformaciones internas; tt denota el tiempo. Las tasas de entrada y salida se definen positivas y sus signos se explicitan en el balance. La generación es un término algebraico que puede corresponder a producción o a consumo.

dXdt=X˙ent−X˙sal+X˙gen.\frac{dX}{dt}=\dot X_{\mathrm{ent}}-\dot X_{\mathrm{sal}}+\dot X_{\mathrm{gen}}.

Los puntos indican tasas por unidad de tiempo. Una reacción química ordinaria no genera masa total, aunque cambien las cantidades de cada especie. Para una especie, el término químico depende del coeficiente estequiométrico y de la velocidad de reacción; el sumidero de un reactivo es compatible con la conservación de masa y elementos en el conjunto de reactivos y productos.

Sea nn la cantidad de especie, VV el volumen y x=n/Vx=n/V su concentración. En una mezcla homogénea, QinQ_{\mathrm{in}} y QoutQ_{\mathrm{out}} son caudales volumétricos, cinc_{\mathrm{in}} la concentración de entrada y rr la velocidad volumétrica neta de producción. La mezcla perfecta hace que la concentración de salida sea xx y el balance queda

d(Vx)dt=Qincin−Qoutx+Vr.\frac{d(Vx)}{dt}=Q_{\mathrm{in}}c_{\mathrm{in}}-Q_{\mathrm{out}}x+Vr.

La regla del producto distingue los balances de cantidad y de concentración; al aumentar VV, la concentración puede bajar sin pérdida de soluto. Esa dilución resulta del almacenamiento geométrico. Para un líquido incompresible y volúmenes aditivos, dV/dtdV/dt es la diferencia de caudales; fuera de esas condiciones, el balance de masa debe incluir densidades.

Vdxdt=Qincin−Qoutx+Vr−xdVdt.V\frac{dx}{dt}=Q_{\mathrm{in}}c_{\mathrm{in}}-Q_{\mathrm{out}}x+Vr-x\frac{dV}{dt}.

A volumen fijo, el último término desaparece y los caudales deben respetar esa restricción, pero persiste el lavado de soluto por la salida. La dilución por crecimiento de volumen y la sustitución continua del contenido pueden dar ecuaciones semejantes, aunque describen experimentos y condiciones de frontera distintos.

En un volumen de control fijo, despreciando las energías cinética y potencial, el balance de energía, que contabiliza calor, trabajo y transporte material, es

dUdt=Q˙h−W˙+∑inm˙h−∑outm˙h,\frac{dU}{dt}=\dot Q_h-\dot W+\sum_{\mathrm{in}}\dot m h-\sum_{\mathrm{out}}\dot m h,

donde UU es la energía interna almacenada, Q˙h\dot Q_h el calor entrante, W˙\dot W el trabajo saliente distinto del trabajo de flujo ya incluido en las entalpías, m˙\dot m los caudales másicos y hh las entalpías específicas de las corrientes.

La reacción redistribuye energía química y térmica sin crear energía total. En una ecuación separada de temperatura, el calor de reacción depende de la elección de variables y debe ser coherente con las entalpías y los intercambios térmicos modelados.

171.3. Equilibrio y estado estacionario

Un estado estacionario mantiene constantes las variables macroscópicas que lo describen bajo condiciones externas fijas. Esta definición atañe a la evolución temporal y depende del conjunto de observables elegido. El equilibrio termodinámico exige, además, la desaparición de las fuerzas impulsoras de los procesos irreversibles permitidos. Distinguir ambas nociones evita atribuir a una concentración constante propiedades termodinámicas que su medición no demuestra.

En un balance estacionario, la acumulación se anula porque los términos restantes se compensan. Una entrada puede ser igual a la suma de una salida y un consumo químico. Por tanto, derivadas observables nulas no implican flujos nulos. El sistema puede renovar continuamente su contenido y producir entropía mientras conserva la misma composición media. Una fotografía del estado no revela esa circulación; hacen falta medidas de corrientes o una perturbación controlada.

dXdt=0⟹X˙ent−X˙sal+X˙gen=0.\frac{dX}{dt}=0 \quad\Longrightarrow\quad \dot X_{\mathrm{ent}}-\dot X_{\mathrm{sal}}+\dot X_{\mathrm{gen}}=0.

En esta expresión, XX es una cantidad extensiva y las tres tasas son la entrada, la salida y la generación interna algebraica. La igualdad permite inferir una compensación, pero deja indeterminados sus términos individuales. Incluso medir varias cantidades constantes puede resultar insuficiente para reconstruir las corrientes. La identificación del mecanismo exige información adicional sobre transporte, reacciones y relaciones constitutivas.

En equilibrio no hay corrientes macroscópicas netas debidas a diferencias persistentes de temperatura, potencial químico o presión que puedan relajarse mediante los intercambios admitidos. Los intercambios microscópicos continúan y las fluctuaciones no desaparecen. La condición de equilibrio es la ausencia de sesgos termodinámicos sostenidos. En modelos estocásticos apropiados, el balance detallado expresa esa compensación entre transiciones elementales y sus reversas.

Un sistema abierto puede estar en equilibrio si los reservorios son compatibles, porque el contacto que permite intercambio no obliga a que exista una corriente neta. Cuando coinciden las temperaturas y los potenciales químicos pertinentes y las restricciones mecánicas son compatibles, la apertura puede coexistir con el equilibrio. La permeabilidad de la frontera determina qué diferencias pueden relajarse; una membrana selectiva puede impedir algunos intercambios y permitir otros.

Recíprocamente, un sistema cerrado puede permanecer fuera del equilibrio durante un transitorio o bajo un impulso energético sostenido. Una diferencia de temperatura impuesta a sus paredes puede mantener conducción de calor sin transferencia de materia. La clasificación abierto o cerrado y la distinción entre equilibrio y no equilibrio responden, pues, a preguntas diferentes. Ambas deben especificarse al describir un experimento, junto con la duración y la resolución de la observación.

Estar abierto tampoco basta para la autoorganización. Un flujo puede conducir a un estado homogéneo estable cuya dinámica solo amortigua perturbaciones. La aparición de patrones, oscilaciones o nuevas estructuras requiere mecanismos capaces de amplificar determinadas fluctuaciones y de limitar después su crecimiento. El análisis pertinente combina condiciones de frontera, cinética, transporte y estabilidad; la presencia de una corriente es solo una parte de esa descripción.

171.4. El reactor agitado de volumen fijo

Consideremos un reactor de volumen constante VV, alimentado y vaciado con el mismo caudal QQ. La agitación mantiene uniforme la concentración x(t)x(t), de modo que el efluente tiene esa misma concentración. La alimentación tiene concentración cinc_{\mathrm{in}}. Supondremos una conversión irreversible de primer orden, con consumo volumétrico kxkx y constante cinética positiva kk.

Definimos la tasa de renovación D=Q/VD=Q/V, que, como kk, tiene dimensiones de inverso del tiempo. El balance de cantidad, dividido por el volumen fijo, distingue aporte, extracción y consumo por reacción; la derivada corresponde a la concentración media.

dxdt=D(cin−x)−kx=Dcin−(D+k)x.\frac{dx}{dt}=D(c_{\mathrm{in}}-x)-kx =Dc_{\mathrm{in}}-(D+k)x.

Anulando la acumulación, para D+k>0D+k>0 se obtiene un único estado estacionario x∗x^*, que compensa alimentación, lavado y conversión. Con alimentación y reacción se cumple x∗<cinx^*<c_{\mathrm{in}}; la diferencia sostiene el aporte neto de reactivo al volumen de control.

x∗=DcinD+k.x^*=\frac{Dc_{\mathrm{in}}}{D+k}.

Para derivar la solución completa, definimos y(t)=x(t)−x∗y(t)=x(t)-x^*. Restar el balance estacionario elimina la alimentación constante; el procedimiento equivale a aplicar un factor integrante y deja la ecuación homogénea separable

dydt=−(D+k)y,y(t)=y(0)e−(D+k)t.\frac{dy}{dt}=-(D+k)y, \qquad y(t)=y(0)e^{-(D+k)t}.

Con x0=x(0)x_0=x(0), la reconstrucción da la concentración para todo t≥0t\geq 0. El tiempo característico de relajación τ\tau es el inverso de la tasa total de amortiguamiento; tras un intervalo τ\tau, la desviación inicial se reduce por un factor e−1e^{-1}.

x(t)=x∗+(x0−x∗)e−t/τ,τ=1D+k.x(t)=x^*+(x_0-x^*)e^{-t/\tau}, \qquad \tau=\frac{1}{D+k}.

El único autovalor, −(D+k)-(D+k), es negativo, de modo que el estado estacionario es estable y el acercamiento, monótono. Con concentración inicial físicamente admisible y alimentación no negativa, xx permanece no negativa. La dinámica, de una variable y linealmente amortiguada, no produce oscilaciones ni múltiples estados estacionarios.

La derivación supone parámetros constantes, mezcla suficientemente rápida y temperatura controlada. Una cinética dependiente de la temperatura o una alimentación variable modifica la ecuación y puede exigir variables adicionales. La estabilidad obtenida es propia de este modelo y del régimen experimental que representa; otros reactores abiertos pueden carecer de ella.

171.5. Un caso numérico didáctico

Tres curvas de concentración frente a horas: reactivo y producto convergen a 5 mmol/L y su suma a 10 mmol/L; todas parten de cero bajo los parámetros especificados.
Figura 171.1. Reactivo, producto e inventario combinado en un reactor abierto con D=k=0,5 por hora y alimentación de 10 mmol/L. El inventario total relaja más lentamente que el reactivo. Ampliar SVG ↗
Tasas de entrada, lavado y reacción de un reactor didáctico y compensación de los flujos con la acumulación de concentración.
Figura 171.2. En un reactor alimentado, la concentración estable coexiste con una entrada de cinco y salidas por lavado y conversión de 2,5 cada una, medidas por litro y por hora. Ampliar SVG ↗

Este ejemplo es didáctico; sus valores no proceden de mediciones experimentales. Tomamos una tasa de renovación D=0,5 h−1D=0{,}5\ \mathrm{h}^{-1}, una constante cinética k=0,5 h−1k=0{,}5\ \mathrm{h}^{-1} y una concentración de entrada cin=10 mmol/Lc_{\mathrm{in}}=10\ \mathrm{mmol/L}. La igualdad de tasas facilita seguir el reparto del reactivo entre salida y conversión; caudal y volumen pueden variar si su cociente conserva el valor especificado.

Las expresiones generales dan una concentración estacionaria x∗=5 mmol/Lx^*=5\ \mathrm{mmol/L} y un tiempo de relajación τ=1 h\tau=1\ \mathrm{h}. El asterisco señala el valor estacionario interior. Por mezcla perfecta, la concentración de salida también es de cinco milimoles por litro; la alimentación, más concentrada, compensa la reacción con su aporte neto.

x∗=0,5×100,5+0,5=5 mmol/L,τ=10,5+0,5=1 h.x^*=\frac{0{,}5\times10}{0{,}5+0{,}5}=5\ \mathrm{mmol/L},\qquad \tau=\frac{1}{0{,}5+0{,}5}=1\ \mathrm h.

Por unidad de volumen, la entrada aporta Dcin=5Dc_{\mathrm{in}}=5, la salida transporta Dx∗=2,5Dx^*=2{,}5 y la conversión consume kx∗=2,5kx^*=2{,}5, todas en mmol L−1 h−1\mathrm{mmol\,L^{-1}\,h^{-1}}. Multiplicarlas por el volumen da tasas extensivas de cantidad de sustancia.

dxdt∣x=x∗=(5−2,5−2,5) mmol L−1 h−1=0.\left.\frac{dx}{dt}\right|_{x=x^*}=(5-2{,}5-2{,}5)\ \mathrm{mmol\,L^{-1}\,h^{-1}}=0.

La desaparición del reactivo corresponde a su conversión en producto. Supongamos transformación uno a uno, sin producto en la alimentación ni pérdidas adicionales. La concentración de producto p(t)p(t), que se forma a la tasa volumétrica de reacción y sale por la misma corriente, cumple

dpdt=kx−Dp,p∗=kx∗D=5 mmol/L.\frac{dp}{dt}=kx-Dp,\qquad p^*=\frac{kx^*}{D}=5\ \mathrm{mmol/L}.

En el estado estacionario, la salida transporta reactivo y producto, cada uno a 2,5 mmol L−1 h−12{,}5\ \mathrm{mmol\,L^{-1}\,h^{-1}}; juntos igualan la entrada. La reacción se cancela al sumar los balances. La conservación de la suma de cantidades molares depende de la estequiometría uno a uno; en una reacción general se conservan la masa y los elementos correspondientes. La suma de concentraciones evoluciona hacia la concentración de entrada porque la reacción no hace desaparecer material y solo el intercambio exterior modifica ese inventario conjunto. Una medida del reactivo aislado atribuiría a la reacción una pérdida aparente; incluir el producto y cualquier corriente omitida la corrige y amplía la frontera contable.

La suma de ambas concentraciones muestra otra escala temporal. Si z=x+pz=x+p, la reacción se cancela y queda z˙=D(cin−z)\dot z=D(c_{\mathrm{in}}-z). Por tanto,

z(t)=cin+(z0−cin)e−Dt,p(t)=z(t)−x(t).z(t)=c_{\mathrm{in}}+(z_0-c_{\mathrm{in}})e^{-Dt},\qquad p(t)=z(t)-x(t).

En el ejemplo, el inventario conjunto relaja en una escala de dos horas, mientras que el reactivo lo hace en una de una hora. La concentración del producto combina en general ambas exponenciales. Si el reactor comienza con z0=cinz_0=c_{\mathrm{in}}, el modo de renovación del inventario no aparece y p(t)=cin−x(t)p(t)=c_{\mathrm{in}}-x(t). Identificar el tiempo de un observable con el de todo el sistema puede ocultar esta dependencia de la preparación. Un registro aislado del reactivo puede no mostrar el modo más lento del balance conjunto.

El tiempo de relajación mide la rapidez con que se alcanza el estado estacionario y es distinto del tiempo que cada molécula permanece en el reactor. Desde cualquier concentración inicial admisible, la desviación del reactivo decae exponencialmente con una hora como escala característica. La composición estable coexiste con renovación y conversión continuas, identificables mediante los flujos calculados.

171.6. Costes y recursos para mantener gradientes

Dos curvas de concentración y velocidad de consumo muestran protocolos de cierre y lavado con la misma condición inicial, pero balances distintos.
Figura 171.3. Partiendo de cinco milimoles por litro, cerrar el reactor y sustituir la alimentación por disolvente son protocolos distintos: el primero relaja en dos horas y el segundo en una. Ampliar SVG ↗

Mantener una diferencia de temperatura, composición o potencial químico requiere compensar los procesos que tienden a relajarla. Esta compensación consume recursos disponibles en el entorno y suele generar disipación. En el reactor, la corriente entrante repone reactivo y la saliente retira material transformado. La persistencia de la composición interior depende de esas condiciones externas, de la cinética y de la capacidad del dispositivo para conservar el volumen y la temperatura supuestos.

El coste físico debe especificarse mediante magnitudes concretas. Puede incluir trabajo de bombeo, energía para agitar, transferencia térmica o consumo de sustancias con potencial químico elevado respecto de los productos. La ecuación de concentración no fija ese coste; para calcularlo se necesitan presiones, temperaturas, propiedades termodinámicas y detalles del equipo. Un balance material correcto es compatible con consumos energéticos muy distintos según la realización experimental.

Un quimiostato ideal mantiene constantes ciertas concentraciones o potenciales químicos externos durante toda la evolución. Esta idealización equivale a suponer reservorios suficientemente grandes o mecanismos de reposición que impiden su agotamiento apreciable. Es útil para formular una dinámica autónoma con parámetros fijos. Sin embargo, desplaza al entorno parte del problema, porque el dispositivo que prepara la alimentación y retira los productos también necesita recursos y una descripción física coherente.

Los reservorios reales tienen capacidades finitas. Su composición puede cambiar, sus temperaturas fluctuar y sus bombas alcanzar límites operativos. Si esas variaciones ocurren en escalas comparables a la respuesta del reactor, la alimentación deja de ser una condición constante y debe incorporarse a la dinámica. El estado estacionario del modelo reducido puede entonces convertirse en una aproximación temporal, válida mientras el agotamiento o la deriva de los parámetros sean pequeños.

Retirar la alimentación conduce a otro régimen, cuya forma depende del protocolo. Si se detienen a la vez entrada y salida y se conserva el volumen, queda un recipiente cerrado con reacción. Si la alimentación se sustituye por disolvente sin reactivo, continúa el lavado junto con la conversión. En ambos casos disminuye el reactivo, aunque las escalas temporales y los balances externos difieren. El protocolo debe describirse antes de interpretar la respuesta posterior.

dxdt=−kx(Q=0),dxdt=−(D+k)x(cin=0, D>0).\frac{dx}{dt}=-kx \quad(Q=0), \qquad \frac{dx}{dt}=-(D+k)x \quad(c_{\mathrm{in}}=0,\ D>0).

Aquí xx es la concentración de reactivo, kk la constante cinética, QQ el caudal y D=Q/VD=Q/V la tasa de renovación del volumen VV. Las ecuaciones suponen mezcla uniforme y temperatura constante. Un sistema térmicamente no controlado necesitaría además un balance de energía, porque detener corrientes también modifica la extracción de calor y, con ella, la velocidad de reacción.

El no equilibrio es una noción distinta del caos y de la complejidad. La alimentación continua del reactor lineal mantiene una corriente y una conversión sin producir dinámica irregular. La complejidad temporal o espacial depende de retroalimentaciones, acoplamientos y condiciones de estabilidad específicas. Cuantificar los recursos que sostienen un gradiente y establecer el mecanismo que genera una organización son problemas conexos, aunque requieren evidencias y ecuaciones distintas.

171.7. Primera y segunda ley

La primera ley expresa la conservación de la energía en un balance de almacenamiento, transporte, calor y trabajo. La segunda restringe la dirección de los procesos y la capacidad de producir trabajo. En sistemas abiertos, una energía almacenada constante puede coexistir con transferencias y producción de entropía; la estacionariedad energética no implica reversibilidad.

Para un sistema con frontera y variables termodinámicas definidas, el balance de entropía separa el intercambio de la irreversibilidad. SS es la entropía del sistema, S˙e\dot S_e la tasa neta recibida del entorno y S˙i\dot S_i la producción interna. La entropía entrante se cuenta como positiva y la producción interna satisface la desigualdad macroscópica.

dSdt=S˙e+S˙i,S˙i≥0.\frac{dS}{dt}=\dot S_e+\dot S_i,\qquad \dot S_i\geq0.

La entropía se intercambia mediante corrientes de materia y de calor. En cada contacto térmico a temperatura de frontera TjT_j, el aporte es Q˙j/Tj\dot Q_j/T_j, con Q˙j\dot Q_j el calor entrante; las corrientes materiales llevan caudales másicos m˙\dot m y entropías específicas ss. El índice jj recorre los contactos térmicos y las demás sumas recorren entradas y salidas.

S˙e=∑jQ˙jTj+∑inm˙s−∑outm˙s.\dot S_e=\sum_j\frac{\dot Q_j}{T_j}+\sum_{\mathrm{in}}\dot m s-\sum_{\mathrm{out}}\dot m s.

En estado estacionario, con inventario de entropía constante, frontera fija y contabilidad completa, la producción interna se compensa con una exportación neta; con disipación, el intercambio es negativo. Una concentración constante no basta, porque también deben ser estacionarias las variables térmicas y las demás contribuciones al almacenamiento.

dSdt=0⟹S˙e=−S˙i.\frac{dS}{dt}=0\quad\Longrightarrow\quad\dot S_e=-\dot S_i.

El sistema puede conservar su organización macroscópica mientras aumenta la entropía del conjunto sistema-entorno. La organización no implica una entropía arbitrariamente pequeña, porque la entropía termodinámica depende de estados, propiedades materiales y resolución física. Un patrón visual regular puede coexistir con temperaturas elevadas, corrientes intensas o una producción de entropía considerable.

En sistemas pequeños, las fluctuaciones exigen precisar el nivel de descripción. La termodinámica estocástica relaciona trayectorias, intercambios y producción de entropía; en ella se distinguen los resultados fluctuantes de las restricciones sobre promedios válidas bajo hipótesis definidas (Seifert, 2012). Su uso requiere identificar reservorios, mecanismos y condiciones experimentales antes de atribuir significado termodinámico a una observación.

No hay un principio universal de máxima producción de entropía que prediga cualquier organización fuera del equilibrio; los principios extremales tienen dominios y restricciones específicos. Predecir estados y estabilidad exige estudiar las ecuaciones dinámicas y las condiciones de contorno. El balance solo restringe las evoluciones admisibles; seleccionar patrones requiere además información constitutiva y análisis de estabilidad.

171.8. Escalas de tiempo y número de Damköhler

Gráficos de fracciones estacionarias de reactivo lavado y convertido, y tiempo de relajación en horas según número de Damköhler.
Figura 171.4. En el reactor ideal, el número de Damköhler determina las fracciones que salen sin convertir y que reaccionan; al crecer Da, disminuye también la escala de relajación con D fijo. Ampliar SVG ↗

La comparación de escalas temporales permite interpretar un balance antes de resolverlo y evaluar sus aproximaciones. En el reactor aparecen el tiempo medio de residencia, el tiempo de reacción y el tiempo de mezcla, que corresponden, respectivamente, al transporte por la corriente, a la transformación química y a la homogeneización del contenido. Su jerarquía determina qué procesos pueden tratarse como rápidos y qué variables necesitan una descripción espacial o temporal explícita.

Para volumen fijo VV y caudal QQ, el tiempo medio de residencia es τres=V/Q=1/D\tau_{\mathrm{res}}=V/Q=1/D, donde DD es la tasa de renovación. En mezcla perfecta no todas las moléculas permanecen ese tiempo, porque la salida extrae una muestra representativa del contenido; la residencia media solo caracteriza el recambio global del reactor.

Un trazador inerte inyectado como pulso permite precisar la distribución de tiempos de residencia. En el reactor ideal cada partícula sale con tasa constante DD, de modo que la supervivencia es e−Dte^{-Dt} y la densidad de tiempos de residencia es E(t)=De−DtE(t)=De^{-Dt}, para t≥0t\geq0. Su integral vale uno y su media es 1/D1/D. Hay partículas que salen casi de inmediato y otras que permanecen varias veces la media. Con la curva de un trazador real se contrasta esa hipótesis de mezcla sin depender de la cinética del reactivo.

Si una molécula de reactivo puede convertirse con tasa kk, compiten dos eventos, la salida sin transformarse y la reacción. La supervivencia como reactivo dentro del tanque es e−(D+k)te^{-(D+k)t}. La probabilidad de convertirse antes de salir resulta

∫0∞ke−(D+k)t dt=kD+k.\int_0^\infty k e^{-(D+k)t}\,dt=\frac{k}{D+k}.

Esa probabilidad coincide con la fracción de conversión del balance estacionario. La duración media hasta cualquiera de esos eventos es 1/(D+k)1/(D+k), pero la materia transformada puede seguir dentro del tanque como producto. Por eso esa duración y la residencia del trazador responden a preguntas distintas.

Para una reacción de primer orden con constante kk, definimos el tiempo de reacción τrxn=1/k\tau_{\mathrm{rxn}}=1/k. El número de Damköhler compara el tiempo de residencia con esa escala química. En este modelo su definición usa solo las dos tasas ya presentes; otras geometrías o leyes cinéticas pueden requerir definiciones diferentes. Al comunicar un valor de este número adimensional debe indicarse siempre qué tiempo de transporte interviene.

Da=τresτrxn=kD,x∗cin=11+Da.\mathrm{Da}=\frac{\tau_{\mathrm{res}}}{\tau_{\mathrm{rxn}}} =\frac{k}{D}, \qquad \frac{x^*}{c_{\mathrm{in}}}=\frac{1}{1+\mathrm{Da}}.

Cuando Da≪1\mathrm{Da}\ll1, la reacción es lenta respecto del recambio. Gran parte del reactivo abandona el recipiente antes de convertirse y la concentración estacionaria se aproxima a cinc_{\mathrm{in}}. Se habla de un régimen limitado por la cinética frente al tiempo de residencia disponible. Esa lentitud se mide frente al lavado; respecto de otras escalas externas, la misma reacción puede ser rápida.

Cuando Da≫1\mathrm{Da}\gg1, la conversión es rápida comparada con la renovación y la concentración interior de reactivo resulta pequeña. La fracción convertida aumenta y la alimentación limita el aporte disponible. Esta conclusión vale para un reactor ideal homogéneo; en un equipo real, una cinética muy rápida puede consumir reactivo antes de que se mezcle y hacer relevantes perfiles locales que el modelo omite.

τmix≪min⁡(τres,τrxn),τ=1D+k.\tau_{\mathrm{mix}}\ll \min(\tau_{\mathrm{res}},\tau_{\mathrm{rxn}}), \qquad \tau=\frac{1}{D+k}.

Aquí τmix\tau_{\mathrm{mix}} es el tiempo de mezcla y τ\tau la relajación de la concentración media. La desigualdad es un criterio práctico suficiente para sostener la aproximación de mezcla rápida. Deben evaluarse además zonas muertas, cortocircuitos y transporte cerca de superficies reactivas. Que haya un agitador no demuestra que todas las escalas relevantes cumplan esta separación.

La relajación combina reacción y lavado porque ambos amortiguan una desviación del reactivo; de ahí que pueda ser más rápida que cada una de las escalas individuales. Medir esa relajación da la suma de tasas, mientras que separar sus contribuciones exige conocer el caudal, el volumen u otra observación independiente. Esta distinción relaciona el análisis dimensional con la identificación experimental de los parámetros.

171.9. Organizaciones sostenidas y transitorios

Una organización sostenida fuera del equilibrio necesita una fuente que mantenga diferencias, procesos disipativos y una dinámica compatible con su persistencia. En el caso de los patrones espontáneos hace falta también un mecanismo que vuelva inestable un estado previo o seleccione una estructura entre varias posibles. La alimentación y la disipación permiten sostener corrientes, pero la forma concreta de la organización depende de acoplamientos, geometría y condiciones de contorno.

El reactor lineal estudiado muestra el alcance limitado de esas condiciones generales. Recibe reactivo, lo transforma y elimina material continuamente, pero todas las perturbaciones de concentración decaen. Con una sola variable, carece de una retroalimentación capaz de producir una inestabilidad oscilatoria. La estacionariedad disipativa está plenamente representada, mientras que la hipótesis de mezcla uniforme y la ausencia de coordenadas espaciales excluyen los patrones espaciales.

La convección térmica ilustra una situación distinta. Una capa de fluido calentada de forma apropiada puede perder la estabilidad del estado conductor cuando el impulso supera un umbral determinado por propiedades materiales y condiciones de frontera. La flotación amplifica los movimientos, mientras que la viscosidad y la difusión térmica los amortiguan. El inicio de la convección depende de esa competencia; una diferencia de temperatura cualquiera no garantiza la formación de celdas.

Las paredes, la profundidad de la capa y las condiciones térmicas y mecánicas de sus límites participan en la selección de modos. El patrón observado tampoco tiene que coincidir con el modo inicialmente más amplificado cuando intervienen efectos no lineales y competencia entre estructuras. La teoría de formación de patrones estudia precisamente estas relaciones entre inestabilidad, simetrías, escalas espaciales y saturación dinámica (Cross y Hohenberg, 1993).

Las oscilaciones químicas requieren mecanismos adicionales, como acoplamientos entre especies y retroalimentaciones con escalas diferentes. Un flujo constante puede mantener disponibles los reactivos, pero la existencia de un ciclo depende de la cinética y de los parámetros. El trabajo de Prigogine y Lefever (1968) es un referente teórico sobre inestabilidades en sistemas químicos disipativos. Su relevancia no hace de la oscilación una consecuencia automática de la apertura.

Un transitorio puede mostrar estructuras llamativas sin alcanzar una organización sostenida. Un frente de reacción, una separación temporal de regiones o una oscilación amortiguada pueden consumir diferencias almacenadas inicialmente. Para distinguir estos casos hay que comparar la duración observada con las escalas de relajación y comprobar el papel de los recursos externos. Una trayectoria prolongada puede parecer estacionaria si el intervalo de medida es demasiado breve respecto de la deriva real.

La estabilidad se establece alrededor de un estado o solución definidos, bajo perturbaciones y restricciones especificadas. En sistemas con varias variables, el análisis lineal identifica el crecimiento o el amortiguamiento inicial; las no linealidades determinan la evolución posterior y pueden limitar la amplitud. La interpretación física necesita relacionar explícitamente fuentes, disipación y mecanismo de inestabilidad. Las categorías abierto, complejo y autoorganizado solo adquieren precisión cuando se acompañan de esa descripción dinámica verificable.

171.10. Identificación experimental y límites del modelo

Identificar un estado sostenido requiere medir variables y corrientes bajo condiciones de frontera documentadas. En el reactor se registran concentración, caudal, volumen, temperatura y composición de entrada, con su incertidumbre y resolución temporal. Una señal estable no demuestra ausencia de actividad interna. Las medidas simultáneas permiten comprobar compensaciones en el balance y distinguir una meseta de una deriva lenta o de un promedio que oculta fluctuaciones.

La respuesta a una perturbación da información que el estado estacionario deja indeterminada. Con alimentación constante, la relajación permite estimar la suma de las tasas de renovación y de reacción; si caudal y volumen se conocen por separado, puede inferirse la constante cinética dentro del modelo. El ajuste debe capturar el transitorio, separar la asíntota de la amplitud inicial y no tratar parámetros correlacionados como estimaciones independientes.

Con un pequeño cambio de la alimentación se contrasta de forma independiente la concentración estacionaria. Si la concentración de entrada cambia bruscamente y después permanece constante, el nuevo estado estacionario se calcula con la misma fracción de conversión, mientras que la diferencia respecto de él decae con la suma de las tasas de renovación y de reacción. Medir a la vez esa suma y la fracción estacionaria permite separar ambos parámetros dentro del modelo. Una sola concentración estacionaria, por el contrario, identifica principalmente su razón.

Este procedimiento requiere saber cuándo llegó realmente la nueva alimentación al reactor. Un tubo de entrada con volumen apreciable introduce un retraso de transporte; una mezcla imperfecta puede añadir varias escalas de relajación. Ajustar toda la señal con una exponencial sin considerar esos elementos puede atribuir a la química un tiempo generado por el aparato. Medir el caudal, el volumen de conexión y un trazador que no reaccione ayuda a evaluar esa explicación. Los balances siguen siendo válidos, pero la descripción de un único compartimento puede necesitar ampliarse para representar la medición y el transporte.

Apagar el impulso externo permite investigar la dependencia de recursos continuos. El protocolo distingue entre detener las corrientes, sustituir la alimentación por disolvente y eliminar el control térmico, porque cada intervención modifica términos diferentes de los balances. Una persistencia finita puede reflejar reservas almacenadas; atribuir autonomía exige compararla con el agotamiento y con las escalas de relajación pertinentes.

Un balance extendido incorpora productos, corrientes secundarias y almacenamiento en superficies o fases omitidas. Los invariantes estequiométricos detectan errores contables, porque las combinaciones de cantidades correspondientes a elementos conservados solo cambian por transporte exterior. Un residuo sistemático puede señalar una especie no medida, un caudal incorrecto o una hipótesis inadecuada; conviene revisar la descripción antes de invocar mecanismos nuevos o excepciones a la conservación material.

El dominio de validez comprende mezcla rápida, volumen fijo, cinética de primer orden y parámetros aproximadamente constantes. Fuera de él pueden aparecer limitaciones de transporte, reversibilidad química, variación térmica o múltiples especies acopladas. Extrapolar exige comprobar estas modificaciones y su efecto sobre la estabilidad, ya que un buen ajuste en una ventana experimental no garantiza capacidad predictiva ante cambios de escala, alimentación o temperatura.

Un modelo físico no autoriza metáforas sociales sobre la entropía. Aplicar vocabulario termodinámico a otro dominio requiere definir variables, unidades, fronteras, mecanismos de intercambio y relaciones contrastables. La regularidad visual o una medida estadística de diversidad no equivalen automáticamente a entropía termodinámica. La precisión conceptual mantiene unidos la definición matemática, la operación de medida y las restricciones físicas, sobre todo si una analogía pretende sostener una explicación causal.

Referencias

Seifert, U. (2012). Stochastic thermodynamics, fluctuation theorems and molecular machines. Reports on Progress in Physics, 75, 126001. DOI.

Cross, M. C. y Hohenberg, P. C. (1993). Pattern formation outside of equilibrium. Reviews of Modern Physics, 65, 851–1112. DOI.

Prigogine, I. y Lefever, R. (1968). Symmetry breaking instabilities in dissipative systems. II. The Journal of Chemical Physics, 48, 1695–1700. DOI.