Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

COMPARATIVE ANALYSIS OF PARAMETERS COUNTERFLOW HEAT EXCHANGER, PREPARED ACCORDING TO METHOD OF COMPUTER FLUID DYNAMICS AND PRACTICAL DESIGN METHOD OF HEAT EXCHANGEAPPARATUS

Zharov A.V. 1 Pavlov A.A. 1 Vlasov V.V. 1 Favstov V.S. 1 Gorshkov R.V. 1
1 FGBOU VPO «The Yaroslavl state technical university»
Calculation of heat exchangers is a necessary element in the design of power equipment. The currently used methods of calculation of heat exchangers based on experimental data and theory of similarity in the form of criterial equations, generalizing this data. Application of technology of computer fluid dynamics modelling of processes in heat exchangers will allow to increase the reliability of the basic design parameters of heat exchange. The purpose of this work is to compare the results of calculation of the heat exchanger, using the technologies of the computer fluid dynamics and practical methodology for the design of heat exchangers. Considered hydrodynamic and heat-exchange processes in heat exchangers using the practical methodology for the design of heat exchangers and technology computational fluid dynamics. The results of calculations discovered the advantages and disadvantages of both methods and degree of reliability of the obtained results.
heat exchanger
computer fluid dynamics
heat transfer
coolant temperature
coolant velocity
the heat transfer coefficients
1. Spravochnik po teploobmennikam. V 2-kh t. T. 1 (Directory on heat exchangers. In 2 vol. Vol. 1). Per. s angl. pod red. O.G. Martynenko i dr. Moscow: Energoatomizdat, 1987. 560 p.
2. Spravochnik po teploobmennikam. V 2-kh t. T. 2 (Directory on heat exchangers. In 2 vol. Vol. 2). Per. s angl. pod red. O.G. Martynenko i dr. Moscow: Energoatomizdat, 1987. 352 p.
3. Spravochnik khimika. V 7-mi t. T. 5. 2-e izd., pererab. i dop. (Directory of the chemist. In 7 vol. Vol. 5. The 2nd edition, processed and added). Pod red. B.P. Nikolskogo. Leningrad: Chemistry, 1968. 976 p.
4. Khauzen Kh. Teploperedacha pri protivotoke, pryamotoke i perekrestnom toke. Per. s nem. (Heat transfer at a countercurrent, a direct flow and cross current. Translation from German). Мoscow: Energoizdat, 1981. 384 p.
5. Fletcher K. Vychislitelnyye metody v dinamike zhidkostey. V 2-kh t. T. 1. Per. s angl. (Computational techniques for fluid dynamics. In 2 vol. Vol. 1. Translation from English). Мoscow: World, 1991. 504 p.
6. Fletcher K. Vychislitelnyye metody v dinamike zhidkostey. V 2-kh t. T. 2. Per. s angl. (Computational techniques for fluid dynamics. In 2 vol. Vol. 2. Translation from English). Мoscow: World, 1991. 552 p.

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

Процесс теплообмена в теплообменниках весьма сложен для математического моделирования и расчета. Широко применяемые до настоящего времени методики расчета теплообменников [1, 2, 3, 4] основаны на использовании средних по всему объему параметров: температурного напора, условного среднего коэффициента теплоотдачи, эффективной площади контакта теплоносителей и др. Кроме того, данные методики базируются целиком на опытных данных и теории подобия в форме критериальных уравнений, обобщающих эти данные. В этих условиях точность расчета и проектирования теплообменных аппаратов целиком и полностью базируется на производственном опыте создания теплообменников, который является весьма специфическим и для переноса на новые объекты требует экспериментальной проверки.

С активным развитием технологий математического моделирования, в основе которых лежат фундаментальные уравнения гидрогазодинамики [5, 6], стало возможно расчетным путем определять параметры потока различных сред (давления, скорости) с определением характеристик конвективного теплообмена между твердой поверхностью теплопередающего элемента и текущего теплоносителя. Вышеуказанное математическое моделирование, в основе которого лежат фундаментальные уравнения гидрогазодинамики, осуществляется с применением численных методов и мощных вычислительных средств. В профессиональном сообществе подобные технологии получили название – «компьютерная гидрогазодинамика» (международная аббревиатура CFD). Применение технологии компьютерной гидрогазодинамики для моделирования процессов в теплообменниках позволит повысить степень достоверности основных расчетных параметров теплообмена и выбрать рациональные конструкции и режимы их эксплуатации. Вместе с тем в практике создания теплообменников достаточно редко используются методы их расчета с применением технологии компьютерной гидрогазодинамики. По нашему мнению, отсутствуют работы по критическому сравнению и анализу возможностей технологии компьютерной гидрогазодинамики применительно к расчету теплообменников и широко применяемой в настоящее время практической методики проектирования теплообменных аппаратов. В связи с этим целью данной работы является проведение сравнительного анализа результатов расчета теплообменника, полученных с помощью технологии компьютерной гидрогазодинамики и практической методики проектирования теплообменных аппаратов.

В качестве объекта исследования выбран рекуперативный теплообменник типа «труба в трубе» с противоточной схемой движения теплоносителей, у которого в качестве материала труб используется нержавеющая сталь. Схема рассматриваемого теплообменника с конструктивными размерами представлена на рис. 1. В качестве граничных условий для расчета процессов в теплообменнике двумя вышеуказанными методами заданы физические свойства холодного и горячего теплоносителей, в качестве которых используется вода (холодная и горячая вода соответственно), их массовые расходы, а также их температуры на входе в соответствующие его полости. Граничные условия для расчета теплообменника представлены в табл. 1.

pic_9.tif

Рис. 1. Расчетная схема теплообменника

Таблица 1

Конструктивные параметры теплообменника и теплофизические параметры теплоносителей

Параметр

Значение

Температура холодной воды на входе, °С

60

Температура горячей воды на входе, °С

90

Расход холодной воды на входе, л/мин

20

Расход горячей воды на входе, л/мин

20

Коэффициент теплопроводности материала труб, Вт/(м·°С)

27

Для исследования процессов в полостях теплообменника методом компьютерной гидрогазодинамики применен программный продукт «Autodesk Simulation CFD». При этом с использованием вышеуказанного инструмента определялись следующие параметры: локальные скорости и давления теплоносителей в теплообменнике; коэффициенты теплоотдачи между поверхностями теплообменника и теплоносителями; тепловые потоки; температуры теплоносителей и стенок теплообменника. Отметим, что задача по исследованию рабочих процессов в теплообменнике методами компьютерной гидрогазодинамики решалась для установившихся гидродинамических и тепловых процессов без моделирования конвективного теплообмена между наружной поверхностью теплообменника с окружающей средой. Практическая методика проектирования теплообменных аппаратов также основывается на том, что гидродинамические и тепловые процессы стационарны.

Применение в качестве инструмента расчета теплообменника программного продукта «Autodesk Simulation CFD» определяет набор необходимых операций: создание модели исследуемой области; описания физико-механических свойств материалов; указания граничных условий; разбивки исследуемой области на конечные элементы; опциональной настройки процедуры решения задачи; выполнение процедуры расчета и анализа полученных результатов. Далее рассмотрим некоторые особенности вышеуказанного порядка действий по процедуре исследования. Создание модели исследуемой области подразумевает разработку трехмерной модели теплообменника, которая выполняется с применением соответствующих инструментов моделирования. В итоге перед исследованием согласно рис. 1 создается трехмерная модель теплообменника. Процедуры описания физико-механических свойств материалов являются стандартными функциональными возможностями программного обеспечения и заключаются в выборе из имеющихся баз данных требуемых свойств конкретного материала или формирование их в ручном режиме по имеющимся параметрам применяемого материала. Указанная процедура задания физико-механических свойств выполняется как для материалов, из которых выполнен теплообменник, так и для параметров теплоносителей. Описание граничных условий заключается в указании на входных сечениях трехмерной модели теплообменника таких параметров, как расходы теплоносителей (имеется возможность указать как массовые, так и объемные расходы) и температуры теплоносителей, а на выходных сечениях трехмерной модели теплообменника – перепады давлений. На основании данных по геометрии исследуемой области, свойствам материалов, граничных условий численным методом решается система дифференциальных уравнений [5, 6], состоящая из уравнений движения, неразрывности и энергии (уравнение, учитывающее перенос теплоты путем конвекции и теплопроводностью), дополненное уравнением теплообмена на границе между твердым телом и потоком. В итоге для всей исследуемой области (трехмерной модели теплообменника) определяются скорости и давления теплоносителей, их температуры, а также температуры стенок теплообменника, коэффициенты конвективной теплоотдачи, тепловые потоки.

Практическая методика проектирования теплообменных аппаратов функционально также базируется на решении фундаментальных уравнений гидрогазодинамики и теплообмена, но в отличие от методов компьютерной гидрогазодинамики решение указанных уравнений осуществляется с применением теории подобия с широким использованием опытных данных. В условиях постановки нашей задачи по определению параметров теплообменника (когда неизвестны температуры теплоносителя на выходе из соответствующих полостей) температуры теплоносителей на выходе из его соответствующих полостей определяются по формулам из работы [4]:

garov01.wmf (1)

garov02.wmf (2)

где t1 – начальная температура горячей воды, °С; t1 – конечная температура горячей воды, °С; garov03.wmf – начальная температура холодной воды, °С; garov04.wmf – конечная температура холодной воды, °С; G – массовый расход горячей воды, кг/с; G′ – массовый расход холодной воды, кг/с; cp – удельная теплоемкость горячей воды, Дж/(кг·°С); garov05.wmf – удельная теплоемкость холодной воды, Дж/(кг·°С); k – коэффициент теплопередачи, Вт/(м2·°С); F – площадь теплопередающей поверхности, м2.

Согласно практической методике проектирования теплообменных аппаратов для определения средних коэффициентов теплоотдачи используются критерии Нуссельта, которые применительно к данной задаче при вынужденном продольном турбулентном течении обоих теплоносителей рассчитываются по эмпирической формуле [3]

garov06.wmf (3)

где εl – поправочный коэффициент, учитывающий влияние геометрии теплообменника; Re, Pr – критерии Рейнольдса и Прандтля; Pr, Prст – критерий Прандтля и пристеночный критерий Прандтля соответственно.

Средний по поверхности коэффициент теплоотдачи, Вт/(м2·°С) [3]

garov07.wmf, (4)

где λ – коэффициент теплопроводности теплоносителя, Вт/(м°С); d – диаметр проходного сечения для теплоносителя, м.

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

Согласно данным, приведенным в работе [4], гидравлическое сопротивление для гладких труб при вынужденном продольном турбулентном течении обоих теплоносителей составляет, Па:

garov08.wmf (5)

где ψ(Re) – коэффициент гидравлического сопротивления; L – длина трубы теплообменника, м; ρ – плотность теплоносителя, кг/м3; w – средняя скорость теплоносителя, м/с.

КПД теплообменника для обеих используемых методик определяется одинаково [4]:

garov09.wmf (6)

Тепловая мощность, кВт [4]:

garov10.wmf (7)

Результаты расчета параметров теплоносителей противоточного теплообменника сведены в табл. 2.

Таблица 2

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

Параметр

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

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

Температура холодной воды на выходе, °С

61,01

61,06

Температура горячей воды на выходе, °С

88,97

88,95

Средняя скорость холодной воды, м/с

0,774

0,693

Средняя скорость горячей воды, м/с

0,469

0,473

Средний коэффициент теплоотдачи для холодной воды, Вт/(м2·°С)

2226,55

4763,93

Средний коэффициент теплоотдачи для горячей воды, Вт/(м2·°С)

3890,93

3961,49

Гидравлическое сопротивление по тракту холодной воды, Па

387,55

442,41

Гидравлическое сопротивление по тракту горячей воды, Па

281,76

228,88

КПД

0,034

0,035

Тепловая мощность, кВт

1,407

1,477

Результаты проведенных расчетов с использованием двух методик показывают, что для обоих теплоносителей выходные температуры совпадают. Расхождение значений средних коэффициентов теплоотдачи для соответствующих теплоносителей можно объяснить приближенными расчетами по эмпирической формуле (3), используемой в практической методике проектирования теплообменных аппаратов, в которую невольно заложены погрешности измерений при экспериментальном исследовании, в то время как с помощью метода компьютерной гидрогазодинамики средний коэффициент теплоотдачи определяется непосредственно на основании вычисляемого теплового потока, площади поверхности теплопередающего элемента и разности средних температур между теплоносителем у стенки и самой стенки. В подтверждение этого был выполнен проверочный расчет выходных температур теплоносителей по формулам (1) и (2) при известных значениях средних коэффициентов теплоотдачи, полученных с применением технологии компьютерной гидрогазодинамики. По результатам данного проверочного расчета выходные температуры изменились на 0,3 °С, что укладывается в пределы погрешности термометров сопротивления. Сравнение гидравлических сопротивлений также выявило расхождение в полученных результатах из-за возможных неточностей в эмпирической формуле (5). Кроме того, стоит отметить, что технология компьютерной гидрогазодинамики позволяет визуально оценить полученные результаты расчета (рис. 2, 3, 4, 5).

pic_10.tif

Рис. 2. Гидравлическое сопротивление по тракту холодной воды

pic_11.tif

Рис. 3. Гидравлическое сопротивление по тракту горячей воды

pic_12.tif

Рис. 4. Картина распределения скоростей потоков теплоносителей

pic_13.tif

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

Отметим, что у используемой практической методики проектирования теплообменных аппаратов, описанной в работах [1, 2, 3, 4], есть один существенный недостаток. Он состоит в том, что расчеты основаны на использовании балансовых уравнений по всему теплообменнику, где потоки теплоты определяются по некоторым средним значениям параметров состояния взаимодействующих теплоносителей. Это приводит к существенным погрешностям расчета, так как разность потенциалов теплообменных процессов существенно меняется по длине процесса, и предсказать заранее характер этих изменений довольно сложно, а иногда и практически невозможно. Но здесь необходимо отметить также и достоинство практической методики проектирования теплообменных аппаратов, которое состоит в том, что она позволяет без использования мощной вычислительной техники приближенно определять средние коэффициенты теплоотдачи в теплообменниках различной конструкции, а на основании их значений вычислить передаваемые тепловые потоки или, напротив, для заданного передаваемого теплового потока определить необходимую площадь теплообмена.

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

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

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

Рецензенты:

Мурашов А.А., д.т.н., заведующий кафедрой математических и естественнонаучных дисциплин Московского финансово-юридического университета МФЮА, г. Ярославль;

Епархин О.М., д.т.н., профессор, директор Ярославского филиала ФГБОУ ВПО «Московский государственный университет путей сообщения», г. Ярославль.

Работа поступила в редакцию 16.12.2013.