Neste trabalho foi estudado o comportamento em altas pressões e temperatura ambiente de dois compostos com estrutura de escuterudita, CoSb3 e LaFe3CoSb12. A determinação dos parâmetros das equações de estado isotérmicas, bem como a verificação das diferenças relativas nas curvas de compressibilidade desses compostos, foram os objetivos principais deste trabalho. Ambos compostos foram sintetizados conforme rota proposta na literatura. O sucesso da síntese foi verificado por meio de difração de raios X, não sendo observada a presença de fases contaminantes. Como meio gerador de altas pressões foi utilizada uma câmara de bigornas de diamante (DAC). A evolução das estruturas destes compostos com a pressão foi acompanhada por meio de difração de raios X por dispersão em energia. Até o limite de pressão estática de cerca de 14 GPa (definido pelas condições experimentais previamente existentes no laboratório), não foram observadas anomalias nas curvas de volume versus pressão desses compostos ou mesmo evidências de transições de fase. Com o objetivo de estender este limite de pressão, foram feitas melhorias na célula de alta pressão existente no laboratório, possibilitando assim a geração de pressões de até 40 GPa. Três equações de estado isotérmicas foram ajustadas aos dados experimentais, utilizando como parâmetros no ajuste o módulo volumétrico isotrópico e sua derivada primeira com relação à pressão. Como resultado deste procedimento, verificou-se que a escuterudita com lantânio apresenta um módulo volumétrico ligeiramente menor que o CoSb3, sendo, portanto, mais compressível, apesar da presença do íon lantânio preenchendo as cavidades deste composto. Foi observada também uma mudança na compressibilidade da escuterudita CoSb3 acima de 20 GPa, devido, talvez, ao congelamento do meio transmissor de pressão ou ao pinçamento do rubi (calibrante de pressão) entre as bigornas da DAC. O CoSb3 foi comprimido até a pressão de 40 GPa, apresentando indícios de um estranho comportamento na seqüência decrescente de pressão. Após a descompressão, o parâmetro de rede desse composto, a pressão ambiente, mostrou-se maior que o da fase original. Este aumento do parâmetro de rede não foi acompanhado de alterações significativas do espectro de difração de raios X. Este efeito anômalo dever ser sujeito a verificação por meio de novos experimentos.

O fator de compressibilidade (Z) de gás natural é utilizado em vários cálculos na engenharia de petróleo (avaliação de formações, perda de carga em tubulações, gradiente de pressão em poços de gás, cálculos de balanço de massa, medição de gás, compressão e processamento de gás). As fontes mais comuns de valores de Z são medições experimentais, caras e demoradas. Essa propriedade também é estimada por correlações empíricas, modelos baseados no princípio dos estados correspondentes ou equações de estado (EOS). Foram avaliadas as capacidades das EOS de Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), Patel-Teja (PT), Patel-Teja-Valderrama (PTV), Schmidt-Wenzel (SW), Lawal-Lake-Silberberg (LLS) e AGA-8 para previsão desta propriedade em aproximadamente 2200 pontos de dados experimentais. Estes pontos foram divididos em quatro grupos: Grupo 1 (Presença de frações C7+, Grupo 2 (temperaturas inferiores a 258,15 K), Grupo 3 (pressões superiores a 10000 kPa) e Grupo 4 (pressões inferiores a 10000 kPa). Os cálculos utilizando as equações de estado sob diferentes esquemas de previsão de coeficientes binários de interação foram cuidadosamente investigados. Os resultados sugerem que a EOS AGA-8 apresenta os menores erros para pressões de até 70000 kPa. Entretanto, observou-se uma tendência de aumento nos desvios médios absolutos em função das concentrações de CO2 e H2S. As EOS PTV e a EOS SW são capazes de predizer o fator de compressibilidade (Z) com desvios médios absolutos entre os valores calculados e experimentais com precisão satisfatória para a maioria das aplicações, para uma variada faixa de temperatura e pressão. Este estudo também apresenta uma avaliação de 224 métodos de cálculo de Z onde foram utilizadas 8 correlações combinadas com 4 regras de mistura para estimativa de temperaturas e pressões pseudorreduzidas das amostras, junto com 7 métodos de caracterização das propriedades críticas da fração C7+, quando presente na composição do gás. Em função dos resultados são sugeridas, para diferentes tipos de sistemas, as melhores combinações de correlações com regras de mistura capazes de predizer fatores de compressibilidade (Z) com os menores erros absolutos médios relativos / The compressibility factor (Z-factor) of natural gases is necessary in many petroleum engineering calculations (evaluation of formation, pressure drop through a pipe, pressure gradient in gas wells, material balance calculations, gas metering, gas compression and processing). The most common sources of Z-factor values are experimental measurements, which are expensive and time consuming. This property is also estimated from empirical correlations, corresponding state models or equations of state (EOS) when experimental data is unavailable. Capabilities of Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), Patel-Teja (PT), Patel-Teja-Valderrama (PTV), Schmidt-Wenzel (SW), Lawal-Lake-Silberberg (LLS) and AGA-8 to predict this property of 2200 data points under various schemes of binary interaction numbers are thoroughly investigated. This database was split on four groups: Group 1 (presence of hydrocarbon plus fraction - C7+), Group 2 (temperatures lower than 258.15 K), Group 3 (pressures higher than 10000 kPa) e Group 4 (pressures lower than 10000 kPa). The best results were obtained from EOS AGA-8 for pressures up to 70000 kPa. However, an increasing trend in average absolute deviations was observed as a function of CO2 e H2S concentrations. EOS PTV and EOS SW are capable to predict the compressibility factor (Z) with average absolute deviation between calculated and experimental values with satisfactory accuracy for most applications for a wide range of temperature and pressure. This study also presents an evaluation of 224 possible methods of calculating the gas compressibility factor, based on 8 correlations and corresponding state models, combined with 4 mixing rule that predict the pseudo-reduced pressure and temperatures of the mixture, that were combined with 7 methods of characterizing the plus fraction critical properties when present in the gas composition. Results suggest for different systems conditions, the best correlation and mixing rule combination capable of predicting Z-factor with the least average absolute relative error

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Atualmente, existem modelos matemáticos capazes de preverem acuradamente as relações entre propriedades de estado; e esta tarefa é extremamente importante no contexto da Engenharia Química, uma vez que estes modelos podem ser empregados para avaliar a performance de processos químicos. Ademais, eles são de fundamental importância para a simulação de reservatórios de petróleo e processos de separação. Estes modelos são conhecidos como equações de estado, e podem ser usados em problemas de equilíbrios de fases, principalmente em equilíbrios líquido-vapor. Recentemente, um teorema matemático foi formulado (Teorema de Redução), fornecendo as condições para a redução de dimensionalidade de problemas de equilíbrios de fases para misturas multicomponentes descritas por equações de estado cúbicas e regras de mistura e combinação clássicas. Este teorema mostra como para uma classe bem definidade de modelos termodinâmicos (equações de estado cúbicas e regras de mistura clássicas), pode-se reduzir a dimensão de vários problemas de equilíbrios de fases. Este método é muito vantajoso para misturas com muitos componentes, promovendo uma redução significativa no tempo de computação e produzindo resultados acurados. Neste trabalho, apresentamos alguns experimentos numéricos com misturas-testes usando a técnica de redução para obter pressões de ponto de orvalho sob especificação de temperaturas. / Nowadays, there are mathematical models capable of predicting accurately the relationships among state properties; and this task is very important in Chemical Engineering context, since these models can be used to evaluate the performance of chemical processes. Moreover, they are of fundamental importance for oil reservoir simulation and separation processes. These models are known as equations of state and can be used in phase equilibrium problems, mainly in vapor-liquid equilibrium. Recently, a mathematical theorem was formulated (Reduction Theorem), which provides the conditions for the reduction in the dimensionality of the phase equilibrium problem for multicomponent mixtures described by cubic equations of state and classical mixing and combination rules. This theorem shows how to a well-defined class of thermodynamic models (cubic equations of state and classical mixing rules), one can reduce the dimensionality of various phase equilibrium problems. This method is very advantageous for mixtures of many components, providing a significant reduction in the computation time and producing accurate results. In this work, we presented some numerical experiments with test-mixtures using the reduction technique to obtain dew point pressures at specified temperatures.

Neste trabalho foi estudado o comportamento em altas pressões e temperatura ambiente de dois compostos com estrutura de escuterudita, CoSb3 e LaFe3CoSb12. A determinação dos parâmetros das equações de estado isotérmicas, bem como a verificação das diferenças relativas nas curvas de compressibilidade desses compostos, foram os objetivos principais deste trabalho. Ambos compostos foram sintetizados conforme rota proposta na literatura. O sucesso da síntese foi verificado por meio de difração de raios X, não sendo observada a presença de fases contaminantes. Como meio gerador de altas pressões foi utilizada uma câmara de bigornas de diamante (DAC). A evolução das estruturas destes compostos com a pressão foi acompanhada por meio de difração de raios X por dispersão em energia. Até o limite de pressão estática de cerca de 14 GPa (definido pelas condições experimentais previamente existentes no laboratório), não foram observadas anomalias nas curvas de volume versus pressão desses compostos ou mesmo evidências de transições de fase. Com o objetivo de estender este limite de pressão, foram feitas melhorias na célula de alta pressão existente no laboratório, possibilitando assim a geração de pressões de até 40 GPa. Três equações de estado isotérmicas foram ajustadas aos dados experimentais, utilizando como parâmetros no ajuste o módulo volumétrico isotrópico e sua derivada primeira com relação à pressão. Como resultado deste procedimento, verificou-se que a escuterudita com lantânio apresenta um módulo volumétrico ligeiramente menor que o CoSb3, sendo, portanto, mais compressível, apesar da presença do íon lantânio preenchendo as cavidades deste composto. Foi observada também uma mudança na compressibilidade da escuterudita CoSb3 acima de 20 GPa, devido, talvez, ao congelamento do meio transmissor de pressão ou ao pinçamento do rubi (calibrante de pressão) entre as bigornas da DAC. O CoSb3 foi comprimido até a pressão de 40 GPa, apresentando indícios de um estranho comportamento na seqüência decrescente de pressão. Após a descompressão, o parâmetro de rede desse composto, a pressão ambiente, mostrou-se maior que o da fase original. Este aumento do parâmetro de rede não foi acompanhado de alterações significativas do espectro de difração de raios X. Este efeito anômalo dever ser sujeito a verificação por meio de novos experimentos.

Orientador: Martin Aznar / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-08T09:32:38Z (GMT). No. of bitstreams: 1 AlvarezAlvarez_VictorHugo_M.pdf: 1355509 bytes, checksum: c9710252adab51dca90a695fc91a5aa6 (MD5) Previous issue date: 2007 / Resumo: Os líquidos iônicos têm recebido considerável interesse devido a seu potencial como solventes projetados, que podem ser adaptados a vários tipos de processos industriais. A razão principal do interesse nos líquidos iônicos como solventes é sua baixa pressão de vapor, a qual minimiza os riscos de exposição e contaminação ambiental. Esta pesquisa visa realizar um estudo completo e sistemático sobre duas regras de mistura utilizando a equação de estado cúbica de Peng-Robinson para correlacionar o equilíbrio líquido-vapor em sistemas contendo líquidos iônicos. O problema principal nesta modelagem é a predição da baixa concentração do líquido iônico na fase vapor.Até o momento são escassas as publicações das propriedades críticas dos líquidos iônicos, motivo pelo qual escolheu-se uma equação com baixa quantidade de parâmetros, como a equação de estado cúbica de Peng-Robinson. Estudou-se o parâmetro dependente da temperatura a(T) e a regra de mistura de van der Waals. Os parâmetros a(T) comparados foram os propostos por Soave (1972) e por Almeida et al. (1991). Nos primeiros testes da modelagem foi utilizada a regra de mistura de van der Waals, mas para melhores resultados utilizou-se a regra de mistura de Wong-Sandler com os modelos UNIQUAC ou NRTL. Desenvolveu-se uma estratégia de modelagem molecular para calcular os parâmetros estruturais de área e volume do modelo UNIQUAC para os líquidos iônicos. Os parâmetros de interação foram calculados independentes da temperatura e concentração. Foram estudadas misturas binárias para descrever o equilíbrio líquido-vapor em altas pressões (CO2 ou CHF3 supercrítico + líquido iônico) e a baixas pressões (hidrocarbonetos + líquidos iônicos). Logo, os parâmetros de interação binária foram utilizados na correlação de sistemas ternários a baixa pressão com a regra de mistura de Wong-Sandler. A estimação dos parâmetros de interação binários foi realizada em um primeiro momento com o método de Levenberg-Marquardt e depois, com melhor sucesso, utilizou-se um algoritmo genético. A função objetivo utilizada contém a pressão do sistema e a composição do líquido iônico na fase gás. Os resultados para a modelagem com a regra de mistura de van der Waals apresentam altos desvios na pressão, mas com a regra de mistura de Wong-Sandler têm-se baixos desvios na pressão e uma baixa concentração do líquido iônico na fase gás até 100 atm. Os resultados mostram uma boa correlação dos sistemas ternários. Foi usado um teste de consistência termodinâmica para o sistema CO2 + hexafluorofosfato de 1-butil-3-metilimidazolio, para o qual há quatro conjuntos de dados conflitantes na literatura. Os resultados mostraram que, um conjunto de dados é termodinamicamente consistente, outro é não inteiramente consistente e os outros dois são termodinamicamente não consistentes / Abstract: The ionic liquids have received considerable interest due to his potential as designer solvents, that can be adapted in several types of industrial processes. The main reason of the interest in the ionic liquid as solvent is his negligible vapor pressure, which decreases the risks of exposition and environmental contamination. This research tries realize a complete and systematic study of two mixing rule in the Peng-Robinson equations of state for correlate and predict the equilibrium liquid-vapor in systems containing ionic liquids. The main problem in modeling the liquid-vapor equilibrium is the prediction of the negligible concentrations of ionic liquid in the phase vapor. Up to the moment the publications of the critical properties of the ionic liquids are scarce, reason for which chose an equation with low amount of parameters, as the cubic equation of state of Peng-Robinson. In this work, was studied the a(T): parameter dependent of the temperature and the mixing rule of van der Waals. The compared a(T): parameters were proposed for Soave 1972) and Almeida et al. (1991). In the first moment of the modeling was used the van der Waals mixing rule, after for good results was used the Wong-Sandler mixing rule with the UNIQUAC or NRTL model. A molecular modeling strategy was used to calculate the volume and surface area parameters of ionic liquids for UNIQUAC. Independent temperature and concentration interaction parameters were calculated. The binary mixtures were studied to describe the liquid-vapor equilibrium included high pressures (CO2 or CHF3 supercritical + ionic liquid) and low pressures (hydrocarbons + ionic liquid). After, the binary interaction parameters were used for correlated the ternary system at low pressure with the Wong-Sandler mixing rule. For evaluating the binary interaction parameters was used in first moment the Levenberg-Marquardt method and after, with best results was used a genetic algorithm. The objective function uses the pressure of the systems and the ionic liquid fraction mol in the gas phase. The results for correlation with van der Waals mixing rule show high deviations in the pressure system, but the Wong-Sandler mixing rule had low deviations in the pressure system and low concentration of the ionic liquid in the gas phase, up to 100 atm. The ternary system can be correlate with acceptable accuracy. A test of thermodynamic consistency was used for a binary system CO2 + 1-n-butyl-3-methylimidazolium hexafluorophosphate, it has four conflicts data set in the literature. The results show one data set thermodynamically consistent, one data set not fully consistent and two data set thermodynamically inconsistent / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Antonio Carlos Bannwart / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica e Instituto de Geociências / Made available in DSpace on 2018-08-27T04:22:37Z (GMT). No. of bitstreams: 1 BarbosaNeto_AntonioMarinho_M.pdf: 2354581 bytes, checksum: dada46fe50ab798d217e21bcad15eeb1 (MD5) Previous issue date: 2015 / Resumo: Fluidos de petróleo, incluindo gás, condensado e óleo leve, são misturas complexas de hidrocarbonetos. O conhecimento sobre o comportamento de fases destas misturas é de interesse para a indústria de petróleo, desde o reservatório, produção, transporte, até o processamento e refino de petróleo. Além disso, pacotes termodinâmicos computacionais para predizer o comportamento de tais misturas são de grande interesse nesta indústria. Neste sentido, uma ferramenta computacional para prever com precisão as propriedades termodinâmicas utilizando cálculos de equilíbrio de fases pode desempenhar um papel importante para aplicações de projeto e operações de campo. Equações de estado cúbicas tem sido amplamente utilizada para predizer o comportamento de fases e propriedades volumétricas de fluidos de petróleo devido a sua simplicidade e alto desempenho computacional. Diante deste contexto, o presente trabalho teve como objetivo principal desenvolver um simulador PVT composicional, de fácil manuseio, porém completo, capaz de avaliar o comportamento volumétrico e de fases de fluidos de petróleo em diferentes condições operacionais (P, T) a partir da modelagem do equilíbrio líquido-vapor e do equilíbrio líquido-líquido-vapor. A ferramenta computacional foi desenvolvida na linguagem de programação Visual Basic Applications (VBA) em aplicativo Excel. Uma interface gráfica foi construída em aplicativo Excel. O programa, denominado PVTpetro, consiste basicamente em três módulos computacionais e suas funções auxiliares. O módulo de cálculo de flash (P, T) 2 fases isotérmico para misturas de hidrocarbonetos consiste no equilíbrio líquido-vapor (ELV) para sistemas óleo-gás multicomponentes. O módulo de cálculo de flash (P, T) 3 fases isotérmico avalia sistemas envolvendo óleo, água e gás (OWG). E, o módulo de cálculo de propriedades volumétricas de fluidos de petróleo que calcula o peso molecular, o volume e a massa específica das fases existentes no sistema. A ferramenta foi validada via comparação entre os resultados de pacotes termodinâmicos comerciais (Wimprop, Hysys, COMThermo) e os do programa desenvolvido. Os resultados gerados pelo PVTpetro foram avaliados através do cálculo do desvio relativo em relação às respostas dos softwares comerciais. Desta forma, verificou-se para o módulo de cálculo flash (P, T) 2 fases que os valores de desvio em relação ao Wimprop e Hysys foram inferiores a 1% para os estudos realizados. Enquanto que, para o módulo de cálculo flash (P, T) 3 fases os valores de erro relativo se concentraram na faixa de 0 ¿ 5% quando comparados com os pacotes termodinâmicos Hysys e COMThermo. De uma maneira geral, os resultados gerados pelo PVTpetro apresentaram-se satisfatórios em todos os cálculos propostos. Na validação do módulo de cálculo das propriedades volumétricas, feita de forma qualitativa, observou-se que o comportamento das propriedades em função da pressão e em diferentes temperaturas, gerado pelo PVTpetro, reproduziram as curvas fornecidas pelo Hysys, em todas as condições investigadas. Portanto, além de se apresentar confiável em seus cálculos o PVTpetro mostrou ser de fácil manuseio e intuitivo para o usuário, facilitando assim o trabalho em vários casos / Abstract: Petroleum fluids, which include gas, condensate and light oil, are complex hydrocarbons mixtures. There is a great deal of interest in the oil industry in computational thermodynamic packages able to predict the phase behavior of such mixtures, from the reservoir, production, transportation until the oil processing and refining. In this sense, a computational tool that accurately predicts thermodynamic properties using phase equilibrium calculations have an important role in designing applications and field operations. Cubic equations of state have been widely used to predict the phase behavior and bulk properties of petroleum fluids due to its simplicity and high computational performance. This study focused on developing a PVT compositional simulator, easy to use, able to assess the bulk properties and phase behavior of petroleum fluids at different operating conditions (P, T) from the modeling of the vapor-liquid equilibrium (VLE) and the vapor-liquid-liquid equilibrium (VLLE). The computational tool was developed in Visual Basic Applications (VBA) programming language, with a graphic interface built in Excel. The program, named PVTpetro, consists in three basic computational modules and their auxiliary functions. The isothermal two phase flash (P, T) module calculation to hydrocarbon mixtures consists of VLE for multi-component oil-gas systems. The isothermal three phase flash (P, T) module calculation evaluates oil, water and gas (OWG) systems. Lastly, the bulk properties of petroleum fluids calculation module calculates the molecular weight, the volume and density of the system phases. The tool has been validated through the comparison between the results of commercial thermodynamic packages (Wimprop, Hysys, COMThermo) with the developed program. The results were evaluated by the relative deviation calculation between these commercial software and PVTpetro results. For the isothermal two-phase flash (P, T) calculation module were observed deviations below 1% for Wimprop and Hysys results. For the isothermal three phase flash (P, T) calculation module the deviations were about 0 and 5% for Hysys and COMThermo. For the bulk properties calculation module the results by PVTpetro and Hysys showed a good agreement in all investigated conditions, ensuring the PVTpetro reliability / Mestrado / Explotação / Mestre em Ciências e Engenharia de Petróleo

Orientador: Martin Aznar / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-05T00:40:21Z (GMT). No. of bitstreams: 1 Arce-Castillo_PedroFelipe_D.pdf: 4607592 bytes, checksum: 216de649e3febd33712f357f832e071f (MD5) Previous issue date: 2005 / Resumo: Durante a presente década, vários resultados experimentais sobre fenômenos críticos no equilíbrio de fases, tais como efeito de co-solvência, janelas de miscibilidade e regiões de duas fases LV em superfícies de três fases LLV em sistemas ternários foram obtidos experimentalmente por vários grupos de pesquisa. Estes fenômenos são de grande influência em aplicações industriais de tecnologia supercrítica ou quase crítica, já que eles são observados dentro do intervalo de pressão, temperatura e concentração onde acontecem estes processos tecnológicos. O presente trabalho de tese visou modelar computacionalmente o equilíbrio multifásico (ELL, ELV, ELLV) e predizer esses fenômenos, que acontecem perto da região crítica, em sistemas simples (CO2 + n-alcano + n-alcanol, CO2 + 2-nitrofenol + n-alcanol) e complexos, tais como polímero + CO2 + cosolvente [polipropileno (iPP) + CO2 + n-pentano; poliestireno (PS) + CO2 + ciclohexano e poli(D,L-lactida) + CO2 + dimetileter)], fazendo uso de equações de estado (EDE) não cúbicas, que possuem algum embasamento teórico, como modelos termodinâmicos. Esses modelos foram a equação da Teoria Estatística do Fluido Associado de Cadeia Perturbada (Perturbed Chain ¿ Statistical Associating Fluid Theory, PC-SAFT) e de Sanchez- Lacombe (SL). Os resultados obtidos na modelagem do equilíbrio de fases a altas pressões desses sistemas foram comparados com os obtidos pela tradicional EDE de Peng-Robinson (PR). Estes modelos termodinâmicos também foram usados na modelagem do ELV e ELL de sistemas binários compostos de copolímeros comuns e biodegradáveis com solventes a baixas e altas pressões. Cada modelo termodinâmico possui parâmetros de componente puro e um parâmetro de interação (na sua regra de cruzamento) para cada sistema binário. Em todos os modelos termodinâmicos foi usada uma regra de cruzamento convencional e uma regra de mistura de primeira ordem de van der Waals. Os parâmetros de componente puro foram obtidos ajustando por regressão os dados de pressão de vapor e do volume molar do líquido saturado para o componente simples e dados de densidade, pressão e temperatura para o polímero na fase líquida. A modelagem de cada sistema binário foi feita a partir de dados de ELL, ELV e ELLV obtendo um parâmetro de interação binária (que leva em conta as interações entre duas moléculas). A predição do comportamento de fases do sistema ternário foi feita usando os parâmetros de interação binária dos três sistemas binários envolvidos e os resultados foram comparados com os dados experimentais. A otimização foi feita usando o método modificado de Máxima Verossimilhança para determinar o ótimo global dos dois tipos de parâmetros. Em todas as modelagens do comportamento do equilíbrio multifásico dos sistemas apresentados acima, a EDE PC-SAFT teve a melhor performance em termos dos desvios relativos na pressão quando comparada à performance das EDEs SL e PR / Abstract: During the present decade, several experimental results on critical phenomena in phase equilibria, such as co-solvency effect, miscibility windows and two-phases LV regions in three phases LLV surfaces in ternary systems were obtained experimentally by several research groups. These phenomena are of great influence in industrial applications of near-critical and supercritical technology, since they are observed within the pressure, temperature and concentration intervals where these technological processes happen. The present thesis aimed for modeling computationally the phase equilibria (LLE, VLE, VLLE) and predict those phenomena that happen near to critical region in simple systems (CO2 + n-alkane + n-alkanol, CO2 + 2-nitrophenol + n-alkanol) and complex systems, such as polymer + CO2 + co-solvent [polypropylene (iPP) + CO2 + n-pentane; polystyrene (PS) + CO2 + cyclohexane and poly(D,L-lactide) + CO2 + dimethyl ether)] using non-cubic equations of state (EoS), with some theoretical base, as thermodynamic models. Those models were the Perturbed Chain - Statistical Associating Fluid Theory (PC-SAFT) and the Sanchez - Lacombe (SL) EoS. The results obtained in modeling of high-pressure phase equilibria of those systems were compared with those obtained by the traditional Peng- Robinson (PR) EoS. These thermodynamic models were also used in modeling the VLE and LLE of binary systems composed of common and biodegradable copolymers with solvents at low and high pressures. Each thermodynamic model has pure component parameters and one interaction parameter (in its combining rule) for each binary system. Conventional combining and van der Waals one-fluid mixing rules were used in all thermodynamic models. Pure component parameters were obtained by regression of liquid saturated vapor pressure and volume molar data for each simple component, and liquid density, pressure and temperature data for each polymer. The modeling of each binary system was made from LLE, VLE and VLLE data, obtaining one binary interaction parameter (which takes into account the interactions between two molecules). The prediction of phase behavior of ternary system was made using the binary interaction parameters of its three binary systems and the results were compared with experimental data. Optimization was made using the modified likelihood maximum method to determine the global optimum of the two types of parameters. In all the modeling of the multiphase equilibria behavior of the systems presented above, the PC-SAFT EoS had the best performance in terms of relative deviations in pressure when compared to performance of SL and PR EoS / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química