Cálculo de Flash2

1
UNIVERSIDADE FEDERAL DE UBERLNDIA

FACULDADE DE ENGENHARIA QUMICA
CURSO DE GRADUAO EM ENGENHARIA QUMICA
Desenvolvimento de uma Biblioteca Eletrnica de

Modelos
Fernando Rosa Gomes
Uberlndia / MG
2009
UNIVERSIDADE FEDERAL DE UBERLNDIA

FACULDADE DE ENGENHARIA QUMICA
CURSO DE GRADUAO EM ENGENHARIA QUMICA
Desenvolvimento de uma Biblioteca Eletrnica de Modelos
Fernando Rosa Gomes
Monografia de graduao apresentada

Universidade Federal de Uberlndia como
parte dos requisitos necessrios para
aprovao
da
disciplina
de Projeto de
Graduao do curso de Engenharia Qumica.
Uberlndia / MG
2010
BANCA EXAMINADORA DA MONOGRAFIA DA DISCIPLINA PROJETO DE

GRADUAO DE FERNANDO ROSA GOMES APRESENTADA FACULDADE DE
ENGENHARIA QUMICA DA UNIVERSIDADE FEDERAL DE UBERLNDIA, EM
12/01/2010.
BANCA EXAMINADORA:
___________________________________
Prof. Lus Cludio Oliveira Lopes
Orientador- FEQUI/UFU
___________________________________
Prof. Dra. Mria Hespanhol Miranda Reis
FEQUI/UFU
___________________________________
Prof. Cludio Roberto Duarte
FEQUI/UFU
AGRADECIMENTOS
Ao orientador Prof. Dr. Lus Cludio Oliveira Lopes que sempre acreditou no meu
potencial e na minha capacidade em desenvolver este projeto de graduao.
s amigas Juliana Maria e Lvia Tizzo, companheiras de reunio, com quem
compartinhei momentos de dvidas e sucessos deste trabalho. Alm disso, foram momentos
que proporcionaram nos conhecer melhor e aumentar os laos de amizade.
Aos amigos, Diego Souza e Karen Valente, que foram meus companheiros de reunio
antes da minha ida para a Frana.
E principalmente aos meus familiares e amigos que sempre estiveram comigo e que
entenderam os momentos que estive desenvolvendo este e que no passei com eles. Este
apoio, mesmo que indiretamente, colaborou para a execuo da obra e meu xito.
O mestre na arte da vida faz pouca distino entre o

seu trabalho e o seu lazer, entre a sua mente e o seu
corpo, entre a sua educao e a sua recreao, entre
o seu amor e a sua religio. Ele dificilmente sabe
distinguir um corpo do outro. Ele simplesmente persegue sua viso de excelncia em tudo que faz,
deixando para os outros a deciso de saber se est
trabalhando ou se divertindo. Ele acha que est
sempre fazendo as duas coisas simultaneamente.
(Texto Zen-Budista)
SUMRIO
INTRODUO .................................................................................................................. 1
REVISO BIBLIOGRFICA............................................................................................ 3
2.1.
IMPORTNCIA DE FONTES DE DADOS TERMODINMICOS ............................ 3
2.1.1.
DIPPR (Design Institute for Physical Properties)........................................................ 3
2.1.2.
FYDAT ........................................................................................................................ 4
2.2.
CAPE-OPEN ................................................................................................................... 5
2.3.
A REGRA DE FASES DE GIBBS E GRAU DE LIBERDADE ................................... 7
2.4.
CLCULOS DE FLASH SISTEMA DE NICO ESTGIO DE EQUILBRIO....... 9
2.4.1.
Derivaes do Modelo Flash ..................................................................................... 10
2.4.1.1.
Estratgias para os clculos de flash ...................................................................... 13
2.4.1.2.
Sequncia de clculo Flash TP ............................................................................... 14
2.4.2.
2.5.
Mtodo Inside-Out para os Clculos de Flash ........................................................... 15

PROJETO DE REATORES NO-ISOTRMICOS .................................................... 18
2.5.1.
Balano de energia e a Primeira Lei da Termodinmica ........................................... 19
2.5.2.
CSTR (Continuous-Stirred-Tank-Reactor) ................................................................ 20
2.5.3.
PFR (Plug-Flow Tubular Reactor) ............................................................................. 22
RESULTADOS E DISCUSSO ...................................................................................... 24
3.1.
FHASH ISOTRMICO MULTICOMPONENTE........................................................ 24
3.2.
FLASH ADIABTICO MULTICOMPONENTE........................................................ 39
3.3.
REATOR CSTR ADIABTICO MULTICOMPONENTE ......................................... 42
3.4.
REATOR PFR ADIABTICO MULTICOMPONENTE ............................................ 44
3.5.
CASO ESPECIAL DE UM SEPARADOR MEMBRANA DE BIOGS................ 45
3.6.
OBTENO DE PARMETROS FSICOS E TERMODINMICOS ...................... 46
CONCLUSO .................................................................................................................. 46
SUGESTES PARA PRXIMOS TRABALHOS ................................................................. 47

REFERNCIAS BIBLIOGRFICAS ..................................................................................... 47
APNDICES ............................................................................................................................ 48
APNDICE A - CDIGO SCILAB UTILIZADO PARA SIMULAR OS CASOS DE FLASH
ISOTRMICO, PONTO DE BOLHA E PONTO DE ORVALHO. ........................................ 49
APNDICE B CDIGO IMPLEMENTADO NO SCILAB QUE UTILIZA AS
FERRAMENTAS DO COCO. ................................................................................................. 55
APNDICE C CDIGO SCILAB DE CLCULO DO FLASH ISOTRMICO

UTILIZANDO O MTODO DE NEWTON E A PLATAFORMA DO COCO. .................... 57
APNDICE D CDIGO SIMPLIFICADO PARA O CLCULO DO FLASH
ISOTRMICO QUE UTILIZA SOMENTE COMANDOS DO COCO. ................................ 60
APNDICE E CDIGO SCILAB UTILIZADO PARA IMPLEMENTAR O CASO DE
FLASH ADIABTICO NO COCO. ........................................................................................ 60
APNDICE F CDIGO SCILAB UTILIZADO PARA A DEFINIO DA UNIDADE DE
CSTR ADIABTICO. ............................................................................................................. 61
APNDICE G CDIGO SCILAB UTILIZADO PARA A DEFINIO DA UNIDADE
PFR ADIABTICO. ................................................................................................................ 63
APNDICE H.1 CDIGO PRINCIPAL SCILAB DA UNIDADE DE SEPARAO DE
BIOGS. .................................................................................................................................. 66
APNDICE H.2 ARQUIVO ADICIONAL SCILAB PARA A IMPLEMENTAO DA
UNIDADE DE SEPARAO DE BIOGS. .......................................................................... 68
APNDICE I ARQUIVO SCILAB DE LEITURA DE PROPRIEDADES DAS
SUBSTNCIAS CONTIDAS NO DYPPR E FYDAT. ......................................................... 69
LISTA DE FIGURAS
Figura 1 - Janela principal do ambiente de simulao do COCO-COFE. .................................. 6
Figura 2 - Separao contnua e com um nico estgio de equilbrio: (a) vaporizao flash
(flash adiabtico com vlvula, flash isotrmico sem vlvula com T especificada); (b)
condensao parcial (anlogo para flash isotrmico quando T especificado). ...................... 10
Figura 3 - Unidade Flash. ......................................................................................................... 11
Figura 4 - Localizao da seo de criao de pacotes termodinmicos (ConfigureTEA) e
caixa de seleo de pacotes termodinmicos criados previamente dos pacotes de substncias
qumicas. ................................................................................................................................... 28
Figura 5 - Adicionando componentes a um pacote termodinmico. ........................................ 29
Figura 6 - Seleo das propriedades e equaes de estado....................................................... 30
Figura 7 - Insero de unidade de operao Scilab no COCO. ................................................ 32
Figura 8 - Edio de uma unidade de Operao no COCO. ..................................................... 33
Figura 9 - Configurao da corrente de alimentao. ............................................................... 34
Figura 10 - Edio da unidade de operacional Scilab, aproximando a aparncia da unidade a
ser simulada. ............................................................................................................................. 35
Figura 11 - cone de teste do cdigo Scilab implementado no COCO..................................... 35
Figura 12 - Tela de simulao do COCO mostrando o caso do Flash isotrmico I- Mtodo de
Newton, para o qual no obteve sucesso. ................................................................................. 36
Figura 13 - Tela de simulao do COCO mostrando o caso do Flash isotrmico II. ............... 38
Figura 14 - Tela de simulao do COCO mostrando o caso do Flash isotrmico. ................... 39
Figura 15 - Impresso da tela principal do COCO mostrando os resultados simulados para o
caso do flash adiabtico. ........................................................................................................... 41
caso de um CSTR adiabtico. ................................................................................................... 43
caso do flash adiabtico. ........................................................................................................... 44
Figura 18 Fluxograma e planilhas de resultados obtidos para o caso de Separao de Biogs.
.................................................................................................................................................. 46
LISTA DE TABELAS
Tabela 1 - Algoritmo para o CSTR. ......................................................................................... 21
Tabela 2 - Algoritmo para um PFR adiabtico, com reao elementar reversvel em fase
gasosa ( A B ) e perda de presso desprezvel. .................................................................... 23
Tabela 3 - Simulao obtida no software HYSIM para o caso flash isotrmico para a
separao de hidrocarbonetos. .................................................................................................. 25
Tabela 4 - Resultado gerado na simulao do cdigo implementado no Scilab de clculo do
flash isotrmico. ....................................................................................................................... 27
Tabela 5 - Relatrio da simulao da unidade Scilab Flash isotrmico II, resgatando os
dados gerados pelo prprio COFE............................................................................................ 37
Tabela 6 - Relatrio da simulao da unidade Flash isotrmico, gerado atravs da utilizao
da unidade operacional Flash do COFE. .................................................................................. 38
Tabela 7 - Valores para a simulao de um flash adiabtico apresentado por Seader & Henley
(2006). ...................................................................................................................................... 40
Tabela 8 - Valores para a simulao do flash adiabtico no COCO utilizando a equao de
estado de Soave Redlich Kwong Modificada. .......................................................................... 42
Tabela 9 - Valores da simulao do CSTR adiabtico no COCO. ........................................... 43
Tabela 10 - Valores da simulao do PFR adiabtico .............................................................. 45
RESUMO
Este trabalho visou o desenvolvimento de uma base (biblioteca) de dados
composta por modelos dinmicos de baixa e mdia complexidade de atividades nas reas
de separao de fluidos e reatores qumicos, utilizando ferramentas computacionais que
fornecessem, numa linguagem computacional livre e de fcil acesso, resultados mais
prximos possveis daqueles que se encontram na realidade industrial e citados na
literatura. Alm disso, este foi conduzido em dois caminhos diferentes, porm
complementares. O primeiro caminho dedicou-se na implementao no Scilab de cdigos
de busca de dados termodinmicos das substncias fsico-qumicas, atravs da
decodificao binria de dois softwares, o FYDAT e o DIPPR, especficos para esta
atividade e disponveis comercialmente. O segundo caminho foi a utilizao das
plataformas Scilab e COCO (que utiliza o padro CAPE-OPEN e Inside-Out),
verificando a viabilidade de se trabalhar simultaneamente com ambos, atravs da
implementao de problemas tpicos da Engenharia Qumica. A vantagem de utilizar
estes softwares est na interface grfica amigvel do COCO e tambm na possibilidade
de implementar problemas utilizando mtodos que produzem os resultados desejveis. Os
resultados mostraram que a gama de possibilidades grande e que nenhum ou pequenos
desvios relativos dos valores apresentados na literatura podem ocorrer, e essa preciso
depende da equao de estado e das correlaes termodinmicas de predio dos
parmetros utilizadas nos clculos.
Palavras-chave: Biblioteca de Modelos, Simulao, CAPE-OPEN, Linguagem Livre,
Inside-Out.
ABSTRACT
This work aimed at the development of a database (library) composed for dynamic
models of low and average complexity of activities in the areas of separation of fluids
and chemical reactors, using computational tools which supply, in a free computational
language and of easy access, results more nearby to those that is found in the industrial
reality and cited in the literature. Moreover, this work was lead in two different ways,
however complementary. The first one dedicated in the implementation of search Scilab
codes of thermodynamic data of physicist-chemistries substances, through the binary
decoding of two softwares, the FYDAT and the DIPPR, specific for this activity and
available commercially. The second way was the use of the Scilab platforms and
COCONUT (that it uses standard CAPE-OPEN and Inside-Out), verifying the viability to
work simultaneously with both, through the implementation of Chemical Engineering
typical problems. The advantage to use these softwares is the COCONUT friendly
graphical interface and also the possibility to implement problems using methods that
produce the results desirable. The results had shown that the gamma of possibilities is
great and none or small relative shunting lines of the values presented in literature can
occur, and this precision depends on the equation of state and the thermodynamic
prediction correlations of the parameters used in the calculations.
Keywords: Library of Models, Simulation, CAPE-OPEN, Free Language, Inside-Out.
1 INTRODUO
A IMPORTNCIA DA ANLISE DE PROCESSOS
Para se obter um bom entendimento de processos em que o Engenheiro Qumico
ir atuar, necessrio analis-los partindo de suas formas mais simples e isoladas do
todo. Ou seja, necessrio que o engenheiro possua uma base bem fundamentada dos
fenmenos que ocorrem no determinado processo para, depois conseguir inferir as
influncias que este pode sofrer dos fatores externos. Mas medida que o processo exige
uma maior complexidade, j que de suma importncia representar da melhor forma o
que se verifica na realidade, torna-se impossvel o desenvolvimento de modelos e
resolues das equaes de forma manual. Nesse sentido, o uso de computadores, mais
precisamente, de simuladores, de grande ajuda, tanto na questo da velocidade de
processamento quanto na criao de um banco de dados permanente, disponvel sempre
que necessrio e que no exige recomear o procedimento toda vez que mudar os dados
analisados.
Historicamente, os processos qumicos evoluram de uma escala e unidades
pequenas onde eram geralmente operadas no modo batelada ou semi-contnuo. Matriaprima e energia eram relativamente abundantes. As regulamentaes governamentais em
segurana, efluentes e emisses eram pouco rgidas. As margens de lucro eram
relativamente altas e atrativas e a tecnologia, em termos de ferramentas analticas, era
primitiva. Nos ltimos anos considerveis mudanas tm sido adotadas. Primeiro pela
competio e depois pelos custos de matria-prima e energia tm criado uma convincente
necessidade por processos eficientes, levando a um maior grau de integrao entre
energia e matria no design de processos. Cada vez que introduzido um reciclo de
material ou um trocador de calor, criam-se novas interaes entre as variveis do
processo, no qual aumenta o nmero de relaes matemticas que podem ser
consideradas simultaneamente.
importante notar que o sistema como um todo pode exibir caractersticas
superiores e alm daquelas que se consideram quando se utiliza apenas seus componentes
individuais. Em Engenharia Qumica, um sistema bem conhecido a associao de um
reator qumico exotrmico como um simples trocador de calor. Ambas as unidades
apresentam estabilidades. No momento em que o trocador de calor conectado ao reator
para pr-aquecer a alimentao e para resfriar o produto, um comportamento bastante

diferente observado para o sistema como um todo, pois vrias solues possveis
podem ocorrer, mas algumas so solues estveis enquanto outras so instveis.
Portanto, quanto mais unidades forem sendo acopladas, mais complexa ser a sua
resoluo.
A gama de simuladores que hoje existe fruto da necessidade de predizer os
fenmenos que ocorrem na natureza e nos diversos segmentos de indstrias, sem correr o
risco de haver desperdcio de matria-prima e gasto desnecessrio de mo-de-obra, alm
de preservar a segurana das pessoas envolvidas. O que difere um simulador do outro ,
na maioria das vezes, a interface grfica e a forma com que o algoritmo foi transcrito.
Esta ltima o ponto de partida de muitos estudos, j que sempre mais vivel trabalhar
com programas que gerem respostas em menor perodo de tempo e tambm com menor
esforo computacional. Porm, no basta um software bem estruturado, se somente o
programador/autor do mesmo consegue entend-lo, necessrio que o simulador esteja
numa linguagem acessvel, para que outras pessoas consigam aprender e capacitar outras
tantas. O desenvolvimento de modelos codificados em uma linguagem computacional de
fcil acesso permite grande impacto dessas avaliaes no ambiente educacional e de
pesquisa, onde esses modelos podem ser usados em laboratrios e ambiente de testes para
implementaes variadas.
Alm disso, necessrio que o engenheiro possua tambm um conhecimento
sobre pacotes termodinmicos responsveis por predizer a energia requerida pelo sistema,
o momento do equilbrio de fases e o tamanho dos equipamentos. Ou seja, a eficincia de
uma simulao depende dos mtodos de predio das propriedades das misturas ideais e
no-ideais, como o volume especfico, entalpia, entropia, fugacidades e atividades
juntamente com seus coeficientes, todos como funo da temperatura, presso e
composio das fases. Outro parmetro importante a constante de Equilbrio, K, que
mostra a composio limite atingvel para uma mistura de dois ou mais componentes, em
outras palavras, quando se atinge o equilbrio qumico.
A modelagem de processos qumicos permite que representaes conceituais de
processos sejam avaliadas em diversas reas de interesse da Engenharia Qumica,
destacando-se a otimizao de processos, a avaliao de segurana, a parada e partida de
plantas e o controle de processos. Com esse intuito, esse projeto visa o desenvolvimento
de uma base (biblioteca) de dados composta por modelos dinmicos de baixa e mdia
complexidade que implementam atividades nas reas de: separao de fluidos, reatores
qumicos, reatores biotecnolgicos etc. e consiste na investigao e desenvolvimento do

modelo, implementao dos mesmos em uma linguagem computacional livre e de fcil
acesso utilizando ferramentas de anlise de processos qumicos e de simulao
computacional.
2 REVISO BIBLIOGRFICA
2.1. IMPORTNCIA DE FONTES DE DADOS TERMODINMICOS

Toda atividade significativa em qumica ou engenharia qumica requer quantificao
dos parmetros e composio envolvidos nos processos. Um simples computador e at
mesmo um complexo simulador impossvel quantificar sem o conhecimento das
propriedades fsicas dos presentes componentes e suas misturas. Para muitos, esta pode ser
uma atividade tediosa e cansativa. Alm disso, complicam-se quando se tem dados de
substncias menos usadas de pssima credibilidade e, tambm quando se comete erros na
transferncia de unidades. Para um clculo em engenharia qumica ter credibilidade
imprescindvel haver preciso nas propriedades dos dados.
Equaes de estado como aquelas de SOAVE (1972), PENG ROBINSON (1976), e
JACOBSEN et al. (1995), so cada vez mais usadas para correlacionar propriedades
termofsicas. Independente das propriedades particulares do modelo usado, o uso efetivo de
propriedades termofsicas na prtica industrial, dependente de algoritmos rpidos e robustos
para a simulao de processos. por esses motivos que programas computacionais so
criados, ficando com o trabalho de calcular as propriedades das substncias e, com isso, a
dedicao fica restrita somente na resoluo dos problemas.
2.1.1. DIPPR (Design Institute for Physical Properties)

O banco de dados DIPPR 801 um extenso banco de dados de componentes puros
rigorosamente avaliados de compostos qumicos importantes industrialmente. Este contm um
total de 1891 compostos at hoje. Valores recomendados podem ser encontrados para 48
propriedades de cada componente no banco de dados, isto , dados crus da literatura assim
como valores recomendados e criticamente avaliados para cada uma das 34 propriedades
constantes. Dados crus da literatura e coeficientes de equaes para 15 propriedades

dependentes da temperatura so tambm includas para cada composto. Alm disso,
referncias, notas e cdigos de qualidade podem ser encontrados para cada propriedade
termodinmica.
2.1.2. FYDAT
O banco de dados FYDAT um software livre o qual pode ser ampliado inserindo
novas substncias. Este no universalmente usado, mas destinado s tcnicas
computacionais para propsitos prticos e pedaggicos. O FYDAT pode ser utilizado como
um direto dilogo interativo para o clculo de propriedades ou como uma biblioteca usada em
programas de aplicao escrita em Visual Basic 5.0.
Este banco de dados fornece as seguintes propriedades fsicas:
Massa molar;
Constantes crticas (presso, temperatura, volume);
Fator acntrico;
ndice de refrao;
Momento dipolo;
Ponto de ebulio normal;
Ponto de fuso;
Calor de fuso;
Ponto de ebulio (dependente da presso);
Ponto de orvalho;
Calor de evaporao;
Densidade (fase vapor e lquida);
Viscosidade dinmica (fase vapor e lquida);
Capacidade trmica (fase vapor e lquida);
Condutividade trmica (fase vapor e lquida);
Tenso superficial;
Presso de vapor;
Funes termodinmicas (entalpia,entropia, energia interna);
Coeficiente de Joule-Thompson;
Logaritmo do coeficiente de fugacidade, entre outras.

Estas propriedades podem ser fornecidas tanto em moles ou quilogramas. E todos os
valores so listados em unidades SI menos para as excees notificadas.
2.2. CAPE-OPEN
Os padres CAPE-OPEN so normas uniformes para interface de componentes de

softwares de modelagem de processos, desenvolvidas especialmente para o design e operao
de processos qumicos. As normas so pblicas, com multiplataformas, uniformes e
livremente disponveis para realizar download. (http://www.colan.org/index-3.html).
O CAPE-OPEN descrito em um conjunto de documentos formais abrangendo reas
como operaes unitrias, propriedades fsicas e resolvedores numricos. Na prtica, este
permite linkar componentes de outras partes, como os nichos de pacotes de propriedades
fsicas ou modelos de operaes unitrias, para ser usada no modo plug and play em
ferramentas de softwares de modelagem de processos comerciais. Com isso, remove a maior
barreira para a construo de melhores modelos de processos clssicos em reas especficas.
O CAPE-OPEN Laboratories Network (CO-LaN) uma reconhecida organizao de
responsveis pelo gerenciamento do padro CAPE-OPEN (CO). A CO-LaN ajuda a facilitar a
implantao de interfaces do padro CO em softwares comerciais tornando uma simulao em
uma realidade prtica.(http://www.colan.org/Communication/COLAN_2003_Leaflet.pdf).
O ambiente de simulao escolhido para se trabalhar nessa monografia foi o COCO
(CAPE-OPEN to CAPE-OPEN), que utiliza o padro CAPE-OPEN. Este foi criado pelo
grupo amsterCHEM, formado por experts no padro citado. O COCO consiste dos seguintes
componentes:
COFE CAPE-OPEN Flowsheet Environment: uma interface intuitiva grfica para
o usurio construir fluxogramas. Exibe as propriedades das correntes, trabalha com
converso de unidades e fornece facilidades durante a plotagem de grficos.
TEA-COCOs Thermodynamics for Engineering Applications, baseado na biblioteca
termodinmica de cdigos e inclue um banco de dados de mais de 150 substncias
qumicas comumente usadas. O pacote exibe mais de 100 mtodos de clculos de
propriedades com suas derivadas numricas e analticas.
COUSCOUS o CAPE-OPEN Unit Operations Simple contm divisores,

misturadores, trocadores de calor, bombas e reatores alm de outras unidades de
operao.
CORN O pacote CAPE-OPEN Reaction Numerics especifica algum tipo de cintica
ou reao em equilbrio. Simples unidades de reao, como reatores de converso,
CSTRs e PFRs podem usar este pacote.
Figura 1 - Janela principal do ambiente de simulao do COCO-COFE.

O padro CAPE-OPEN pode interagir com outros softwares e com isso, ser possvel
utilizar a linguagem de programao especfica daquele software. Isso possvel com o
Microsoft Office Excel, Matlab e Scilab. Na pgina do grupo amsterCHEM,
http://www.amsterchem.com, existem guias de como operar o CAPE-OPEN nos trs
softwares acima citados e vice-versa. Mais especificamente para este trabalho, a interao
utilizada foi entre o CAPE-OPEN e o Scilab.
As duas formas disponveis para realizar a interao so atravs da utilizao dos
pacotes Scilab CAPE-OPEN Thermo Import e Scilab CAPE-OPEN Unit Operation. O
primeiro permite importar do CAPE-OPEN pacotes de propriedades termodinmicas e fsicas
para o Scilab, atravs de um pacote de comandos, e ento apresentar os clculos das mesmas
propriedades assim como clculos de equilbrio termodinmico. O pacote define todas as
interfaces de comunicao com os pacotes de propriedades CAPE-OPEN e este no
implementa um equipamento termodinmico sozinho. Para utilizar estes pacotes, so
necessrias as verses 5.02 (ou maior) do Scilab assim como a verso 1.1 do CAPE-OPEN.
(http://www.amsterchem.com/scilabthermo.html)
A segunda forma de interao, o Scilab CAPE-OPEN Unit Operation, uma
implementao de uma unidade de operao na qual os clculos so realizados na linguagem
Scilab. Para inici-lo, necessrio iniciar o ambiente de simulao escolhido que utiliza o
padro CAPE-OPEN, no caso deste, o ambiente de simulao COCO. Aps, inserir um
Scilab CAPE-OPEN Unit Operation, definir as portas de alimentao e produto. Para
editar a unidade de operao, insere-se o cdigo Scilab na seo corresponde este, e as
correntes de entrada e sada devem ser especificadas de acordo com a linguagem especfica de
resgate das mesmas no padro CAPE-OPEN.
(http://www.amsterchem.com/scilabunitop.html)
2.3. A REGRA DE FASES DE GIBBS E GRAU DE LIBERDADE
A descrio de um sistema de estgio nico no equilbrio fsico envolve variveis

intensivas, na qual so independentes do tamanho do sistema, e variveis extensivas, no qual
dependem do tamanho do sistema. Variveis intensivas so temperatura, presso, e
composies de fase (fraes molares, fraes mssicas, concentraes, etc.). Variveis
extensivas incluem massa ou moles e energia para um sistema batelada, e vazo mssica ou
molar e taxa de transferncia de energia para um sistema contnuo.
Independente se somente variveis intensivas ou ambas variveis intensivas e
extensivas so consideradas, somente uma pequena parcela das variveis so independentes;
quando essas so especificadas, todas as outras se tornam fixas. O nmero de variveis
independentes a varincia, ou o nmero de grau de liberdade, F, para o sistema.
A regra de fase de J. Willard Gibbs, a qual se aplica somente para as variveis
intensivas no equilbrio, afirma que:
F = CP +2
(1)
onde C o nmero de componentes e P o nmero de fases no equilbrio. A equao (1)

derivada contando, no equilbrio fsico, o nmero de variveis intensivas e o nmero de
equaes independentes que relata essas variveis. O nmero de variveis intensivas, V,
V = C .P + 2
(2)
onde o 2 refere-se a presso e temperatura do equilbrio, enquanto o termo CP o nmero

total de variveis de composio (por exemplo, fraes molares) para os componentes
distribuos ao longo das fases de equilbrio P. O nmero de equaes independentes, E,
relacionando as variveis intensivas
E = P + C (P 1)
(3)
no qual o primeiro termo, P, refere necessidade do somatrio das fraes molares ou

mssicas ser 1 para cada fase e o segundo termo, C(P
P-1) refere ao nmero de equaes de K
independentes da forma geral
Ki =
fraao molar de i na fase (1)

fraao molar de i na fase ( 2)
onde (1) e (2) refere s fases em equilbrio. Para duas fases, h C expresses independentes
deste tipo, trs fases, 2C; para quatro, 3C; e assim por diante.
Assim, o Grau de Liberdade, que o nmero de variveis intensivas (V) menos o
nmero de equaes (E), pode ser definido, como combinao das equaes (2) e (3):
F = V E = (CP + 2 ) - [P + C(P-1)] = C-P + 2

Usando a Regra de Fases de Gibbs, pode-se notar que Ks no so variveis, mas so
funes termodinmicas que dependem de variveis intensivas. A Regra limitada porque ela
no trata de correntes enviadas para o estgio de equilbrio nem de variveis extensivas
usadas no design ou anlise em operaes de separao. No entanto, a regra de fases pode ser
estendida para aplicaes em processos, adicionando a corrente de alimentao e variveis
extensivas, e equaes independentes adicionais relatando as variveis de alimentao,
extensivas, e intensivas j consideradas pela regra. (SEADER; HENLEY (2006))
2.4. CLCULOS DE FLASH SISTEMA DE NICO ESTGIO DE EQUILBRIO
Em programas de simulao de processos, a determinao de um sistema de separao

flash representa o mais frequentemente envolvido e o mais bsico dos conjuntos de clculos.
Este requerido para determinar o estado de alguma corrente do processo aps uma
transformao fsica ou qumica. Isso ocorre aps a adio ou remoo de calor, uma
mudana na presso ou uma carga na composio devido reao. Os clculos do flash
podem resultar em um comportamento de fase complexo, no qual inclui mltiplas fases e
tambm fases slidas.
Um flash um estgio nico de destilao de equilbrio na qual uma alimentao
parcialmente vaporizada para dar um vapor mais rico em componentes mais volteis, do que o
lquido remanescente. Na Figura 2a, uma alimentao lquida aquecida sobre presso e
flasheada adiabaticamente atravs da vlvula para um vaso de menor presso, resultando na
criao de uma fase vapor que separada do lquido remanescente no mesmo vaso. Se a
vlvula omitida, um lquido de baixa presso pode ser parcialmente vaporizado no
aquecedor e ento separado em duas fases no vaso de flash. Alternativamente, uma
alimentao de vapor pode ser resfriada e parcialmente condensada, com separao de fases,
como na Figura 2b, para dar um lquido que mais rico no componente menos voltil. Em
ambos os casos, se o equipamento bem projetado, o vapor e o lquido que deixam o vaso
esto em equilbrio.
10
Figura 2 - Separao contnua e com um nico estgio de equilbrio: (a) vaporizao flash
(flash adiabtico com vlvula, flash isotrmico sem vlvula com T especificada); (b)
condensao parcial (anlogo para flash isotrmico quando T especificado).
A menos que a volatilidade relativa seja muito grande, o grau de separao alcanado
entre os dois componentes em um estgio nico de equilbrio pobre. Portanto, um flash
(parcial vaporizao ou condensao parcial) usualmente uma operao auxiliar usada para
preparar correntes para um processo seguinte. Tipicamente, a fase vapor enviada para o
sistema de separao de vapor, enquanto a fase lquida enviada para o sistema de separao
de lquidos. (SEADER; HENLEY (2006))
2.4.1. Derivaes do Modelo Flash
Para ilustrar a discusso dessa seo, considere a operao de separao de fase

representada na Figura 3.
11
Vapor
Feed
Liquid
Flash isotrmico
Figura 3 - Unidade Flash.
O modelo desenvolvido para esta unidade de separao baseado nos fluxos molares
para os NC componentes i nas correntes de alimentao, vapor e lquida, fi, vi, li,
respectivamente. Assume-se que o estado da alimentao completamente definido, ou seja,
conhece-se a vazo, frao molar (zi) e a entalpia de entrada. Definindo, as fraes molares
como xi=li/(ili) e yi=vi/(ivi) obtm-se o mnimo conjunto de balano de massa:
f i = vi + l i
i=1,...NC
(4)
equaes de equilbrio:
i (l,T)f i 0 (T, P)x i = (v,T)Py i
i=1,...NC
(5)
e um balano de entalpia:
FH f ( f ,T , P ) + Q = VH v ( v ,T , P ) + LH l ( l ,T , P )
(6)
o qual nos d (2 NC+1) equaes para as (2 NC+3) variveis, vi, li, T, P e Q. Tem-se portanto
um grau de liberdade igual dois para especificar para o problema flash.
No entanto, quando uma fase desaparece, um modelo derivado dos balanos de massa
em fluxos molares conduz a composies indefinidas para as condies de ponto de orvalho e
de bolha. Alm disso, j que as relaes de equilbrio de fase no-linear so dependentes da
composio, ser desenvolvido um modelo de flash um pouco diferente em termos das vazes
e fraes molares. Seguindo a mnima descrio acima, o balano de massa para a unidade
dada por:
zi F = Vyi + Lxi
i=1,...NC
(7)
E o balano de entalpia porduz
FH f + Q = VH v ( y ,T , P ) + LH l ( x ,T , P )
(8)
12
Definindo ainda o valor K:

0
K i = i (x,T)f i (T, P)xi / (y,T)P i=1,...NC
(9)
Portanto tem-se (3 NC+5) variveis (yi, xi, L, T, P, Q, V) e somente (3 NC+1)

equaes at agora. Notar que no foi especificado nenhuma condio ainda das condies de
xi e yi (por exemplo, que a soma das fraes molares 1). Interessantemente, esta escolha
necessita ser feita cuidadosamente j que falsas razes so introduzidas mesmo com algumas
escolhas bvias.
J que nem a frao molar de lquida nem de vapor especificada pode-se incluir um
balano de massa global:
F =V + L
i=1,...NC
(10)
Agora considere que o mais simples dos casos onde T e P so especificados. Este
dissocia o balano de entalpia e permite que Q seja calculado uma vez que o balano de massa
resolvido. Combinando a equao do balano total de massa com as expresses de balano
de massa dos componentes e de equilbrio chega-se nas seguintes relaes para as fraes
molares:
xi =
zi
1 + (K i 1)
V
F
yi =
K i zi
1 + (K i 1)
V
F
(11)
Neste momento, necessita-se de uma especificao adicional em qualquer um dos

conjuntos de fraes molares para obter um modelo com o requerido grau de liberdade igual
dois.
Considere as duas escolhas bvias: xi =1 ou yi=1. Par a primeira escolha tem-se:
Fz i
x = F + (K
i
=1
- 1)V
(12)
e o modelo flash trivialmente satisfeito para todo problema flash se for estabelecido que yi =
zi e V=0. Similarmente, se for usado yi=1 encontra-se:
F.K i .z i
=1
i - 1)V
y = F + (K
i
(13)
e o modelo flash trivialmente satisfeito para todo problema flash se for estabelecido que yi
= zi e V=F.
13
Claramente, qualquer uma das equaes conduz pseudo-solues (na composio de

alimentao) que no esto relacionadas com a verdadeira soluo do problema flash. Para
eliminar as solues triviais considera-se uma especificao alternativa de Rachord e Rice
(1952). Tomando a diferena de xi=1 e yi=1 tem-se:
x y =0
i
(14)
Nota-se que esta nova especificao, junto com o balano de massa total, ainda conduz
s corretas especificaes nas fraes molares. Aplicando estas condioes s relaes das
fraes molares tem-se:
F (K i - 1)z i
=0
i
)
1
V
i
y x = F + (K
i
(15)
e v-se que as pseudo-razes no podem resolver esta equao. De fato, xi=zi ou xi=zi so
solues permitidas somente se a (bem apropriada) condio que Ki = 1 e se for uma mistura
azeotrpica. (BIEGLER et al. (1999))
2.4.1.1.
Estratgias para os clculos de flash
O modelo flash pode ser dado concisamente por:
z i F = Vyi + Lxi
i = 1,...NC
y i = K i xi
i = 1,...NC
0
K i = i (x,T)f i (T, P)xi / (y,T)P i = 1,...NC

(16)
F =V + L
x y =0
i
FH f + Q = VH v ( y, T , P ) + LH l ( x, T , P )
e agora dois graus de liberdade podem ser especificados. Enquanto muitas alternativas so
possveis para os clculos do projeto, os clculos de flash so geralmente resolvidos pelos
graus de liberdade escolhidos entre as variveis V/F, Q, P e T, como pode ser visto abaixo,
nas combinaes possveis dessas variveis:
T, P -
Flash isotrmico
V/F = 0, P - Temperatura do ponto de bolha

V/F = 1, P - Temperatura do ponto de orvalho
T, V/F=0 - Presso do ponto de bolha
14
T, V/F=1 - Presso do ponto de orvalho

Q=0, P
Flash adiabtico
Q, P
Flash no adiabtico
V/F, P
Percentagem de vaporizao no flash
O mais simples dos casos dado pelo flash (P,T) j que este no requer iterao para o
balano de energia. Para este caso, o problema flash pode ser resolvido pela sequncia de
clculos, apresentada na seo 2.4.1.2.
2.4.1.2.
Sequncia de clculo Flash TP
1. Para um fixo zi (estando certo que zi=1) e F, especifica-se T e P. Continuar se estiver

entre o ponto de bolha e ponto de orvalho. (Para no-idealidade dependente da
composio, fornecer uma estimativa inicial para xi e yi);
2. Chute V/F;
3. Calcule K i = i (x,T)f i 0 (T, P)x i / (y,T)P ;
4. Calcular x i = z i /(1 + (K i - 1)V/F) e yi =Kixi;
5. Avaliar a relao implcita (V/F) =
xi
yi. Se (V/F) zero (ou com uma
pequena tolerncia), PARE. Siga para o passo 6;

6. Reavalie o chute para V/F e volte ao passo 3.
Especificaes TP so mais comuns para misturas em que todos os componentes tm
pontos de ebulio numa faixa bem estreita (valores bem prximos), como o benzeno e o
tolueno. Aqui V/F pode variar enre zero e um, com uma pequena faixa de temperatura. Este
caso comum para misturas separadas por destilao. Por outro lado, para misturas com uma
faixa larga de ebulio (como a gua e ar) a especificao TP na sequncia de clculos do
flash funciona pobremente porque a temperatura de equilbrio varia bastante para pequenas
mudanas em V/F. Estas misturas comumente separadas por absoro e as especificaes
(V/F, T) e (V/F, P) so usadas. Seno, o algoritmo similar sequncia de clculos apresenta
anteriormente. Alm disso, para estes casos, note que o balano de entalpia no necessrio
no loop de iterao.
Finalmente, quando a especificao no calor de entrada, Q, feita (como em um flash
adiabtico), ento um balano de entalpia imposto e necessita ser incorporado ao algoritmo
15
do flash. Geralmente o balano de entalpia tratado chutando a temperatura, e resolvendo o

flash TP no loop interno. A entalpia ento calculada, combinada ao calor de entrada, e a
temperatura reavaliada no loop externo. Esta sequncia de clculos consome muito mais
tempo. Alternativamente, todas as equaes do modelo flash podem ser resolvidas
simultaneamente usando o mtodo de Newton ou Broyden. Com esta simultnea
aproximao, ambas as misturas de pontos de ebulio prximos e mais distantes podem ser
conduzidos por um caminho simples. No entanto, para todos estes mtodos, rotinas
termodinmicas no-ideias necessitam ser chamadas frequentemente e isso aumenta os gastos
computacionais.
2.4.2. Mtodo Inside-Out para os Clculos de Flash
As sequncias de clculo flash desenvolvidas acima sofrem de duas desvantagens:

Elas so projetadas para misturas de altos ou baixos pontos de ebulio e se
apresentam pobres para os casos opostos;
Elas requerem frequentes chamadas para avaliar as funes termodinmicas noideais, especialmente quando o balano de energia necessita ser incorporada nos
clculos do flash.
Preocupados com esses problemas, Boston e Britt (1978) desenvolveram um algoritmo
inside-out que acelera enormemente a soluo dos problemas de flash. No loop externo,
esta aproximao combina as equaes das propriedades fsicas no-ideais e expresses
simplificadas para K e entalpias e ento usa estas expresses para resolver as equaes de
flash no loop interno. A soluo destas equaes ento usada para atualizar as expresses
simplificadas e o procedimento termina uma vez que as expresses combinam aqueles atuais
no-ideais no loop externo.
Para ilustrar as vantagens do algoritmo inside-out, considere o tipo de flash PQ com as
equaes de flash dadas anteriormente. Boston (1980) alm disso sugeriu as seguintes
simplificaoes para o loop interno:
16
Ki = i Kb
ln(K b ) = A + B (1/T - 1/T * )
H' v = C + D (T - T * )
(17)
H' l = E + F (T - T * )
Onde os parmetros A, B, C, D, E, F e i so livres para combinarem com as expresses noideais de K e entalpias (Hv e Hl computadas em base mssica). Kb uma mdia de K que
baseada na importncia do K componente i. i representa as volatilidades relativas, e Hv e Hl
so as entalpias de gs ideal com temperatura de referncia T*. Para manusear ambas as
misturas com substncias com pontos de ebulio estreitos e largos, Boston e Britt definiram
uma varivel artificial de iterao, R = Kb/(Kb+L/V). Esta varivel toma o domnio de
temperatura ou V/F para misturas mencionadas, respectivamente, e elimina a necessidade de
separar algoritmos para estes sistemas. Este ocorre porque R pode agora variar bastante ambas
para grandes mudanas em T (wide boiling) e L/V (narrow boiling). Neste ponto, uma vez
fixados os parmetros (A, B, C, D, E, F e i) para o loop externo, pode-se derivar as seguintes
relaes atravs da substituio das equaes de flash e expresses simplificadas:
z i F = f i = Vy i + Lxi
(18)
Usando yi =Kixi e definindo Ki =iKi, tem-se:
f i = (VK i + L )xi = ( iVK b + L )xi
(19)
Dividindo por (VKb+L) e substituindo por R chega em:
f i / (VK b + L ) = ( iVK b + L )xi / (VK b + L )

f i / (VK b + L ) = ( i R + 1 - R )xi
(20)
Define ento o novo conjunto de variveis:
pi = xi (VK b + L ) = xi L / (1 - R ) = f i / ( i R + 1 - R )xi
(21)
Note que pis so determinados somente de R e quantidades especificadas no loop

externo. Do somatrio e das equaes de equilbrio obtem-se:
17
L = (1 - R ) pi
V =F-L
K b = ( pi
p )
i
(22)
xi = pi / (VK b + L )
y i = i K b xi
1
T = ((ln K b A) / B + 1 / T * )
Usando R como a varivel de iterao, o clculo do flash completado checando o

balano de entalpia simplificado. O algoritmo de Boston-Britt pode ser resumido pela
seguinte sequncia de clculos:
1. Iniciar A, B, C, D, E, F e i;
2. Chutar R;
3. Resolver para pi, Kb,T, L, V, xi e yi usando as equaoes acima;
4. Converter os fluxos para a base mssica e avaliar as entalpias mssicas simplificadas
para a equao de balano;
(R ) = H ' f +Q / F' +(L' / F' )(H' l (x,T, P ) H ' v ( y,T, P )) H ' v ( y,T, P )
5. Se (R) est dentro da tolerncia estipulada, ir para o passo 6. Seno, modifique o
chute para R e volte para o passo 3;
6. Na primeira passagem, obtem-se os novos valores de A, B, C, D, E, F e i comparando
com as expresses no-ideais. Para os prximos ajuste somente A, C, E e i usando o
mtodo de Broyden para conciliar estes parmetros com as expresses no-ideais.
Boston e Britt preferem a base mssica para o balano de entalpia para evitar
insensibilidade de R (atravs de L/F) quando (Hv-Hl) prximo de zero em termos molares.
Este algoritmo converge muito mais rapidamenteque os algoritmos desenvolvidos
anteriormente e tem sido incorporado como o algoritmo padrao para o flash em simuladores
de processos comerciais.
Para demonstrar este algoritmo, Boston e Britt resolveram uma grande variedade de
sistemas no-ideais incluindo sistemas de ponto de ebulio prximos e afastados, e com
Wilson, UNIQUAC, NRTL, e equao de estado. Tpica experincia nestes exemplos tem
menos que seis interaes externas (onde avaliaes das propriedades fsicas so requeridas).
Finalmente, experimentos numricos tm mostrado que este algoritmo geralmente pode ser
18
resolvido com a composio dependente de K embora as expresses simplificadas no sejam

funes de xi. Para casos altamente no-ideais, no entanto, Boston (1980) sugere uma
modificao que leva em conta o K simplificado dependente da composio e torna o
algoritmo mais robusto. (BIEGLER et al. (1999))
2.5. PROJETO DE REATORES NO-ISOTRMICOS
Reatores qumicos, particularmente para processos contnuos, so geralmente

projetados para envolver mltiplas fases (por exemplo, vapor, lquido, slido reagente e
catalizador slido), diferentes geometrias (tanques agitados, tubulares, transporte atravs da
membrana), e vrios regimes de transferncia de momento, calor e massa (fluxos viscosos,
turbulentos, conduo, radiao, difuso e disperso). So vrias configuraes, envolvendo
diferentes combinaes de seus atributos, que na tentativa de desenvolver modelos de reatores
generalizados tem-se deparado com um sucesso limitado.
A maioria dos simuladores de processo fornece quatro tipos de modelos de reatores,
incluindo:
(1) Um modelo estequiomtrico que permite a especificao das converses e exteno da
reao para uma ou mais reaes especficas;
(2) Um modelo de mltiplas fases no equilbrio qumico, onde a aproximao para o
equilbrio para reaes individuais podem ser feitas;
(3) Um modelo cintico para um reator CSTR que assume perfeita mistura das fases
homogneas (lquido e vapor); e
(4) Um modelo cintico para um reator PFR, para fases homogneas e assumindo
nenhuma disperso.
Estes modelos ideais so usados nos primeiros estgios da sntese de processos,
quando os detalhes do projeto do reator so menos importante, porm os efluentes do reator e
a conduo de calor so necessrios.
Os modelos de reatores ideais so substituidos por modelos comerciais de
significantes detalhes. Para este propsito, todos os simuladores fornecem facilidades para a
insero de modelos gerados pelo usurio e tambm facilitam a estimao dos parmetros dos
modelos cinticos por regresso no-linear. (SEIDER et al. (2003))
19
2.5.1. Balano de energia e a Primeira Lei da Termodinmica
Um sistema qualquer poro limitada do universo, estacionrio ou em movimento,

que escolhida para a aplicao de vrias equaes termodinmicas. Para um sistema
fechado, no qual nenhuma massa cruza as fronteiras do sistema, a variao na energia total do
sistema, d, igual ao fluxode calor para o sistema, Q, menos o trabalho realizado pelo
sistema sobre as vizinhanas, W:
dE = Q W
(23)
no qual os s signficam que Q e W no so diferenciais exatas de uma funo de estado.

J para os sistemas abertos, ou seja, onde massa cruza as fronteiras do sistema, como
nos casos dos reatores de escoamento contnuo citados anteriormente, o balano de energia
pode ser representado pela expresso genrica, representada pela equao (24).
taxa de
taxa de transferncia taxa de transferncia

taxa de
taxa de variao de de energia adicionada de energia perdida
transfern
cia
acmulo

= de calor das trabalho realizado + ao sistema atravs pelo sistema atravs

de energia
pelo sistema sobre

vizinhanas
da massa transferida da massa que sai
no sistema para o sistema suas vizinhanas para o sistema
do sistema
(24)
O balano de energia para o regime no estacionrio para um sistema aberto que tenha
n espcies, cada qual entrando ou saindo do sistema, com suas respectivas vazes molares Fi
(moles de i por tempo), e com suas respectivas energias Ei (joules por mol de i),
n
dE sist & &
= Q W + E i Fi
dt
i =1
E i Fi
entrada
i =1
(25)
saida
De forma mais especfica, no caso do regime estacionrio:

n
& W
& F T C dT H (T ) + T C dT F X = 0
Q
s
A0
Ti 0 i p i Rx R TR p A 0
i =1
(26)
onde:
FA0 fluxo molar da espcie de referncia A;
razo do nmero de moles iniciais (que entram) da espcie i, pelo nmero de moles
iniciais (que entram) da espcie A;
20
Cpi calor especfico da espcie i na temperatura T, podendo o termo
i =1
Ti 0
i C p i dT ser
representado por:
n
i =1
Ti 0
~
i C p i dT = i C p i (T Ti 0 )
i =1
~
Cpi
Ti 0
C p i dT
(Ti 0 T )
H Rx (TR ) + C p dT - a entalpia de reao temperatura T, H Rx (T ) , e H Rx (TR ) a

TR
entalpia de reao temperatura de referncia TR; e

X a converso. (FOGLER, (1999))
2.5.2. CSTR (Continuous-Stirred-Tank-Reactor)
O mais simples dos reatores cinticos o CSTR, no qual o contedo assumido estar
perfeitamente agitado. Com isso, a composio e a temperatura so assumidas uniformes em
todo o volume do reator e igual composio e temperatura do efluente do reator. Um reator
perfeitamente agitado usado geralmente para reaes homogneas em fase lquida. O
modelo CSTR adequado para este caso, desde que a reao ocorra sob condies
isotrmicas ou adiabticas. Embora os clculos envolvam somente equaes algbricas, estas
podem ser no-lineares. Consequentemente, uma possvel complicao que deve ser
considerado a existncia de mltiplas solues, duas ou mais podendo ser estveis.
(SEIDER et al. (2003))
A equao de projeto para um CSTR no qual no h variao espacial na velocidade
de reao
V=
FA 0 X
rA
(27)
Acoplando a equao (27) na equaao (26) e fazendo alguns rearranjos e simplificao

obtm a equao (28)
n
& W
&
Q
s
(T T ) = iC (T T )
X HRx (TR ) + C
p
R
pi
i0
FA0
i=1
(28)
sendo C p os calores especficos mdios ou constantes. Esta equao utilizada para obter o
volume do reator ou a temperatura de operao. Se necessrio, o CSTR ou aquecido ou
21
resfriado por uma camisa de aquecimento ou resfriamento ou por uma serpentina colocada no
interior do reator.
As reaes so frequentemente conduzidas adiabaticamente e usualmente com
aquecimento ou resfriamento das correntes de entrada e sada do vaso de reao. Com a
exceo de processos altamente viscosos, o trabalho feito pelo agitador geralmente pode ser
& na equao (28), o balano de energia pode ser escrito

desprezado. Depois de substituir Q
como
n
UA(Ta T )
~
(T T ) = C
X H Rx (TR ) + C
p
R
i pi (T Ti 0 )
FA 0
i =1
(29)
sendo Ta a temperatura ambiente.

Sob condies de operao adiabtica e trabalho de agitao desprezvel, o balano se
torna
n
~
(T T ) = C
X H Rx (TR ) + C
i pi (T Ti0 )
p
R
(30)
i =1
Abaixo ilustrado um algoritmo para o projeto de um reator no isotrmico

considerando a reao irreversvel de primeira ordem, em fase lquida:
Tabela 1 - Algoritmo para o CSTR.
1. Equao de projeto do CSTR:
V=
FA 0 X
rA
(T.1.1)
2. Lei de velocidade de reao: rA = kC A

com k = Ae
RT
(T.1.2)
(T.1.3)
3. Estequiometria (fase lquida, v=v0):

C A = C A 0 (1 X )
(T.1.4)
4. Combinando, resulta
V=
v0
Ae
RT
1 X
(T.1.5)
CASO A: As variveis X, v0, CA0 e FA0 so especificadas, e o volume do reator, V,

deve ser determinado:
5.a. Encontrar a temperatura, T, para a alimentao. Para o caso adiabtico, resolver
22
~
T + C
X H Rx (TR ) + C
i pi Ti0
p R
i =1
T=
(T.1.6)
~
i C pi + XC p
i =1
Para o caso no adiabtico,

UA(Ta )
T=
~
T + C
X H Rx (TR ) + C
i pi Ti0
p R
FA 0
i =1
UA
FA 0
~
+ C
+ XC
i pi
p
(T.1.7)
i =1
6.a.Calcular k pela equao de Arrhenius

7.a.Calcular o volume do reator, V, a partir da equao (T.1.5).
CASO B: As variveis v0, CA0, V e FA0 so especificadas, e a temperatura de sada,

T, e a converso, X, so quantidades desconhecidas. O procedimento :
5.b.Resolver o balano de energia (adiabtico) para X como uma funo de T:
n
~
C (T T )
i
X EB =
[H
pi
i0
i =1
Rx
(T.1.8)
(TR ) + C p (T TR )]
Para o caso no adiabtico,
X EB
UA(Ta T ) n
~
i C pi (T Ti 0 )
FA 0
i =1
=
(T T )
H (T ) + C
Rx
(T.1.9)
6.b.Resolver a equao (T.1.5) para X como uma funo de T:

X MB =
Ae
1 + Ae
RT
RT
onde =
V
v0
(T.1.10)
7.b.Encontrar os valores de X e T que satisfaam tanto o balano de energia quanto o

balano molar.
(FOGLER, (1999))
2.5.3. PFR (Plug-Flow Tubular Reactor)

Mais complexo, no qual a composio do fluido, escoando como um pisto,
gradualmente muda ao longo do comprimento do reator, sem gradiente de composio ou
23
temperatura na direo do escoamento. Alm disso, as taxas de transferncia de calor e massa

so negligenciveis na direo axial. Ento, o PFR completamente no-agitado, com todos
os elementos tendo o mesmo tempo de residncia no reator. Se o reator opera sob condies
adiabticas ou no-isotrmicas, a temperatura do fluido muda gradualmente ao longo do
reator. (SEIDER et al. (2003))
O balano de energia dado pela equao (26) relaciona a converso em qualquer ponto
temperatura da mistura de reao no mesmo ponto (isto , fornece X como uma funo de
T). Normalmente existe uma quantidade desprezvel de trabalho realizado sobre a, ou pela,
mistura reagente, de forma que, em geral, o termo de trabalho pode ser desprezado no projeto
de reator tubular. Contudo, a menos que a reao seja conduzida adiabaticamente, a equao
(26) ainda difcil de ser resolvida porque em reatores no adiabticos o calor adicionado ou
removido do sistema varia ao longo do comprimento do reator. Este problema no ocorre em
& e
reatores adiabticos, os quais so frequentemente encontrados na indstria. Como o Q
& so iguais a zero, a equao (26) reduz-se a
W
s
X[ H Rx (T )] =
Ti 0
C
i
i =1
pi
dT
(31)
Esta equao pode ser combinada com o balano molar diferencial

dX
FA 0
= rA (X, T ) para obter a temperatura, a converso e os perfis de concentrao ao
dV
longo do comprimento do reator. O algoritmo para resolver PFRs operados adiabaticamente
mostrado na Tabela 2.
Tabela 2 - Algoritmo para um PFR adiabtico, com reao elementar reversvel em fase
gasosa ( A B ) e perda de presso desprezvel.
dX rA
=
(T.2.1)
Balano molar:
dV FA 0
C
Lei de velocidade de reao: rA = k C A B
KC
com k = k 1 exp
E 1 1

R T1 T
K C = K C (T2 )exp
Estequiometria:
H Rx 1 1

R T2 T
Gs, =0, P=P0
(T.2.2)
(T.2.3)
(T.2.4)
24
C A = C A 0 (1 X )
C B = C A0 X
T0
T
T0
T
Combinando,
X T0
rA = kC A 0 (1 X )

K C T
X T0
kC A 0 (1 X )

K C T
dX
=
e
dV
FA 0
(T.2.5)
(T.2.6)
(T.2.7)
(T.2.8)
Balano de energia:
n
~
T + C
X H Rx (TR ) + C
i pi Ti0
p R
T=
i =1
n
~
i C pi + XC p
(T.2.9)
i =1
~
Entrar com os valores dos parmetros k1, E, R, KC(T2), HRx(TR), C pi , C p , CA0, T0, T1, T2 e
P.
Entrar com os valores iniciais X=0, V=0 e os valores iniciais X=Xf e V=Vf.
Definir as vazes iniciais de entrada dos componentes, Fi0.
3 RESULTADOS E DISCUSSO
O objetivo dessa anlise implementar um simples processo qumico usando
diferentes plataformas de simulao e comparando os resultados obtidos entre si e, se
possvel, com aqueles encontrados na literatura, obtidos de dados experimentais previamente
testados ou que so realizados em plantas qumicas reais. No caso deste trabalho foram
implementados casos de flashs isotrmico e adiabtico; um caso de CSTR adiabtico e um
PFR adiabtico. Alm disso, ao final ser mostrado um exemplo de um sistema de separao
de biogs por um sistema de membrana (a unidade implementada) que representa apenas uma
parte de um conjunto de outros equipamentos.
3.1. FHASH ISOTRMICO MULTICOMPONENTE
25
Para o clculo do flash isotrmico multicomponente, vrias formas diferentes de

representar um problema foram testadas de maneira que ao final fosse possvel verificar a
eficincia dos cdigos implementados e tambm das plataformas utilizadas.
Um exemplo de vaso flash isotrmico retirado da literatura, mais precisamente aquele
da simulao no software HYSIM e apresentados no seu Guia do Usurio, pgina 7-110, foi
utilizado para a validao dessa seo: uma corrente de 144 kmol/h a 60F e 600 Psia
alimentada a um vaso flash. Esta formada por uma mistura de hidrocarbonetos, onde os
componentes e suas fraes molares so apresentados na Tabela 3. Os resultados aps
simulao tambm foram descritos na mesma tabela. Alm disso, a alimentao descrita
anteriormente ser a mesma para todos os casos apresentados nessa seo.
Tabela 3 - Simulao obtida no software HYSIM para o caso flash isotrmico para a
separao de hidrocarbonetos.
Corrente
Alimentao
Lquido
Vapor
Unidade
Presso
600
600
600
Psia
Temperatura
60
60
60
Vazo
144
74,3626
69,6374
lbmol/h
Frao molar metano
0,4861
0,1933
0,7988
Frao molar etano
0,1389
0,1435
0,1339
Frao molar propano
0,0694
0,1040
0,0325
Frao molar isobutano
0,0625
0,1067
0,0153
Frao molar n-butano
0,0556
0,0979
0,0104
Frao molar isopentano
0,0486
0,0900
0,0044
Frao molar n-pentano
0,0417
0,0779
0,0030
Frao molar n-hexano
0,0486
0,0929
0,0013
Frao molar n-heptano
0,0278
0,0535
0,0003
Frao molar n-octano
0,0208
0,0403
0,0001
Entalpia
231900,2465
-41325,5303
273225,7809
Btu/h
Primeiramente trabalhou-se com a implementao de um cdigo no Scilab 5.1.1. Este,

que pode ser visto no Apndice A, foi concebido para ser capaz de calcular 3 casos diferentes,
dependendo da necessidade do problema. Estes casos so: (1) clculo de flash isotrmico, (2)
condio de ponto de bolha de uma mistura e (3) ponto de orvalho. A escolha feita inserindo
o nmero do caso desejvel a ser resolvido via teclado, na tela de simulao do Scilab
26
Console, aps a validao do cdigo e clicando em Execute. A partir deste ponto

necessrio inserir os dados solicitados na mesma tela de simulao Console, para que o
problema tenha grau de liberdade zero e possa convergir a uma resposta coerente. Neste
momento se est interessado no cas (1), ou seja, no clculo de um flash isotrmico. O cdigo
implementado considera algumas hipteses que so importantes para o clculo das
composies de equlibrio, como:
Considerou-se que o sistema de mistura comporta-se como um sistema ideal,

por ser composto de espcies da mesma famlia (hidrocarbonetos), embora a
presso no sendo baixa (600 psia ou 41,4 bar);
O modelo termodinmico segue a Lei de Raoult;
So necessrios calcular as variveis K e estes so encontrados dividindo a

presso de saturao de cada componente puro pela presso do sistema;
Para o clculo das presses parciais dos componentes, faz-se uso da Equao
de Antoine e consequentemente, foi necessrio conhecer os valores dos
parmetros dessa equao para cada componente. Os valores dos parmetros
foram obtidos do apndice A do livro The properties of gases and liquids de
Reid e Prausnitz (2001).
Para avaliar a convergncia dos clculos, lanou-se mo do Mtodo de Newton, que

avalia o quanto este mtodo preciso na estimativa das variveis. Para isso, definiu o limite
mximo de iteraes e a tolerncia do desvio entre o valor calculado na iterao atual e aquele
da iterao imediatamente anterior. Foram adotados como limte mximo de iteraes e
tolerncia do desvio os valores 1.107 e 1.10-8, respectivamente. Alm disso, definou-se
tambm qual seria a estimativa inicial para a razo entre a vazo de Vapor e a da
Alimentao, representada pela varivel VF. Esta ser a varivel de avaliao nas iteraes do
Mtodo de Newton.
Com isso, inseriu-se os dados necessrios para a resoluo e gerou-se os seguintes
resultados:
27
Tabela 4 - Resultado gerado na simulao do cdigo implementado no Scilab de clculo do

flash isotrmico.
iterao
V/F
funo
derivada
1
0,5000
0,0318
2,1985
2
0,4855
0,0002
2,1737
3
0,4855
0,0000
2,1736
***************************************************************************
No de iteraes= 3
Valor de VF= 0,485452
A vazo da fase vapor e lquida so, respectivamente (em lbmol/h):
Composio a P= 600 Psia
69,905148
74,094852
*******************************************************************
* componente
composio
vazo
*
*
frao molar
(lbmol/h)
*
*-------------------------------------------------------------------------------------------------*
*
|
x
y
|
liquida
vapor
*
* [ metano]
| 0.149866
0.842486 | 11.104326
58.894074 *
* [ etano]
| 0.154405
0.122466 | 11.440620
8.560980
*
* [ propano]
| 0.115503
0.020534 | 8.558175
1.435425
*
* [isobutano]
| 0.114602
0.007275 | 8.491419
0.508581
*
* [n-butano]
| 0.103812
0.004499 | 7.691922
0.314478
*
* [isopentano] | 0.093073
0.001461 | 6.896259
0.102141
*
* [n-pentano]
| 0.080174
0.000921 | 5.940449
0.064351
*
* [n-hexano]
| 0.094169
0.000300 | 6.977416
0.020984
*
* [n-heptano]
| 0.053982
0.000049 | 3.999795
0.003405
*
* [ octano]
| 0.040414
0.000010 | 2.994472
0.000728
*
****************************************************************
Verifica-se que o Mtodo de Newton foi bastante eficaz, visto que foram necessrios
apenas 3 iteraes para se obter os valores do Estado Estacionrio. O valor final da razo
entre a vazo de Vapor produzida sobre aquela da Alimentao foi 0,485, ou seja, 48,5% da
alimentao vapor e 51,5% lquido. Comparando os valores obtidos nesta simulao com
aqueles apresentados na Tabela 3 verifica-se que a proporcionalidade de cada frao molar de
cada componente nas correntes vapor e lquida se mantm, porm a grandeza de cada um
difere. Este resultado j era esperado, pois se trata de uma simplificao do flash e mtodos
utilizados naquele implementado no HYSIM.
Para que valores mais precisos fossem obtidos seria necessrio, obviamente, aumentar
o grau de detalhamento do cdigo, em outras palavras, detalhar o comportamento
termodinmico da mistura, pois se trata de uma mistura no-ideal (devido condio de
presso e tambm pela variedade de tamanhos das molculas e por haver compostos com
ramificao). Este maior detalhamento viria com um aumento dos esforos tanto mental
quanto computacional.
28
Neste sentido, buscaram-se novos caminhos que diminussem estes esforos e que
pudessem oferecer resultados ainda melhores. As alternativas so aquelas quando se utiliza a
plataforma CAPE-OPEN no Scilab e vice-versa. E o COCO 1.15, foi o software escolhido
para este trabalho que utiliza a plataforma CAPE-OPEN.
Antes de descrever os resultados obtidos da interao entre a plataforma CAPE-OPEN
com o Scilab necessrio inserir um pacote contendo todas as substncias da mistura e o
pacote de propriedades termodinmicas responsveis pelos clculos. extremamente
importante que estes pacotes estejam totalmente criados para que as simulaes feitas com os
mesmos possam ser vlidas. Para isso, os passos a seguir apresentam como so criados e
modificados estes pacotes termodinmicos:
Abre-se a janela de interface ConfigureTEA, que utiliza a verso 1.1 do CAPEOPEN.
Se o pacote de substncias qumicas que se deseja avaliar j existe, clique em

done.
Figura 4 - Localizao da seo de criao de pacotes termodinmicos (ConfigureTEA) e

caixa de seleo de pacotes termodinmicos criados previamente dos pacotes de substncias
qumicas.
29
Caso este no exista, clique em Create template (ver Figura 4). Na janela que
se abre, insere-se um nome para o pacote e uma descrio (opcional) na aba
General. Na aba Compounds, clica-se em Add e selecionam-se as
substncias
desejadas.
Verifica-se
em
Property
Calculations
quais
propriedades termodinmicas e qual modelo das mesmas deseja-se calcular,

como por exemplo, para a Equao de Estado das fases vapor e lquida se o
modelo mais bem representado pelo modelo de Peng Robinson ou Redlich
Kwong enre outras (ver Figuta 5).
Alm disso, pode-se definir a tolerncia do erro, a quantidade mxima de

iteraes a relao das pertubaes da temperatura e presso no sistema, etc.,
inseridas na aba Options.
Existem outras abas que outros dados podem ser alterados de acordo com a
necessidade do problema.
Figura 5 - Adicionando componentes a um pacote termodinmico.
30
Figura 6 - Seleo das propriedades e equaes de estado.
Criado o pacote necessrio escrever um cdigo Scilab que possua alguns

comandos especficos de busca dos dados termodinmicos contidos no pacote
de propriedades citado anteriormente. Os principais comandos utilizados nesta
implementao so aqueles mostrados logo abaixo e todo o cdigo foi descrito
no Apndice B.
Quando utilizada a funo capeOpenPackageManagers( ) esta buscar os Pacotes

Mestres do COCO onde pacotes especficos de substncias qumicas foram definidas pelo
usurio e estes sero apresentados na tela de interface, podendo serem selecionados para se
trabalhar no Scilab. Com este pacote selecionado, pode-se detalhar quais pacotes especficos
de substncias e propriedades sero utilizados, atravs do comando capeOpenPackages(x),
onde x representa o nome recebido para o Pacote Mestre. O cdigo aqui implementado e
aqueles que sero implementados futuramente utilizam os pacotes criados no 'TEA (CAPEOPEN 1.1 e sua busca feita utilizando capeOpenGetPackage( ). E para obter os resultados
simulados utiliza-se a linha de comando
31
[phases,phaseFractions,phaseCompositions,T,P]=capeOpenEquilibrium(handle,X,prop1,val1,
prop2,val2,type), onde:
handle: o pacote escolhido entre os demais do pacote 'TEA (CAPE-OPEN

1.1;
X: composio global;
prop1: a propriedade 1 definida para a realizao dos clculos;
val1: valor da prop1;
prop2: a propriedade 2 definida para a realizao dos clculos;
val2: valor da prop1;
type: string indicando o tipo de soluo desejada. opcional;
phases: quais fases esto presents no equilbrio;
phaseFractions: a frao que cada fase ocupa no equilbrio;
phaseCompositions: a composio molar dos componentes em cada uma das

fases;
T e P so a temperatura e presso no equilbrio, respectivamente.
Ao final da simulao do cdigo no foi possvel obter qualquer resultado, pois

existem algumas inconsistncias no identificadas. Por se tratar de um trabalho inovador erros
como este podem acontecer e, somente com um dispndio de tempo maior so possveis de
serem resolvidos. Mas em momento nenhum se tira o mrito de se construir uma biblioteca de
modelos. O desafio pode ser repassado para futuros pesquisadores e tambm aos usurios dos
cdigos.
A ltima forma de anlise do mesmo processo foi efetuar o caminho inverso
apresentado no segundo teste: utilizou a linguagem Scilab para projetar um vaso Flash
isotrmico no ambiente de simulao do COCO. Neste ambiente j existem simples unidades
operacionais prontas para serem utilizadas, somente necessitando definir o(s) pacote(s)
termodinmico(s) (apresentado anteriormente), as portas de entrada e sada de dados, e
consequentemente, especificarem os campos intrnsecos contidos na porta de entrada. Porm
o objetivo desse trabalho utilizar a interface do COCO somente como apoio para obteno
de propriedades termodinmicas e tambm como uma interface grfica amigvel. A inteno
conseguir atingir um grau de detalhamento e complexidade de um processo, representando-o
da maneira que se aproxime o mximo possvel da realidade, sem a necessidade de utilizar
32
modelos matemticos intrnsecos do software, que podem no ser os mais indicados para o
caso. Porm, como se necessita de um parmetro de comparao, configurou-se em uma
mesma janela de simulao uma unidade operacional com um cdigo Scilab e outra gerada
pelo COCO. As especificaes das entradas e pacotes termodinmicos foram os mesmos para
ambas e, ao final, compararam-se os resultados, analizando os desvios gerados.
Para isso, a unidade operacional que ser implementada com o cdigo Scilab,
inicialmente funciona como uma caixa-preta na qual se pode escrever um cdigo seja qual
for a linguagem escolhida (Scilab, Matlab), que manipula as informaes de entrada e
transforma-as nas informaes de sada do processo.
Os passos para se chegar nesta configurao so os seguintes:
Abre-se a janela de interface COFE 1.15. Seleciona na barra de ferramentas o

cone Insert, Unit operation, Custom e finalmente Scilab Unit
Operation. Confirme em ok.
Figura 7 - Insero de unidade de operao Scilab no COCO.
Para modificar a unidade, clica-se duas vezes sobre esta e na janela que se abrir
seleciona-se Show GUI, na parte debaixo da mesma. Uma nova janela se abre
e nesta possvel definir as portas de entrada e sada (no caso do Flash, uma
entrada e duas sadas), os parmetros e o prprio cdigo. A linguagem utilizada
na seo Scilab no difere muito daquela implementada software de mesmo
nome, como pode ser comprovada no Apndice C.
33
Figura 8 - Edio de uma unidade de Operao no COCO.
Com as portas definidas, a unidade caixa-preta aceitar inserir setas que

representam as correntes, na configurao definida na seo Ports. As
magnitudes e unidades das variveis do processo podem ser especificadas, que
podem ser acessadas atravs do duplo-clique do boto esquerdo do mouse, no
qual uma nova caixa de dilogo se abre, passvel de modificao.
34
Figura 9 - Configurao da corrente de alimentao.
Como j dito, o COCO possui uma interface grfica bastante amigvel. Nesse
sentido, possvel modificar a aparncia das unidades operacionais,
escolhendo entre vrias opes de unidades operacionais j existentes no
software.
35
Figura 10 - Edio da unidade operacional Scilab, aproximando a aparncia da unidade a ser

simulada.
Retornando seo de modificao da unidade operacional, Scilab, pode-se

verificar se a implementao foi bem feita, clicando sobre o cone Test e, se
houver qualquer erro, na seo Output mostrar o erro e onde este est, de
maneira anloga ao Scilab.
Figura 11 - cone de teste do cdigo Scilab implementado no COCO.
36
Se no houver mais nenhum problema, feche as janelas abertas e para resolver

o caso clica-se em Sove, na barra de ferramentas.
Um relatrio da simulao pode ser visualizado, inserindo um Stream report:

definem-se quais portas e variveis so pertinentes para a anlise. Alm deste,
pode-se gerar o relatrio da unidade, o Unit report.
Duas unidades diferentes cada uma com seu cdigo Scilab foram implementadas, uma
utilizando o mesmo mtodo de Newton para convergncia apresentado anteriormente e a
outra usando somente alguns comandos de manipulao da corrente de alimentao (atravs
da linguagem prpria CAPE-OPEN - COCO). Estas recebero as denominaes Scilab
isotrmico I e Scilab isotrmico II, respectivamente. A unidade operacional gerada com os
recursos j definidos pelo COFE ser tratada como Flash isotrmico. Os cdigos foram
exibidos nos Apndices C e D, respectivamente. Porm, no caso da unidade Scilab
isotrmico I o cdigo implementado no forneceu a convergncia esperada, apresentando
uma mensagem de erro, assim como ocorreu para o caso apresentado anteriormente de
utilizao do COCO no Scilab. A Figura 12 mostra a tela principal do COCO onde foi
simulada a unidade Scilab isotrmico I com sua respectiva mensagem de erro.
Figura 12 - Tela de simulao do COCO mostrando o caso do Flash isotrmico I- Mtodo de

Newton, para o qual no obteve sucesso.
37
A Tabela 5 mostra os resultados obtidos na simulao da unidade Scilab Flash

isotrmico II. O pacote termodinmico escolhido foi implementado utilizando a equao de
estado de Peng Robinson e verfica-se que estes apresentaram um desvio relativo menor em
comparao aos valores citados na literatura e tambm bem diferentes daqueles simulados no
Scilab, onde se considerou uma mistura ideal e seguindo a Lei de Raoult.
Tabela 5 - Relatrio da simulao da unidade Scilab Flash isotrmico II, resgatando os
dados gerados pelo prprio COFE.
Corrente
Alimentao
Lquido
Vapor
Unidade
Presso
600
600
600
Psia
Temperatura
60
60
60
Vazo
144
74,5248
69,4752
lbmol/h
Frao molar metano
0,4861
0,191051
0,802593
Frao molar etano
0,1389
0,147447
0,129732
Frao molar propano
0,0694
0,103713
0,032593
0,0625
0,106841
0,0149367
Frao molar n-butano
0,0556
0,0976466
0,0104974
0,0486
0,0895533
0,00467011
0,0417
0,0774261
0,00337729
Frao molar n-hexano
0,0486
0,0927611
0,0012292
0,0278
0,05345
0,000285674
Frao molar n-octano
0,0208
0,0401113
0,0000851717
Entalpia
-693202
-642962
-50239,6
Btu/h
Os mesmos resultados foram obtidos quando se utilizou a unidade operacional Flash

contida na biblioteca interna de modelos do COCO, como pode ser comprovado observando a
Tabela 6 e a Figura 14. Este resultado era esperado, pois no caso do Scilab Flash isotrmico II
trata-se de uma forma mais entendvel e visvel de obter os dados da simulao que feita
internamente pelo caso Flash isotrmico.
38
Figura 13 - Tela de simulao do COCO mostrando o caso do Flash isotrmico II.
Tabela 6 - Relatrio da simulao da unidade Flash isotrmico, gerado atravs da utilizao

da unidade operacional Flash do COFE.
Corrente
Alimentao
Lquido
Vapor
Unidade
Presso
600
600
600
Psia
Temperatura
60
60
60
Vazo
144
74,5248
69,4752
lbmol/h
Frao molar metano
0,4861
0,191051
0,802593
Frao molar etano
0,1389
0,147447
0,129732
Frao molar propano
0,0694
0,103713
0,032593
0,0625
0,106841
0,0149367
Frao molar n-butano
0,0556
0,0976466
0,0104974
0,0486
0,0895533
0,00467011
0,0417
0,0774261
0,00337729
Frao molar n-hexano
0,0486
0,0927611
0,0012292
0,0278
0,05345
0,000285674
Frao molar n-octano
0,0208
0,0401113
0,0000851717
Entalpia
-693202
-642962
-50239,6
Btu/h
39
Figura 14 - Tela de simulao do COCO mostrando o caso do Flash isotrmico.
3.2. FLASH ADIABTICO MULTICOMPONENTE

Prosseguindo com exemplos de unidades operacionais que trabalham simultaneamente
com Scilab e COCO, ser apresentado nessa seo o caso de um flash adiabtico com
alimentao multicomponente. O exemplo escolhido como referncia foi aquele apresentado
no livro Separation Process Principles, 2.ed., de SEADER & HENLEY (2006), pp. 130-131.
Este diz que uma corrente de alimentao contendo hidrognio, metano, benzeno e tolueno
so alimentados a um vaso flash a 120F e 485 psia. Porm, a corrente lquida em equilbrio
saindo do vaso alimentar uma torre de destilao para remover o restante de metano e
hidrognio na mesma. E a presso no prato de alimentao desta torre de 165 psia. Portanto,
uma vlvula que aplica uma queda de presso instalada na entrada do vaso flash,
diminuindo a presso da corrente alimentada de 485 a 165 psia, realizando uma operao
adiabtica. Na resoluo do mesmo, o livro utiliza a equao de estado de Peng-Robinson
para a estimativa termodinmica dos parmetros. A Tabela 7 mostra as quantidades de cada
corrente deste vaso flash.
40
Tabela 7 - Valores para a simulao de um flash adiabtico apresentado por Seader & Henley
(2006).
Corrente
Alimentao
Lquido
Vapor
Unidade
Presso
485
165
165
Psia
Temperatura
120
112
112
Vazo
487,4
471,06
16,34
kmol/h
Vazo molar hidrognio
1,0
0,3
0,7
kmol/h
Vazo molar metano
27,9
12,7
15,2
kmol/h
Vazo molar benzeno
345,1
344,7
0,4
kmol/h
Vazo molar tolueno
113,4
113,36
0,04
kmol/h
Entalpia
-1089,0
-1451,0
362,0
kJ/h
A partir deste exemplo, foi implementado no COCO duas unidades flash e

compararam-se seus resultados. Da mesma forma realizada na seo anterior, foi criada uma
unidade que utiliza um cdigo Scilab (Flash adiabtico Scilab) e outra a partir do modelo
Inside-Out, intrnseco do software COCO (Flash adiabtico COCO). A Figura 15 mostra a
interface obtida aps a simulao de ambas as unidades operacionais no COCO. O cdigo
Scilab mostrado no Apndice E. O mesmo procedimento realizado na seo anterior, de
criao e modificao das unidades foi realizado nessa, com a diferena que a queda de
presso na unidade Flash adiabtico COCO foi instalada uma vlvula antes da alimentao
no vaso.
41

caso do flash adiabtico.
Verifica-se que os resultados para ambas as unidades operacionais simuladas (que
apresentaram os mesmos valores) apresentam valores prximos daqueles da literatura, mas
no idnticos. Devido a este problema foi testado a influncia do pacote termodinmico nos
resultados, mais precisamente das equaes de estado utilizadas, e concluiu-se que estes
influenciam na resposta final. Entre as equaes de estado disponveis no COCO, a que
melhor obteve resultados mais prximos daqueles apresentados na literatura foi aquela de
Soave Redlich Kwong Modificada. A Tabela 8 mostra os resultados da simulao quando se
utiliza esta equao de estado.
42
Tabela 8 - Valores para a simulao do flash adiabtico no COCO utilizando a equao de

estado de Soave Redlich Kwong Modificada.
Corrente
Alimentao
Lquido
Vapor
Unidade
Presso
485
165
165
Psia
Temperatura
120
112
112
Vazo
487,4
471,696
15,7035
kmol/h
Vazo molar hidrognio
1,0
0,0793177
0,920682
kmol/h
Vazo molar metano
27,9
13,5176
14,3824
kmol/h
Vazo molar benzeno
345,1
344,739
0,360552
kmol/h
Vazo molar tolueno
113,4
113,36
0,0398851
kmol/h
Entalpia
-1,41837.107
-1,44959.107
8888,79
kJ/h
3.3. REATOR CSTR ADIABTICO MULTICOMPONENTE

Neste exemplo de CSTR adiabtico foram implementadas duas unidades operacionais,
uma utilizando a unidade CSTR da biblioteca do COCO e a outra inserindo o cdigo Scilab,
assim como nas demais sees vistas anteriormente. Este problema foi retirado da pgina do
grupo Amsterchem, http://www.amsterchem.com, que possui alguns exemplos de utilizao
da linguagem Scilab na criao de unidades. Utilizou-se deste exemplo como ilustrao do
caso de um CSTR adiabtico, pois no foi encontrado na literatura um caso no to complexo
e que tambm trabalhasse com substncias qumicas contidas no banco de dados do COCO
que, na verso gratuita no possui todas as existentes. O cdigo implementado foi descrito no
Apndice F.
O problema testado o seguinte: deseja-se produzir etano, em fase gasosa, a partir da
reao entre o etileno e hidrognio. A alimentao de 10 kmol/h contendo etileno, hidrognio
e nitrognio, nas propores molares mostradas na Tabela 9, alimentada a um CSTR que
opera adiabaticamente. A reao isotrmica tem uma taxa de reao igual a -rA = 4,16667e-6
*C(Etileno)*C(Hidrognio) [mol/s/m3] e um calor de reao de -136330 J/mol. O volume do
reator 0,09025 m3.
Os resultados obtidos por essa simulao, iguais para ambas as unidades, foram
discretizados tambm na Tabela 9, e a Figura 16 mostra a tela de simulao do COCO com as
duas unidades simuladas.
43
Tabela 9 - Valores da simulao do CSTR adiabtico no COCO.

Corrente
Alimentao
Produto
Unidade
Presso
1466560
1466560
Pa
Temperatura
300
324,583
Vazo
10
9,93534
kmol/h
Frao molar etano
0,00650798
Frao molar etileno
0,4
0,396095
Frao molar hidrognio
0,4
0,396095
Frao molar nitrognio
0,2
0,201302

caso de um CSTR adiabtico.
Existe a possibilidade de modificar as condies de operaes das unidades, na
tentativa de obter uma melhor converso na sada do reator. Porm, modificando a presso e a
temperatura, por exemplo, as mudanas so pequenas ou piora-se a produo de etano.Esta
uma vantagem de utilizao do COCO, onde as simulaes so realizadas com rapidez e de
uma maneira organizada na forma de planilhas.
44
3.4. REATOR PFR ADIABTICO MULTICOMPONENTE

Para este caso foi considerado como parmetro de comparao o exemplo mostrado
por FOGLER (1999), no captulo 8, sobre a isomerizao em fase lquida do n-butano em um
PFR a alta presso, usando essencialmente traos de um catalisador lquido. Este catalisador
d reao uma taxa especfica de 31,1 h-1 a 360 K. Deseja-se calcular o volume necessrio
para processar 163 kmol/h de uma mistura de 90 mol % de n-butano e 10 mol% de isopentano
(considerado inerte). O HRx igual a -6900J/molbutano e a energia de ativao 65,7 kJ/mol. A
constante da taxa de reao segue a Equao de Arrhenius e a concentrao inicial de nbutano, CA0, 9,3 kmol/m3.
Neste exemplo, foi implementado a unidade operacional Scilab e seu cdigo foi
mostrado no Apndice G. A tela contendo a simulao realizada aquela mostrada na Figura
17.

caso do flash adiabtico.
A planilha contendo os resultados da simulao mostrada na Tabela 10. Verifica-se
que os valores obtidos so muito prximos daqueles exibidos na literatura. Porm, o cdigo
ainda necessita de melhorias, como por exemplo, no momento de integrar a equao do
balano de massa para obter o volume do reator. No cdigo implementado foi utilizado o
45
mtodo numrica de discretizao da integral, o que gera um pequeno desvio do valor real. O
valor do volume desse reator PFR citado na literatura para chegar converso de equilbrio
igual a 3,4 m3, enquanto que na simulao foi encontrado um valor de 2,58 m3.
Tabela 10 - Valores da simulao do PFR adiabtico
Corrente
Alimentao
Produto
Unidade
Presso
600
600
Psi
Temperatura
330
359,469
Vazo
163
163
kmol/h
Frao molar n-butano
0,9
0,27
0,1
0,1
0,63
3.5. CASO ESPECIAL DE UM SEPARADOR MEMBRANA DE BIOGS

Um caso especial, retirado tambm do site do Amsterchem aquele sobre a separao
de compostos contidos num biogs atravs de um sistema de membranas. Esta ilustrao
feita apenas para mostrar que possvel implementar uma unidade operacional que utiliza o
cdigo Scilab interligada a outros equipamentos e montar um fluxograma de um processo
industrial. A Figura 18 mostra o fluxograma e as planilhas de resultados aps simulao.
Outra novidade neste caso a utilizao de um arquivo adicional, o gpcross.sci que calcula a
transferncia de massa dos compostos entre as fibras e a calandra. Este cdigo, por ser escrito
na planilha de simulao do Scilab, pode ser simulado separadamente neste ambiente ou
mesmo servir como arquivo-base para outros problemas.
Os dois cdigos, o principal simulado no COCO e o gpcross.sci foram descritos nos
Apndices H.1 e H.2.
46
Figura 18 Fluxograma e planilhas de resultados obtidos para o caso de Separao de Biogs.

3.6. OBTENO DE PARMETROS FSICOS E TERMODINMICOS
A outra vertente desse trabalho foi a criao de um cdigo Scilab (e que trabalha neste
ambiente) que l e busca dados termodinmicos e fsicos das substncias requeridas e so
apresentadas na tela CONSOLE do mesmo software estes dados reqsueridos. Este cdigo,
descrito no Apndice I, procura dados de duas listas diferentes de compostos simultaneamente
gerados a partir da decodificao binria de dois softwares lvres disponveis atualmente, o
FYDAT e o DIPPR. O FYDAT possui uma lista de aproximadamente 700 substncias
qumicas e o DIPPR, 1023 compostos. Este bastante til e rpido na determinao das
propriedades, pois evita a necessidade de recorrer literatura, o que pode custar bastante
tempo.
4 CONCLUSO
Conclui-se com este trabalho que h uma infinidade de possibilidades que podem ser
utilizadas para a simulao de unidades operacionais, servindo como uma base de modelos
dinmicos para a utilizao por pesquisadores ou at mesmo como material de aprendizado
em sala de aula. A comunicao entre softwares muito vantajosa, visto que cada um
contribui com que melhor possui para facilitar e melhorar o entendimento de processos
47
qumicos, como o caso da ligao entre o Scilab e o COCO. Quanto mais ferramentas de
modelagem e simulao so disponveis s pessoas mais tempo pode-se dedicar ao
aperfeioamento dos cdigos, sem se preocupar se a propriedades fsico-qumicas e
termodinmicas dos compostos esto corretas e se as fontes so confiveis.
O objetivo deste trabalho bastante audacioso, e nestas pginas foram descritas
somente o incio de um projeto que ainda tem um potencial enorme de contribuio.
SUGESTES PARA PRXIMOS TRABALHOS

Futuros trabalhos que completem este trabalho sero bem-vindos e sugestes para
estes seriam: a verificao dos cdigos implementados e possveis melhoramentos e solues
para os problemas de no-convergncia encontrados; simulao de novas unidades
operacionais, como reatores biolgicos, colunas de destilao, absoro e flash trifsico. Alm
disso, criar fluxogramas de processos consagrados e verificar suas eficincias outra forma
bastante interessante para se realizar.
REFERNCIAS BIBLIOGRFICAS
BIEGLER, L.T. et al. Sistematic Methods of Chemical Process Design. Rio de

Janeiro, Prentice-Hall do Brasil Ltda, 1999.
BOSTON, J. F. & H. I. BRITT, A Radically Different Formulation and Solution of the
Single-Stage Flash Problem. Comp. Chem.Eng., 2, 109, 1978.
FOGLER, H.S. Elementos de Engenharia das Reaes Qumicas 3.ed. Rio de Janeiro,
LTC Editora, 1999.
PENG, D. Y. & ROBINSON, D. B. A new two-constant equation of state. Ind. Eng.
Chem. Fundam.. v. 15, p. 59-64, 1976.
REID, R.C.; PRAUSNITZ, J.M. The properties of gases and liquids. 5.ed. New York
McGraw-Hill, 2001.
SEADER, J.D.; HENLEY, E.J. Separation process principles. 2.ed. EUA, John Wiley
& Sons, 2006.
SEIDER, W.D. et al. Product & Process Design Principles. 2.ed. EUA, John Wiley &
Sons, 2003.
Soave, G., Equilibrium Constants for a Modified Redlich-Kwong Equation of State.
Chem. Eng. Sci., 27, 1186, 1972.
48
APNDICES
49
APNDICE A - CDIGO SCILAB UTILIZADO PARA SIMULAR OS CASOS DE

FLASH ISOTRMICO, PONTO DE BOLHA E PONTO DE ORVALHO.
50
51
52
53
54
55
APNDICE B CDIGO IMPLEMENTADO NO SCILAB QUE UTILIZA AS

FERRAMENTAS DO COCO.
56
57
APNDICE C CDIGO SCILAB DE CLCULO DO FLASH ISOTRMICO

UTILIZANDO O MTODO DE NEWTON E A PLATAFORMA DO COCO.
//Flash isotrmico utilizando o Mtodo de Newton
//Caso para uma separao entre metano, etano, propano, n-butano, isobutano,
//n-pentano, i-pentano, n-hexano, n-heptano, n-octano.
CAS_Methane="74-82-8";
CAS_Ethane="74-84-0";
CAS_Propane="74-98-6";
CAS_Butane="106-97-8";
CAS_Isobutane="75-28-5";
CAS_Pentane="109-66-0";
CAS_Isopentane="78-78-4";
CAS_Hexane="110-54-3";
CAS_Heptane="142-82-5";
CAS_Octane="111-65-9";
//Obtendo os indices dos componentes
cas=getCompoundConstant("CASRegistryNumber");
indexMethane=0;
indexEthane=0;
indexPropane=0;
indexButane=0;
indexIsobutane=0;
indexPentane=0;
indexIsopentane=0;
indexHexane=0;
indexHeptane=0;
indexOctane=0;
for i=1:nComp
casNumber=cas(i,:);
if (indexMethane==0) then if (strcmp(casNumber,CAS_Methane)==0) then indexMethane=i; end; end;
if (indexEthane==0) then if (strcmp(casNumber,CAS_Ethane)==0) then indexEthane=i; end; end
if (indexPropane==0) then if (strcmp(casNumber,CAS_Propane)==0) then indexPropane=i; end; end
if (indexButane==0) then if (strcmp(casNumber,CAS_Butane)==0) then indexButane=i; end; end
if (indexIsobutane==0) then if (strcmp(casNumber,CAS_Isobutane)==0) then indexIsobutane=i; end; end
if (indexPentane==0) then if (strcmp(casNumber,CAS_Pentane)==0) then indexPentane=i; end; end
if (indexIsopentane==0) then if (strcmp(casNumber,CAS_Isopentane)==0) then indexIsopentane=i; end; end
58
if (indexHexane==0) then if (strcmp(casNumber,CAS_Hexane)==0) then indexHexane=i; end; end

if (indexHeptane==0) then if (strcmp(casNumber,CAS_Heptane)==0) then indexHeptane=i; end; end
if (indexOctane==0) then if (strcmp(casNumber,CAS_Octane)==0) then indexOctane=i; end; end
end
if (indexMethane==0) then error("Componente Metano nao encontrado"); end
if (indexEthane==0) then error("Componente Etano nao encontrado"); end
if (indexPropane==0) then error("Componente Propano nao encontrado"); end
if (indexButane==0) then error("Componente Butano nao encontrado"); end
if (indexIsobutane==0) then error("Componente Isobutano nao encontrado"); end
if (indexPentane==0) then error("Componente Pentano nao encontrado"); end
if (indexIsopentane==0) then error("Componente Isopentano nao encontrado"); end
if (indexHexane==0) then error("Componente Hexano nao encontrado"); end
if (indexHeptane==0) then error("Componente Heptano nao encontrado"); end
if (indexOctane==0) then error("Componente Octano nao encontrado"); end
// Busca dos dados de entrada

T=getFeedProp(1,"temperature")
P=getFeedProp(1,"pressure")
imax=getParameter("max_int") // numero maximo de iteracoes para Metodo de Newton
tol=getParameter("tolerance") //precisao desejada para o Metodo de Newton
zf=getFeedProp(1,"fraction")
F=getFeedProp(1,"flow")
totalFlow=sum(F);
VF=getParameter("vf")
// Calculo dos valores de K
Ac_coef=getSinglePhaseProperty("activityCoefficient","liquid",T,P,[zf])
Fug_coef=getSinglePhaseProperty("FugacityCoefficient","vapor",T,P,[zf])
Fug=getSinglePhaseProperty("Fugacity","vapor",T,P,[zf])
for r=1:nComp
k(r)=(Ac_coef(r)*Fug(r))/(P*Fug_coef(r));
end
disp('Valores de k:');
disp (k);
//Metodo de Newton
iter=0;
toler=1e10;
//iteracao alfa funcao derivada
while (toler>tol) & (iter<imax) then
for d=1:nComp
x(d)=zf(d)/(1+VF*(k(d)-1));
y(d)=k(d)*x(d);
end
disp('Xs sao'); disp(x);
disp('Ys sao'); disp(y);
sum_x=sum(x);disp('somatorio de x eh:'); disp(sum_x);
sum_y=sum(y);disp('somatorio de y eh:');disp(sum_y);
Ac_coef=getSinglePhaseProperty("activityCoefficient","liquid",T,P,[x]);
Fug_coef=getSinglePhaseProperty("FugacityCoefficient","vapor",T,P,[y])
59
for s=1:nComp
k(s)=(Ac_coef(s)*Fug(s))/(P*Fug_coef(s));
end
disp('valores de k eh'); disp(k);
if abs(sum_x-sum_y)>1e-5 then
f=0;
for j=1:nComp
f=f+zf(j)*(1-k(j))/(1+VF*(k(j)-1));
end
disp("f eh"); disp(f);
df=0;
for g=1:nComp
df=df+zf(g)*(1-k(g))^2/(1+VF*(k(g)-1))^2;
end
disp("df eh"); disp(df); VF_old=VF;
if (VF>0) & (VF<1) then
if abs(df)>1e-8 then
VF= VF - f/df;
disp("vf eh"); disp(VF);
else
error ('Derivada nula');
end
else
error('Sistema nao fornece duas fases em equilibrio nessas condicoes');
end
else
break;
end
disp('tolerancia eh');
toler=abs(VF-VF_old);
disp (toler);
disp('iteracao:'); iter=iter+1; disp(iter);
end
disp('Composicao da fase liquida')
disp(x);
disp('Composicao da fase vapor')
disp(y);
V=totalFlow*VF;
L=totalFlow-V;
disp("A vazao de vapor eh:");
disp(V);
disp("A vazao de liquido eh:");
disp(L);
disp('setproduct y');
setProduct(1,V,y,"pressure",P,"temperature",T);
disp('setproduct x');
setProduct(2,L,x,"pressure",P,"temperature",T);
60
APNDICE D CDIGO SIMPLIFICADO PARA O CLCULO DO FLASH

ISOTRMICO QUE UTILIZA SOMENTE COMANDOS DO COCO.
T=getFeedProp(1,"temperature")
P=getFeedProp(1,"pressure")
F=getFeedProp(1,"flow")
totalFlow=sum(F)
[phases,phasefractions,compositions,T,P]=getEquilibrium(zf,"pressure",P,"temperature",T)
vaporFraction=phasefractions(1);
vaporPhase=0;
phaseCount=size(phasefractions,1);
for i=1:phaseCount
if (getAggregationState(phases(i,:))==AGGSTATE_VAPOR) then
vaporPhase=i;
break
end
end
liquidPhase=0;
for j=1:phaseCount
if (getAggregationState(phases(j,:))==AGGSTATE_LIQUID) then
liquidPhase=j;
break
end
end
fv=totalFlow*vaporFraction;
fl=totalFlow-fv; disp(fl);
setProduct(1,fv,compositions(vaporPhase,:),"pressure",P,"temperature",T);
setProduct(2,fl,compositions(liquidPhase,:),"pressure",P,"temperature",T);
APNDICE E CDIGO SCILAB UTILIZADO PARA IMPLEMENTAR O CASO

DE FLASH ADIABTICO NO COCO.
// Exemplo de Flash Adiabatico
F=getFeedProp(1,"totalFlow")
delta_p=getParameter("Pressure drop")
delta_T=getParameter("Temp drop")
delta_h=getParameter("Heat duty")
if ((delta_p==0)&(delta_h==0)) then
//usando a alimentao de equilibrio
[phases,phasefractions,compositions,T,P]=getFeedEquilibrium(1)
else
//Calculando o equilibrio depois da queda de presso e trabalho
if (delta_h~=0) then if (F==0) then delta_h=0, end; end
if (delta_h~=0) then delta_h=delta_h/F; end
p=getFeedProp(1,"pressure")-delta_p
t1=getFeedProp(1,"temperature")-delta_T
h=getFeedProp(1,"enthalpy")+delta_h
[phases,phasefractions,compositions,T,P]=getEquilibrium(zf,"pressure",p,"temperature",t1)
61
end
//Encontrando a fase vapor
vaporPhase=0;
for i=1:phaseCount
if (getAggregationState(phases(i,:))==AGGSTATE_VAPOR) then
vaporPhase=i;
break
end
end
//Checando se temos uma fase vapor
if (vaporPhase>0) then
vaporFraction=phasefractions(vaporPhase)
if (phaseCount==1) then
//Tem-se somente a fase vapor
setProduct(1,F,zf,"pressure",P,"temperature",T)
setProduct(2,0,zf,"pressure",P,"temperature",T)
else
//Fase vapor
setProduct(1,F*vaporFraction,compositions(vaporPhase,:),"pressure",P,"temperature",T)
//Todas as outras fases por diferena
f=zf*F-(F*vaporFraction*compositions(vaporPhase,:));
ftot=sum(f);
if (ftot>0) then x=f/ftot; end;
setProduct(2,ftot,x,"pressure",P,"temperature",T);
end
else
disp(No tem fase vapor);
setProduct(1,0,zf,"pressure",P,"temperature",T)
setProduct(2,F,zf,"pressure",P,"temperature",T)
vaporFraction=0
end
APNDICE F CDIGO SCILAB UTILIZADO PARA A DEFINIO

DA UNIDADE DE CSTR ADIABTICO.
disp('Este eh um exemplo de um CSTR adiabatico')
disp('Reao: Etileno + Hidrogenio -> Etano')
//
//O sistema termico apresentado abaixo deve incluir os trs componentes, e
// expor CASRegistryNumber, enthalpyF e volume
//Pressumir que a reao acorre na fase vapor
CAS_Ethane="74-84-0";
CAS_Ethylene="74-85-1";
CAS_Hydrogen="1333-74-0";
//Estequiometria
stoi_Ethane=1;
stoi_Ethylene=-1;
stoi_Hydrogen=-1;
62
compCount=size(cas,1)
indexEthane=0;
indexEthylene=0;
indexHydrogen=0;
for i=1:compCount
casNumber=cas(i,:);
if (indexEthane==0) then if (strcmp(casNumber,CAS_Ethane)==0) then indexEthane=i; end; end
if (indexEthylene==0) then if (strcmp(casNumber,CAS_Ethylene)==0) then indexEthylene=i; end; end
if (indexHydrogen==0) then if (strcmp(casNumber,CAS_Hydrogen)==0) then indexHydrogen=i; end; end
end
if (indexEthane==0) then error("Compound Ethane not found"); end
if (indexEthylene==0) then error("Compound Ethylene not found"); end
if (indexHydrogen==0) then error("Compound Hydrogen not found"); end
//A taxa da reao (mol/m3/s) dada por
//r*C(Ethylene)*C(Hydrogen)
r=4.16667e-006;
//Nome da fase vapor
vaporPhaseName=[];
for i=1:nPhase
if (phaseAggregationStates(i)==AGGSTATE_VAPOR) then
vaporPhaseName=phaseNames(i,:)
break;
end
end
if (isempty(vaporPhaseName)) then error("Nenhuma definicao para a fase vapor"); end
//Checando que alimentao contm somente vapor
[phases]=getFeedEquilibrium(1);
if (size(phases,1)~=1) then error("Alimentao deve ser Feed must be in vapor state"); end
if (strcmpi(phases(1,:),vaporPhaseName)~=0) then error("Feed must be in vapor state"); end
//heat of reaction, J/mol
// Hreac=-136330;
//using heat of reaction assumes that formation terms are not included in
// default enthalpy; instead we will use enthalpy, for which formation
// terms are included and we will not use heat of reaction
//get reactor volume, m3
vol=getParameter("Reactor volume")
//get pressure drop, Pa
dP=getParameter("Pressure drop")
//get feed pressure
p=getFeedProp(1,"pressure");
//reactor pressure
p=p-dP
//get total and component flow in
flow=getFeedProp(1,"flow")
totalFlow=sum(flow)
if (totalFlow==0) then error("Reactor calculations cannot be performed for zero flow"); end
//calculate total reactor H
hfeed=totalFlow*getFeedProp(1,"enthalpyF");
63
//we solve for reactor extent, mol/s

function d=ErrorFunc(extent)
//at the current extent, calculate the compound flows
f=flow;
f(indexEthane)=f(indexEthane)+extent*stoi_Ethane;
if (f(indexEthane)<0) then f(indexEthane)=0; end
f(indexEthylene)=f(indexEthylene)+extent*stoi_Ethylene;
if (f(indexEthylene)<0) then f(indexEthylene)=0; end
f(indexHydrogen)=f(indexHydrogen)+extent*stoi_Hydrogen;
if (f(indexHydrogen)<0) then f(indexHydrogen)=0; end
//calc X
ft=sum(f);
x=f/ft;
//calculate reactor temperature from enthalpy flash
[phases,phasefractions,compositions,T,P]=getEquilibrium(x,"enthalpyF",hfeed/ft,"pressure",p);
//we ignore the fact that non-vapor phases can be returned, we will check at the solution
//calc vapor volume at T and p
vapVol=getSinglePhaseProperty("volume",vaporPhaseName,T,p,x);
//total concentration is 1 / volume
Ctotal=1/vapVol;
//r*C(Ethylene)*C(Hydrogen)
CEthylene=x(indexEthylene)*Ctotal;
CHydrogen=x(indexHydrogen)*Ctotal;
rate=vol*r*CEthylene*CHydrogen; //mol / m3 / s * m3 = mol / s
//return error for current extent guess
d=extent-rate;
endfunction
//solve for the reaction extent
extent=fsolve(0,ErrorFunc)
//calculate the resulting compound flows
f=flow;
f(indexEthane)=f(indexEthane)+extent*stoi_Ethane;
f(indexEthylene)=f(indexEthylene)+extent*stoi_Ethylene;
f(indexHydrogen)=f(indexHydrogen)+extent*stoi_Hydrogen;
//calc X
ft=sum(f);
x=f/ft;
[phases,phasefractions,compositions,T,P]=getEquilibrium(x,"enthalpyF",hfeed/ft,"pressure",p);
//check that only vapor results
if (size(phases,1)~=1) then error("Reaction product includes multiple phases"); end
if (strcmpi(phases(1,:),vaporPhaseName)~=0) then error("Reaction product is not vapor"); end
//set the product (no need to perform enthalpy flash, we have done so above; use the resulting T, more efficient)
setProduct(1,ft,x,"pressure",p,"temperature",T)
APNDICE G CDIGO SCILAB UTILIZADO PARA A DEFINIO

DA UNIDADE PFR ADIABTICO.
disp('Isomerizacao do n-Butano em fase liquida');
disp('.........n-C4H10<-->i-C4H10..............');
disp('A reacao eh conduzida abiabaticamente em fase liquida sob alta pressao,');
disp('usando essencialmente tracos de um catalisador liquido (iso-pentano)');
CAS_Butane="106-97-8";
CAS_Isopentane="78-78-4";
64
CAS_i_Butane="75-28-5";
//Estequiometria
stoi_i_Butane=1;
stoi_Butane=-1;
stoi_Isopentane=0;
compCount=size(cas,1)
indexButane=0;
indexIsopentane=0;
indexi_Butane=0;
for i=1:compCount
casNumber=cas(i,:);
if (indexButane==0) then if (strcmp(casNumber,CAS_Butane)==0) then indexButane=i; end; end
if (indexIsopentane==0) then if (strcmp(casNumber,CAS_Isopentane)==0) then indexIsopentane=i; end; end
if (indexi_Butane==0) then if (strcmp(casNumber,CAS_i_Butane)==0) then indexi_Butane=i; end; end
end
if (indexButane==0) then error("Componente Butano nao encontrado"); end
if (indexIsopentane==0) then error("Componente Isopentano nao encontrado"); end
if (indexi_Butane==0) then error("Componente i_Butano nao encontrado"); end
//Caracterizacao da fase liquida
liquidPhaseName=[];
for i=1:nPhase
if (phaseAggregationStates(i)==AGGSTATE_LIQUID) then
liquidPhaseName=phaseNames(i,:)
break;
end
end
if (isempty(liquidPhaseName)) then error("Nao ha formacao de fase liquida"); end
//Checando que alimentao contm somente liquida
[phases]=getFeedEquilibrium(1);
if (size(phases,1)~=1) then error("Alimentao deve estar no estado liquido"); end
if (strcmpi(phases(1,:),liquidPhaseName)~=0) then error("Alimentao deve estar no estado liquido"); end
//Obtendo a queda de pressao, Pa
dP=getParameter("Pressure drop")
//Obtendo a pressao e temperatura de alimentacao
p=getFeedProp(1,"pressure");
To=getFeedProp(1,"temperature");
//pressao do reator
p=p-dP;
//Obtendo o fluxo de entrada por componente e total
flow=getFeedProp(1,"flow");disp(flow)
totalFlow=sum(flow);
if (totalFlow==0) then error("Os calculos para o reator nao podem ser desenvolvidos para fluxo zero"); end
zf=getFeedProp(1,"fraction");
//Obtendo os parametros da reacao
dHrx1=getParameter("dHrx");//J/mol
R=getParameter("Constgas");
E=getParameter("Eativ");
Kc1=getParameter("Const_C");
65
k1=getParameter("Const_reac");
v=getParameter("vazao_vol");disp(v);
Ca0=totalFlow/v*zf(indexButane); disp('ca0 eh'); disp(Ca0);
Cp(indexButane)=getParameter("CpnB");
Cp(indexi_Butane)=getParameter("CpiB");
Cp(indexIsopentane)=getParameter("CpiP");
theta(indexButane)=flow(indexButane)/flow(indexButane);
theta(indexi_Butane)=flow(indexi_Butane)/flow(indexButane);
theta(indexIsopentane)=flow(indexIsopentane)/flow(indexButane);
sum_thetaCp=
Cp(indexButane)*theta(indexButane)+Cp(indexi_Butane)*theta(indexi_Butane)+Cp(indexIsopentane)*theta(ind
exIsopentane);
disp(sum_thetaCp);disp(flow(indexButane));
deltaCp=Cp(indexButane)-Cp(indexi_Butane);
x=0:0.05:1;
for z=1:21
T(z)=To-(dHrx1*x(z))/sum_thetaCp;
k(z)=k1*exp(E/R*(1/360-1/T(z)));
Kc(z)=Kc1*exp(dHrx1/R*(1/333-1/T(z)));
Xe(z)=Kc(z)/(1+Kc(z));
rA(z)=k(z)*(Ca0)*(1-(1+1/Kc(z))*x(z));
t(z)=flow(indexButane)/rA(z);
end
disp (T);
disp(k);
disp(Kc);
disp(Xe);
disp(rA);disp(t);
V=0;
h1=(x(15)-x(1))/14;
V=V+3/8*h1*(t(1)+3*t(2)+3*t(3)+2*t(4)+3*t(5)+3*t(6)+2*t(7)+3*t(8)+3*t(9)+2*t(10)+3*t(11)+3*t(12)+2*t(1
3)+3*t(14)+t(15));
disp('Volume do PFR eh:');
disp(V);
//calculate total reactor H
hfeed=totalFlow*getFeedProp(1,"enthalpyF");
Ctotal=Ca0/zf(indexButane);
Ci_Butane=zf(indexi_Butane)*Ctotal;disp(Ci_Butane);
CButane=Ca0;
CIsopentane=zf(indexIsopentane)*Ctotal;
//calculate the resulting compound flows
disp(x(15));
flow(indexi_Butane)=flow(indexi_Butane)+x(15)*stoi_i_Butane*flow(indexButane);disp(flow(indexi_Butane));
flow(indexButane)=flow(indexButane)+x(15)*stoi_Butane*flow(indexButane);
flow(indexIsopentane)=flow(indexIsopentane)+x(15)*stoi_Isopentane*flow(indexButane);
//calc X
disp('flow eh:'); disp(flow);
ft=sum(flow);disp('ft eh:');disp(ft)
b=flow/ft; disp('b eh:');disp(b);
[phases,phasefractions,compositions,T,P]=getEquilibrium(b,"enthalpyF",hfeed/ft,"pressure",p);
66
disp(phases);
//set the product (no need to perform enthalpy flash, we have done so above; use the resulting T, more efficient)
setProduct(1,ft,b,"pressure",p,"temperature",T)
APNDICE H.1 CDIGO PRINCIPAL SCILAB DA UNIDADE DE

SEPARAO DE BIOGS.
/ Membrane properties library.
// Mem.Comp - chemical formula
// Mem.CAS - CAS number for each component (used to identify components in the feed stream)
// Mem.alfa membrane selectivity relative to a reference component ref
//
alfa(i) = Q(i)/Q(ref) where Q(i) represents membrane permeability [mol/(s*Pa*m^2)]
// Mem.Qref membrane permeability of a reference component [mol/(s*Pa*m^2)]
Mem.Qref=1*(3.348E-9);
// Permeance [mol/(s*Pa*m2)] of reference component
// 1 Barrer = 3.348E-9 [mol/(s*Pa*m2)] with a membrane thikness of 1000.
Mem.Comp=[];
Mem.CAS=[];
Mem.alfa=[];
Mem.Comp($+1)="CO2";Mem.CAS($+1)="124-38-9";
Mem.alfa($+1)=15;
Mem.Comp($+1)="H2O";
Mem.CAS($+1)="7732-18-5";
Mem.alfa($+1)=1000;
Mem.Comp($+1)="CH4";
Mem.CAS($+1)="74-82-8";
Mem.alfa($+1)=1;
Mem.Comp($+1)="C2H6";
Mem.CAS($+1)="74-84-0";
Mem.alfa($+1)=1.2;
Mem.Comp($+1)="H2S";
Mem.CAS($+1)="7783-06-4";
Mem.alfa($+1)=750;
Mem.Comp($+1)="N2"; Mem.CAS($+1)="7727-37-9";
Mem.alfa($+1)=1;
Mem.Comp($+1)="O2"; Mem.CAS($+1)="7782-44-7";
Mem.alfa($+1)=5;
// Verify that membrane properties are defined for all components.
alfa=zeros(nComp,1); Comp=[];
CompNames=getCompoundConstant("chemicalFormula");
CASRegNum=getCompoundConstant("casRegistryNumber");
for i=1:nComp,
ind=find(Mem.CAS==CASRegNum(i));
if isempty(ind) then
error('Missing permeability data for '+CompNames(i));
else
alfa(i)=Mem.alfa(ind); Comp(i)=Mem.Comp(ind);
end
end
// Feed stream properties.
nf=getFeedProp(1,"totalFlow"); // mol/s
xf=getFeedProp(1,"fraction"); // Tf=getFeedProp(1,"temperature"); // K
Pf=getFeedProp(1,"pressure"); // Pa
// Get input parameters
g=getParameter("Trans membrane pressure ratio"); // feed to permeate pressure ratio [-]
A=getParameter("Membrane area"); // membrane area [m2]
R=Pf*A*Mem.Qref/nf; // dimensionless membrane area (permeation factor)
// Simulate membrane module
getf('gpcross.sci'); // load the simulation model (cross flow model)
67
[xr,xp,th,rp,rr]=gpcross(xf,alfa,R,g); // run simulation

nr=(1-th)*nf;// retentate
np=th*nf;// permeate
Pp=Pf/g; // permeate pressure
Pr=Pf; // retentate pressure
// Report product streams
setProduct(1,nr,xr,"temperature",Tf,"pressure",Pr);
setProduct(2,np,xp,"temperature",Tf,"pressure",Pp);
// Set output parameters
setParameter("Dimensionless membrane area",R);
// Generate a report
txt=[
'|===================================================================|'
'| Membrane module for gas separation
|'
'|
Feed side
|'
'| Feed
|--------------------------| Retentate
|'
'| nf,z --->|
plug flow
|-----> xr, nr
|'
'|
|==========================|
|'
'| np,y <---|
crossflow
|
|'
'| Permeate |--------------------------|
|'
'|
Permeate side
|'
'|===================================================================|'
'| Comp. |Q(i)/Q(c)| Feed | Retentate | Permeate |'
'|
| alfa(i) | xf | xr | rr | xp | rp |'
'|===================================================================|'];
for j=1:nComp,
txt=[txt;sprintf('|%7s| %7.1f | %0.3f | %0.3f | %0.3f | %0.3f | %0.3f
|\n',Comp(j),alfa(j),xf(j),xr(j),rr(j),xp(j),rp(j))];
end
txt=[txt;
sprintf('|===================================================================|'
);
sprintf('|Module cut rate
th = np/nf
= %6.3f
sprintf('|Pressure ratio
g = Pf/Pp
= %6.1f
sprintf('|Permeation factor R = Q(c)*Pf*A/nf = %6.3f
|\n',th);
|\n',g);
|\n',R);
sprintf('|===================================================================|\
n')];
setReport("Membrane module",txt);
68
APNDICE
H.2
ARQUIVO
ADICIONAL
SCILAB
PARA
IMPLEMENTAO DA UNIDADE DE SEPARAO DE BIOGS.
69
APNDICE I ARQUIVO SCILAB DE LEITURA DE PROPRIEDADES

DAS SUBSTNCIAS CONTIDAS NO DYPPR E FYDAT.
70
71
72

Cálculo de Flash2

Enviado por

Direitos autorais:

Formatos disponíveis

Cálculo de Flash2

Enviado por

Dados do documento

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Cálculo de Flash2

Enviado por

Direitos autorais:

Formatos disponíveis

1

UNIVERSIDADE FEDERAL DE UBERLNDIA

Desenvolvimento de uma Biblioteca Eletrnica de

Fernando Rosa Gomes

UNIVERSIDADE FEDERAL DE UBERLNDIA

Desenvolvimento de uma Biblioteca Eletrnica de Modelos

Fernando Rosa Gomes

Monografia de graduao apresentada

Graduao do curso de Engenharia Qumica.

BANCA EXAMINADORA DA MONOGRAFIA DA DISCIPLINA PROJETO DE

O mestre na arte da vida faz pouca distino entre o

IMPORTNCIA DE FONTES DE DADOS TERMODINMICOS ............................ 3

DIPPR (Design Institute for Physical Properties)........................................................ 3

A REGRA DE FASES DE GIBBS E GRAU DE LIBERDADE ................................... 7

CLCULOS DE FLASH SISTEMA DE NICO ESTGIO DE EQUILBRIO....... 9

Derivaes do Modelo Flash ..................................................................................... 10

Estratgias para os clculos de flash ...................................................................... 13

Sequncia de clculo Flash TP ............................................................................... 14

Mtodo Inside-Out para os Clculos de Flash ........................................................... 15

Balano de energia e a Primeira Lei da Termodinmica ........................................... 19

CSTR (Continuous-Stirred-Tank-Reactor) ................................................................ 20

PFR (Plug-Flow Tubular Reactor) ............................................................................. 22

RESULTADOS E DISCUSSO ...................................................................................... 24

FHASH ISOTRMICO MULTICOMPONENTE........................................................ 24

FLASH ADIABTICO MULTICOMPONENTE........................................................ 39

REATOR CSTR ADIABTICO MULTICOMPONENTE ......................................... 42

REATOR PFR ADIABTICO MULTICOMPONENTE ............................................ 44

CASO ESPECIAL DE UM SEPARADOR MEMBRANA DE BIOGS................ 45

OBTENO DE PARMETROS FSICOS E TERMODINMICOS ...................... 46

SUGESTES PARA PRXIMOS TRABALHOS ................................................................. 47

APNDICE C CDIGO SCILAB DE CLCULO DO FLASH ISOTRMICO

para pr-aquecer a alimentao e para resfriar o produto, um comportamento bastante

qumicos, reatores biotecnolgicos etc. e consiste na investigao e desenvolvimento do

2.1. IMPORTNCIA DE FONTES DE DADOS TERMODINMICOS

2.1.1. DIPPR (Design Institute for Physical Properties)

constantes. Dados crus da literatura e coeficientes de equaes para 15 propriedades

Logaritmo do coeficiente de fugacidade, entre outras.

Os padres CAPE-OPEN so normas uniformes para interface de componentes de

COUSCOUS o CAPE-OPEN Unit Operations Simple contm divisores,

Figura 1 - Janela principal do ambiente de simulao do COCO-COFE.

2.3. A REGRA DE FASES DE GIBBS E GRAU DE LIBERDADE

A descrio de um sistema de estgio nico no equilbrio fsico envolve variveis

onde C o nmero de componentes e P o nmero de fases no equilbrio. A equao (1)

onde o 2 refere-se a presso e temperatura do equilbrio, enquanto o termo CP o nmero

no qual o primeiro termo, P, refere necessidade do somatrio das fraes molares ou

fraao molar de i na fase (1)

F = V E = (CP + 2 ) - [P + C(P-1)] = C-P + 2

2.4. CLCULOS DE FLASH SISTEMA DE NICO ESTGIO DE EQUILBRIO

Em programas de simulao de processos, a determinao de um sistema de separao

2.4.1. Derivaes do Modelo Flash

Para ilustrar a discusso dessa seo, considere a operao de separao de fase

i (l,T)f i 0 (T, P)x i = (v,T)Py i

E o balano de entalpia porduz

Definindo ainda o valor K:

K i = i (x,T)f i (T, P)xi / (y,T)P i=1,...NC

Portanto tem-se (3 NC+5) variveis (yi, xi, L, T, P, Q, V) e somente (3 NC+1)

Neste momento, necessita-se de uma especificao adicional em qualquer um dos

Claramente, qualquer uma das equaes conduz pseudo-solues (na composio de

Estratgias para os clculos de flash

O modelo flash pode ser dado concisamente por: