Энерготехнология_1 / ENERGY

{$N+}
uses crt;

type vect5 = array[1..5] of real;
type vect8 = array[1..8] of real;
type vect10 = array[1..10] of real;
type vect15 = array[1..15] of real;

(***** €б室­лҐ ¤ ­­лҐ *****)
const
ni:vect5 = (0.0001, 0.0138, 0.00078, 0.0008,0.01122);{kmol/s}
nj:vect10 = (0.0034,0.0084,0.0046,0.00045,0.00019,0.00008,0.0001,0.0001,0,0);
n1:vect8= (0.858,0.002,0.001,0.001, 0 ,0.002,0.137,0.001);{¬®«м/¬®«м}
Tc = 518; Te = 880; Tf = 1130; Td = 673; Tper=999; Pf=2.5*100000;
d=0.086; dnar = 0.102;
alpha=1.01;
ha=1519.39; hb=2647.32; h1=3455; h2=2740;

(***** ”Ё§ЁзҐбЄЁҐ бў®©бвў  *****)
R=8.31451;
Tos = 298.15; { [K] }
Vm = 22.4141;
Pat = 101.325;{ [kPa] }
M:vect15=(16.04,28.05,30.07,42.08,44.09,54.09,58.12,78.12,28.01,44.01,2.02,28.01,32,18.02,17.01);
H:vect15=(-74.85,52.3,-84.67,70.41,-103.85,110.16,-126.15,82.93,-110.53,-393.51,0,0,0,-241.81,42.1);
Ex298:vect15=(830,1360,1494,2001,2150,2500,2803,3432,275.4,20.1,235,0.7,3.95,8.6,157);
S:vect15=(186.27,219.45,229.49,266.94,269.91,278.74,310.12,269.2,197.55,213.66,130.52,191.5,205.04,188.72,0);
a:vect15 =(14.32,11.32,5.72,12.44,1.72,8.08,18.23,-21.09,28.41,44.14,27.28,27.88,31.46,30,29);
b:vect15
=(0.07466,0.12201,0.17511,0.18838,0.27075,0.27322,0.30356,0.40012,0.0041,0.00904,0.00326,0.00427,0.00339,0.01071,0.004);

(***** ЏҐаҐ¬Ґ­­лҐ *****)
var
Input,Out, Qef, Qfd, Qab, Qce, Pe, Pef,
NEx_pir, NEx_gor, NEx_xr,NEx_to ,niExi,njExj,n1Ex1,n2Ex2,
x,y,nO2,nN2,Tad_gor, Snj,nB,nL,nV,T,n1I1,n2I2, T3: real;
j: word;
yi:vect5;
yj,n2,y1,y2:vect10;
ntp:longint;

function pwr(n,t:real):real;
(***** t ў б⥯Ґ­Ё n *****)
begin
pwr:=exp(n*ln(t));
end;

function Cph(j:byte; T:real):real;
var cphh:real;
begin
Cphh:=a[j]+0.5*b[j]*(T+Tos);
cph:=cphh;
end;

function Cps(j:byte; T:real):real;
var cpss: real;
begin
Cpss:=a[j]+b[j]*(T-Tos)/ln(T/Tos);
Cps:=Cpss
end;

function I( j:byte; T:real):real;
(*****ќ­в «мЇЁп ў Є„¦/Є¬®«м****)
var ii:real;
begin
II:=H[j]*1000+Cph(j,T)*(T-Tos);
I:=ii;
end;

function Ex( j:word; T,y:real):real;
var exx:real;
begin
if T = Tos
then Exx:=Ex298[j]*1000
else Exx:=Ex298[j]*1000+Cph(j,T)*(T-Tos)-Tos*Cps(j,T)*ln(T/Tos)+Tos*R*ln(y);
Ex:=Exx;
end;

procedure Ad_Gor(T0:real;var T:real);
var q,Ii,Ij:real;
begin
T:=T0;
Ii:=n1[1]*I(1,Tos)+n1[2]*I(2,Tos)+n1[3]*I(5,Tos)+n1[4]*I(7,Tos)+n1[5]*I(7,Tos);
Ii:=Ii+n1[6]*I(11,Tos)+(n1[7]+nN2)*I(12,Tos)+n1[8]*I(10,Tos)+nO2*I(13,Tos);
repeat
T:=T+0.1;
Ij:=n2[1]*I(10,T)+n2[2]*I(14,T)+n2[3]*I(12,T)+n2[4]*I(13,T)+n2[5]*I(9,T);
Ij:=Ij+n2[6]*I(15,T)+n2[7]*I(11,T);
q:=Abs(Abs(Ii)-Abs(Ij));
until q < 100;
end;

Procedure I1;
begin
(*I.1 Њ вҐаЁ «м­л© Ў « ­б *)
Writeln('I.1. Њ вҐаЁ «м­л© Ў « ­б');
nj[9]:=ni[5]-nj[8];
nj[10]:=2*ni[1]+3*ni[2]+3*ni[3]+4*ni[4]+ni[5]-2*nj[1]-2*nj[2]-3*nj[3]-3*nj[4]-3*nj[5]-3*nj[6]-5*nj[7]-nj[9];
Input:=M[2]*ni[1]+M[3]*ni[2]+M[4]*ni[3]+M[5]*ni[4]+M[14]*ni[5];
Out:=M[1]*nj[1]+M[2]*nj[2]+M[3]*nj[3]+M[4]*nj[4]+M[6]*nj[5]+M[8]*nj[6]+M[7]*nj[7]+M[9]*nj[8]+M[14]*nj[9]+M[11]*nj[10];
writeln('n(H2O) = ',nj[9]:2:4, ' kmol/s');
writeln('n(H2) = ',nj[10]:2:4, ' kmol/s');
Writeln (Input:2:5 ,' = ', Out:2:5,' kg/s');
end;

Procedure I2;
(*I.2 ќ­ҐаЈҐвЁзҐбЄЁ© Ў « ­б *)
begin
Writeln;
Writeln('I.2. ќ­ҐаЈҐвЁзҐбЄЁ© Ў « ­б');
Qef:=I(1,Tf)*nj[1]+I(2,Tf)*nj[2]+I(3,Tf)*nj[3]+I(4,Tf)*nj[4];
Qef:=Qef+I(6,Tf)*nj[5]+I(8,Tf)*nj[6]+I(7,Tf)*nj[7]+I(9,Tf)*nj[8]+I(14,Tf)*nj[9];
Qef:=Qef+I(11,Tf)*nj[10]-(I(2,Te)*ni[1]+I(3,Te)*ni[2]+I(4,Te)*ni[3]+I(5,Te)*ni[4]+I(14,Te)*ni[5]);
Writeln('Qef = ',Qef/1000:3:4, ' MW');
end;

Procedure I3;
var
j:integer;
f,n,Roe,Rof,We,Wf,PSIef,dP,Me,Mf,sni,snj,Re,Rf:real;
(*Ѓ « ­б ¬Ґе. н­ҐаЈЁЁ*)
begin
writeln;
Writeln('I.3. ѓЁ¤а ў«ЁзҐбЄЁ© а б祥⠧¬ҐҐўЁЄ ');
Snj:=0;
sni:=0;
for j:=1 to 5 do sni:=sni+ni[j];
for j:=1 to 10 do snj:=snj+nj[j];
for j:=1 to 5 do yi[j]:=ni[j]/sni;
for j:=1 to 10 do yj[j]:=nj[j]/snj;
Me:=M[2]*yi[1]+M[3]*yi[2]+M[4]*yi[3]+M[5]*yi[4]+M[14]*yi[5];
Mf:=M[1]*yj[1]+M[2]*yj[2]+M[3]*yj[3]+M[4]*yj[4]+M[6]*yj[5]+
M[8]*yj[6]+M[7]*yj[7]+M[9]*yj[8]+M[14]*yj[9]+M[11]*yj[10];
Rf:=R*1000/Mf;
Rof:=Pf/(Rf*Tf);
f:=Pi*d*d/4;
Wf:=Out/(Rof*f);
j:=1;
Pe:=Pf+320000;
repeat
Pe:=Pe+1;
Re:=R*1000/Me;
Roe:=Pe/(Re*Te);
n:=ln(Pe/Pf)/ln(Roe/Rof);
We:=Input/(Roe*f);
dP:=0.032*75*(Roe+Rof)*(We+Wf)*(We+Wf)/(d*16);
PSIef:= dP*2/(Roe+Rof);
until abs(n/(n-1)*(Pe/Roe-Pf/Rof)-((Wf*Wf-We*We)/2+PSIef))<1;
Writeln(abs(n/(n-1)*(Pe/Roe-Pf/Rof)):5:4,' = ',((Wf*Wf-We*We)/2+PSIef):5:4);
writeln('dP = ',dp:4:2);
Writeln('Pe = ',Pe:9:2, ' Pa') ;
end;

Procedure I4;
var ExQef:real;
begin
writeln;
writeln('I.4. ќЄбҐаЈҐвЁзҐбЄЁ©  ­ «Ё§ Їа®жҐбб  ЇЁа®«Ё§ ');
niExi:=Ex(2,Te,yi[1])*ni[1]+Ex(3,Te,yi[2])*ni[2]+Ex(4,Te,yi[3])*ni[3]+Ex(5,Te,yi[4])*ni[4]+Ex(14,Te,yi[5])*ni[5];
njExj:=Ex(1,Tf,yj[1])*nj[1]+Ex(2,Tf,yj[2])*nj[2]+Ex(3,Tf,yj[3])*nj[3]+Ex(4,Tf,yj[4])*nj[4]+Ex(6,Tf,yj[5])*nj[5];
njExj:=njExj+Ex(8,Tf,yj[6])*nj[6]+Ex(7,Tf,yj[7])*nj[7]+Ex(9,Tf,yj[8])*nj[8]+Ex(14,Tf,yj[9])*nj[9]+Ex(11,Tf,yj[10])*nj[10];
ExQef:=Qef*(1-Tos/((Tf-Te)/ln(Tf/Te)));
NEx_pir:=(njExj-niExi)/ExQef;
Writeln('ќЄбҐаЈҐвЁзҐбЄЁ© ЉЏ„ ЇЁа®«Ё§  = ',NEx_pir*100:3:4,' %');
end;

Procedure II1(x,y,T0:real;out: boolean);
var sni:real;
j:word;
(* Њ в. Ў « ­б Ј®аҐ­Ёп *)
begin
{CO2}n2[1]:=(n1[1]+n1[2]*2+n1[3]*3+n1[4]*4+n1[5]*5+n1[8]);
{H2O}n2[2]:=(n1[1]*2+n1[2]*2+n1[3]*4+n1[4]*5+n1[5]*6+n1[6]);
{O2} nO2:=alpha*(n2[1]+n2[2]/2-n1[8]);{} n2[4]:=(alpha-1)*nO2+0.5*n2[1]*x;
{N2} nN2:=3.7619*nO2; {} n2[3]:=nN2;
{CO} n2[5]:=n2[1]*x;
{OH} n2[6]:=n2[2]*y;
n2[1]:=n2[1]*(1-x);
n2[2]:=n2[2]*(1-y);
{H2} n2[7]:=0.5*n2[2]*y;
Sni:=0;
snj:=0;
for j:=1 to 8 do sni:=sni+n1[j];
sni:=sni+nO2+nN2;
for j:=1 to 7 do snj:=snj+n2[j];
for j:=1 to 8 do y1[j]:=n1[j]/sni;
for j:=1 to 7 do y2[j]:=n2[j]/snj;
y1[9]:=nO2/sni;
y1[7]:=nN2/sni;
if out = true then
begin
Ad_Gor(T0,Tad_gor);
writeln('T ¤.Ј®а. = ',Tad_gor:4:2,' K');
Writeln('P* CO2 = ',y2[1]*Pat:1:3,' kPa');
Writeln('P* H2O = ',y2[2]*Pat:1:3,' kPa');
end;
end;

procedure II3;
var T:real;
begin
writeln('II.3.');
T:=Tos;
n1Ex1:=0;
n1Ex1:=n1[1]*Ex(1,T,y1[1])+ n1[2]*Ex(2,T,y1[2])+ n1[3]*Ex(5,T,y1[3])+
n1[4]*Ex(7,T,y1[4])+ n1[6]*Ex(11,T,y1[6])+ (n1[7]+nN2)*Ex(12,T,y1[6])+
n1[8]*Ex(10,T,y1[7])+ nO2*Ex(13,T,y1[8]);
T:=Tad_gor;
n2ex2:=0;
n2Ex2:=n2[1]*Ex(10,T,y2[1])+ n2[2]*Ex(14,T,y2[2])+ n2[3]*Ex(12,T,y2[3])+ n2[4]*Ex(13,T,y2[4])+ n2[5]*Ex(9,T,y2[5])
+n2[6]*Ex(15,T,y2[6])+ n2[7]*Ex(11,T,y2[7]);
NEx_gor:=n2Ex2/n1Ex1;
Writeln('ќЄбҐаЈҐвЁзҐбЄЁ© ЉЏ„ Ј®аҐ­Ёп = ',NEx_gor*100:1:4,' %');
end;

procedure III1;
begin
writeln('III.1. P бзҐв ЇҐзЁ ЇЁа®«Ё§ ');
T:=Tos;
n1I1:=n1[1]*I(1,T)+ n1[2]*I(2,T)+ n1[3]*I(5,T)+ n1[4]*I(7,T)+
n1[6]*I(11,T)+ n1[7]*I(12,T)+ n1[8]*I(10,T)+ nO2*I(13,T)+
nN2*I(12,T);
T:=Tper;
n2I2:=n2[1]*I(10,T)+ n2[2]*I(14,T)+ n2[3]*I(12,T)+ n2[4]*I(13,T)+
n2[5]*I(9,T)+ n2[6]*I(15,T)+ n2[7]*I(11,T);
nB:=-Qef/(n2I2-n1I1);
Writeln('ђ б室 в®Ї«Ёў  nB = ',nB:2:4, ' kmol/s');
nL:=nB*(nO2+nN2);

II1(0,0,2000,false);

Writeln('ђ б室 ў®§¤ге  nL = ',nL:2:4, ' kmol/s');
nV:=Snj*nB;
Writeln('ђ б室 в®Ї®з­ле Ј §®ў nV = ',nV:2:4, ' kmol/s');
T:=Tper;
n2ex2:=0;
n2Ex2:=n2[1]*Ex(10,T,y2[1])+ n2[2]*Ex(14,T,y2[2])+ n2[3]*Ex(12,T,y2[3])+ n2[4]*Ex(13,T,y2[4]);

NEx_xr:=(njExj-niExi)/(nb*n1Ex1-nb*n2Ex2);
Writeln('ќЄбҐаЈҐвЁзҐбЄЁ© ЉЏ„ ЇҐзЁ ЇЁа®«Ё§  ',NEx_xr*100:3:4,' %');
end;

procedure IV;
begin
writeln('IV.');
NEx_to:=NEx_xr/(NEx_gor*NEx_pir);
Writeln('ќЄбҐаЈҐвЁзҐбЄЁ© ЉЏ„ ⥯«®®Ў¬Ґ­  ',NEx_to*100:3:4,' %');
end;

procedure V1;
var nBdop,Qper,Lmex,mH2O:real;
begin
Writeln('V.1. ќ­ҐаЈҐвЁзҐбЄЁ© Ў « ­б ‡€Ђ');
Qfd:=I(1,Td)*nj[1]+I(2,Td)*nj[2]+I(3,Td)*nj[3]+I(4,Td)*nj[4]+I(6,Td)*nj[5];
Qfd:=Qfd+I(8,Td)*nj[6]+I(7,Td)*nj[7]+I(9,Td)*nj[8]+I(14,Td)*nj[9]+I(11,Td)*nj[10];
Qfd:=Qfd-(I(1,Tf)*nj[1]+I(2,Tf)*nj[2]+I(3,Tf)*nj[3]+I(4,Tf)*nj[4]+I(6,Tf)*nj[5]);
Qfd:=Qfd-(I(8,Tf)*nj[6]+I(7,Tf)*nj[7]+I(9,Tf)*nj[8]+I(14,Tf)*nj[9]+I(11,Tf)*nj[10]);
Writeln('Qfd = ',Qfd/1000:3:4, 'MW');
mh2o:=Qfd/(ha-hb);
Writeln('ђ б室 ЈаҐо饣® Ї а  mH2O =',mH2O:2:4,' kg/s');
Writeln;
Writeln('V.2. ђ бзҐв ¬Ґе. ¬®й­®бвЁ вгаЎЁ­л');
Lmex:=4*mH2O*(h1-h2);
Writeln('Lmex = ',Lmex/1000:3:4,' MW');
Qper:=mH2O*(h1-hb);
Writeln('Qper = ',Qper/1000:3:4,' MW');
nBdop:=-Qper/(n2I2-n1I1);
writeln;
Writeln('„®Ї. ђ б室 в®Ї«Ёў  nBdop = ', nBdop:2:4,' kmol/s');
nB:=nB+nBdop;
Writeln('ЋЎйЁ© а б室 в®Ї«Ёў  nB = ',nB:2:4, ' kmol/s');
nL:=nL+nB*(nO2+nN2);
II1(0,0,2000,false);
Writeln('ЋЎйЁ© а б室 ў®§¤ге  nL = ',nL:2:4, ' kmol/s');
nV:=nV+Snj*nB;
Writeln('ЋЎйЁ© а б室 в®Ї®з­ле Ј §®ў nV = ',nV:2:4, ' kmol/s');

end;

procedure VI1;
var teta, dt, at_obm, ev,ef,eh,SH,delta,H,HL,HF,HS,hm,sF,beta,Rofh,
Qp, Qc, Qk :real;

begin
teta:=(Te+Tf)/2+30;
dT:=abs(Tper-teta);
ev:=2/(1+2.15*alpha);
ef:=0.9;
eh:=0.9;
SH:=2*dnar;
delta:=2+SH*0.866;
ntp:=round(206/(dnar*dnar*Pi*10/4));
H:=(SH*(ntp-1)+dnar);
HL:=H*1.72;
hm:=SH*ntp+2*dnar;
sF:=delta*10+2*hm*9.5;
HF:=sF-HL;
if (HF/HL >=0) or (HF/HL<=2) then ROfh:=HL/(HF+HL) else ROfh:=HL/HF;
beta:=1/(1+ev/((1-ev)*eh*Rofh));
HS:=ev*(eh*HL+beta*ef*HF)/0.85;
Qp:=HS*5.66*(pwr(4,Tper/100)-pwr(4,teta/100));
at_obm:=2.09*pwr(0.25,dT);
Qc:=at_obm*206*(Tper-teta);
Qk:=Abs(Qp+Qc);
Writeln('VI.1. Џ®ўҐа®з­л© а бзҐв а ¤Ё ­в­®© §®­л ЇҐзЁ');
Writeln('4*Qef = ',4*Qef/1000:2:3,' MW');
Writeln('Qk = ',Qk/1000000:2:3, ' MW');
end;

procedure VI21;
var q,Ii,Ij:real;
begin
Qce:=(I(2,Te)*ni[1]+I(3,Te)*ni[2]+I(4,Te)*ni[3]+I(5,Te)*ni[4]+I(14,Te)*ni[5]);
Qce:=Qce-(I(2,Tc)*ni[1]+I(3,Tc)*ni[2]+I(4,Tc)*ni[3]+I(5,Tc)*ni[4]+I(14,Tc)*ni[5]);
{Qce:=Qce*}
writeln;
writeln('VI.2.1 ќ­ҐаЈҐвЁзҐбЄЁ© Ў « ­б Є®­ўҐЄжЁ®­­®© §®­л ЇҐзЁ');
Writeln('Qce = ',Qce/1000:2:4,' MW');
Ii:=nB*(n2[1]*I(10,Tper)+n2[2]*I(14,Tper)+n2[3]*I(12,Tper)+n2[4]*I(13,Tper));
T3:=Tper-452.8;
q:=Qce-Ii;
repeat
T3:=T3-0.05;
Ij:=nB*(n2[1]*I(10,T3)+n2[2]*I(14,T3)+n2[3]*I(12,T3)+n2[4]*I(13,T3));
until abs(abs(Ij)-abs(q)) < 0.5;
Writeln('T3 = ',T3:3:2,' K');
end;

procedure VI22;
var Mm,E,f,ak,W,MV,Tcp,
dTcp,Seff,PCO2, PH2O,aCO2,aH2O,ar,ateor,areal:real;
begin
writeln;
writeln('VI.2.2. Џ®ўҐа®з­л© а бзҐв Є®­ўҐЄжЁ®­­®© §®­л ЇҐзЁ');
Tcp:=(Tper+T3)/2;
{Write('Tcp = ',T-273:3:0, ' C; E = ');readln(E);}
E:=23.9;
MV:=M[10]*y2[1]+M[14]*y2[2]+M[12]*y2[3]+M[13]*y2[4];
Mm:=(2*dnar)*(4-1)+dnar+0.05;
f:=(Mm-4*dnar)*10;
W:= nV*MV/f;
ak:=0.34*E*pwr(0.6,W)/pwr(0.4,dnar);
Writeln('ak = ',ak:2:4);
dTcp:=((Tper-Te)-(T3-Tc))/ln((Tper-Te)/(T3-Tc));
Seff:=(3.49*2-4.1)*dnar;
PCO2:=Seff*y2[1]*Pe/100000;
{Write('P''CO2 = ',PCO2:3:2,' bar; aCO2 = ');readln(aCO2);}
aCO2:=46;
PH2O:=Seff*y2[2]*Pe/100000;
{Write('P''H2O = ',PH2O:3:2,' bar; aH2O = ');readln(aH2O);}
aH2O:=54;
ar:=aCO2+aH2O;
ateor:=1.1*(ak+ar);
areal:=Qce*3600/(225*dTcp);
Writeln('ateor. = ',ateor:2:4,' areal = ',areal:2:4);
end;

begin
clrscr;
Writeln('I. ЏЁа®«Ё§ нв ­®ў®© да ЄжЁЁ');
I1;I2;I3;I4;
Writeln('_______________________________________________________________________________');
writeln('II.Џа®жҐбб Ј®аҐ­Ёп в®Ї«Ёў ');
writeln('II.1.');
{‘⥯Ґ­м ¤Ёбб®жЁ жЁЁ CO2=}x:=8/100;
{‘⥯Ґ­м ¤Ёбб®жЁ жЁЁ H2O=}y:=2.4/100;
II1(x,y,2232,true);
II3;
READLN;
clrscr;

III1;
Writeln('_______________________________________________________________________________');
IV;
Writeln('_______________________________________________________________________________');
Writeln('V. ђ §а Ў®вЄ  н­ҐаЈ®бЎҐаҐЈ о饩 вҐе­®«®ЈЁЁ');
V1;

readln;
clrscr;
Writeln('VI. Џ®ўҐа®з­л© а бзҐв ЇҐзЁ ЇЁа®«Ё§ ');
VI1;
VI21;
VI22;
readln;
end.