Belgorod, Belgorod, Russian Federation
GRNTI 67.53 Инженерное обеспечение объектов строительства
BBK 313 Теплоэнергетика. Теплотехника
The use and prevalence of heat exchangers in Russian heat supply systems are considered. Attention is paid to the improvement of serial heat exchangers with smooth tubes – the increasing of heat transfer coefficient. One of the ways to increase the heat transfer coefficient is considered: it is the turbuliza-tion of the fluid flow on the heat exchange surface. The original design of the heat exchange surface for shell and tube devices of heat supply systems is presented. The dynamics of the heated fluid in the annular space of a shell and tube heat exchanger when flowing around the heat exchange surface with a modified geometry is studied (RF Patent 149737). A feature of the dynamics is a circular edge (element of the surface of heat exchange), which contributes to the creation of turbulence in the flow of the heated liquid on the plate and on the surface of the next edge. Emphasis is placed on heat ex-change processes between the solid surface of the edges and the heated fluid. For a circular cross sec-tion, the equation of thermal conductivity in cylindrical coordinates is compiled, taking into account the stationary heat exchange process, with an internal source of thermal energy. Solution of equation makes possible to determine the change in temperature on the surface and the average temperature of the edge. This value allows determining the Prandtl number to calculate the heat transfer coefficient.
heat exchange surface, turbulization, temperature, heat transfer coefficient.
Введение. Теплообменные аппараты поверхностного типа, в которых теплота передается от греющего к нагреваемому теплоносителю через разделяющую их твердую теплопроводную стенку, широко применяются в стационарных и транспортных энергетических установках, в тепло- и электроэнергетике, а также в различного рода теплотехнологических установках. Эти аппараты составляют значительную долю объема производства энергетического машиностроения. В системах теплоснабжения и теплофикации такими аппаратами являются сетевые и регенеративные подогреватели воды, охладители конденсата пара, охладители различных технологических продуктов, конденсаторы и переохладители тепловых насосов [1, 2].
В 2017 календарном году в Российской Федерации было отпущено потребителям 4,94 ·108 Гкал тепловой энергии. Это на 13 % больше, чем за предыдущий год [3].
Применение серийных теплообменных аппаратов, в частности с гладкими трубками (ГОСТ ГОСТ 27590), во многих случаях приводит к тому, что такие аппараты имеют, как правило, невысокий коэффициент теплопередачи. Это приводит к значительной площади поверхности теплообмена, повышенным металлоемкости и габаритам [4]. В термодинамических циклах, процессах производства при преобразовании энергии удельный вес теплообменных аппаратов как оборудования по сравнению с общим паротурбинной установки достигает 30–35 % [5, 6].
Таким образом, из-за широкого распространения кожухотрубных теплообменных аппаратов актуален вопрос об интенсификации тепловых процессов.
Методика. В БГТУ им. В.Г. Шухова разработана оригинальная конструкция кожухотрубного теплообменного аппарата. Особенностью аппарата являются теплообменные трубки 1, оснащенные пластинами 2 с расположенными на них ребрами 3 цилиндрической формы [6]. Это позволяет увеличить теплообменную поверхность и создать дополнительную турбулизацию потока нагреваемой жидкости при обтекании этих ребер [7] (рис. 1).
В классических работах Жукаускаса А.А., Шлихтинга Г., Хинце И.О. [8 - 10] установлено, что при обтекании потоком жидкости цилиндра, начиная с числа Re =60, за кормовой частью образуется зона завихрений, а при значении Re >5000 и выше за кормовой частью – полное турбулентное перемешивание. Это приводит к уменьшению толщины ламинарного подслоя жидкости у пластины, а, следовательно, увеличению передачи теплоты через этот слой.
|
|
Рис. 1. Элемент теплообменной поверхности оригинального кожухотрубного теплообменного аппарата: 1 – трубка, 2 – пластина, 3 – ребра цилиндрической формы
|
На основании СП 41-01-95 «Проектирование тепловых пунктов» коэффициент теплопередачи К, Вт/(м2К) для аппаратов с гладкими трубками рассчитывается:
(1)
где Ψ – коэффициент эффективности теплообмена, β – коэффициент, учитывающий загрязнение поверхности труб [11], – коэффициент теплоотдачи от внешней поверхности нагретой теплообменной трубки к нагреваемой жидкости межтрубного пространства, Вт/(м² К), δ – толщина стенки теплообменной трубки, м, λ – коэффициент теплопроводности материала трубки, Вт/(м·К),
– коэффициент теплоотдачи от горячей жидкости, протекающей в трубке, к внутренней поверхности этой трубки, Вт/(м² К).
Так как в нашем случае изменена геометрия внешней поверхности трубки, то необходимо определить коэффициент .
Основная часть. Для этого рассмотрим обтекание потоком нагреваемой жидкости теплообменной поверхности, состоящей из пластины и ребер цилиндрической формы (рис. 2). За первым ребром образуется зона повышенной турбулизации потока (с многочисленными вихрями), которая распространяется на определенное расстояние. Следующее ребро установлено на таком расстоянии, при котором лобовая часть этого ребра находится еще интенсивном вихревом потоке. Таким образом, по поверхности следующего ребра будут так же распространяться многочисленные вихри жидкости. Турбулизированный поток жидкости будет способствовать уменьшению толщины ламинарного подслоя жидкости у поверхности, а,
следовательно, – увеличению передачи теплоты (повышению коэффициента теплоотдачи) через этот слой [3].
Рис. 2. Динамика жидкости приобтекании поверхности теплообмена
1 – ребро; 2 – зона турбулизации; 3 – пристенная область; 4 – промежуточный слой, 5 – вязкий подслой;
6 - пластина, 7 – зоны вихрей, 8 – внешняя граница вихрей
Для нахождения изменения температуры внутри цилиндрического тела радиусом и длиной L с вынужденным источником тепла
выполненного в виде тонкой плоской пластинки, расположенной вдоль оси цилиндра (рис. 2.3) воспользуемся стационарным уравнением теплопроводности в цилиндрических координатах, учитывая, что в силу аксиальной симметрии функция распределения температуры внутри цилиндра зависит от координаты φ [12];
(2)
где λ – коэффициент теплопроводности цилиндрического тела, Вт/(м·К),
Рис. 3 Расчетная схема ребра цилиндрической формы для расчета температуры внутри цилиндрического тела с источником теплоты в виде плоской пластины
1 – вихри, 2 – пластина, 3 – ребро
В состоянии термодинамического равновесия количество тепла, выделенного в единицу времени в единице объема цилиндрического тела определяется соотношением:
(3)
где V – объем цилиндра, м3, – количество тепла, выделенного плоской пластиной, Вт.
Учитывая что на основании [12] количество тепла, выделенного пластиной определяется по формуле:
(4)
где – коэффициент теплопроводности пластины, Вт/(м К), S – площадь поверхности пластины, м2,
– градиент температуры по оси z, К/м.
Объем цилиндра м3, определяется по формуле:
(5)
где - радиус ребра, м, L – высота ребра, м (см. рис. 3). Тогда площадь поверхности цилиндра равна:
(6)
И формула (4) принимает вид:
(7)
Подстановка (7) в (2) приводит к следующему уравнению:
(8)
где введено следующее обозначение ( - коэффициент теплопроводности стали,
- коэффициент теплопроводности меди, Вт/(м К)):
(9)
Для нахождения решения (8) применим метод Фурье [13], согласно которому решение ищем в виде:
(10)
Из-за громоздкости математических выкладок запишем конечное уравнение изменения температуры , К, по поверхности ребра, которое имеет вид:
(11)
где Т0 – температура в точке соприкосновения пластины и трубки с координатой z=0, Т1 – температура на конце пластины с координатой z=L.
Таким образом, изменение температуры по поверхности ребра зависит от температур Т0, Т1, радиуса ребра R0 и высоты ребра L.
Нами получена объемная картина изменения температуры по поверхности ребра. График распределения температуры по поверхности ребра (11) представлен на рис. 4, где температура меняется от 350,3 К до 352,4 К; координата z (длина ребра) меняется от 0 до 4 мм, координата r (радиус) варьируется от 0 до 4 мм. График изменения температуры носит линейный характер, причем температура убывает с увеличением радиуса ребра r, что хорошо согласуется с физической интерпретацией процессов теплообмена и теплопроводности в поверхности кожухотрубного аппарата.
Полученная зависимость (11) используется для определения температуры в любой точке ребра. Это позволит определить среднюю температуру на поверхности ребра по формуле:
(12)
Рис. 4 Распределения температуры по
поверхности ребра
Выводы. Зависимость (12) позволит определить число Прандтля Pr, характеризующее отношение скоростных и температурных полей в потоке жидкости. Число Pr используется для расчета коэффициента теплоотдачи , который является важной величиной для определения коэффициента теплопередачи К, Вт/(м2 К).
Таким образом нами решена задача аналитического определения коэффициента , Вт/(м2 К), и коэффициента теплопередачи К, Вт/(м2 К) для измененной геометрии поверхности теплообмена (см. рис. 1) в интенсивном кожухотрубном теплообменном аппарате.
Источник финансирования. Программа развития опорного университета на базе БГТУ им. В.Г. Шухова.
1. Case V.M., London A.L. Compact heat exchangers [Kompaktnye teploobmenniki]. M.: State Energy Publishing House, 1962. 160 p. (rus)
2. Kasatkin K.A., Barochkin A.E., Zhukov V.P. Development of a mathematical model of multi-flow heat exchangers taking into account the phase transition in coolants [Razrabotka matematicheskoj modeli mnogopotochnyh teploobmennyh apparatov s uchetom fazovogo perekhoda v teplonositelyah]. Vestnik IPEU. 2018. No. 5. Pp. 61-66. (rus)
4. Kuschev LA, Nikulin N.Yu., Feoktistov A.Yu., Yakovlev E.A. Intensification of thermal processes in a shell-and-tube heat exchanger [Intensifikaciya teplovyh processov v kozhuhotrubnom teploobmennom apparate]. Scientific Bulletin of the Voronezh GASU. Construction and architecture. 2016. No. 3. Pp. 9-17. (rus)
5. Chicherin S.V. Analysis of the structure of the existing and future consumption of the centralized heat supply system [Analiz struktury sushchestvuyushchego i perspektivnogo potrebleniya sistemy centralizovannogo teplosnabzheniya]. Engineering Bulletin of the Don. 2017. No. 1. Pp. 91-102. (rus)
6. Malyshev A.A., Mamchenko V.O., Mizin V.M. Modern heat exchangers in low-temperature technology and processing industries [Sovremennye teploobmennye apparaty v nizkotemperaturnoj tekhniki i pererabatyvayushchih proizvodstvah]. Scientific journal NRU ITMO. Series "Refrigeration and air conditioning". 2015. No. 1. Pp. 38-44. (rus)
7. Nikulin N.Yu., Kuschev L.A., Suslov D.Yu. Shell-and-tube heat exchanger; applicant and patentee Pat. RF no. 149737, 2014. (rus)
8. Nikulin N.Yu., Kushchev LA, Semenok V.S. Development of a shell-and-tube heat exchanger with a modified geometric surface [Razrabotka kozhuhotrubnogo teploobmennogo apparata s izmenennoj geometricheskoj poverhnost'yu]. Systems Requirements: AdobeAcrobatReader. URL: http://www.scienceforum.ru/2015/pdf/17842.pdf (request date 18.03.19). (rus)
9. Kruglov G.A., Bakunin V.V., Andreeva M.V. Theoretical studies of the degree of interrelation of flow turbulence with heat transfer coefficient [Teoreticheskie issledovaniya stepeni vzaimosvyazi turbulizacii potoka s koehfficientom teplootdachi]. Herald KrasGASU. 2015. No. 6. Pp. 67-73. (rus)
10. Zhukauskas A.A. Convective transfer in heat exchangers [Konvektivnyj perenos v teploobmennikah]. M.: Science, 1982. 472 p. (rus)
11. Hinze J.O. An introduction to its mechanism and theory. London: McGraw-Hill book company, 1959. 680 p.
12. Schlichting G. Theory of the boundary layer [Teoriya pogranichnogo sloya]. M.: Science, 1974. 712 p. (rus)
13. Galkovsky V.A., Chupova M.V. Analysis of the reduction in heat transfer coefficient of heat exchangers due to surface contamination [Analiz snizheniya koehfficienta teploperedachi teploobmennyh apparatov vsledstvie zagryazneniya poverhnosti]. Systems. Requirements: AdobeAcrobatReader. URL: naukovedenie.ru/PDF/41TVN217.pdf (request date 18.03.19). (rus)
14. Nashchokin V.V. Technical thermodynamics [Tekhnicheskaya termodinamika]. M.: Higher School, 1995. 496 p. (rus)
15. Karpovich D.S., Susha O.N., Korovina N.P. Analytical and numerical methods for solving the heat equation [Analiticheskij i chislennyj metody resheniya uravneniya teploprovodnosti]. Proceedings of the Belarusian State Technical University. 2015. No. 6. Pp. 122-127. (rus)
