Программа расчёта теплообменника

Обозначения, принятые в программе:

d – диаметр внутренней трубы теплообменника;

pi – ;

v1, v2 – объёмные скорости теплоносителя и хладоагента;

ro1, ro2 – плотности теплоносителя и хладоагента;

Cp1, Cp2 – теплоёмкости теплоносителя и хладоагента;

h – шаг по длине теплообменника для вывода результатов;

hl – шаг по длине теплообменника для расчёта температур и вывода соответствующих точек на графике;

K – коэффициент теплопередачи;

L – текущая длина теплообменника;

lt – координата текущей длины теплообменника на графике;

t1, t2 – температуры теплоносителя и хладоагента;

tt1, tt2 – текущие координаты температур теплоносителя и хладоагента;

n – число вычислений между выводом результатов в раздел;

gd, gm – параметры графического режима;

f1 – файловая переменная, соответствующая разделу с результатом.

Program Teplo;

Const d=0.03; v1=2.3e-4; v2=5.1e-4; pi=3.14; ro1=900.0; ro2=1000.0;

cp1=3.35e3; cp2=4.19e3; h=0.2; hl=0.01; k=4900.0;

Var n:integer;

t1,t2,lk,l,b1,b2:real;

f1:text;

Begin

assign(f1,'rteplo.pas');

rewrite(f1);

t1:=200.0; t2:=35.0; lk:=3.0;

b1:=k*pi*d/(ro1*cp1*v1);

b2:=k*pi*d/(ro2*cp2*v2);

l:=0.0;

writeln(f1,'Расчёт процесса теплообмена ');

writeln(f1,' Длина, Т1, Т2,');

writeln(f1,' м град.С град.С');

writeln(f1,l:5:1,t1:10:2,t2:10:2);

n:=0;

while l<=lk do

begin l:=l+hl;

n:=n+1;

t1:=t1-b1*(t1-t2)*hl;

t2:=t2+b2*(t1-t2)*hl;

if n=round(h/hl) then

begin

writeln(f1,l:5:1,t1:10:2,t2:10:2);

n:=0;

end; end;

Close(f1);

end.