Prévia do material em texto
Outline 6. Equilíbrio líquido-vapor (ELV) e a “lei” de Raoult 6.1. Terminologia e representação 6.1.1. Representação gráfica 6.2. A lei de Raoult 6.3. Cálculo de estados de equilíbrio usando a lei de Raoult 6.4. Separação flash 6.4.1. Separação flash multicomponente (c > 2) 6. Equilíbrio líquido-vapor (ELV) e a “lei” de Raoult 6.1. Terminologia e representação Antes de iniciar o estudo específico da termodinâmica do equilíbrio líquido- vapor, especifiquemos a terminologia envolvida em sua descrição e sua representação gráfica. Primeiro, façamos uma análise dos graus de liberdade envolvidos. Um sistema com c componentes e duas fases em equilíbrio (líquida e vapor) tem c graus de liberdade, conforme a regra das fases de Gibbs: 𝐿 = 𝑐 − 𝜋 + 2 = 𝑐 ( 6-1 ) Por facilidade, consideremos inicialmente um sistema formado por dois componentes que coexistam em fases líquida e vapor. Esse sistema tem dois graus de liberdade, ou seja, são necessárias duas variáveis para determinar completamente seu estado. A nomenclatura e as análises são as mesmas para um sistema com mais componentes: apenas a representação gráfica é mais difícil caso haja mais de dois componentes. Normalmente, uma das condições (pressão ou temperatura) e uma das composições (da fase líquida ou da fase vapor) é definida, e o cálculo das condições restantes é realizado. Caso a composição da fase líquida seja conhecida, e se calcule a composição da fase vapor, o cálculo é denominado cálculo de ponto de bolha. Essa denominação se dá por analogia do que aconteceria fisicamente: partindo-se de uma fase líquida, aumentando-se a temperatura (em pressão constante) ou diminuindo-se a pressão (em temperatura constante), uma bolha seria formada, indicando que a região de duas fases foi atingida. Em cálculos de ponto de bolha, deve-se determinar simultaneamente a condição desconhecida (temperatura ou pressão) e a composição da fase vapor. Esta composição, como regra, é diferente da composição da fase líquida. O caso inverso, em que a composição da fase vapor é conhecida, e se calcula a composição da fase líquida, é chamado de cálculo do ponto de orvalho. Novamente, a denominação ocorre por analogia do que ocorreria fisicamente: diminuindo-se a temperatura (em pressão constante) ou aumentando-se a pressão (em temperatura constante), uma gota de líquido seria formada, indicando que a região de duas fases foi atingida. Alguns pontos devem ser esclarecidos. Primeiro, o equilíbrio está completamente determinado ao se conhecerem as variáveis necessárias. Não há relação entre quantidades nesse tipo de cálculo. Uma fase líquida, de composição determinada, a uma certa temperatura, estará em equilíbrio com uma fase vapor de certa composição em certa pressão, independentemente das quantidades dessas fases. Um segundo ponto é que a análise dos graus de liberdade nos diz quantas variáveis independentes devem ser definidas para que o sistema esteja corretamente descrito, mas nada diz sobre quais são essas variáveis ou seus valores. É perfeitamente possível que o problema físico não possa ser resolvido com as variáveis fornecidas: o conjunto de dados pode ser inconsistente, por exemplo, ou efetivamente o sistema pode não existir no estado considerado. 6.1.1. Representação gráfica É muito difícil representar graficamente uma propriedade como função de outras duas (por exemplo, pressão de equilíbrio em função de temperatura e composição do líquido); mais difícil ainda seria representar simultaneamente duas propriedades (pressão e composição do vapor) em função de outras duas. Por isso, dados de equilíbrio líquido-vapor são usualmente apresentados em condições de pressão constante ou de temperatura constante; além de uma razão gráfica, há também uma razão prática: conjuntos de dados experimentais são normalmente obtidos em uma dessas condições. É costume apresentar as curvas de temperatura ou pressão de equilíbrio como função da composição, e não o contrário: dessa maneira, podem-se representar simultaneamente os dados de composição de líquido e de vapor. Vamos considerar a Figura 6-1 seguinte. Ela apresenta dados de ELV para o sistema acetona (1) e acetonitrila (2) em temperatura constante de 323,15 K. Nessa Figura, observamos dados experimentais (pontos) e estados de equilíbrio calculados com auxílio de modelos (curvas). Os dados experimentais são encontrados na literatura, e são resultado de uma meticulosa investigação com experimentos feitos em condições muito controladas. O recurso aos experimentos é sempre necessário – nós sempre necessitamos de dados experimentais para conhecer o comportamento de uma mistura. Por outro lado, os modelos também são importantes. Não é possível obter dados experimentais em todas as situações possíveis: precisamos de modelos que nos forneçam informações sobre o comportamento de modo contínuo. Modelos são fundamentais, por exemplo, no cálculo de processos de separação, em que precisamos determinar o estado de equilíbrio em inúmeras condições não estudadas experimentalmente. Na Figura 6-1, losangos são dados experimentais de pressão de bolha (ou seja, pressão de equilíbrio em função da composição do líquido), enquanto os triângulos são os dados experimentais de pressão de orvalho (ou seja, pressão de equilíbrio em função da composição do vapor). A curva contínua é a curva de pressão de bolha calculada, enquanto a curva pontilhada é a curva de orvalho calculada. Figura 6-1. Diagrama de fases do sistema acetona (1) e acetonitrila (2) em temperatura constante de 323,15 K. Dados experimentais (Elsi, Patel, Abbott & van Ness, J. Chem. Eng. Data 23, 242-245, 1978; van Ness & Kochar, J. Chem. Eng. Data, 12, 38-39, 1967) e modelagem. Por convenção, em sistemas binários o componente mais volátil, ou seja, aquele com maior pressão de saturação em uma mesma temperatura, é identificado como componente (1). Entretanto, como o próprio nome diz, isso é somente uma convenção. Também por convenção, identificamos a composição do líquido pela letra x, e a composição do vapor pela letra y. Vamos analisar as regiões desse gráfico. Acima da curva de bolha temos a região de líquido. Uma mistura nessas condições não apresenta duas fases. Por exemplo, uma mistura em que 𝑥1 = 0,5 e 𝑃 = 0,7 𝑏𝑎𝑟 (ponto a na Figura 6.2) é uma mistura líquida. Não há fase vapor em equilíbrio. Por essa razão, além da temperatura, também podemos estabelecer a pressão e a composição do líquido – como estamos uma região monofásica, há um grau de liberdade a mais. Se abaixarmos a pressão, com composição constante, atingimos a curva de bolha, que é a condição de equilíbrio. No caso, para 𝑥1 = 0,5, isso corresponde à 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.0 0.2 0.4 0.6 0.8 1.0 P / b ar x1, y1 pressão 𝑃 = 0,58 𝑏𝑎𝑟, identificada como ponto a’ na mesma Figura. Da mesma forma, abaixo da curva de orvalho temos a região de vapor. Uma mistura em que 𝑦1 = 0,5 e 𝑃 = 0,35 𝑏𝑎𝑟 (ponto b na Figura 6.2) é uma mistura em fase vapor – não há fase líquida em equilíbrio (e novamente temos um grau de liberdade a mais, por estarmos em região monofásica). Se aumentarmos a pressão, atingimos a curva de orvalho, que também é a condição de equilíbrio. Para que 𝑦1 = 0,5, a pressão correspondente é 𝑃 = 0,48 𝑏𝑎𝑟 (ponto b’). Duas fases em equilíbrio estão necessariamente à mesma pressão – este é um dos critérios de equilíbrio. Desse modo, as composições de duas fases em equilíbrio estão na mesma linha horizontal. Os pontos c e d da Figura 6.2 correspondem a essa situação: o líquido em c (𝑥1 = 0,40) está em equilíbrio com o vapor em d (𝑦1 = 0,62); ambos estão em pressão 𝑃 = 0,53 𝑏𝑎𝑟. A linha que une duas composições em equilíbrio é chamada de linha de amarração: no caso do ELV de sistemas binários, alinha de amarração é a linha horizontal que une as duas composições (a linha pontilhada cd). Um diagrama de equilíbrio nos diz quais são os estados de equilíbrio, mas nada diz sobre a proporção entre as fases em equilíbrio. Assim, os pontos c e d estão em equilíbrio, mas esse fato nada diz sobre suas quantidades: podemos ter apenas uma gota na composição c em equilíbrio com uma fase vapor extensa na composição d, ou apenas uma bolha de composição d em equilíbrio com uma fase líquida extensa na composição c. Na região “interna” às curvas de bolha de orvalho, um sistema não existe como uma única fase estável: ele se separa em uma fase líquida e uma fase vapor. Uma situação assim pode ser encontrada, por exemplo, ao se abaixar subitamente a pressão exercida sobre uma fase líquida. Figura 6-2. Diagrama de fases do sistema acetona (1) e acetonitrila (2) em temperatura constante de 323,15 K. Figura 6-3. Diagrama de fases do sistema acetona (1) e acetonitrila (2) em pressão constante de 1,013 bar. Dados experimentais (Pratt, Trans. Inst. Chem. Engr., 25, 43-67, 1947) e modelagem. 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.0 0.2 0.4 0.6 0.8 1.0 P / b ar x1, y1 a b c d a' b' 325.0 330.0 335.0 340.0 345.0 350.0 355.0 360.0 0.0 0.2 0.4 0.6 0.8 1.0 T / K x1, y1 O mesmo sistema pode ser analisado como auxílio de um diagrama em pressão constante. Nesse caso, a Figura 6-3 anterior seria obtida. A curva em pressão constante é, em certo sentido, invertida: a região de líquido está em baixas temperaturas, e a região de vapor, em altas temperaturas. O composto mais volátil é o que tem menor temperatura de saturação. A convenção de símbolos e tracejados é a mesma da curva em temperatura constante. 6.2. A lei de Raoult Em sistemas em baixas pressões duas simplificações são bastante comuns: A fase vapor pode ser considerada ideal. Excetuando-se alguns poucos componentes que se associam fortemente em fase vapor (como o ácido acético e o ácido fluorídrico), em geral pode-se considerar que a fase vapor é uma mistura de gases ideais para pressões da ordem da pressão atmosférica ou superiores. As propriedades da fase líquida dependem muito pouco da pressão. Também nessa situação pode-se considerar que a fase líquida é praticamente independente da pressão. Uma dependência considerável ocorre somente em sistemas a altas pressões ou próximos de condições críticas. Esses dois fatos permitem que se estabeleçam algumas relações úteis, válidas para misturas ideais em fase líquida. Por exemplo, sabe-se que em uma mistura ideal, valerá: 𝑓𝑖 𝐿[𝑇, 𝑃, 𝒙] = 𝑥𝑖𝑓𝑖 𝐿[𝑇, 𝑃] ( 6-2 ) ou seja, existe uma proporcionalidade entre a fugacidade do componente e a fração desse mesmo componente na mistura. Em uma situação de equilíbrio vale também a igualdade da fugacidade: 𝑓𝑖 𝐿[𝑇, 𝑃, 𝒙] = 𝑓𝑖 𝑉[𝑇, 𝑃, 𝒚] ( 6-3 ) Como a fase vapor é constituída por uma mistura de gases ideais: 𝑓𝑖 𝑉[𝑇, 𝑃, 𝒚] = 𝑦𝑖𝑃 ( 6-4 ) Resta descobrir o valor da fugacidade do composto puro em fase líquida, para estabelecer a relação entre as duas fases. Para isso, considera-se a situação de equilíbrio, em que a fugacidade da fase vapor será igual à pressão de saturação, ou seja: 𝑓𝑖 𝐿[𝑇, 𝑃𝑖 𝑠𝑎𝑡] = 𝑃𝑖 𝑠𝑎𝑡[𝑇] ( 6-5 ) Considerando que a fugacidade da fase líquida não se altere com a pressão, tem-se: 𝑓𝑖 𝐿[𝑇, 𝑃, 𝒙] = 𝑥𝑖𝑃𝑖 𝑠𝑎𝑡[𝑇] ( 6-6 ) Igualando a fugacidade da fase líquida à da fase vapor, tem-se para a situação de equilíbrio: 𝑦𝑖𝑃 = 𝑥𝑖𝑃𝑖 𝑠𝑎𝑡[𝑇] ( 6-7 ) Essa é a chamada lei de Raoult. Deve-se notar que ela não é uma lei, como as leis da Termodinâmica – é uma simplificação válida em alguns poucos casos especiais, baixas pressões e fase líquida ideal. Apesar de todas as considerações feitas, ela é muito útil no estudo do equilíbrio líquido-vapor de sistemas formados por compostos semelhantes. As curvas calculadas para o sistema acetona (1) e acetonitrila (2) anteriormente apresentadas foram calculadas com auxílio da lei de Raoult. Embora a lei de Raoult seja usada principalmente para sistemas de compostos pouco polares, como misturas de hidrocarbonentos, ela fornece bons resultados até mesmo para misturas de compostos auto-associativos (capazes de formar pontes de hidrogênio), desde que esses compostos sejam semelhantes entre si e as pressões sejam baixas, como se pode ver na Figura seguinte. Figura 6-4. Diagrama de fases do sistema metanol (1) e etanol (2) em pressão constante de 1,013 bar. Dados experimentais (Amer et al., Ind. Eng. Chem., 48, 142- 146, 1956) e modelagem usando a lei de Raoult. 336.0 338.0 340.0 342.0 344.0 346.0 348.0 350.0 352.0 354.0 0 0.2 0.4 0.6 0.8 1 T / K x1 , y1 6.3. Cálculo de estados de equilíbrio usando a lei de Raoult De maneira geral, cálculos de equilíbrio líquido vapor (como outros tipos de equilíbrio) são iterativos: parte-se de uma estimativa inicial e corrige-se a estimativa de modo a buscar uma condição em que os critérios de equilíbrio sejam satisfeitos. Para pressões baixas, em que a fase vapor pode ser considerada ideal e a fase líquida independente da pressão, o cálculo é bastante simplificado. Como regra, cálculos de condição de bolha são mais simples que de orvalho, e cálculos de pressão são mais simples que de temperatura – em situações extremas, são necessários algoritmos robustos específicos. A melhor maneira de aprendermos a fazer o cálculo de estados de equilíbrio… é fazendo cálculos de estados de equilíbrio. Os exercícios resolvidos a seguir mostram os tipos mais comuns de cálculo de ELV usando a lei de Raoult. Exemplo 6-1. Seja o sistema formado por benzeno (1) e tolueno (2). A pressão de vapor dos compostos puros é dada pela seguinte equação: 𝑙𝑛 𝑃𝑖 𝑣𝑎𝑝 𝑏𝑎𝑟 = 𝐴𝑖 − 𝐵𝑖 (𝑇/𝐾) + 𝐶𝑖 ( 6-8 ) válida para 𝑇𝑖 𝑚𝑖𝑛 ≤ 𝑇 ≤ 𝑇𝑖 𝑚𝑎𝑥, conforme a Tabela: Composto 𝐴𝑖 𝐵𝑖 𝐶𝑖 𝑇𝑖 𝑚𝑖𝑛/𝐾 𝑇𝑖 𝑚𝑎𝑥/𝐾 1a 10,8816 3823,79 -1,461 333,4 373,5 2b 9,3905 3094,54 -53,773 308,5 384,7 a – Eon, C., Pommier, C. Guiochon, G., J. Chem. Eng. Data, 16, 408-410, 1971. b – Willingham, C. B., Taylor, W. J., Pignocco, J. M., Rossini, F.D., J. Res. Natl. Bur. Stand., 35, 219-244, 1945. Calcular a pressão de bolha e a composição da fase vapor para T = 360 K. Resolução. A primeira questão que precisamos resolver é: vale a lei de Raoult? Observamos que nessa temperatura, 𝑃1 𝑣𝑎𝑝 = 1,242 𝑏𝑎𝑟 e 𝑃2 𝑣𝑎𝑝 = 0,489 𝑏𝑎𝑟. Portanto, as pressões de equilíbrio devem ser baixas, o que justifica considerar que a fase vapor seja uma mistura de gases ideais, e que a fase líquida independa da pressão. Ademais, benzeno e tolueno são suficientemente semelhantes para que possamos considerar que a fase líquida seja uma mistura ideal. Nesse caso, temos um problema em que são conhecidas a temperatura e a composição da fase líquida, e devem ser determinadas a pressão e a composição da fase vapor. Trata-se de um cálculo de pressão de bolha. O sistema de equações que precisamos resolver é: 𝑦1𝑃 = 𝑥1𝑃1 𝑣𝑎𝑝[𝑇] ( 6-9 ) 𝑦2𝑃 = 𝑥2𝑃2 𝑣𝑎𝑝[𝑇] ( 6-10 ) As incógnitas são 𝑦1, 𝑦2e P. De fato, 𝑦2 não é uma incógnita propriamente dita, pois está ligada a 𝑦1 por meio de: 𝑦1 + 𝑦2 = 1,0 ( 6-11 ) A maneira mais simples de resolver essas equações é acoplá-las, somando- as: 𝑃 = 𝑥1𝑃1 𝑣𝑎𝑝 + 𝑥2𝑃2 𝑣𝑎𝑝 ( 6-12 ) O valor de 𝑦1 é calculado a partir da equação correspondente: 𝑦1 = 𝑥1𝑃1 𝑣𝑎𝑝 𝑥1𝑃1 𝑣𝑎𝑝 + 𝑥2𝑃2 𝑣𝑎𝑝 ( 6-13 ) Com isso, pode-se construir a tabela de resultados: 𝑥1 𝑃 / 𝑏𝑎𝑟 𝑦1 0,0 0,489 0,0 0,1 0,565 0,220 0,2 0,640 0,388 0,3 0,715 0,521 0,4 0,790 0,629 0,5 0,866 0,717 0,6 0,941 0,792 0,7 1,016 0,856 0,8 1,091 0,910 0,9 1,167 0,958 1,0 1,242 1,0 Algumas coisas são interessantes de seremnotadas. Primeiro, a pressão se comporta linearmente com respeito a 𝑥1. Isso é consequência direta da forma da equação da pressão. Em segundo lugar, a fase vapor é proporcionalmente mais rica no componente mais volátil (ou seja, temos que 𝑦1 ≥ 𝑥1 e 𝑦2 ≤ 𝑥2, sendo que as igualdades só ocorrem no limite dos compostos puros). Isso não é verdadeiro sempre em problemas de equilíbrio líquido vapor, mas é verdadeiro sempre que houver um problema que use a lei de Raoult. Exemplo 6-2. Seja o mesmo sistema formado por benzeno (1) e tolueno (2), com as informações dadas no exercício anterior. Calcular a pressão de orvalho e a composição da fase líquida para T = 350,0 K. Resolução. Nesse caso, temos um problema em que são conhecidas a temperatura e a composição da fase vapor, e devem ser determinadas a pressão e a composição da fase líquida. Trata-se de um cálculo de pressão de orvalho. O sistema de equações que precisamos resolver é o mesmo. A diferença é que, nesse caso, as incógnitas são 𝑥1, 𝑥2e P, considerando também a seguinte equação: 𝑥1 + 𝑥2 = 1,0 ( 6-14 ) Nesse caso, isolando os valores da composição da fase líquida nas expressões da lei de Raoult, e substituindo nessa equação, temos: 𝑦1𝑃 𝑃1 𝑣𝑎𝑝 + 𝑦2𝑃 𝑃2 𝑣𝑎𝑝 = 1,0 ( 6-15 ) que resulta em: 𝑃 = ( 𝑦1 𝑃1 𝑣𝑎𝑝 + 𝑦2 𝑃2 𝑣𝑎𝑝) −1 ( 6-16 ) Observamos que nessa temperatura, 𝑃1 𝑣𝑎𝑝 = 0,915 𝑏𝑎𝑟 e 𝑃2 𝑣𝑎𝑝 = 0,348 𝑏𝑎𝑟. Com isso, pode-se construir a tabela de resultados seguinte. Deve-se notar que a relação da composição da fase vapor e da fase líquida se mantém (𝑦1 ≥ 𝑥1 e 𝑦2 ≤ 𝑥2). Obviamente, P não é uma função linear de 𝑦1. 𝑦1 𝑃 / 𝑏𝑎𝑟 𝑥1 0,0 0,348 0,0 0,1 0,371 0,041 0,2 0,397 0,086 0,3 0,427 0,140 0,4 0,462 0,202 0,5 0,504 0,276 0,6 0,554 0,363 0,7 0,614 0,470 0,8 0,690 0,603 0,9 0,786 0,774 1,0 0,915 1,0 Nota: um conceito que aparece bastante no cálculo de colunas de destilação é a volatilidade relativa. Ela é definida exatamente por: 𝛼12 = 𝑦1 𝑥1 ⁄ 𝑦2 𝑥2 ⁄ ( 6-17 ) Sendo válida a lei de Raoult, a volatilidade relativa torna-se: 𝛼12 𝑅𝑎𝑜𝑢𝑙𝑡 = 𝑃1 𝑣𝑎𝑝 𝑃2 𝑣𝑎𝑝 ( 6-18 ) Nota-se que 𝛼12 𝑅𝑎𝑜𝑢𝑙𝑡[360 𝐾] = 2,54 e 𝛼12 𝑅𝑎𝑜𝑢𝑙𝑡[350 𝐾] = 2,63. Exemplo 6-3. Seja o mesmo sistema formado por benzeno (1) e tolueno (2), com as informações dadas no exercício anterior. Qual a temperatura em que a pressão de equilíbrio de uma mistura líquida com 𝑥1 = 0,5 é igual a 1,0 bar. Resolução. Nesse caso, temos um problema em que são conhecidas a pressão e a composição da fase líquida, e devem ser determinadas a temperatura e a composição da fase vapor. Trata-se de um cálculo de temperatura de bolha. As equações disponíveis são as mesmas que para o cálculo da pressão de bolha. Entretanto, o cálculo, nesse caso, é iterativo, pois considerando que: 𝑃 = 𝑥1𝑃1 𝑣𝑎𝑝 + 𝑥2𝑃2 𝑣𝑎𝑝 ( 6-19 ) e as equações para pressão de vapor, não é possível isolar a temperatura na equação acima. Isso acontece mesmo considerando a lei de Raoult, e deve-se à forma da equação para a pressão de vapor. A resposta ao problema é 𝑇 = 364,77 𝐾. O cálculo é ligeiramente mais difícil do que os cálculos de pressão de equilíbrio, mas ainda é factível, e pode ser facilmente resolvido com o auxílio de uma planilha de cálculo. 6.4. Separação flash Consideremos o seguinte problema, uma continuação dos exercícios anteriores: Exemplo 6-4. Uma mistura de benzeno (1) e tolueno (2) cuja composição é 𝑧1 = 0,5, atinge uma condição em que 𝑇 = 370 𝐾 e 𝑃 = 1,0 𝑏𝑎𝑟. Pergunta-se: qual a composição das fases em equilíbrio nessa condição de temperatura e pressão? Qual a fração da mistura original que estará em fase vapor no equilíbrio? Esse tipo de cálculo é chamado de cálculo flash. Normalmente, ele ocorre quando uma corrente inicialmente líquida atinge uma condição em que a pressão é menor que sua pressão de bolha. Um esquema geral pode ser visto na Figura seguinte. Figura 6-5. Esquema geral de um vaso de separação flash. Para um sistema binário em que vale a lei de Raoult, como é o caso, a resolução é bastante simples, e pode ser feita analiticamente. Notamos, inicialmente, que na condição de temperatura do sistema, tem-se 𝑃1 𝑣𝑎𝑝 = 1,659 𝑏𝑎𝑟 e 𝑃2 𝑣𝑎𝑝 = 0,6735 𝑏𝑎𝑟. Note-se que nada é dito sobre o estado da mistura inicial – nem é necessário, pois o estado de equilíbrio não depende de como ele é atingido. A pressão de bolha calculada para a composição de 𝑧1 = 0,5 é igual a 1,166 bar, e a pressão de orvalho, 0,958 bar. Ou seja, a pressão é inferior à pressão de bolha, mas superior à pressão de orvalho, o que garante que estejamos em uma situação de equilíbrio líquido-vapor. Nesse caso, usando a relação: 𝑃 = 𝑥1𝑃1 𝑣𝑎𝑝 + 𝑥2𝑃2 𝑣𝑎𝑝 = 𝑥1𝑃1 𝑣𝑎𝑝 + (1 − 𝑥1)𝑃2 𝑣𝑎𝑝 ( 6-20 ) pode-se isolar 𝑥1 e obter: 𝑥1 = 𝑃 − 𝑃2 𝑣𝑎𝑝 𝑃1 𝑣𝑎𝑝 − 𝑃2 𝑣𝑎𝑝 ( 6-21 ) Pelo balanço de massa, considerando-se uma alimentação (F), uma fase vapor (V) e uma fase líquida (L), tem-se: 𝐹 = 𝑉 + 𝐿 ( 6-22 ) 𝑧1𝐹 = 𝑦1𝑉 + 𝑥1𝐿 ( 6-23 ) Tomando-se uma base de cálculo 𝐹 = 1, chamando-se 𝑉 = 𝜓, tem-se: 𝑧1 = 𝑦1𝜓 + 𝑥1(1 − 𝜓) ( 6-24 ) que resulta em: 𝜓 = 𝑧1 − 𝑥1 𝑦1 − 𝑥1 ( 6-25 ) No caso específico, temos que 𝑥1 = 0,331 e 𝑦1 = 0,549. Essas são as composições das fases líquida e vapor em equilíbrio a 𝑇 = 370 𝐾 e 𝑃 = 1,0 𝑏𝑎𝑟. Com esses resultados, obtém-se 𝜓 = 0,772, ou seja, 77,2% da mistura original estará na fase vapor quando ela atingir o equilíbrio nas condições especificadas (considerando uma base molar). Uma pergunta: o que ocorreria anterior se 𝑧1 = 0,3? E se 𝑧1 = 0,7? Pelas equações anteriores, a composição das fases em equilíbrio (𝑥1 = 0,331 e 𝑦1 = 0,549) não se altera. Fazendo-se os cálculos, observa-se que para 𝑧1 = 0,3, tem-se 𝜓 = −0,144. Esse resultado aparentemente absurdo significa que somente há fase líquida no sistema em equilíbrio, e não há fase vapor; tal fase líquida é muitas vezes conhecida por líquido sub-resfriado, por estar em temperatura abaixo da temperatura de bolha. Por outro lado, para 𝑧1 = 0,7, tem-se 𝜓 = 1,690. Nesse caso, inversamente, só há fase vapor no sistema, não havendo fase líquida no equilíbrio; tal fase vapor é conhecida por vapor superaquecido, por razões análogas. 6.4.1. Separação flash multicomponente (c > 2) A análise de uma separação flash fica mais complicada quando a mistura tiver mais de dois componentes. Nesse caso, necessariamente o cálculo do equilíbrio e o balanço de massa devem ser acoplados: não é possível deduzir, nesse caso, qual seria a composição das fases em equilíbrio conhecendo-se somente a pressão e a temperatura, em virtude do número de graus de liberdade envolvidos. Nesse caso específico, define-se a variável 𝐾𝑖: 𝐾𝑖 = 𝑦𝑖 𝑥𝑖 ( 6-26 ) que, para um sistema que siga a lei de Raoult, torna-se: 𝐾𝑖 = 𝑃𝑖 𝑠𝑎𝑡(𝑇) 𝑃 ( 6-27 ) O balanço de massa para o componente i pode ser escrito como: 𝑧𝑖 = (1 − 𝜓)𝑥𝑖 + 𝜓𝑦𝑖 ( 6-28 ) Substituindo-se o valor de 𝑦𝑖: 𝑧𝑖 = (1 − 𝜓)𝑥𝑖 + 𝜓𝐾𝑖𝑥𝑖 ( 6-29 ) e isolando-se 𝑥𝑖: 𝑥𝑖 = 𝑧𝑖 1 − 𝜓(1 − 𝐾𝑖) ( 6-30 ) O valor de 𝜓 poderia em princípio ser calculado sabendo-se que a soma de todos os valores de x deve ser igual à unidade, ou seja: ∑ 𝑥𝑗 𝑗 = 1 ⇒ ∑ 𝑧𝑗 1 − 𝜓(1 − 𝐾𝑗) 𝑗 = 1 ( 6-31 ) Entretanto, essa equação tem uma raiz trivial, a saber 𝜓 = 0. Poderíamos escrever uma equação análoga para as frações molares em fase vapor: ∑ 𝑦𝑗 𝑗 = 1 ⇒ ∑ 𝐾𝑗𝑥𝑗 𝑗 = 1 ⇒ ∑ 𝐾𝑗𝑧𝑗 1 − 𝜓(1 − 𝐾𝑗) 𝑗 = 1 ( 6-32 ) mas nesse caso teríamos a raiz trivial 𝜓 = 1. A solução para esse problema está em considerar que: ∑(𝑥𝑗 − 𝑦𝑗) 𝑗 = 0 ⇒ ∑ (1 − 𝐾𝑗)𝑧𝑗 1 − 𝜓(1 − 𝐾𝑗) 𝑗 = 0 ( 6-33 ) Essa equaçãoé conhecida por equação de Rachford-Rice. Para um sistema com c componentes, ela é na verdade uma equação polinomial de grau 𝑐 − 1. Com isso, já podemos prever que haverá raízes que serão simplesmente raízes matemáticas, mas não físicas, do problema. Vejamos um caso no exemplo seguinte. Exemplo 6-5. Seja o sistema formado por n-hexano (1), n-heptano (2) e n-octano (3). A pressão de vapor dos compostos puros é dada pela seguinte equação: 𝑙𝑛 𝑃𝑖 𝑣𝑎𝑝 𝑏𝑎𝑟 = 𝐴𝑖 − 𝐵𝑖 (𝑇/𝐾) + 𝐶𝑖 ( 6-34 ) válida para 𝑇𝑖 𝑚𝑖𝑛 ≤ 𝑇 ≤ 𝑇𝑖 𝑚𝑎𝑥, conforme a Tabela: Composto 𝐴𝑖 𝐵𝑖 𝐶𝑖 𝑇𝑖 𝑚𝑖𝑛/𝐾 𝑇𝑖 𝑚𝑎𝑥/𝐾 1a 9,52214 2869,781 -40,162 178,2 510,2 2b 9,27555 2921,142 -56,199 299,1 372,4 3b 9,32241 3120,293 -63,633 326,7 399,7 a – dbonline.ddbst.de/AntoineCalculation/AntoineCalculationCGI.exe?component=Hexane b – Willingham, C. B., Taylor, W. J., Pignocco, J. M., Rossini, F.D., J. Res. Natl. Bur. Stand., 35, 219-244, 1945. Uma corrente contendo esses três compostos (𝑧1 = 0,35 e 𝑧2 = 0,35) atinge uma condição em que 𝑇 = 370 𝐾 e 𝑃 = 1,0 𝑏𝑎𝑟. Qual o estado da mistura? Caso a mistura esteja em equilíbrio líquido-vapor, qual a composição das fases em equilíbrio? Resolução. O cálculo das pressões de equilíbrio, nesta temperatura, fornece 𝑃1 𝑣𝑎𝑝 = 2,274 𝑏𝑎𝑟, 𝑃2 𝑣𝑎𝑝 = 0,967 𝑏𝑎𝑟 e 𝑃3 𝑣𝑎𝑝 = 0,422 𝑏𝑎𝑟. Com a composição fornecida, calcularíamos uma pressão de bolha de 1,261 bar e uma pressão de orvalho de 0,815 bar. Ou seja, a mistura em questão está em equilíbrio líquido- vapor, pois a pressão se encontra entre esses dois limites. Pela equação de Rachford-Rice: ∑ (1 − 𝐾𝑗)𝑧𝑗 1 − 𝜓(1 − 𝐾𝑗) 3 𝑗=1 = 0 ( 6-35 ) No caso específico, o fato de a pressão do sistema ser igual a 1,0 bar faz com que os valores de K sejam numericamente iguais às pressões de vapor. Calculando- se, obtemos 𝜓 = 0,5413, com as composições de fase líquida e vapor iguais a: 𝑥1 = 0,2071 𝑦1 = 0,4711 𝑥2 = 0,3563 𝑦2 = 0,3446 𝑥3 = 0,4366 𝑦3 = 0,1843 Note-se que a equação de Rachford-Rice, nesse caso, é uma equação do segundo grau. A segunda raiz (𝜓 = 19,98) é espúria, pois fornece 𝑥1 = 0,0132, 𝑥2 = 1,0156 e 𝑥3 = −0,0285, o que é um completo absurdo. Note-se que não é o fato de 𝜓 ser negativo ou maior que 1,0 que torna o resultado necessariamente espúrio, como visto anteriormente para um sistema binário: tal conclusão é obtida da análise das composições das fases em equilíbrio.