TechnicralArticle

ogo
Q

Desenvolvimento de uma

metodolo g is p era analise e

determine ao dos coeficientes

globais de transferencia de
calor de uma unidade de

evapora de
ao licor negro
determine and
Development of a methodology to analyze global coeficients of heat transfer of a black liquor evaporatingplant

Autores Alberto Avelar Costa Marcelo Cardoso Song Won Park Hamilton Saulo Moreira Brandao3
George

chave
Palavras ao eva
simula
aocoeficiente global
pora de transfe Abstract
rencia de calor In order to evaluate the influence of the black liquor viscosity in the
evaporation recovery units of a Kraft process more speclfzcally in the heat
it data and process parame
RESUMO transference phenomenon was
accomplished a

ters collecting of the Unit

ofEvaporation II of Cenibra During a period of
de analisar ainflu
Com o objetivo four months between May and September of 2003 information such as
encia da viscosidade do licor negro em vapor consumption liquor feeding rate concentrated liquor production

unidade de evapora temperatures andpressures at each effect among others was collected hour
uma aomais es

ly Based on these data an analysis of the performance of operation was

pecificamente no fenomeno de trans elaborated with assistance of the commercial simulator nGEMS 5 3 The
ferencia de calor foi efetuado um le parameters not evaluated directly such as the vapor consumption the con
centration of solids in the liquor and the global heat transfer coefficient
vantamento de dados e parametros de

naUnidade de U for each body of evaporation have been determined from process simu
processo aoII
Evapora lation and analyzea nGEMS evaluation models material and energy
da Cenibra Durante periodo de
um balances have displayed an average relative error of 2 4for the con
alive and the
quatro mews entre maio e setembro sumption of vapor 2 5 for vapor economy which are compa
rable to the experimental error fi
om the measurements and therefore the
de oesreferentes
2003 informa a con
simulations have been considered consistent About the heat transfer coe

de vapor taxa de aoe

alimenta
sumo
cient modeling the analysis led to a set of equations capable to describe
de the phenomenon of heat transference inside the plate fallingfzlm evapora
ao
produ licor temperatures e pres
tors Among the
soes a cede efeito entre outras foram equations cited on the literature Kandiklar 1990 was the
one that leads to the best results This model displays values of average
coletadas de hora em hora De posse relative errorfor the effects hI to IB of the industrial evaporation system of
desses dados foi elaborada uma ana approximately rangingfi om 2 to 11 and about 44for the first unit IA
These values indicate that the
lise da performance de opera
ao da
seem to
proposal of this researcher can be
used to describe the phenomenon of thermal exchange in the
unidade auxilio do simulador evaporation
com o
bodies ofKraft recovery processes
comercial WinGEMS 3
5 no qual os

Referencias dos autores

Departamento de Engenharia Quimica da Universidade Federal de Minas Gerais UFMG Belo Horizonte MG
Escola Politecnica da Universidade de Sao Paulo USP Sao Paulo SP
3Celulose Nipo Brasileira S
A Belo Oriente MG
parametros nao mensurados direta vez acarreta ao da eflciencia
redu da Conforme esquematizada na Figu
caldeira de unidade II de evapora
mente Como por exemplo o consu ao
recupera com ao
redu re 1 a aoda
de vapor aode soli da de vapor da eflciencia industrial da Cenibra e consti
mo a concentra ao
produ e
planta
dos coeficiente de do sulfato de sodio sul tuida de seis efeitos contendo
sews e o proprio glo conversao a um to

bal de transferencia de calor feto de sodio assim Como aumento das tal de oito corpos de evapora
U para aodo
cada corpo de aopuderam
evapora emissoes de et Al fllme descendente de places
ser
SOZ e
HZS Adams tipo sen

determinados e analisados Os mode 1997 do oprimeiro efeito constituido por

los de ao do
avalia simulador Win Uma das alternatives para se con tres corpos No primeiro efeito a ope
GEMS propostos mostraram
se ade tornar o problema imposto pelo acres ra de
ao lavagem continua e conduzi

quados uma vez que forneceram um cimo da viscosidade e a aode tec

ado da em um dos corpos enquanto os

Erro Relativo Medio ERM de 2

4 nicas de aode
redu viscosidade tra outros doffs operam normalmente Ao

os calculos do consumo de va tamento termico ou ao

adi de sais flm de determinado periodo de tempo
para

por vivo e de 5
2 para a economic que conduzam a ganhos reais na con a lavagem e alternada de modo a pro
de vapor Quanto modelo do coefi aofinal do licor negro Toda
ao centra porcionar uma
limpeza adequada em
ciente de transferencia de todos
global ca via para precisar esses ganhos tor os tres equipamentos Salienta
lor analise efetuada levou se necessario relacionar viscosi e efetuada normal
a a umcon na a se que a
lavagem
de oes capazes
equa de descre dade do licor da uni licor negro fraco
junto no
desempenho mente com
aproxi
ofenomeno de troca termica dade de evapora madamente
ver na ao 20m3h e o produto da
aofor
convec ada dentro dos evapo A complexidade dos os que
balan lavagem e enviado para concentra
ao
radores de filme descendente de pla envolvem sistemas de evaporadores de no terceiro efeito onde a corrente se

ces Entre as varies oes

proposi da multiplo efeito inibem oes
avalia une ao produto do quarto efeito Se

literature testadas neste de presen de corpos pa

trabalho a quantitativas A guindo ofluxo normal contra
corren
a
Kandiklar 1990 foi a que conduziu ralelos tanques
flash pre e aque
pos te com o vapor
cos melhores resultados As cedores reciclos entre outros elemen Seguindo o fluxograma da Figure
oes
equa
obtidas renderam valores de erro re tos torna os calculos manuais muito 1 se que
tem o licor negro fraco isto

lativo medio para os corpos do Siste diflceis A aodos

determina coeflci e com baixa concentra
aode solidos
de evapora
aoda unidade indus entes de transferencia de calor e alimentado nos sexto
ma
globais quinto e quar
trial da ordem de 11 de evapora
IB ao VI a 2 e a capacidade aoSao da to efeitos nas oes aproxima
propor
de ate 44 para dos de vital na busca des das de 43 50 7 respectiva
e o
primeiro corpo importancia e

de evapora
aoIA condi otimas de
oes ao
opera mente O licor segue contra
Seguindo a mo corrente

ao efetuada
diflca na aode
proposi ao fluxo do vapor e ao sair do primei
Kandiklar 1990 Esses valores efeito expandido
ace aodo Sistema de Evapora
Avalia ao ro e em um tanque
nam que a aodesse pesqui
proposi O primeiro passo da avalia da
ao flash onde atinge a aofi
concentra
sador pode ser utilizada na ao
descri unidade foi o levantamento do fluxo nal antes de ser enviado cos tanques
do fenomeno de troca termica para grama da de evapora
aoII da
os planta de estocagem de licor concentrado

corpos de evapora
aode IB ao VI Cenibra S
A Celulose Brasi
Nipo Quanto as oes
adi de vapor nos

leira A
S localizada a 230 km de Belo evaporadores de multiplo efeito ocor

INTRODUfAO Horizonte no municipio de Belo Ori rem em tres entradas independentes

Um dos problemas enfrentados pe ente MG De posse do fluxograma e localizadas no primeiro efeito sendo
las unidades de recupera
aodo pro entrada para cada corpo de
acompanhando a ao na
opera planta uma eva

Kraft refere
se a viscosidade do seas
estabeleceram ao O vapor
industrial efeito
cesso cone pora gerado nesse

licor negro du da uni cinco efeitos restantes

e seu comportamento xoes entre os equipamentos segue para os

rante a sue concentra nos evapora

ao dade tais Como os corpos de evapo O condensado obtido em cada efeito

dores de multiplo efeito O acrescimo ao

ra e os tanques de evapora
ao e por sue vez encaminhado a tan

no teor de solidos sews ocorrido com flash Com a sequencia estabelecida ques de expansao flash e o vapor ob
a evapora da
ao ague altera as carac se um
elaborou fluxograma no simu tido e reaproveitado juntamente com

teristicas do licor omais

tornando lador WinGEMS versao 3 3 uma
5 o evaporado Como meio de aqueci
viscoso e pegajoso bem Como afetan vez a complexidade dos calculos mento do efeito seguinte
que
do diretamente transferencia de que envolvem sistemas de evapora adi de
a ca Outra importante ao calor
ao
lor de evapora
aono Sistema
e a taxa
multiplo efeito inibe qualquer avali ao Sistema ocorre com a ao de
adi
de 5oSoderhjelm 1986
concentra aoquantitative manual Edwards coluna de do
a vapor na
lavagem con

Um licor menos concentrado por sue Baldus 1978 densado contaminado com Gases

Pacific Simulation Inc 2004

1999 WinGEMS programa de aopara
simula industrais de papel e celulose versao 3www
5com
pacsim
tem
da de
Figura 1 Fluxograma planta aoda
evapora unidade II da Cenibra

Nao Condensaveis GNC Essa co ao do

densa vapor que retorna para finidos pertencentes a biblioteca do

luna recebe todo o condensado de a coluna de separa

ao e o GNC ob programa WinGEMS semelhante a

cada corpo de evapora

aoe efetua a tido e enviado para queima na cal abordagem de Edwards Baldus

separa dos gases atraves de lava deira de biomassa 1978 e Cardoso et al 1998 As
ao

gem por contato direto com vapor de A Tabela 1 apresenta os dados de oestomadas
simpliflca envolveram a

baixa pressao A aodo sistema de evapora

aoda aodas
elimina linhas de baixa
GNC
misturavapor opera vazao

com grande quantidade de energia unidade II da Cenibra que nao interferem nos os
balan de

de e aodo
elimina ci
passivel aproveitamento enca massa e
energia e a

minhada para de aodo Sistema de Evapora

Simula ao clo de do corpo sob lava
um
conjunto pla lavagem e
exclusivo no terceiro onde A
cas efeito partir do fluxograma levantado gem no primeiro efeito Na sequen
ocorre uma troca de calor com o li se um
criou diagrama de blows no si cia de lavagem do primeiro efeito
cor do proprio efeito Ali se da a con mulador que envolviam blows prede como os tres corpos A B e C Sao

iguais nao se justifica a ado des

ao
sa alternancia no simulador Assim
Va or Licor
o fluxograma foi adotado com apenas
Vazao 74600 h
kg 600 m3
483 h doffs corpos no primeiro efeito um

7C
151 1C
86 concentrando produto do Segundo
Temperatura o

efeito e o outro enviando o licor con

Pressao 6 kgf
3 2
cm centrado aos tanques de estocagem
Teor de solidos dis 4
15
Assim foi criada virtual
uma planta
Perdas Termicas capaz de avaliar com detalhes a ope

1 e 2 efeitos ra da unidade de evapora

ao aoA
0
1
Figura 2 apresenta a aoda
representa
do 3 ao 6 efeito 5
0
planta no simulador WinGEMS

Tabela 1 Dados de opera do sistema de

ao evapora da unidade II
ao Foram especificados os dados das

da Cenibra correntes de entradas de licor

 peratura vazao e aoas
concentra
Variaveis
temperatures do evaporado de cada

vazao do vapor
efeito e o teor de solidos sews do li Temperature pressao e
cor concentrado final se
obtendo vivo no primeiro efeito 03

Como resposta as oesva

concentra Temperature e aodo licor fraco
concentra 02
zoes pressoes coeficientes
e
globais
de transferencia de calor de cada cor ao
Distribui do licor fraco entre os quarto quinto e sexto efeitos 03

po de evapora
aoalem do consumo
Temperature concentra
aoe vazao de licor concentrado 03
total de vapor e a economic do siste
Temperature do licor temperature do vapor e
pressao
ma A aofoi
especifica conduzida
cada corpo de evapora 21
dados
em ao
se os
compartilhando previa
mente formatados planilha ele Vazao de vapor de baixa vazao de concensado 02
numa
pressao
tronica do Excel com o WinGEMS
Pressao na coluna de 01
stripper
atraves de um
arquivo programado
Temperature e vazao da do condensador de superficie 02
em linguagem Visual Basic Opro ague
grama envia as oes da pla
informa aodo licor
Concentra na Saida do Segundo efeito 01
nilha ao WinGEMS e epos o forne
38
cimento de todos os dados necessa

rios inicia a ao pare

itera conver

gencia do calculo O criterio de con

vergencia e o alcance da concentra

a

ao de solidos sews no licor concen

1

trado final com o valor fornecido

t TI
durante a entrada de dados Antes de

finalizar o programa testa a conver

do WinGEMS a
alcan
gencia e se
1

da retorna os valores calculados a

do Excel Esse
uma outs planilha
1

e
procedimento repetido varies ve I
1

xes ate que todos os conjuntos de

1

dados sejam avaliados

ao do
Valida Modelo de Analise

Para a ao dos
valida modelos

contidos no simulador WinGEMS

foram coletados os dados horarios

de opera
ao da unidade II da Ceni

bra
pressoes
se temperatures
incluindo

e vazoes de vapor e licor

1

4
9 Y
a N ss
E

bem Como as oes de

informa con

aode solidos
centra sews do licor
a

Saida do efeito li
na Segundo e no

N cor negro concentrado final forne

0
cidas pelos refratometros da unida

de O periodo de analise estendeu

c se de maio a setembro de 2003 com

z
um volume de aoem
informa torno
a de 2
800 conjuntos de dados em
W que
cada
i
conjunto corresponde as me

dics horarias
O registradas por equi
pamento somando um total de 38 Figure 2 Fluxograma do primeiro ao quinto efeitos da plants de evapora

variaveis a saber ao da unidade II da Cenibra representado no simulador WinGEMS

 Os dados foram obtidos a partir Portanto se necessario
faz o co Assim procurou
serelacionar nes

dos relatorios de 115 nhecimento do termo referente lado te trabalho os dados da viscosidade do
impressos cerca ao

mil Em
dadosseguida foram tabu dolicor ja que para o lado do vapor licor das variaveis operacionais e das
lados em uma planilha eletronica do condensante diversos trabalhos tem variaveis de projeto com o U atraves
Excel e preparados para a utiliza sido na literatura para des do
ao publicados ajuste dos parametros das equa
oes
pelo programa WinGEMS Os dados cri do fenomeno
ao Collier Tho desenvolvidas por Chen 1966 Shah

relativos a woes de inicio de opera rne 1994 e Chen 1961 Para o coefl 1976 Gungor Winterton 1986
ao bem Como os de oes do
interrup ciente de transferencia de calor no lado Steiner Taborek 1992 Chun

processo por estarem ligados a pe do licor seComo

considerara um fato Seban 1971 e Alhusseini et al 1998

riodos transientes foram seleciona estabelecido no presente momento que conforme as oesapresentadas na
equa
dos e eliminados do conjunto com o doffs mecanismos principais participam Tabela 2 Os parametros coeficientes
intuito de evitar transmissao de calor das
erros nos os
balan na em escoamen e expoentes oes foram
equa de

de massa e energia to bifasico em dutos fechados Kli terminados por analise de regressao

menko 1988 O primeiro e a ebuli nao linear avaliados pelo metodo dos

Equacionamento do Coeficiente Glo ao nucleada caracterizada pela pre minimos quadrados

bal de Tranfferencia de Calor a de sitios ativos de
sen aona
nuclea
Uma vez validado o modelo de si parede aquecida bem Como por uma RESULTADOS

ao
mula os valores das variaveis que forte dependencia do Coeficiente con ao do Simulador
Valida WinGEMS

nao Sao avaliados na planta industrial vectivo de transferencia de calor h com Corn base em trabalhos recentes

puderam ser obtidos entre os quaffs as o fluxo de calor q semelhante ao pro publicados por Almeida et al 2003 e

vazoes dos evaporados e os teores de cesso de aoestagnada

ebuli No outro 2005 a aodo modelo proposto
valida
solidos sews assim como os proprios mecanismo denominado ao
vaporiza foi conduzida a partir da analise dos

valores dos coeficientes globais de ao for

por convec ada bifasica os erros relativos medios e dos coeficien

transferencia U de calor em cada sitios ativos Sao completamente supri tes de correla R obtidos
ao na analise

corpo de evapora
ao De posse dos midos e o calor e transferido por eva de regressao O erro relativo medio
valores dos coeficientes de ao na interface
pora fllme li dado e calculado
globais entre o
pela equa
ao2 a par
transferencia de calor obtidos ovapor Na de tir da diferen do valor de
a uma varia
a
partir quido e
regiao vapori
do Simulador ao a transferencia de calor
WinGEMS vel obtida da leitura
se ao
passou za depen a partir na planta
equacionamento para aode uma
obten de principalmente do fluxo massico e industrial com a mesma variavel cal

aoque relaciona
equa as
propriedades da qualidade do vaporz culada pelo Simulador WinGEMS Para

flsicas do licor os dados operacionais Diversos trabalhos tern sido publi o consumo de vapor vivo no primeiro
de projeto U dos evaporado cados tentativa de formular efeito Figura 3 ERM foi de
e com o na se uma o 4
2
res de fllme descendente de placas ao generica para
equa o calculo do

Sendo U o Coeficiente global de Coeficiente convectivo h em sistemas

transferencia de calor um bifasicos Algumas das maiores difl

que para rl
r
Y

equipamento em opera sem a

ao pre culdades envolvidas nessa ao
formula
Figura 4 exibe
sen de
a oes
deposi e oes e
incrusta Sao crescimento e desprendimento das A uma compara
ao
dado inverso do somatorio das bolhas Campo de escoamento da valores de economic de vapor
pelo no entre os

resistencias taffy i aodas dugs

misturabifasica distribui reais unidade industrial obti
como pela con na e os

lado do vapor termodinamico dos apartir do Simulador WinGEMS

vec no
ao condensante fries equilibrio nay

ii pela parede metalica da superflcie oeslocais

condi caracteristicas e ge se as
Analisando Figuras 3 e 4
de aquecimento e iii pela convec
ao ometric da superflcie de aquecimento se que
percebe de forma geral os da

dos do de
no lado do licor em aquecimento Sen e efeitos das propriedades do fluido no comparativos consumo va

do do pro sistema de evapora economic de vapor estao ali

assim o equacionamento aoUm modelo por e a

blema somente e possivel realistico que englobe today nhados

como co essay com proximidade ou mesmo se
na

nhecimento da contribui de cada

ao plexidades ainda nao foi desenvolvi sobrepondo em pontos na linha dia
termo da equa
ao 1 do Kandiklar 1990 A Tabela 2 apre gonal Portanto os valores calculados
senta resumidamente 10 importantes pelo Simulador estao proximos cos
i

oes disponiveis
contribui na literatu encontrados na unidade industrial
1 1
J re no
equacionamento dos coeficien Essa aopode
aflrma ainda ser corro

tes
globais de transferencia de calor borada pelos baixos valores do ERM

para diversas oes

solu encontrados que foram de 4 para
2

z
Na ao e
ebuli na
aoe frequentemente conveniente o use da fra
condensa ao
de vapor em
ao ao fluxo massico total detinida como qualidade de vapor
rela M
escoamento
r tYatiY3
lnwesliymJr rNa rrg es
c
4

py sve y

trm rlgrl Uc
r fir hy E i
lr
n
r nral
lr
r tlt
fi
inl
Y R
r u
ItIY

5
yr
G1GkYr nP3el ued
na 1
6avcirn
Inl
I 4 1
Irdlrii
r
x
C ilt
L
l SdeY
ulif
S llTdl
d Leo r
l rim 6t
rtiiu
Lam tl
9
YiL
l
et lixn
c reiJu er
e
Jsxs
p
Iluic7ix r
r
h
esiy
P

rl lac i g
r
liar

ry
de
niml fi Fig T
rrl ar
hlmea rin I
um rrra
nPx7
h
Titi sn5 Pt
p
FaUll
n d1
rnrt
8r Tyi
B n
ua tuba rrrrlc
r s kE
A Per
ea rrPl I
InpP rC
I eu h
arel aeale nathva
ra
r
4
r
ilrI
1 E nlc
l
rr
s5

Pr rA rJU PrUrrtill c

J1rAd5d SBFtY
i 13
Jridrralit
riuS rr
s
n

D e
s dx3 I
rati

jnrwirl
rivna nVNfC
C
QrSU en 9r
1
h 19 E
a
I re4 c
q 9f
1
F

nam
C
nl Jt k
Sri FUrrle
Cel
e
era rtrr
a
n
r
O
I
l
l 1
r rY I
7 ri
h
b
l I t
rel A
1
7 l151P
vi
rxKi dnsv
s
11 12 R
113
8 l
CC
duf
1a I t
r 41 rs ajur rJ
r
I
VI l
6 5
hIG
tae rta F
a
r
1 P x
I r1
7
inxnlaa
rxpr
y

gl9o
I 11
lld111
qu 5
1
d

ade
i
Rc hrvx nualdade I R
tlS
Y
esvro mwJ
fJ a a1

lral
nerrlasr
to
1i s
fr p l ra1
eer s
r
3r prapoGla
rel
er
rZ 3 IAaAI
nror cr h
f n
n idfo 3
r dodos r
C hAnr
1rar
sY11ty1k5Ul
r
k7 itit
u
gl7rpna mrryr d8
I
r
r
ralr sesinix S
r
ii IrI
rlr 4Q l
ku1
dd
I
rsla

Qua FaC 415 ETrfS

C 117 Ki tl iD trig
q 9k7uli
Jrr rJrr
1
rir
rer
tl
rk errxuivislQ

II r I
I P A w
7J e f r 1 I 5Pr
I
r i 1

en7e uertae
E al e

ia
h
a
r l anal
n 1 1a e efah 3d
a1 mi
a4 gy ntsa
3
ton
au X1
rr 7 I lel I1 1
J vela rarrrmla4m d
9 tip Irk L
91Ca rr
11
1 p

114 7
r
n
nr e re wn h
ry ma
r
r7
ranl iaxrsl eaieg rrarrv
irr
1 I
IG Frrrtir
rl
ii
a
nxyRrrlm
v sYi Il i
l
A
7
2 l7 dr

Tabela 2 Algumas correla

oes importantes para a aosob
ebuli
ao
Jer
4r7are5lili 9
rt
Qe4rJ ei
4e
rt dC
r
Glr

k
F
om
r
E
cal
verl
narllr rt
ry
1 e ti
i
4
hciriYmial
ua FI
1 I9x A
N vianirdia
s
de
n ne ral
r
1I
FI p 91 A 7
3
1 pans Jocu
1
dI tid F
Kd
r
aGFnno Ir ersladau
t cal lr
nsl r6i
fisi
a

a
panl r
tidn flr rh
sa
re
n
n
rx
ena l 1l
lran druli id
E
S
i di r5d
l

ri

ti LF
9
lii
IIP lIL

rsp ryial u
d l nr
rdC
r
ICIGLngx d
iraursl fa
fc5r
c
rr
rlc 1l
r3 rn
l 2
lho
tS
L
1
irrloria5
@ti wl
lrow
yr
xeip s
nis
acaP IrarsP
Fk n
cakartissrwaxh
fie rr91 IIll a 51r11
r
EI
Jf
ia r53 Pr ur

7
ralriartev
rg lnr

IrrFrflclua a
en
par
a CInr4 lab
lYa 9
r
t7 r 5 p R11 YS r
dYytnf
nl
t II l
fd9
Frr 1
iP 7
Y 7l rl e r K
d
Ci inrrrsa da
jr na

oarnrnlavrvtical a
Esr hul raaclma
rasrzmial
I xs rn
u
vraa rl
na 4
r
agdia H I 1 Ii IH xrnra us d
15
s1 21 a 281 s

snSY
c
rlrs
97
1 FI tf l l rl I I
wpars i os

hrt rceel ir
nsl rr
iPeysir
I
3
13
FI aya r
o 1 Irdr
3
t Y iJkJrvu
o
v
Y tarri nrpn
rxrrrx orFdo a
C
lrr
r
aranlr
dM da
ionaao
r Ir r I
Ir
w riad r
smh a
queaiW rn
p
s rd abuli
a rmnr
h
p
r
n
5rrl
kal d
I i9i pnrtlas
r kSpirla
Ja S
17
E n do Frauds Rru h eeahnh
ac
rcia
l I R
lP
11 rr4
G
p fddQBfi7l J
5lluidr
rnlaYa hew
prtla r

4r xi
ul flnl
Y
n
I I
lr
rrt
xvrrf
x
eal
1re
er
ns e
I s 1

Canderrlc
d
Ikl t fi t it f fri
irk

ii
s
pordf r
p

pare raualru riw ip

r l4r
Iq 1lll rr11
9 i
11

rrlriprr I R
f 1 1 Fl alrry
G
I f
Pri sp
J
ar
eal
ryrrrF
c
Rc
yarla
r2 P1
11 4i
G
r r n ue
r Ii
iier R h
C rnf
lerhrJrx
3
a
faels
lr
ur 4eid 5 e
r is
i rp5nr dy
r3 B
g fF
ul
7
i9
r6
a rt
rtdrt
tiL darra
K ssnerrlg esln
e9Go sle
r ut7lira
F
I
d h p
tiiI C r7d e d119 89
1
an4 rr
Id slarxrl 44r1r
d
J
l
If
l
3
t
@
3
1
yr li 1
rl
vlan a L
II Imo r I za v
b
fe
S vfl
Ir UC
a
j 19i C4n1
fr1 5 rxi
d r
l
Jf
piBi
iie hejrm
nn
rcr
li

b
14 u I t Sri
zanyanha dilrna

gua
de a

rar
s i I
htol
pr4piFrrn itireldllvu rtzl
VE
h rx IrarasF ca4ar rroareanh5ro
pr drr
s
rna
rtt
lir
aq s a rrrstanl 1
iuta vcrraesrl
7
ti c Pi 8
rnes farninar
Fl e

ie
7urhular
irIGarmrGa
r
rwmax
da
1 F prvrlrx
s

0
r
apl1f gr
4
r5C

Tabela 2 Algumas correla

oes importantes pars a ao sob
ebuli escoamento Continua
 0
10D
0
90

a 0
95

c L
0
90 0
75
a
a D
ss Q
o
p o
0
60
c
A
O
SD
O 8 gj

0
45
0 0
75 a

m
0
70
a o
a 0
30
0
65
p
0 O a o

0
15
0
6D
c o

c
V 0
55
0

0
50 D 5D0 1D00 1500 2D00 25D0 3D00

0
50 0
60 0
70 0
80 0
90 0
100

consumo de vaporvivo real tlh e Vapor real o

Vapor predito

3 aoentre
Compara valores reais a preditos pelo simulador WinGEMS para o consumo de vapor vivo
Figura os

o
aD
D
s a
U
c 0
78
a

8
5

o
u 0
76

o w
R 6
5 r
o

a
o

o ao w a 0
74
o

C 4
5
o o a
w

0
72
R o0 0 o ac
pp O o o
w U
p 2
5 b
o

o 0
7D
m o
U
o
O
O
O O
41 O
W 0
5
0
68
R

8
4 i
a 0
66
v
c

v
0
64
6
4
0
64 0
66 0
68 0
70 0
72 0
74 0
76 0
78 0
80
6
4 8
4 0
5 2
5 4
5 6
5 8
5 0
6

aoreal
Concentra de sdlidos secos da licor conce ntrado

Economic real

Figura 4 aoentre
Compara os valores reais a estima Figura 5 aoentre
Compara os valores reais a es

dos para a economia de vapor do sistema de ao

evapora timados para o teor de solidos secos no licor ne

gro concentrado

o consumo de vapor e ERM de 5

2 Isso e importante aa va
refor os erros foram de 8e 4
7 9 e para
porque
economia de vapor lida do modelo utilizado
ao efeito ERM
para a o
Segundo Figura 7 o

Pela analise da Para efeito de foi de Isso mostra que balan

Figura 5 percebe aoprocurou
valida se 1
2 o

WinGEMS tambem de pressao dos o efetuado WinGEMS foi satis

se que o convergiu apro produzir os perfis pelo
priadamente para a aode
concentra longo do sistema de eva
evaporados ao fatorio pois conseguiu reproduzir ade
0
solidos
N secos
estipuladana entrada de ao Na Figura 6 e representada a
pora quadamente os resultados esperados
0
dados Cabe mencionar novamente aodos perfis
compara de pressao predi Ainda a fim de se verificar a valida

na metodologia adotada o simu tos pelo simulador WinGEMS e deter de do modelo empregado foi realizacla
que
c
lador foi configurado a efetuar os ba minados na planta industrial para calla uma analise de variancia Anova para
Z
os
lan de modo a atingir determina corpo de aoPode
evapora se observar os pares reais e preditos dal variaveis
a
da concentra
aode solidos secos es nessa figura os valores preditos se chave de valida
aoda simula do pro
ao
W que

pecificada na entrada de dados Sen encontram muito pr6ximos cos obtidos cesso Entretanto os testes baseados na

do essa aoatingida
concentra em to na unidade industrial studentrealizados paracom
estatisticat
O
das resultado apre dos
as oes o
simula Para as pressoes evaporados ao dal
para medial reais e preditas para
da ao
simula dos corpos A B do as variaveis economia de vapor con
senta a convergencia e primeiro efeito
 1800 D
1150
s
1600 O
e E
E
1400 D
1050 g
c e8
1200
A 88
o

1000 D
950
R
v 0
Boa 0 0
a

soD
D
sso a
v

aoD a
o
A
200 a
D
750 ao
0 a
IA IB II III N V VI

D
650

Predito Real D
650 0
750 0
85D 0
950 0
1050 0
1150
o e

PressSo real da evaporado II mmHg

6 Perfis de pressao do evaporado real ecal 7

Figura Figura aoentre
Compara os valores reais a estima

culado ao da unidade de ao
evapora
longo dos para a
pressao do evaporado no segundo efeito

de vapor vivo no do evapora

sumo
primeiro efeito pressao
do do efeito concentra do licor do
ao
segundo e
segundo 80
efeito nao foram conclusivos Para a aodesses testes
realiza 70 0
rn
e aoessencial que seja satisfeita a premissa
condi de distribui
60
1
normal dal dos residu 0
ao oes e consequentemente
observa 50
N
os residuo yr a ynr r Todavia
a a analise desses residuos
a a 40
demonstrou o contrario Os residuos obtidos nao apresenta 0
i
n o
N
30 0 p
vam aonormal
distribui afetados talvez pela descalibra
ao a o 20
nos instrumentos de medi ao longo
ao da campanha de opera
a 10
ao da unidade Dessa forma sem a aoda condi
satisfa ao D
basica na qual baseia a analise de variancia dal
se a qualidade IA IB I III N V VI
estatisticas de teste foram
comprometidas de modo a compro
meter tambem a conclusao quanto a rej ei
ao ou nao da hipote o Calalado o Real
se de igualdade dal medial obtidas pelas leituras dos instru

mentos de ao da unidade
medi com as mesmas
preditas pelo
Perfil de concentra
aode solidos secos dissol
modelo do WinGEMS Figura 8

Alem disso destaque ao

simula vidos no licor negro ao da unidade de ao
evapora
outro ponto merece a longo
foi efetuada em um processo industrial e o modelo dese

nhado foi concebido para ser capaz de predizer as varias Saida de cada corpo de evapora
aoNa unidade estudada

oesdesse
flutua processo A aoda unidade
opera nao ocor existem apenas dois refratometros que medem a concentra

unico patamar sendo composta

re em um na verdade por ao de solidos sendo um na Saida do segundo efeito e outro

varios conjuntos de estados operacionais sobrepostos Cada na Saida do licor concentrado final apos o tanque de con

tem media variancia aoflash

centra do efeito Considerando
se que
patamar sua
propria e e mesmo as
primeiro nas

leituras dos instrumentos de medi da unidade fornecem

ao leituras fornecidas dois refratometros da unidade ha
pelos
valores de erro que dependem da grandeza da variavel lida despreziveis devido aos cuidados e ao rigor da cali
erros

O processo tem varias distribui

oeso que acaba por preju ao efetuada pela equipe de manuten
bra aoindustrial ve
dicar a realiza do teste
ao t
s tudent que compara dugs situ se
riflca mail uma a qualidade do modelo elabora
vez

oes
a em que cada uma dal oes comparadas
popula tem do visto que os valores medidos na Saida do segundo efeito

uma unica media e aonormal

distribui e no tanque flash de produto final concordaram com os va

Assim se
tomando Como valido o modelo do WinGE lores preditos pelo simulador

etapa seguinte envolveu analise do desempenho da A da analise da compara

MS a a partir e aodos dados obti
unidade industrial A Figura 8 ilustra o perfll da concentra dos na planta industrial com os preditos pelo simulador
ao de solidos secos dissolvidos presentee no licor negro na se que
concluiu o programa WinGEMS se mostrou ade
 4000

3000

s
2000

Y
1000

0 500 1000 1500 2000 2500 3000

3000 1500

I I

2500 0 00
o 0

lU
08 0
0 0 0

U
1250

O O g c
O
2000
O
pOaO
L O
g O

S
O
O
o
O CQ o0
O 00 O
L
1000
ooap 00
0

R 1500 m
Y Y

o 0
0
8
0 0 00 0 750
1000

500 500

0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000

3000 3000

2500 2500

U 7
t

t s

2000 2000

R R
V V
Y Y

1500 7 1500

1000 1000

0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000

o Real Media
2000

0
17so
0
N

0
F 1500

t
1250
o
9 Coeficiente de transferen
Z Figura global
Y 1000
cia de calor do efeito
a primeiro ao sexto

W
750

500

0 500 1000 1500 2000 2500 3000

quado ao desenvolvimento dos ba

os
lan de masse e energia e conse
t

quentemente no levantamento das W

I
variaveis de da unida n
performance d y a aFf n r n

de industrial de evapora
ao
A da ao do
valida simu Equa 3
ao aoempirica proposta
Modifica a
equa de Kandiklar 1990
ao
partir para

lador e da ao das
obten variaveis

de performance do sistema de eva

so

pora foram obtidos

ao os coefici

entes globais de transferencia de 40

c

calor para cada corpo de evapora

ao do sistema conforme ilustrado c 30

na Figura 9
20
i
i
A valo
primeira observa
aoos o

w
de U obtidos para
res os evaporado o

de filme descendente de
res
places
o
concordam com os descritos na li
IA IB II III IV V VI
Venka
terature Segundo Nguyen Efeita
tesh 1992 e Foust 1982 A segun
Chen1966 Shah 1982 Gungor e on1986
Winter

da a dispersao dos pontos aumenta Kandiklar 1990 Steiner e Taborek 1992 Chun e Seban 1971
Alhusseini et al 1998
com o acrescimo do teor de solidos

dissolvidos O acrescimo de solidos

10 Erros relativos medios avaliados conforme a ao 2 obti
equa
modifica o comportamento reologi Figura
dos para todas as equa testadas
oes
co do fluido e aumenta a incidencia

de aoque
deposi sue vez afeta
por
o processoprovocando oscila
oes com as oes contidas
equa na Tabe Uma ao na Figura
inspe 10 e na

U A la 2 foram obtidos parametros Tabela 3 revela que para efeitos

periodical propria sistema
no os os

dal de numero III modifica efe

tica de lavagens periodical no pri oes coeficientes
equa eexpo a VI a ao
meiro efeito tambem e uma fonte entes ATabela 3 e Figura 10 ilus tuada na ao de
proposi Kandiklar

de ela tram os valores dos Errol Relati 1990 se mostrou promissora equa
periodica oes
oscila e a se

Medios dos coeficien

credita a cause da large dispersao vos ERM ao 3 Para esses corpos de evapora
observada para B do Upara cada corpo de evapora ao os ERMs encontrados foram bai
os corpos A e tes

efeito e na Figura 11 Sao apresenta da ordem de 2 a 4

primeiro ao xos enquan
dos valores dos coeficientes de
os to para os corpos IB e II os ERMs

Equacionamento do Coeficiente U correla R obtidos

ao na analise de foram mail elevados de respectiva
Para cada uma dal equa tes
oes regressao para os pares de U pre mente 11 e 7 Isso significa que

tadas dito U avaliado via WinGEMS efeitos

combina da equa
ao ao 1 e para os em que a ao
opera
ocorre em oes de temperature
condi
Efeito
e teor de solidos dissolvidos menos
Pesquisador
IA IB II III IV V VI elevados e o licor apresenta um com

36 13 9 6 5 4 4 portamento reologico newtoniano a

Chen 1966
ao efetuada
proposi se mostrou apta
Shah 1982 36 14 8 8 5 4 4
a descrever o fenomeno de transferen

Winterton cia de calor Ja para primeiro corpo

Gungor 1986 44 13 8 5 27 23 4 o

2 3 3 lA a aodo
predi Coeficiente U nao
Kandiklar 1990 44 11 7 4
foi eficaz de Arazao dis
ERM 44
Steiner Taborek 95 437 11 7 5 4 4
1992 so esta ligada ao fato de que nesses

Chun Seban 1971 31 13 12 5 5 5 5 corpos a ao e mail drastica

opera

Alhusseini et al 1998 35 13 9 5 27 4 5 maiores temperatures e concentra

oes de solidos dissolvidos e o licor
Tabela 3 Erros relativos medios avaliados conforme a ao 2
equa processado tem comportamento tipi
obtidos para as equa testadas
oes camente pseudoplastico
 Os valores de R apresentados na
100

90 Figura 11 tambem corroboram as in

oes
forma da Tabela 3 da
so e
Figura
o 10 relativas aos valores de ERM Da
so
mesma forma a analise revela que a
50
melhor obtida e basea
40 ao
proposi a

s
30 da no modelo de Kandiklar 1990
20
cujos valores de R foram da ordem
10
de 56 a 88 para os efeitos de nu
o
mero III a aode
VI cuja concentra
IA IB II III IV V VI

licor e mais baixa e

Efeito consequente
mente o escoamento do fluido e mais
Chen 1966 Shah 1982 Gungor e Wirnerton 1986
Kandiklar 1990 Steiner e Taborek 1992 Chun e Seban 1971 bem comportado Ou seja para os
Alhusseini et al 1998
efeitos de numero III a VI de 56 a

88 dos valores da variavel U obti

Figura 11 Coeficientes de correla R2 obtidos para todas as equa

ao dos pelo WinGEMS puderam ser ex
oes testadas aoefetuada
plicados pela modifica
no modelo de Kandiklar Esses va

lores podem ser considerados bons

L

c 8000 em se tratando de um fenomeno de

U transferencia de calor avaliado a par
a

tir de uma ao empirica

equa
6000
Na analise dos altos valores de
c
d

ERM para os corpos IA e IB 44 e

U
11 respectivamente baixos
R 4000 e va

lores de R alguns pontos merecem

s
d

O refere
destaque primeiro se ao ca

2000
s v ster altamente transiente da opera
0

ao nessa parte da evapora

aoA ins
d

0 ao da
pe Figura 12 revela o compor

tamento altamente oscilatorio da

6
24 6
24 6
25 6
25 6
26 6
26 6
27
do
00 00
12 00 00
12 00 00
12 00 ao e
opera consequentemente
coeficiente global de transferencia de

calor nesse estagio

Na analise da ao de U
oscila o
L

4000 efeito da deposi

aode sais de burqui
ta precisa ser levado em considera
d

ao A ao de
deposi solidos afeta mui
3000
to ofenomeno de transferencia de

U calor uma vez que cria uma resisten

ul o

2000 cia adicional a troca termica Como

i efeito
s
esse perturba a aonor
opera
0 d

N mal principalmente no primeiro efei

0 v 1000
to nesse a qualidade dos dados obti
0

C7 dos pelo WinGEMS e afetada e por

d

c
0 consequencia e afetado tambem o
z d

valor da ao de
predi U pelo modelo
6
24 6
24 6
25 6
25 6
26 6
26 6
27
a d
matematico Alem disso a
W 00 00
12 00 00
12 00 00
12 00 ao
equa
matematica utilizada na aode
predi
U nao apresenta o termo relativo a
O
12 ao do
Oscila coeficiente de transferencia de calor du resistencia oferecida filme de
Figura global pelo
rante a
ao para
opera os corpos de evapora IA
ao a a IB b deposito A medida que o
equipamen
to opera as placas de troca de calor a do comportamento da viscosida rencia de calor de cada corpo de eva

resistencia adicional a de aoda concentra

altera obtidos
se sujam e a com a ao ao puderam
pora ser

troca de calor aumenta de solidos Quanto modelo do coeficien

gradativa secos ao

mente Com a opera

aode lavagem Na analise do erro obtido para o te global de transferencia de calor
o filme de deposito e removido e a fenomeno oscila analise efetuada levou
Segundo efeito o a a um con

resistencia por sua vez e elimina torio proporcionado pelas lavagens junto de oes capazes
equa de des

da Apesar de ser um fenomeno os nao pode ser considerado pois nes crever ofenomeno de troca termica

cilatorio constante de efeito existe nenhuma ope ada dentro dos

o processo se nao na aofor
convec

ao
deposi e altamente dependente ra de
ao lavagem e tampouco o flui evaporadores de filme descendente

das oes de polpa

condi aoda madei do apresenta comportamento estri de placas Entre as varies proposi
ra aodo
composi licor e de opera tamente pseudoplastico Todavia oes da literature testadas neste tra
ao da planta de evapora com nesse efeito pelo nivel de solidos balho a de Kandiklar 1990 foi a
aoe

isso a cinetica de aoe

incrusta com secos dissolvidos ados pela
alcan que conduziu cos melhores resulta

plexa eprecisa ser

equacionada para aoo fluido
concentra se encontra dos As oes obtidos
equa renderam

termo in de ao reologi
transi valores de R para ultimos efei
que esse posse entao ser em uma zone os

a equa
ao de U Nos niveis de TDS T ali tos do sistema de evapora
ao da
corporado McCabe ca e en

unidade industrial
Robinson 1924 McCabe 1926 contrados surge novamente o pro III ao VI da
blema aoprecisa
determina do ordem de 56 87
Epstein 1979 na a Segundo a pro
O Segundo ponto diz respeito a valor da viscosidade do fluido apar ao
posi de Kandiklar 1990 Tais

viscosidade do fluido nesses doffs tir de fonte utilizada Zvinakevicius valores evidenciam que a proposi
corpos de evapora
aoA et al ao e
primeiros 1979 algo acon adequada a aodo
descri feno
que passa

proposi inicial
ao de Kandiklar se Como a principal
figurar fonte de meno de troca termica nesses cor

1990 foi testada inicialmente en erro ao de

responsavel pela obten pos de evapora
aoPara os primei
tre outros fluidos com agua cujo um ERM de 7 ros efeitos IAI em virtude das ca

comportamento reologico e tipica Para o entendimento dos erros racteristicas particulares do licor

mente newtoniano Ja fluido obtidos de forma

o em
geral
uma um negro com seu comportamento
e agua ponto tambem omodelo adotado
questao nao pure e
possui pode ser negli
nao
pseudoplastico
de solidos a
uma carga crescente genciado a geometric abordadanos nao foi capaz de descrever adequa
medida que avan Kandiklar de damente ofenomeno
a pelos efeitos experimentos 1990
principalmente nos corpos IA e IB senvolveu sua ao em
equa experi Alem de municiar para o equa
Essa ano perfil
mudan de solidos mentos realizados em evaporadores cionamento do problema de trans

conforme ja mencionado anterior de tubos verticais lisos diferente ferencia de calor nos evaporadores
mente altera as propriedades liga mente das placas observadas nos da unidade industrial o modelo de

das ao escoamento do fluido em equipamentos da Cenibra Essa di simula do WinGEMS revelou

ao
que difi apode tambem
especial a viscosidade o feren ser obj eto fon uma ao direta
aplica e de extrema

culta ainda mais o equacionamento te dos desvios entre os dados expe importancia no ambiente fabril

do problema rimentais e os avaliados pela equa Como a opera de

ao lavagem dos

O licor negro de eucalipto utili ao empirica corpos do primeiro efeito e condu

zado fluido caracteristi zida
e um com a partir de uma premissa tem
cas
particulares Em baixas concen CONCLUSAO poral a ao do
inser modelo do

oes de solidos secos

tra seu compor O modelo de aodo
avalia simu WinGEMS no ambiente operacio
tamento reologico e semelhante ao lador WinGEMS proposto mostrou nal e de extrema valia na determi

da agua mas em niveis maiores de se satisfatorio uma vez que forne ao do ponto otimo para
na lavagem
solidos dissolvidos a viscosidade ceu erros relativos medios de 4
2 dos equipamentos Com o modelo

a depender da tensao de cisa os calculos do consumo de va elaborado e possivel a ao

condu
passa para
lhamento Essa a de
mudan com por vivo e de 5 para
2 a economic dos os
balan de masse e energia
portamento certamente nao esta en de vapor economic essa com media para aodos
determina valores de

globada fator dependente do de 3

5 padraoigual
desvio 08
0 U tempo real para cada corpo
no
Ffl e a em

fluido proposto por Kandiklar O modelo proposto permitiu uma de ao

evapora Analisando
seadi

Dessa forma analise mais detalhada de toda ao do

minui coeficiente U
1990 possivel que a e o

a ao necessite da inclusao de
equa unidade de evapora
aoe os valores acrescimo do diferencial de tempe

fator que considere mudan dos coeficientes de transfe rature determiner

algum a
globais se
pode com
mais acertividade o memento ide possivel e submeter o licor concen
TDS Teor de solidos secos dissolvidos

al para alternancia dos corpos trado obtido efeito

em no Segundo ou
Coeficiente de transferencia
U global
ao reduzindo assim inter corpo intermediario do
opera em um
pri
de
oes
ven lavagem desnecessari meiro efeito a um tratamento para

as no processo as quaffs prejudicam ao de viscosidade

redu tratamento ad Dado de processo obtido
eu
Yea
estabilidade da opera termico
a ao esse que pode ser ou mes
pelo programs WinGEMS
Em virtude dos resultados alcan mo por meio da tecnica de ultracen
yreal Dado de processo obtido
ados com a analise de regressao a aoou adi
trifuga ao de sais Frede
partir dos dados do simulador Win rick 1988 Ferreira et al 2003 na unidade industrial

GEMS foi observado que e possi Com isso operando com licor de

vel descrever o fenomeno de trans viscosidade mais baixa o primeiro

GREGAS
ferencia de calor nos evaporadores efeito pode lever a ao de
obten um

de licor negro de eucalipto operan produto mais concentrado que por

do filme descendente sobre vez teria direto
com
pla sua
impacto e
posi aP
4H Entalpia de aoda ague KJ Kg
vaporiza
se necessario tive
ces Para tanto terns na performance da caldeira de
Condutividade termite C
m
W
um conhecimento detalhado do ao Obviamente
recupera para

comportamento reologico do liter mensurar esses

ganhos o conheci Viscosidade dinamica s
m
Kg
se a
incluindo de mento do fenomeno de transferen
empregado curve v Viscosidade cinematica p

viscosidade para faixas de licores cia de calor na aoe condi

evapora
p Densidade m3
kg
com altas oes de
concentra solidos ao essential

secos dissolvidos e em altas tempe

6 Tensao superficial N
m

ratures de processamento NOMENCLATURA

De posse do modelo do Coefici
ADIMENSIONAIS
ente global de transferencia de ca
VARIAVEIS
associado modelos do simu
Bo Numero de ao
ebuli pm
x
q
lor aos

lador se
fez Co Numero ao x
deconvec P
y
possivel umaavalia
ao
dos e das Vazao massica
impactos consequencias de m mz
kg s Fr Numero de Froude m a
g
i
oes nas propriedades fisico
altera x Qualidade do vapor fra
ao massica
Ka Numero de Kapitza y
R
quimicas do liter negro ou mesmo de em ao a
rela masse total do Fluido
vapor
da configura Pr Numero de Prandtl p
ao do processo no fe
cp Calor especifico C
Kg
J
nomeno de transferencia de calor Re Numero de Reynolds para
economic da evapora Diametro hidraulico 4
na aoe con dh m
e
S
o escoamento laminar
na queima do liter
sequentemente E Fator de ganho
e na germao de vapor na caldeira
Ffi Parametro dependente do fluido
de recupera SUBSCRITOS
ao
g aoda gravidade m
Acelera sz
Uma aplica direta da meto
ao
desenvolvida h Coeficiente de transferencia
dologia na ao
otimiza
Bi Valor observado no ponto de inicio
de plantas de evapora
aopode in

cluir por exemplo a analise do im do regime de ao complete

ebuli
Espessura da parede metalica
de tecnicas de redu de vis
ao
pacto c Termo ligado ao efeito convective

cosidade
da superficie de aquecimento m
no processo de evapora
ao
f Relative a parte liquida do fluido ri
z
i
de li Numero de dados operacionais
Como apresentado a queima n

de alts concentra Licor 0

cores aode soli q Calor transferido do vapor condensante N

dos acima de 80 leva a um au LO Fluido escoa totalmente

0

liter que
para o h
J
mento das eficiencias energeticas e
come uma unite fase liquida
NB
q Calor critico minimo necessario ao inicio y
de reaproveitamento de quimicos na 0

NcB Termo ligado ao efeito de ao

nuclea
caldeira de recupera
aoA limita da ao J
nuclea h

ao de
obten tais licores Calor ambiente Nu Nusselt
a
ao para Perdido para o W
qP
reside hoje na dificuldade de con

loa medida que

centra sua viscosi
nos evaporadores J
h tp Two phase fluido bifasico
a
S Fator de O
dade supressao vapor liquido
aumenta com a ao do
eleva
teor de solidos Uma alternative T Temperature C vap Vapor
LITERATURA CITADA 1966 S G An improved
velopment 5 3 p 329
322 KANDIKLAR
ADAMS T
N FREDERICK CHEN MM An analytical stu correlation for predicting two
phase
J GRACE T
Wm M HUPA M dy oflaminar film condensation part flow boiling heattransfer coefficient in

LISA K JONES A
K TRAN H 1 flat Transactions of the horizontal and vertical tubes em Heat
plates
a
Kraft Recovery Boilers 1 ed Ta
ASME 2 54
48 1961 exchangers for phase
two flow appli
p

ppi Press New York 1997

CHUN K R SEBAN R A cations ASME Nova Iorque 1983

ALHUSSEINI A A TUZLA K Heat transfer to evaporating liquid KLIMENKO V V A generali

CHEN J C film
Falling evaporati films Journal of Heat Transfer 11 zedcorrelation for phase
two forced

of Int
on
single component liquids p 396
391 1971 flow heat transfer Int Journal of

Journal of Heat and Mass Transfer COLLIER J G THOME J R Heat and Mass Transfer 31 3 p

41 12 p 1623
1632 1998 Convective and condesati
boiling 552
541 1988

ALMEIDA G M CARDOSO on 3a ed Clarendon New York MCCABE W L ROBINSON

M OLIVEIRA E D PARK S W
596p 1994 C S Evaporator Scale Formation
Predi da germao de vapor
ao e da
EDWARDS L L BALDUS R Industrial and Engineering Chemis

razao tonelada de vapor gerado por F effect evaporation

Multiple sys try 16 5 p 478 479 1924

tonelada de solidos sews em caldei tems for kraft black Evalua MCCABE W L Economic Side
liquor
ran de recupera
aovia modelos de tion and design Tappi 61 3 of Evaporator Scale Formation
pp
redes neurais O Papel 64 10 p 78
75 1978 Chemical and Metallurgical Engine
105
99 2003
EPSTEIN N Optimum Evapora ering 33 2 p 86 87 1926

ALMEIDA G M BARBOSA tor with Scale Formation The NGUYEN X N VENKATESH

Cycles
C M L CARDOSO M OLIVEI Canadian Journal of Chemical Engi V Evaporation and concentration of

RA E D PARK S W Analise de
neering 57 8 p 659 661 1979 black liquor em Chemical recovery

dados operacionais da caldeira de re das in the alkaline pulping ed

FOUST A S Principios process 3a

cupera
aopor meio das tecnicas de
oes Unitarian
Opera 2 ed Guana 1992
Tappi Press Atlanta pp 35
5
visualiza de
ao dados aode
sele va bara Dois Rio de Janeiro 670p 1982 SHAH M M Chart correlation

riaveis de redes neurais O

e
Papel GUNGOR K E WINTERTON for natured boiling heat transfer equa

66 8 p 73 2005
66 R H S A general correlation for tions and further study ASHRAE

BJORGE R W HALL G R flow in tubes and annuli Trans 88 1 196

185 1982
boiling p

ROHSENOW W M Correlation of Int Journal of Heat and Mass Trans L of

SODERHJELM Viscosity
forced convection boiling heat trans 1986 Black Puu
fer 29 3 p 358
351 Strong Liquor Paperi ja
fer data Int Journal of Heat and Mass och Tra
GUNGOR K E WINTERTON Paper 9 pp 652
642 1986

Transfer 25 6 757
753 1982 R H S J
p Simplified general corre STEINER D TABOREK
CARDOSO M PASSOS M L lation for flow and Flow heat transfer in verti
boiling compari boiling
A CARAZZA F SANTOS S M sonsof correlation The Cana cal tubes correlated
data by an asympto

aopreliminar
Avalia do processo de dian Journal of Chemical tic Heat Transfer N
Enginee model Enginee 0
Y

aodo licor de
recupera eucalipto ring 65 1 p 156
148 1987 ring 13 2 p 69 1992
43
para sua aoO Papel 10
otimiza KANDIKLAR S G A general ZVINAKEVICIUS C SAN y
0

pp 84
80 1998 correlation for natured F C E B Ca
phase
two TOS P FOEKEL
a
CHEN J C Correlation for boi flow heat transfer inside racteristicas de licores kraft
boiling pretos W

ling heat transfer to natured fluids in horizontal and vertical de Relatorio interno se
tubes eucalipto
O
convective flow Industrial Engine Journal of Heat rie processamentos n
Transfer 2 112 102 Cenibra
Process De
ering Chemistry Design p 228
219 1990 Belo Oriente llp 1979