Расчет теплообменика

Содержание

 

                                                                                                                     лист

Задание на проектирование  3

Введение  4

1  Тепловой расчет кожухотрубчатого аппарата  7

2  Гидравлический расчет кожухотрубчатого аппарата 15

3  Определение толщины тепловой изоляции 19

4  Механический расчет  кожухотрубчатого аппарата 21

5  Техника безопасности 24

Список использованной литературы 25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЗАДАНИЕ  НА  ПРОЕКТИРОВАНИЕ

 

Спроектировать горизонтальный кожухотрубчатый теплообменник для конденсации D, кг/с, перегретого пара с температурой tпп, ºС, и давлением Р, МПа. Охлаждающая жидкость поступает в трубки с начальной температурой tн, ºС.

 

Исходные данные:

– вид пара                                                                            водяной;

– охлаждающая жидкость                                                  морская вода;

– производительность                                                         D = 1,2 кг/с;

– температура перегретого  пара                                       tпп = 90ºС;

– давление перегретого  пара                                              Р = 0,04 МПа;

– начальная температура  охлаждающей жидкости          tн = 10ºС.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ВВЕДЕНИЕ

 

Теплообменники –  аппараты, в которых происходит теплообмен между рабочими средами независимо от их технологического или энергетического назначения (подогреватели, выпарные аппараты, конденсаторы, пастеризаторы, испарители и др.).

Классификация теплообменников  возможна по нескольким признакам.

По способу передачи тепла различают теплообменники смешения, в которых рабочие среды непосредственно соприкасаются или перемешиваются, и поверхностные теплообменники – рекуператоры, в которых тепло передается через поверхность нагрева – твердую (металлическую) стенку, разделяющую эти среды.

По основному назначению различают подогреватели, испарители, холодильники, конденсаторы.

В зависимости от вида рабочих сред различают теплообменники:

– жидкостно-жидкостные – при теплообмене между двумя  жидкими средами;

– паро-жидкостные –  при теплообмене между паром  и жидкостью;

– газо-жидкостные –  при теплообмене между газом  и жидкостью.

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

В качестве теплоносителя  наиболее широко применяется насыщенный или слегка перегретый водяной пар. В смесительных аппаратах пар  обычно впускают под уровень жидкости, при этом конденсат пара смешивается с продуктом, что не всегда допустимо. В поверхностных аппаратах пар конденсируется на поверхности теплообмена, и конденсат удаляется отдельно от продукта с помощью водоотводчиков.

Водяной пар, как теплоноситель, обладает рядом преимуществ: легкостью транспортировки по трубам и регулирования температуры, высокой интенсивностью теплоотдачи.

Общим недостатком парового и водяного обогрева является быстрый  рост давления при повышении температуры.

Для нагревания и охлаждения сред разработаны теплообменники разнообразных конструкций.

1) Теплообменники с рубашками. Обладают низкими коэффициентами теплопередачи, обусловленными малой скоростью движения жидкого теплоносителя в сечении рубашки.

2) Элементные теплообменники. Примером такого теплообменника является теплообменник типа «труба в трубе». При подборе соответствующих параметров можно сообщить рабочим телам желательную скорость и достичь высоких значений коэффициента теплопередачи. Достоинством таких теплообменников является соблюдение противотока, что обеспечивает наиболее полное использование теплоносителя. Недостатком является громоздкость и большой расход материала.

3) Погружные трубчатые теплообменники. Коэффициент теплопередачи таких теплообменников невелик, так как жидкость снаружи змеевика движется с малой скоростью или вовсе не движется.

4) Оросительные теплообменники. Снаружи трубы омываются жидкостью, которая вытекает из регулируемой щели на верхний виток трубы и стекает на нижние витки. Недостатком является то, что при большом количестве стекающей воды часть её сливается мимо трубок и не участвует в теплообмене.

5) Теплообменники с плоскими поверхностями нагрева. К этим теплообменника относятся различные ребристые, пластинчатые и другие теплообменники. Оребрение поверхности производится с целью увеличения поверхности соприкосновения с рабочей средой.

6) Кожухотрубчатые теплообменники. Эти теплообменники представляют собой пучок труб, помещенных в цилиндрической камере (кожухе). Полость камеры является межтрубным пространством. Трубы развальцованы в трубных решетках, которые ограничивают камеру с двух сторон. К трубным решеткам крепятся распределительные коробки с патрубками для впуска рабочей жидкости, протекающей внутри трубок. Камера снабжена также патрубками для подвода и отвода второго рабочего тела. Трубные решетки могут быть наглухо приварены или прикреплены к корпусу. Одна из решеток может быть не соединена с камерой. В этом случае уплотнение достигается резиновыми прокладками, зажимающими щель между корпусом и решеткой. Такой способ крепления обеспечивает свободное удлинение трубок при нагревании и предохраняет от нарушения крепления трубок в решетках. Кожух теплообменника обычно стальной, цилиндрической формы.

Большим недостатком  одноходовых теплообменников является несоответствие между пропускной способностью пучка труб и площадью поверхности теплообмена. Поэтому по трубкам приходится пропускать меньше жидкости, уменьшая её скорость движения.

Улучшение трубчатых  теплообменников достигается путем  группировки труб в отдельные пучки (ходы). Рабочая жидкость проходит через трубное пространство в несколько ходов, протекая последовательно через все пучки трубок. При этом скорость жидкости при данном расходе за счет уменьшения сечения проточной части возрастает, и коэффициент теплопередачи увеличивается.

Все трубчатые теплообменники выполняются  двух типов – вертикальные и горизонтальные.

 

 

 

 

 

1  ТЕПЛОВОЙ  РАСЧЕТ  КОЖУХОТРУБЧАТОГО  АППАРАТА

 

Температура конденсации  водяного пара при заданном давлении

tкп = 75ºС.

Удельная теплота конденсации  пара при температуре конденсации

r = 2320 кДж/кг.

Тепловая нагрузка Q в соответствии с заданными условиями

Q = Dr,


Q = 1,2 ∙ 2320 = 2784 кВт.

Во избежание значительного образования накипи конечную температуру охлаждающей воды ([3], с. 213) не рекомендуется принимать выше (45÷50)ºС. Принимаем:                                                

tк = 50ºС.

Из уравнения теплового  баланса находим расход охлаждающей воды

G = Q/[cв(tк – tн)],


где  св – теплоемкость воды, св = 4,2 кДж/(кг∙ºС);

G = 2784/[4,2∙(50 – 10)] = 16,6 кг/с.

Средняя температура охлаждающей воды

tв = 0,5(tк + tн),


tв = 0,5∙(50 + 10) = 30ºС.

Разность температур потоков продукта и охладителя –  средний температурный напор – определяется как среднее логарифмическое между большей и меньшей разностями температур теплоносителей на концах аппарата:

Δtcр = ,


Δtcр =

= 47,3ºС.

Задаемся числом Рейнольдса для потока воды

Re = 2300.

 

Принимаем для охладителя трубки с размером

20 × 2,0 мм.

Определяем соотношение числа труб, n, заданного размера к числу ходов, z, теплообменника

n/z = ,


где  n – общее число труб охладителя;

z – число ходов по трубному пространству;

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

μв – динамическая вязкость воды при средней температуре,

μв = 802∙10–6 Па/с;

n/z =

= 717.

Ориентировочно принимаем ([1], табл. 2.1) коэффициент теплопередачи для воды при вынужденном движении

К = 1000 Вт/(м2∙К).

Требуемая поверхность  теплообмена при принятом коэффициенте теплопередачи

Fтр = Q/(KΔtcp);


Fтр = 2784/(1,0 ∙ 47) = 59,2 м2.

Предварительно принимаем ([1], табл. 2.3) стандартный теплообменник с ближайшим значением отношения n/z:

– диаметр кожуха  D = 800 мм,

– число ходов   z = 1;

– общее число труб  n = 717;

– длина труб  L = 2,0 м;

– поверхность теплообмена  F = 90 м2;

– площадь сечения  одного хода по трубам Fтр = 0,144 м2;

– площадь самого узкого сечения потока в межтрубном

пространстве Fмт = 0,069 м2.

Рисунок 1 – Распределение температур на стенке

 

Расчет действительного коэффициента теплопередачи К производим для стационарных условий теплообмена, когда тепловые потоки от конденсирующегося пара к стенке, через стенку и от стенки к воде равны, т.е.

q1 = qст = q2.

Так как коэффициенты теплоотдачи α1 и α2, определяющие интенсивность теплообмена между теплоносителями и стенкой, зависят от неизвестных нам температур стенки tcт1 и tcт2, определяем их методом последовательных приближений.

С целью определения  коэффициента теплоотдачи от горячей стенки к охлаждающей воде α2 задаемся в первом приближении

Δ't2 = 28 град.

Находим термодинамические  критерии:

– Рейнольдса                        

Re = ,


Re =

= 2298;

– Прандтля

Pr = свμв/λв,


где  λв – теплопроводность воды при средней температуре,

λв = 0,618 Вт/(м∙К);

Pr = 4,2∙103 ∙ 802∙10–6/0,618 = 5;

– Грасгофа

Gr = gdв3ρв2βΔ't2/μв2,


где  g – ускорение свободного падения, g = 10 м/с2;

ρв – плотность воды при средней температуре, ρв = 996 кг/м3;

β – коэффициент объемного  расширения воды, β = 3,21∙10–4;

Gr = 10 ∙ 0,0163 ∙ 9962 ∙ 3,21∙10–4 ∙ 28/(802∙10–6)2 = 567797;

При движении воды в горизонтальных трубах ламинарно (Re < 2300) критерий Нуссельта равен:

Nu = C(RePr)0,2(GrPr)0,1,


где  С – опытный коэффициент, в горизонтальных трубах С = 0,74;

Nu = 0,74∙(2298 ∙ 5)0,2∙(567797 ∙ 5)0,1 = 21.

Коэффициент теплоотдачи  к воде равен:

α'2 = λвNu/dв,


α'2 = 0,618 ∙ 21/0,016 = 811 Вт/(м2∙К).

Коэффициент теплопроводности материала труб (нержавеющей стали)

λст = 17,5 Вт/(м∙К).

Термическое сопротивление  загрязнений со стороны воды ([1], табл. 2.2)

rзв = 18∙10–5 м2∙К/Вт.

Термическое сопротивление  загрязнений со стороны пара ([1], табл. 2.2)

rзп = 18∙10–5 м2∙К/Вт.

Сумма термических сопротивлений  стенок труб из нержавеющей стали и загрязнений со стороны воды и пара

Σ(δ/λ) = δ/λст + rзв + rзп,


Σ(δ/λ) = 0,002/17,5 + (18 + 18)∙10–5 = 47∙10–5 м2∙К/Вт.

Перепад температур на стенке из условия установившегося процесса

Δ'tст = α'2Δ't2Σ(δ/λ),


Δ'tст = 811 ∙ 28 ∙ 47∙10–5 = 10,7ºC.

Разность температур конденсирующегося пара и стенки со стороны пара

Δ't1 = Δtcp – Δ't2 – Δ'tст,


Δ't1 = 47,3 – 28,0 – 10,7 = 8,6ºС.

Коэффициент теплоотдачи  пара, конденсирующегося на пучке  горизонтально расположенных труб

α'1 = аελк ,


где  а – коэффициент, для горизонтальных труб

а = 0,645;

ε – коэффициент, при числе труб n > 100 

ε = 0,6;

λк – теплопроводность конденсата при температуре tкп,

λк = 0,671 Вт/(м∙К);

ρк – плотность конденсата при температуре tкп,

ρк = 975 кг/м3;

r – удельная теплота конденсации при температуре tкп,

r = 2320 кДж/кг;

μк – вязкость конденсата при температуре tкп,

μк = 381∙10–6 Па∙с;

α'1 = 0,645 ∙ 0,6 ∙ 0,671∙

= 13068 Вт/(м2∙К).

Проверим правильность первого приближения по равенству удельных тепловых нагрузок

q = αΔ't,


q'1 = 13068 ∙ 8,6 = 112385 Вт/м2,

q'2 = 811 ∙ 28 = 22708 Вт/м2.

Равенство удельных тепловых нагрузок не соблюдается:

q'1 ¹ q'2.

Для второго приближения  принимаем 

Δ''t2 = 33 град

и повторяем расчет по тем же формулам:

Re = 2298;

Pr = 5;

Gr = 10 ∙ 0,0163 ∙ 9962 ∙ 3,21∙10–4 ∙ 33/(802∙10–6)2 = 669189;

Nu = 0,74∙(2298 ∙ 5)0,2∙(669189 ∙ 5)0,1 = 22;

α''2 = 0,618 ∙ 22/0,016 = 850 Вт/(м2∙К);

Δ''tст = 850 ∙ 33 ∙ 47∙10–5 = 13,2ºС;

Δ''t1 = 47,3 – 33 – 13,2 = 1,1ºС;

α''1 = 0,645 ∙ 0,6 ∙ 0,671∙

= 12180 Вт/(м2∙К);

q''1 = 12180 ∙ 0,8 = 13398 Вт/м2;

q''2 = 850 ∙ 33 = 28050 Вт/м2;

q''1 ¹ q''2.

Для расчета в третьем  приближении строим графическую  зависимость (рисунок 2) удельной тепловой нагрузки q от разности температур Δt2 между стенкой и охлаждающей водой.

Рисунок 2 Графическая зависимость q = f(Δt2)

 

Пересечение двух прямых дает примерное значение той разности температур Δt2, при которой возможно равенство q1 = q2. Так как действительные зависимости q = f(Δt) не являются прямолинейными, для «истинной» разности температур Δt2и, взятой с графика, проводим третий расчет, аналогичный первым двум.

Для третьего приближения принимаем

Δt2 = 32,2ºС;

Re = 2298;

Pr = 5;

Gr = 10 ∙ 0,0163 ∙ 9962 ∙ 3,21∙10–4 ∙ 32,2/(802∙10–6)2 = 652966;

Nu = 0,74∙(2298 ∙ 5)0,2∙(652966 ∙ 5)0,1 = 22;

α2 = 0,618 ∙ 22/0,016 = 850 Вт/(м2∙К);

Δtст = 850 ∙ 32,2 ∙ 47∙10–5 = 12,9ºС;

Δt1 = 47,3 – 32,2 – 12,9 = 2,2ºС;

α1 = 0,645 ∙ 0,6 ∙ 0,671∙

= 12280 Вт/(м2∙К);

q1 = 12280 ∙ 2,2 = 27016 Вт/м2;

q2 = 850 ∙ 32,2 = 27370 Вт/м2;

q1 ¹ q2.

Определим величину расхождения  между тепловыми нагрузками

Δq = 100∙(q1 – q2)/q1,


Δq = 100∙(27016 – 27370)/27016 = –1,3% < ±5%.

Расчет коэффициентов α1 и α2 заканчиваем.

Коэффициент теплопередачи

К = 1/(1/α1 + Σ(δ/λ) + 1/α2),


К = 1/(1/12376 + 47∙10–5 + 1/850) = 579 Вт/(м2∙К).

Требуемая поверхность  охлаждения

F = Q/(KΔtср),


F = 2784∙103/(579 ∙ 47,3) = 102 м2.

Окончательно выбираем ([1], табл. 2.3) стандартный теплообменник:

– размер труб   20 × 2 мм;

– диаметр кожуха                            D = 800 мм;

– число ходов           z = 6;

– общее число труб                       n = 618;

– длина труб                                        L = 3,0 м;

– поверхность теплообмена                            F = 116 м2;

– площадь сечения  одного хода по трубам Fтр = 0,020 м2;

– площадь самого узкого сечения в межтрубном

пространстве    Fмт = 0,065 м2;

При этом запас поверхности теплообмена

ΔF = 100∙(116 – 102)/102 = 14%.


 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2  ГИДРАВЛИЧЕСКИЙ  РАСЧЕТ  КОЖУХОТРУБЧАТОГО  АППАРАТА

 

Гидравлическое сопротивление  обусловлено сопротивлением трения и местными сопротивлениями, возникающими при изменениях скорости потока. Потери давления ΔР на преодоление сопротивления трения и местных сопротивлений в трубопроводах определяются по формуле:

∆Р = 0,5(λL/dэ + Σξ)ρω2,


где  λ – коэффициент  потерь на трение;

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

Σξ – сумма коэффициентов местных сопротивлений.

При ламинарном режиме (Re < 2300) шероховатость практически не влияет на коэффициент трения:

λ = 64/Re,


λ = 64/2298 = 0,028.

Скорость охлаждающей  воды в трубках

ωтр = ,


ωтр =

= 0,81 м/с.

Скорость воды в штуцерах

ωвш = ,


где  dвш – диаметр штуцера в распредколлекторе ([1], табл. 2.5), 

dвш = 0,25 м;

ωвш =

= 0,34 м/с.

В трубном пространстве имеют место следующие местные  сопротивления:

– входная и выходная камеры (ξ1 = 1,5);

– поворот между ходами (ξ2 = 2,5);

– вход в трубы и выход из них (ξ3 = 1,0).

Гидравлическое сопротивление трубного пространства

ΔРтр = λ(Lz/d)(ωтр2ρв/2) + 2ξ1(ωвш2ρв/2) + [ξ2(z – 1) + 2ξ3z](ωтр2ρв/2),

ΔРтр = 0,028∙(3 ∙ 6/0,016)∙(0,812 ∙ 996/2) + 2 ∙ 1,5∙(0,342 ∙ 996/2) +

+ [2,5∙(6 – 1) + 2 ∙ 1 ∙ 6]∙(0,812 ∙ 996/2) = 18470 Па.

Скорость пара в штуцере  на входе в межтрубное пространство

ωпш = ,


где  dпш – диаметр штуцера на входе в межтрубное пространство ([1], табл. 2.5),

dпш = 0,25 м;

ρп – плотность насыщенного пара при заданном давлении Р,

ρп = 0,25 кг/м3;

ωпш =

= 97,83 м/с.

Скорость конденсата в штуцере на выходе из межтрубного  пространства

ωкш = ,


где  dмш – диаметр штуцера на выходе из межтрубного пространства ([1], табл. 2.5),

dмш = 0,25 м;

ρк – плотность конденсата при температуре tкп,

ρк = 975 кг/м3;

ωкш =

= 0,03 м/с.

Скорость конденсата в наиболее узком сечении межтрубного пространства

ωмтр = D/(Fмтρк),


ωмтр = 1,2/(0,065 ∙ 975) = 0,02 м/с.

 

 

 

Число Рейнольдса при  течении пара в межтрубном пространстве

Reмтр = Ddн/(Fмтμк),


где  μк – вязкость конденсата при температуре tкп,

μк = 381∙10–6 Па∙с;

Reмтр = 1,2 ∙ 0,020/(0,065 ∙ 381∙10–6) = 969.

Число рядов труб, омываемых  конденсатом в межтрубном пространстве (округляем до целого числа в большую сторону)

m = ,


m =

= 18.

Число сегментных перегородок  при принятых размерах теплообменника ([1], табл. 2.6)

Х = 6.

В межтрубном пространстве имеют место следующие местные  сопротивления:

– вход пара и выход  конденсата через штуцера (ξ1 = 1,5);

– повороты через сегментные перегородки (ξ2 = 1,5);

– сопротивления трубного пучка при его поперечном омывании ξ3:

ξ3 = 3m/Reмтр0,2,


ξ3 = 3 ∙ 18/9690,2 = 13,6.

В межтрубном пространстве теплообменника потери на преодоление  сопротивлений трения малы по сравнению с потерями на преодоление местных сопротивлений, поэтому они не учитываются.

Гидравлическое сопротивление межтрубного пространства

ΔРмтр = ξ1(ωпш2ρп/2) + ξ1(ωкш2ρк/2) + ξ2Х(ωмтр2ρк/2) + ξ3(Х – 1)(ωмтр2ρк/2),


ΔРмтр = 1,5∙(97,832 ∙ 0,25/2) + 1,5∙(0,032 ∙ 975/2) + 1,5 ∙ 6∙(0,022 ∙ 975/2) +

+ 13,6∙(6 – 1)∙(0,022 ∙ 975/2) = 1810 Па.

Для обеспечения заданной производительности по конденсации пара его давление перед входом в теплообменник должно быть не менее

Рп.тр = Р + Рмтр,


Рп.тр = 0,04∙106 + 1810 = 41810 Па.

Для обеспечения заданной производительности по конденсации  пара давление охлаждающей воды перед входом в теплообменник должно быть не менее

Рв.тр = ΔРтр = 18470 Па.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3  ОПРЕДЕЛЕНИЕ  ТОЛЩИНЫ   ТЕПЛОВОЙ  ИЗОЛЯЦИИ

 

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

Толщина изоляции должна быть такой, чтобы температура на её поверхности была не более 40÷50ºС.

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

Выбираем в качестве материала для тепловой изоляции совелит, имеющий коэффициент теплопроводности

λи = 0,1 Вт/(м·К).

Ввиду незначительного  термического сопротивления стенки аппарата по сравнению с термическим сопротивлением изоляции, температуру изоляции со стороны теплообменника tст1 принимаем равной температуре перегретого пара в межтрубном пространстве

tст1 = tпп = 90ºС.

Температура окружающей среды (воздуха)

tв = 20ºС.

Температуру изоляции со стороны окружающей среды (воздуха) принимаем

tст2 = 30ºС.

Коэффициент теплоотдачи  от внешней поверхности изоляционного  материала в окружающую среду

αв = 9,3 + 0,058tст2,


αв = 9,3 + 0,058 ∙ 30 = 11,0 Вт/(м2·К).

Удельные тепловые потери в окружающую среду

qп = αв(tст2 – tв),


qп = 11,0∙(30 – 20) = 110,0 Вт/м2.

Толщина тепловой изоляции

δи = ,


δи =

= 0,055 м = 55 мм.

Принимаем толщину изоляции теплообменника

δи = 56 мм.

Тепловые потери в окружающую среду

Qп = π(D + 2δи)Lqп,


Qп = 3,14∙(0,8 + 2 ∙ 0,056)∙3 ∙ 110 = 945 Вт.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4  МЕХАНИЧЕСКИЙ  РАСЧЕТ  КОЖУХОТРУБЧАТОГО  АППАРАТА

 

Толщина стенки корпуса воздухоохладителя

δс = РD/(2[σ]φ) + с,


где  Р – внутреннее давление в корпусе, Р = 0,04 МПа;

D – диаметр корпуса аппарата, D = 0,8 м;

[σ] – допускаемое напряжение, для стали Х18Н10Т  [σ] = 140 МПа;

φ – коэффициент прочности сварного шва, φ = 0,8;

с – прибавка на коррозию, овальность и т.д., с = 0,001;

δс = 0,04 ∙ 0,8/(2 ∙ 140 ∙ 0,8) + 0,001 = 0,0011 м.

Принимаем толщину стенки корпуса равной:

δс = 2 мм.

Толщина эллиптической крышки

δк = DR/(2[σ]φ – 0,5Р) + с,


где  R – радиус кривизны в вершине днища,

R = 0,25Dв2/h, ;


Dв – внутренний диаметр корпуса,

Dв = D – 2δс,


Dв = 0,8 – 2 ∙ 0,002 = 0,796 м;

h – высота выпуклой части наружной поверхности крышки,

h = 0,2Dв,


h = 0,2 ∙ 0,796 = 0,159 м;

R = 0,25 ∙ 0,7962/0,159 = 0,996 м;

δк = 0,8 ∙ 0,996/(2 ∙ 140 ∙ 0,8 – 0,5 ∙ 0,04) + 0,001 = 0,0046 м.

Принимаем толщину эллиптической  крышки равной:

δк = 5 мм.

Общая сила нормального  давления на прокладку

N = π(D2 – Dв2)σу/4,


где  σу – удельное давление на прокладку, для резины σу = 3,5 МПа;

N = 0,785∙(0,82 – 0,7962)∙3,5 = 0,018 МН.

Усилие, отрывающее крышку от фланца,

Q = PFкр,


Расчет теплообменика