O presente trabalho consiste em testar a consistência termodinâmica dos dados experimentais obtidos da literatura para o equilíbrio isotérmico (T=373,15 K) do sistema Metanol – CH3OH (1) e Polietileno glicol dimetil éter (PEGDME 250) - CH3O(C2H4O)nCH3 (2), os quais foram submetidos aos testes de área e dos desvios. Além disso, fez-se uso da modelagem e simulação, com aplicação de modelos termodinâmicos, para predição dos dados experimentais com o auxílio de simuladores e ferramentas matemáticas como o CONSIST, ANTOINE e SPECS. Para calcular as propriedades críticas do polímero PEGDME 250, utilizou-se o método de contribuição de grupos de Constantinou e Gani (1994). Em referência aos modelos, utilizou-se a abordagem gama-phi com o modelo de coeficiente de atividade UNIFAC aliado a equação de estado de Soave-Redlich-Kwong (SRK) e a regra de mistura QMR, além da abordagem phi-phi, ao empregar a equação de estado de Soave-Redlich-Kwong com associação CPA (Cubic-Plus-Association) aliado a regra de mistura QMR. Por fim, fez-se uso do modelo de cálculo do coeficiente de atividade com a equação de Hildebrand & Scatchard. Em relação ao teste da área, os dados isotérmicos do sistema trabalhado não foram consistentes, uma vez que o logaritmo da razão entre os coeficientes de atividade resultou em um valor acima de 0,01. O modelo que melhor ajustou-se aos dados experimentais foi através do UNIFAC aliado a equação de estado de Soave-Redlich-Kwong (SRK). Por fim, observou-se um erro médio com a equação de Hildebrand & Scatchard para o cálculo da fase vapor de 0,58% e 31,54% no cálculo da pressão.