УДК 621.321
РАЗРАБОТКА ЭМПИРИЧЕСКОГО ОБЕСПЕЧЕНИЯ МАТРИЧНОЙ МОДЕЛИ НАГРЕВА И ДЕАЭРАЦИИ ВОДЫ В СТРУЙНЫХ ОТСЕКАХ АТМОСФЕРНЫХ
ДЕАЭРАТОРОВ
А.В. МОШКАРИН, Е.В. БАРОЧКИН, В.П. ЖУКОВ, Г.В. ЛЕДУХОВСКИЙ, А.А. КОРОТКОВ ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
Иваново, Российская Федерация E-mail: admin@tes.ispu.ru
Авторское резюме
Состояние вопроса: Преимуществом существующих матричных математических моделей процессов тепломассообмена в многоступенчатых системах со сложной конфигурацией потоков является их универсальность - применение их к конкретным условиям требует лишь наличия значений параметров идентификации. Использованию этих моделей для расчета характеристик процессов термической деаэрации воды в струйных отсеках деаэраторов атмосферного давления препятствует отсутствие эмпирического обеспечения, разработка которого и является целью настоящей работы.
Материалы и методы: Идентификация модели выполнена по результатам экспериментальных исследований процессов деаэрации в двух аппаратах, включающих в общей сложности четыре различающихся по конструкции и условиям работы струйных отсека. Эмпирическое обеспечение включает метод определения эффективной площади поверхности контакта фаз в струйном потоке и записанные в безразмерном виде уравнения для расчета коэффициентов теплопередачи и массопередачи.
Результаты: Получены значения параметров идентификации матричной математической модели термической деаэрации воды в струйных отсеках деаэраторов атмосферного давления, описывающей процессы нагрева воды в струйном потоке при конденсации водяного пара из движущейся парогазовой смеси и десорбции растворенного кислорода. Показано, что матричная модель совместно с полученными параметрами ее идентификации позволяет рассчитать температуру воды и концентрацию растворенного в ней кислорода за струйными отсеками деаэраторов атмосферного давления с точностью 10,8 % при расчете теплообмена и 17,2 % при расчете десорбции кислорода. Выводы: Полученная модель может быть использована для расчета процессов тепломассообмена при термической деаэрации воды в задачах конструирования струйно-барботажных деаэраторов атмосферного давления и режимной наладки эксплуатируемых аппаратов этого типа.
Ключевые слова: термический деаэратор, струйный отсек, теплообмен при конденсации, десорбция кислорода.
DEVELOPMENT OF EMPIRICAL SUPPORT FOR MATRIX MODEL OF WATER HEATING AND DEAERATION IN SPRAY SECTIONS OF ATMOSPHERIC DEAERATORS
A.V. MOSHKARIN, E.V. BAROCHKIN, V.P. ZHUKOV, G.V. LEDUHOVSKY, A.A. KOROTKOV Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: admin@tes.ispu.ru
Abstract
Background: The authors describe universality of the existing matrix mathematical models of heat-mass exchange processes in multistage systems with complex stream configuration as their advantage. Their usage in specific conditions requires only identification parameters. The lack of empirical support prevents from applying these models for calculation of process characteristics of thermal water deaeration in spray sections of deaerators of atmospheric pressure. The authors have the aim to develop such an empirical support in the article.
Materials and methods: Model identification is made according to the results of experimental researches of deaeration processes in two devices, which include four distinguished in construction and operation conditions spray sections. Empirical support includes the determination method of effective area of phase contact surface in spray stream as well as written equalizations in non-dimential aspect for calculating the coefficient of heat and mass transfers. Results: The authors provide the values of identification parameters of the matrix mathematical models of thermal water deaeration in spray sections of deaerators of atmospheric pressure. This mathematical model describes the processes of water heating in spray stream with the water vapour condensation from the moving gas-vapor mixture and desorption of dissolved oxygen. The authors prove that the matrix model combined with the received parameters of its identification allows to calculate the water temperature and concentration of dissolved oxygen outside the spray sections of deaerators of atmospheric pressure with accuracy of 10,8 % of heat exchange calculation, and 17,2% of oxygen desorption calculation. Conclusions: The developed model can be used for calculating heat-mass exchange processes with thermal water deaeration in designing problems of spray deaerators of atmospheric pressure and security adjustment of operated devices of this type.
Key words: thermal deaerators, spray section, heat transfer in condensation process, oxygen desorption.
Эффективность работы тепломассооб-менного аппарата зависит от его конструктивного исполнения и режима работы. Наладка режима работы деаэратора заданной конструкции на практике производится в ходе натурных испытаний, требующих определенных материальных и трудовых ресурсов. Затраты этих ресурсов могут быть сокращены при использовании компьютерных программ, обеспечивающих возможность расчета показателей работы деаэраторов с требуемой точностью, в том числе, при изменении конструкции.
Создание математической модели деаэратора в целом предполагает наличие моделей отдельных струйных и барботажных элементов. Количество таких отдельных элементов в деаэраторах современных конструкций может достигать десяти, а схема их соединения по потокам воды и греющего пара предусматривать рециркуляцию, объединение и разделение потоков. Однако даже при наличии адекватных моделей каждого элемента в отдельности при изменении конструкции аппарата каждый раз приходится «вручную» корректировать программные коды, описывая сложные взаимосвязи между элементами по потокам теплоносителей.
Указанная проблема разработки прикладных расчетных программ может быть решена при переходе к универсальным математическим моделям элементов, разработанным в рамках матричной формализации расчета тепломасообменных систем со сложной конфигурацией потоков [1].
Математическая модель теплообмена и десорбции растворенного кислорода в деаэраторе представлена системой дифференциальных уравнений, которые описывают изменение температуры воды, расхода пара и расхода воды, концентрации кислорода в водяной и паровой фазах вдоль поверхности контакта фаз Р [1]:
С т СР
ССС-СР Св2
кт к т2
1 _
О^СС2 кт,
Г1, ссС1
ГА
СР СР '
Сод1 _ кт (кдОд2 - °д1)
(1)
дР
дО
д 2
6,
кт(кдОд2 - Од1)
дР
где т - температурный напор; к -коэффициент теплопередачи; с* - удельная теплоемкость; г - удельная теплота парообразования; сд -концентрация газа; кт - коэффициент массо-передачи по рассматриваемому газовому ком-
поненту; кд -коэффициент фазового равновесия, определяющий связь между концентрацией газа в воде и равновесной концентрацией газа в паровой фазе; нижний индекс 1 относится к горячему теплоносителю; 2 - к холодному; п - к состоянию насыщения.
Баланс аддитивных характеристик потоков позволяет представить расчетное уравнение процесса для систем со сложной конфигурацией потоков в матричном виде [1]:
(
-I К21В1
К12В2 -I
V Кп1В1 кп 2В2 •• ([VI.^ Г-[Хех]1 ^
-[ Хех ]2
КщВп ^ К2пВп
-I
[ V ]1 [ V ]2
(2)
[Х ]п
-[ Хех ]п
где В - матрица процесса в ступени; К - матрица коммутации; X - вектор аддитивных параметров; I - единичная матрица; индекс «ех» указывает на внешний поток, подаваемый на вход ступени.
При известных матрицах В и К решение системы (2), которая для целого ряда практически важных случаев является системой линейных уравнений, позволяет определить значения параметров теплоносителей в любой точке установки. Авторы подхода дают простые рекомендации относительно порядка решения задачи в рамках матричной формализации [1]. При этом не только вычисления, но и предшествующее составление всех матриц хорошо поддаются автоматизации (формализации), что и требуется для решения поставленной выше задачи автоматизированного синтеза математической модели деаэратора на основе моделей его отдельных элементов.
Рассмотренная математическая модель требует идентификации - для элементов остаются неизвестными площадь поверхности контакта фаз Р, коэффициент теплопередачи к и коэффициент массопередачи кт. Проблема связана с тем, что опубликованные модели струйных и барботажных элементов деаэраторов [2] традиционно представлены в статистическом виде и обычно устанавливают зависимость непосредственно изменения температуры воды или концентрации растворенного в ней кислорода в элементе от каких-либо параметров. То есть такие модели не позволяют рассчитать собственно значения коэффициентов тепло- и массопередачи.
Анализ опубликованных данных [2] показывает, что наиболее сложными с точки зрения математического моделирования являются струйные отсеки деаэраторов. Ранее нами проведены испытания деаэраторов ДСА-300 и
2
ДА-300м с отбором проб воды из внутренних элементов [3]. По результатам работы предложены критериальные уравнения, обеспечивающие расчет процессов нагрева воды и десорбции из нее кислорода для условий работы струйных отсеков при струйно-капельном гидродинамическом режиме [4]:
lg
^ | = 0,0137(H
tu -12) I d
<Fr
н '2 0,054
0,982
Lap
0,163 ,
pr-0,832 k-0,290
lg(С | = 0,8910 f H
0,964
Lap =
К =
Рж^п 2d
, Fr =
gd
Lap0264Sc-1K
-, Рг — —,
0,882
(3)
(4)
Sc — , Wж —
ж
4Gl
С ж (t2 - " 0ж 'ж nd 2првжх'
при Wnex/ Wneba > 1,7,
W вх - Wвых вх вь,х
W — п п 1Л/ вх /|Л/вых-
п Wвх
2,3lg Wn—
Wпвых
Wn — ^
Wn вх + Wn вых
2
при Wn* / Wnba < 1,7,
I-) вх вх I-) вых вых
W вх = ип ип W вых = ип ип
п Q ' п Q '
где t1, t2 - температура воды на входе и выходе струйного отсека соответственно, оС; tH -температура насыщения пара в отсеке, оС; с1 и с2 - концентрации растворенного в воде кислорода на входе и выходе струйного отсека соответственно (проба при 25 оС), мкг/дм3; Lap, Fr, Pr, K и Sc - критерии Лапласа, Фруда, Прандтля, Кутателадзе и Шмидта соответственно; Н - высота отсека от нижней поверхности струеобразующей тарелки до борта нижней тарелки, м; d - диаметр отверстий струе-образующей тарелки, м; рж - средняя плотность воды, кг/м3; Wп - средняя скорость парового потока в отсеке, м/с; стж - средний коэффициент поверхностного натяжения воды, Н/м; W^c - средняя скорость истечения воды из отверстий струеобразующей тарелки, м/с; g - ускорение свободного падения, м/с2; vж - коэффициент кинематической вязкости воды, м2/с; аж - коэффициент температуропроводности воды, м2/с; г - скрытая теплота парообразования, кДж/кг; сж - истинная изобарная теплоемкость воды, кДж/(кг оС); иж - коэффициент молекулярной диффузии кислорода в воде, м2/с; Gжвх - массовый расход воды на входе в отсек, кг/с; п - число отверстий струеобразующей тарелки, шт.; ржвх - плотность воды на входе в струйный отсек, кг/м3; Wпвх, Wпвых - скорости парового потока во входном и выходном сечениях отсека соответственно, м/с; ипвх и ипвых -
удельные объемы пара во входном и выходном сечениях отсека соответственно, м /кг; &вх и Овых - площади проходного сечения на входе и выходе пара из отсека соответственно, рассчитанные через высоту отсека Н, за вычетом площади, занимаемой условными струями воды постоянного по высоте диаметра С, м2.
Критериальные уравнения (3) и (4) характеризуются высокими показателями точности по сравнению с ранее опубликованными моделями струйных отсеков, что позволяет рекомендовать их для расчетов [3, 4]. Однако полученные в таком виде уравнения не могут быть использованы в модели (2), поскольку не позволяют рассчитать в явном виде значения коэффициентов теплопередачи к и массопе-редачи кт.
В целях получения расчетных зависимостей для коэффициентов к и кт проведена специальная обработка результатов испытаний деаэраторов ДСА-300 и ДА-300м. Порядок определения расходов, теплофизических и химических характеристик теплоносителей на входе и выходе каждого струйного отсека подробно описан ранее [3, 4] и поэтому здесь не приводится.
Значения искомых коэффициентов к и кт определяются по минимальному рассогласованию между экспериментальными и рассчитанными согласно (2) значениями выходных параметров. При определении коэффициента теплопередачи таким выходным параметром являлась температура воды на выходе из отсека, а при определении коэффициента мас-сопередачи - концентрация компонента в воде на выходе из струйного отсека.
При идентификации модели (2) с определением коэффициентов к и кт площадь поверхности контакта фаз в струйном отсеке Р должна быть известна. Анализ опубликованных данных по гидродинамике в струйных отсеках теплообменных аппаратов [2, 5] позволил выбрать базовое расчетное выражение для определения искомой площади Рстр, м2, справедливое для чисто струйного режима течения воды в отсеке:
2ndW2 Fcmp 3ц3/2д
W2
\3/4
- 1
L —
вх^ 3
(5)
(6)
где ц - коэффициент расхода, принимаемый равным 0,75 при диаметре отверстий тарелки 5-8 мм и толщине этой тарелки 4-6 мм [2, 5]; £ - средняя высота зоны чисто струйного режима течения воды в отсеке, м.
Все режимы работы струйных отсеков рассматриваемых деаэраторов являются
ж
ж
ж
струйно-капельными [3, 4]. В таких условиях выражение (5), примененное в целом к высоте отсека Н, дает заниженные значения Р. При определении площади поверхности контакта фаз предложен следующий подход. Площадь поверхности контакта определяется как сумма площадей поверхности струй и поверхности капель. Длина и площадь поверхности струй определяются согласно (5), (6). Для определения площади поверхности капель в выражение (5) вводится поправочный коэффициент, диапазон значений которого по результатам специального анализа фотографических изображений составил от 1,4 до 1,6, а вместо длины струй £ используется высота зоны капельного режима (Н - £).
Полученные в результате обработки данные позволили записать критериальные уравнения, обеспечивающие расчет искомых параметров идентификации к и кт:
N11 = ехр(15,021) | Н
<рг0,288 рг-1,833 к-0,865
1,594
Ьар'
0,046 ,
БИ = ехр(5,044)
0,897
Ьар'
0,318 ,
(7)
(8)
хрг°,468Бс -0,413к-а802,
Nu ■
кб
К ж
Б И ■
ктб °ж Рж
где Nu, БИ - критерии Нуссельта и Шервуда соответственно; к - средний по поверхности контакта фаз коэффициент теплопередачи в отсеке, Вт/(м2К); "ж - средний в отсеке коэффициент теплопроводности воды, Вт/(м-К); кт -средний по поверхности контакта фаз коэффициент массопередачи по растворенному кислороду в отсеке, кг/(м2с).
Выраженные в процентах среднеквадра-тические отклонения, характеризующие точность предложенных уравнений (7) и (8), составили: для модели теплообмена ± 10,8 % (рис. 1); для модели десорбции кислорода ± 17,2 % (рис. 2).
Для полученных ранее уравнений (3) и (4) эти показатели составили ±1,8 и ±8,2 % соответственно. Сравнение показателей точности моделей (3), (4) и (7), (8) позволяет сделать вывод, что внесение в модель дополнительной неопределенности в виде площади поверхности контакта фаз Р привело к заметному ухудшению точности. Однако полученные показатели точности уравнений (7), (8) в целом характерны для критериальных уравнений, описывающих процессы конденсации и испарения [5].
ззо
280
230
180
130
80
30
Ми» А* / ■
А */
А А //Г ■ /' о / • ■'0 о
й * . * 1,' / 6 л
N4'
30
80
130
180
230
280
330
Рис. 1. Сопоставление экспериментальных и расчетных значений критерия Нуссельта: Nuэ и Nuр -экспериментальные и расчетные значения критерия Нуссельта соответственно; сплошная линия - совпадение экспериментальных и расчетных значений; пунктирные линии - границы 95 %-го доверительного интервала, обусловленного погрешностью измерения параметров в условиях каждого опыта; точки - результаты расчета по разработанному критериальному уравнению: о и • - первая и вторая по ходу воды зоны верхнего струйного отсека деаэратора ДА-300м соответственно; ▲ и Л - верхний и нижний струйные отсеки деаэратора ДСА-300 соответственно
14000
12000
10000
8000
6000
4000
2000
/ ■ •
• , / / •
/
Л .•'/.О
Ж' 1 т Л
8 Ь3
0 2000 4000 6000 8000 10000 12000 14000
Рис. 2. Сопоставление экспериментальных и расчетных значений критерия Шурвуда: Биэ и БИр - экспериментальные и расчетные значения критерия Шервуда соответственно; прочие обозначения аналогичны обозначениям на рис. 1
Полученные уравнения (7), (8), как и уравнения (3), (4), справедливы при следующих параметрах, обеспеченных испытаниями: диаметр отверстий струеобразующей тарелки от 0,006 до 0,01 м; высота струйного отсека от
0,3 до 0,95 м; среднее абсолютное давление пара в отсеке от 109 до 137 кПа; средняя скорость истечения воды из отверстий струеобра-зующей тарелки от 0,2 до 3,0 м/с; средняя скорость пара в отсеке от 0,8 до 48,2 м/с.
Нужно отметить, что определяемые по уравнениям (7), (8) коэффициенты теплопередачи k и массопередачи km не следует рассматривать как истинные, поскольку их значения в условиях опытов определялись не только характеристиками собственно процесса, но и значениями площади поверхности контакта фаз Р, которые являются условными.
Несмотря на ухудшение прогностических свойств уравнений (7), (8) по сравнению с уравнениями (3), (4), обеспечена возможность их использования в модели (2), что, в свою очередь, открывает возможность автоматического учета изменений конструкции деаэратора в компьютерной программе, то есть создания систем автоматизированного проектирования деаэраторов.
Список литературы
1. Жуков В.П., Барочкин Е.В. Системный анализ энергетических тепломассообменных установок / ГОУВПО «Ивановский государственный энергетический университет им. В.И. Ленина». - Иваново, 2009. - 176 с.
2. Шарапов В.И., Цюра Д.В. Термические деаэраторы. - Ульяновск: Ульян. гос. техн. ун-т., 2003. - 560 с.
3. Экспериментальные исследования и моделирование технологических процессов атмосферной струй-но-барботажной деаэрации воды / А.В. Мошкарин, В.Н. Виноградов, Г.В. Ледуховский и др. // Теплоэнергетика. -2010. - № 18. - С. 21-25.
4. Ледуховский Г.В. Совершенствование технологии десорбции кислорода в струйно-барботажных деаэраторах атмосферного давления: дис. ... канд. техн. наук: 05.14.14. - Иваново, 2008. - 226 с.
5. Кутателадзе С.С. Теплопередача при конденсации и кипении. - 2-е изд., перераб. и доп. - М.: Гос. науч.-техн. изд-во машиностроит. лит., 1952. - 231 с.
References
1. Zhukov, V.P., Barochkin, E.V. Sistemnyy analiz en-ergeticheskikh teplomassoobmennykh ustanovok [Systems Analysis of power engineering heat and mass exchange units], Ivanovo, 2009, 176 p.
2. Sharapov, V.l., Tsyura, D.V. Termicheskie deaeratory [Thermal Deaerators], Ul'yanovsk: UGTU, 2003, 560 p.
3. Moshkarin, A.V., Vinogradov, V.N., Leduhovsky, G.V., Korotkov, A.A., Barochkin, A.E. Teploenergetika, 2010, issue 8, pp. 21-25.
4. Ledukhovskiy, G.V. Sovershenstvovanie tekhnologii desorbtsii kisloroda v struyno-barbotazhnykh deaeratorakh atmosfernogo davleniya [Technological Development of Oxygen Desorption in Spray-buble Atmospheric Deaerators], dis. ... kand. tekhn. nauk, 05.14.14., Ivanovo, 2008, 226 p.
5. Kutateladze, S.S. Teploperedacha pri kondensatsii i kipenii [Heat Transfer in Condensation and Boiling], Moscow: Gosudarstvennoe nauchno-tekhnicheskoe izdatel'stvo mashi-nostroitel'noy literatury, 1952, 231 p.
\Мошкарин Андрей Васильевич,
ФГБОУВПО «Ивановский государственный энергетический университет им. В. И. Ленина»,
Заслуженный деятель науки Российской Федерации, Почетный работник высшего профессионального образования России, доктор технических наук, профессор, заведующий кафедрой тепловых электрических станций
Барочкин Евгений Витальевич,
ФГБОУВПО «Ивановский государственный энергетический университет им. В.И. Ленина», доктор технических наук, профессор кафедры тепловых электрических станций, адрес: г. Иваново, ул. Рабфаковская, д. 34, кор. В, ауд. 408, телефон (4932) 41-60-56, 26-99-31, e-mail: admin@tes.ispu.ru
Жуков Владимир Павлович,
ФГБОУВПО «Ивановский государственный энергетический университет им. В.И. Ленина»,
доктор технических наук, профессор кафедры прикладной математики,
адрес: г. Иваново, ул. Рабфаковская, д. 34, кор. А, ауд. 202,
телефон (4932) 26-97-45,
e-mail: zhukov@home.ivanovo.ru
Ледуховский Гоигорий Васильевич,
ФГБОУВПО «Ивановский государственный энергетический университет им. В. И. Ленина»,
кандидат технических наук, доцент кафедры тепловых электрических станций,
адрес: г. Иваново, ул. Рабфаковская, д. 34, кор. В, ауд. 408,
телефон (4932) 41-60-56, 26-99-31,
e-mail: lgv83@yandex.ru
Коротков Александр Александрович,
ФГБОУВПО «Ивановский государственный энергетический университет им. В.И. Ленина»,
старший преподаватель кафедры тепловых электрических станций,
адрес: г. Иваново, ул. Рабфаковская, д. 34, кор. В, ауд. 408,
телефон (4932) 41-60-56, 26-99-31,
e-mail: admin@tes.ispu.ru