CC BY
38
ВЕСТНИК ПНИПУ
2019 Химическая технология и биотехнология № 3
DOI: 10.15593/2224-9400/2019.3.08 УДК 66-97
Г.И. Ефремов1, В.И. Мартыненко2
1 Московский финансово-юридический университет (МФЮА),
Москва, Россия 2 ООО «Фармсистемы», Москва, Россия
УТОЧНЕННЫЙ РАСЧЕТ РЕКУПЕРАТИВНЫХ ТЕПЛООБМЕННЫХ АППАРАТОВ
Рассмотрен уточненный расчет при моделировании работы рекуперативных теплообменников с использованием сплайн-интерполяции. Рассмотрен как прямой (конструктивный), так и обратный (проверочный) расчеты теплообменников. В расчетах использовали усредненные по длине аппаратов значения температур и ряда параметров теплофизических свойств воды: среднюю теплоемкость, среднюю плотность, вязкость и др., в том числе и коэффициенты теплопередачи.
Выбрана противоточная, как наиболее эффективная, схема движения теплоносителей. Она предпочтительнее прямоточной, так как позволяет получить больший средний температурный напор и, тем самым, уменьшить требуемую площадь поверхности теплообмена. Циклы компьютерных расчетов повторяются до тех пор, пока не будет получена требуемая точность сходимости результатов.
Расчеты усложняются еще тем, что на каждой стадии цикла при изменении температуры надо вновь рассчитывать значения теплофизических свойств теплоносителей. При выполнении компьютерного расчета теплообменников предложено вносить на каждом цикле итерационных вычислений корректировку значений теп-лофизических свойств теплоносителей методом сплайн-интерполяции в среде Mathcad по табличным данным зависимости критерия Прандтля для обрабатываемой среды от температуры. В работе приведен конкретный пример расчета для водо-водяного рекуперативного теплообменника в среде Mathcad.
Ключевые слова: моделирование теплообмена, рекуперативный теплообмен, уточненный итерационный расчет, сплайн-интерполяция в среде Mathcad.
G.I. Efremov1, V.I. Martynenko2
1 Moscow University of Finance and law (MFUA), Moscow, Russian Federation 2 LLC "Pharmsystems", Moscow, Russian Federation
REFINED CALCULATION OF RECUPERATIVE HEAT EXCHANGERS
The refined calculation in the simulation of recuperative heat exchangers using spline interpolation is considered. Both direct (constructive) and reverse (verification) calculations of heat exchangers are considered. The calculations used averaged over the
length of the apparatus temperature values and a number of parameters of the thermal properties of water - the average heat capacity, average density, viscosity, etc., including heat transfer coefficients.
The countercurrent, as the most effective, scheme of movement of heat fluids is chosen. It is preferable to direct flow, as it allows to obtain a larger average temperature head and thus reduce the required heat exchange surface area. Computer calculation cycles are repeated until the required precision of convergence of the results is obtained.
Calculations are complicated by the fact that at each stage of the cycle, when the temperature changes, the values of the thermophysical properties of the heat carriers must be calculated again. For performing the computer-aided analysis of heat exchangers is proposed to be paid on each cycle of the iterative calculation, the correction values of thermophysical properties of fluids by the method of spline interpolation in Mathcad using tabular data dependency of the Prandtl number for the treated medium temperature. The paper presents a specific example of calculation for water-water recuperative heat exchanger in Mathcad.
Keywords: modeling of heat exchange, recuperative heat exchange, refined iterative calculation, spline interpolation in Mathcad.
При проектировании теплообменных аппаратов выполняют как прямой (конструктивный), так и обратный (проверочный) расчеты. При конструктивном расчете, как правило, известны: расход нагреваемого теплоносителя G2, а также начальная и конечная температуры обоих теплоносителей. В процессе расчета выбирают тип теплообменника и определяют его тепловую мощность Q, расход горячего теплоносителя Gi и площадь поверхности теплообмена F, по которой устанавливают габаритные размеры аппарата.
При проверочном расчете теплообменных аппаратов известны тип и размеры теплообменника, а следовательно, и поверхность теплообмена, расходы теплоносителей, их начальная температура, удельные теплоемкости и коэффициент теплопередачи. И нужно определить тепловую мощность аппарата и конечную температуру теплоносителей, т. е. проверить применимость имеющегося теплообменника для конкретных условий [1-8].
Для стационарного процесса нагрева (охлаждения) различных сред рекомендуется [3] использовать в расчетах усредненные по длине аппаратов значения температур и ряда параметров, в том числе и коэффициентов теплопередачи. Так, по средней температуре теплоносителей t1 и t2 по таблицам теплофизических свойств воды находят среднюю теплоемкость, среднюю плотность, вязкость и другие параметры теплоносителей.
Выбирают обычно противоточную, как наиболее эффективную, схему движения теплоносителей. Она предпочтительнее прямоточной, так как позволяет получить больший средний температурный напор и, тем самым, уменьшить требуемую площадь поверхности теплообмена.
На начальном этапе ни коэффициент теплопередачи K в теплообменнике, ни площадь поверхности теплообмена F не известны. Поэтому величиной коэффициента теплопередачи приходится задаваться. В дальнейшем расчете, когда уже известна геометрия проточной части теплообменника, выполняется уточненный расчет коэффициента теплопередачи K, в результате чего уточняется и величина площади поверхности теплообмена F. Для определения конструкции теплообменника задаются рекомендуемой скоростью воды в трубках. Так, для напорного движения жидкостей рекомендуются [3] скорости 0,5-2,5 м/с.
По каталогу размеров труб выбирают их размер dнx5, т.е. наружный диаметр и толщину стенки 5, мм. Требуемое число трубок п можно определить из уравнения неразрывности. Приняв расчетное число трубок в пучке п, уточняют значение скорости нагреваемого теплоносителя в трубках, рассчитывают площадь поверхности теплообмена К
При приближенном расчете находят среднюю температуру горячего и нагреваемого теплоносителя, по средней температуре теплоносителей по таблицам теплофизических свойств воды находят средние теплоемкости и передаваемый в теплообменнике тепловой поток. Вычисляют средний логарифмический температурный напор и по ориентировочному значению коэффициента теплопередачи К в конструируемом рекуператоре и находят оценочное значение площади поверхности теплообмена F. Массовый расход греющего теплоносителя определяют, исходя из уравнения теплового баланса (предварительно задают скорость воды в трубках) и рассчитывают геометрию трубного пучка теплообменника.
Уточненный расчет рекуперативного теплообменника включает расчет коэффициента теплопередачи по рекомендациям [3] как и для плоской стенки. Для расчета коэффициента теплоотдачи от греющего теплоносителя к стенке и коэффициента теплоотдачи от стенки к нагреваемому теплоносителю обычно применяют [3] формулу М.А. Ми-хеева для турбулентного режима потока:
Ми = 0,021- Ие0,8- РГ
0,8 тэ 0,43
\ 0,25 рг ^
РГ
Vх ст У
(1)
Поскольку температуры поверхности теплообмена 01 и 02 заранее неизвестны, коэффициенты теплоотдачи а1 и а2 приходится рассчитывать методом последовательных приближений. Задаваясь температурой поверхности теплообмена, определяем численные значения коэффициентов теплоотдачи, зная которые можно уточнить температуру поверхности теплообмена F по общему уравнению теплообмена
Q = ^(11 - 12) = а1 F(tl - 01) = а2 F(02 - г 2) = ^ F(01 - 0 2). (2)
о
Потом вновь повторяют расчет коэффициентов теплоотдачи. Цикл расчетов повторяют до тех пор, пока не будет получена требуемая точность сходимости результатов. Расчеты усложняются еще тем, что на каждой стадии цикла при изменении температуры надо вновь рассчитывать значения теплофизических свойств теплоносителей. Они выражаются значением критерия Рг от температуры. Для ускорения вычислений и повышения точности лучше выполнять расчеты на ПК [7, 9].
Нами рекомендовано вносить на каждом цикле вычислений корректировку значений теплофизических свойств теплоносителей методом сплайн-интерполяции на основе введенных табличных данных зависимости критерия Прандтля для воды от температуры. Ниже приведен такой пример расчета для водо-водяного рекуперативного теплообменника в среде Ма&саё [10-15].
На рис. 1 показан выполненный в среде МаЛсаё листинг ввода табличных теплофизических свойств теплоносителей (горячая и холодная вода) по средней температуре и расчет средней логарифмической разности температур для конкретного примера. Индекс 1 относится к горячему теплоносителю, а 2 - к холодному. Параметры теплоносителей можно определить по таблице методом интерполяции [7] .
На рис. 2 показано продолжение листинга, согласно изложенной выше методике - расчет теплового потока, ориентировочной поверхности теплообмена и числа трубок теплообменника.
На рис. 3 показан дальнейший расчет размеров сечения кожухо-трубного теплообменника и критериев Ие для теплоносителей в среде МаЛсаё. Минимальный зазор между крайними трубками и корпусом теплообменника обычно принимается равным к > 5 мм. Принимаем к = 5 мм.
Уточненный расчет теплообменника Заданы начальная и конечная температуры горячего теплоносителя - 90 и 78 оС Средняя температура горячего теплоносителя Н := 05-(90 + 78) 11 = 34 оС Параметры горячего теплоносителя
теплоемкость разность тем-р плотность теплопроводность вязкость число Рг
с 1 := 4190 ЛИ := 90 - 75 р1 := Ш2 11 := 0..676 VI := 0Л49 10~ 6 Рг1 2.1 Пример линейной интерполяции по табл. 1 для липкости при 84 оС
г := \ 30 : р := : ? 13 : 1т1ыр(1.р.Ш = 9692 \_90у ^965 ■
Заданы начальная и конечная температуры нагреваемого теплоносителя -20 и 37 оС Средняя температура нагреваемого теплоносителя: 12 := 0_5-(37 +20) Й = 25.5 оС
Параметры нагреваемого теплоносителя аналогично:
с2 := 4175 АО := 37 - 20 р2 := 996.03 12 := 0.615 \'2 := 0.335-Ю~ 6 Рг2 := 5.66
Большая и меньшая разность температур Л1б := 73- 20 Айн := 90- оС
Средняя логарифмическая разность температур: ^^
Д1 := -г-г- А1 = 55.452 оС
_ЧАЬпу_
Рис. 1. Расчет теплофизических свойств теплоносителей и средней логарифмической разности температур
Теплофизические свойства воды
°С Р ■ 10-5, Па кг/м3 й, кДж/кг кДж/(кгград) Ы0-2, Вт/(мград) а ■ 106, м2/с V ■ 10-6, м2/с Рг
0 1,013 999,9 0,00 4,212 55,1 13,1 1,789 13,67
10 1,013 999,7 42,04 4,191 57,4 13,7 1,306 9,52
20 1,013 998,2 83,91 4,183 59,9 14,3 1,006 7,02
30 1,013 995,7 125,7 4,174 61,8 14,9 0,805 5,42
40 1,013 992,2 167,5 4,174 63,5 15,3 0,659 4,31
50 1,013 988,1 209,3 4,174 64,8 15,7 0,556 3,54
60 1,013 983,2 251,1 4,179 65,9 16,0 0,478 2,98
70 1,013 977,8 293,0 4,187 66,8 16,3 0,415 2,55
80 1,013 971,8 335,0 4,195 67,4 16,6 0,365 2,21
90 1,013 965,3 377,0 4,208 68,0 16,8 0,326 1,95
100 1,013 958,4 419,1 4,220 68,3 16,9 0,295 1,75
Задан массовый расход нагреваемого теплоносителя: 02 := 1.3 кг/с
Тепловой поток нагреваемого теплоносителя:
Л
(3 := 02-с2-Ш <3 = 5.227 х 10 Приближенное значение коэф. теплообмена:
Вт
К := 1700
Поверхность теплообмена в первом приближении:
м2К
Р - —— Р = 0 979 Принимаем скорость холодного потока: та-2 = 1 5 — КМ с
Расход горячего потока: 01 :
с1-ЛИ
01 = 1.335 кг/с
Размеры трубок.мм: с1н:=1б 6 := 1 с! := с!н - 2 6 ¿=14 с1с := 0.5 (с)н + с!) с1с = 15 Число трубок: 4-в2
лр2 «'2 ^10 ^
п = 5.652 Примем п:=6
Уточненная скорость холодного потока: :=
4 02
л р2 п й 10
^2= 1.413 -2 .„-6 с
Рис. 2. Расчет теплового потока поверхности теплообмена и числа трубок
Длина трубок теплообменника Р
Ь:=
-3
Ь = 3.461 м Примем Ь := 3 м Термосопротивление стенки:
лп с1с-10
Размеры сечения теплообменника: Шаг: й := 1.25-йн з = 20 к:=5 Бе := 2 е + йн + 2-к Бе = 66 мм
0.001 -б
Ег ---Ег = 9.346 х 10
107
Сечение горячего потока: :=- '.Ов^ - п с1н^' ю ^ f = 2.215x10 ^ м^
4
Скорость горячего потока: те-1 = «-1 = 0 355 — Эквивалентный диаметр:
с
- з
Ке2 := —- Йе2 = 2.3692 х 10^ Режим турбулентный
^ := 0.0174 м
Ее1 := ■
v2 те1! ■ с!е
VI
БЫ = 4.262 х 10 Режим турбулентный
Рис. 3. Расчет размеров сечения теплообменника и критериев Ие
теплоносителей
На рис. 4 показаны введенные табличные значения критерия РГ при различных температурах [6], график, построенный по этим значениям, и расчет по ним сплайн-интерполяция в среде МаЛсаё.
Рис. 4. Сплайн-интерполяция в среде Mathcad для зависимости Pr
от температуры
На рис. 5 показано решение системы нелинейных уравнений (1) и (2) для обоих теплоносителей методом итераций в среде Mathcad. Заданы для начала итераций ориентировочные значения температур на стенке трубы 91 и 02 и коэффициенты теплоотдачи а1 и а2. Запуск итераций производится по кодовому слову Given (дано).
Рис. 5. Решение в МаШсаё системы нелинейных уравнений методом итераций
В результате расчета на ПК для заданных параметров потоков получены значения коэффициентов теплоотдачи ai = 5025 Вт/м-К и a2 = 6944 Вт/м2-К, значения температур стенки труб 91 = 52,039 °С и 02 = 51,183 °С, поверхность теплообмена F = 1,007 м2 и коэффициент теплопередачи К = 2838 Вт/м2-К.
Список литературы
1. Боровко В.М., Калютик А.А., Сергеев В.В. Теплотехническое оборудование: учеб. - 2-е изд., испр. - М.: Академия, 2013. - 192 с.
2. Васильченко Ю.В. Теплогенерирующие установки. - Белгород: Изд-во БГТУ им. Шухова, 2008. - 162 с.
3. Вафин Д.Б. Источники производства теплоты / Нижнекам. хим.-технол. ин-т (филиал) ФГБОУ ВПО «КНИТУ». - Нижнекамск, 2014. - 242 с.
4. Дегтяренко А.В. Теплоснабжение. - Томск: Изд-во Том. гос. архит.-строит. ун-та, 2010. - 185 с.
5. Евсеева Т.А., Ластовец Н.В. Отопление: конспект лекций / Харьк. нац. акад. гор. хоз-ва. - Харьков, 2012. - 84 с.
6. Касаткин А.Г. Основные процесс и аппараты химической технологии. - М.: Химия, 1971. - 784 с.
7. Павлов К.Ф., Романков П.Г., Носков А.А. Примеры и задачи по курсу процессов и аппаратов химической технологии. - Л.: Химия, 1976. - 550 с.
8. Основные процессы и аппараты химической технологии: пособие по проектированию / под ред. Ю.И. Дытнерского. - М.: Химия, 1983. - 272 с.
9. Ефремов Г.И. Моделирование химико-технологических процессов: учеб. - М.: ИНФРА-М, 2016. - 254 с.
10. Амосова О.А., Вестфальский А.Е. Применение пакета Mathcad к решению вычислительных задач: метод. пособие. - М.: Изд. дом МЭИ, 2007. - 30 с.
11. Бундаев В.В. Решение задач линейной оптимизации с использованием Mathcad и ExceL: метод. пособие и контр. задания. - Улан-Удэ, 2006. -31 с.
12. Гурьяшова Р.Н., Шеянов А.В. Информатика. Пакет Mathcad: учеб. пособие. - Н. Новгород, 2005. - 140 с.
13. Mathcad 6.0 Plus. Финансовые, инженерные и научные расчеты в среде Windows 95. - М.: Филинъ, 1996. - 712 с.
14. Очков В.Ф. Mathcad 14 для студентов, инженеров и конструкторов. - СПб.: БХВ- Петербург, 2007. - 368 с.
15. Дьяконов В .П. Справочник по Mathcad PLUS 6.0 PRO. - М.: СК Пресс, 1997. - 336 с.
References
1. Borovko V.M., Kaliutik A.A., Sergeev V.V. Teplotekhnicheskoe oborudovanie [Heat engineering equipment]. 2nd ed. Moscow, Izdatel'skii tsentr «Akademiia», 2013, 192 p.
2. Vasil'chenko Iu.V. Teplogeneriruiushchie ustanovki [Heat Generating plants]. Belgorod, Izd-vo BGTU im. Shukhova, 2008, 162 p.
3. Vafin D.B. Istochniki proizvodstva teploty [Sources of heat production]. Nizhnekamsk, Nizhnekamskii khimiko-tekhnologicheskii institut, 2014, 242 p.
4. Degtiarenko A.V. Teplosnabzhenie [Heat]. Tomsk, Izd-vo Tom. gos. arkhit.-stroit. un-ta, 2010, 185 p.
5. Evseeva T. A., Lastovets N. V. Konspekt lektsii po kursu «Otoplenie» [Lecture Notes on the course "heating"]. Kharkov, KhNAGKh, 2012, 84 p.
6. Kasatkin A.G. Osnovnye protsess i apparaty khimicheskoi tekhnologii [Basic process and apparatus of chemical technology]. Moscow, Khimiia, I971, 784 p.
7. Pavlov K. F., Romankov P. G., Noskov A. A. Primery i zadachi po kursu prot-sessov i apparatov khimicheskoi tekhnologii [Examples and tasks on the course of processes and apparatus of chemical technology]. Leningrad, Khimiia, 1976, 550 p.
8. Dytnerskii Iu.I.Osnovnye protsessy i apparaty khimicheskoi tekhnologii [Basic processes and apparatus of chemical technology]. Moscow, Khimiia, 1983, 272 p.
9. Efremov G.I. Modelirovanie khimiko-tekhnologicheskikh protsessov [Modeling of chemical and technological processes]. Moscow, NITs INFRA-M, 2016, 254 p.
10. Amosova O.A., Vestfal'skii A.E. Primenenie paketa Mathcad k resheniiu vychis-litel'nykh zadach [The Use of Mathcad to the solution of computational problems]. Moscow, Izdatel'skii dom MEI, 2007, 30 p.
11. Bundaev V.V. Reshenie zadach lineinoi optimizatsii s ispol'zlvaniem MATHCAD i EXCEL [the Solution of problems of linear optimization ispolzuyet MATHCAD and EXCEL]. Ulan-Ude, 2006, 31 p.
12. Gur'iashova R.N., Sheianov A.V. Informatika. Paket Mathcad [Computer Science. Mathcad]. Nizhny Novgorod, Izd-vo FGOU VPO VGAVT, 2005, 140 p.
13. Mathcad 6.0 Plus. Finansovye, inzhenernye i nauchnye raschety v srede Windows 95 [Mathcad 6.0 Plus. Financial, engineering and scientific calculations in Windows 95]. Moscow, Inform.-izdat. dom «Filin"», 1996, 712 p.
14. Ochkov V.F. Mathcad 14 dlia studentov, inzhenerov i konstruktorov BKhV [Mathcad 14 for students, engineers and designers BHV]. Saint Petersburg, 2007, 368 p.
15. D'iakonov V. P. Spravochnik po Mathcad PLUS 6.0 PRO [Handbook of Mathcad PLUS 6.0 PRO]. Mosscow, «SK Press», 1997, 336 p.
Получено 31.07.2019
Об авторах
Ефремов Герман Иванович (Москва, Россия) - доктор технических наук, профессор кафедры строительства Московского финансово-юридического университета (МФЮА) (117342, г. Москва, ул. Введенского, д. 1А; e-mail: [email protected]).
Мартыненко Валерий Иванович (Москва, Россия) - директор по производству ООО «Фармсистемы» (142715, Московская обл., Ленинский р-н, с. Беседы, Промышленный проезд, д. 1, стр. 1, e-mail: [email protected]).
About the authors
German I. Efremov (Moscow, Russian Federation) - Doctor of Technical Sciences, Professor, Department of construction of the Moscow University of Finance and Law (MFUA) (1A, Vvedensky st., Moscow, 117342, e-mail: efremov_german@ mail.ru).
Valery I. Martynenko (Moscow, Russian Federation) - Director of production of LLC "Pharmsystems" (1, b. 1, Promyshlennyy proezd, s. Besedi, Leninskiy r-n, Moscow region, 142715, e-mail: [email protected]).
