Flash PT Flash PT Arquivos MATLAB: ex1.m fncex1.m Flash PT Tv=Tf Pv=Pf F Tf Pf z1, z2, z3, z4 Tl=Tf Pl=Pf Flash PT Composição da corrente de entrada: Etano (substancia 48) Propano (substancia 64) n-butano (substancia 87) n-pentano (substancia 104) 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 ? Atenção: nunca use a regra “liquido – vapor = = 0” 1 1 1 1 0 0 0 Ponto de cruzamento 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 selecionada no flash PT deve ficar entre a Tbolha e a Torvalho.” Temperatura de bolha Temperatura de bolha Arquivos MATLAB: ex2.m fncex2.m Temperatura de bolha curva de bolha pressão fixa Temperatura de bolha Zoom: Temperatura de bolha de 214.4 F para uma pressão de 420 psia. Temperatura de orvalho Temperatura de orvalho Arquivos MATLAB: ex3.m fncex3.m Temperatura de orvalho curva de orvalho pressão fixa Temperatura de orvalho Zoom: Temperatura de orvalho de 307.8 F para uma pressão de 420 psia. Temperatura de bolha e orvalho Temperatura de bolha e orvalho Arquivos MATLAB: ex4.m manipulab4.m manipulaa4.m teste8b.m orvalho.m teste6b.m bulbo.m Temperatura de bolha e orvalho Flash PT Tbolha Torvalho 307F 214F Exemplo 1 230F 300F Exemplo 2 100F 400F Temperatura de bolha (influência da composição) Temperatura de bolha (influência da composição) Arquivos MATLAB: ex5.m teste5b.m bulbo.m Temperatura de bolha (influência da composição) z1=0.10; % composicao da alimentacao z1=0.00; % composicao da alimentacao z2=0.20; % composicao da alimentacao z2=0.40; % composicao da alimentacao z3=0.30; % composicao da alimentacao z3=0.00; % composicao da alimentacao cor='-r'; cor='-k'; z1=0.00; % composicao da alimentacao z1=0.00; % composicao da alimentacao z2=0.40; % composicao da alimentacao z2=0.00; % composicao da alimentacao z3=0.20; % composicao da alimentacao z3=0.50; % composicao da alimentacao cor='-y'; cor='-c'; Temperatura de bolha (influência da composição) Temperatura de orvalho (influência da composição) Temperatura de orvalho (influência da composição) Arquivos MATLAB: ex6.m teste7b.m orvalho.m Temperatura de orvalho (influência da composição) z1=0.10; % composicao da alimentacao z1=0.00; % composicao da alimentacao z2=0.20; % composicao da alimentacao z2=0.40; % composicao da alimentacao z3=0.30; % composicao da alimentacao z3=0.00; % composicao da alimentacao cor='-r'; cor='-k'; z1=0.00; % composicao da alimentacao z1=0.00; % composicao da alimentacao z2=0.40; % composicao da alimentacao z2=0.00; % composicao da alimentacao z3=0.20; % composicao da alimentacao z3=0.50; % composicao da alimentacao cor='-y'; cor='-c'; Temperatura de orvalho (influência da composição) Flash T Beta (variando beta) Flash T Beta (variando beta) Arquivos MATLAB: ex7.m 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 temperatura) Flash T Beta (variando temperatura) Arquivos MATLAB: ex8.m 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 e beta) Flash T Beta (variando beta e temperatura) Arquivos MATLAB: ex9.m 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) x = -1:0.01:1 y = -1:0.01:1 [X,Y] [X,Y] = meshgrid(x,y); www.scicomput.com 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) Flash P Beta (variando beta) Arquivos MATLAB: ex10.m fnc4.m Flash P Beta (variando beta) Flash P Beta (variando beta) Flash P Beta (variando pressão) Flash P Beta (variando P) Arquivos MATLAB: ex11.m fnc4.m Flash P Beta (variando P) Flash P Beta (variando P) Flash P Beta (variando pressão e beta) Flash P Beta (variando P e beta) Arquivos MATLAB: ex12.m fnc4.m hvfcn hlfcn Flash P Beta (variando P e beta) Flash P Beta (variando P e beta) Flash P Beta (variando P e beta)