FlashPT

ex1.m

fncex1.m

Arquivos MATLAB:

Flash PT

FTfPf

z1, z2, z3, z4

Tv=TfPv=Pf

Tl=TfPl=Pf

Flash PT

Etano (substancia 48)

Propano (substancia 64)

n-butano (substancia 87)

n-pentano (substancia 104)

Composição da corrente de entrada:

Flash PT

%Valores especificados:

F=100; % vazão de alimentação (kgmole/hr) Tfvetor=230:1:300; % temperatura (Fahrenheit) Pf=420; % pressão (psia) z1=0.10; % composição da alimentação z2=0.20; % composição da alimentação z3=0.30; % composição da alimentação% z4=0.40; Tv=350; % temperatura da saída vapor (Fahrenheit) Pv=450; % pressão da saída vapor (psia)

Flash PT

Flash PT

Como localizar o ponto de cruzamento?

Flash PT

>> ginput(1)

ans = 277.9839 49.8538

Flash PT

Líquido – Vapor > 0 ? 1111000

Ponto de cruzamentoAtenção:nunca use a regra

“liquido – vapor = = 0”

Flash PT

Ampliando a faixa de simulação de Tf:

%Valores especificados:

F=100; % vazão de alimentação (kgmole/hr) Tfvetor=100:1:400; % temperatura (Fahrenheit) Pf=420; % pressão (psia) z1=0.10; % composição da alimentação z2=0.20; % composição da alimentação z3=0.30; % composição da alimentação% z4=0.40; Tv=350; % temperatura da saída vapor (Fahrenheit) Pv=450; % pressão da saída vapor (psia)

Flash PT

?!

?!

Flash PT

Flash PT

Flash PT

“Para uma dada pressão, a temperatura selecionadano flash PT deve ficar entre a Tbolha e a Torvalho.”

Flash PT

Temperaturade bolha

ex2.m

Arquivos MATLAB:

fncex2.m

Temperatura de bolha

Temperatura de bolha

curva debolha

pressãofixa

Zoom:

Temperatura de bolha de 214.4 F para uma pressão de 420 psia.

Temperatura de bolha

Temperaturade orvalho

ex3.m

Arquivos MATLAB:

fncex3.m

Temperatura de orvalho

curva deorvalho

pressãofixa

Temperatura de orvalho

Zoom:

Temperatura de orvalho de 307.8 F para uma pressão de 420 psia.

Temperatura de orvalho

Temperaturade bolha e orvalho

ex4.m

Arquivos MATLAB:

manipulab4.m

manipulaa4.mteste8b.m

teste6b.m

Temperatura de bolha e orvalho

orvalho.m

bulbo.m

Temperatura de bolha e orvalho

Flash PT

Tbolha Torvalho214F 307F

230F 300F

100F 400F

Exemplo 1

Exemplo 2

Temperatura de bolha (influência da composição)

ex5.m

Arquivos MATLAB:

teste5b.m

bulbo.m

Temperatura de bolha (influência da composição)

Temperatura de bolha (influência da composição)

z1=0.10; % composicao da alimentacao

z2=0.20; % composicao da alimentacao

z3=0.30; % composicao da alimentacao

cor='-r';

z1=0.00; % composicao da alimentacao

z2=0.40; % composicao da alimentacao

z3=0.20; % composicao da alimentacao

cor='-y';

z1=0.00; % composicao da alimentacao

z2=0.40; % composicao da alimentacao

z3=0.00; % composicao da alimentacao

cor='-k';

z1=0.00; % composicao da alimentacao

z2=0.00; % composicao da alimentacao

z3=0.50; % composicao da alimentacao

cor='-c';

Temperatura de bolha (influência da composição)

Temperatura de orvalho (influência da composição)

ex6.m

Arquivos MATLAB:

teste7b.m

orvalho.m

Temperatura de orvalho (influência da composição)

z1=0.10; % composicao da alimentacao

z2=0.20; % composicao da alimentacao

z3=0.30; % composicao da alimentacao

cor='-r';

z1=0.00; % composicao da alimentacao

z2=0.40; % composicao da alimentacao

z3=0.20; % composicao da alimentacao

cor='-y';

z1=0.00; % composicao da alimentacao

z2=0.40; % composicao da alimentacao

z3=0.00; % composicao da alimentacao

cor='-k';

z1=0.00; % composicao da alimentacao

z2=0.00; % composicao da alimentacao

z3=0.50; % composicao da alimentacao

cor='-c';

Temperatura de orvalho (influência da composição)

Temperatura de orvalho (influência da composição)

Flash T Beta (variando beta)

ex7.m

Arquivos MATLAB:

fnc3.m

Flash T Beta (variando beta)

%Valores especificados:

F=100; % vazão de alimentação (kgmole/hr)

Tf=230; % temperatura de alimentação (Fahrenheit)

fivetor=0.1:0.01:0.99; %fi

z1=0.10; % composição da alimentação

z2=0.20; % composição da alimentação

z3=0.30; % composição da alimentação

Flash T Beta (variando beta)

Flash T Beta (variando beta)

Flash T Beta (variando beta)

Flash T Beta (variando temperatura)

ex8.m

Arquivos MATLAB:

fnc3.m

Flash T Beta (variando temperatura)

%Valores especificados:

F=100; % vazao de alimentacao (kgmole/hr)

Tfvetor=230:1:300; % temperatura de alimentacao (Farenheit)

fi=0.6; %fi

z1=0.10; % composicao da alimentacao

z2=0.20; % composicao da alimentacao

z3=0.30; % composicao da alimentacao

Flash T Beta (variando temperatura)

Flash T Beta (variando temperatura)

Flash T Beta (variando temperatura)

Flash T Beta (variando temperatura

e beta)

ex9.m

Arquivos MATLAB:

fnc3.m

Flash T Beta (variando beta e temperatura)

%Valores especificados:

F=100; % vazao de alimentacao (kgmole/hr)

Tfextensao=230:1:300; % temperatura de alimentacao (Farenheit)

fiextensao=0.1:0.01:0.99; %fi

[Tfvetor,fivetor]=meshgrid(Tfextensao,fiextensao); %cria a malha

z1=0.10; % composicao da alimentacao

z2=0.20; % composicao da alimentacao

z3=0.30; % composicao da alimentacao

Flash T Beta (variando beta e temperatura)

www.scicomput.com

x = -1:0.01:1

y = -1:0.01:1

[X,Y] = meshgrid(x,y);[X,Y]

Flash T Beta (variando beta e temperatura)

Flash T Beta (variando beta e temperatura)

Flash T Beta (variando beta e temperatura)

Flash T Beta (variando beta e temperatura)

Flash P Beta (variando beta)

ex10.m

Arquivos MATLAB:

fnc4.m

Flash P Beta (variando beta)

Flash P Beta (variando beta)

Flash P Beta (variando beta)

Flash P Beta (variando pressão)

ex11.m

Arquivos MATLAB:

fnc4.m

Flash P Beta (variando P)

Flash P Beta (variando P)

Flash P Beta (variando P)

Flash P Beta (variando pressão e beta)

ex12.m

Arquivos MATLAB:

fnc4.m

Flash P Beta (variando P e beta)

hvfcnhlfcn

Flash P Beta (variando P e beta)

Flash P Beta (variando P e beta)

Flash P Beta (variando P e beta)