Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Gelperin_N_I_-_Osnovnye_protsessy_i_apparaty_h.rtf
Скачиваний:
71
Добавлен:
08.09.2019
Размер:
10.7 Mб
Скачать

5. Теплообмен в трубках Фильда

Возможны два основных варианта теплообмена в трубках Фильда, отличающиеся тем, что в одном из них теплоноситель входит в узкую трубку и выходит из кольцевого пространства, а во втором — наоборот. При этом второй теплоноситель, омы­вающий наружную поверхность широкой трубки, движется про­тивотоком или прямотоком относительно потока в кольцевом про­странстве (рис. VI1-23, а).

Вариант 1. Внешний теплоноситель движется противотоком относительно теплоносителя в кольцевом пространстве (рис. VII-23, б). Обозначив коэффициенты теплопередачи через стенку узкой трубки Kt и через стенку широкой трубки Kit можно написать следующие уравнения:

Wt dt2 = Ki (в — t2) dx (a); —WidQ = Ka(l1 e)'dx (6); —W2 dtt = Wt dt2 Wx dO (в)

Дифференцируем уравнения (а) и (б):

Подставляя в уравнение (д) значение dtx из уравнения (в) и dQ/dx из уравнения (г), получаем:

dx2 ~ Kt \ Wt Wa ) dx2 ~*~ Wt ' dx Продифференцировав еще раз уравнение (г) и подставив най­денное значение d2Q/dx2, найдем:

л_ (Jk. л. Л* К2 \ d2t2 KLK2 dt2 _Q

dx3 \ Wt Wa / dx2 w2 ' dx

откуда 1 1 i

4 = Aeax -f- Bebx + С (e)

Здесь А, В, С — постоянные величины, а а и b — являются

корнями характеристического уравнения:

(а, Ь)1 + (KlfWl + K2/W2 - Kjj/W,) (а, Ь) - КХК21^\ = 0 (ж)

Подставив в уравнение (а) значения 4 и dt2/dx из уравне­ния (е), мы найдем температурный профиль теплоносителя в коль­цевом пространстве:

6 = Аеах + ВеЬх + С + (№t/Ki) (Ааеах + ВЬеЬх) (з)

= 4 = 4 + (Wy/Ci) + ВЪ)\ при д: = I имеем 4 = 9, поэтому Aaeal + Bbebl + С.

Путем совместного решения последних трех уравнений полу­чим:

A = [K1(f2- Q е! WyiWi (11 <6~а> «)]; В A (a/b) e~l <ь-°>;

С = (А+В)

После подстановки значений А, В и С в уравнения (е), (з) и (и) находим в окончательном виде уравнения температурных про­филей обоих теплоносителей:

„Д-А [<'--« -(4)"~>]

(VII. 18а)

_^(е«еМь-°)-е^] (VI 1.186)

+('+^)('-t)«"-4«'»-41+^)(1-^)'"]

(VII. 18в?

Полагая в уравнении (VII.186) х = I, находим температуру теплоносителя в точке его перехода из узкой трубки в кольцевое пространство:

Если по условию задачи требуется определить длину соосных труб, необходимую для достижения заданных конечных темпера­тур теплоносителей, то удобнее всего воспользоваться уравне­нием (VII.18b); положив в нем х = 0, после простых алгебраиче­ских действий найдем:

I = [\l(b - a)) In {[К2 (<; - О + (bWx + /(,) (?2 - ЩЦк2 (*[ - Q +

+ (aW1+Kl)(f2-Q-\} (УИ.18д)

dx2 W\ ' dx

В случае прямотока теплоносителей в кольцевом пространстве и вне труб (пунктирная линия на рис. VI1-23, б) уравнения (а) и (б) сохраняются, а уравнение (в) будет иметь следующий вид: W% dtx = Wx dt2 Wt dQ. Используя тот же метод решения си­стемы уравнений, что и в предыдущем случае, получаем:

dx3 ^KWx W2 Wx I

откуда t2 =A1ea**+ BlebtX + CL.

Величины ах и bx являются корнями характеристического урав­нения:

(«ь bif + (Ki/^i - K2/W2 - K2/Wx) (av *t) - K,K2IW\ = 0

Уравнения двух других температурных профилей напишутся следующим образом:

6 = Aiea'x + BLeb*x + Cj + (WVKi) + BXVM)

t1 = A1ea^(l+a1W1/K1)(l-a1Wl/K2) +

+ Bie»" (I + bxWJKt) (1 -+ ct

Легко убедиться, что при переходе к граничным условиям мы получим Ах = А, В, = В и С4 = С, если заменить в их выра­жениях а и b на ах и Таким образом, уравнения (VI 1.18а), (VII. 186), (VII.18в) и (VII. 18г) сохраняются и в рассматриваемом случае, если везде заменить а и b величинами аг и Ьх. Изменится лишь выражение для требуемой длины труб:

1 = - а011п {[*2 (1 - Q+(Ki+biwi) - 2 C'i - Q +

-f-^+a^) (<;-<;)]} (VH.lSe)

В частном случае, когда наружная поверхность широкой трубы омывается теплоносителем постоянной температуры (tx const), исчезает разница между противотоком и прямотоком. Тогда Wx = = оо, а = at = K2/Wx, b = bx= —K\l^\- Подставляя эти зна­чения а и b в предыдущие уравнения, получим применительно к tx = const следующие выражения:

+ [Kll(Kl + K2)\(tl-t"2)(e-(KtIW1)X_e(K,IWi)x^ (VII.l&K)

e^fj-ft — QeCKt/*!)* (VI 1.18з)

= Q e^"™' (VII.18H)

' = + Ч\U*2 (<1 - Q + *1 (<2 - <2)l/[*2 ('l - Q])

(VI I. 18k)

Располагая выведенными выше уравнениями для темпера­туры f2 и обозначив диаметры узкой и широкой труб через dx и d2, воспользуемся для расчета общим уравнением теплопередачи:

Q = w, (t\ - Q = wx {}\ -1'2) + wl (f2 - Q =

где Дер и Аср — средние разности температур для узкой и широкой труб, вы­числяемые обычным способом.

Разумеется все выведенные уравнения справедливы как при нагревании внутреннего потока теплоносителя (t[ *> t2), так и при его охлаждении (t[ < ин­вариант 2. В данном варианте один теплоноситель входит в кольцевое пространство и уходит из узкой трубы, а второй, как и в предыдущем случае, омывает наружную поверхность широкой

363

трубы. Если теплоноситель в кольцевом пространстве движется противотоком относительно второго (внешнего) теплоносителя (рис. VI1-23, в), то процесс теплообмена можно описать следую­щей системой уравнений:

= Kl(t2 Q)dx (a); Wtd% = Ка (t± — Q)dx-\- Кх (t2 Q) dx (6);

W2dtx + Wxdt2 = WtdQ (в)

Сохраняя тот же ход рассуждений, что и в предыдущем вари­анте, получим:

А-к. = Jk. ( Ah. dQ\. " dx2 wx \ dx ИГ)'

^ = _Ix(Jk-_Jk- d^ Кг dt2

dx2 Kx \W2 Wt Wx) dx2 ~~W1"~dT

J<3_ k2_\ M2 KxK2 dt2 _n

dx-3 ^KWx W2) dx2 ~~Щ dT~° t2 = Ae(r)

ax+Bebx + C (д)

Величины а и b являются корнями характеристического урав­нения:

(я, 6)а + (K2/Wt - K2/W2) (a, b) - KtKJWl = О

С помощью уравнений (а), (б), (в) и (д) мы находим уравнения двух других температурных профилей:

6 = Аеах + ВеЬх + С - (WyKJ (Ааеах + ВЬе**) (е)

^('--fc-'^H+^O-'-u-'^B+'w

Для определения значений постоянных коэффициентов А, В и С воспользуемся граничными условиями:

t2=t"2 = A + B + C и е = /2-(1^11)(Ла + В6) при х = 0

Aaeal = Bbebl = 0 при х=1

Из граничных условий вытекает, что величины А, В к С имеют те же значения, что и в первом варианте. После подстановки зна­чений Л, В и С в уравнения (д), (е) и (ж) находим в окончательном виде искомые уравнения температурных профилей теплоносителей:

t2—t2-


bwx(l—el <b-a>)

(VII.19a)


в = С + *2> Uebx - I) - el '6-fl> (eax -1)4-


bKL(el (b-a)+aX _ ^jtv] (VI1.196)

t f , Kl^-fl Г/ & fMb-g) A

-TV Kl KtK2)e +V Ki KiK2)e J

(VI 1.19b)

Для определения температуры теплоносителя в точке его пере­хода из кольцевого пространства в узкую трубу воспользуемся уравнением (VI 1.19а), положив х = h

с=с+ K;(4~l;)an fы-0-4-0ebl] <VII-19r>

bWxil el (6_a>) J

Приняв в уравнении (VII.19b) x = 0, находим требуемую длину труб для заданного процесса теплообмена;

I = [ 1/(6 - a)] In {2 {t\- Q + bWx (t'2 - g ]/[К2 (t'[ - f2) +

+ 0^,(^-/2)]} (VII.19д)

В случае прямотока теплоносителей вне широкой трубы и в кольцевом пространстве (рис. VI1-23, в — пунктирная кривая) мы располагаем следующей исходной системой уравнений:

Wx dt2 = Kt (t2 - Q) dx; Wx dQ = K2 (tx -Q)dx + Ki (t2 - Q) dx; W2 dix + Wx dt2 = Wx dQ

Решение этой системы тем же методом, что и в предыдущих случаях, приводит к следующему уравнению:

d% , / К, . Кщ\ d42 KtK2 dt^_ dx3 ~r \ W2 ~*~ Wx j dx2 W\ dx ~ Полученное уравнение отличается от уравнения (г) только коэффициентом при производной второго порядка и имеет то же решение при новых значениях а и Ь, обозначаемых ах и Ъх и опре­деляемых из характеристического уравнения:

х, Ьх)* + (K2IWX + K2/W2) (ox, bt) - КХК2 = 0

Таким образом, уравнения (VII. 19а)—(VII. 19г), выведенные для противотока, применимы также к прямотоку, если значе­ния а и b заменить значениями ах и Ьх. Выражение же для требуе­мой длины труб будет в данном случае иметь следующий вид:

u

1 = [ чь1 - in {lb ei - q+bi ri & - t:MK* -0+

+ axW1(t;-f2)}\ (VII.19e)

Наконец, при tx = const исчезает разница между прямотоком и противотоком; приведенные уравнения одинаково справедливы (W-i " 00, а = ах, Ь = 60, причем

а = -K2/2WX + (rf2/2Wx) Vl + Ki/Kv b = i/^i - (^/2^1) Vl + Kx/K2

365

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

при противотоке

«с;-<;)+»,<<;-;)>» «. -£<4ff(frf+lH

при прямотоке

*.(,;-,:)+wl((:_0>, „„ i<|ri(|5|_i_,)

Ж. НЕСТАЦИОНАРНЫЙ ПРОЦЕСС КОНВЕКТИВНОГО ТЕПЛООБМЕНА

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

Допустим, что в аппарате (рис. VI1-24) содержится несменяе­мое количество жидкости, водяной эквивалент которой равен ЦРДж/(е-К) и остается постоянным при ее нагревании или охла­ждении от начальной температуры 4 до конечной 4.. Второй тепло­носитель (греющий или охлаждающий) омывает с постоянной скоростью поверхность теплообмена, имея постоянную началь­ную температуру 6i и водяной эквивалент Wx Дж/(с-К). Совер­шенно очевидно, что с течением времени будут изменяться тем­пературы жидкости внутри аппарата (t) и теплоносителя, покида­ющего поверхность теплообмена (6), осуществляемой в виде по­груженных пучков труб, змеевиков или наружных рубашек. Жидкость внутри аппарата обычно интенсивно перемешивается, поэтому можно считать ее температуру в каждый момент времени постоянной по всему объему. Как и во всех ранее рассмотренных процессах теплообмена, примем также К const.

Допустим, что жидкость в аппарате нагревается при полном перемешивании, а греющий теплоноситель в конце некоторого участка / поверхности теплообмена понизил свою температуру от 9Х до 6. Для элемента поверхности df можно написать следую­щее уравнение теплопередачи: dQ = —Wi dQ = К (9 — f) df, от­куда d9/(9 — 0 = —{KIWL)df. Интегрируя полученное уравне­ние в пределах от 9, до 9 и от 0 до f, находим:

ln[(0_<)/(e1-f)]=-(A:/B7i)/ или в-t =(di-t)e-Kf/w' (а)

Среднее значение 9 — t для любого момента времени по всей поверхности теплообмена F вы­разится так:

-Kfivldf =


t)e


(в-0сР = -]г |(в-0#

4W

= (9i

Рнс. VII-24. К расчету нестационарного конвективного теплообмена.

Количество тепла, передаваемого через всю поверхность тепло­обмена F за время dx, составляет:

dQ = W dt = KF (6 - Оср dx = Wt (8X -1) (1 - (Г**!**) dx

x t, л

f . W f dt

откуда j dx= (l_e-jcf/y.) j 6^=10 " xx- -

Следовательно, продолжительность нагревания жидкости в ап­парате от начальной температуры 4 до конечной температуры 4 выразится так:

т = \ W/[Wt (1 - e-1

KF<w>)]\ in - ^/(Oi - 4)1 (VI1.20) Выражение (VII.20) позволяет решить также обратную за­дачу, т. е. определить температуру, до которой может быть на­грета жидкость в аппарате за любое время т:

t3 = Ql^((>l-tl)e^x/W) (.-(№/*.>_,) (VII.20a)

Располагая выражением (VI 1.20а), можно определить конеч­ную температуру греющего теплоносителя 92 в любой момент вре­мени. В самом деле, из уравнения (а) следует:

Х(1

йа = в-^' +4(1- *-KF/Wl) = 6х-(в1 - 4) е<^1^ (r-«F'W')- »> X

(VI 1.206)

. е-КГ ГШ г)

Средняя конечная температура греющего теплоносителя может быть найдена из общего Уравнения теплового баланса: Wxt (8Х — — 0ср) = W (4 к)- Дли определения средней разности темпе­ратур теплоносителей на протяжении рассматриваемого про-

цесса ЛСп воспользуемся общим уравнением теплопередачи: Q = = KAcpFr = W (t2 tx). После подстановки найденного значе­ния т [выражение (VII.20) ] получим:

AcP = (Wi//CF)(l-e-KF/u70(4-4)/{ln[ei-fl)/(ei- tt)}} (VI 1.20b)

Если греющий теплоноситель сохраняет постоянную темпера­туру (конденсирующийся насыщенный пар или кипящая инди­видуальная жидкость), то Wt = оо, и после раскрытия неопре­деленности [lim (WjKF) (1 — e-KF/Wl)-* 1 получим обычное вы­ражение:

Дср = (*, - *i)/{ln Wt - Швх - 4)]}

Заметим, что выведенные формулы справедливы не только в случае нагревания, но и в случае охлаждения жидкости в аппа­рате.

3. ОПРЕДЕЛЕНИЕ ОСНОВНЫХ РАЗМЕРОВ ПОВЕРХНОСТИ ТЕПЛООБМЕНА РЕКУПЕРАТИВНЫХ АППАРАТОВ

Поверхность теплообмена, или рабочая поверхность, тепло-обменного аппарата F определяется из основного уравнения тепло­передачи:

F = Q/KAcp (а)

Необходимо, однако, учесть, что рассчитанная поверхность теплообмена F будет удовлетворять заданию только при геометри­ческих ее размерах, обеспечивающих условия, принятые в расчете коэффициента теплопередачи /(. Так, в случае кожухотрубных аппаратов эти условия сводятся к скоростям потоков теплоно­сителей ад, диаметру труб d, их числу п, длине I и шагу, коэффи­циенту теплопроводности материала труб %ст и толщине их сте­нок s. Материал труб (следовательно, и Хст) диктуется физико-химическими свойствами теплоносителей (коррозия, температура); это могут быть чугуны, углеродистые и легированные стали, цвет­ные металлы и- металлические сплавы, различные неметаллические материалы. Диаметры труб и толщины стенок регламентированы государственными стандартами и выбираются соответственно ра­бочему давлению и требованиям технологии машиностроения; для обеспечения компактности аппарата стремятся к минимальному диаметру труб.

Число труб, их длина и скорость потока^ w связаны следую­щими зависимостями: F = ndtn и G = 3600 (jtd2/4) wnp, где G — часовой расход теплоносителя во всех трубах пучка, ар — его плотность.

В тех случаях, когда теплоносители сохраняют свое агрегат­ное состояние в процессе теплообмена, по известным значениям G, р, d, F и ад (величина ад была уже выбрана при расчете коэффи­циента теплоотдачи в трубах а2) из приведенных двух зависимостей

однозначно определяют пи/. Напомним, что при значениях /, превышающих стандартную длину труб, в зависимости от их ко­личества п применяют либо многоходовые аппараты, либо последо­вательное соединение нескольких аппаратов. Скорость теплоно­сителя в межтрубном пространстве предопределена найденным числом труб п и должна, разумеется, соответствовать принятой в расчете коэффициента теплоотдачи а2.

Если же один или оба теплоносителя изменяют свое агрегатное состояние (конденсирующийся пар, кипящая жид­кость), то усложняется расчет самого коэффициента теплопередачи. В самом деле

К= l/(l/«i+sACT+ 12) (а)

Однако при конденсации пара и ки­пении жидкости величины ах и а2 сами являются функциями температуры по­верхности теплообмена, с которой они

Рис. V11-25. К расчету коэффициентов теплопередачи при изменении агрегатного состояния одного или обоих теплоносителей.

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

Нагревание жидкости или газа паром, конденсирующимся в межтрубиом пространстве аппарата. Определение коэффициента теплоотдачи а2 от поверхности стенки (средняя температура 92) к потоку жидкости или газа по приведенным в главе VI критери­альным формулам не встречает затруднений. Для расчета коэф­фициента теплоотдачи а,х от конденсирующегося пара (темпера­тура насыщения tH = const) к противоположной поверхности стенки (средняя температура 9Х) мы располагаем формулой (VI.68), конкретный вид которой зависит от расположения теплообменных труб. Для пучка горизонтальных труб диаметром d фор­мула (VI.68), если придерживаться обозначений на рис. VI1-25,

будет иметь следующий вид: ах = 0,728<р YРж^ж^/рж d (tH — 0i), где ф —коэффициент, учитывающий число расположенных друг над другом рядов труб в горизонтальном пучке.

При выбранном диаметре труб d и ориентировочной оценке ф, учитывая зависимость рж, Хж, г и рж для любого насыщенного пара

только от его температуры, величина 0,728ф V Рж^жП^Рж^ — = А = const, поэтому

<*4 = Л/Сн-6Ч)1/4 (б)

Так как ах itu — 9i) = К Дср> т0 после подстановки в выра­жение (б) значения (tB — 9,) находим:

а1=Л4/3/х1д(в)

Выражение (а) можно написать теперь следующим образом:

ЦК = I/a, + s/%„ + I/o2 = KV3A\'3/A4'3 + s/\„ + 12

Для удобства определения искомого коэффициента теплопере­дачи К графическим способом перепишем последнее уравнение в окончательном виде:

1 - (K/Af3 Д1/3 = К (s/Kt + 1/<ч) (VI1 '21)

Применительно к вертикальным трубам высотой Нм формула (VI.68) имеет следующий вид: ах = 0,943 у'£ржА.ж/7рж/г (/„—6])= = AJ[ft1/4 (4 - во1'*], где Л = 0,943 (§рж^3жг/рж)1/4= const -для любого насыщенного пара заданной температуры tH.

Выражение (в) будет в данном случае иметь следующий вид:

аг = А\<31/3&1р1/а (г)

Если теплоноситель, движущийся внутри труб со скоростью до, имея водяной эквивалент W и нагреваясь от t'c до совершает в аппарате i ходов, содержащих каждый п труб, то

KAcpndhni = W (t"c - Q = 3600 (nd2/4) пшРжсж (t"c - Q

откуда h = [900до фжсж (t"c QVK Acpi.

Таким образом, at = Л1/3г1/3/[900шфжсж QY/3 й

I//C = [900шфжсж (fe - QyP/AiW3 + sACT + l/oj (VII .21a) где еж — удельная теплоемкость нагреваемой жидкости.

Испарение жидкостей при температуре кипения. Поверхность теплообмена наиболее распространенных испарителей представ­ляет собою пучок вертикальных труб, внутри которых кипит испаряемая жидкость, а межтрубное пространство занято грею­щим теплоносителем. Последним служит чаще всего конденси­рующийся пар, а в ряде случаев — высоконагретые газы или жидкости. Высота и диаметр труб, как мы увидим в главе VIII, выбираются исходя из технологических факторов. Искомой вели­чиной в инженерных расчетах является требуемая поверхность теплообмена (нагрева), для определения которой необходимо знать коэффициент теплопередачи К-

В случае обогрева вертикального испарителя (высота труб л) конденсирующимся паром можно воспользоваться выражением (г):

а, = АПКГ AcpV'3.

Для определения коэффициента теплоотдачи к кипящей жид­кости применим формулу (VI.58), которую напишем в следующем виде: а2 = 9ХЖ (92 - ШожужТ„ = В (62 - tc)\ где величина В = = Ь3%жж\жТн постоянна для каждой жидкости и зависит от ее свойств и температуры кипения tc. Но а224) = К Аср и 02 — ta = К Аср2, поэтому

а2==В1/3^/Зд2/3 (д)

Соответственно выражению (а), придерживаясь обозначений на рис. VII-25, получаем следующее уравнение для расчета иско­мого коэффициента теплопередачи:

1 - (KIA)"3Д»р3 = К(sAct + l/BV3A%3) (VII .216)

Заметим, что в испарителях непрерывного действия темпера­тура кипения жидкости tc обычно постоянна, а так как темпера­тура конденсации пара ^ тоже постоянна, то Дср = tH tc. Если греющим теплоносителем является газ или жидкость, то величина аг определяется по соответствующей критериальной формуле (см. главу VI), а для аг мы располагаем выражением (д). Таким образом

ЦК = I/a, +sACT + \/B^3K2/3A2J3 (VII.21B)

Гидравлическое сопротивление аппаратов. Гидравлическое сопротивление теплообменного аппарата Ар слагается из потерь давления вследствие трения, возникающего при движении теплоносителей в трубах (или межтрубном простран­стве) Арт, и потерь давления на преодоление местных сопротивлений Дрм. К числу последних относятся нижняя и верхняя камеры (пространства между трубными решетками и крышками аппарата), вход в трубы и в межтрубное пространство и выходы из них, перетекание жидкости из одной секции многоходового аппарата в другую, огибание перегородок в межтрубном пространстве. Как известно из главы I, величина Ар может быть выражена следующим образом:

Ар = Арт + Ары = (lid) + £ nMU] (о>2/2) р (е)

где X — коэффициент гидравлического сопротивления для прямых участков; пм — число одноименных местных сопротивлений; £„ — коэффициент местных сопротивлений; тир - скорость и плотность жидкости (газа). Значения £м приведены в справочниках и руководствах по теплопередаче. При определении значений Арт следует помнить, что в коротких трубах (//d<0,0575Re при лами­нарном течении и lid < 60 при турбулентном) сопротивление трення больше, чем иа участке гидравлически стабилизированного потока.

Заметим, наконец, что выражение (е) является, строго говоря, приближен­ным, так как оно игнорирует потери давления, вызванные изменением плотностей теплоносителей при их нагревании или охлаждении, а также подъемную силу в вертикальных теплообменных трубах. Эти потери, однако, существенны при значительных перепадах температур и больших высотах труб, но не играют заметной роли в большинстве теплообменных аппаратов на химических пред­приятиях.

Гидравлическое сопротивление пластинчатых и спиральных аппаратов опре­деляется по частным формулам, предложенным для каждой конструктивной мо­дификации.

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

Оптимальные рабочие режимы аппаратов. Оптимальный рабочий режим теплообменного аппарата определяется не только скоростями теплоносителей (или величинами Др), ио также перепадами их температур иа концах поверхности теплообмена. Оба эти фактора сопряжены с коэффициентом теплопередачи и средней разностью температур, предопределяя поверхность теплообмена н, следовательно, затраты на эксплуатацию аппарата. Последние слагаются, во-первых, нз амортизационных отчислений на капиталовложения (приобретение и монтаж аппарата, насосов, труб и т. д.), расходов на текущий ремонт н обслу­живание. С некоторым приближением можно принять, что эти эксплуатационные затраты пропорциональны поверхности теплообмена F, составляя а руб./(м2-год). Во-вторых, к эксплуатационным затратам относится стоимость энергии, расхо­дуемой на создание принятых скоростей потоков теплоносителей. Если мощность насосов, нагнетающих теплоносители, равна N кВт, число рабочих часов в год п и цена электроэнергии Ь руб./(кВт-ч), то второе слагаемое эксплуатационных затрат составит bnN руб./год. Суммарные эксплуатационные затраты (а/г+ bnN) противостоят стоимости рекуперированного тепла (холода) cQn руб./год, где с—цена единицы тепла (холода), a Q—тепловой поток.

Таким образом, экономическую эффективность тепло-обменного аппарата Э можно выразить следующим уравнением:

3 = cQn—(aF + bnN) (а)

При заданных водяных эквивалентах теплоносителей (Wx н W2), их на­чальных температурах (t[ и t'2) и известных значениях a, b и с величина Э яв­ляется функцией скоростей теплоносителей {wb w2) н их конечных температур (f; и Q:

Э = / (a»j, w2, tv .,) (С)

Функциональная зависимость (б) имеет очень сложное алгебраическое вы­ражение, поскольку К зависит от щ и а2, а последние, как и N — от w1 и ш2. Оптимальные значения wit w2, t\ н t2, соответствующие максимуму, должны

дЭ пдЭ_пдЭ_пдЭ__г. удовлетворять условиям: = 0; — 0; — U; — и; расчет этих

значений целесообразно производить при помощи ЭЦВМ.

Очень часто, однако, скорости wt и w2 предопределены ходом технологиче­ского процесса или возможные пределы их варьирования очень малы. Тогда уравнение (а) упрощается:

3 = cQn aF (в)

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

Q=W, ft - Q = W2 (t"2 - Q = KF J(Aj - A2)/lln (r)

причем в случае прямотока Дх = t[ t'2 и Д2 = 4 — t2, а в случае противотока д, = f ^ и А2 = t"x t'2.

Подставляя в уравнение (в) значения Q ъ F, дифференцируя его по 4 и г2,

приравнивая первые производные нулю (jffr = °> -Щ~ ~ °)' находим оп­тимальные конечные температуры теплоносителей: для прямотока

л fWt^2) (tj + Q + a/Kcn tx + {W2IW{)i2-alKcn

'i = 1 + wt/w2 2 ~ i + «V^i

для противотока

. - в ± у в*+4 (цуг2) {<; (g + g + [(<; - g//rCB] a} ^

~ 2 (Г!/Г2) '

'2 = '2 + (Wi/^2)Pi-0 (VI 1.226)

где В = 1 — 24 + 4 (1 — НУ Г*).

Если одним нз теплоносителей является конденсирующийся насыщенный пар, то, как уже отмечалось, разница между прямотоком и противотоком исче­зает. В этом случае Wx = оо, t[ = t'{=tv и из предыдущих уравнений находим:

4 = <, —a//(cn (VI 1.22 в)

И. СТАЦИОНАРНЫЙ ТЕПЛООБМЕН МЕЖДУ ПОТОКОМ ГАЗА (ЖИДКОСТИ) И СЛОЕМ ЗЕРНИСТОГО МАТЕРИАЛА ПРИ НАЛИЧИИ ВНУТРЕННЕГО ИСТОЧНИКА ТЕПЛА

Рассматриваемый процесс теплообмена часто встречается в ап­паратах для осуществления гетерогенного газового катализа. На пути движения газовой смеси через слой зернистого катализа­тора требуется в этих случаях отводить (или подводить) тепло, выделяющееся (или поглощаемое) в результате химического пре­вращения. Мы ограничимся здесь определением температурных перепадов в поперечных сечениях слоя катализатора, полагая температурный профиль по высоте слоя известным, а внутренний источник тепла равномерно распределенным по объему.-

Для упрощенного решения задачи представим себе такое одно­родное твердое тело, которое при тождественных температурных условиях и геометрических размерах проводит в единицу времени столько же тепла, сколько отдает (или поглощает) слой зернистого катализатора, пронизываемый газовым потоком, ограничиваю­щей его поверхности в результате одновременного действия тепло­проводности, конвекции и лучеиспускания. Коэффициент тепло проводности такого тела э) можно представить суммой: Кэ = — Хт + + Хл, где Хт — коэффициент теплопроводности зер­нистого материала (катализатора); — коэффициент теплопро­водности, эквивалентный коэффициенту конвективной теплоот­дачи газового потока в зернистом слое; Хл — коэффициент тепло­проводности, эквивалентный теплоотдаче излучением.

Величина Хэ, называемая, как уже отмечалось, эквивалентным коэффициентом теплопроводности, определяется в каждом кон­кретном случае экспериментальным путем. После введения Хэ к рассматриваемому зернистому слою применимо уравнение (VI.13):

а2е ffle д»в q

~W + ду2 + дг* + ТГ ~ где q удельное количество выделяющегося (поглощаемого) тепла, Вт/м(а)

3.

373

Ниже приводятся решения уравнения (а) применительно к наи­более употребительным формам поперечного сечения зернистого слоя (рис. VI1-26).

Если зернистый слой помещен в круговом цилиндре (рис. VII-26, а), омываемом снаружи охлаждающим или греющим потоком, то с учетом изменения температуры в радиальном На-

рве. VII-26. К расчету теплообмена с потоком газа в слое зернистого материала при наличии внутреннего источника тепла.

(6)

правлении целесообразно написать уравнение (а) в цилиндриче­ских координатах:

д2в

дг* г dr + Кэ rf8

Полагая —и и умножая все члены уравнения (б) на rdr,

JL

(в)

2 v

+ г

dr

получаем: d (иг) = —(<q/K) rdr. Двукратно интегрируя последнее уравнение, находим:

e=—ij[-r» + C1lnr + C1

В центре рассматриваемого сечения = 0) мы имеем = 0,

поэтому Сх = 0.

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

где R — наружный радиус цилиндра; 8д — температура наружной поверх­ности цилиндра в данном его сеченин; t — температура греющего или охлаж­дающего агента; а — коэффициент теплоотдачи от наружной поверхности ци­линдра к потоку охлаждающего (или греющего) агента.

Относя уравнения (в) к поверхности цилиндра R, 0 =

= 0Д) и подставляя из них значения 0# и (-jjjr)r=J? в уравнение (г),

получим: С2t + (qRV4k3) (1 + 2K3/aR).

С найденными значениями Сх и С2 уравнение температурного профиля в сечении зернистого слоя принимает следующий окон­чательный вид:

9 = 1 + (qR2/4%3) [1 + 2),3/aR (r/R)2] (VI 1.23)

Полагая в последнем уравнении г — 0, находим температуру в центре сечения (0О), а при г = R — температуру на его пери­ферии (0Л):

ео = * + (iRV^K) 0 + 2K/aR): е« = * + <7#/ (VI1 -23a)

Отсюда

в0-6я = Й&14Ь3 (VI 1.24а)

Если зернистый материал расположен в кольцевом простран­стве между двумя соосными цилиндрами (рис. VII-26, б), то на ограничивающих поверхностях (радиусы Rt и R2) могут устанав­ливаться разные температуры (0Х и 02). Температура будет макси­мальна на промежуточной поверхности радиуса R, т. е. Rt < < R <R2. Для определения температурного профиля в попереч­ном сечении зернистого слоя можно воспользоваться уравне­ниями (в). А так как (-^г) = 0, то С, = qR2/2X3.

Полагая, что коэффициенты теплоотдачи различны для пото­ков теплоносителя, омывающих ограничивающие поверхности слоя (а, и а2), мы будем располагать следующими граничными условиями:

*• ("STL,, = *^ - <>: ~ *• ("STL,. = «2 <9* ~ <> Отсюда при помощи уравнений (в), зная величину Си находим: С2 = - (9/2а,) (Я, - Я2/Я,) + /2\з) 2/2 - R2 In Я,) +1 = = (<?/2a2) 2- Я22) + (q/2\) (R2/2 -R2 In R2) + t

После подстановки полученного значения С2 в уравнение (в) находим искомое уравнение температурного профиля в двух вари­антах:

9 = t + q/iK3 (R2 - г2) + (qR2/2l9) In (г/Я,) -- (#2a,) (Я, - Я2/Я,) = * + (<7/4Хэ-) (Я2 - л2) +

+ (?Я2/2^э) In (г2) + (<7/2а2) 2 - Я4/R2) (VI1.25)

Температура на внутренней поверхности цилиндра (соответ­ственно г = Ri) выразится так:

в, = t + ((//2а,) 2/Я, - Я) = i + (t;/*A,) 22 - Я2) +

+ (<?Яа/2^) 1п х2) + (?/2а2) 2 - Я22) (VII .25а)

375

Температура на поверхности наружного цилиндра (соответ» ственно г = R2) будет:

92 = t + (?/2о2) (Я2 - R2/R2) = t + (д/41э) (R\ - Rl) +

+ ((^/2^) In (RJRi) + (qftaj (R*/Ri ~ Ri) (VII .256)

Аналогично получим максимальную температуру на поверх­ности радиуса R:

еМакс=t+(Л) (r\ И72ч).ш те>+

+ (?/2Л1) (/?2//?! = (,/4X3) (#2 - #2) +

+ (qR*/2k3) In (Я/Я,) + (q/2a2) (Я, - R*/R2) (VII .25в)

Из написанных выражений находим перепады температур в се-

чении слоя:

?2 D2

макс

еМакс-е2 = №)(^-/?2) + (^2/2^)1п(/?//?2) Для определения радиуса R достаточно решить уравне­ние (VII.25B):

= |Л('^э) (Rl - «О + *i/«i + %a2]/l((V1I25r)

I/4)1n(/?2/^i)+ !/<*!«, + 1/а2«21

(VII.25fl)

Наконец, если зернистый слой расположен между пластинами, омываемыми снаружи потоком охлаждающего или греющего аген­та (рис. VI1-26, в), то для определения температурного про­филя в сечении слоя исходное уравнение будет иметь следующий

вид: -^г + 17 = 0, откуда

В точке максимальной температуры, т. е. при я = 0, имеем: ^=0иС1 = 0.

Для нахождения С2 воспользуемся уравнениями (д) и гранич­ным условием: —к3 (-^j) _6 = a (0i — г)» гДе ^ — половина тол­щины зернистого слоя. Получаем: С2 = t + (q82/2k3) (1 + 2Хэ/аб). ■ С этим значением С2 находим искомое уравнение температурного профиля:

6 = t + (qbV2k3) [1 + 2k/afi - (*/°)21 (VI1 -26)

Полагая л = О, получаем выражение для максимальной тем­пературы в центре сечения слоя:

Bume = t + (q6*/2\a){l + 2X3/a6) (VII.26а)

Я =

Температура на поверхности слоя (соответственно х — 6):

q1 = t + qb/a (VI 1.266)

Максимальная разность температур в сечении слоя:

BH№e bi = q&№s (VII.26B)

К. ПРИБЛИЖЕННЫЙ МЕТОД РАСЧЕТА РЕГЕНЕРАТОРОВ

Регенераторы, как уже известно, в отлнчне от рекуперативных теплообмен­ных аппаратов, работают в нестационарном режиме. Точный расчет нх весьма сложен, поэтому в инженерной практике обычно пользуются приближенными методами, оперируя средними параметрами процесса и относя тепловые потоки не к единице времени, а к продолжительности цикла тц с. Периоды нагревания (тн) и охлаждения (т0) регенераторной насадкн, составляющие цикл (тц = тн + + т0), могут быть равными н неравными.

Обозначив средние температуры потоков теплоносителей в периоды нагре­вания и охлаждения насадкн через ?н и t0, а соответствующие средние темпера­туры насадкн — через 0Н и 90, можно с некоторым приближением применить к рассматриваемому процессу уравнения стационарной конвективной тепло­передачи:

Ян = Ки. Сн — 'о) = «н^н Сн — Вн) = аоЧ 0 — 'о) (а)

где ан н а0 — коэффициенты теплоотдачи в периоды нагревания и охлаждения насадкн; Кц — коэффициент теплопередачи, отнесенный к одному циклу. Из уравнений (а) следует:

'н — бн = Яц/авхы 90 — 'о = WaoTo:

tH to flh + 90 = Яи. (1/«нТн + 1/а0То) Подставляя в последнее уравнение значение qa = Кц (*н — 'о)> находим:

Кц= [1/(1/а„гн+ 1/аоТо)] [1-(вн—в0)/(<н —<»)] (б)

Прн 0Н = 80 регенератор называется идеальным; в этом случае

< = '/(Wh + VVo) (в)

В реальных регенераторах используется не вся теплоаккумулнрую-щая способность регенераторной насадкн, а лишь доля ее, выражаемая коэф­фициентом использования насадкн регенератора т|. Последний зависит от значения критерия Fo = aT4/0,25s2 (где s — толщина элемента на­садкн) и может быть рассчитан по формуле:

т|= 1/(1+1/Fo) (г)

Легко видеть, что величина Ц тем ближе к единице, чем меньше толщина элемента насадкн и больше ее коэффициент температуропроводности (а). Заме­тим, что формула (г) справедлива прн синусоидальном изменении температур теплоносителей на протяжении каждого периода и является приближенной прн отклонении от этой закономерности изменения этих температур.

В практических расчетах оперируют обычно всей располагаемой поверх­ностью насадки F, учитывая при этом степень ее использования при оценке коэф­фициента теплопередачи Кц- Поэтому последний выражают произведением коэф­фициента теплопередачи идеального регенератора К'ц на коэффициент полез­ного действия поверхности теплообмена т)т, т. е. Кц = г|тКц. Величина % за­висит не только от значения г), но и от продолжительности цикла и составляю­щих его периодов, теплофнзических свойств и геометрических размеров элемен­тов насадки, а также коэффициента теплоотдачи от греющего потока к поверх-

н

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

ости насадки. Теоретический анализ в сочетании с опытными данными приво­дит к следующей зависимости:

<Пт = / [(aHTH/0,5cHpfls) rj] (д)

где % и рн — удельная теплоемкость и плотность греющего потока.

З ависимость (д) представлена на рис. VI1-27, откуда можно заимствовать численные значения г)т. Средние температуры ta и t0, фигурирующие в урав­нении {а), как показывает опыт, могут быть найдены исходя из ли­нейного изменения выходных тем­ператур теплоносителей (£, Q во времени за каждый период рабо­чего цикла регенератора. Следова­тельно, обозначив начальные тем­пературы теплоносителей через ^ я £, получим: f„ = + Qll;

Рис. V1I-27. К расчету регенераторов.

где WH — водяной эквивалент греющего потока за время одного цикла; АСр — среднелогарифмическая разность температур при противотоке или прямотоке, определяемая по ранее приведенным формулам для рекуперативного теплооб­менного аппарата.

Л. ПРОМЫШЛЕННЫЕ ТЕПЛОНОСИТЕЛИ

Наиболее экономичным и доступным теплоносителем (греющим агентом) является водяной пар, широко приме­няемый на химических предприятиях благодаря большому удель­ному, теплосодержанию (скрытая теплота испарения при нормаль­ном давлении составляет 2256,8 кДж/кг), постоянной температуре и высокому коэффициенту теплоотдачи при конденсации (~ 10 ООО Вт/(м*-К)). Нагревание водяным паром становится, од­нако, невыгодным для достижения температур выше 180—200 °С из-за высокого давления пара {температура насыщенного пара 200 °С соответствует давлению ~2МПа), так как это сопряжено с удорожанием аппарата. Эти же недостатки присущи воде при высоких температурах, уступающей к тому же водяному пару по значению коэффициента теплоотдачи.

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

^ Для обогрева аппаратов с рабочими температурами выше 200 °С применяются высококипящие органиче­ские и неорганические теплоносители в жидком и парообразном агрегатном состояниях при атмосферном или не­больших избыточных давлениях. К числу основных требований, предъявляемых к промышленным теплоносителям, относятся: воз­можно большая рабочая температура, большая объемная тепло­емкость, низкая вязкость, термическая и химическая стойкость, огне- и взрывобезопасность, нетоксичность, невысокая стоимость и низкие эксплуатационные затраты. Поиски веществ с таким сложным сочетанием физических и химических свойств обусло­вили появление многочисленного ряда теплоносителей, каждый из которых либо только частично удовлетворяет предъявляемым требованиям, либо в полной мере, но в ограниченном диапазоне рабочих условий. В связи с этим выбор оптимального теплоноси­теля в каждом конкретном случае представляет важную практи­ческую задачу.

Чрезвычайно разнообразные свойства и большое число пред­ложенных высококипящих жидких теплоносителей вызывают не­обходимость их классификации. Если базировать последнюю на принципе термодинамического подобия (молекулярная структура, тип слабейшей химической связи и критический коэффициент одинаковы), то можно разделить известные высококипящие тепло­носители на три основные группы: 1) органические (с остаточной связью); 2) ионные (с ионной связью); 3) жидкометаллические.

Органические теплоносители. Эта группа насчитывает большое число циклических, ациклических и смешанных соединений с тем­пературами кипения (при нормальном давлении) до 400 °С. Они применяются в виде индивидуальных веществ, а также в виде би­нарных и многокомпонентных смесей — эвтектических и неэвтек­тических. Свойства смесей данной группы, как и двух других, подробно освещены в специальной литературе (см. список литера­туры в конце книги). В табл. VII-1 приведены основные свойства теплоносителей, представляющих наибольший практический ин­терес для химической технологии.

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

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

п

Основные физические свойства ряда органических высококипящих теплоносителей при нормальном давлении

ТАБЛИЦА VIM

Температура, °С

Удельная теплоем-

Плот-

Энталь-

Теплоноситель

ность при 20 "С,

плавле-

кипения

кость.

пия жид­кости,

кг/м'

ния

(кг-К)

кДж/кг

Этиленгликоль

1113

—15,6

197,3

2,35

812,2

Глицерин

1260

—19

290,0

2,345

774,6

Нафталин

1150

80,2

218,0

2,010

340,8

Дифениловый эфир

1080

27,0

258,5

1,671

275,1

Дифенил бензол

1019

56,3

337,5

1,834

264,0

Дифенилметан

1006

25,4

264,3

1,733

376,8

Дитолилметан

982,6

—33,0

296,0

1,553

334,9

Дикумилметан

947,5

—23,0

335,0

1,436

332,9

Диксилилметан

989,9

—42,0

324,0

1,490

330,8

Тетраизопропилдифенилме-

927,8

—8,0

384,0

1,331

287,2

танМоноизопропилдифенил

979,0

—47,0

290,0

1,410

431,2

Изопропилтерфенил

1013,0

—15,0

346,0

1,158

Эвтектические смеси

26,5% дифенила+73,5%

1060,0

12,3

258,0

1,591

дифенилового эфира

28% дифеиила +72%

1007,0

12,3

263,0

1,400

дифенилметана

0,3% дифенила+11,7%

760

60,0

371,0

2,512

ортотерфенила

(316 °С)

(316 °С)

+60% ж-терфенила+

+28% л-терфенила

Неэвтектические смеси

34,4% дифенила+

1070,0

-7,0

290,0

1,616

+48,9% о-терфени-

ла + 11,5% л-терфе-

1006,0

—10,0

240,0

1,884

нила + 5,2% фенан-

трена + метилнафта-

лины

метилнафталины +10%

976,0

—10,0

240,0

1,842

парафинов

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

Органические теплоносители, как правило, горючи и взрыво­опасны, но малотоксичны и химически не агрессивны по отноше­нию к обычным конструкционным материалам (углеродистые и другие стали). Исключение составляют хлоропроизводные соеди­нения, которые отличаются высокой токсичностью и коррозион­ным действием на цветные металлы (медь, латунь, бронзы, алю­миний).

К числу факторов, определяющих в каждом случае выбор теп­лоносителя, относятся: требуемая рабочая температура, плот­ность, вязкость, удельная теплоемкость и коэффициент теплопро­водности. Верхний предел рабочей температуры ограничен нача­лом разложения теплоносителя, а нижний предел — его вяз­костью, возрастающей с понижением температуры и практически неприемлемой при кинематической вязкости выше 4'Ю-4 м2/с из-за большого расхода энергии на циркуляцию. С увеличением удельной объемной теплоемкости (рс) теплоносителя уменьшается его расход, необходимый для переноса требуемого количества тепла при заданном перепаде температур, и, следовательно, рас­ход энергии на его циркуляцию. Напомним также, что с уменьше­нием вязкости и увеличением теплопроводности теплоносителя возрастает его коэффициент теплоотдачи. Легко видеть, что наи­более выгодным является тот теплоноситель, который обеспечи­вает перенос требуемого количества тепла при минимальном рас­ходе энергии, наибольшем коэффициенте теплоотдачи и наимень­шем термическом разложении.

Процесс термического разложения органических теплоносите­лей достаточно изучен и освещен в литературе. Скорость разло­жения увеличивается экспоненциально с ростом температуры и в среднем примерно удваивается на каждые 10°. Допустимая ско­рость разложения измеряется несколькими процентами в год.

Среди применяемых в промышленности органических тепло­носителей первое место (около 40 % установок) занимает эвтек­тическая смесь дифенила и дифенилового эфира («даутерм»), используемая как в жидком (до температуры 330 °С), так и в парообразном состояниях (до температуры 400 °С). При использовании жидкой смеси с температурой выше 257 РС требуется ее подпрессовка сжатым азотом. Степень разложения данной смеси при температурах до 400 °С невелика; она несколько темнеет, но без существенного изменения своих физических свойств. При более высоких температурах происходит полимери­зация; образующиеся вещества растворимы, и смесь легко реге­нерируется путем простой перегонки. Допустимое накопление продуктов полимеризации равно 10—15%. В этом случае срок службы теплоносителя составляет 45—60 месяцев при рабочей температуре 350 °С и лишь 1,5—2 месяца при температуре 410 °С.

Самыми дешевыми, но и наименее термически устойчивыми теплоносителями являются минеральные масла. При использовании масел необходимо учесть, что вблизи температуры вспышки, которая обычно ниже 300 °С, наблюдается их термиче­ское разложение, сопровождающееся отложением кокса на по­верхности нагрева и выделением газообразных веществ. Послед­ние образуют с кислородом воздуха взрывоопасные смеси, а отло­жение кокса приводит к уменьшению коэффициента теплопере-

Основные физические свойства ряда высококипящих ионных теплоносителей при нормальном давлении

Температура, °С

135,9 230,0 263,3

200,0

250,0 254,0


Теплоноситель


Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]