7. Estabilidad de sistemas termodinámicos. Principio de le Chatelier * Hasta ahora hemos trabajado frecuentemente con la condición de equilibrio dS = 0 a U = cte, o dU = 0 a S =cte. Similarmente mediante otras funciones termodinámicas. *En este capítulo p q queremos explorar p las consecuencias de exigir g q que el equilibrio sea estable: d2S < 0, a U =cte o d2U > 0 a S = cte Adelantamos las p principales p conclusiones,, que q son: 1) cP cV 0 T S 0 2) Principio de le Chatelier: Un sistema estable, al aplicarle cualquier perturbación (o por una fluctuación espontánea) éste reacciona oponiéndose a la perturbación inicial e intentando restaurar el equilibrio equilibrio. Un sistema inestable reacciona apartándolo más todavía del equilibrio. Estabilidad en Mecánica Energía potencial U(x) F r r F Fr x x x0 Sistema estable x x Sistema inestable Sistema localmente estable o metaestable Desplazamiento producido por medios externos: x (se supone eq eq. en x0 =0) Fuerza aplicada por el sistema como respuesta a la perturbación exterior Equilibrio: Estabilidad: Fr 0 Fr dU dx dU 0 dx F r ( x 0) 0 dF d 2U 0 0 dx dx 2 F r ( x 0) 0 U cóncava: estable Localmente estable Estabilidad intrínseca de los sistemas termodinámicos. Volumen Consideremos dos sistemas iguales cuya ecuación fundamental S(U,V,N) es como la figura y separados por una pared móvil. S(U) cóncava debe ser convexa Por simetría, simetría la posición con la pared en el centro y V1 = V2 = V es de equilibrio. Nos preguntamos si la pared puede moverse a q con los volúmenes otra situación de equilibrio desiguales. Si S(V) fuera cóncava la mover la pared hacia la izquierda el sistema 1 tendrá un volumen V-V y el 2 V+V La respuesta es que ocurrirá si la entropía es mayor que la inicial. Veamos, la entropía inicial es: Y la final: V 0(V pequeño) : 1 2S S S (V V ) S (V ) V V 2 2 V 2! V S 1 2S V 2 S (V V ) S (V ) V 2 V 2! V 1 2S S (V V ) S (V V ) 2 S (V ) V 2 2 2 V S i 2 S (U , V , N ) S f S (U ,V V , N ) S (U ,V V , N ) 2 S (U ,V , N ) El sistema se descompondría en 2 partes de diferente v (o densidad) Al contrario, si S = convexa: S (U , V V , N ) S (U , V V , N ) 2 S (U ,V , N ) 2S 2 0 V U , N Estabilidad global Estabilidad local Estabilidad intrínseca de los sistemas termodinámicos. Energía S(U) cóncava debe ser convexa Consideremos dos sistemas iguales cuya ecuación fundamental S(U,V,N) es como la figura g y separados p p por una p pared diatérmana. U ? Si S(U) fuera cóncava pasaría una cantidad de energía U (en forma de calor) de uno a otro aumentando la entropía total. Si 2 S (U , V , N ); S f S (U U , V , N ) S (U U , V , N ) 2 S (U , V , N ) Si U 0(U pequeño) : El sistema se descompondría en dos con diferente U S 1 S U U 2 2 U 2! U S 1 2S S (U U ) S (U ) U U 2 2 U 2! U 1 2S S (U U ) S (U U ) 2 S (U ) U 2 2 2 U S (U U ) S (U ) 2 Lo importante sobre todo es que serían diferentes u y s (densidades de energía y de entropía) S (U U , V , N ) S (U U , V , N ) 2 S (U , V , N ) 2S 0 2 U V , N Estabilidad global Estabilidad local R l de Regla d lla palanca l Sea x un parámetro extensivo molar cualquiera Ec fundamental subyacente s(x) deducida de la Mecánica Estadística imponiendo sólo la cond. de equilibrio Sistema en zona roja (inestable) se descompone en dos partes, en AyB B, con distinto valor de x Sistema en zona ámbar (estado C, localmente estable) puede permanecer en C, pero con una perturbación (o fluctuación) suficientemente fuerte se descompondrá en 2partes y ya no volverá a C. Sistema en zona verde: estable Semáforo de la estabilidad La línea recta tangente en A y en B representa la entropía molar media del sistema más estable pero descompuesto en dos partes con diferente valor del parámetro intensivo molar x Es decir, si x = xC = XC/N, el estado más estable no es el C, homogéneo, sino que el sistema se descompone en dos, una parte en A y otra en B, de modo que: NA xB xC N xB x A NA Y la entropía no es NsC sino: xC x A N xB x A S s A N A s B N B NsC ' NsC V i variables. Varias i bl Ej. Ej U y V S(U,V,N) CONVEXA exige que la matriz “hessiana": Sea definida negativa, para lo cual : 2S 2 U 2 S UV 2S UV 2S V 2 2S 1 T 1 1 / T 0 2 2 2 U U T U T Nc V ,N V ,N V , N V cV 0 Y además, como los dos valores propios son negativos, el determinante es positivo: 2S U 2 2S UV También 2S 2 2 2 UV S S S U 2 V 2 2S V , N V , N UV V 2 2S P / T 2 0 V V U ,N U , N 2 0 2S P / T 2 0 0 V U , N V U , N , pero no añade nada nuevo C di i Condiciones de d estabilidad t bilid d para llos potenciales t i l ttermodinámicos di á i Energía: Si S(U,V,N) S(U V N) es convexa convexa, U(S,V,N) U(S V N) debe ser CÓNCAVA 2U T T 2 0 cV S S Nc N V ,N V , N V 0 2U S 2 2U SV 2U 2 S2 U SV 2U SV 2U V 2 Sea definida p positiva 2U 2 2 2 SV U U U S 2 V 2 2U V , N V , N SV V 2 2U 1 P 2 0 V S , N Vk S V S , N Transformaciones de Legendre: Sea X un parámetro extensivo y P U P U P X X P X 2U P 2U 2U 2U P 2 ; 2 0 0 P P 2 P 2 X X 2 0 kS 0 U X 1 Por tanto, para F, H y G: Función de Helmholtz: 2F 2F Nc S 2 0 2 V 0 T T V , N T V , N T V , N cV 0 2F 2F 1 P 2 0 2 0 V V V Vk T , N T , N T ,N T kT 0 Entalpía: p 2H 2 S T T 0 0 S P , N NcP P,N cP 0 2H 2 P V 0 k SV 0 P S , N S ,N kS 0 Energía libre de Gibbs: 2G Nc S 2 0 P 0 T T P , N T P , N cP 0 Por otra parte se ha visto que: 2G V 2 0 kT V 0 P P T , N T , N kT 0 Tv 2 cP cV cP cV 0 kT k S cV kT k S 0 kT cP Principio p de le Chatelier y le Chatelier-Braun: Ejemplo introductorio Consideremos en primer lugar un sistema (derecha) que puede intercambiar calor y trabajo con un reservorio (Figura ) Si aumentamos el volumen en Vf (externamente o por una fluctuación) la presión disminuye. El reservorio actúa moviendo la pared intermedia hacia la derecha hasta restablecer la igualdad de presiones. 1) Produce una disminución de volumen: Vr de signo opuesto a Vf (= -Vr si el reservorio es muy grande) Pr. de le Chatelier Veamos qué pasa con otro parámetro que no sea el que hemos modificado ni su conjugado de Legendre. 2) Si > 0 all aumentar t ell volumen l su ttemperatura t h ha di disminuido i id y ell sistema i t recibe ib calor del reservorio aumenta la presión: También en sentido contrario de la perturbación inicial, que disminuía la presión. Pr. de le Chatelier-Braun Si <0 el sistema da calor al reservorio reservorio, pero también aumenta la presión presión. Pero…un sistema inestable aumenta la perturbación inicial Otro ejemplo: Una reacción química se encuentra en su situación de equilibrio. Si añadimos más cantidad de un n componente (Ni f > 0) el eq equilibrio ilibrio se ve e alterado y parte de lo que hemos añadido reacciona con los demás componentes, disminuyendo la cantidad de ese componente añadido. (Ni r > 0) DEMOSTRACIÓN: (recordemos que en formulación entrópica los parámetros intensivos son 1/T y P/T, ponemos con superíndice res magnitudes del reservorio y sin i nada d llas d dell sistema) i t ) Sean Vf y Uf los cambios pequeños en V y U del sistema por la acción inicial ((en general g el superíndice p f se refiere a cambios p producidos p por la acción inicial o fluctuación). P1 y T1 la temperatura y presión del sistema alcanzadas debido a la acción inicial. Sean Vr y Ur los cambios producidos en el sistema como respuesta del reservorio (y en general el superíndicer r se refiere a la respuesta). La respuesta del reservorio se produce de modo que Utotal = cte y Vtotal =cte de modo que Stotal aumenta S total S S res res 1 1 P P U r 1 V r res U r res V r T1 T1 T T 1 P1 P res r 1 1 P r res U res V U r V r T T T1 T T1 T f f Stotal 1 0 U r 0 (A) T f f y (? , no es obvio) P V r 0 T (será cierto si Ur y Vr son independientes, pero no necesariamente) (B) T Tomemos la l 2ª desigualdad d i ld d (B): (B) f d P P f r f r V r V V 0 V V 0 dV T T y también: d P 2S 0(convexidad ) dV T V 2 U , N Multipliquemos por d dV V f V r 0 le Chatelier P 0 T d P P / T P P f d P r f P / T r f r V V V V 0 P P 0 f r dV T dV T V V T T f r f r T Tomemos la l 1ª desigualdad d i ld d (B): (B) f d 1 U r dU T 1 f r f r U U 0 U U 0 T d 1 S 0(convexidad ) dU T U 2 le Chatelier-Braun 2 Multipliquemos por r d 1 (1 / T f ) r f (1 / T ) f d 1 r U U U U dU T dU T U f U r 1 1 f r 0 T f T r 0 T T T f 1 1 f ; r 2 T T Tf T r r 2 T Principio p de le Chatelier y le Chatelier-Braun: Formulación general en representación entrópica Consideremos un sistema en contacto con un reservorio. Supongamos que por una acción exterior producimos un cambio dX1f en el parámetro extensivo X1 del sistema. Entonces el correspondiente parámetro intensivo P1 cambia : dP1 f P1 dX 1f X 1 La fluctuación también afecta a otro parámetro intensivo P2 : Entonces se cumple: Y también: dP1 f dP1r 0 dP2 f dP2r 0 dP2f P2 dX 1f X 1 LE CHATELIER LE CHATELIER-BRAUN Se demuestra (?) exactamente igual que en el caso particular de X1 = U y X2 = V usando la condición de convexidad de S(U, X1,X2,…) P i i i de Principio d le l Chatelier Ch t li y le l Chatelier-Braun: Ch t li B Fl Fluctuaciones t i Observar: Las fluctuaciones son la esencia de la Termodinámica Ej: en un gas la densidad (volumen específico) es el promedio del número de moléculas que hay en un volumen infinitesimal. * Las fluctuaciones aparecen en pequeñas regiones de un sistema aislado, manteniéndose la energía constante. Un región de esas se puede considerar "subsistema" subsistema y el resto (mucho mayor) "reservorio" reservorio . En un sistema estable la respuesta del resto del sistema es anular las *En fluctuaciones espontáneas que se producen. * En un sistema inestable la respuesta del resto tiende a aumentar la fluctuación y termina por apartar todo el sistema de la situación de equilibrio inicial. C íti del Crítica d l principio i i i de d le l Chatelier Ch t li •No hay duda de su validez si el reservorio es infinitamente grande y la perturbación infinitesimal infinitesimal. De hecho los parámetros intensivos del reservorio no varían y en el equilibrio final los del sistema son iguales a los iniciales, o sea que la respuesta no sólo es de signo contrario sino igual en valor: Pr = -Pf. Esto en particular indica la estabilidad de un sistema respecto de fluctuaciones en pequeñas partes del mismo. •En cambio no siempre se cumple cuando el reservorio es finito o el cambio grande aunque sí la mayoría de las veces. •En particular si P es la presión no necesariamente PfPr <0 •En E todos t d estos t casos excepcionales i l ell comportamiento t i t obedece b d all postulado II, pero de ahí no se deduce el principio de le Chatelier Ejemplo 1 El sistema es un sólido con gran T y pequeña (es de goma) goma). El reservorio es un gas ideal a P = P1 y T = T1. Se suministra una energía U al sistema en forma de calor, a V =cte. Luego la pared se hace móvil y diatérmana. Cambios iniciales: U N s CV s T f T f U 0 N s CV s U P P f 0 U (reducir ) N C U V T s Vs Vamos a poner datos numéricos razonables: V1s = 20 cm3, P1 = 105 Pa (1 atm) , U = 1000 J, Ns=Ng = 1 mol, T1 = 300 K, CVg = 3R/2, CVs = 3R, = 10-5K T = 10-7 Pa-1 (y V1g = 24.94 24 94 litros) Se tiene pues: Tf = 1000/(138.3145)= 40.1 K Pf = (10-5/10-7)1000/(138.3145)= 4009 Pa Es decir el sólido (sistema) ha alcanzado una temperatura T2 = 340.1 K presión P2 = 104009 Pa Ahora hacemos la p pared móvil y diatérmana. En la nueva situación de equilibrio la temperatura T3 y presión P3 del gas y del sólido son iguales y la energía se conserva respecto de la situación. El cambio de energía del sólido es: U sr NC Ps PV T r P T T VP r NC Ps T3 T2 Y ell d dell gas U gr N g CVg T3 T1 Debe ser U gr U sr 0 Luego: 1 2 N g CVg T3 T1 N s C Ps (T2 T3 ) T3 T1 T2 326.7 K 3 3 r Ts T3 T2 13.4 K Luego: Tsr Ts f 13.4 K 40.1 K 0 (LeCh OK) Pero también debe ser: Pero…también Vgr Vsr 0 Vg 3 Vg1 N RT T Y P3 g 3 3 P1 108900 Pa Vg 3 T1 r f r P P3 P2 108900 104009 4891 Pa 0 P P 0 !!! En cambio sí q que ocurre q que: P3 P2 P2 P1 P 2 Pa P P 27.5 27.5 2 0 K T T T3 T2 T2 T1 r f (Chemical Engineering Education, 2003) Ejemplo 2 (Contra-ejemplo) Sea la reacción entre gases: N2+3H23NH3 Consideremos el N2 como “sistema” y el resto de componentes como reservorios además de una fuente de trabajo y de calor que mantienen T, reservorios, T y P =ctes. La Termodinámica predice (y así sucede en efecto, como se verá más tarde) que a P y T =ctes, si la fracción molar de nitrógeno x(N2) <0.5 al añadir una pequeña cantidad de N2 el equilibrio se desplaza hacia la derecha, consumiendo N2 (de acuerdo con Le Chat.), pero si x(N2) >0.5, all añadir ñ di un poco de d N2 la l reacción ió se d desplaza l h hacia i lla iizquierda, i d produciendo más N2. La razón es que las x(N2), x(H2) y x(NH3) no son independientes (están ligadas por la ecuación de reacción) y la situación de mínima G no siempre corresponde con la predicción de Le Chatelier. En el ejemplo anterior añadir N2 implica aumentar x(N2) pero disminuir x(H2). Este segundo efecto prevalece si x(N2) es grande (reservorio pequeño). La reacción se desplaza hacia la izquierda produciendo más H2, a costa de d aumentar más á todavía d í x(N (N2). )
© Copyright 2024