Система уравнений для расчета скорости естественного упаривания, страница 9

read(tt);writeln;

Qcp:=(1-Bt)*Af1[4,j]*Af1[7,j]*Af1[2,j];

tp:=Af1[3,j]+Qcp/(Hp*Pp*Cp*Af1[2,j])-0.5;

clrscr;

writeln('   Mesyats ',j);writeln;

1:

if tp<(Af1[3,j]) then begin writeln(' Viparivanie nevozmojno !');

goto 2; end;

Sm:=0; Scp:=0;

for y:=1 to a do begin

Fi:=0; Qi:=0;

Fi:=z[7,y]/1000000+z[8,y]/exp(8*ln(10))*tp+z[9,y]/exp(9*ln(10))*tp*tp;

Qi:=z[10,y]/exp(8*ln(10))+z[11,y]/exp(8*ln(10))*tp+z[12,y]/exp(11*ln(10))*tp*tp;

Sm:=Sm+Fi*Ni[y];

Scp:=Scp+Qi*Ni[y]*Ns;

end;

ep:=(7.0185471/10000*tp*tp*tp*tp+0.0024391*tp*tp*tp+2.0528933*tp*tp+37.4777559*tp+

+636.01788989)*exp((Sm+Scp)*ln(10));

ПРИЛОЖЕНИЕ А (продолжение)

E:=(1.4/1000000)*Af1[2,j]*(ep-Af1[6,j])*(1+0.72*Af1[5,j]);

Qnp:=Hp*Pp*Cp*(tp-Af1[3,j])*Af1[2,j];

Qicp:=2.26*exp(9*ln(10))*E;

if abs((Qnp+Qicp-Qcp)/Qcp)<=0.05  then begin writeln

('   Raschyot okonchen:');

goto 2;

end

else  begin

tp:=tp-tt;

goto  1;

end;

2:

for y:=1 to a do begin write('  Xicp[',y,']=',Xicp[y]:3:6,'    Ni[',y,']=',Ni[y]:3:8);

writeln;

writeln('  Xccp=',Xccp:3:6);

writeln('  Xvcp=',Xvcp:3:6);

writeln('Ns=',Ns:3:8);

writeln('Cp=',Cp:7:4);

writeln;

writeln(' tp=',tp:2:4,';  Qcp=',Qcp:11:4,';  Qicp=',Qicp:11:4,';   Qnp=',Qnp:11:4);

writeln(' ep=',ep:7:4,  ';  E=',E:5:6, ‘ ;  t=',Af1[3,j]:2:2,';  B=',Af1[4,j]:3:1,                              ';  w=',Af1[5,j]:1:2,  '; e=',Af1[6,j]:4:1););

end;

end.


ПРИЛОЖЕНИЕ Б

Примечание: знак  ##  в результатах расчетов свидетельствует о том, что проведение расчетов в пределах первой итерации некорректно, т.к. расчетное значение давления паров воды над рассолом меньше давления паров воды в окружающем воздухе.

Состав 1:

NOMER        1          2          3          4

  KOD     NaCl  MgCl2     KCl     CaCl2

    M        58.44    95.21    74.55   110.99

   Cp        864.1    746.7    688.0    654.4

   Sv           35.9     54.8     34.4     74.5

                      1         1         1          1

  1 - est dannie dlya raschyotov koeffitsientov aktivnosti

  0 - net dannih dlya raschyotov koeffitsientov aktivnosti

Xicp[1]=0.092638                    Ni[1]=38.75487962

Xicp[2]=0.026375                    Ni[2]=6.77276720

Xicp[3]=0.007767                    Ni[3]=2.54715104

Xicp[4]=0.136974                    Ni[4]=30.17200837

Xccp=0.263754

Xvcp=0.736246

Ns=78.24680622

Cp=3279.5927

Нр=0.2м

4

tp=12.2860   Qnp=140724294.9600   Qicp=60100970.8850   Qcp=191799360.0000

ep=861.4620   E=0.026593   ;  t=6.50;  B=158.0;  w=3.74; e=690.0

5

tp=19.6449   Qnp=146894597.7000  Qicp=136828587.2300  Qcp=272555472.0000

ep=1393.8938   E=0.060544  t=13.80;  B=191.0;  w=3.53; e=1000.0

6

tp=23.8255  Qnp=168437586.5700  Qicp=149170929.2700  Qcp=304637760.0000

ep=1809.8314  E=0.066005  t=16.90;  B=207.0;  w=3.16; e=1330.0

7

tp=24.4827  Qnp=150357028.5600  Qicp=108793885.8800  Qcp=250885728.0000

ep=1884.2383  E=0.048139  t=18.50;  B=171.0;  w=2.96; e=1530.0

8

tp=23.2803  Qnp=152812068.9100  Qicp=83248715.3540  Qcp=225695376.0000

ep=1750.0783  E=0.036836  t=17.20;  B=171.0;  w=2.82; e=1470.0

ПРИЛОЖЕНИЕ Б (продолжение)

9

tp=17.5123  Qnp=124338338.3200  Qicp=22622948.1270  Qcp=141363360.0000

ep=1215.9422  E=0.010010  t=12.40;  B=129.0;  w=2.97; e=1140.0

Нр=0.3м

4

tp=10.9573;  Qnp=162613608.5500 ; Qicp=34191464.7630;  Qcp=191799360.0000;

ep=787.5448;  E=0.015129 ;  t=6.50;  B=158.0;  w=3.74; e=690.0

5

tp=18.5299  Qnp=178309816.2800  Qicp=103590483.6900  Qcp=272555472.0000

ep=1298.2100  E=0.045836  t=13.80;  B=191.0;  w=3.53; e=1000.0

6

tp=22.6503  Qnp=209784067.7900  Qicp=109797404.1200  Qcp=304637760.0000

ep=1683.1803  E=0.048583   t=16.90;  B=207.0;  w=3.16; e=1330.0

7

tp=23.4551  Qnp=186798682.1100  Qicp=73415106.6030  Qcp=250885728.0000

ep=1769.0432  E=0.032485   t=18.50;  B=171.0;  w=2.96; e=1530.0

8

tp=22.0869  Qnp=184227287.4800  Qicp=46149087.8110  Qcp=225695376.0000

ep=1625.2619  E=0.020420   t=17.20;  B=171.0;  w=2.82; e=1470.0

## 9

tp=12.3749  Qnp=2731040.6152  Qicp=-79698598.8640  Qcp=141363360.0000

ep=872.4625  E=-0.035265  t=12.40;  B=129.0;  w=2.97; e=1140.0

Нр=0.4м

4

tp=9.9430  Qnp=167477900.4600  Qicp=15800122.0520  Qcp=191799360.0000