Neste trabalho estudou-se a modelagem termodinâmica do equilíbrio de fases de sistemas aquosos bifásicos formados por um polímero neutro, especificamente o poli(etileno glicol), e um sal. O modelo implementado é baseado na equação de Pitzer para soluções de eletrólitos, dele diferindo ao considerar a influência das propriedades do polímero no cálculo da constante dielétrica e da densidade do meio, presentes no cálculo do termo relativo às interações de longo alcance da energia de Gibbs excedente. Além disso, procurou-se também considerar a dependência dos parâmetros relativos ao polímero com respeito ao tamanho de sua cadeia. Dados experimentais obtidos da literatura foram utilizados para a obtenção dos parâmetros ajustáveis do modelo, em cuja implementação foi utilizada a linguagem FORTRAN. A alteração nas propriedades do meio modificou o estado de equilíbrio calculado, fato que pôde ser observado por meio da análise de fatores como o desvio entre a composição calculada pelos modelos e os valores experimentais em cada fase e a inclinação da linha de amarração. Em alguns sistemas o desempenho do modelo modificado foi superior nesses parâmetros, mas em termos gerais a modificação no cálculo da constante dielétrica e da densidade não resultou em melhora no cálculo do equilíbrio de fases, por vezes piorando-o, de modo que investigações específicas sobre o termo de longo alcance foram conduzidas para verificar a razão deste comportamento. A inclusão de parâmetros ternários aumentou a capacidade de correlação de ambos os modelos. No que diz respeito à possibilidade de correlação dos parâmetros de interação como função do tamanho da cadeia polimérica, que poderia levar a um modelo preditivo, observou-se que somente para tamanhos grandes de cadeia foi possível estabelecer uma dependência direta e inequívoca.

The thermodynamic modeling of phase equilibrium in aqueous two phase systems containing the neutral polymer poly(ethylene glycol) and a salt was studied in this work. The implemented model is based on the Pitzer equation for electrolytes solutions, modified in order to account for the influence of polymer properties in both the medium dielectric constant and the medium density, values which are present in long range term of the excess Gibbs energy. The dependence of the adjustable parameters on the polymer molecule size was also investigated. Experimental data from literature were used to obtain the adjustable parameters of the model, in whose implementation the FORTRAN computer language was used. The changes introduced into the long range term resulted in a shift of the calculated equilibrium compositions, which could be observed by analyzing the deviation between calculated and experimental compositions and tie-line slopes. For some systems the performance of the modified model was superior, but in most cases the changes introduced did not result in a significant improvement in phase equilibrium calculations, and even worsened them. Specific investigations on the long range term were carried out to verify the reason for such behavior. The insertion of ternary parameters increased the correlation capacity of both models. Concerning the hypothesis that the interaction parameters can be directly related to the polymer chain size, which might eventually lead to a predictive model, it was noticed that a direct dependency holds only for larger chain sizes.