python ciclo rankine

•
UTFPR

André Cruz
27/10/2021
Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados
16 milhões de materiais de várias disciplinas
Impressão de materiais
Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes
E aí, curtiu este material?
Ajude a incentivar outros estudantes a melhorar o conteúdo
Gostou desse material? Compartilhe! 🧡
Termodinâmica

28.559 Materiais compartilhados
Baixe o app para aproveitar ainda mais
Leia os materiais offline, sem usar a internet. Além de vários outros recursos!
Prévia do material em texto
Entrega 1 - GDV
Alunos: André F. da Cruz, Reinaldo de Vaz Lima
Programa Ciclo Rankine - cogeração
Tabela 1 - Dados operacionais da planta termoelétrica
Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título
1 485 6495 27.9 -
2 - 900 - -
3 - 250 - -
4 51 Saturação - -
5 51 Saturação - 0
6 - 250 - -
7 110 250 - -
8 - 6495 27.9 -
9 110 250 - -
10 - 900 - -
11 Saturação 900 ≤ 6.94 1
12 Saturação 900 ≤ 6.94 0
 
 ---------- ESTADO 1 - Entrada da Turbina ---------- 
 
 
Temperatura no ponto 1 -> T1 = 485.00 °C 
 
Pressão no ponto 1 -> P1 = 6495.00 kPa 
 
Entalpia no ponto 1 -> h1 = 3381.20 Kj/kg 
 
Entropia no ponto 1 -> s1 = 6.79 Kj/kg 
 
Densidade no ponto 1 -> ρ1 = 19.65 Kg/m^3 
 
Título no ponto 1 -> X1 = -1.00 
 
 
 
 ---------- ESTADO 2 - Primeira Extração da Turbina ---------- 
 
 
Temperatura no ponto 2 -> T2 = 516.42 K 
 
Entalpia no ponto 2 -> h2 = 2932.08 Kj/kg 
 
Entropia no ponto 2 -> s2 = 6.95 Kj/kg 
 
Densidade no ponto 2 -> ρ2 = 3.91 Kg/m^3 
 
Título no ponto 2 -> X2 = -1.00 
 
 
 
 ---------- ESTADO 3 - Segunda Extração da Turbina ---------- 
 
 
Temperatura no ponto 3 -> T3 = 127.41 °C 
 
Entalpia no ponto 3 -> h3 = 2714.68 Kj/kg 
 
Entropia no ponto 3 -> s3 = 7.05 Kj/kg 
 
Densidade no ponto 3 -> ρ3 = 1.39 Kg/m^3 
 
Título no ponto 3 -> X3 = 1.00 
 
 
 
 ---------- ESTADO 4 - Última Extração da Turbina ---------- 
 
 
Pressão no ponto 4 -> P4 = 12.98 kPa 
 
Entalpia no ponto 4 -> h4 = 2333.12 Kj/kg 
 
Entropia no ponto 4 -> s4 = 7.26 Kj/kg 
 
Densidade no ponto 4 -> ρ4 = 0.10 Kg/m^3 
 
Título no ponto 4 -> X4 = 0.89 
 
 
 
 ---------- ESTADO 5 - Saída do Condensador ---------- 
 
 
Temperatura no ponto 5 -> T5 = 51.00 °C 
 
Pressão no ponto 5 -> P5 = 12.98 kPa 
 
Entalpia no ponto 5 -> h5 = 213.52 Kj/kg 
 
Entropia no ponto 5 -> s5 = 0.72 Kj/kg 
 
Densidade no ponto 5 -> ρ5 = 987.54 Kg/m^3 
 
Título no ponto 5 -> X5 = 0.00 
 
 
 
 ---------- ESTADO 6- Saída da Bomba e Entrada do Desareador ---------- 
 
 
Temperatura no ponto 6 -> T6 = 51.02 °C 
 
Pressão no ponto 6 -> P6 = 250.00 kPa 
 
Entalpia no ponto 6 -> h6 = 213.81 Kj/kg 
 
Entropia no ponto 6 -> s6 = 0.72 Kj/kg 
 
Densidade no ponto 6 -> ρ6 = 987.64 Kg/m^3 
 
Título no ponto 6 -> X6 = -1.00 
 
 
 
 ---------- ESTADO 7 - Saída do Desaerador ---------- 
 
 
Temperatura no ponto 7 -> T7 = 110.00 °C 
 
Pressão no ponto 7 -> P7 = 250.00 kPa 
 
Entalpia no ponto 7 -> h7 = 461.49 Kj/kg 
 
Entropia no ponto 7 -> s7 = 1.42 Kj/kg 
 
Densidade no ponto 7 -> ρ7 = 951.00 Kg/m^3 
 
Título no ponto 7 -> X7 = -1.00 
 
 
 
 ---------- ESTADO 8 - Entrada da Caldeira ---------- 
 
 
Temperatura no ponto 8 -> T8 = 110.75 °C 
 
Pressão no ponto 8 -> P8 = 6495.00 kPa 
 
Entalpia no ponto 8 -> h8 = 469.21 Kj/kg 
 
Entropia no ponto 8 -> s8 = 1.42 Kj/kg 
 
Densidade no ponto 8 -> ρ8 = 953.43 Kg/m^3 
 
Título no ponto 8 -> X8 = -1.00 
 
 
 
 ---------- ESTADO 9 - Entrada da bomba para HX ---------- 
 
 
Temperatura no ponto 9 -> T9 = 110.00 °C 
 
Pressão no ponto 9 -> P9 = 250.00 kPa 
 
Entalpia no ponto 9 -> h9 = 461.49 Kj/kg 
 
Entropia no ponto 9 -> s9 = 1.42 Kj/kg 
 
Densidade no ponto 9 -> ρ9 = 951.00 Kg/m^3 
 
Título no ponto 9 -> X9 = -1.00 
 
 
 
 ---------- ESTADO 10 - Saída da Turbina ---------- 
 
 
Temperatura no ponto 10 -> T10 = 110.08 °C 
 
Pressão no ponto 10 -> P10 = 900.00 kPa 
 
Entalpia no ponto 10 -> h10 = 462.30 Kj/kg 
 
Entropia no ponto 10 -> s10 = 1.42 Kj/kg 
 
Densidade no ponto 10 -> ρ10 = 951.25 Kg/m^3 
 
Título no ponto 10 -> X10 = -1.00 
 
 
 
 ---------- ESTADO 11 - Entrada do Processo Industrial Vizinho ---------- 
 
 
Temperatura no ponto 11 -> T11 = 175.35 °C 
 
Pressão no ponto 11 -> P11 = 900.00 kPa 
 
Entalpia no ponto 11 -> h11 = 2773.03 Kj/kg 
 
Entropia no ponto 11 -> s11 = 6.62 Kj/kg 
 
Densidade no ponto 11 -> ρ11 = 4.65 Kg/m^3 
 
Título no ponto 11 -> X11 = 1.00 
 
 
 
 ---------- ESTADO 12 - Saída do Processo Industrial Vizinho ---------- 
 
 
Temperatura no ponto 12 -> T12 = 175.35 °C 
 
Pressão no ponto 12 -> P12 = 900.00 kPa 
 
Entalpia no ponto 12 -> h12 = 742.56 Kj/kg 
 
Entropia no ponto 12 -> s12 = 2.09 Kj/kg 
 
Densidade no ponto 12 -> ρ12 = 891.92 Kg/m^3 
 
Título no ponto 12 -> X12 = 0.00 
 
 
Tabela 1 Atualizada - Dados operacionais da planta termoelétrica
Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título
1 485.00 6495.00 27.9 -1.000000
2 243.27 900.00 - -1.000000
3 127.41 250.00 - 0.999172
4 51.00 12.98 - 0.890762
5 51.00 12.98 - 0.000000
6 51.02 250.00 - -1.000000
7 110.00 250.00 - -1.000000
8 110.75 6495.00 27.9 -1.000000
9 110.00 250.00 - -1.000000
10 110.08 900.00 - -1.000000
11 175.35 900.00 ≤ 6.94 1.000000
12 175.35 900.00 ≤ 6.94 0.000000
In [85]: from sympy import * 
from math import pi 
from IPython.display import display, Math, Latex 
x , y = symbols('x y') 
init_printing() 
from CoolProp.CoolProp import PropsSI as ps 
from CoolProp.CoolProp import State as st 
import CoolProp.CoolProp as coolprop 
import numpy as np 
import pandas as pd 
In [86]: i = [[485,6495, 27.9, '-'],['-',900,'-','-'], 
 ['-',250,'-','-'],[51,'Saturação','-','-'], 
 [51,'Saturação','-',0],['-',250,'-','-'], 
 [110,250,'-','-'],['-',6495,27.9,'-'], 
 [110,250,'-','-'],['-',900,'-','-'], 
 ['Saturação',900,'≤ 6.94',1],['Saturação',900,'≤ 6.94',0]] 
df=pd.DataFrame(i,index='1 2 3 4 5 6 7 8 9 10 11 12'.split(), 
 columns='Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título '.split()) 
In [87]: df 
Out[87]:
In [123… fl = 'Water' 
 
#Estado 1 - Entrada da Turbina 
P1 = 6495 #kPa 
m1 = 27.9 #kg/s 
T1 = 485+273.15 #K 
 
st1 = st(fl,{'P':P1,'T':T1}) 
h1 = st1.h 
s1 = st1.s 
ρ1 = st1.rho 
tit1 = st1.Q 
 
print('\n\n','-'*10,'ESTADO 1 - Entrada da Turbina','-'*10) 
print ('\n\nTemperatura no ponto 1 -> T1 = {:.2f}'.format(T1-273.15),'°C') 
print ('\nPressão no ponto 1 -> P1 = {:.2f}'.format(P1),'kPa') 
print ('\nEntalpia no ponto 1 -> h1 = {:.2f}'.format(st1.h),'Kj/kg') 
print ('\nEntropia no ponto 1 -> s1 = {:.2f}'.format(st1.s),'Kj/kg') 
print ('\nDensidade no ponto 1 -> ρ1 = {:.2f}'.format(ρ1),'Kg/m^3') 
print('\nTítulo no ponto 1 -> X1 = {:.2f}\n\n'.format(tit1)) 
In [105… #Estado 2 - Primeira Extração da Turbina 
 
P2 = 900 #kPA 
 
 #Estado ideal 
st2s = st(fl,{'P':P2,'S':st1.s}) 
 
 #Eficiência da Turbina 
ηTurb = 0.85 
 
 #Estado real 
h2 = st1.h - (st1.h - st2s.h)*ηTurb 
st2 = st(fl,{'P':P2,'H':h2}) 
T2 = st2.T 
s2 = st2.s 
ρ2 = st2.rho 
tit2 = st2.Q 
 
print('\n\n','-'*10,'ESTADO 2 - Primeira Extração da Turbina','-'*10) 
print ('\n\nTemperatura no ponto 2 -> T2 = {:.2f}'.format(st2.T),'K') 
print ('\nEntalpia no ponto 2 -> h2 = {:.2f}'.format(st2.h),'Kj/kg') 
print ('\nEntropia no ponto 2 -> s2 = {:.2f}'.format(st2.s),'Kj/kg') 
print ('\nDensidade no ponto 2 -> ρ2 = {:.2f}'.format(ρ2),'Kg/m^3') 
print('\nTítulo no ponto 2 -> X2 = {:.2f}\n\n'.format(tit2)) 
In [125… #Estado 3 - Segunda Extração da Turbina 
 
P3 = 250 #kPa 
 
 #Estado ideal 
st3s = st(fl,{'P':P3, 'S':st2.s}) 
 
 #Eficiência da Turbina 
ηTurb = 0.85 
 
 #Estado real 
h3 = st2.h - (st2.h - st3s.h)*ηTurb 
st3 = st(fl,{'P':P3,'H':h3}) 
T3 = st3.T 
s3 = st3.s 
ρ3 = st3.rho 
tit3 = st3.Q 
 
print('\n\n','-'*10,'ESTADO 3 - Segunda Extração da Turbina','-'*10) 
print ('\n\nTemperatura no ponto 3 -> T3 = {:.2f}'.format(T3-273.15),'°C') 
print ('\nEntalpia no ponto 3 -> h3 = {:.2f}'.format(st3.h),'Kj/kg') 
print ('\nEntropia no ponto 3 -> s3 = {:.2f}'.format(st3.s),'Kj/kg') 
print ('\nDensidade no ponto 3 -> ρ3 = {:.2f}'.format(ρ3),'Kg/m^3') 
print('\nTítulo no ponto 3 -> X3 = {:.2f}\n\n'.format(tit3)) 
In [107… # Estado 4 - Última Extração da Turbina 
 
T4 = 51 + 273.15 #K 
 
 #Estado ideal 
st4s = st(fl,{'T':T4, 'S':st3.s}) 
 
 #Eficiência da Turbina 
ηTurb = 0.85 
 
 #Estado real 
h4 = st3.h - (st3.h - st4s.h)*ηTurb 
P4 = ps('P','T',T4,'Q',1 , fl)/1000 
st4 = st(fl,{'P':P4,'H':h4}) 
s4 = st4.s 
ρ4 = st4.rho 
tit4 = st4.Q 
 
print('\n\n','-'*10,'ESTADO 4 - Última Extração da Turbina','-'*10) 
print ('\n\nPressão no ponto 4 -> P4 = {:.2f}'.format(P4),'kPa')print ('\nEntalpia no ponto 4 -> h4 = {:.2f}'.format(st4.h),'Kj/kg') 
print ('\nEntropia no ponto 4 -> s4 = {:.2f}'.format(st4.s),'Kj/kg') 
print ('\nDensidade no ponto 4 -> ρ4 = {:.2f}'.format(ρ4),'Kg/m^3') 
print('\nTítulo no ponto 4 -> X4 = {:.2f}\n\n'.format(tit4)) 
In [126… #Estado 5 - Saída do Condensador 
P5 = P4 
st5 = st(fl,{'P':P5,'Q':0}) 
T5 = st5.T 
s5 = st5.s 
ρ5 = st5.rho 
tit5 = st5.Q 
 
print('\n\n','-'*10,'ESTADO 5 - Saída do Condensador','-'*10) 
print ('\n\nTemperatura no ponto 5 -> T5 = {:.2f}'.format(T5-273.15),'°C') 
print ('\nPressão no ponto 5 -> P5 = {:.2f}'.format(P5),'kPa') 
print ('\nEntalpia no ponto 5 -> h5 = {:.2f}'.format(st5.h),'Kj/kg') 
print ('\nEntropia no ponto 5 -> s5 = {:.2f}'.format(st5.s),'Kj/kg') 
print ('\nDensidade no ponto 5 -> ρ5 = {:.2f}'.format(ρ5),'Kg/m^3') 
print('\nTítulo no ponto 5 -> X5 = {:.2f}\n\n'.format(tit5)) 
In [128… # Estado 6 - Saída da Bomba e Entrada do Desareador 
 
P6 = P3 
 
 #Estado ideal 
st6s = st(fl,{'P':P6, 'S':st5.s}) 
 
 #Eficiência da Bomba 
ηBomb = 0.85 
 
 #Estado real 
h6 = ((st6s.h - st5.h)/ηBomb)+st5.h 
st6 = st(fl,{'P':P6,'H':h6}) 
T6 = st6.T 
s6 = st6.s 
ρ6 = st6.rho 
tit6 = st6.Q 
 
print('\n\n','-'*10,'ESTADO 6- Saída da Bomba e Entrada do Desareador','-'*10) 
print ('\n\nTemperatura no ponto 6 -> T6 = {:.2f}'.format(T6-273.15),'°C') 
print ('\nPressão no ponto 6 -> P6 = {:.2f}'.format(P6),'kPa') 
print ('\nEntalpia no ponto 6 -> h6 = {:.2f}'.format(st6.h),'Kj/kg') 
print ('\nEntropia no ponto 6 -> s6 = {:.2f}'.format(st6.s),'Kj/kg') 
print ('\nDensidade no ponto 6 -> ρ6 = {:.2f}'.format(ρ6),'Kg/m^3') 
print('\nTítulo no ponto 6 -> X6 = {:.2f}\n\n'.format(tit6)) 
In [129… #Estado 7 - Saída do Desaerador 
 
P7 = P3 
T7 = 110 + 273.15 #kPa 
st7 = st(fl,{'P':P7,'T':T7}) 
h7 = st7.h 
s7 = st7.s 
ρ7 = st7.rho 
tit7 = st7.Q 
 
print('\n\n','-'*10,'ESTADO 7 - Saída do Desaerador ','-'*10) 
print ('\n\nTemperatura no ponto 7 -> T7 = {:.2f}'.format(T7-273.15),'°C') 
print ('\nPressão no ponto 7 -> P7 = {:.2f}'.format(P7),'kPa') 
print ('\nEntalpia no ponto 7 -> h7 = {:.2f}'.format(st7.h),'Kj/kg') 
print ('\nEntropia no ponto 7 -> s7 = {:.2f}'.format(st7.s),'Kj/kg') 
print ('\nDensidade no ponto 7 -> ρ7 = {:.2f}'.format(ρ7),'Kg/m^3') 
print('\nTítulo no ponto 7 -> X7 = {:.2f}\n\n'.format(tit7)) 
In [130… #Estado 8 - Entrada da Caldeira 
P8 = P1 
 
 #Estado ideal 
st8s = st(fl,{'P': P8, 'S': st7.s}) 
 
 #Eficiência da Bomba 
ηBomb = 0.85 
 
 #Estado real 
h8 = ((st8s.h-h7)/ηBomb)+h7 
st8 =st(fl,{'P': P8, 'H': h8}) 
T8 = st8.T 
s8 = st8.s 
ρ8 = st8.rho 
tit8 = st8.Q 
 
print('\n\n','-'*10,'ESTADO 8 - Entrada da Caldeira ','-'*10) 
print ('\n\nTemperatura no ponto 8 -> T8 = {:.2f}'.format(T8-273.15),'°C') 
print ('\nPressão no ponto 8 -> P8 = {:.2f}'.format(P8),'kPa') 
print ('\nEntalpia no ponto 8 -> h8 = {:.2f}'.format(st8.h),'Kj/kg') 
print ('\nEntropia no ponto 8 -> s8 = {:.2f}'.format(st8.s),'Kj/kg') 
print ('\nDensidade no ponto 8 -> ρ8 = {:.2f}'.format(ρ8),'Kg/m^3') 
print('\nTítulo no ponto 8 -> X8 = {:.2f}\n\n'.format(tit8)) 
In [131… #Estado 9 - Entrada da bomba para HX 
 
st9 = st7 
P9 = P7 
T9 = st9.T 
h9 = st9.h 
s9 = st9.s 
ρ9 = st9.rho 
tit9 = st9.Q 
 
print('\n\n','-'*10,'ESTADO 9 - Entrada da bomba para HX ','-'*10) 
print ('\n\nTemperatura no ponto 9 -> T9 = {:.2f}'.format(T9-273.15),'°C') 
print ('\nPressão no ponto 9 -> P9 = {:.2f}'.format(P9),'kPa') 
print ('\nEntalpia no ponto 9 -> h9 = {:.2f}'.format(st9.h),'Kj/kg') 
print ('\nEntropia no ponto 9 -> s9 = {:.2f}'.format(st9.s),'Kj/kg') 
print ('\nDensidade no ponto 9 -> ρ9 = {:.2f}'.format(ρ9),'Kg/m^3') 
print('\nTítulo no ponto 9 -> X9 = {:.2f}\n\n'.format(tit9)) 
In [132… #Estado 10 - Saída da Turbina 
P10 = P2 
 
 #Estado ideal 
st10s = st(fl,{'P': P10, 'S': st9.s}) 
 
 #Eficiência da Bomba 
ηBomb = 0.85 
 
 #Estado real 
h10 = ((st10s.h-h9)/ηBomb)+h9 
st10 =st(fl,{'P': P10, 'H': h10}) 
T10 = st10.T 
s10 = st10.s 
ρ10 = st10.rho 
tit10 = st10.Q 
 
print('\n\n','-'*10,'ESTADO 10 - Saída da Turbina ','-'*10) 
print ('\n\nTemperatura no ponto 10 -> T10 = {:.2f}'.format(T10-273.15),'°C') 
print ('\nPressão no ponto 10 -> P10 = {:.2f}'.format(P10),'kPa') 
print ('\nEntalpia no ponto 10 -> h10 = {:.2f}'.format(st10.h),'Kj/kg') 
print ('\nEntropia no ponto 10 -> s10 = {:.2f}'.format(st10.s),'Kj/kg') 
print ('\nDensidade no ponto 10 -> ρ10 = {:.2f}'.format(ρ10),'Kg/m^3') 
print('\nTítulo no ponto 10 -> X10 = {:.2f}\n\n'.format(tit10)) 
In [133… #Estado 11 - Entrada do Processo Industrial Vizinho 
Q11 = 1 
P11 = P10 
st11 = st(fl,{'P': P11, 'Q': Q11}) 
T11 = st11.T 
s11 = st11.s 
ρ11 = st11.rho 
tit11 = st11.Q 
 
print('\n\n','-'*10,'ESTADO 11 - Entrada do Processo Industrial Vizinho ','-'*10) 
print ('\n\nTemperatura no ponto 11 -> T11 = {:.2f}'.format(T11-273.15),'°C') 
print ('\nPressão no ponto 11 -> P11 = {:.2f}'.format(P11),'kPa') 
print ('\nEntalpia no ponto 11 -> h11 = {:.2f}'.format(st11.h),'Kj/kg') 
print ('\nEntropia no ponto 11 -> s11 = {:.2f}'.format(st11.s),'Kj/kg') 
print ('\nDensidade no ponto 11 -> ρ11 = {:.2f}'.format(ρ11),'Kg/m^3') 
print('\nTítulo no ponto 11 -> X11 = {:.2f}\n\n'.format(tit11)) 
In [134… #Estado 12 - Saída do Processo Industrial Vizinho 
Q12 = 0 
P12 = P10 
st12 = st(fl,{'P': P12, 'Q': Q12}) 
T12 = st12.T 
s12 = st12.s 
ρ12 = st12.rho 
tit12 = st12.Q 
 
print('\n\n','-'*10,'ESTADO 12 - Saída do Processo Industrial Vizinho ','-'*10) 
print ('\n\nTemperatura no ponto 12 -> T12 = {:.2f}'.format(T12-273.15),'°C') 
print ('\nPressão no ponto 12 -> P12 = {:.2f}'.format(P12),'kPa') 
print ('\nEntalpia no ponto 12 -> h12 = {:.2f}'.format(st12.h),'Kj/kg') 
print ('\nEntropia no ponto 12 -> s12 = {:.2f}'.format(st12.s),'Kj/kg') 
print ('\nDensidade no ponto 12 -> ρ12 = {:.2f}'.format(ρ12),'Kg/m^3') 
print('\nTítulo no ponto 12 -> X12 = {:.2f}\n\n'.format(tit12)) 
In [142… i1 = [[485,6495, 27.9, tit1], 
 [round(T2-273.15,2),900,'-',tit2],[round(T3-273.15,2),250,'-',tit3],[51,round(P4,2),'-',tit4], 
 [51,round(P5,2),'-',0],[round(T6-273.15,2),250,'-',tit6], 
 [110,250,'-',tit7],[round(T8-273.15,2),6495,27.9,tit8], 
 [110,250,'-',tit9],[round(T10-273.15,2),900,'-',tit10], 
 [round(T11-273.15,2),900,'≤ 6.94',1],[round(T12-273.15,2),900,'≤ 6.94',0]] 
df1=pd.DataFrame(i1,index='1 2 3 4 5 6 7 8 9 10 11 12'.split(), 
 columns='Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título '.split()) 
In [143… df1 
Out[143…
In [ ]: 
	Entrega 1 - GDV
	∙ Alunos: André F. da Cruz, Reinaldo de Vaz Lima
	∙ Programa Ciclo Rankine - cogeração
	∙ Tabela 1 - Dados operacionais da planta termoelétrica
	∙ Tabela 1 Atualizada - Dados operacionais da planta termoelétrica
{
 "cells": [
 {
 "cell_type": "markdown",
 "id": "052916da",
 "metadata": {},
 "source": [
 " # Entrega 1 - GDV \n",
 " \n",
 " ## Alunos: André F. da Cruz, Reinaldo de Vaz Lima\n",
 " \n",
 " ## Programa Ciclo Rankine - cogeração"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 85,
 "id": "8490ce73",
 "metadata": {},
 "outputs": [],
 "source": [
 "from sympy import *\n",
 "from math import pi\n",
 "from IPython.display import display, Math, Latex\n",
 "x , y = symbols('x y')\n",
 "init_printing()\n",
 "from CoolProp.CoolProp import PropsSI as ps\n",
 "from CoolProp.CoolProp import State as st\n",
 "import CoolProp.CoolProp as coolprop\n",
 "import numpy as np\n",
 "import pandas as pd"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 86,
 "id": "d99b17f6",
 "metadata": {},
 "outputs": [],
 "source": [
 "i = [[485,6495, 27.9, '-'],['-',900,'-','-'],\n",
 " ['-',250,'-','-'],[51,'Saturação','-','-'],\n",
 " [51,'Saturação','-',0],['-',250,'-','-'],\n",
 " [110,250,'-','-'],['-',6495,27.9,'-'],\n",
 " [110,250,'-','-'],['-',900,'-','-'],\n",
 " ['Saturação',900,'≤ 6.94',1],['Saturação',900,'≤ 6.94',0]]\n",
 "df=pd.DataFrame(i,index='1 2 3 4 5 6 7 8 9 10 11 12'.split(),\n",
 "columns='Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título '.split())\n"
 ]
 },
 {
 "attachments": {
 "trab1.JPG": {
 "image/jpeg": "/9j/4AAQSkZJRgABAQEAYABgAAD/4RDsRXhpZgAATU0AKgAAAAgABAE7AAIAAAAMAAAISodpAAQAAAABAAAIVpydAAEAAAAWAAAQzuocAAcAAAgMAAAAPgAAAAAc6gAAAAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAEFuZHLDqSBDcnV6AAAFkAMAAgAAABQAABCkkAQAAgAAABQAABC4kpEAAgAAAAM2NQAAkpIAAgAAAAM2NQAA6hwABwAACAwAAAiYAAAAABzqAAAACAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMjAyMToxMDowMiAxODoxNDozMgAyMDIxOjEwOjAyIDE4OjE0OjMyAAAAQQBuAGQAcgDpACAAQwByAHUAegAAAP/hCx5odHRwOi8vbnMuYWRvYmUuY29tL3hhcC8xLjAvADw/eHBhY2tldCBiZWdpbj0n77u/JyBpZD0nVzVNME1wQ2VoaUh6cmVTek5UY3prYzlkJz8+DQo8eDp4bXBtZXRhIHhtbG5zOng9ImFkb2JlOm5zOm1ldGEvIj48cmRmOlJERiB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiPjxyZGY6RGVzY3JpcHRpb24gcmRmOmFib3V0PSJ1dWlkOmZhZjViZGQ1LWJhM2QtMTFkYS1hZDMxLWQzM2Q3NTE4MmYxYiIgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIi8+PHJkZjpEZXNjcmlwdGlvbiByZGY6YWJvdXQ9InV1aWQ6ZmFmNWJkZDUtYmEzZC0xMWRhLWFkMzEtZDMzZDc1MTgyZjFiIiB4bWxuczp4bXA9Imh0dHA6Ly9ucy5hZG9iZS5jb20veGFwLzEuMC8iPjx4bXA6Q3JlYXRlRGF0ZT4yMDIxLTEwLTAyVDE4OjE0OjMyLjY1MjwveG1wOkNyZWF0ZURhdGU+PC9yZGY6RGVzY3JpcHRpb24+PHJkZjpEZXNjcmlwdGlvbiByZGY6YWJvdXQ9InV1aWQ6ZmFmNWJkZDUtYmEzZC0xMWRhLWFkMzEtZDMzZDc1MTgyZjFiIiB4bWxuczpkYz0iaHR0cDovL3B1cmwub3JnL2RjL2VsZW1lbnRzLzEuMS8iPjxkYzpjcmVhdG9yPjxyZGY6U2VxIHhtbG5zOnJkZj0iaHR0cDovL3d3dy53My5vcmcvMTk5OS8wMi8yMi1yZGYtc3ludGF4LW5zIyI+PHJkZjpsaT5BbmRyw6kgQ3J1ejwvcmRmOmxpPjwvcmRmOlNlcT4NCgkJCTwvZGM6Y3JlYXRvcj48L3JkZjpEZXNjcmlwdGlvbj48L3JkZjpSREY+PC94OnhtcG1ldGE+DQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIDw/eHBhY2tldCBlbmQ9J3cnPz7/2wBDAAcFBQYFBAcGBQYIBwcIChELCgkJChUPEAwRGBUaGRgVGBcbHichGx0lHRcYIi4iJSgpKywrGiAvMy8qMicqKyr/2wBDAQcICAoJChQLCxQqHBgcKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKir/wAARCAGyAbkDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD6RooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACisLR9MhvNDsbm4nvmmmto5JGF/OMsVBJwHwOT2q5/Ydp/wA9b/8A8GM//wAXQBo0Vnf2Haf89b//AMGM/wD8XR/Ydp/z1v8A/wAGM/8A8XQBo0Vnf2Haf89b/wD8GM//AMXR/Ydp/wA9b/8A8GM//wAXQBo0Vnf2Haf89b//AMGM/wD8XR/Ydp/z1v8A/wAGM/8A8XQBo0Vnf2Haf89b/wD8GM//AMXR/Ydp/wA9b/8A8GM//wAXQBo0Vnf2Haf89b//AMGM/wD8XR/Ydp/z1v8A/wAGM/8A8XQBo0Vnf2Haf89b/wD8GM//AMXR/Ydp/wA9b/8A8GM//wAXQBo0Vnf2Haf89b//AMGM/wD8XR/Ydp/z1v8A/wAGM/8A8XQBo0Vnf2Haf89b/wD8GM//AMXR/Ydp/wA9b/8A8GM//wAXQBo0Vnf2Haf89b//AMGM/wD8XR/Ydp/z1v8A/wAGM/8A8XQBo0VmaOhhn1OASTPHDdBY/NlaQqDDG2MsSerE/jWnQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRUVxdW9nF5t3PHBH/flcKPzNAEtFZbeItOP/HvJJd+htoWlU/8CUbf1qNtZvJP+PXSpF/2rqZYx/47uP6UrpAbFFZel6heXOoXVtexwKYoopF8kk/fLjBJ6/c9B1rUp7gFFFFABRRRQAUUUUAFFFU9YuZLPQ765gIEsNtJIhIzghSR/KgC5RWF9mvv+g3e/wDfuD/43R9mvv8AoN3v/fuD/wCN1n7SJp7ORu0VhfZr7/oN3v8A37g/+N0fZr7/AKDd7/37g/8AjdHtIh7ORc8Pf8ixpf8A15w/+gCtGudtdOurOzhtoNavhFDGsaApAcADA/5Z+1S/Zr7/AKDd7/37g/8AjdHtIh7ORu0VhfZr7/oN3v8A37g/+N0fZr7/AKDd7/37g/8AjdHtIh7ORu0Vm6DPPcaWWupmnkSeaPzGVQWCysoztAHQDtWlWhmFFFFABRRRQAUUUUAFFFFABRRRQAVzfj6yju/BWpvK0o8i1kkVUkKhmCnG4Drjrg8V0lUdV0ey1q1+zaikkkJzlEmeMMCMEHaRkexoAtW//HrF/uD+VSVBZWcNhaJbW3meWmdvmStIeufvMSf1qegDO0z/AJCGsf8AX4v/AKIirRrO0z/kIax/1+L/AOiIq0aACiiigAooooAKKKKACiikJCgliAB1JoAWis1/EOkqxVL6Kdx1S2JmYf8AAUyarXfiZLe0muY9NvpYoUaRnKLGMAZPDkN29KV0Bt0UA5ANY2vhpJtOg82aOOWdg4ilaMsBG5AypBxkChuyuBqXF1b2cJlu544Ix1eVwo/M1QbxFpp/495ZLv0NrC0qn/gSjb+tVINLsLaYSw2cKy/89dgLn6t1NW6z9oK5G2s3kn/HrpUi+91MsY/8d3H9KjafWJvvXFrbL6RQl2H/AAJiB/47Viip52K5Tawkl/4+9Rvp/TE/lY/797f1pYNLsLaXzYLOBJf+evljefq3U1bopXbAKKKKQEWmf8jFf/8AXpb/APoc1bNY2mf8jFf/APXpb/8Aoc1bNbx2KCiiiqAKKKKACiiigArO8Rf8ivqn/XnN/wCgGtGs7xF/yK+qf9ec3/oBoASiiiuQ6wooooAKKKKACiiigBnh3/kFyf8AX3c/+jnrVrK8O/8AILk/6+7n/wBHPWrXUtjle4UUUUxBRRRQAUUUUAFFFFABRRRQAUUUUAFVb/U7HS4kl1O8gs45JBGjzyBFZz0XJ4ycdKtVjeLvDVp4w8JahoWoAeVeQlA2M+W/VXHuGAP4UAR2+r6bp82rXN9qFrbQNeLtlmmVFP7iHoSea3etfIXwS+GF3efGG5i122Kw+GZS9yrDhpgSI19wSC3uF96+tbvUbKwUG+vILYHoZpQmfzNAFiispvENq3FpDd3Z7eVAyqR6h32qfzqNtU1Ob/UWENup6NcT5Yf8BUEf+PUuZIDZorBZdTn/ANfqhjH921hVB+bbj/KoZNItJEY3SyXfHS6laUZ/3WJA/AVPOhXOkpkz+VBJJ/dUt+QqnoJz4c03P/PpF/6AKl1NtmkXj/3YHP8A46asZh2iXt/YW897qd0TLErtHCViUEjOAVG7/wAeqT+xtPZt01styw6NckzMPoXyRU9kuywt1/uxKP0qeudtkiABVAUAAcADtVDX/wDkW9T/AOvSX/0A1oVn6/8A8i3qf/XpL/6AaQHRr9xfpWRrf/H9pP8A18P/AOinrXX7i/SsjW/+P7Sf+vh//RT1vL4RklFFFYCCiiigAooooAKKKKAItM/5GK//AOvS3/8AQ5q2axtM/wCRiv8A/r0t/wD0Oatmt47FBRRRVAFFFFABRRRQAVneIv8AkV9U/wCvOb/0A1o1neIv+RX1T/rzm/8AQDQAlFFFch1hRRRQAUUUUAFFFFADPDv/ACC5P+vu5/8ARz1q1leHf+QXJ/193P8A6OetWupbHK9wooopiCiiigAooooAKKKKACiiigAooooAKwPtup3s1wIbmG1hjmeIBYd0h2nGdxOOev3a365+x4uNQX+7eP8Arg/1qJtpCZWh8PWsd1dXLy3Ms946vcP5pj81goUFgm0HCqByKuW2n2dmxa0tYYWbq0cYBP1PerNFZXbEFFFFIApGGVI9RS0UAUdPudYsdNtrQ2Fi/kQpHu+2uN20AZx5XtS311rF5p1zaiwsYzPE0e/7a525BGceVz1q7RVc7C41V2oFHYYp1FFSAVV1O2a90m7tYyFeeB41LdAWUjn86skgAknAHUmkVldAyMGVhkEHIIoAjGpawFA/s2x4H/P8/wD8aqCZtQv720e7tra3jt3Z8x3LSFiUK4wUX+9nOe1XKWq5mAUUUVIBRRRQAUUUUAFFFFAEWmf8jFf/APXpb/8Aoc1bNY2mf8jFf/8AXpb/APoc1bNbx2KCiiiqAKKKKACiiigArO8Rf8ivqn/XnN/6Aa0azvEX/Ir6p/15zf8AoBoASs6TxBo8Wqrpkmq2SX7EBbVrhRIc9PlzmtBs7Tt644rgPh2+inwRbPqpsv7S+2P9v+0bfMF55x6553bsbe/TFch1noFVrDULXVLGO8sJlnt5M7JF6HBIP6g157pnia7uPHOifYrzUJtN1W5u4W+2yQbZBHHIwKRoNyAMgGTjI6jJqh4WvLmew8N6MdXm0ezbT7y8aaEorTSJcBQu5gRhQxYjvx2FPlFzHqVjf2up2a3VhMs8DFlDp0JUlT+RBFOe8tkvY7N7iJbmVGkjhLgO6rjJA6kDIyfcV494f1bVV0HwxpOnz3DQXcd/PJPazRQvcOlxgAM4IAw5YgcnjsDXQf2vr9tcaTDqdyn2o6JqcsrRlG3PG8QjYlRjcAeccZJo5Q5j0eivOLe9vYfCOgNfa1rF7qeuKkyrbNBFk+TvZQzABEA56liR35roPh5qt3rXgm1vNRlM1x51xCZCVJYRzOikleCcKMkcGiwXOk8O/wDILk/6+7n/ANHPWrWV4d/5Bcn/AF93P/o561a6Vscz3CiiimIKKKKACiiigAooooAKKKKACiiigArn7TjU9WX+7eD9Yoz/AFroK5+2/wCQvrH/AF+L/wCiIqiewmXKKKKxEFFFFABRRRQAUUUUAef+IPC+jz/EjQkktCVvo7yS5USuBKyiMqSAe2T+dWtQ8WXum6D4quLeKDOh3SW9srKcFfLib5ueT+8Pp2rqp9MtLnVLTUZoy1zZrIsL7iNocANx0Odo61kap4D0HWLy6uL6CdvthVriJLmRI5WUAK7IGALAAc47D0p3AzLrVta1PVvFFkrWcOl6UhiYGJmmm32wfruAXBbrg5BxgYycvwrrGs6bovhaw1Q2NxY6rpJWFII2V4DHAHG5ix3gqDnAXBxXdJotjHLqUiQkNqZBujvPz4QJ+HygDis7SvA+h6NPHLZQz/uYGt4EmuZJFt42xuWNWJC5wM49KLoDmNF8Q65e2OhaboSabYiXw/HqDNNFJIqEELsVd4OORyTxjvU1547vpvCum6nY3NrbXdzpn297L7BNdu3Gf+WZGxM5G5v6V1lh4Z0rTJLZ7O3MbW1iLCI+Yx2wAgheTzyBz1qhL4B0GS3toFhuIY7e0FkFgu5I/MgHIjchhuGSevqfWi6A2NHv/wC1dDsdQ2eX9rt459mc7dyg4/WrtQWVnDp9hb2douyC3jWKNck4VRgDJ9hU9IAooooAKKKKACiiigCLTP8AkYr/AP69Lf8A9DmrZrG0z/kYr/8A69Lf/wBDmrZreOxQUUUVQBRRRQAUUUUAFZ3iL/kV9U/685v/AEA1o1neIv8AkV9U/wCvOb/0A0AJWc3h/R31Uam+lWTX45F0bdTIPfdjNaNFch1mfBoGkWt817baXZxXTOZDOkChyxBBO4DOSCfzNFzoGj3lrDbXel2c0Fu2+KKSBWWNvUAjArQooAz7nQdJvLFbK70yzmtUcusDwKUViSSQMYByTz71LHpWnxRxJFY2yJDE0MarEoCI2NyjjgHAyO+BVuigCld6Npl9YxWV7p9rcWsOPLhlhVkTAwMAjAwOKntbS3srcQWcEdvCpJEcSBVBJyeB6kk1NRQAzw7/AMguT/r7uf8A0c9atZXh3/kFyf8AX3c/+jnrVrqWxyvcKKKKYgooooAKKKKACiiigAooooAKKKKACuftv+QvrP8A1+L/AOiIq6CsKaz1C11G9uLe2juoLmUS4WbbIpEaJjaRtP3M53Dr0qJJtCZZoqm2oeT/AMflneWx77oC6j/gSbl/Wn2+oWd2cW11DKw6qjgkfhWVmhFmiiikAUUUhOFJ9BQAtFUbR9ZvbGC6isbEJPGsihr18gEZGf3XvRdPrNnZzXMtjYlIY2kYLeuSQBk4/ddeKrlYWL1FIp3KCOhGaWpAKKKKACiqFnJrF9YwXcFjYiK4jWVA964YBhkZxEeeae02oW13bxX9rbRpcMUVoLlpCCFLcgovGAe9VysC5RRRUgFFFFABRRRQAUUUUARaZ/yMV/8A9elv/wChzVs1jaZ/yMV//wBelv8A+hzVs1vHYoKKKKoAooooAKKKKACs7xF/yK+qf9ec3/oBrRrO8Rf8ivqn/XnN/wCgGgBKKKK5DrCiiigAooooAKKKKAGeHf8AkFyf9fdz/wCjnrVrK8O/8guT/r7uf/Rz1q11LY5XuFFFFMQUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAVXutPs70YvLWGcdvNjDY/OrFFAGU3h6zXm1kubQ9vJnbav0Rsr+lRtpWpRf8AHvqUcy9cXVvlj7bkKgf98mtmilZMDnJbrUbW6it7jTPNeYMyfZJ1fhcZzv2eo6Zok1W3iRvtSz2vHJuIWRf++iNv61pXHPiWwHpa3Df+PRD+taNTyIVihoQK+HdNBGCLSIEHt8gqbUk8zSrtP70Lj81NWaRlDqVbkEYNWMxLF/M062f+9Ep/QVYrz7xd4ysvAvjPw14Zi1OQf2k/lyC5VZEtYsbYycBWOXwMlugJruWttZh/5Z2d2o6skjQsfopDD/x6seRk2LFFUZtSazUNqFheWyllUN5YlBJIAH7st1JAp0WrWEsgjW7iEn/PN22v/wB8nmps0Be8O/8AIr6V/wBecP8A6AKh1v8A4/tJ/wCvh/8A0U9TeHf+RX0r/rzh/wDQBUOt/wDH9pP/AF8P/wCinraWwySiiisBBRRRQAUUUUAFFFFAEWmf8jFf/wDXpb/+hzVs1jaZ/wAjFf8A/Xpb/wDoc1bNbx2KCiiiqAKKKKACiiigArO8Rf8AIr6p/wBec3/oBrRqlrUElzoGoQQLvlltZERR/ExUgCgCKiqf2+T/AKBt/wD9+f8A69H2+T/oG3//AH5/+vXLyvsdXMu5coqn9vk/6Bt//wB+f/r0jai6KWbTr4KBkkw9B+dHK+wcy7l2isux12LU7CC90+0u7i2uEEkUscWVdT0I5qx9vk/6Bt//AN+f/r0cr7BzLuXKKp/b5P8AoG3/AP35/wDr0fb5P+gbf/8Afn/69HK+wcy7lnw7/wAguT/r7uf/AEc9atZnh+OWPSj58Lws9xPIEkGGAaViM/gQa066Vscz3CiiimIKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDOn/wCRnsv+vO4/9DhrRrH1G+tLHxJYPe3UNsjWlwA00gQE74eOan/4SHRf+gvYf+BKf40AaNB4B4z7Vnf8JDov/QXsP/AlP8aP+Eh0X/oL2H/gSn+NAHzt43+B/jnxn4yl1/U9S0u1n1K5MdtbtM7GBFRmRSQpHCp2zya+g/Ci6xF4WsIfEyxjVIYhFctE+9ZGXjeD/tAA/jUdxqdhf6ppMdjfW1y63TOywzK5C+RKMkA9MkD8a26AM7XP+QfF/wBflr/6Pjq7Nbw3MZS4iSVD1V1DD9aoeIJEi0tJJWVEW7tizMcADz4+Sam/tnS/+glaf9/1/wAaALcaJFGscSqiIAqqowAB0AFVtQ0621COMXXmDyW3o8crRlTgjOVI7E+1N/tnS/8AoJWn/f8AX/GsHxrdS6t4QvtN8N6zp1tfXieQtxNcgLCrcM3GSTtzj3xQByvgDx1B458QeIdN0nVWV9LuStuLmJZVnt+F8wbdh+8G7ngrXctDrEH3re1ulH8UUpRj/wABYYH/AH1XiHw7+CXiDwf4vi8QeH/Een332G4EFxDtZVuImRWcBgSOj8Z7qK+i6nlQHPtftD/x+WF7b+5h8wfnHuFOg1KyupNlvdwvJ/cDjcPqOoreqG5s7W8j2XltDOv92VAw/Wp5EKxRoobw7YD/AI9hNaH+EW87oq/8Azt/MVG2kX8X/HtqvmDuLu3Vz+BTZj8jU8jCxJRVZk1eH/WWEM49ba4+Y/g4UfrVeHWoJIy8tvd24V3QmS3YqCrFT8y5XqD3qeVoRc0z/kYr/wD69Lf/ANDmrZrA0S7t7zX9Qe0njnQWtuC0bhgDvm44rfraOxQUUUVQBRRRQAUUUUAFFFFABRRRQAV5Z+0D42PhP4cy2VnJt1LWSbSAKfmVMfvG/wC+Tt+rCvU6geytZLxbuS2he5RSizNGC6qewPXHtQB4H+zD46a40+68F6m5E1qDc2O/q0RPzoPoTkezH0r6DrH0KxtJdB0W4ltYXngtIjFK0YLRkxgHaeoyDjitigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKK5K6gmt/HGmxaZqF9NNI8k9/FLcM8SQFTtGz7q/NtC4APB680Abk//ACM9l/153H/ocNaNZ0//ACM9l/153H/ocNaNABRRRQAUUVBfR+bYTIGnXKHm3OJP+AnsaAKmugNpsYIyDeWuQf8ArvHV77PD/wA8Y/8AvkVwul3MhutUtWkvooY72waGz1GVpJowZl3PliflYjjk9D0rvqAI/s8P/PGP/vkUfZ4f+eMf/fIqSigBFVUXCKFHoBiloooAKKKKACiiigArO0P/AJB8v/X5df8Ao+SuJ8ca/aT6ibZ9TjtV0u7tT5XnbGlkMqFiR1KqhPtlj6V2nh+RJdKaSJg6Pd3LKynIIM74NAGkFAJIABPXjrS0UUAFFFFABRRRQAUUUUAFFFFABRRRQAVxWn+KNSuvE0Vi00JnN3LFc6b9nIe2hUNtl355z8nbB38dK7WuPs/CN9DNYxTz2Zgsb1rtbpI2FzKSSdrHoM5wx5yB0oA3/D3/ACLGl/8AXnD/AOgCtGs7w9/yLGl/9ecP/oArRoAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACsOx8MLp+oz3dvquoZuJzPNGxiZZD6ElN2AOAM8DpW5SPu2NswGx8uRxmgDPn/5Gey/687j/ANDhrRr581T9pJNH8Tva6z4TuY9S04T2c0MVyCrOXTkErnHyencV7rolzfXmh2lzq1otneTRCSW2Vt3kk87M9yOh980AXqKKKACobu3N3ayQiaWAuMCWFsOvuDzU1FAHNXOhRacj3sl1cXt5cXVmj3FwV3bFuEwoCgADJJ4HeulrO1z/AJB8X/X5a/8Ao+OtGgAooooAKKKKACiiigAooooA5Tx94m8P+GdLtZfFqTpp1xcIPtEUZZYpFIdN235hkr6HpzXP+GPjF4O1K+tdF0a+m1LUL6+nEUNvbONqtK772LAAKFO498dq67xp4esPFfg3UtF1Z0jt7qEr5rkDyn6q/PowB/CvHP2bvh8uhXWra3qzQNqKyvY28aSBiiI2HcY7MRgH0HvQB9AUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAZ3h7/kWNL/684f8A0AVo1neHv+RY0v8A684f/QBWjQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUVm6h4i0jSdRs7HVNQgs7i+3fZlnbYJSuMgE8E8jjOaALN/fR6dZtczK7qGVAsYyzFmCgD8SKo/8JB/1C7//AL5j/wDi6d4k/wCQMv8A1923/o+OqN7qNrp7WwvJREbqdbeH5Sd8hBIXjpwD1qopMzlJp6HnfiD4e2Wt/GrS/GUmlXf2SGLzLuArHmWePAiON3TGM/7g9a9Q/wCEg/6hd/8A98x//F1k6h4k0rS2vFvrvyjY26XNz+7ZvLiYkBjgHjKt9MZOBVhdVsn1VdNS4V7toPtIjUE/u843EjgZPTJ55x0NVyoj2kjY07U49RMypDNC8LAOkoAPIyDwSKu1i6J/yFNT/wB6L/0GtqoejNou6uFFFFIZna5/yD4v+vy1/wDR8daNZ2uf8g+L/r8tf/R8daNABRRXP+IPF0Giaja6XbWF3quq3aNJFY2YXcI1ODI7MQqLkgZJ5PAzQB0FFZuj6ydU017q60+70p4nMcsF8qqykYOQQSrLzwwJFaCyo0YdXUo3RgeDQA6impLHKCYnV9pwdpzg+lYuo+JorHxLoukxxLP/AGq86GVZB+5MUe85HfPTtigDcqhLcT3c721gwjWM7ZbnAO0/3VB4Le54Hv0p+qXMlrp7tBjz3ZYosjIDuwUEjuATk+wqa2tI7WyS2i3bEXbnPzH1JPqeufWgCtFotgkglkt1uJxz59wPMf8AAnp9Bge1T3On2V4oW8tIJwOgljDY/OsU6DPZHKT3+oQ/3W1GZJl+h3hW+h2/U1PaWemXjMkVxqCzIMvDJf3CyJ9VL5x79D2oAtfY59P+fT3kmhH3rWV93H+wx5B9icduOtXbe4jurdJoW3Iw44wR6gjsQeCKpf2Haf8APW//APBjP/8AF1UtNOg0nxEBC9yY7yNnAlupXUSLjPDMQSykH/gBoA3KKKKACiiigAooooAKKKKACiiigAooooAzvD3/ACLGl/8AXnD/AOgCtGs7w9/yLGl/9ecP/oArRoAKKKKACiiigAooooAKKKKACiiigAooooAxrMalfxyzjURCouJo1RYFOAkrIOT7LVj7FqX/AEFv/JZaND/5B8v/AF+XX/o+StGgDO+xal/0Fv8AyWWj7FqX/QW/8llrRooAzvsWpf8AQW/8llo+xal/0Fv/ACWWtGigDO+xal/0Fv8AyWWvEv2lNL1LU9I8PaTaNNquo3F67wW0Fv8APtCYY4XnGWWvfqb5aeb5uxfM27d+OcemfSgDw34deCviF4Y8NxN4u17Fh9othHpDkTvH+/jx+8/gx/dUkV1vxAhmuI/D0drO1tK2twbZkQMU+STnB4P412XiT/kDj/r7tv8A0fHUdaR2Mam550Jf+Ef8XeJ7rxFeyapBHotsXLwIpkBknHlhVABJJwB3LUnwzsLnw3fXmi67Ht1KeCK6gkLl/wDRwoUQBj/zyPy++4HvXo1FVYzuJon/ACFNT/3ov/Qa2qxdE/5Cmp/70X/oNbVZS3OiPwoKKKKRRna5/wAg+L/r8tf/AEfHWjWdrn/IPi/6/LX/ANHx1o0AFedalqdp4T+NEuqeIp1s9O1XSIrW1vpzthjljldmiZzwpIcMM4Bwa9FpssUc0ZjmjWRD1VxkH8KAPMPF2v8Ah/xBrfhe5vby3u/CIu7mK8nl/wCPR7gRgw72PysuS2Dyu4dciuZ1KCyn0PxRF4fwPDE2vaUlkbVisO8yxC48kjou7HK8ZzivdGgiaHyWiQxYxsKjbj0xQIo1jWNY1CLjaoUYGOnFAHkPiPST4b1vxdZeCbT7AZ/CouFgsl2bplkkXeoH8e3jPU4FQaXF4JX4k+Bn8DC1J8m4+0tZ8jBtzt80jjzOG+983XNez7V3bto3YxnHOKZHbwxf6qKNOS3yqByepoAxBrOn69Z6Xf6Tcrc2bX5TzFBALKHQjkdnGK36p6hZfaNNeC2CRyAiSLjCiRWDKTjtuAz+NSWl2t3ZLcRowJHzRnG5WHVT7g8UAWKzNUk0+R0guImubofNFFB/rl9wQQUHbcSB71Qa+1q+YibSL+yt8/cikgMrD3bzML9Bk/7VW7W4ayjMdr4fvo1Jy2Hgyx9SfNyT7mgCxpUV/FC/9oyhwW/dKSGdF9GYAAn6D8T1pup/8hHR8dftbH8PIl/+tR/ad1/0BL//AL7g/wDjlVLG8m1TXzM+nXMEFrG0KSSNGVMhILfdc9NoGR3LDtQBuUUUUAFFFFABRRRQAUUUUAFFFFABRRRQBneHv+RY0v8A684f/QBWjWd4e/5FjS/+vOH/ANAFaNABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAGdof/ACD5f+vy6/8AR8laNZ2h/wDIPl/6/Lr/ANHyVo0AFFFFABRRRQAUUUUARXNtBeW7QXcEc8LY3RyoGU4ORkHjqM1S/wCEc0P/AKA2n/8AgKn+FXlniadoVlQyoAzRhhuUHoSKFnhed4UlRpUALxhgWUHoSO1AGDN4f0YeI7OMaRYhGtJ2K/ZkwSHhwcY9z+ZrQ/4RzQ/+gNp//gKn+FLP/wAjPZf9edx/6HDWjQBXtNPs9PRksLSC1VzuZYYggY+pwKsUUUAFFFFAGdrn/IPi/wCvy1/9Hx1o1na5/wAg+L/r8tf/AEfHWjQAUUUUAFFFFABRRRQAVRntZobhrrT9nmPjzYXOEl7Zzg7Wx3wc4wexF6igDO/tyzi4vy9g/dbpdg/B/ut+BNH/AAkOjt/qtStZm7JDKJGPsFXJJ+laNZ2h/wDIPl/6/Lr/ANHyUAHm3epfLBHJZ2x6zSDbK49FU8r9WwfbvV6GGO3hSGFQkaDCqOwp9FABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAGd4e/5FjS/wDrzh/9AFaNZ3h7/kWNL/684f8A0AVo0AFFFFABRRRQAUUUUAFFFFABRRRQBylgsun+L9fH2m4uttjBMBNJuwxabgDoBwBgDtWXpcRs7fwlq0NzNJearIq3heZmE4khdzlScDayjGBxjFdNbeGbS21qTVEub5riQYcPdMyMozhSvTA3HA7Umn+FdN028iuLfz2Fvu+zRSTM0dvu67FPA449hwKALOh/8g+X/r8uv/R8laNZ2h/8g+X/AK/Lr/0fJWjQAUUUUAFFFFABRRRQBy+l2lrZfETU0tESPzNPgkfByXcyy5JPUmk0m0tLH4harHZxpGJLC3d8Hl3Mk2ST1JrXi8O6NDqH26LTLVLvcX89YgH3Hqc/iaSPw/otnenUItNs4bhSXNwIlDDOcnd+JoAdP/yM9l/153H/AKHDWjWbKwbxLYspBU2VwQQevzw1oJIki7o2VxkjKnPNADqKKKACiiigDO1z/kHxf9flr/6PjrRrO1z/AJB8X/X5a/8Ao+OtGgAooooAKKKKACiiigAooooA42/jOka/dapq9hcTWcl7F5V1FeHEAKxoMxbhxvznr16V0Oh/8g+X/r8uv/R8lQ3Hhqzu74XF1NeTIJRMLaS5Ywhwcg7OnBGcdPaptD/5B8v/AF+XX/o+SgDRooooAKKKKACiiigAooooAKKKKACiiigAooooAzvD3/IsaX/15w/+gCtGs7w9/wAixpf/AF5w/wDoArRoAKKKKACiiigAooooAKKKKACiiigAooooAztD/wCQfL/1+XX/AKPkrRrO0P8A5B8v/X5df+j5K0aACiiigAooooAKK5DxNqepwa99msdQltIltkk2xxxtuZmcEkup7KKyv7V1z/oO3X/fm3/+N1rGlKSujmniYQlys9EpGVXQq6hlYYIIyCK88/tXXP8AoO3X/fm3/wDjdH9q65/0Hbr/AL82/wD8bqvYTI+uUvM+e/FvgXxLpHxnHgrQtQvobTUZg1htncIlu5yeAeiYIP8AuV9c6Fo1p4e0Gy0jTkK21nEsSZ6nHUn1JOST6mvPZUvJ9at9Xl1Odr+2ieGGcwQbkRyCwH7vvtH6+pq7/auuf9B26/782/8A8bo9hMPrlLzPRKK87/tXXP8AoO3X/fm3/wDjdH9q65/0Hbr/AL82/wD8bo9hMPrlLzPRKK4PSdW1c+JNMguNWnuIJ5nSSOSKEBgIZGHKoD1Ud67yspRcXZnRTqRqR5omdrn/ACD4v+vy1/8AR8daNZ2uf8g+L/r8tf8A0fHWjUmgUUUUAFFFFABRRRQAUUUUAFZ2h/8AIPl/6/Lr/wBHyVzvxV0fV9U8B3k3hm9urLWLAG6tZLWQozlR80fHXcuRj1xXiPwHufGPxA8XG513XtQn0TSH+0SwmUqk07MWVSBjPzZc59B60AfUNFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAGd4e/5FjS/+vOH/ANAFaNZ3h7/kWNL/AOvOH/0AVo0AFFFFABRRRQAUUUUAFFFczrnje20XX10ddK1TUrv7L9rZbCASbI9xXJywPUdBmgDpqK5r/hN7O50aw1LRNP1HWYr9GeJbK35ULw28uVCkHjBOcg8Vo+HfEFj4o0OHVdMMnkSll2ypsdGVirKw7EEEUAX7q6hs4DNcPtQcdMknsAO5PpWLeatqItpLhVhsYEGcSoZpG9sBgAfQZbOaf5n2/VJp25htXMMC9tw4d/rnK+2D6mq8n+n6oIutvZkO/o8uMqP+Ag7vqV9KpIlsp6LJrll5lpcXls0xZrkJJb8MHYswBVhjDMQeD2PfFdBZap58wt7yH7NckZVd25ZPXa3f6EA+1UNRgkeNLi2Gbm3O+MdN4/iT6EcfXB7U4iLU9PV4nIWQB45BwyN2I9CDTsFzdoqnpV419p0c0oCygtHKB0DqSrY9sg1cqCgooooA4fxT/wAjU3/XlF/6HLXIXvi7S9Pvrm2ujcL9kKC4mW3doodwBG5wMDgiuv8AFP8AyNTf9eUX/octeay6Bd614l8TW0l5JaaddNAkyLAC06+SoIVz0HY4Brtpt8iseRXUXWlzf1sdsDkZHIpa4tdIkl8Q+J7yCBze26IunO5O1HNuBlR064BPtWXpccT3GknwvBcRahHayjUmljdSW8o4EpYYZ/Nwe54PatOYw9mn1PSKK81hjt28B3KaBb3Sa59ijF8xilEhbcvmhicbn+/0OfSn29jaS6Tqxi1SGC0mSAMlnp8yQxsr5yyknO4YVgMcdaXMHs/M9HorzO8kim8PWNzJZ20BtnuBDYNbStbXoyPmQAZRjj5cg4ycA9al1X7LNd6ufEFldC6ltYv7KjVHdk/dfdjKjAcSZyeD07Ucw/ZHpul/8jVov/XzJ/6Ty16LXmHhcXK6n4ZXUDm7DYnPq/2WTd+ua9PrlrfEejg/4fzM7XP+QfF/1+Wv/o+OtGs7XP8AkHxf9flr/wCj460axOwKKKKACiiigAorH1rxb4f8OzRQ67rFnp8kyl40uJghYDgkZ7Val1vS4NGOry6japp3l+Z9rMy+Vt9d2cYoAvVHcXEVrA01w4jjXqT/AJ5PtSWt1BfWcN3ZyrNbzxrLFIhyrqwyCD6EHNcvql1fXniONYJoIrKFmgjMsJkzOBknhl7EqPQq3rTQGo+o6ldf8ecENpH2e6Bdz/wBSMfi2fasTwz4cl8H2dzb6INPEVzdSXcsQt3j3u5yRu3tgDgD5TgAVo+RrH/P/Zf+ATf/AB2jyNY/5/7L/wAAm/8AjtVZE3ZqWWqLcS/Z7mI21zjPlschx3Kt3H5H1Aq/XHatFq62Py31mbhnAt9tmwYS/wAJB8zt1PXgHIIzXTaVeNfabHNKAsvKSqvQOpKtj2yDj2qWrDTLdFFFIYUUUUAFFFFABRRRQAUUUUAFFFFAGd4e/wCRY0v/AK84f/QBWjWd4e/5FjS/+vOH/wBAFaNABRRRQAUUUUAFFFFABXnGsQa3cfG1k8P31nYynw8vmS3Vq0+F+0N91Q6859cj2r0emeTF5/neWnm7dnmbRu25zjPpQBxEmlWHhPwrYeGrbxl/YNwBJL9rl8jzLkklpGxICv3nzxyOKPg/cCXwAsSRReTa3txBFdRbtt6qyH/SBuJPzkknk85xxXZXdhZ38YS+tYLlFOQs0YcA/jUyRpFGqRKqIowqqMAD0AoA5e3ujaaRIVXfO11NHHHn78hmfA+nc+gyav2VqLO0WLdvblncjl2PJb8TWfqunXlr4iW7t7qCK0myUE0BkEc7YBOQ643ADHuWH8QFWPJ1j/n/ALH/AMAn/wDjtaIhmhWdF/oGpmE8W92xeL/Zk6sv4jLD3De1L5Osf8/9j/4BP/8AHao6rb6vNafZ0vrJriVh9nVbNw28HIb/AFvAB5J9KANvQOYr5h9xrx9n4AA/+PBq1qradZiw0+K23b2UZd8Y3sTlm/Ekn8as1mWFFFFAHD+Kf+Rqb/ryi/8AQ5azK6rXPDNxquqC8tb+K2/crEySWxkztZiCCHXH3qz/APhC9R/6DFr/AOADf/Ha66dSMYpM8yvh6k6jkkYtFbX/AAheo/8AQYtf/ABv/jtH/CF6j/0GLX/wAb/47V+2gY/VKvYxaK0X8LakmqQWf9rWpEsMku/7C3GxkGMeb33/AKVZ/wCEL1H/AKDFr/4AN/8AHaPbQD6pV7GLRW1/wheo/wDQYtf/AAAb/wCO0f8ACF6j/wBBi1/8AG/+O0e2gH1Sr2MvS/8AkatF/wCvmT/0nlr0WuV0/wAIXdprNle3GqQTJayM/lpaFCxMbJ94yHH389O1dUSAMngCuarJSldHoYanKnC0jxv4hfG638GeIp/D/iLQLwNHNBc21zbSKyzwiRXzhsYPylcc8iu4+Hfjd/iB4fk1yPSpNOsXmMdr50gZ5gvDOQBgDOQOT0NcR8cfBtj8RdHsH0eQy6pY3CqssNvJKrQuwEgLIpHy8N1/hI716R4XTR9L0Wy0PRpAkdlAsUcMilJNqjG4qwB56k46msjpNuiiigAooooA818UXN/a/GrS20vSBq0zaFOpgM6RADzk+YluPw680+HwTHpHwhvrPXre0uLqCK9vkjRd0VrJIJHCx57KGwDj34run0ewfXY9ZeDOoRW7WyTbm4jZgxXGcdQDnGai1zw/p/iOxFnqyTyQZJKQ3UsG7IIIJjZSRgng8UAUfAH/ACTXwz/2CLX/ANErUFrareaVcwOSpa7uSHXqjCdyGHuDg/hWpoPh3TfDOnix0eOaK2GNsctzLNsAAACmRmKgADgcVXnj/s/VZC3FveuHRuyy4wV9s4BHvu9qa3Exun3TXVrmYBZ42Mcyjorjrj2PBHsRVqs65/0HUY7scQzlYZ/Y5wjfmdp+o9KfqU0m2O0tmKz3JKhh1jQfef8AAHj3IqyRlt/p2oPdnmGAmKD3PR3/ADG0fQ+tXtB6agR937Y23/vlc/rmoJGSws44baLLYEUEKn7xxwo/Lk9gCa09Ns/sGnxwM29xlpHA+87Esx/MmpkNFqiiipKCiiigAooooAKKKKACiiigAooooAzvD3/IsaX/ANecP/oArRrO8Pf8ixpf/XnD/wCgCtGgAooooAKKKKACiiigAooooAKKKKAGyxRzxNFMiyRuMMrDIYehFZMmk3NsCbC/CxAZ8u7QyBR7NkMPxzWxXF/FY+IZ/Ad3png6xlu9U1IfZVZGCiGNh87liQB8uQOc5YelAEHgvxOfH2l3N7o99aRw211JbSbY2dsqeGGSAAy4YdetdfZabDZFnDPNO4w88pyzD09APYACvnr4H+EfHvw88WM2paQz6HqTNa3TQTpJ5MqMVWQqDnAYMpPo2e1fSNO4BRRRSAKKKKACiiigDFttaupvFk2kz2IghW2M8UrSAtIA+3OBwAeo5z9Khl1vVLLxFY2eoWVqtrqE0kUDQzlpEKqWBYFQMEL26ZqMQawfHI1A6bGLMWxtPM+1DdjzN2/bj07VDYWutt4pk1DVtMikG5obaRbobbaEnqFxks2ASfoOgoA2Z/8AkZ7L/rzuP/Q4a0azp/8AkZ7L/rzuP/Q4a0aACiiigArJ1352063fmG4uikqHo6iKRgD7ZUHHtWtWTrf/AB96R/1+t/6IloAnAAGAMAdhVTU4lk0+ZmGHiRnjcfeRgOGB7GrdV7//AJBtz/1xf+RoA0bOVprGCV/vPGrHHqRmpqrad/yC7X/rin/oIqzQAUUUUAFFFFABTJoY7iFop41kjcYZWGQRT6KAMO50S68iSC2uUntpFKmC7BJweMeYDnH1BPvWZodhrlwktxdTWQl3vbCY75CFicp935eSQSTnv7V19Z2h/wDIPl/6/Lr/ANHyU7sViSy0yO0czPI9xcsMGaTGQPRQOFHsPxzV2iikMKKKKACiiigAooooAKKKKACiiigAooooAzvD3/IsaX/15w/+gCtGs7w9/wAixpf/AF5w/wDoArRoAKKKKACiiigAooooAKKKKACiiigBk88VtC01zKkMSDLPIwVVHuTSW9xDdQLNazRzROMrJGwZW+hFY/ifTLnUYbCSzSGdrO8W4a2nbak4CsMZwcEFgw46rVLwX5pn15ngit4jqJ2RQvvRT5ab8HA/iznjrmgDY0P/AJB8v/X5df8Ao+StGs7Q/wDkHy/9fl1/6PkrRoAKKKKACiiigAooooAKKKKAM6+t73+1La8sYreXyoZYnSaZo/vNGQQQrf3D+dHn61/0D7D/AMDn/wDjNaNFAGYL/UIb21hvrK2jS5kMQeG6ZypCM/QxrxhD3rTrO1P/AJCGj/8AX43/AKIlrRoAKydb/wCPvSP+v1v/AERLWtWTrf8Ax96R/wBfrf8AoiWgCeq9/wD8g25/64v/ACNWKr3/APyDbn/ri/8AI0AXtO/5Bdr/ANcU/wDQRVmq2nf8gu1/64p/6CKs0AFFFFABRRRQAUUUUAFZ2h/8g+X/AK/Lr/0fJXIXGpZ1XVtV1exu7iy0+/ECzR3exbZFCciMMCxy2SfTjmuv0P8A5B8v/X5df+j5KANGiiigAooooAKKKKACiiigAooooAKKKKACiiigDO8Pf8ixpf8A15w/+gCtGs7w9/yLGl/9ecP/AKAK0aACiiigArznxDqWpzfEXU9PTxa3h6xsdHgu0PlwFDI0koYv5ikkYReARXo1chP4IttR+JV5r+safp99aNptvbWwuIlkeORJJGYgMMAYdeQe1AGfpVz4g8XfD7SdeufEEnhsvYme5NvbRnc3OJCZA2E2jdgAdetbXw71jUvEHw70bVdcQLfXVvvkKptDjJCvjtuUBse9ZnxA0TxJr0thYaTa2Nzoa/Pf2s141u1yQfljJVG/d9yO/Tp16nRm1BtLj/teytrG5XKmC1mMsaqOmGKr27YoAv0UUUAFFFFAFa/02z1S2+z6jbRXMQYMElXIBHf681JbWsFlbJb2cMcEMYwkcahVUewFS0UAZ2h/8g+X/r8uv/R8laNZ2h/8g+X/AK/Lr/0fJWjQAUUUUAFFFFABRRRQAUUUUAFHXrRWfc6vFFM0FrG93cLwyRYwn+8x4H06+xoA+W/ifc/EDwP8SE0DQ9f1RrG9mE2kR+aXK+ZlPLUnJBUsV+hFfTvhTR59B8L2On315NfXcUQNzdTSF2llPLtk84yTj0GBWRqugy6/q+lanfw2EVzpMzTWhMbTFGZSp5yuex6dQD2rU8vVe+qLn2tlx/OgDarJ1v8A4+9I/wCv1v8A0RLUfl6p/wBBT/yWWqd8l9FcWN1d3n2i3tpzJIBCF2Axuu7jnALDPtz2oA16r3//ACDbn/ri/wDI1NFLHPGskLrIjDKshyCPrVPVLqKO0kt9wa4mjZYoVPzOSMcD055PQd6ANTTv+QXa/wDXFP8A0EVZrAtrbVYLSGI6qMxoqnFuvYYqXy9U/wCgp/5LLQBtUVi7dXX7mpQt/wBdbXP8mFPW/wBTt/8Aj5tYbtO7WzbH/wC+G4/8eoA16Kr2d9b30Za2k3FTh0YFWQ+jKeQfrVigAooooAxbnwppt1fy3Mv2gLPIss9uszCKZ1xhmToT8o+uBnNWdD/5B8v/AF+XX/o+StGs7Q/+QfL/ANfl1/6PkoA0aKKKACiiigAooooAKKKKACiiigAooooAKKKKAM7w9/yLGl/9ecP/AKAK0azvD3/IsaX/ANecP/oArRoAKKKKACqt5qVpYbRdTBGf7kYBZ3+ijJP4Co9TvXtUjhtwDc3BKx55CgDlj7D9SQO9UEht7COSeVwGI3TXEpG5vcn09ug7U0ribLH/AAkVp1MF8E7P9il/lt3fpV2zv7W/RmtJ1k2nDqOGQ+jA8g+xrE+1XV9xYJ5EB/5eZl5b/cQ8/icD0BFURpUUXiVZIJJFu/shP2otlyd46+o/2enoBT5RXOxoqnpt6b23bzFCTwt5cyDoGHOR7EEEex9auVJQUUUUAFFFFAGdof8AyD5f+vy6/wDR8laNZ2h/8g+X/r8uv/R8laNABRRRQAUUUUAFFFFABRRVbUbz7Bps9yF3tGvyJnG9jwq/iSB+NAFS8upLq5eytXKJHj7RMp5GRnYp7HHJPYEdzkLFFFbQrHEqxxr0A4ApllbfZLRImbzH5aSTGN7k5ZvxJJrj/i5Gs3gFoZBmOXULNHXONym4jBH5GgDtwQehz9KK8513w7pvgfUfD+peFom0559Wt7G4gjlYx3EUxKsGQkjI+8D1GKoSfFuUyT39u1g9jBeNb/2d5Upu5Y1k2NIHHyA8FgmOg6g0Aeq0V55dePdSg8ZTaXM+l2AjvVt4bPUEliku4jt/exzk+WTycIATxjOTUHiz4j3/AIevNRdLrRSlhIANNy81zNH8uWLIcRHk4DKenJ5oA72XRdNmkaR7KESMcs6LtLfUjrU1rYWlkGFnbQwbvveWgXd9cda89174l3Nrr2q2enTabbrpO1TFexytJeSFA5VSnCDDAZO7ntitGw8Xa3r/AIkhtdFtbODTxp1pqM73gfzQsxfMYAIAbCdTwMdDngA7ignHXivNbD4j30niLSbS4utFu49SuzayWun+ZI9oSrsMzZKOflwQAPbpVLXte8QeIfC9pq6RWMOi3Gs2qRRfOLhY1vEUSFs7TuK/dwMA9T0oA9XorgLzxl4hbT9a1/TLXT20bRriaJreXf59ysJxKyuDtTkNtBBzjkjNaGleItb1vxpqFpZxWMejae0O+WQOZphLAsgCgHAILdT2IoA6ee18yRZ4G8m5QYSUDt/dYd19vyweau2F6LyFty+XNE2yWPOdrfXuCMEH0P4VFVSVvserW10vCTEW034/cP4Nx/wM0AbVFFFABWdof/IPl/6/Lr/0fJWjWdof/IPl/wCvy6/9HyUAaNFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAGDo2qw2ug2FvPb36SxW0aOv9nznDBQCPuVd/ty0/55X/8A4Lp//iK0aKAM7+3LT/nlf/8Agun/APiKP7ctP+eV/wD+C6f/AOIrRooA5C41xLnxbJb2kMr3C2cflLcxtAFDO+8/OASPlToD07davRadulWfUJPtUynK5XEcZ/2V/qcn3p2r2cLa1G91EskN3EIcsPuyIWZcehIZ+f8AZqLy76w/1Ja+tx/yzdgJV+jHhvo2D7npVrYl7mjWef8AkZU/69G/9DFWLW9gvFbyHyynDowKsh9Cp5FVz/yMqf8AXo3/AKGKYjN8QeLdO8DzXOq6uLlrZ7UFltoGlYsjYHTgZ39WIHHWvKdB+PerePvirpGj6fbtpGhtM0kqxIZridURnAO0HAJUZCjpnk17lpkQu9UvLplDQogtkyOGIJL/AFGSB9VNYEvwj8Lx+MrHxRo9p/ZOqWc3mFrQBY5gQQysnTkEjIweah7lLY6f+3LT/nlf/wDgun/+Io/ty0/55X//AILp/wD4itGikMzv7ctP+eV//wCC6f8A+Io/ty0/55X/AP4Lp/8A4itGigDA0jVoILKRJrfUEY3Vw4B06fo0zsD9zuCDV7+3LT/nlf8A/gun/wDiK0aKAM7+3LT/AJ5X/wD4Lp//AIij+3LT/nlf/wDgun/+IrRooAzv7ctP+eV//wCC6f8A+Io/ty0/55X/AP4Lp/8A4itGigDO/ty0/wCeV/8A+C6f/wCIo/ty0/55X/8A4Lp//iKk1KTU41Q6VBZynnzPtU7RgDtjajZ/Sq/hvVrjW9Ejv7m3jgMjuEEchdXQMQHBIBwcZHHQigCT+3LT/nlf/wDgun/+IrM1zW7WSG0i8q+2yXce7NhOPuneP4PVRXSVl+IBt01bj/n2njmJ9FDAOfwUsaAKD+JNOjYq5ulYdQbKYH/0CsfxMNK8X6FJpbXOpW4MsUwmtrGUujRuHU8xkdVHUVsaxpRvEElpHEJ93zFhGNw92aN/5VNpenLY2w3qhnYfOwVB+GVRcj8KAOVs9Fthq9pqOuavr+tS2TGS1jutPdI4nII37Y4Vy2CQCc4zxUK+HbSGaSCy1jxBaaTNdG7fTobCRV3l97Ksnl71QtyVB7nscV31FAHA6p4ds9XnuIbvVtdfSbq5W5m059OkdSwIbCuYi6qSoOAfXGKrX3g+wvbXV9PGr69b6Zq0z3E1pFpr8SOQSQ5iLbcgHbn26cV6PRQBwt7ocEuoXtzpes69pQ1JVF9Hbac5EpVdu9S0RKOVABK+g7ir9tY6Ra6jeXS/2s/2uwhsHR7Oc4SPzMHdsyWPmHJJ7CurooA86sfClnaR6NDNrOvXFroc6S2NudMZFQKpUBtsQL/KSMn19aWXwnp7qtpFquvQ6THfLfxacmnPsSQSiTG7yt2zdk7c8Z+leiUUAcBe+GrC7fULeLUtcttJ1SZpr3To9Ok2Ss2N+HMW5Q2PmAPc4xmtzTF0zStT1O9tk1HOoyRu8Z0+bbHsiWMBfk6YUV0dFAGf/bdr/wA8r7/wXz//ABFUda1q2/se4dYr0NEBKubCcDKEMOdnqK3qo6uPMsVthybmaOHHqCw3f+O7j+FAHjXxF/aYTSrm40nwbpsj3cRKPeX8TRqh/wBmI4Y/8Cx9DXqPgPxBJc+ANEutVfULy9urOO4nn+xSHc8i7zgquMDdgY4wBVrxf8O/C/ji2MfiHSop5cYS5QbJk+jjn8Dke1aPhfQV8MeGLHRIrqS6isY/JillADlAflBxwSBgZ74oAl/tu2/597//AMAJv/iaNDVxphLxvHvubiQLIpVtrTOwJB5GQQa0aKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAIbq1ivLZoJxlG9Dgg9iD2IPOayJDd6edt3DJcxD7txboXJH+0g5B9wCD146Vu0U07COVurvSrtlZpXW4ThHiRhKvsABn8CMH0pllYaxfaoZZT5NosJiFzIhjmcFgeI+x4+8cey11tFFwsRwQR21ukMCBI0GFUdhUlFFIYUUUUAFFFFABRRRQAUUUUAFFFFAGX4itr6+0Oez0whJrnELSFseWjHDsPcLnHvir9rbRWdpDbW6BIYUEaKOiqBgCpaKACmSxJPC8Uyh45FKspHDA8EU+igDFsWkhDWNyxae2wu9usqfwv+I6+4NW6lvrEXaq8b+TcR5McoGceoI7qe4/kcGs7+0BbyLDqai0lJwpZv3ch/wBl+h+hwfagC5RRRQAUUUUAFFFFABRRRQAUUVHPcQ2sJluZUhjHV5GCj8zQBJVazT7fqn2nrb2hZIz/AH5Twx/4CMr9S3pSRifVOIlkt7Q/emYFXkHog6gf7Rx7dcjWiiSCFIoUCRoAqqowAKAH0UUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRUVxcwWkXmXU0cCZxukcKM/U0AS0VUi1fTZ5Vig1C1kkY4VEmUk/QA1La3treo72VzDcLHI0TtFIHCupwynHQg9R2oAmooooAKKKKACimSzRQqpmkSMMwUF2AyScAfU0+gAooqrNqmn28xhuL62ilHVHmVWH4E0AWqKr3N/Z2Vus95dwW8LEKJJZAqknoMnjmrFABRRRQAUUUUAFFQJe2sl5JaR3MLXMSh5IFkBdFPQleoBqegArL8RosmjbJFDq1zbgqwyCPPStSq99Zpf2jW8jugLKwdMblZWDAjII4IHUUAZP9i2a8Qefbr/cguHjUf8AAQcfpUWlaTHdLdGe7vnMdyyKftcgwAB6H3q9/Ys//Qav/wDviD/41VvT7BdPgaNZZJmdzI8ku3cxP+6APyFAFX/hH7X/AJ+L7/wNl/8AiqP+Eftf+fi+/wDA2X/4qtSigDL/AOEftf8An4vv/A2X/wCKo/4R+1/5+L7/AMDZf/iq1KKAMv8A4R+1/wCfi+/8DZf/AIqj/hH7X/n4vv8AwNl/+KrUooA5i30xX1DUYmu74pDcKiL9rk4HlRt1znqx/Ok1DSrOCK3lWIyTC8tsSzO0jj9+nRmJIrWm0bzLuaeDULq1M7B5EiEZUsFC5+ZCeijv2qP+wi7xm41S9nSORJPLcRAMVYMM7YweoHegDWooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigArnPiIiyfDPxKrqGH9l3BwRn/lm1dHWV4n0ibxB4X1HSLe6jtGvrd7dpnhMuxXUqSFDLzz61FRNwaRcGlJNnGeCEhv9D8K2tz4VaARaZFdRahMIcGWNYwu0ozMM7ifmA47dRU9l4+isvBt/qzaLBaSx60+mi2t5fkeYzCPzGfaOCxyTtz9TXWeG9JuNC8MWGk3F1HdPZW6W6zJCYwwVQoJUs3PHrXn3iXQpPCXw91DTtR1iB4tY1kzvePppMFsJZPMcSrvb5Pl27sr1HI6jWclzt9P82v0MoRfIl1/4D/U09X+IOs6UfFMLaVYtJ4cs0u3l+0vtuFdSyqF2fKflYHLHGBwc8a58WX3/AAl2m6QlnbumoaTJfxsZGUo6FBsJweDv644x0rl/D1jN4l0/V9CfVNC1bStVs5EuNT0WB1ZJCAq73aR1c7ScDORt9K3dP8Cavba9pmr3viZLu60+wksVUaeI43VtuDgOSD8gzzz221Oq/ryf62K6af1t/wAEqaL4+8R6x4Xj19PDFstlOmI8ahlkfzvLYvlABGq5Yt1+XoM1t+D/ABcfEt1rNpJHAJNKuVhM9tKXinVkDq65Axwcd+nBNZEPw2vIvhvZeF/7eTzbG7W6hu1ssI22UyBZIjId65PTcOgrX8NeE7vQvEOtardax9uOrtFJJELURBJFQISCCeMDgdu5brTVrv8Art/wQfl/X9aGT8UbWCR/CVy8KNPF4js1SQqNygk5APoa73p1rnPF/hm68TLpS2uow2Q0/UIr795ambzGjJIXh1wDn3rohnaN2Ccc4HWktItef6L/ACB/Ffy/VmbpviXQtZuXt9H1rTr+eNdzxWt2krKM4yQpJAzXB+NEsh8bvCj39j9sjOmXoaNbbzi33MfKASep/OvRLXStOsZWlsbC1tpGGGeGFUJHoSBWDqfhK9v/AIh6T4mj1SCGPTYJYFtGsyxkEmNxL+YMH5Rj5fzo+0n6/k0Po/l+aPMrzw9qGj/CbWTqdm1laTeI4bnSrCYhms4GuY8LgEhM8naDxmu88XfEC68Ky30s2mwfZbNoQomudst4rkB2iUA8Ju5J6kEcda2fG3hmfxb4cOlW9/HYFp4ZjK9uZv8AVyK4GA69So5zXPa58Mr7WW8QoPEnkQa75LyobEO0UkYUDaxf/V/LnZ1Gfvdci7f1tFfoxef9btlzVfGeuQ+KtV0PSNEs7iSx09L9Jp75kV1YsMECMkN8hwOR6kVmWvjPWdf8WeDX0vyLbTdY0uW+e3lLFtwCZBI643nHHXk9sbSeDdSXxRqGsnWrd3vdMTTyj2J+XaWO/IkGeXPGB2571S0n4cXekr4VaHXl8/w/ayWbOLLAuIX25GC52N8g5+Ye1EdLX/r4v/tQfl/Xw/8ABHSfEG6tvEunaffabBAmoai9gsP2ndcQgBikrqBgK+zIGc4IOTXQ+JvESeHbG1cQ/aLq+u47O0gL7BJK54y2DgAAknB4HQ1ycPwv1CFbCJfFBMWn6u+p2+bBTIxcuWWRi/zt85AbAx3B4x0vjHwoni3SYLdb2TT7yzuo7yzu41DmGZD8rbTww5II75oXwq/f8NL/AKjfxO39PX/gHDRay/hj4n+LtV1mzt0ePSbLC2bEi4Z5ZFXquQSzBe/TPtXZaD4k1XUPEVxpupaJLBAtuJ4b+OOZYmOcGI+YiEMOvGQR6YxWNd/DG41q81a48R699qbVdPis5VtLMW4jaNy6SJl2IIY5xk8j04G/4Z0LXNL+bxD4mk1t0j8uLFotuoHGWYKTvfgc59cDk047Wfn+bJe7a8vyX/BHal4huF8R/wBgaNbwz6gtkb2RriQoiJu2qOASSzA+mACfY0LbxfqV9qFjoyaQllrc1g19dW11cBktUD7ANyA7yzZxjHAyfSl8R+DLzU/E9p4i8P66+i6pBbm0lc2y3EdxCW3bWQkcg8g5pJPBE0Ou6frumaw8erWts9pcT3cHnLeRs28h1VkwQ/I2kAdMYqVsr+f62/S435f1tf8AWxRtPiWL22sIE07yNWudWk0iWCSXMdvNGpZzvA+ZdoyvAJyOnNFt8RrgWmqvqOlRRS2OqrpESwXDyfaJ2IwcCPKrhgeNx6jGRzg+NNKs/DWnaXp8mr29hJf6rLqd3qmo2ha2afBOGIZfLOWXZ8w4Tqcc2tJ0Gfxh4XvtGuNQ0w2scsV3YazoMDxKtyHLFvnZhIwIBLBudxB5prVX/rpf9bX8gen9etv0v8y3qvxG1vRtF1m8u/DLt/ZrQmKZmlhiuUkYL8vmRg7lJ5U8HrmtK98U+ItOVV1PRtPsS5mZbiS/LQ7VC+WmdoYyOSRtA42k89Kr6r4D17XfB95pGs+LRd3V0Y1+0nTgkcao4fiJZB8xIGWLH2ArR1fwjqGpeJdO1q31tbWa1s5LSaP7IJEdXILNGC37t+Op3ccEGk720/rT/P8Aqw15lTRvH7+JLTRU0exjS/1TTm1Bo7iU7LeNWCckDLEscDgcAnjoc/xPJca/oGkp4l0H+zZZvNMzz6kY47GRPuEPGQWdiPkxg4J+hxtY8C3ngbwhperaV4hn/tbQIWtYpY7ONvtVu7gmExM43EdRhsk5wOgHWXPgq/udS0TVLbxDMt3p6TCZ7q1WXz/OxuYLlRGwxgEZAHGCKcrPbv8A52/T5iV1/X3mb4W8awyeFvC1tplrK2oa4ZxHHe3bSlPKLGWR5CNzcjjgZ3Acdui8NeKH1nVNY0fULRbTU9HlRJ0jk3xyI67kkUkA4IzwRwRjnrWBp3wuk0zSNFit9dY6joV1NNY3j2g2hJSd8TxhhuBB5IYHIBGOlbEfhK+tv7UvbHWltta1S4hkuL5bNWUJHgCNY2Y4G0EZJJyxPoKd7u/9dP6/pC2Wn9b/ANf0zqaKKKQwooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAP/2Q=="}
 },
 "cell_type": "markdown",
 "id": "857dd879",
 "metadata": {},
 "source": [
 "![trab1.JPG](attachment:trab1.JPG)"
 ]
 },
 {
 "cell_type": "markdown",
 "id": "990f1ea5",
 "metadata": {},
 "source": [
 " ## Tabela 1 - Dados operacionais da planta termoelétrica\n"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 87,
 "id": "dd559476",
 "metadata": {},
 "outputs": [
 {
 "data": {
 "text/html": [
 "<div>\n",
 "<style scoped>\n",
 " .dataframe tbody tr th:only-of-type {\n",
 " vertical-align: middle;\n",
 " }\n",
 "\n",
 " .dataframe tbody tr th {\n",
 " vertical-align: top;\n",
 " }\n",
 "\n",
 " .dataframe thead th {\n",
 " text-align: right;\n",
 " }\n",
 "</style>\n",
 "<table border=\"1\" class=\"dataframe\">\n",
 " <thead>\n",
 " <tr style=\"text-align: right;\">\n",
 " <th></th>\n",
 " <th>Temperatura(°C)</th>\n",
 " <th>Pressão(kPa)</th>\n",
 " <th>Vazão(kg/s)</th>\n",
 " <th>Título</th>\n",
 " </tr>\n",
 " </thead>\n",
 " <tbody>\n",
 " <tr>\n",
 " <th>1</th>\n",
 " <td>485</td>\n",
 " <td>6495</td>\n",
 " <td>27.9</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>2</th>\n",
 " <td>-</td>\n",
 " <td>900</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>3</th>\n",
 " <td>-</td>\n",
 " <td>250</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>4</th>\n",
 " <td>51</td>\n",
 " <td>Saturação</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>5</th>\n",
 " <td>51</td>\n",
 " <td>Saturação</td>\n",
 " <td>-</td>\n",
 " <td>0</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>6</th>\n",
 " <td>-</td>\n",
 " <td>250</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>7</th>\n",
 " <td>110</td>\n",
 " <td>250</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>8</th>\n",
 " <td>-</td>\n",
 " <td>6495</td>\n",
 " <td>27.9</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>9</th>\n",
 " <td>110</td>\n",
 " <td>250</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>10</th>\n",
 " <td>-</td>\n",
 " <td>900</td>\n",
 " <td>-</td>\n",
 " <td>-</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>11</th>\n",
 " <td>Saturação</td>\n",
 " <td>900</td>\n",
 " <td>≤ 6.94</td>\n",
 " <td>1</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>12</th>\n",
 " <td>Saturação</td>\n",
 " <td>900</td>\n",
 " <td>≤ 6.94</td>\n",
 " <td>0</td>\n",
 " </tr>\n",
 " </tbody>\n",
 "</table>\n",
 "</div>"
 ],
 "text/plain": [
 " Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título\n",
 "1 485 6495 27.9 -\n",
 "2 - 900 - -\n",
 "3 - 250 - -\n",
 "4 51 Saturação - -\n",
 "5 51 Saturação - 0\n",
 "6 - 250 - -\n",
 "7 110 250 - -\n",
 "8 - 6495 27.9 -\n",
 "9 110 250 - -\n",
 "10 - 900 - -\n",
 "11 Saturação 900 ≤ 6.94 1\n",
 "12 Saturação 900 ≤ 6.94 0"
 ]
 },
 "execution_count": 87,
 "metadata": {},
 "output_type": "execute_result"
 }
 ],
 "source": [
 "df"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 123,
 "id": "f54ee48f",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 1 - Entrada da Turbina ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 1 -> T1 = 485.00 °C\n",
 "\n",
 "Pressão no ponto 1 -> P1 = 6495.00 kPa\n",
 "\n",
 "Entalpia no ponto 1 -> h1 = 3381.20 Kj/kg\n",
 "\n",
 "Entropia no ponto 1 -> s1 = 6.79 Kj/kg\n",
 "\n",
 "Densidade no ponto 1 -> ρ1 = 19.65 Kg/m^3\n",
 "\n",
 "Título no ponto 1 -> X1 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "fl = 'Water'\n",
 "\n",
 "#Estado 1 - Entrada da Turbina\n",
 "P1 = 6495 #kPa\n",
 "m1 = 27.9 #kg/s\n",
 "T1 = 485+273.15 #K\n",
 "\n",
 "st1 = st(fl,{'P':P1,'T':T1})\n",
 "h1 = st1.h\n",
 "s1 = st1.s\n",
 "ρ1 = st1.rho\n",
 "tit1 = st1.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 1 - Entrada da Turbina','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 1 -> T1 = {:.2f}'.format(T1-273.15),'°C')\n",
 "print ('\\nPressão no ponto 1 -> P1 = {:.2f}'.format(P1),'kPa')\n",
 "print ('\\nEntalpia no ponto 1 -> h1 = {:.2f}'.format(st1.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 1 -> s1 = {:.2f}'.format(st1.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 1 -> ρ1 = {:.2f}'.format(ρ1),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 1 -> X1 = {:.2f}\\n\\n'.format(tit1))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 105,
 "id": "c4dc29ac",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 2 - Primeira Extração da Turbina ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 2 -> T2 = 516.42 K\n",
 "\n",
 "Entalpia no ponto 2 -> h2 = 2932.08 Kj/kg\n",
 "\n",
 "Entropia no ponto 2 -> s2 = 6.95 Kj/kg\n",
 "\n",
 "Densidade no ponto 2 -> ρ2 = 3.91 Kg/m^3\n",
 "\n",
 "Título no ponto 2 -> X2 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 2 - Primeira Extração da Turbina\n",
 "\n",
 "P2 = 900 #kPA\n",
 "\n",
 " #Estado ideal\n",
 "st2s = st(fl,{'P':P2,'S':st1.s})\n",
 "\n",
 " #Eficiência da Turbina\n",
 "ηTurb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h2 = st1.h - (st1.h - st2s.h)*ηTurb\n",
 "st2 = st(fl,{'P':P2,'H':h2})\n",
 "T2 = st2.T\n",
 "s2 = st2.s\n",
 "ρ2 = st2.rho\n",
 "tit2 = st2.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 2 - Primeira Extração da Turbina','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 2 -> T2 = {:.2f}'.format(st2.T),'K')\n",
 "print ('\\nEntalpia no ponto 2 -> h2 = {:.2f}'.format(st2.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 2 -> s2 = {:.2f}'.format(st2.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 2 -> ρ2 = {:.2f}'.format(ρ2),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 2 -> X2 = {:.2f}\\n\\n'.format(tit2))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 125,
 "id": "a7bf6095",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 3 - Segunda Extração da Turbina ----------\n","\n",
 "\n",
 "Temperatura no ponto 3 -> T3 = 127.41 °C\n",
 "\n",
 "Entalpia no ponto 3 -> h3 = 2714.68 Kj/kg\n",
 "\n",
 "Entropia no ponto 3 -> s3 = 7.05 Kj/kg\n",
 "\n",
 "Densidade no ponto 3 -> ρ3 = 1.39 Kg/m^3\n",
 "\n",
 "Título no ponto 3 -> X3 = 1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 3 - Segunda Extração da Turbina\n",
 "\n",
 "P3 = 250 #kPa\n",
 "\n",
 " #Estado ideal\n",
 "st3s = st(fl,{'P':P3, 'S':st2.s})\n",
 "\n",
 " #Eficiência da Turbina\n",
 "ηTurb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h3 = st2.h - (st2.h - st3s.h)*ηTurb\n",
 "st3 = st(fl,{'P':P3,'H':h3})\n",
 "T3 = st3.T\n",
 "s3 = st3.s\n",
 "ρ3 = st3.rho\n",
 "tit3 = st3.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 3 - Segunda Extração da Turbina','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 3 -> T3 = {:.2f}'.format(T3-273.15),'°C')\n",
 "print ('\\nEntalpia no ponto 3 -> h3 = {:.2f}'.format(st3.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 3 -> s3 = {:.2f}'.format(st3.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 3 -> ρ3 = {:.2f}'.format(ρ3),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 3 -> X3 = {:.2f}\\n\\n'.format(tit3))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 107,
 "id": "920a3a89",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 4 - Última Extração da Turbina ----------\n",
 "\n",
 "\n",
 "Pressão no ponto 4 -> P4 = 12.98 kPa\n",
 "\n",
 "Entalpia no ponto 4 -> h4 = 2333.12 Kj/kg\n",
 "\n",
 "Entropia no ponto 4 -> s4 = 7.26 Kj/kg\n",
 "\n",
 "Densidade no ponto 4 -> ρ4 = 0.10 Kg/m^3\n",
 "\n",
 "Título no ponto 4 -> X4 = 0.89\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "# Estado 4 - Última Extração da Turbina\n",
 "\n",
 "T4 = 51 + 273.15 #K\n",
 "\n",
 " #Estado ideal\n",
 "st4s = st(fl,{'T':T4, 'S':st3.s})\n",
 "\n",
 " #Eficiência da Turbina\n",
 "ηTurb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h4 = st3.h - (st3.h - st4s.h)*ηTurb\n",
 "P4 = ps('P','T',T4,'Q',1 , fl)/1000\n",
 "st4 = st(fl,{'P':P4,'H':h4})\n",
 "s4 = st4.s\n",
 "ρ4 = st4.rho\n",
 "tit4 = st4.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 4 - Última Extração da Turbina','-'*10)\n",
 "print ('\\n\\nPressão no ponto 4 -> P4 = {:.2f}'.format(P4),'kPa')\n",
 "print ('\\nEntalpia no ponto 4 -> h4 = {:.2f}'.format(st4.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 4 -> s4 = {:.2f}'.format(st4.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 4 -> ρ4 = {:.2f}'.format(ρ4),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 4 -> X4 = {:.2f}\\n\\n'.format(tit4))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 126,
 "id": "b1c33dac",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 5 - Saída do Condensador ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 5 -> T5 = 51.00 °C\n",
 "\n",
 "Pressão no ponto 5 -> P5 = 12.98 kPa\n",
 "\n",
 "Entalpia no ponto 5 -> h5 = 213.52 Kj/kg\n",
 "\n",
 "Entropia no ponto 5 -> s5 = 0.72 Kj/kg\n",
 "\n",
 "Densidade no ponto 5 -> ρ5 = 987.54 Kg/m^3\n",
 "\n",
 "Título no ponto 5 -> X5 = 0.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 5 - Saída do Condensador\n",
 "P5 = P4 \n",
 "st5 = st(fl,{'P':P5,'Q':0})\n",
 "T5 = st5.T\n",
 "s5 = st5.s\n",
 "ρ5 = st5.rho\n",
 "tit5 = st5.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 5 - Saída do Condensador','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 5 -> T5 = {:.2f}'.format(T5-273.15),'°C')\n",
 "print ('\\nPressão no ponto 5 -> P5 = {:.2f}'.format(P5),'kPa')\n",
 "print ('\\nEntalpia no ponto 5 -> h5 = {:.2f}'.format(st5.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 5 -> s5 = {:.2f}'.format(st5.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 5 -> ρ5 = {:.2f}'.format(ρ5),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 5 -> X5 = {:.2f}\\n\\n'.format(tit5))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 128,
 "id": "b4cf3ec1",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 6- Saída da Bomba e Entrada do Desareador ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 6 -> T6 = 51.02 °C\n",
 "\n",
 "Pressão no ponto 6 -> P6 = 250.00 kPa\n",
 "\n",
 "Entalpia no ponto 6 -> h6 = 213.81 Kj/kg\n",
 "\n",
 "Entropia no ponto 6 -> s6 = 0.72 Kj/kg\n",
 "\n",
 "Densidade no ponto 6 -> ρ6 = 987.64 Kg/m^3\n",
 "\n",
 "Título no ponto 6 -> X6 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "# Estado 6 - Saída da Bomba e Entrada do Desareador\n",
 "\n",
 "P6 = P3\n",
 "\n",
 " #Estado ideal\n",
 "st6s = st(fl,{'P':P6, 'S':st5.s})\n",
 "\n",
 " #Eficiência da Bomba\n",
 "ηBomb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h6 = ((st6s.h - st5.h)/ηBomb)+st5.h\n",
 "st6 = st(fl,{'P':P6,'H':h6})\n",
 "T6 = st6.T\n",
 "s6 = st6.s\n",
 "ρ6 = st6.rho\n",
 "tit6 = st6.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 6- Saída da Bomba e Entrada do Desareador','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 6 -> T6 = {:.2f}'.format(T6-273.15),'°C')\n",
 "print ('\\nPressão no ponto 6 -> P6 = {:.2f}'.format(P6),'kPa')\n",
 "print ('\\nEntalpia no ponto 6 -> h6 = {:.2f}'.format(st6.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 6 -> s6 = {:.2f}'.format(st6.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 6 -> ρ6 = {:.2f}'.format(ρ6),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 6 -> X6 = {:.2f}\\n\\n'.format(tit6))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 129,
 "id": "fcf294bc",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 7 - Saída do Desaerador ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 7 -> T7 = 110.00 °C\n",
 "\n",
 "Pressão no ponto 7 -> P7 = 250.00 kPa\n",
 "\n",
 "Entalpia no ponto 7 -> h7 = 461.49 Kj/kg\n",
 "\n",
 "Entropia no ponto 7 -> s7 = 1.42 Kj/kg\n",
 "\n",
 "Densidade no ponto 7 -> ρ7 = 951.00 Kg/m^3\n",
 "\n",
 "Título no ponto 7 -> X7 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 7 - Saída do Desaerador \n",
 "\n",
 "P7 = P3 \n",
 "T7 = 110 + 273.15 #kPa\n",
 "st7 = st(fl,{'P':P7,'T':T7})\n",
 "h7 = st7.h\n",
 "s7 = st7.s\n",
 "ρ7 = st7.rho\n",
 "tit7 = st7.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 7 - Saída do Desaerador ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 7 -> T7 = {:.2f}'.format(T7-273.15),'°C')\n",
 "print ('\\nPressão no ponto 7 -> P7 = {:.2f}'.format(P7),'kPa')\n",
 "print ('\\nEntalpia no ponto 7 -> h7 = {:.2f}'.format(st7.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 7 -> s7 = {:.2f}'.format(st7.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 7 -> ρ7 = {:.2f}'.format(ρ7),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 7 -> X7 = {:.2f}\\n\\n'.format(tit7))\n"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 130,
 "id": "6ece0c74",
 "metadata": {
 "scrolled": true
 },
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 8 - Entrada da Caldeira ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 8 -> T8 = 110.75 °C\n",
 "\n",
 "Pressão no ponto 8 -> P8 = 6495.00 kPa\n",
 "\n","Entalpia no ponto 8 -> h8 = 469.21 Kj/kg\n",
 "\n",
 "Entropia no ponto 8 -> s8 = 1.42 Kj/kg\n",
 "\n",
 "Densidade no ponto 8 -> ρ8 = 953.43 Kg/m^3\n",
 "\n",
 "Título no ponto 8 -> X8 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 8 - Entrada da Caldeira \n",
 "P8 = P1 \n",
 "\n",
 " #Estado ideal\n",
 "st8s = st(fl,{'P': P8, 'S': st7.s})\n",
 "\n",
 " #Eficiência da Bomba\n",
 "ηBomb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h8 = ((st8s.h-h7)/ηBomb)+h7\n",
 "st8 =st(fl,{'P': P8, 'H': h8})\n",
 "T8 = st8.T\n",
 "s8 = st8.s\n",
 "ρ8 = st8.rho\n",
 "tit8 = st8.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 8 - Entrada da Caldeira ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 8 -> T8 = {:.2f}'.format(T8-273.15),'°C')\n",
 "print ('\\nPressão no ponto 8 -> P8 = {:.2f}'.format(P8),'kPa')\n",
 "print ('\\nEntalpia no ponto 8 -> h8 = {:.2f}'.format(st8.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 8 -> s8 = {:.2f}'.format(st8.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 8 -> ρ8 = {:.2f}'.format(ρ8),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 8 -> X8 = {:.2f}\\n\\n'.format(tit8))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 131,
 "id": "04a9b58d",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 9 - Entrada da bomba para HX ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 9 -> T9 = 110.00 °C\n",
 "\n",
 "Pressão no ponto 9 -> P9 = 250.00 kPa\n",
 "\n",
 "Entalpia no ponto 9 -> h9 = 461.49 Kj/kg\n",
 "\n",
 "Entropia no ponto 9 -> s9 = 1.42 Kj/kg\n",
 "\n",
 "Densidade no ponto 9 -> ρ9 = 951.00 Kg/m^3\n",
 "\n",
 "Título no ponto 9 -> X9 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 9 - Entrada da bomba para HX\n",
 "\n",
 "st9 = st7\n",
 "P9 = P7\n",
 "T9 = st9.T\n",
 "h9 = st9.h\n",
 "s9 = st9.s\n",
 "ρ9 = st9.rho\n",
 "tit9 = st9.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 9 - Entrada da bomba para HX ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 9 -> T9 = {:.2f}'.format(T9-273.15),'°C')\n",
 "print ('\\nPressão no ponto 9 -> P9 = {:.2f}'.format(P9),'kPa')\n",
 "print ('\\nEntalpia no ponto 9 -> h9 = {:.2f}'.format(st9.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 9 -> s9 = {:.2f}'.format(st9.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 9 -> ρ9 = {:.2f}'.format(ρ9),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 9 -> X9 = {:.2f}\\n\\n'.format(tit9))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 132,
 "id": "d2e64558",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 10 - Saída da Turbina ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 10 -> T10 = 110.08 °C\n",
 "\n",
 "Pressão no ponto 10 -> P10 = 900.00 kPa\n",
 "\n",
 "Entalpia no ponto 10 -> h10 = 462.30 Kj/kg\n",
 "\n",
 "Entropia no ponto 10 -> s10 = 1.42 Kj/kg\n",
 "\n",
 "Densidade no ponto 10 -> ρ10 = 951.25 Kg/m^3\n",
 "\n",
 "Título no ponto 10 -> X10 = -1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 10 - Saída da Turbina \n",
 "P10 = P2\n",
 "\n",
 " #Estado ideal\n",
 "st10s = st(fl,{'P': P10, 'S': st9.s})\n",
 "\n",
 " #Eficiência da Bomba\n",
 "ηBomb = 0.85\n",
 "\n",
 " #Estado real\n",
 "h10 = ((st10s.h-h9)/ηBomb)+h9\n",
 "st10 =st(fl,{'P': P10, 'H': h10})\n",
 "T10 = st10.T\n",
 "s10 = st10.s\n",
 "ρ10 = st10.rho\n",
 "tit10 = st10.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 10 - Saída da Turbina ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 10 -> T10 = {:.2f}'.format(T10-273.15),'°C')\n",
 "print ('\\nPressão no ponto 10 -> P10 = {:.2f}'.format(P10),'kPa')\n",
 "print ('\\nEntalpia no ponto 10 -> h10 = {:.2f}'.format(st10.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 10 -> s10 = {:.2f}'.format(st10.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 10 -> ρ10 = {:.2f}'.format(ρ10),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 10 -> X10 = {:.2f}\\n\\n'.format(tit10))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 133,
 "id": "0bb9078e",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 11 - Entrada do Processo Industrial Vizinho ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 11 -> T11 = 175.35 °C\n",
 "\n",
 "Pressão no ponto 11 -> P11 = 900.00 kPa\n",
 "\n",
 "Entalpia no ponto 11 -> h11 = 2773.03 Kj/kg\n",
 "\n",
 "Entropia no ponto 11 -> s11 = 6.62 Kj/kg\n",
 "\n",
 "Densidade no ponto 11 -> ρ11 = 4.65 Kg/m^3\n",
 "\n",
 "Título no ponto 11 -> X11 = 1.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 11 - Entrada do Processo Industrial Vizinho\n",
 "Q11 = 1\n",
 "P11 = P10\n",
 "st11 = st(fl,{'P': P11, 'Q': Q11})\n",
 "T11 = st11.T\n",
 "s11 = st11.s\n",
 "ρ11 = st11.rho\n",
 "tit11 = st11.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 11 - Entrada do Processo Industrial Vizinho ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 11 -> T11 = {:.2f}'.format(T11-273.15),'°C')\n",
 "print ('\\nPressão no ponto 11 -> P11 = {:.2f}'.format(P11),'kPa')\n",
 "print ('\\nEntalpia no ponto 11 -> h11 = {:.2f}'.format(st11.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 11 -> s11 = {:.2f}'.format(st11.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 11 -> ρ11 = {:.2f}'.format(ρ11),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 11 -> X11 = {:.2f}\\n\\n'.format(tit11))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 134,
 "id": "3bb6bdde",
 "metadata": {},
 "outputs": [
 {
 "name": "stdout",
 "output_type": "stream",
 "text": [
 "\n",
 "\n",
 " ---------- ESTADO 12 - Saída do Processo Industrial Vizinho ----------\n",
 "\n",
 "\n",
 "Temperatura no ponto 12 -> T12 = 175.35 °C\n",
 "\n",
 "Pressão no ponto 12 -> P12 = 900.00 kPa\n",
 "\n",
 "Entalpia no ponto 12 -> h12 = 742.56 Kj/kg\n",
 "\n",
 "Entropia no ponto 12 -> s12 = 2.09 Kj/kg\n",
 "\n",
 "Densidade no ponto 12 -> ρ12 = 891.92 Kg/m^3\n",
 "\n",
 "Título no ponto 12 -> X12 = 0.00\n",
 "\n",
 "\n"
 ]
 }
 ],
 "source": [
 "#Estado 12 - Saída do Processo Industrial Vizinho\n",
 "Q12 = 0 \n",
 "P12 = P10\n",
 "st12 = st(fl,{'P': P12, 'Q': Q12})\n",
 "T12 = st12.T\n",
 "s12 = st12.s\n",
 "ρ12 = st12.rho\n",
 "tit12 = st12.Q\n",
 "\n",
 "print('\\n\\n','-'*10,'ESTADO 12 - Saída do Processo Industrial Vizinho ','-'*10)\n",
 "print ('\\n\\nTemperatura no ponto 12 -> T12 = {:.2f}'.format(T12-273.15),'°C')\n",
 "print ('\\nPressão no ponto 12 -> P12 = {:.2f}'.format(P12),'kPa')\n",
 "print ('\\nEntalpia no ponto 12 -> h12 = {:.2f}'.format(st12.h),'Kj/kg')\n",
 "print ('\\nEntropia no ponto 12 -> s12 = {:.2f}'.format(st12.s),'Kj/kg')\n",
 "print ('\\nDensidade no ponto 12 -> ρ12 = {:.2f}'.format(ρ12),'Kg/m^3')\n",
 "print('\\nTítulo no ponto 12 -> X12 = {:.2f}\\n\\n'.format(tit12))"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 142,
 "id": "11c18191",
 "metadata": {},
 "outputs": [],
 "source": [
 "i1 = [[485,6495, 27.9, tit1],\n",
 " [round(T2-273.15,2),900,'-',tit2],[round(T3-273.15,2),250,'-',tit3],[51,round(P4,2),'-',tit4],\n",
 " [51,round(P5,2),'-',0],[round(T6-273.15,2),250,'-',tit6],\n",
 " [110,250,'-',tit7],[round(T8-273.15,2),6495,27.9,tit8],\n",
 " [110,250,'-',tit9],[round(T10-273.15,2),900,'-',tit10],\n"," [round(T11-273.15,2),900,'≤ 6.94',1],[round(T12-273.15,2),900,'≤ 6.94',0]]\n",
 "df1=pd.DataFrame(i1,index='1 2 3 4 5 6 7 8 9 10 11 12'.split(),\n",
 " columns='Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título '.split())\n"
 ]
 },
 {
 "cell_type": "markdown",
 "id": "77ee8c5c",
 "metadata": {},
 "source": [
 " ## Tabela 1 Atualizada - Dados operacionais da planta termoelétrica\n"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": 143,
 "id": "de6603ba",
 "metadata": {},
 "outputs": [
 {
 "data": {
 "text/html": [
 "<div>\n",
 "<style scoped>\n",
 " .dataframe tbody tr th:only-of-type {\n",
 " vertical-align: middle;\n",
 " }\n",
 "\n",
 " .dataframe tbody tr th {\n",
 " vertical-align: top;\n",
 " }\n",
 "\n",
 " .dataframe thead th {\n",
 " text-align: right;\n",
 " }\n",
 "</style>\n",
 "<table border=\"1\" class=\"dataframe\">\n",
 " <thead>\n",
 " <tr style=\"text-align: right;\">\n",
 " <th></th>\n",
 " <th>Temperatura(°C)</th>\n",
 " <th>Pressão(kPa)</th>\n",
 " <th>Vazão(kg/s)</th>\n",
 " <th>Título</th>\n",
 " </tr>\n",
 " </thead>\n",
 " <tbody>\n",
 " <tr>\n",
 " <th>1</th>\n",
 " <td>485.00</td>\n",
 " <td>6495.00</td>\n",
 " <td>27.9</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>2</th>\n",
 " <td>243.27</td>\n",
 " <td>900.00</td>\n",
 " <td>-</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>3</th>\n",
 " <td>127.41</td>\n",
 " <td>250.00</td>\n",
 " <td>-</td>\n",
 " <td>0.999172</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>4</th>\n",
 " <td>51.00</td>\n",
 " <td>12.98</td>\n",
 " <td>-</td>\n",
 " <td>0.890762</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>5</th>\n",
 " <td>51.00</td>\n",
 " <td>12.98</td>\n",
 " <td>-</td>\n",
 " <td>0.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>6</th>\n",
 " <td>51.02</td>\n",
 " <td>250.00</td>\n",
 " <td>-</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>7</th>\n",
 " <td>110.00</td>\n",
 " <td>250.00</td>\n",
 " <td>-</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>8</th>\n",
 " <td>110.75</td>\n",
 " <td>6495.00</td>\n",
 " <td>27.9</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>9</th>\n",
 " <td>110.00</td>\n",
 " <td>250.00</td>\n",
 " <td>-</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>10</th>\n",
 " <td>110.08</td>\n",
 " <td>900.00</td>\n",
 " <td>-</td>\n",
 " <td>-1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>11</th>\n",
 " <td>175.35</td>\n",
 " <td>900.00</td>\n",
 " <td>≤ 6.94</td>\n",
 " <td>1.000000</td>\n",
 " </tr>\n",
 " <tr>\n",
 " <th>12</th>\n",
 " <td>175.35</td>\n",
 " <td>900.00</td>\n",
 " <td>≤ 6.94</td>\n",
 " <td>0.000000</td>\n",
 " </tr>\n",
 " </tbody>\n",
 "</table>\n",
 "</div>"
 ],
 "text/plain": [
 " Temperatura(°C) Pressão(kPa) Vazão(kg/s) Título\n",
 "1 485.00 6495.00 27.9 -1.000000\n",
 "2 243.27 900.00 - -1.000000\n",
 "3 127.41 250.00 - 0.999172\n",
 "4 51.00 12.98 - 0.890762\n",
 "5 51.00 12.98 - 0.000000\n",
 "6 51.02 250.00 - -1.000000\n",
 "7 110.00 250.00 - -1.000000\n",
 "8 110.75 6495.00 27.9 -1.000000\n",
 "9 110.00 250.00 - -1.000000\n",
 "10 110.08 900.00 - -1.000000\n",
 "11 175.35 900.00 ≤ 6.94 1.000000\n",
 "12 175.35 900.00 ≤ 6.94 0.000000"
 ]
 },
 "execution_count": 143,
 "metadata": {},
 "output_type": "execute_result"
 }
 ],
 "source": [
 "df1"
 ]
 },
 {
 "cell_type": "code",
 "execution_count": null,
 "id": "b7107280",
 "metadata": {},
 "outputs": [],
 "source": []
 }
 ],
 "metadata": {
 "kernelspec": {
 "display_name": "Python 3",
 "language": "python",
 "name": "python3"
 },
 "language_info": {
 "codemirror_mode": {
 "name": "ipython",
 "version": 3
 },
 "file_extension": ".py",
 "mimetype": "text/x-python",
 "name": "python",
 "nbconvert_exporter": "python",
 "pygments_lexer": "ipython3",
 "version": "3.9.1"
 }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}