La Teor´ eor´ıa de Flory–Huggins Flory–Hu ggins en la Ingenier´ Ingenier´ıa de Soluciones Solucio nes y Mezclas Mezc las de Pol´ Pol´ımeros. ımer os.
Recibido: 15 de enero de 2008. Aceptado: 15 de mayo de 2008.
tria de pol´ımeros ımeros sino que tambi´en en refleja el atractivo que representa un tema te´orico or ico dif´ıcil. ıci l.
Resumen
Es bien conocida la carencia de una teor´ teor´ıa satisfactoria para explicar todas las propiedades de los sistemas pol p ol´´ımero–disolvente, ımero–dis olvente, a pesar del aluvi´ a luvi´on de trabajos te´oricos oricos en las ´ultimas ultimas d´ecadas. ecadas. Sin embargo, es indudable que se ha obtenido un progreso considerable e ininterrumpido, tanto en la propuesta de nuevos enfoques como enmendando los defectos de las teor´ te or´ıas ıa s m´as as antiguas.
El ´exito exito en la modelaci´on on de procesos industriales depende de forma decisiva de una descripci´on eficaz de las propiedades termodin´amicas amicas y del comportamiento de fase de los sistemas qu´ımicos ımicos asociados. En particular, la industria de los pol´ pol´ımeros, tanto por su relevancia econ´omica omica como por su dificultad te´orica orica inherente, requiere de modelos termodin´amicos amicos apropiado apropiadoss para la pr´actica actica ingenieril. En este trabajo, se ofrece una perspectiva introductoria del enfoque cl´asico asico de Flory–Huggins para describir la termodin´amica amica de soluciones y mezclas polim´ ericas. ericas. Se muestran aplicaciones aplicacion es de inter´es es desde el punto de vista industrial.
El primer modelo de soluciones soluc iones polim´ p olim´ericas, ericas, debido a Paul Flory y Maurice Huggins, es todo un cl´asico en el ´area area y la raz´on on de esto es simple: el criterio b´asico para p ara juzgar el ´exito exito de una teor´ıa ıa dada, d ada, es la concordancia o no con las observaciones experimentales. Al mismo tiempo, se espera que las teor´ teor´ıas exitosas proporcionen informaci´on on fundamental acerca del sistema en cuesti´on, on, tal como el tama˜no no y la forma del soluto en presencia del disolvente, en todo el intervalo de concentraci´on. on. Sin embargo, incluso si una teor´ıa ıa posee p osee todo to do lo anterior, muchas veces puede no ser del todo ´util util desde el punto de vista pr´actiactico. Por ejemplo, en aplicaciones industriales, si se trata de encontrar el disolvente m´as apropiado para un pol´ pol´ımero particular, lo m´as as probable es que se elija un m´ etodo etodo de selecci´on on basado en par´ameametros de solubilidad, antes que utilizar una teor´ teor´ıa con fundamento riguroso en la mec´anica anic a estad´ e stad´ıstica. ısti ca. En este sentido, la teor´ teor´ıa de Flory–Huggins (FH) ofrece un buen balance entre la profundidad de sus bases te´oricas oricas y su tratabilidad en aplicaciones de ingeni ge nier´ er´ıa. ıa .
Palabras Palabras clave: clave: Termodin´ amica,, soluci amica solucione oness po-
lim´ ericas, ericas, Flory-–Huggins, Flory-–Hu ggins, aplicaciones. aplicacion es. Abstract
The success in industrial process modeling is dependent of decisive form of a good description of thermodynamical properties and phase behavior of chemical systems associated. Particularly, the polymer industry industry,, for its economic economic releva relevance nce and his inherent theoretical difficulty, require of thermodynamical models proper for engineeri engineering ng calculati calculations. ons. In this paper, we offer an introductory perspective of Flory–Hug Flory–Huggins gins classical approach approach in order to describe the polymer solutions and blends thermodynamics. We show interesting applications of industrial point of view.
Espec´ıficamente, ıficamente, en el estudio de sistemas polim´eriericos en disolu disoluci´ ci´ on, on, tres tres son las pregun preguntas tas imporimportantes: ¿cu´ales ales disolventes pueden disolver a cu´ales ales pol po l´ımeros ıme ros?, ?, ¿c´ ¿ c´omo omo influyen las interacciones pol´ımeımero–disolvente en las propiedades de la soluci´on?, y finalmente: ¿qu´ ¿q u´e aplicaciones aplica ciones se tienen tie nen a partir de las propiedades de las soluciones p olim´ ericas? ericas? En concreto, la solubilidad o no de pol´ pol´ımeros en determi-
Introducci´ on on
En el mundo, existe un amplio inter´es es en la industria de pol´ pol´ımeros. ımeros. Incluso Incluso se ha llegado llegado a afirmar, afirmar, antes antes del boom de las computado computadoras, ras, que se viv´ viv´ıa en la era del pl´astico, astico, por la gran presencia que estos productos ten´ ten´ıan y a´un un tienen en las actividades humanas. En realidad, este gran inter´ inter´es es refle ja no s´olo olo la fuerte relevancia econ´omica omica de la indus54
La Teor´ eor´ıa de Flory–Huggins Flory–Hu ggins en la ingenier´ıa ıa de. . . pol´ımeros. ımeros. Eduardo Macario Moctezuma N.
nados nados disolve disolvente ntes, s, depende depende de las condicion condiciones es termodin´amicas amicas de cada sistema. Por ende, se requiere un buen marco te´orico orico para describir tales condiciones termodin´amicas. amicas. As´ As´ı, aunque se han desarrollado numerosos enfoques para aquellas soluciones en las cuales el soluto est´a conform con formado ado por p or mol´ mo l´eculas ecul as de pol´ımero, ımero, en realidad, la teor´ teor´ıa de Flory–Huggins Flory–Hu ggins consti constituy tuyee uno de los pilares pilares sobre los que se ha construido la interpretaci´ interpretaci´on on termodin´ termodin´ amica amica vigenvigente del comportamiento comportamiento de las soluciones polim´eriericas, de manera que se ha vuelto punto de partida para los estudiosos en el tema, y por ello la estudiaremos en este trabajo. De la teor´ teor´ıa de soluciones, sabemos que su comportamiento tamiento est´a determ determina inado do por las variaci ariacione oness de energ´ energ´ıa libre del sistema. Adem´as, as, el cambio en la energ energ´ıa libre libre de Gibbs Gibbs para el proceso proceso de mezclamezclado est´a dada por ∆Gm = ∆H m − T ∆ T ∆S m
(1)
donde H donde H,, S , se s e refier r efieren en a la l a entalp´ e ntalp´ıa ıa y entrop´ e ntrop´ıa, ıa, y T es la temperatura del sistema. En particular, sabemos que si ∆G ∆Gm < 0 entonces la disoluci´on on ser´a un proceso favorable y espont´aneo. aneo. A continuaci´ continuaci´on, on, veamos el enfoque FH para describir matem´aticamenaticamente estos esto s t´erminos. ermi nos. Teor´ ıa ıa de Flory– Flo ry–Hug Huggin ginss
Las soluciones sol uciones polim´ericas, ericas, tienen un comportamiencomp ortamiento muy diferente al de las soluciones ideales, b´asicamente por causa de las fuertes interacciones entre soluto polim´ erico erico y disolvente; disolvente; contrario a lo que se presenta pres enta entre mol´eculas ecul as m´ m ´as as peque˜ nas. nas. En contraste con lo que ocurre en una soluci´on polim´ pol im´erica, eric a, las interacciones entre mol´eculas eculas peque˜ peque ˜nas, nas, son de menor intensidad. La teor´ teor´ıa de Flory—Huggins modela las interacciones interacciones pol´ pol´ımero–disolvente ımero–disolvente conduciendo a resultados cualitativos importantes, como lo es el comportamiento de la solubilidad y las condiciones para que ´esta esta se presente, dando lugar a un me jor entendimiento entendimiento de –por ejemplo–, los procesos de separaci´ on on a nivel industrial. Recordemos que, para el an´alisis alisis termodin´amico amico de las soluciones polim´ ericas, ericas, se requiere el c´alculo de la entrop´ entrop´ıa y la entalp´ entalp´ıa de mezclado, a partir de cuyos resultados se deduce la energ´ energ´ıa libre asociada aso ciada y luego la actividad de cada especie. esp ecie. A partir de ah´ ah´ı, podemos vincular la teor´ teor´ıa con propiedades medibles. Utilizando un modelo de malla, veamos a continuaci´on o n c´omo omo se hizo esto, en el enfoque pionero de FH.
55
Entrop´ıa de mezcla Flory y Huggins (en las referencias 2-–4, 7 y 8), derivaron de forma independiente una expresi´on para la entrop´ entrop´ıa de mezclado en soluciones solucione s polim´eriericas, ba jo la idea general de que la entrop´ entrop´ıa est´a asociada al acomodo de mol´ eculas eculas en una malla o red tridimensional, y concretamente, al c´alculo de configuraciones p osibles que las mol´ eculas eculas pueden ocupar en la red. FH partieron de los siguientes supuestos: a) el sistema tema est´ est´a repres represen entad tadoo por un modelo modelo de mal malla la o red; red; b) cada cada siti sitioo en la ma mall llaa tien tienee un voluolumen igual al de las mol mol´´eculas eculas de disolven disolvente; te; c) la mol´ecula ecul a de pol po l´ımero ımer o se compone comp one de m unidades monom´ nom´ericas ericas o segmento segmentos, s, cada uno definido definido a partir tir de los los vol´ umenes umenes de la mol´ecula ecula de pol´ımero ımero (v p ) y de la mol mol´´ecula ecula de disolven disolvente te ( vs ), como sigue: m = v p /vs ; d) cada segmento tiene un n´umeumero z de sitios vecinos disponibles (a z se le llama n´ umero umero de coordinaci´on); on); e) las interacciones moleculares no contribuyen a la entrop´ entrop´ıa pero s´ı al calor de mezclado. Como ya dijimos para el c´alculo alculo de la entrop´ entrop´ıa, FH hacen uso de un modelo de malla o reticular para representar la disoluci´on on de pol´ımero ımero en disolvente (v´ease ease la Figura 1). Adem´as, as, se propone que la disoluci´ on se lleva a cabo mediante la ocupaci´on on o n al azar de los sitios o lugares disponibles en la malla, por parte de las mol´eculas eculas de disolvente y de pol´ımeımero. A partir de esta ocupaci´on on o llenado, se calcula Ω, es decir, el n´umero umero de formas fo rmas en que las mol´ mo l´ecueculas de pol´ımero ımero podr po dr´´ıan llenar la malla. Esta Ω, permite determinar determina r la entrop´ entrop´ıa del sistema, a trav´es es de la ley de Boltzmann. Tengamos en cuenta que, por la ley de Boltzmann, el n´umero umero de estados microsc´opiopicos (o microestados) de un sistema, asociados a un estado macrosc´opico opico particular, est´a directamente ligado con la entrop´ entrop´ıa de dicho sistema. sis tema. Matem´ M atem´aticamente, esta ley se escribe S = k = k ln Ω
(2)
donde k es una consta constant ntee unive universa rsall tal que k = −23 1.380662 × 10 J/K. Con un poco de trabajo adicional, Flory y Huggins obtuvieron la siguiente expresi´ on para el mezclado en una soluci´on polim´ on po lim´erieri ca:
∆S m =
−k (n1 ln φ1 +
n2 ln φ2 )
(3)
56
ContactoS ContactoS 68, 54–62 (2008)
donde z donde z es es el n´ umero umero de coordinaci´on on (n´ umero umero de sitios que son vecinos m´as as cercanos a una celda en un espacio tridimensional), ∆w ∆ w12 se refiere a la energ´ıa ıa de contacto asociada asocia da con la mezcla de pol po l´ımero y disolvente, y solvente, y es es el cociente co ciente del de l volumen de una un a mol´ecuecula de pol´ımero ımero entre el volumen de una mol´ecula ecula de disolvente (y (y = m); y finalmente, T T es la temperatura. Con la introducci´on on de este par´ametro, amet ro, la entalp ental p´ıa de la soluci´on on queda: ∆H m = kTχn1 φ2 Figura 1. Llenado aleatorio de sitios en una malla o red.
donde φ1
=
n n + n m nm 1
1
φ2
=
2
2
n1 + n2 m
(4)
∆S m =
−k (n1 ln x1 +
n2 ln x2 )
donde χ es el denomi denominad nadoo par´ par´ametro ametro de interacinteracci´ on on de Flory-–Huggins. Flory-–Hu ggins. En cuanto al significado signi ficado f´ısiısico concreto de χ de χ,, ´este este se interpreta i nterpreta como una cantidad adimensional ad imensional que caracteriza car acteriza la energ ene rg´´ıa de interacci´ on on de mol´eculas ecul as de pol po l´ımero ımer o con c on disolvente diso lvente.. Por ende, ende , el t´ermino ermi no kT k T χ modela los efectos de mezclado: energ´ ene rg´ıa ıa del
son las fracciones–volumen de los componentes de la disoluci´on on (n1 y n2 son las moles de disolvente y de pol´ pol´ımero, respectivamente). respectivamente). En contraste, la entrop´ trop´ıa de mezclado en una soluci´on on ideal viene dada por (5)
donde x donde x1 y x y x 2 son las fracciones–molares de las componentes de la disoluci´on. on. Podemos concluir entonces que en solucione solucioness polim´ polim´ericas, ericas, no debe usarse la fracci´on on molar como variable independiente para describir descr ibir la entrop e ntrop´´ıa de mezclado. mez clado. En su lugar, lugar , debe utilizarse la fracci´on on volumen. Entalp´ıa de mezcla Para determinar la entalp´ entalp´ıa del sistema, Flory propuso la idea id ea de energ ener g´ıa de contacto conta cto entre las mol´ecueculas como causa de la entalp´ entalp´ıa de mezcla. As´ As´ı, optando por p or el enfoque de van Laar de energ´ıa ıa de contacto entre soluto y disolvente, utiliz´o espec´ esp ec´ıficamenıfica mente la expres expresi´ i´ on on de aqu´ aqu´el el para el calor calor de mezclamezclado: ∆H m = z∆ z ∆w12 yn 1 φ2 (6)
(8)
kT χ =
disolvente en soluci´on on
energ´ ene rg´ıa ıa del −
disolvente en su estado puro
(9)
As´ As´ı, con las ecuaciones (2) y (7), ya podemos completar el c´alculo alculo de la energ´ energ´ıa libre de la soluci´on. on. ´ Energ´ıa libre para soluciones polimericas ericas Con los anteriores resultados podemos obtener una expresi´ on on para la energ´ energ´ıa libre en el sistema. Sustituyendo (2) y (7) en (1), concluimos que ∆Gm = kT [( kT [(n n1 ln φ1 + n2 ln φ2 ) + χn1 φ2 ]
(10)
donde (n (n1 ln φ1 + n2 ln φ2 ) es un t´ ermino ermino combinatorio y χn y χn 1 φ2 es un t´ ermino ermino de contacto. contac to. Debe De be destacarse que el par´ametro ametro χ no se conoce a priori y se obtiene experimentalmente. Sin embargo, a partir de la energ´ energ´ıa libre podemos deducir las dem´as as propiedades termodin´amicas, amicas, como hacemos en el siguiente apartado.
En este este punto, punto, Flory Flory defini´ defini´ o el siguie siguient ntee par´ par´ameametro de interacci´on on (caracterizado por las interacciones entre pares de segmentos polim´ p olim´ ericos, ericos, entre pares de mol´ eculas eculas de disolvente disolvente y entre un segmento de pol´ımero ımero y una mol´ecula ecula de disolvente):
Actividad y potencial qu´ımico Para vincular las expresiones derivadas por Flory y Huggins con cantidades medibles, es pertinente diferenciar la ecuaci´on on de d e la energ´ıa ıa libre libr e de mezclado con respecto a n a n 1 para obtener el potencial potenc ial qu´ qu´ımico y luego aprovechar la relaci´on entre ent re poten p otencial cial qu´ımiımico y actividad. As´ As´ı, partir de (10) nos conduce a la siguiente expresi´on on para la actividad del disolvente
z ∆w12 y χ = kT
(7)
ln a1 = ln(1 − φ2 ) + 1 −
1 m
φ2 + χφ22 (11)
La Teor´ eor´ıa de Flory–Huggins Flory–Hu ggins en la ingenier´ıa ıa de. . . pol´ımeros. ımeros. Eduardo Macario Moctezuma N.
y an´alogamente alogamente para el pol´ımero: ımero: ln a2 = ln φ2 − (m − 1)(1 − φ2 ) + χm(1 χm(1 − φ2 )2 (12) ´ lculo de χ Calculo a de χ Hildebrand y Scott (1950, pp. 75-–92) mostraron que χ puede relacionarse con los par´ametros ametros de solubilidad, δ dad, δ ,, a trav´ tra v´es es de: χ =
V 1 (δ (δ 1 − δ 2 )2 RT
(13)
donde V donde V 1 es el volumen molar del disolvente y R es R es la constante universal universal de los gases. En realidad, la ecuaci´ on (13) no ofrece una buena descripci´on cuantion tativa de soluciones solucion es polim´ericas, ericas, pero s´ı provee un buen criterio para el an´alisis alisis cualitativo de la solubilidad de tales sistemas (como veremos en la secci´ on on de aplicaciones). M´as as bien, los resultados experimentales, sugieren que χ = χ = χs +
V 1 (δ (δ 1 − δ 2 )2 RT
(14)
donde dond e el t´ermino ermi no χ s mide la contribuci´on on entr´opica opica y suele asociarse al n´umero umero de coordinaci´on, z on, z,, de la siguiente manera: χs = 1/z (a χs se le suele llamar constante de malla). Valores usuales son 0.3 < χs < 0 < 0..4 o bien χ bien χ s = 0.34. Por otro lado y para fines de completar el c´alculo, alculo, se destaca que los par´ametros ametros de solubilidad, δ , se encuentran tabulados para muy diversas sustancias. Como puede verse, las ecuaciones (13) y (14) permiten considerar el efecto que la temperatura tiene en χ en χ.. Por otro lado, aunque en la derivaci´on original de Flory, χ Flory, χ es es constante respecto a la composici´ on, en realidad, la evidencia experimental muestra on, que χ que χ no no s´olo olo var´ıa ıa con la temperatu tem peratura ra sino si no que qu e tambi´en en var´ var´ıa ante cambios cambi os en la concentra conc entraci´ ci´on de la soluci´on. on. En este sentido, existen manuales con una lista muy amplia de valores y ecuaciones de variaci´ on on de χ de χ (po (porr ejemplo, v´ease ease Brandrup Brandru p et al., 1989, p. 173-–182). Aplicaciones ingenieriles
Tras la breve revisi´on on de la teor´ teor´ıa, en esta secci´on on examinaremos varios varios sistemas pol´ pol´ımero–disolvente ımero–disolvente de inter´ int er´es es en ingen in genier´ ier´ıa. ´ n de aditivos 1. Seleccion o En la industria industria de los pl´asticos asticos es frecuente el uso de aditivos para optimizar o mejorar determinadas propiedades de los productos (por ejemplo, aumentar la
57
resistencia, impedir la depositaci´on de polvos, me jorar la impermeabilidad, retardar la degradaci´on, on, etc.). Algunos problemas asociados a la mejora por agregaci´ on de aditivos, es que estas substancias pueon den migrar, y que a corto plazo puede haber hab er p´erdierdida o desincorporaci´on on por medios f´ısicos, entre otros inconvenientes. El an´alisis alisis termodin´amico amico del sistema pol p ol´´ımero–aditivo ımero–aditi vo puede proporcionar propor cionar una primeprime ra gu´ gu´ıa en la selecci´on de aditivos. Un criterio criterio inicial para distinguir distinguir solubilid solubilidad, ad, hace uso de la ecuaci´on on (13) y afirma que la solubilidad se presentar´a si χ ≤ 0.5. En caso caso contra contrario rio (χ > 0. 0 .5), no habr´a miscibilidad y el aditivo presentar´ a alguno de los problemas que mencionamos en l´ıneas ınea s anteriore anter iores. s. Haciendo uso de los datos de la Tabla 1 y considerando temperatura ambiente, indique cu´ales ales de las posibles combinaciones c ombinaciones pol´ımero–aditivo ımero–adit ivo deben deb en seleccioselec cionarse (y cu´ales ales no). Tabla 1. Par´ ametros de solubilidad, ametros tomados de las Referencias 5, 9 y 11. Sustancia PBMA Poli( oli(me meta tacr cril ilat ato o de n–bu n–buti tilo lo)) PEMA Pol Poli( i(me mettacri acrila lato to de etil tilo) PMMA P Pol oli( i(me meta tacr cril ilat ato o de meti metilo lo)) PVAC Po Poli(acetato de vinilo) TiO2 Bi´ oxido de Titanio Al Aluminio ´ xido de hierro III Fe2 O3 O
δ en (cal/cm3 ) 8.75 8.75 8.95 8.95 9.50 9.50 9.40 11.00 9.00 5.50
1 2
Los pesos moleculares y densidades de los aditivos son: PM = 80 g/gmol y ρ = ρ = 3.3g 3g//cm3 para el TiO2 , PM = 27 g/gmol y ρ = 2.7g 7g//cm3 para el Al, y por u ultimo, ´ ltimo, PM = 159.7 g/gmol y ρ = 5.2g 2g//cm3 para el Fe2 O3 . ´n Solucion o Veamos el par PBMA–TiO PBMA–TiO2 . El volumen molar del aditivo aditivo es V 1 = P M 1 /ρ1 = 24 24..24cm3 /gmol, y el par´ ametro ametro de interacci´on on se calcula a partir de (13), como sigue: V 1 (δ (δ 1 − δ 2 )2 RT 24 24..24 = (11 − 8.75)2 (1. (1.9872)298 = 0.207
χ =
donde R donde R = = 1.9872cal/ 9872cal/gmol gmolK K es la constante constante univeruniversal de los gases. As´ As´ı, como co mo χ χ = = 0.207 ≤ 0.5, entonces
58
ContactoS ContactoS 68, 54–62 (2008)
el aditivo se mantendr´a unido u nido al pol p ol´´ımero, pues presentan solubilidad. Con el mismo procedimiento, elaboramos la Tabla 2. Tabla 2. Par´ ametros ametros de interacci´on, on, χ, del problema 1. Aditivo P PB BMA TiO2 0. 0.207 Al 0.085 Fe2 O3 0. 0.548
Pol´ Pol´ımer ım ero o PEMA PMMA 0.172 0.092 0.071 0.038 0.617 0.830
PVAC 0.105 0.043 0.789
Podemos ver que los tres pares asociados a Fe 2 O3 , verifican χ > 0.5, por tanto no habr´a miscibilidad con ninguno de los pol´ımeros ımeros de inter´es. es. Por otro lado, como todos los pares vinculados a TiO 2 y Al, presentan χ ≤ 0.5, entonces podemos concluir que estos aditivos s´ı tendr´an an un comportamiento termodin´ amico propic amico propicio, io, puesto puesto que son solubl solubles es en cada uno de los cuatro cuatro pol´ pol´ımeros. ımeros. De hecho, hecho, la solubilidad se favorece en la medida que χ es m´as a s peque˜ no. no. As´ı, ı, seg´ se g´un un este criterio, el aluminio es el aditivo con el mejor grado de adopci´on en los pol´ımeros ımer os de inter´ i nter´es, es, pues da lugar a par´ p ar´ametros ametros χ menores que los correspondientes al TiO 2 y Fe2 O3 . En particular, la adici´on on de aluminio a poli(metacrilato de metilo) es la combinaci´on m´as as favorecida desde el punto de vista termodin´amico*. amico*. Tabla 3. Caracterizaci´on on de aditivos a ditivos en pol´ımeros ımeros Aditivo PBMA TiO2 Mi Miscible Al Misc. Fe2 O3 Inmi Inmisc scib ible le
Pol´ Pol´ımer ım ero o PEMA PMMA Misc. Misc. Misc. Misc.* Inmi Inmisc sc.. Inmi Inmisc sc..
PVAC Misc. Misc. Inmi Inmisc sc..
´ n sanitaria 2. Legislacion o Una aplicaci´on on interesante interesante de la teor´ teor´ıa de soluciones polim´ericas ericas es la descrita p or Prausnitz et al. (2000, p. 408, problema propuesto), en el marco del cumplimiento cumplimiento de la normatividad sanitaria para ciertos productos pro ductos de inter´es es comercial. comercia l. Resulta Result a que el p oli(acetato de vinilo) –en adelante PVAC–, es utilizado como pegamento para madera y otros materiales; por ejemplo, se suele emplear como adhesivo de encuadernaci´ on, de bolsas de papel, de cartones de leon, che, sobres, etc. Pues bien, en la ´ultima etapa del proceso de producci´on on de pel´ pel´ıculas de PVAC, PVAC, trazas de alcohol a lcohol isoprop´ılico ılico acompa˜nan el producto final. Por razones sanitarias, las autoridades gubernamentales exigen exig en un contenido de d e alcohol alcoho l en la pel´ıcuıcu−4 la que verifique: φ 1 < 10 < 10 .
Para eliminar el alcohol residual, se ha propuesto evaporarlo a 125o C. Sin embargo, se ignora cu´al es el valor de la baja presi´on que debe mantenerse en el evaporador para alcanzar la pureza requerida en la pel pe l´ıcula. ıcu la. En el marco de la teor´ teor´ıa de FH, ¿podr´ ¿podr´ıa usted dar una respuesta al problema de la presi´on de operaci´ on del evaporador? Como datos, tiene usted que a on 125o C, χ C, χ = = 0.44 y la presi´on on de vapor del alcohol isoprop´ılico ılico en estado es tado puro es 4.49 4 .49 bar. b ar. ´n Solucion o En una mezcla mezcla binaria binaria de pol´ pol´ımero ımero y disolven disolvente te a presiones ordinarias, ordinaria s, el pol´ımero ımero no es vol´atil, atil, s´olo olo lo es disolvente. Por ende, la presi´on on total del sistema es igual a la presi´on on parcial del disolvente. Por otro lado, el enunciado del problema sugiere una “baja presi´ on on en el evaporador”, as´ as´ı que podemos suponer que la fase gaseosa del sistema se comporta como gas ideal, por tanto P 1 a1 = 0 (15) P 1 donde P donde P 1 es la presi´on on parcial y P y P 10 es la presi´on on de vapor en estado puro, ambos del disolvente. Por la teor´ teor´ıa de FH, sabemos que
1 m
ln a1 = ln(1 − φ2 ) + 1 −
φ2 + χφ22 (16)
y si trabajamos con los datos que ya tenemos, resolveremos el problema (v´ease ease el inciso b) del siguiente problema). Sin embargo hay una manera m´as f´acil acil de resoluci´on on por las condiciones del sistema. En disoluciones de pol po l´ımeros, es conveniente definir d efinir el coeficiente de actividad Γ 1 = a1 /φ1 donde don de el sub´ındiındi ce 1 se refiere al componente m´as vol´atil. atil. Utilizando (16), tenemos que ln Γ1 =
1−
1 m
φ2 + χφ22
(17)
tal que si φ2 → 1 como en nuestro caso (pues solo hay trazas traz as de d e alcohol alc ohol isoprop´ılico), ılico), entonces ∞
ln Γ1 = 1 −
1 m
φ2 + χ
(18)
donde dond e el e l sup s uper´ er´ındice ındi ce ∞ denota que φ que φ 1 → 0 (es decir, φ2 → 1). Adem´as, as, como m ≫ 1 para pol´ımeros ımeros con alto peso molecular, (18) se simplifica a: ln Γ∞ 1 = 1 + χ
(19)
La Teor´ eor´ıa de Flory–Huggins Flory–Hu ggins en la ingenier´ıa ıa de. . . pol´ımeros. ımeros. Eduardo Macario Moctezuma N.
59
la fase gas, podemos decir que P 1 = a1 P 10 donde P 10 (25◦C ) = 0.0492 049211 bar. bar . As´ı, ı, P 1 ≈ 0.0437 bar. Como en la fase gaseosa s´olo olo hay disolve disolvente nte,, se tiene que la presi´on on total del sistema est´a dada por la presi´ presi´on on de vapor del toluen toluenoo en la disolu disoluci´ ci´on: on: P = P 1 = 0.0437 bar. b) ¿Podr´ıa ıa obtener o btener una expresi´on para el c´alculo alculo de la presi´on on total de la disoluci´on, on, a lo largo de todo el intervalo intervalo de concentraci´ on? Procure que su resultado on? est´ es t´e s´olo olo en funci´on on de las fracci fraccione oness volum volumen en de pol po l´ımero ımer o y disol d isolvente. vente. Figura 2. Disoluci´ on on Tolueno-–Poliestireno olueno-–Poliestireno a temperatura ambiente.
As´ As´ı, la presi´on o n de opera operaci ci´´on on del evapor evaporado adorr deber´a corresponder a la de la presi´on on parcial del alcohol isopoprop´ isopo prop´ılico, ılico, dada por: P 1 P 1
0 = a1 P 10 = φ 1 Γ∞ 1 P 1 = φ1 (1 + χ1 )P 10
(20) (21)
Finalmente, por los datos del problema, concluimos que la baja presi´on on en el evaporador, deber´a ser de P 1 = (1×10−4 )(1+0. )(1+0.44)(4. 44)(4.49 bar) = 6. 6 .465×10−4 bar ´ n de vapor y seguridad industrial 3. Presion o En una planta piloto, se tiene un recipiente que contiene una disoluci´on on compues compuesta ta por un 50 % en masa de tolueno y el resto de poliestireno, a una temperatura de 25o C. La densidad del poliestireno es 903 kg/m3 y la del tolueno es de 867 kg/m 3 . a) ¿Cu´al al es la presi´on on total que ejerce la disoluci´on? on? ´n Solucion o En primer lugar, calculemos las fracciones volumen de los elementos de la disoluci´on: on: φ1 φ2
=
w1 /ρ1 = 0. 0 .51 w1 /ρ1 + w2 /ρ2
´n Solucion o Sustituyendo χ = χ0 + χ 1 φ2 + χ 2 φ22 en la ecuaecuaci´ on (11) y luego despejando a1 para a su vez suson tituirlo en P 1 = a1 P 10 , obtenemos la siguiente expresi´ on, on, v´alida alida para todo el intervalo de concentraci´ on: on: P 1 = φ = φ 1 P 10 exp φ2 + χ0 φ22 + χ1 φ32 + χ2 φ42 = P = P (23) c) Para evitar una fuga, las consideraciones de seguridad establecen que a 25 o C, la mayor presi´on on que puede soportar el recipiente es de 0.0485 bar. Determine la masa m´axima axima de tolueno y de poliestireno que puede guardarse en este contenedor particularmente fr´agil. agil. Suponga que el volumen del recipiente es de 0.2 m 3 .
´n Solucion o Si P max = 0.0485 bar = P 1max , entonces amax = 1 max 0 0 0 P 1 /P 1 siempre siempre que f 1 /f 1 ≈ P 1 /P 1 . Por ende, sustituyendo amax = 0.985 98555 en (11) (11) y resolv resolvien ien-1 do para φ para φ 2 , encontramos que φ 2 = 0.215, por tanto: tanto: φmax = 0. 0 .785. Por relaciones elementales entre frac1 ci´ on on masa y fracci´on on volumen, tenemos el siguiente sistema de ecuaciones lineales: w1 + w2 = 1
(22)
= 1 − φ1 = 0.49
w1 − w2
ρ φ 1
1
ρ2 φ2
(24) = 0
De la literatura (Brandrup et al., 1989, p. 175) se tiene que para el sistema en cuesti´on y a una temperatura de 25o C, χ C, χ var var´´ıa con la concentra conc entraci´ ci´on on de acuerdo con la siguiente expresi´on: on:
Resolviendo (24), vemos que w1 = 0.7781 y w2 = 0.2219. Con esto, podemos calcular la densidad de la soluci´on, ρ on, ρ,, pues :
χ = χ = χ0 + χ1 φ2 + χ2 φ22
1 w1 w 2 = + ρ ρ1 ρ2
tal que χ0 = 0.431 431,, χ1 = −0.311 y χ2 = −0.036. As´ı, ı, con co n φ2 = 0.49 entonces entonces χ = 0.27. UtiliUtilizando zando (11), (11), calcul calculemo emoss la activi actividad dad del disolv disolvenente (tolue (tolueno) no) en la fase fase gaseos gaseosa: a: a1 = 0.8882. Co0 0 mo a1 = f 1 /f 1 ≈ P 1 /P 1 si asumimos idealidad en
(25)
As´ As´ı, tenemos que la masa m´axima axima de tolueno que puede puede guarda guardarse rse en este este conten contenedo edor, r, es de m1 = w1 m = w1 ρV = (0. (0.7781)(874. 7781)(874.7)(0. 7)(0.2) = 136. 136.1 kg. An´ alogamente, la masa m´axima alogamente, axima de poliestireno es
60
ContactoS ContactoS 68, 54–62 (2008)
de m2 = w2 m = w2 ρV = (0. (0.2219)(874. 2219)(874.7)(0. 7)(0.2) = 38 38..8 kg. ´ n de miscibilidad en mezclas 4. Identificacion o ´ polimericas ericas (Polyblends ) Es bien sabido que la mezcla mezcla f´ısica ısica de dos o m´as as pol´ımeros ımeros qu´ qu´ımicamente diferentes suele dar lugar a ma mater terial iales es con propie propiedad dades es superio superiores res a las de los pol´ pol´ımeros ımeros por separado, separado, lo cual genera un mayor valor valor agregado. agregado. Por esta raz´on, on, las mezclas o aleaciones polim´ericas ericas ( polyblends o polymer alloys , en ingl´ ingl´es) es) son objeto de considerab considerables les esfuerzos esfuerzos de investig investigaci´ aci´ on on cient cient´ıfica y desarroll desarrolloo comercial comercial de productos. El punto de inter´es es para par a este trabajo, traba jo, est´ est ´a en que las mezclas polim´ericas ericas no solo pueden formar una ´uniunica fase, sino que tambi´ tambi´en en pueden dar lugar a dos o m´ as fases distintas, consistentes principalmente de as los componentes comp onentes individuale in dividuales. s. As´ı, ı, lo que qu e interesa es determinar las condiciones bajo las cuales se presenta la miscibilidad y mantener el sistema en ese estado. Sin embargo, m´as as que una miscibilidad permanente, lo que a menudo ocurre es que se observan las denominadas “ventanas de miscibilidad”, es decir, un intervalo de composici´on on en e n el cual cua l un pol´ımeımero es miscible en el otro. En este este contex contexto, to, una pregun pregunta ta import importan ante te es: ¿c´ omo saber si la mezcla omo mezcla f´ısica ısica de dos pol´ pol´ımeros ser´ ser´a esta establ blee (es (es deci decir, r, que que no habr habr´´a separa separa-ci´ on de fases)? Bien, el efecto de la estructura moon lecular en la miscibilidad miscibilida d pol´ımero–pol ımero–p ol´´ımero ha sido abordado por Paul y Barlow en un trabajo pionero (Ref. 14). Ellos derivaron derivaron una expresi´ expresi´on o n para la ental entalp p´ıa de mezcla mezclado do de sistem sistemas as homohomopol po l´ımero– ıme ro–hom homop opol´ ol´ımero ıme ro,, copol´ cop ol´ımero ıme ro–ho –homop mopol´ ol´ımeıme ro y copol´ımero–copol ımero–cop ol´´ımero. La expresi´on obteniobtenida est´a en funci´on on de par´ametros ametros de interacci´on on binaria de Flory y de la composici´on on de los elementos de la mezcla, de manera que para el caso homopol´ımero–cop ımero –copol´ ol´ımero ımer o la ecuaci´ ecua ci´on del calor de mezclado resulta ser: ¯ A φB ∆H m = V = V Bφ
(26)
tal que V = V A + V B (27) ′ ′ ′ ′ ¯ = B13 φ1 + B23 φ2 − B12 φ1 φ2 (28) B ¯ es un par´amedonde V donde V es es el volumen de la mezcla, B ametro de interacci´on on global, φA y φB son las fracciones volumen volume n del de l copol´ cop ol´ımero ımer o y del d el homop h omopol´ ol´ımero ımer o resre s-
Mon´omero omero 1
Mon´omero omero 2
φ1′
φ2′
A: Copol´ Cop ol´ımero ıme ro φA Mezcla Mezc la polim´ pol im´erica eric a φ1 + φ2 + φ3 = 1
Mon´omero omero 3 φ1′ + φ2′ = 1
φB B: Homo H omopo poll´ımero ıme ro φ1 + φ2 = φ = φ A
φA + φB = 1 φ3 = φ = φ B
φ1′ =
φ1 φA
φ2′ =
φ2 φA
Figura 3. Representaci´on on esquem´ atica atica de la mezcla homop hom opol´ ol´ımero ım ero–c –cop opol´ ol´ımero ıme ro..
pectivamente (φ ( φA + φ B = 1), Bij son los par´ameametros de interacci´on on entre cada par de mon´omeros, omeros, y por ultimo, u ´ltimo, φ1′ y φ 2′ son fracciones volumen auxiliares, tales que se verifican las relaciones de la Figura 3. En el l´ımite ımite de alto peso p eso molecular, Paul y Barlow encontraron que para que exista miscibilidad se requiere cumplir: ∆H m < 0 < 0
⇐⇒
¯ <0 B
(29)
lo cual significa s ignifica que q ue la miscibilidad misc ibilidad de d e pol´ımeros ımeros (con (c on su alto peso molecular caracter´ caracter´ıstico), ´unicamente unicamente puede asegurars as egurarsee cuando cuand o el mezclado mez clado es exot´ermico. ermico. Con la informaci´on on anterior y utilizando los datos de la Tabla abla 4, consid considere ere todas todas las mezcl mezclas as homohomopol´ımero—copol ımero—cop ol´´ımero posibles posibl es e identifique alguna “venta “ventana na de miscibil miscibilidad” idad” en el l´ımite ımite de alto peso molecular. Tabla 4. Par´ ametros amet ros energ´ ener g´eticos etico s de interacci´on on molecular. Par monom mo nom´ ´erico eri co Estireno—Acrilonitrilo Estir stiren enoo-–M –Meetac tacrilat ilato o de metil tilo Acri Acrilo loni nitr tril ilo— o—Me Meta tacr cril ilat ato o de meti metilo lo
Bij en cal/cm3
6.740 0.18 0.181 1 4.11 4.110 0
´n Solucion o Con los tres mon´omeros omeros indicados y sustituyendo los par´ ametros ametros de interacci´on on correspondientes en (28), podemos pod emos formar for mar tres combinaciones homopol´ımero– ımero– cop co pol´ımer ımero: o: La 1ra. combinaci´on, on, es la de la mezcla Poliestireno / Poli(Metacrilato de Metilo–co–Acrilonitrilo), con
La Teor´ eor´ıa de Flory–Huggins Flory–Hu ggins en la ingenier´ıa ıa de. . . pol´ımeros. ımeros. Eduardo Macario Moctezuma N.
61
un par´ ametro ametro de interacci´on on global dado por ¯ = 0.181(φ B 181(φ1′ )2 + 2. 2.449 449φ φ1′ + 4. 4 .11
(30)
La 2da. combinaci´on, on, corresponde a la mezcla Poliacri liacrilon lonit itril riloo / Poli( Poli(Met Metacr acrila ilato to de Metilo Metilo–co –co–– Estireno), y tiene el siguiente par´ametro ametro de interacci´ on: on: ¯ = 4.11(φ B 11(φ1′ )2 + 2. 2.449 449φ φ1′ + 0. 0 .181 (31) La 3ra. combinaci´on, on, es la mezcla Poli(Metacrilato de Metilo) Metilo) / Poli(Est Poli(Estiren ireno–co– o–co–Acri Acrilonit lonitrilo rilo), ), cuyo cuyo par´ ametro ametro de interacci´on on global es: ¯ = 6.74(φ B 74(φ1′ )2 − 10 10..669 669φ φ1′ + 4. 4 .11
(32)
Para determinar la existencia de “ventanas de miscibilidad” asociadas a las expresiones (30)–(32), po¯ y buscar demos tabular tabular φ1′ vs. B buscar region regiones es donde donde ¯ B < 0. La represent representaci´ aci´ on on gr´ afica afica de este procedimiento se muestra en las Figuras 4–7.
Figura 5. Efecto de la l a composici´ compos ici´on on del copol´ımero ımero en ¯ el par´ ametro ametro global B para la mezcla
PAN/P(MMA–co–S).
Figura 6. Efecto de la composici´on on del copol´ımero ımero en ¯ el par´ ametro ametro global B para la mezcla
PMMA/P(S–co–AN).
Figura 4. Efecto de la composici´on on del copol´ cop ol´ımero ımer o ¯ en el par´ ametro ametro global B para la mezcla PS/P(MMA–co–AN). Limitaciones
De las Figuras 4 y 5 se observa que las dos primeras mezclas polim´ ericas ericas no presentan un cambio ent´ alpico alpico de d e tipo tip o exot´ermico, ermico, de manera que son sistemas inmiscibles. Por el contrario, en la mezcla PMMA / P(S–co–AN) habr´a una regi´on on de miscibilidad dad bien bien defin definid ida, a, debi debido do a que que el par´ par´ametro ametro de ¯ interacci´on on glo global bal B alcanza alcanza valores valores negativ negativos os en 0.67 ≤ φ1′ ≤ 0.91 (v´ease ease Figura 7). As´ As´ı, por el criterio dado en la expresi´on on (29), existe en ese intervalo una regi´on on o ventana de miscibilidad, con un m´aximo aximo ¯ = B ¯ de miscibilidad en φ en φ1′ = 0.79 (ya que qu e ah´ a h´ı B min ).
La validez alidez (y limita limitaci´ ci´on) on) de la teor´ teor´ıa de Flory– Flory– Huggins se ha comprobado v´ıa ıa la determinaci´on on experimen perimental tal de medida medidass de presi´ presi´ on o n de vapor, apor, prepresi´on on osm´otica, otica, diagramas diagramas de fase, fase, etc. Sin embarembargo, si bien las aplicaci aplicaciones ones de la teor´ teor´ıa de Flory– Flory– Huggins son abundantes y muy variadas, existen casos importantes en los que su eficacia es restringida, en concreto: a) cuando la soluci´on es muy diluida (menor al 1 % en peso), b) cuando cuando el pol´ pol´ımero es r´ıgido, c) cuando se consideran sistemas polares, d) para c´alculos alculos de equilibrio de fases a presiones altas.
Tabla 5. Resultados del Problema 4. # Mezcla polim´erica Conclusion o´n 1 PS/P(MMA–co–AN) Inmiscible 2 PAN/P PAN/P (MMA–co–S) Inmiscible 3 PMMA/P(S–co–AN) Parcialmente Parcialmente miscible
A favor favor del formalismo formalismo Flory–Hugg Flory–Huggins, ins, se asevera asevera que “esta teor´ teor´ıa es relativamente relativamente simple y contiene casi todas las caracter´ caracter´ısticas esenciales que diferencian las soluciones de mol´ eculas eculas grandes de las
62
Figura 7. Localizaci´ on de la “ventana de miscibilidad” on en la mezcla PMMA/P(S–co–AN).
que contienen solamente mol´ eculas eculas de tama˜no no normal” (Prausnitz et al., 2000, p. 348). Conclusi´ on on
La carencia de modelos termodin´amicos amicos eficaces para soluciones y mezclas polim´ ericas ericas puede resultar en consecuencias indeseables en la operaci´on de las plantas industriales. Por tanto, para los ingenieros de proceso, es importante contar con modelos que incorporen corporen la mejor mejor descripci´ descripci´on on fisicoqu fisico qu´´ımica posible. posib le. As´ As´ı, un buen b uen modelo mo delo ayuda a los ingenieros i ngenieros a identificar las condiciones que dar´an an lugar o no a un proceso industrial eficiente. Para sistemas macromoleculares como los vistos en este trabajo, el cl´asico asico modelo de malla de Flory– Huggins, Huggins, captura captura dos caracter caracter´´ısticas ısticas clave clave de los pol´ pol´ımeros: el efecto del tama˜no no molecular sobre la entrop´ entrop´ıa de mezcla y el efecto de la interacci´ interacci´on molecular sobre el t´ ermino ermino de entalp´ entalp´ıa. Si bien el modelo FH tiene limitaciones importantes, su alcance descriptivo cualitativo lo sit´ua ua como un pilar y un antecedente de revisi´on on obligada en la frontera actual del conocimiento termodin´amico amico de soluci solucione oness y mezmezclas polim´ pol im´ericas. eric as. Agradecimientos
A la Academia Mexicana de Ciencias por su apoyo para participar en el X Verano de la Investigaci´ on on Cient´ Cient´ıfica (en Julio–Agosto del 2000), y a los profesores Dr. Eduardo Vivaldo Lima (Facultad de Qu´ Qu´ımica de la UNAM) y Dr. Enrique Sald´ıvar ıvar GueG uerra (Centro de Investigaci´on on y Desarrollo Tecnol´ogiogico de GIRSA), as´ as´ı como a la M. C. Martha Roa Luna, en ese entonces, alumna de postgrado de ambos; por haberme ofrecido un ambiente propicio para aprender de ellos sobre ciencia de pol´ pol´ımeros.
ContactoS ContactoS 68, 54–62 (2008)
Bibl Bi blio iogr graf´ af´ ıa ıa 1. Brandrup, Brandrup, J. y E. H. Immergut Immergut (editors), (editors), Polymer Handbook , John Wiley & Sons, New York, 1989 (3rd edition), pp. 173–182. 2. Flory, Flory, Paul Paul J., “Thermodynamic “Thermodynamicss of High Polymer Polymer Solutions”, J. Chem. Phys., 9[8], pp. 660–661, 1941. 3. Flory, Flory, Paul Paul J., “Thermodynamic “Thermodynamicss of High Polymer Polymer Solutions”, J. Chem. Phys., 10[1], pp. 51–61, 1942. Principles of Polymer Polymer Chemistry Chemistry , 4. Flory Flory,, Paul Paul J., Principles Cornell University Press, Ithaca, 1953, pp. 495–593. 5. Hansen, Hansen, C. M., “50 years years with solubility solubility parameparameProg. Org. Coat., 51[1], pp. ters–past and future”, Prog. 77–84, 2004. 6. Hilde Hildebr bran and, d, J. H. y R. L. Scot Scott, t, “Sol “Solut utio ions ns of Non–electrolytes”, Ann. Rev. Phys. Chem., 1, pp. 75–92, 1950. 7. Huggins, Huggins, Mauric Mauricee L., “Solut “Solution ionss of Long Long Chain Chain Compounds”, J. Chem. Phys., 9[5], pp. 440, 1941. 8. Huggin Huggins, s, Mauric Mauricee L., “Some “Some Proper Propertie tiess of Some Some Phys. Soluti Solutions ons of Long Long Chain Chain Compoun Compounds” ds”,, J. Phys. Chem., 46[1], 151–158, 1942. 9. Lindvig, Lindvig, T., M. L. Michelsen Michelsen y G. M. KontogeorKontogeorgis, “A Flory–Huggins model based on the Hansen solubility parameters”, Fluid Phase Equilib., 203[2], pp. 247–260, 2002. 10. Luna—B´ Luna—B´ arce a rcena nas, s, G., G., Nota Notass de curs curso: o: Termodin´ amica de pol´ pol´ımeros , XIX XIX Semina Seminario rio de IngeIngenier´ıa ıa Qu´ Qu´ımica, Instituto Tecnol´ogic o gico o de Cela Cela-ya, Enero de 2000, p. 1–22. Plastics cs PolyPoly11. 11. Mah Mahen endr dra a D. Baij Baijal al (edi (edito tor) r),, Plasti mer Science and Technology echnology , John Wiley & Sons, 1982, pp. 29–59. 12. Moctez Moctezuma uma,, E. M. y C. Hansmann Hansmann,, Sesi´ Sesi´ on on Oral: on de propi propied edade adess cin´ cin´ etica e ticass de polipoli“Estimaci´ merizaci´ on de estir estireno eno a partir artir de la teor teor´ ´ıa de Flory-–Huggins ”, ”, XV Congre Congreso so Nacion Nacional al de Ter-
modin´ amica, amica, Instituto Tecnol´ ogico ogico de Oaxaca, Septiembre de 2000. 13. Nishimoto, Nishimoto, M., H. Keskkula y D. R. Paul, “Misci“Miscibility of blends of polymers based on styrene, acrylonitrile and methyl methacrylate”, Polymer , 30[7], pp. 1279–1286, 1279–1286, 1989. 14. Paul, Paul, D. R., y J. W. Barlo Barlow, w, “A binary binary intera interacction model for miscibility of copolymers in blends”, Polymer , 25[4], pp. 487–494, 1984. 15. Perry, Perry, R. H., D. Green y J. Maloney (editors), Chemical mical Engineer´s Engineer´s Handbo Handbook ok , McGrawMcGraw-–Hill –Hill,, New York, 1999 (7th edition), pp. 2–120, 2–304, 4–9 y 4–15 a 4–21. 16. Prausnitz, J. M, R. N. Lichtenthale Lichtenthalerr y E. G. de Azeamica molecular de los equilibrios de vedo, Termodin´ fases , Prentice Hall, Espa˜ na, na, 2000 (3ra. edici´ on), on), p. 337–409. cs