Resumos
No processo de extração do óleo de soja forma-se a Miscela, que é uma mistura de óleo e hexano. O processo de separação destes dois componentes é denominado de Destilação da Miscela. Esta é composta basicamente da operação de evaporação e da operação de stripping do hexano. Neste trabalho, desenvolveu-se uma modelagem matemática para o stripping do hexano. O modelo é baseado em balanços de massa e energia e em relações de equilíbrio; e, a sua validação foi feita comparando-se os resultados obtidos com os dados operacionais da indústria da COAMO (Cooperativa Agrícola Mouraoense). Em seguida, com o modelo elaborado, fez-se as simulações para estudar: as influências da concentração da alimentação e do modo de injeção de vapor direto nos resultados operacionais do stripper.
modelagem; simulação; miscela; stripping
In the soybean oil extraction process make the Miscella that is a mixture of oil and hexane. The separation process in both components is named Miscella Distillation. It is composed mainly of the hexane evaporation and of the hexane stripping. In this work, it was developed a mathematical modeling of the hexane stripping. This model is based in mass and energy balances and in equilibrium relations. The validity of the model was made comparing the obtained results with data from a industry of the COAMO (Cooperativa Agrícola Mouraoense). In following, with the elaborated model, it was made the simulation of the stripping to study: the influences of the feed flow and the injection form of direct steam in the operation results of the stripper.
modeling; simulation; miscella; stripping
Destilação da miscela II: modelagem e simulação do stripping do hexano
Distillation of the miscella II: modeling and simulation of the hexane stripping
Paulo R. ParaísoI, * * A quem a correspondência deve ser enviada. ; Cid M. G. AndradeI; Roger J. ZempII
IDepartamento de Engenharia Química, UEM- Maringá. E-mail: paulo@deq.uem.br, cid@deq.uem.br
IIDepartamento de Sistemas Químicos-FEQ/UNICAMP-Campinas. E-mail: zemp@desq.feq.unicamp.br
RESUMO
No processo de extração do óleo de soja forma-se a Miscela, que é uma mistura de óleo e hexano. O processo de separação destes dois componentes é denominado de Destilação da Miscela. Esta é composta basicamente da operação de evaporação e da operação de stripping do hexano. Neste trabalho, desenvolveu-se uma modelagem matemática para o stripping do hexano. O modelo é baseado em balanços de massa e energia e em relações de equilíbrio; e, a sua validação foi feita comparando-se os resultados obtidos com os dados operacionais da indústria da COAMO (Cooperativa Agrícola Mouraoense). Em seguida, com o modelo elaborado, fez-se as simulações para estudar: as influências da concentração da alimentação e do modo de injeção de vapor direto nos resultados operacionais do stripper.
Palavras-chave: modelagem; simulação; miscela; stripping.
SUMMARY
In the soybean oil extraction process make the Miscella that is a mixture of oil and hexane. The separation process in both components is named Miscella Distillation. It is composed mainly of the hexane evaporation and of the hexane stripping. In this work, it was developed a mathematical modeling of the hexane stripping. This model is based in mass and energy balances and in equilibrium relations. The validity of the model was made comparing the obtained results with data from a industry of the COAMO (Cooperativa Agrícola Mouraoense). In following, with the elaborated model, it was made the simulation of the stripping to study: the influences of the feed flow and the injection form of direct steam in the operation results of the stripper.
Keywords: modeling; simulation; miscella; stripping.
1 _ INTRODUÇÃO
A industrialização de oleaginosas constitui-se num importante setor do sistema agroindustrial, pela importância de seus produtos nas indústrias siderúrgicas, de cosméticos e como matéria-prima na indústia de alimentos para o consumo animal e humano, segundo BARBOSA FREITAS & FRANÇA [1]. Os produtos gerados atendem o mercado interno e externo, e para serem competitivos, necessitam sempre buscar a redução de custos.
Uma estratégia visando a redução de custos, bem como a proteção do meio ambiente e do consumidor, nas indústrias já instaladas, é procurar sempre estabelecer as melhores condições de operação do processo. Tais condições podem propiciar um consumo menor de energia e um nível mais elevado de recuperação de hexano que podem ser importantes tanto na redução de custos quanto na proteção do meio ambiente e do consumidor. As indústrias instaladas no Paraná/Brasil utilizam o processo contínuo de extração por solvente, sendo o solvente utilizado o hexano. As operações do processo que deverão ser as principais numa estratégia de redução de custos operacionais são: a extração do óleo por solvente, a destilação da miscela e a dessolventização-tostagem do farelo uma vez que estas consomem muita energia e a manipulam com o hexano.
A miscela é uma mistura entre óleo de soja e o hexano formada no extrator. A separação destes componentes é denominada de Destilação da Miscela, que se compõe: da evaporação e do stripping do hexano. Nestas operações, quantidades expressivas de energia são utilizadas, na forma de vapor.
Este trabalho apresenta a modelagem matemática e simulação da etapa de evaporação do hexano. A validação do modelo desenvolvido sera feito através da comparação dos dados gerados com os dados operacionais da indústria da COAMO (Cooperativa Agrícola Mouraoense). Na seqüência, um conjunto de simulações serão executadas a fim de estudar: influência da concentração da alimentação e do modo de injeção de vapor direto nos resultados operacionais do stripper.
2 _ MATERIAL E MÉTODOS
2.1 _ O stripper
Uma etapa de extrema importância na destilação da miscela é a operação de stripping do óleo. Esta tem a função de finalizar a destilação da miscela, ou seja, deve haver a separação do baixo teor de hexano que ainda permanece no óleo após a operação de evaporação. O equipamento utilizado neste tipo de operação é denominado de stripper e é mostrado no esquema da Figura 1.
A Figura 1 mostra o esquema de um stripper típico utilizado numa planta de obtenção de óleo bruto de soja, que é uma coluna dividida em compartimentos denominados de estágios e com estruturas específicas para propiciar o fluxo e o contato dos materiais. A miscela concentrada em óleo forma a corrente líquida que entra na parte superior e passa por uma estrutura que facilita sua distribuição dentro do equipamento. Esta corrente ao atingir o estágio forma um pequeno reservatório de líquido, que possibilita o contato entre o vapor direto e a miscela. Este vapor entra no fundo do estágio, que é o distribuir no seio do líquido de modo a ocorrer o stripping do hexano presente no óleo. Após o reservatório atingir um certo nível, o líquido flui ao estágio inferior através de uma tubulação instalada externamente, conforme mostra a Figura 1. Este fluxo da corrente líquida ocorre até que o óleo bruto saia no fundo do stripper. A corrente de vapor, composta de vapor d'água e hexano, flui de maneira contracorrente em relação à corrente líquida e passa de um estágio inferior para o superior através de tubulações também instaladas externamente, conforme mostra a Figura 1. Assim, pode-se observar que o fluxo de correntes, de estágio a estágio, neste equipamento se diferencia um pouco da forma tradicional. Ele se diferencia na medida em que a corrente de vapor que sobe e a corrente de líquido que desce não se cruzam no fundo do estágio uma vez que este é fechado, como mostra o esquema da Figura 1.
O stripper é um equipamento bem específico e para maximizar a recuperação do hexano, KEMPER [6] recomenda que a miscela na entrada tenha uma concentração em massa de óleo de 97% ou mais e uma temperatura superior a 107ºC; o vapor deve ser superaquecido com temperatura de 177ºC e uma pressão em torno de 4kgf/cm2.
O óleo destilado no stripper deve ser testado a fim de verificar se ainda contém traços de hexano. O teste é feito pelo "ponto de fulgor", que é definido como o ponto mais baixo da temperatura em que os vapores saídos do líquido, misturados com o ar, e em contato com uma chama, queimam na superfície desse líquido. A Tabela 1 mostra o ponto de fulgor (ºC) do hexano no óleo de soja e a sua respectiva concentração (ppm).
2.2 _ Modelo matemático do stripper
A Figura 2 mostra o esquema do stripper que se pretende modelar. Este se apresenta como uma coluna dividida em estágios onde o vapor entra em contato direto e em contracorrente com a miscela concentrada, arrastando, assim, o hexano e formando uma corrente de gases que sai no topo; no fundo sai a corrente de óleo bruto praticamente livre do hexano.
Os modelos matemáticos que representam o stripper são baseados em balanços de massa e de energia, equações de equilíbrio e relações empíricas.
Existem na literatura, vários algorítmos para a solução do sistema de equações resultante da modelagem prato a prato de processos de separação. Neste trabalho, pretende-se utilizar o algoritmo do método SR (Sum-Rates) de Burningham-Otto, segundo HENLEY & SEADER [2]. A escolha se justifica uma vez que ele é um método de solução rigoroso e se adequa muito bem ao sistema, pois existe uma diferença significativa de volatilidade entre os componentes do sistema proposto.
O stripper é dividido em estágios e opera continuamente com as correntes de líquido e vapor fluindo, no seu interior, de forma contracorrente. As seguintes hipóteses são assumidas: I) contato entre a fase líquida e a fase de vapor dentro da coluna é perfeito; II) o vapor que sai de cada prato está em equilíbrio com o líquido de saída; III) o vapor d'água é considerado o inerte no processo; IV) o vapor d'água usado dentro da coluna não se condensa; e V) não existem perdas de calor para o meio ambiente.
A modelagem matemática que será utilizada na análise desse extrator é baseada na modelagem geral teórica para separações, em estágio de equilíbrio, já apresentada na literatura por vários autores HOLLAND [4], HENLEY & SEADER [2]. Evidentemente, que no decorrer do desenvolvimento do modelo matemático específico para o caso em estudo serão necessárias algumas adaptações e elas serão feitas.
O esquema do estágio de equilíbrio j é mostrado na Figura 3.
Como se observa pela Figura 3, entrando no estágio j tem-se quatro correntes especificadas de corrente Fj, de corrente Lj-1, de corrente Vj+1 e de corrente Gj. A corrente Fj é a vazão molar de alimentação da mistura com fração molar zh,j do hexano, com temperatura TFj, com pressão PFj e a correspondente entalpia molar HFj; a pressão de alimentação é assumida ser maior ou igual à pressão do estágio Pj. A corrente Lj-1 é a vazão molar do líquido que sai do estágio situado acima do estágio j com fração molar xh,j-1 do hexano, com temperatura Tj-1, com pressão Pj-1 e a correspondente entalpia molar Hj-1. Assumido-se que a pressão Pj-1 é menor ou igual à pressão do estágio j. A corrente Vj+1 é a vazão molar dos gases que saem do estágio situado imediatamente abaixo do estágio j com fração molar yh,j+1, com temperatura Tj+1, com pressão Pj+1 e com a correspondente entalpia Hvj+1 A corrente Gj é a vazão do vapor d'água injetado direto no estágio, com temperatura TGj e com pressão PGj. Saindo do estágio j têm-se duas correntes especificadas como corrente Lj e como corrente Vj. A corrente Lj é a vazão molar de líquido com composição xh,j, com temperatura Tj, com pressão Pj e com a correspondente entalpia molar HLj; a corrente Vj é a vazão molar do gás com fração molar yh,j de hexano, com temperatura Tj, com entalpia molar HVj e pressão Pj a qual está em equilíbrio com a pressão do líquido.
As equações que relacionam todas as variáveis as-sociadas com o estágio teórico j são explicitadas a seguir.
Equação do balanço de massa para o hexano por estágio:
Equação do balanço de energia por estágio:
Equação de equilíbrio para o hexano por estágio:
sendo Kh,j a constante de equilíbrio.
O somatório das frações molares dos componentes em cada estágio
sendo c o número de componentes da mistura.
A equação do balanço material total pode ser usada no lugar das equações (2.4) e (2.5). Ela pode ser resultante da combinação destas duas equações e a com a equação (2.1), somando-se para todos c componentes a partir do estágio 1 até o estágio j. Assim,
A Figura 4 esquematiza um modelo físico de N estágios de uma configuração semelhante ao stripper, a fim de aplicar das equações dos balanços para se obter uma solução para o problema da separação do hexano do óleo de soja, conforme o algoritmo de Burningham-Otto.
Uma coluna de N estágios, contracorrente, como mostrado na Figura 4 é representada matematicamente por N(2c + 3) equações, as quais relacionam N(3c + 9) variáveis. Se N e todos Fj, zh,j, TFj, PFj, Pj, Gj, TGj, PGj são especificados, o modelo é representado por N(2c + 1) equações algébricas simultâneas com N(2c + 1) incógnitas, as quais são todos xh,j, yh,j, Lj, Vj e Tj onde as equações do balanço de massa, de energia e as equações de equilíbrio são equações não lineares. Evidentemente, outras variáveis podem ser especificadas substituindo as incógnitas mencionadas acima. Independentemente das variáveis especificadas, forma-se um conjunto de equações não lineares que deve ser resolvido mediante técnicas de caráter iterativo.
Existe na literatura diversos métodos que podem ser empregados na solução de balanços de massa por componente, relações de equilíbrio e balanços de energia que ocorrem na modelagem extratores tipo stripper. Em geral, estes procedimentos fazem uso do particionamento de equações em conjunto com a equação de iteração e/ou técnicas de linearização de Newton-Raphson. Entretanto, a escolha do método depende da diferença de volatilidade entre os componentes presentes na mistura a ser separada. Então, uma vez que a volatilidade entre os componentes da mistura em estudo é grande, opta-se pela utilização do método SR (Sum-Rates) para a análise do stripper, em estudo.
Este método foi proposto inicialmente por SUJATA [9] e melhorado por Burningham e Otto, segundo HENLEY & SEADER [2] mediante a introdução no método do algoritmo da matriz tridiagonal. Esta resulta da forma modificada das equações (2.1) e (2.2) quando estas são transformadas em outras equações utilizando-se de Tj e Vj como variáveis de iteração. As equações modificadas são lineares em relação às frações molares das correntes líquidas. Este conjunto de equações pode ser resolvido pelo algoritmo de Thomas, segundo HENLEY & SEADER [2].
A obtenção das equações modificadas resulta da substituição da equação (2.2) na equação (2.1) com o objetivo de eliminar y, assim como substituindo a equação (2.6) na equação (2.1) para eliminar L. O resultado destas substituições é o conjunto de equações para o hexano e para cada estágio como está escrito a seguir.
sendo:
com xi,0 = 0 e VN+1 = 0 como indicado em Figura 4. O valor de Kh,j, para sistemas que operam a baixa pressão, pode ser estimado pela seguinte equação:
sendo gh,j o coeficiente de atividade do hexano na solução líquida no estágio j e P0h,j a pressão de vapor do hexano puro no estágio j. Este coeficiente é obtido conforme PARAÍSO, ANDRADE & ZEMP [8] e seu valor, considerado constante ao longo do equipamento, consta na Tabela 2, a seguir.
O valor desta pressão de vapor pode ser estimado pela equação de Antoine, conforme HIMMELBLAU [3],
Sendo a pressão Pºh,j expressa em mmHg e Tj em K.
Se as equações do balanço de massa forem agrupadas por componente, elas podem ser particionadas e colocadas como uma série de Cm matrizes tridiagonais onde a variável de saída para cada matriz equação é xh,j em todos os N estágios do stripper.
As constantes Bj e Cj para cada componente depende somente das variáveis de iteração T e V, uma vez que se considera que o valor da constante de equilíbrio, K, não varie significativamente com a composição ao longo da coluna.
A resolução da matriz tridiagonal gera os valores das composições do componente i em cada estágio j, ou seja, os valores de xh,j. Porém, os valores obtidos não são normalizados nesta etapa, mas utilizados diretamente para gerar novos valores de Lj através da aplicação da seguinte equação
sendo que a equação (2.14) recebe a denominação de equação sum-rates e k é o número de iterações realizadas.
Os valores de Lj(k) são obtidos a partir dos valores de Vj(k) calculados através da equação (2.6). Os correspondentes valores de Vj(k+1) são obtidos do balanço material global, o qual é derivado do somatório da equação (2.1) sobre os c componentes, combinando o resultado com as equações (2.3) e (2.4) e somando aquele resultado sobre todos estágios j até N para se obter a seguinte equação:
Os valores de xi,j necessitam ser normalizados uma vez que eles deverão ser utilizados posteriormente durante a iteração. Isto é feito através da seguinte equação:
Os correspondentes valores de yh,j são calculados através da equação (2.3)
Neste ponto encerra-se o cálculo das composições de saída de cada estágio para o hexano tanto na fase líquida como na fase de vapor, xh,j e yh,j na iteração k, passando-se, então, para o cálculo do novos valores das temperaturas dos estágios, Tj.
Os valores de Tj são obtidos pela solução simultânea do conjunto de equações do balanço de energia para os N estágios. As temperaturas estão embutidas nas correspondentes entalpias específicas das correspondentes correntes de líquido e vapor. Em geral, estas entalpias são não-lineares em relação à temperatura. Portanto, um procedimento de solução iterativa é necessário. Neste caso, o método de Newton-Raphson, segundo HENLEY & SEADER [2], pode ser aplicado com sucesso.
A aplicação desse método para determinar os novos valores de Tj é baseado na equação (2.2), cuja equação recursiva para o método de Newton-Raphson é a seguinte:
sendo:
As derivadas parciais dependem das correlações de entalpia que são utilizadas. Neste trabalho considera-se que a entalpia é independente da composição, portanto, usar-se-á para os cálculos de HLj de HVj e suas respectivas derivadas, as seguintes equações:
sendo cph e cpo os calores específicos do hexano líquido e do óleo de soja, respectivamente, e Tr a temperatura de referência.
sendo th o calor latente de vaporização do hexano na temperatura de referência Tr, e cphy e cpya os calores específicos do hexano e da água no estado de vapor.
As N equações formadas a partir da equação (2.17) formam uma equação matriz tridiagonal que é linear em relação a DTj. A solução da equação matriz gera o conjunto de correções DTj(k). Novos supostos valores de Tj são determinados da seguinte forma:
sendo f um fator escalar de atenuação que é útil quando os valores supostos inicialmente não estão razoavelmente próximos dos valores verdadeiros. Geralmente, f é tomado como unitário. Quando todas correções DTj(k) estão próximas de zero, deve-se procurar um critério para determinar se a convergência foi alcançada. Neste caso, uma maneira de verificar tal convergência é através do seguinte critério:
Substituição sucessiva é geralmente usada no processo de iteração, ou seja, valores de Vj e de Tjgerados a partir das equações (2.15) e (2.24), respectivamente, durante uma iteração são usados diretamente para iniciar a próxima iteração.
2.3 _ Simulação da operação de stripping
2.3.1 _ Desenvolvimento do algoritmo
Aqui, apresentam-se os procedimentos adequados para análise de desempenho da operação de stripping da miscela.
Como já foi abordado anteriormente, fixando-se todas as vazões de alimentação Fj ao lado de sua concentração zh,j e condições termodinâmicas de pressão PFj e temperatura TFj, as vazões de vapor direto Gj e suas condições termodinâmicas TGj e PGj, a pressão de estágio Pj e o número de estágios N; determina-se os valores de xh,j, yh,j, Lj, Vj e Tj. Para tanto, as seguintes etapas são necessárias: I) supõe-se um valor inicial para Tj e Vj; II) calcula-se a pressão de vapor do hexano e a constante de equilíbrio; III) calcula-se os valores das concentrações do hexano, xh,j, resolvendo-se equação matriz (2.7); IV) calcula-se novos valores de Lj a partir da equação sum-rates (2.14) e dos novos valores de Vj utilizando a equação (2.15); V) normaliza-se os valores de x para cada estágio através da equação (2.16) e calcula-se os valores de y correspondentes; VII) calcula-se os novos valores de Tj a partir da equação (2.2) em conjunto com as equações do método de Newton-Raphson; e, VIII) se necessário, faz-se a iteração, conforme mostra a Figura 5.
O programa computacional derivado do algoritmo da Figura 5 é composto de três elementos: entrada de dados, programa principal e saída de dados.
A entrada de dados se compõe basicamente do coeficiente de atividade do hexano na solução e daqueles valores considerados fixos na condição de funcionamento do stripper. Esses dados são: vazão de alimentação do estágio, Fj em moles/hora; fração molar do hexano na alimentação, zh,jem moles de hexano/moles de miscela; e, pressão do estágio, Pj em mmHg.
O programa principal é o elemento que realiza seqüencialmente os cálculos e as iterações necessárias.
A saída de dados compõe-se dos seguintes dados: as composições molares do hexano no topo e no fundo do stripper, yh e xh respectivamente, em moles de hexano/moles de miscela; a temperatura de cada estágio, Tj em ºC; e, a vazão de vapor, Vjem moles/hora.
2.3.2 _ Dados operacionais do stripper
A fim de cumprir o objetivo proposto inicialmente nesta seção, fez-se um levantamento de dados operacionais na indústria da COAMO. Utilizando-se ora somente os dados operacionais ora em conjunto com os dados simulados, em diferentes situações, foi possível obter vários resultados que podem servir de base para um programa de redução de custos numa indústria de obtenção de óleo de soja. O stripper da COAMO apresenta uma configuração geral semelhante à da Figura 6 mostrada a seguir.
Conforme mostra a Figura 6, o stripper da COAMO é composto de dois estágios e opera sob pressão na faixa de 450 a 600 mmHg. Este equipamento recebe a miscela, que vem do sistema de evaporação com uma concentração, em termos de fração mássica, em torno de 98% de óleo e temperatura em torno de 120ºC, com a função de tornar o óleo de soja praticamente livre de hexano. A miscela entra no seu topo de forma pulverizada e flui para baixo entrando em contato direto com o vapor d'água. Este contato possibilita a desorção do hexano, de tal forma, que o óleo sai no fundo com uma concentração mássica de hexano em torno de 1300ppm. O vapor d'água entra, em cada estágio, superaquecido a uma pressão e a uma temperatura de aproximadamente 4,0kgf/cm2 e de 180ºC, respectivamente.
É importante também ressaltar, como já foi mostrado acima, que o stripper da fábrica da COAMO possui dois estágios cada qual com sua entrada de vapor d'água. No entanto, na atualidade a entrada de vapor está sendo feita apenas no primeiro estágio. O argumento dos operadores da fábrica, para tal, é que a injeção de vapor no segundo estágio possibilita a saída de água juntamente com o óleo no fundo do extrator; ainda segundo os operadores isto deve ser evitado uma vez que irá prejudicar o desempenho da operação seguinte, denominada de secagem de óleo. Entretanto, a não injeção de vapor nesse segundo estágio aumenta significativamente o consumo de vapor e diminui a recuperação do solvente como será mostrado, a seguir, através dos resultados simulados.
A Tabela 2 apresenta os operacionais obtidos na indústria da cooperativa COAMO, os quais serão comparados com os resultados obtidos através do processo de simulação a ser realizado.
Os dados de vazão e de concentração presentes na Tabela 2 coletados na COAMO estão expressos em unidades de massa. Para transformá-los em unidades molares utiliza-se o peso molecular do hexano puro, que é igual a 86,17g/mol e o peso molecular do óleo como sendo a média dos pesos moleculares dos ácidos graxos presentes no óleo, ou seja, 795,6g/mol.
3 _ RESULTADOS E DISCUSSÃO
Os resultados obtidos com as simulações do stripper serão apresentados de duas maneiras. A primeira delas se apresenta como a mais relevante: a comparação entre os valores de saída das vazões, das concentrações e das temperaturas obtidos através do modelo e aqueles coletados na indústria da cooperativa COAMO. A segunda maneira, se refere à verificação do comportamento de determinadas variáveis provenientes de alterações impostas nas condições de operação do processo.
3.1 _ Comparação dos resultados do modelo com os resultados operacionais do stripper.
Os resultados obtidos através do modelo são comparados com os dados operacionais obtidos na indústria. A fim de realizar esta análise da forma mais fiel possível, a simulação foi realizada introduzindo vapor direto apenas no primeiro estágio e utilizando-se os dados de entrada conforme é especificado no algoritmo da Figura 5. A Tabela 3 mostra esses dados.
A Tabela 3 mostra que embora, a COAMO não efetua medidas da vazão de vapores que sai no topo do stripper e, nem do consumo de vapor d'água direto utilizado, a empresa avalia o desempenho de seu processo somente através da vazão mássica de óleo bruto, do teor de hexano que sai junto com o óleo e da temperatura do óleo, que saem no fundo do stripper. Assim, somente estes puderam ser comparados diretamente com aqueles obtidos pelo modelo. A comparação mostra uma boa aproximação entre os dados que foram coletados e os valores obtidos através do modelo aplicado. Pode-se observar que o valor da fração mássica do hexano é de 0,00136 ou 1360ppm na indústria e de 0,001413 ou 1413ppm através do modelo, ou seja, a diferença entre o valor real e o valor simulado apresenta um erro relativo de apenas 3,9%; além disso, esses valores estão de acordo com os padrões aceitos para o hexano residual, conforme mostra a Tabela 2. Com relação à temperatura do óleo na saída, pode-se observar que o seu valor real na indústria é de 114,0ºC na indústria, e de 113,6ºC através da simulação, portanto, com um erro de apenas 0,4%. Com relação ao dado de consumo de vapor d'água direto no stripper, apesar de não existir o dado real da indústria, alguns comentários podem ser ditos. O consumo obtido através do modelo foi de 24kg de vapor por tonelada de soja esmagada; este valor é 60% do valor estimado por JONGENEELEN [5] é de 40kg de vapor por tonelada de soja processada em todas opearações de Destilação da Miscela. Isto indica que a COAMO está gastando mais vapor que o recomendado, uma vez que a operação de evaporação, que também compõe a destilação da miscela, consome 34kg de vapor por tonelada de soja esmagada, conforme PARAÍSO [7]. Cabe ressaltar, ainda que a simulação foi efetuada conforme está operando atualmente a coluna da indústria da COAMO com a introdução de vapor somente no primeiro estágio; como será visto posteriormente, no tópico 4.3, esta forma de operar traz ineficiência para o sistema tanto na recuperação de hexano como no consumo de vapor.
3.2 _ Influência da concentração da alimentação nos resultados operacionais do stripper
O objetivo principal aqui é mostrar o efeito que a concentração de hexano na entrada afeta os resultados de saída do equipamento. A fim de realizar esta análise utilizou-se os dados da Tabela 2 e aplicou-se o modelo para dois valores de concentração, em massa: 2% e 1,2%. Os resultados obtidos são mostrados na Tabela 4.
A Tabela 4 mostra que a redução na concentração do hexano na alimentação provoca os seguintes efeitos relevantes nos dados de saída do stripper:
-
uma redução significativa do consumo de vapor direto na operação com a redução da concentração do hexano na alimentação passando de 960kg/h para 522kg/h, e mesmo com esta redução, a quantidade de hexano residual no óleo, nos dois casos, se manteve próximos, em torno de 1400ppm;
-
as vazões mássicas no topo do
stripper diminuiu, passando de 1100kg/h para 601kg/h, assim como a vazão no fundo que passou de 7358kg/h para 7346kg/h. Estas reduções se explicam devido à menor quantidade de hexano no sistema, bem como o menor consumo de vapor direto; e
-
a temperatura no
stripper aumentou, passando de 113,6ºC para 115,7ºC. Este aumento pode ser justificado, uma vez que a quantidade de hexano evaporada no segundo caso foi inferior ao primeiro caso, com isso o resfriamento do sistema foi menor.
3.3 _ Efeito do modo de injeção de vapor direto nos resultados operacionais do stripper
Aqui, pretende-se verificar os efeitos produzidos nos resultados operacionais desse equipamento quando se introduz o vapor direto nos dois estágios. Para tal, introduz-se vapor direto nos dois estágios sendo que as vazões G1 e G2 variam cada uma de 0 a 960,0 kg/h, os demais dados solicitados pelo modelo são como os da Tabela 2. Os resultados obtidos através das simulações são mostrados na Tabela 5.
Uma análise dos resultados apresentados na Tabela 5 mostra que a melhor condição operacional do stripper é a introdução integral do vapor necessário apenas no último estágio desse equipamento. A recuperação do hexano foi superior a todas as outras situações simuladas, especialmente àquela em que se introduz o vapor apenas no primeiro estágio. Além disso, a temperatura no segundo estágio manteve-se um valor próximo ao valor real de 113ºC, e até apresentou um certo declínio para o valor de 112,9ºC, o que é benéfico para a qualidade do óleo produzido.
4 _ CONCLUSÕES
O objetivo deste trabalho foi apresentar um modelo capaz de analisar as condições de operação stripping na tarefa de recuperação do hexano ainda presente no óleo de soja após a evaporação.
Os resultados obtidos mostraram que o modelo proposto é bastante coerente e pode ser perfeitamente utilizado como ferramenta no processo de análise e otimização da operação de stripping do hexano presente no óleo.
5 _ REFERÊNCIAS BIBLIOGRÁFICAS
6 _ AGRADECIMENTOS
À COAMO, que disponibilizou os dados operacionais, situada no município de Campo Mourão-PR.
Recebido para publicação em 28/06/2002. Aceito para publicação em 28/12/2004 (000915).
- [1] BARBOSA, M.Z.; FREITAS, S.M. de; FRANCA, T.J.F. Considerações sobre os desafios da cadeia de produção de óleo de soja no Brasil, Óleos & Grãos, p. 44, Setembro/Outubro, 1998.
- [2] HENLEY, J.E.; SEADER, D.J. Equilibrium-Stage Separation Operations in Chemical Engineering, John Wiley & Sons, New York, 1981.
- [3] HIMMELBLAU, D.M. Engenharia Química: Princípios e Cálculos. Trad. Jossyl de Souza Peixoto, Sindicato Nacional dos Editores de Livros, Rio de Janeiro, 1984.
- [4] HOLLAND, C.D. Fundamentals and Modeling of Separation Processes, Prentice-Hall, New Jersey, 1975.
- [5] JONGENEELEN, H.P.J. Energy Conservation in Solvent Extration Plants, JAOCS Journal of the American Oil Chemists Society, v. 53, p. 291, June, 1976.
- [6] KEMPER, G.T. Minimização da perda de solvente, Grãos & Óleos, p. 22, Setembro/Outubro, 1998.
- [7] PARAÍSO, P.R. Modelagem e Análise do Processo de Obtenção do Óleo de Soja, 2001, Campinas. Tese de Doutorado _ FEQ/UNICAMP, 212 p.
- [8] PARAÍSO, P.R.; ANDRADE, C.M.G.; ZEMP, R.J. Destilação da Miscela I: Modelagem e Simulação da Evaporação do Hexano. Ciênc. Tecnol. Aliment., v. 23, n. 3, p. 459-467, set.-dez., 2003.
- [9] SUJATA, A.D. Fundamentals and Modeling of Separation Processes in: HOLLAND, C. D. Prentice-Hall, New Jersey, 1975.
Datas de Publicação
-
Publicação nesta coleção
08 Set 2005 -
Data do Fascículo
Mar 2005
Histórico
-
Aceito
28 Dez 2004 -
Recebido
28 Jun 2002