CC BY
27
УДК 69.027.1
Володин Ю.Г. - кандидат технических наук, доцент
E-mail: [email protected]
Ханнанов Р.Р. - аспирант
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Математическое моделирование рабочего процесса в газовоздушном тракте котлового агрегата ТГМ-84
Аннотация
В статье приведена математическая модель турбулентного течения дымовых газов в газовоздушном тракте котла ТГМ-84. Расчеты выполнены с использованием программного комплекса ANSYS. При математическом моделировании рабочего процесса в газовоздушном тракте реализована идея дросселирования потока дымовых газов в районе отводящего газового патрубка за РВП. Получены оптимальные размеры дроссельной заслонки. Выполненные измерения присосов воздуха на работающем котле до и после установки заслонки показали заметное снижение присосов воздуха в газовый тракт.
Ключевые слова: присосы воздуха, математическая модель, дроссельная заслонка, котел ТГМ-84, воздухоподогреватель РВП-54, газовоздушный тракт.
Научно-технический прогресс стимулирует развитие и совершенствование различных отраслей, в том числе и энергоемких, в промышленном комплексе. Освоение новых месторождений полезных ископаемых, расширение масштабов добычи на действующих площадках требует формирования соответствующей инфраструктуры для жизнедеятельности персонала, что в итоге приводит к строительству новых или расширению площадей уже существующих населенных пунктов. Все вместе взятое нуждается в немалых количествах электрической и тепловой энергии. Этим объясняется не только введение новых мощностей в энергетическом комплексе страны, но модернизация действующих агрегатов и совершенствование технологических процессов. Одним из актуальных направлений в теплоэнергопроизводящей промышленности является энергоресурсосбережение при эксплуатации технологического оборудования.
к турбине Hit
S 9 Ю
t t тагпипа
Рис. 1. Технологическая схема ТЭС (фрагмент). Основные элементы парового котла:
I - газовые горелки; 2 - экраны труб, для подогрева питательной воды;
3, 4, 5 - потолочный, ширмовый, конвективный пароперегреватель; 6, 7 - экономайзер;
8 - линия рециркуляции; 9 - дутьевой вентилятор; 10 - дымосос;
II - регенеративный воздухоподогреватель (РВП); 12 - дымовая труба
Тепловую электрическую станцию (ТЭС) образует комплекс оборудования и устройств, основными из которых являются котельная установка, производящая пар высоких параметров, турбинная или паротурбинная установка, преобразующая теплоту пара в механическую энергию ротора турбоагрегата, и электрические устройства
(электрогенератор, трансформатор и т.д.), обеспечивающие выработку электроэнергии (рис. 1). Основным элементом в котельной установке, производящей пар, является котел. Котел выглядит в виде П-образной конструкции с газоходами прямоугольного сечения. В левую часть, называемую топкой, подается смешанный с подогретым воздухом газ. Воздух подогревается путем рециркуляции дымовых газов, вышедших из котла. Подогревается воздух в регенеративном воздухоподогревателе (РВП) дымовыми газами. Частично охлажденные дымовые газы дымососом ДС подаются в дымовую трубу. Через дымовую трубу продукты сгорания (дымовые газы) выталкиваются в верхние слои атмосферы и рассеиваются. Кроме того, дымовая труба необходима для создания разрежения в топке и газоходах котла.
Одним из паровых котлов, которые в настоящее время эксплуатируются на ТЭС, является ТГМ-84 (ТГМ-84А, ТГМ-84Б), в газовоздушном тракте которого устанавливается регенеративный воздухоподогреватель РВП-54 (или РВП-68). Во время функционирования происходят присосы (перетоки) воздуха (рис. 2) из воздушного тракта в газовый, что снижает полезную работу вращающегося РВП. Эксплуатационный опыт показывает, что более 20 % воздуха попадает в газовый тракт, не совершая работы. Соответственно ухудшается показатель по присосам воздуха. Данное обстоятельство приводит к увеличению нагрузки на дымососы и дутьевые вентиляторы, заметному ухудшению по величине такого показателя как удельный расход электроэнергии на тягу и дутье и в итоге к перерасходу топлива.
I иты
Рис. 2. Иллюстрация движения дымовых газов, воздуха и присосов воздуха в РВП
Для снижения присосов воздуха в газовый тракт котла имеется большое количество разнообразных конструкторских решений: установка промежуточных радиальных полос, войлочных пластин на прижимные планки, максимальное уплотнение щелей и зазоров в конструкции РВП и т.п. Реализация таких решений является эффективной и приводит к снижению присосов воздуха в РВП до 10 %, вместо 15-20 %. Однако, все они обладают одним недостатком - механическим износом.
Вопрос снижения присосов воздуха в газовый тракт котла можно решить путем дросселирования дымовых газов (рис. 3) в хвостовой части газохода РВП. Для этого необходимо установить дросселирующее устройство. Известно [1, 2], что при прохождении потока через дросселирующее устройство, происходит местное сжатие потока. Согласно закону сохранения энергии суммарная энергия движущейся среды может уменьшиться только на величину потерь, следовательно, часть потенциальной энергии потока преобразуется в кинетическую энергию. Давление потока на стенки трубопровода несколько возрастает за счет подпора перед сужающим устройством и далее понижается до некоторого минимума за дросселем. Таким образом, разогнав поток и снизив давление на стенки газовоздушного тракта, можно достичь снижения присосов воздуха.
и ■¿.■гил -л.иьл ь.-зиь-
ЗЛТ ^.ат» 1-ТА'л ъ.ии
Рис. 3. Дросселирование газового потока
Течение неизотермического газового потока, в том числе, и при его ускорении исследовалось в работах [3-6]. Общеизвестно, что при сужении канала формируется отрицательный продольный градиент давления, под воздействием которого профиль скоростей становится более заполненным, относительный коэффициент трения по величине увеличивается, увеличивается и скорость течения газового потока. Неизотермичность при тепловом потоке, направленном от более нагретых дымовых газов к стенкам газовоздушного тракта, также способствует большей заполненности профиля скоростей, увеличению относительного коэффициента трения и скорости течения дымовых газов.
При турбулентном режиме течения рабочего тела поля актуальных значений скоростей, давлений и температур в потоке представляют собой сложную структуру. Причиной такой структуры является нерегулярность или случайность характера изменений названных параметров в пространстве и во времени. Поэтому при математическом описании турбулентных течений все подходы, включая прямое моделирование турбулентности, прямое моделирование крупномасштабных структур и так называемый моментный подход Рейнольдса, основаны на нестационарных трехмерных уравнениях динамики вязкого газа - система уравнений Навье-Стокса, а также на выполнении операций статистического осреднения, которое реализуется на одном из этапов при формировании постановки задачи либо при ее решении. Когда выполняют прямое моделирование турбулентности на основе трехмерных нестационарных уравнений Навье-Стокса, то статистическое осреднение делают после интегрирования этих уравнений. Осреднение преследует основную цель, заключающуюся в удобстве традиционной формы практического анализа.
Когда используется моментный подход Рейнольдса, то на начальном этапе постановки задачи выполняется операция статистического осреднения. При таком подходе сначала осредняются статистически уравнения динамики вязкого газа, затем проводится интегрирование после предварительного полуэмпирического замыкания. Операция осреднения выполняется по всем масштабам одновременно. Поэтому здесь утрачивается универсальность и эти методы остаются востребованными при анализе узкого диапазона гидродинамических процессов.
Когда проводится прямое моделирование, то статистическое осреднение уравнений выполняется на этапе постановки задачи. При этом охватываются не все масштабы, так как полуэмпирическое моделирование турбулентности необходимо проводить только на подсеточных масштабах, т.е. на масштабах явлений, которые меньше чем размеры статистической ячейки. Таким образом, во всех случаях получается система уравнений Навье-Стокса, осредненная по Рейнольдсу. Это означает, что любая газодинамическая величина в турбулентном потоке представляется в виде суммы осредненных и пульсационных составляющих. Причем, при определении средних значений осреднение выполняется по некоторому интервалу времени, т.е. временное осреднение. Кроме такого осреднения существуют и другие способы - статистического осреднения по объему (пространственное осреднение), статистическое осреднение по ансамблю.
Классический подход статистического осреднения по некоторому интервалу времени (по Рейнольдсу) в основном применяют при решении уравнений, которые описывают течение рабочего тела с постоянной плотностью - несжимаемые жидкости. Когда речь идет об описании турбулентных течений сжимаемого газа, то чаще пользуются комбинированным способом осреднения. Здесь давление и плотность газа
осредняют по времени (по Рейнольдсу), а другие параметры потока вводят по так называемым средневзвешенным значениям. Этот способом осреднения называют методом средневзвешенных параметров или осреднением по Фавру. Его главным преимуществом является то, что система уравнений турбулентного движения оказывается более компактной по сравнению другими подходами.
При описании турбулентного течения газового потока, в том числе, и в дымовых трубах ТЭС, система уравнений включает в себя уравнения количества движения, уравнение неразрывности, уравнение энергии и уравнения к-е модели турбулентности, т.е. при трехмерной постановке в декартовой системе координат имеем: - уравнение неразрывности:
дpu + дpv + дpw _ о
дx су дz
ду
Уравнение количества движения в проекции на ось X:
д ( ди Л д ( ду
дри + д (рии) + д (риу) + д(puw) дг дх ду дх
дх I т дх
ди Л д ( ди Л др
ц— |—| т— \+ _о.
ду 0 дх Г дх 0 дх
Уравнение количества движения в проекции на ось У:
дру + д(риу) + д(руу) + д(рvw) д ( ду\ д дг дх ду дх дх I дх 0 ду
(
ду Л д ( ду Л др
т— |—| ц— \+—_ о.
ду 0 дх Г дх 0 ду
Уравнение количества движения в проекции на ось Z:
дpw + д (puw) + д( pуw) + д( pww)
дг
дх
ду
дх
д ( дw —| Ц—
дх I дх
_д_ ду
(
Ц
дн> ду
д ( ды Л др
--1 Ц— \+ _ о,
дх I дх 0 дх
(1)
(2)
(3)
(4)
где х, у, х - декартовы координаты, г - время, и, у, ы - проекции вектора скорости на оси х, у, х соответственно, р - давление, ц - коэффициент эффективной вязкости, р -плотность рабочего тела. Уравнение энергии:
д(рк) д(рки) д(рку) д{ркм>) _ + Ц | дк + ^ + ^^дк! + + (5)
дг + дх + ду + дх ~ дх [ Рг Ргт05х ду [ Рг Ргт05у дх [ Рг Ргт05х' ( )
где Рг - число Прандтля, Ргт - турбулентное число Прандтля, к - энтальпия:
к _ X Уа^^О Сра (Т)йТ)
(6)
где Уа - массовая доля компонента а, а Т - температура.
Уравнения переноса кинетической энергии турбулентности и скорости ее диссипации:
д(рк) д(рки) д(рку) д(ркм>) д
дг
дх
ду
дх
дх
цт\дк д (
Цмол + ^ цмол + ^ Цмол + ^ + -ре,
дх ду Г
дк д
ду дх
Цт
дк
дх
(7)
+ д[риг) + д(руг) + д(р^) _ .ЭГ + цЛ®. + + Ь^ + + цЛ® + р(рло - се2 в) £. (8)
дг дх ду дх дх Г пв)дх ду Г пв)ду дх Г пв)дх к у '
Для замыкания системы осредненных уравнений используется к-е модель, генерация турбулентности за счет сдвиговых напряжений:
О _Цт-р
(
„ди ду (ди Л „ди ды (ди Л2 „ду ды ( ду Л2 ( ду Л2 ( ды Л2 ( ды
2--+ | —| + 2--+ | —| + 2--+ | —| +| —| +| — | +| —
ду дх ГСу 0 дх дх \дх 0 дх ду Гдх 0 \дх 0 Гдх 0 Г ду
2
2 (ди ду ды -| — + — + —
3 Г дх ду дх
, цт I ди ду ды к + — | — + — + —
р Г дх ду дх
Л
Турбулентная вязкость определялась по формуле Колмогорова-Прандтля:
Г к2
Цт _ Гцр-,
(9)
(10) (11)
а эффективная вязкость по формуле:
Ц_Цт +Цмол,
где цмал - молекулярная вязкость. Значения констант принимались равными Г=0,09, Ге1=1,30, ГВ2=1,92, Ге1=1,30, ок=1,0, ае=1,3, Ргт=0,7, 8от=0,7.
а
П
П
П
к
к
к
Плотность р определялась из уравнения состояния:
рД0Т
Р = >
М
(12)
где И - универсальная газовая постоянная, а давление Р полагали равным атмосферному.
Рис. 4. Вращающийся регенеративный воздухоподогреватель
Таблица
Давление на воздушной стороне перед РВП Р(В), кПа Давление на газовой стороне за РВП Р(Г), кПа Разность давлений АР, кПа Объем присосов V, м3 Процент снижения присосов, %
До После До После До После
1 7 -0,29 0,74 7,29 6,26 0,0729772 0,06763 7,333362939
2 6,86 -0,37 0,54 7,23 6,32 0,0726762 0,06795 6,504783741
3 6,84 -0,4 0,373 7,24 6,467 0,0727265 0,06873 5,489045915
4 6,94 -0,3 0,62 7,24 6,32 0,727265 0,06795 6,569374572
5 7,05 -0,34 0,4 7,39 6,65 0,073476 0,0697 5,138802348
Рис. 5. Изменение объема присосов
Моделирование течения дымовых газов в хвостовой части газохода РВП выполнено с использованием программного комплекса ЛЫБУБ. Все размеры газохода соответствуют технической документации на конструкцию РВП, а именно: высота шахты 5,0 м, глубина 5,076 м, ширина 2,493 м. Сужающее устройство представляет собой
дроссельную заслонку [7, 8]. В результате математического моделирования течения дымовых газов в хвостовой части газохода определены оптимальные параметры - высота выступа (рис. 4) дроссельной заслонки (0,1-0,17)-а?экв (эквивалентный диаметр отводящего газового патрубка) и место ее расположения в газовом тракте. Увеличение давления на 1-2 мм вод. ст. приводит к выравниванию давлений за РВП потока дымовых газов с давлением воздуха в подводящем воздушном патрубке.
Дроссельная заслонка была установлена на РВП-54 котлоагрегата ТГМ-84Б. В процессе эксплуатации были проведены измерения присосов воздуха. Сравнительные результаты приведены на рис. 5.
Вывод. Путем математического моделирования спрогнозирован рабочий процесс в газовоздушном тракте за РВП-54. Дросселирование дымовых газов за счет перераспределения потенциальной и кинетической энергии потока позволяет получить эффект снижения присосов воздуха. Реализация данного подхода в производственных условиях позволила реально снизить объем присосов воздуха в газовый тракт.
Список библиотечных ссылок
1. Володин Ю.Г., Богданов А.Н. Технические измерения и приборы. Расчет сужающих устройств. Учебное пособие. - Казань: Казан. гос. энерг. ун-т, 2007. - 63 с.
2. Кремлевский П.П. Расходомеры и счетчики количества. - Л.: Машиностроение, 1989. - 701 с.
3. Володин Ю.Г., Марфина О.П. Математическое моделирование пусковых режимов энергетических установок. - СПб.: Инфо-да, 2007. - 128 с.
4. Володин Ю.Г., Марфина О.П. Расчет коэффициентов трения и теплоотдачи при нестационарном неизотермическом течении несжимаемого газа в осесимметричных каналах // Известия ВУЗов. «Машиностроение», 2007, № 3. - С. 21-26.
5. Володин Ю.Г., Марфина О.П. Расчет нестационарного неизотермического течения несжимаемого газа в осесимметричных каналах энергетических установок // Известия ВУЗов. «Машиностроение», 2007, № 5. - С. 39-42.
6. Володин Ю.Г., Марфина О.П., Цветкович М.С. Газодинамические характеристики отводящих стволов дымовых труб // Известия КГАСУ, 2015, № 4 (34). - С. 236-242.
7. Вращающийся регенеративный воздухоподогреватель: пат. 2324114 Рос. Федерация. № 2006127340/06; заявл. 17.07.06; опубл. 10.05.08. Бюл. № 13. - 4 с.
8. Володин Ю.Г., Ханнанов Р.Р. Повышение эффективности работы котлоагрегатов типа ТГМ-84 путем дросселирования газового потока в регенеративный воздухоподогреватель // Приволжский научный вестник, 2014, № 6. — С. 21-23.
Volodin Yu.G. - candidate of technical science, reader
E-mail: [email protected]
Hannanov R.R. - post-graduate student
Kazan State University of Architecture and Engeneering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1
Mathematical modeling of the worker process in gas canal for the boiler aggregate TGM-84
Resume
One of the current trends in the energy-producing industry is the heat energy and resource saving in the operation of the technological equipment.
The article presents the mathematical model of the turbulent flow of the flue gases in the gas canal of the boiler TGM-84. The calculations realization at utilization of the program complex ANSYS showed. Mathematical modeling of the worker process in the gas canal bases on the three-measuring unsteady-state equations Navier-Stokes. At modeling realized an idea of the throttling flow of the smoke gases in the seat for the deferent gas piping over RVP. By the
redistribution of potential and kinetic energy flow received the local increase of the pressure in deferent gas piping. The optimal position of the throttle screen and its optimal size identified. The measured air suction on the working boiler before and after installing the throttle screen showed a marked reduction of the inflow air into the gas canal.
Keywords: the inflow air, the mathematical model, the throttling screen, the boiler TGM-84, the economizer of air RVP-54, the gas canal.
Reference list
1. Volodin Yu.G., Bogdanov A.N. Technical measurement and instrumentation. Calculation of narrowing devices. Study book. - Kazan: Kazan. Gos. Energy. Un-t, 2007. - 63 p.
2. Krevlevskii P.P. Flowmeters and counters of amounts. - L.: Mashinostroenie, 1989. -701 p.
3. Volodin Yu.G., Marfina O. P. Mathematical modeling of starting modes of power plants. - SPb.: Info-da, 2007. - 128 p
4. Volodin Yu.G., Marfina O. P. Calculation of the coefficient of friction and heat transfer in unsteady flow of an incompressible non-isothermal gas in axisymmetric channels // Izvestija VUZov «Mashinostroenie», 2007, № 3. - P. 21-26.
5. Volodin Yu.G., Marfina O. P. Calculation of non-stationary non-isothermal flow of incompressible gas in axisymmetric channels of power plants // Izvestija VUZov «Mashinostroenie», 2007, № 5. - P. 39-42.
6. Volodin Yu.G., Marfina O.P., Tsvetkovich M.S. Gas dynamic characteristics efferent trunks chimneys // Izvestija KGASU, 2015, № 4. - P. 236-242.
7. The rotary regenerative air heater: - pat. 2324114 Ros. Federatsia. № 2006127340/06; zaiavl. 17.07.06; publ. 10.05.08. Bul. № 13. - 4 p.
8. Volodin Yu.G., Hannanov R.R. Improving the efficiency of boiler units TGM-84 by throttling the gas flow in the regenerative air heater // Privolzhskii nauchnii vestnik, 2014, № 6. - P. 21-23.
