CC BY
146
МАШИНОСТРОЕНИЕ • ГИДРАВЛИЧЕСКИЕ МАШИНЫ, ГИДРОПНЕВМОАГРЕГАТЫ
УДК 621.6
Ф. Г. БАКИРОВ, Ю. М. АХМЕТОВ, А. А. СОЛОВЬЕВ, С.В.ГУРИН,
А. Ю. ПАРХИМОВИЧ
ОПЫТ РЕАЛИЗАЦИИ КВАЗИИЗОТЕРМИЧЕСКОГО РЕДУЦИРОВАНИЯ В ВИХРЕВЫХ РЕГУЛЯТОРАХ ДАВЛЕНИЯ ЭНЕРГЕТИЧЕСКИХ СИСТЕМ
Рассматривается актуальная проблема обмерзания магистральных газопроводов и выпадения кристаллогидратов при дросселировании природного газа на газораспределительных пунктах. Представлены обобщенные результаты работ по созданию квазиизотермических вихревых регуляторов давления газа. Приведены результаты численного моделирования внутренних гидравлических и термодинамических процессов. Природный газ; вихревая труба; дросселирование; математическая модель; численное моделирование; газодинамика;
термодинамика
В настоящее время транспортировка и хранение газовых сред осуществляется под большими давлениями, что позволяет в ограниченном объеме трубопровода или баллонной системы пропускать или хранить большее количество газа за счет его сжатия. Однако по условиям технологического процесса, как правило, требуются более низкие давления газа, что вызывает необходимость установки специальных устройств для редуцирования давления до потребительского уровня. Такие устройства, предназначенные для снижения давления в газовых сетях и поддержания входного или выходного давления газа на заданном уровне независимо от отбора газа потребителями, называются регуляторами давления газа [1].
В 2001 г. в УГАТУ в рамках выполнения работ по федеральной целевой программе «Интеграция» были начаты работы по разработке теоретических и методологических основ создания энергосырьевых комплексов на основе вихревого эффекта для утилизации энергии давления магистральных газопроводов. По заказу газоснабжающих организаций (ООО «Баштрансгаз» и ОАО «Тюменьмеж-райгаз») в отдельное направление работ были выделены разработка и создание вихревых регуляторов давления газа (ВРДГ), предназначенных для квазиизотермического редуцирования давления газа. Заинтересованность газоснабжающих организаций к созданию ВРДГ связана со снижением температуры природного газа при его редуцировании на газораспределительных станциях (ГРС) и га-
зорегуляторных пунктах (ГРП) при использовании существующих регуляторов давления, в которых для снижения давления газа создается переменное регулируемое гидравлическое сопротивление потоку газа. Наблюдаемое снижение температуры при дросселировании газа связано с проявлением дроссель-эффекта Джоуля-Томсона. Это приводит к переохлаждению газового потока, обмерзанию трубопроводов, арматуры и распределительных устройств, а также к образованию кристаллогидратов. Образовавшиеся гидраты создают пробки в трубопроводах и отлагаются на деталях регуляторов, что приводит к примерзанию затворов к седлам, закупориванию проходных сечений регуляторов давления, образованию гидратных пробок в импульсных линиях обвязки регуляторов, что снижает надежность работы и может привести к возникновению аварийных ситуаций.
Проведенный анализ характеристик получаемой температуры смешения результирующих потоков на выходах вихревой трубы по результатам экспериментов [3,4] показал, что при определенных условиях вполне вероятно получение процесса дросселирования, близкого к изотермическому. При проведении серии экспериментов по исследованию газодинамических процессов в вихревой трубе было подтверждено, что при смешении «холодного» и «горячего» потоков, получаемых после температурного и массового разделения, возможно получение более высокой температуры результирующего потока, чем при дроссе-
лировании, в соответствии с эффектом Джоуля-Томсона.
Полученные результаты легли в основу разработки первого экспериментального образца ВРДГ на базе вихревой трубы (ВТ) с полным внутренним смешением потока. Закручивающее устройство представляет собой регулируемый двухсопловой ввод прямоугольного сечения, максимальное открытие которого соответствует рекомендации [2], по выбору отношения ширины к высоте как 2:1.
Рис. 1. Экспериментальный образец вихревого регулятора давления газа
На рис. 1 представлен экспериментальный образец ВРДГ, а на рис. 2 — поперечный разрез данного регулятора. Описание регулятора, принцип его работы и полученные результаты представлены в ранее опубликованных работах авторов [5,7,8,10]. Все эксперименты проводились на газодинамическом стенде УГАТУ с использованием воздуха.
Обработка полученных результатов показала, что АТ — значение разницы между температурой потока на выходе регулятора и на входе в регулятор большее или равное нулю, наблюдается для данной конструкции регулятора только при входном давлении газа до 1,2 МПа. Также была выявлена зависимость температуры газа на выходе из регулятора от геометрии соплового ввода. С увеличением площади критического сечения соплового ввода наблюдалось снижение температуры газа на выходе из регулятора.
б
Рис. 3. Продольный и поперечный разрез ВРДГ для ГРП
Рис. 4. Экспериментальный образец ВРДГ для ГРП
С учетом полученных результатов и выявления характера влияния различных факторов на термогазодинамические характеристики регулятора, специально для использования на ГРП с входным давлением до МПа, был спроектирован ВРДГ с фазоотделени-ем [9]. Особенностью данного регулятора является возможность отделения несанкционированной влаги в поступающем газе, а также перепуск горячей части газа по трубкам для подогрева газа на входе и обогрева корпуса регулятора. Внешняя труба регулятора выполнена пластиковой, закручивающее устройство ВТ представляет собой односопловой ввод прямоугольного сечения, причем за максимально открытое положение было выбрано соотношение его ширины к высоте как 8:1. На рис. 3,а представлен продольный, а на рис. 3,б — поперечный разрез регу-
лятора. На рис. 4 представлен внешний вид регулятора.
Данный экспериментальный образец ВРДГ с фазоотделением прошел испытания на сжатом воздухе на газодинамичесокм стенде УГАТУ и на природном газе в ОАО «Тю-меньмежрайгаз» на резервной нитке ГРП «Ушаково». Результаты испытаний регулятора на воздухе и природном газе представлены на рис. 5.
Анализ показывает, что при работе ИВРД на сжатом воздухе наблюдается повышение температуры, при избыточном давлении
МПа, на С,
а на природном газе повышение температуры составляет С, от-
носительно температуры на входе в регулятор. Разница в степени подогрева на воздухе и на природном газе объясняется теплофизическими свойствами газа, что отмечается и другими исследователями [2, 4]. Увеличение степени открытия соплового ввода приводит к уменьшению значения , а рост давления от
0,1 до 0,4 МПа ведет к росту АТ.
Наряду с экспериментальными работами проведены работы по численному моделированию газодинамических процессов в вихревых регуляторах с помощью CFD пакета ^-mosFloWorks. Пакет CosmosFloWorks предназначен для моделирования трехмерных течений жидкости и газа в технических и природных объектах, а также визуализации этих течений методами компьютерной графики.
! !- ш -» - -П ’н КГ V— 1* Г“ " і -и-*я ■f
и ■ СІ. *. р_, t ■» » ч U' : — l_ # 4 * *
І —ч і* E r* h » "¥—
0 1 2 3 4 5 ДР-10МПа
* 2 «б. рег-ра - Qmin (бяоеьй ачапог) < об. рег-ра - tjma* (6ajo№ft -аналог)
—•—gw-Тіяич -----реї^'лиюра
• НіПмґіні'* фаій01АЄЛн1ЄЛЧ — •— ДвґІуіГнган *рі.й<?нь Гй ТЗ fwi
—— р —Гз.-ч.іс it»r- а;грыгп9 палії* (,'іг iv —Ц «ГішінезсІьа птррьнси паллхі'мни
— — Кіінии*г*ио сироїде положение Qinri -.*--Миі*іма^ЬНЄ открыто? попр-тїпмв (Jnin
■ Экс ічрг"іе«і4п&н0е данное для природного rasa
Моделируемые течения включают в себя стационарные и нестационарные, сжимаемые и несжимаемые потоки многокомпонентных текучих сред в каналах и/или вокруг тел, с учетом гравитации, пограничного слоя, в том числе с учетом шероховатости стенок, с теплообменом между текучей средой и твердым телом. Использование различных моделей турбулентности и адаптивной расчетной сетки позволяет моделировать сложные движения жидкости, включая течения с сильной закруткой и теплообмен между потоком и стенкой.
В COSMOSFloWorks движение и теплообмен текучей среды моделируется с помощью уравнений Навье-Стокса, описывающих в нестационарной постановке законы сохранения массы, импульса и энергии этой среды. Для моделирования турбулентных течений упомянутые уравнения Навье-Стокса осредня-ются по Рейнольдсу, т. е. используется усредненное по малому масштабу времени влияние турбулентности на параметры потока, а крупномасштабные временные изменения усредненных по малому масштабу времени составляющих газодинамических параметров потока (давления, скоростей, температуры) учитываются введением соответствующих производных по времени. В результате уравнения имеют дополнительные члены — напряжения по Рейнольдсу, а для замыкания этой системы уравнений используются уравнения переноса кинетической энергии турбулентности и ее диссипации в рамках модели турбулентности.
Уравнения сохранения массы, количества движения и энергии в Декартовой системе координат могут быть записаны следующим образом:
др
~dt
д(р-щ)
дхи
др ■ щ
т
щ
= 0:
дР
Tik) + ^ = S;: дх.;
д(р-Е) д .. о m
—ОІ— + ~дхї-'\\Р'Е+Р)"ик + (1к^тік-щ) =
= Sk • щ + Qu.
где — скорость жидкости, — плотность, — внешняя массовая сила на единицу массы, — полная энергия на единицу массы,
— тепловой источник на единицу объема,
— тензор вязкостных касательных напря-
жений, — утечка тепла, и индекс — указывает на суммирование выражений по всем трем координатным осям.
Следуя допущению Буссинеска, касательные напряжения для Ньютоновских жидкостей определяются в данной модели следующим образом:
дщ дій
т‘> = " ■ W +
2 дщ
3 ' п.п '
где ;
Sij — дельта-функция Кронекера, равная единице при , или нулю в любом дру-
гом случае;
— динамическая вязкость;
— турбулентная вязкость;
— турбулентная кинетическая энергия.
Для замыкания системы уравнений применяются уравнения переноса турбулентной кинетической энергии и уровня ее диссипации согласно модели. В модели турбулентности определяется с использованием двух основных свойств турбулентности — турбулентной кинетической энергии и турбулентной диссипации
cft-P■ к2
Здесь — показатель турбулентной вязкости, определяемый из выражения:
(}—0.025-К,,
)*■(!
20,5
Tlj
где RT = Rt = “
Vk
14
— расстояние от
стенки.
Два дополнительных уравнения переноса используются для описания турбулентной кинетической энергии и диссипации:
др ■ к д(р ■ iih • к)
Ot
дхк
д
Ох к
,, . fit ч дк .
Sk,
др ■ є д(р ■ iik • є)
dt
дхк
JL (( + BL\ дє ■
дхк'ШН crj' дхк
SF
где члены уравнения и определяются как
Sk = т§ • - p • є + щ • Рв,
Pi ■ С в ■ Рв)
Cs2-f2
р ■ Є
Здесь Рв представляет возникновение турбулентности вследствие воздействия выталкивающих сил и может быть записан как
рв = _ ^!_ 1 °р
а в р дх{'
где ^ — проекция вектора ускорения свободного падения в направлении ж*, а в — константа, равная .
Св — константа, определяемая следующим образом: при Рв>0, Св = 1, в других случаях .
/і = 1 + /2 = 1
0,05
//*
Константы ССе\, Се2, (Jk, &е определены эмпирически.
СУ1 = 0,09: Се1 = 1,44: Се2 = 1,92:
а к — 1: е — 1,3.
как
[Н,
°С
ОТ
Охк
к = 1,2,3.
Здесь стс — константа, равная ас = 0,9;
Рг — число Прандтля;
— теплоемкость при постоянном давлении.
Эти уравнения описывают и ламинарный и турбулентный потоки. Параметры и равны нулю для ламинарного течения. Теплообмен в жидкостях описывается уравнением сохранения энергии, где утечка тепла определяется в модели турбулентности.
Для описания потока в пристеночных областях используется модель пограничного слоя. Модель основана на так называемом MWF (модифицированные функции стенки) приближении. Модель используется для того, чтобы охарактеризовать ламинарный и турбулентный потоки близ стенок и описать переходы из ламинарного к турбулентному течению и наоборот. Модель предоставляет значения скоростей и температур в пограничных областях уравнениям сохранения.
Система уравнений Навье-Стокса должна быть дополнена выражениями, определяющими теплофизические свойства жидкости. Предусматривается моделирование потоков газа с плотностью, вязкостью, теплопроводностью и теплоемкостью как функций температуры.
Уравнение состояния: р = р(р, Т, у).
Для газовых потоков применяется уравнение состояния идеального газа:
Р
Р =
R ■ Т
Для учета сжимаемости газа, плотность может определяться по одной из следующих двух зависимостей:
1. Степенная зависимость
р = ро/ 1 - С • In
Диффузионная утечка тепла определяется р0;
В + Р В + Pq
где Ро — плотность жидкости при давлении
С, В — коэффициенты;
Ро,С,В зависят от температуры;
— расчетное давление.
2. Логарифмическая зависимость
Р = Ро/
В + Р В + Pq
где п — показатель степенной зависимости плотности от температуры.
Для моделирования использована твердотельная модель проточной части экспериментального образца регулятора, представленная на рис. 6. На стенках задается граничное условие «Адиабатическая стенка». На входе и выходе задавались граничные условия по статическим давлениям и температуре, соответствующие экспериментальным данным. Турбулентность моделировалась с помощью стандартной модели для высоких чисел Рейнольдса.
Рис.6. Расчетная модель проточной части регулятора
Особенностью представленного регулятора является наличие в вихревой трубе крестовины, состоящей из профилированного лопаточного венца, устройства разворота, и трубки возврата. Профилированные лопатки обеспечивают безударное выпрямление закрученного потока и его плавный ввод в устройство разворота, обеспечивающее торможение и разворот газа. Далее часть горячего потока по центральной трубке возврата перетекает к диафрагме, повышая температуру в сопловом сечении. Вторая часть потока возвращается к сопловому сечению по внешнему каналу, обтекая снаружи поверхность трубки возврата, и в отличие от внутреннего потока участвует в процессе теплообмена со свободным вихрем.
На рис. 7 с помощью векторов представлена картина течения газа в завихрителе, которая показывает наличие прямого свободного вихря, обратного вынужденного вихря и потока газа через трубку.
Рис. 7. Вектора скоростей и изолинии давлений на выходе из трубки возврата
Представленная картина распределения скоростей реализуется распределением давлений, представленным на рис. 8 изолиниями. Данный рисунок также подтверждает предположение многих исследователей [2,3,4,6,7] о наличии центральной зоны разряжения. Рис. 8 демонстрирует направление движения газа и изолинии давлений. На рис. 9 с помощью линий тока можно увидеть дробление вихревого движения в крестовине по межло-паточным секторам.
На рис. 10 представлено поле распределения полных температур по осевому сечению регулятора. Интерес представляет повышение температуры практически на 25 градусов в камере торможения регулятора за устройством разворота.
Рис. 8. Вектора скоростей и изолинии давлений в крестовине
25
20
15
10
-5
■ Г
+ 0. —■ ■■ 2 0. 4 0. 6 0 8
Поле температур в сечении после трубки возврата
Поле температур на выходе из диафрагмы
Рис. 11. Распределение полных температур как функция от безразмерного радиуса (г = 1 соответствует центру, г = 0 — стенке)
Вторая кривая представляет эпюру потока на выходе из диафрагмы, после полного перемешивания потоков. Данный график подтверждает возможность получения изотермического снижения давления газа, так как на входе в регулятор температура задавалась равной 0оС, а среднемассовая температура газа на выходе, судя по второй кривой, будет составлять порядка 0 С.
Рис. 10. Поле полной температуры
Рис. 12. Распределение числа Маха по длине канала
На рис. 11 представлена эпюра полной температуры газа. Первая кривая соответствует сечению на выходе из трубки возврата, в котором потоки еще не успели перемешаться.
Другой задачей было определение скоростей в сопловом сечении по длине подводящего канала. На рис. 12 представлен график изменения числа Маха по длине канала. На длине 0,015 м регулируемые клинья начина-
ют сужать подводной канал, разгоняя поток до критических скоростей. Самое узкое геометрическое сечение (срез сопла) находится на расстоянии 0,03 метра, однако, критика, за счет нарастания пограничного слоя, «садится» несколько раньше. После среза сопла поток расширяется до скоростей порядка 1,4 М, после чего практически сразу скорость потока падает, через серию косых скачков, и устанавливается на уровне 0,9 М. На данный момент изучается характер течения при различных положениях регулирующих клиньев, и влияние сверхзвуковой области течения на процесс энергетического разделения.
По результатам анализа экспериментальных данных с учетом полученных результатов численного моделирования разработан вихревой регулятор давления газа, для использования в составе газораспределительного пункта ВРДГ-3 [6]. В регуляторе имеется фазоотборное устройство для отбора несанкционированной влаги и трубка перепуска горячего потока на обогрев регулятора.
3D-вид разработанного регулятора представлен на рис. 13 и 14. Для сохранения постоянного положения критического сечения сопловой ввод выполнен двухсопловым с регулируемыми поворотными лопатками.
-ч. *
Рис. 13. Общий 3D-вид вихревого регулятора в сборе
Рис. 14. 3D-вид внутреннего устройства вихревого регулятора
Разрабатываемый вихревой регулятор представляет собой гидропневмоагрегат, в котором для поворота допаток предусмотрена пневмогидравлическая система регулирования. Она позволяет значительно упростить кинематику синхронизации движения поворотных лопаток и снизить влияние трения в механизме синхронизации.
На регулятор подготовлен полный комплект конструкторской документации. Разработанная конструкция регулятора предположительно позволит осуществлять квазиизо-термический процесс снижения давления газа на любых режимах его работы и обеспечит необходимую точность и надежность регулирования.
СПИСОК ЛИТЕРАТУРЫ
1. Плотников, В. М. Регуляторы давления газа / В. М. Плотников, В. А. Подрешетников,
A. П. Дроздов [и др.]. Л.: Недра, 1982. 125 с.
2. Меркулов, А. П. Вихревой эффект и его применение в технике / А. П. Меркулов. М.: Машиностроение, 1969. 292 с.
3. Мухутдинов, Р. Х. Эффективность внедрения вихревых аппаратов (применительно к нефтехимическим производствам) / Р. Х. Мухутдинов, Р. Я. Амиров, Л. Э. Альмеев [и др.]; под ред. Я. С. Амирова. Уфа: Реактив, 2001.347 с.
4. Пиралишвили, Ш. А. Вихревой эффект. Эксперимент, теория, технические решения / Ш. А. Пиралишвили, В. М. Поляев, М. Н. Сергеев ; под ред. А. И. Леонтьева. М. : УНПЦ «Энергомаш», 2000. 412 с.
5. Патент на изобр. № 2237918 РФ. Регулятор давления газа с положительной обратной связью (варианты) / Институт технологии и организации производства. Приор. от 19.05.2003.
6. Патент на изобр. № 2282285 РФ. Вихревой регулятор давления газа. Опубл. 27.08.2006. Бюл. № 24.
7. Русак, А. М. Редуцирование давления природного газа без устройств предварительного подогрева / А. М. Русак, В. А. Целищев,
B. Л. Юрьев [и др.] // Сб. науч. тр. IV Между-нар. совещ. по проблемам энергоаккумулирования и экологии в машиностроении, энергетике и на транспорте. М. : ИМАШ РАН, 2004.
C. 133-144.
8. Русак, А. М. Разработка изотермического регулятора давления для редуцирования магистрального газа без подогрева на ГРС / А. М. Русак, В. Л. Юрьев, Ю. М. Ахметов [и др.] ; под ред. В. Л. Юрьева // Наука — производству. 45 лет на пути технического прогресса : ежегодн. науч.-техн. сб. Уфа : Гилем, 2003. С. 78-83.
9. Русак, А. М. Использование особенностей вихревых течений для получения квазии-зотермического процесса дроселирования давления природного газа / А. М. Русак,
B. Л. Юрьев, Ю. М. Ахметов [и др.] // Проблемы машиноведения и критических технологий в машиностроительном комплексе РБ : сб. науч. тр. АН Рб. Уфа: Гилем, 2005.
C. 37-49.
10. Соловьев, А. А. Исследование возможности получения изотермического процесса дросселировании в вихревом регуляторе давления газа / А. А. Соловьев, С. В. Гурин // Вестник УГАТУ. 2006. Т. 8, № 1 (17). С. 3-6.
ОБ АВТОРАХ
Бакиров Федор Гайфулло-вич, проф., зав. каф. теории авиац. и ракетн. двиг. Дипл. инж.-мех. (УАИ, 1969). Др техн. наук по тепл. двиг. ЛА (УГАТУ, 1995). Иссл. в обл. авиадвигателестроения, ракетн. двиг., энергоустановок, проблем энергосбережения.
Ахметов Юрий Мавлюто-вич, доц. каф. прикл. гидромех., зам. ген. дир. НИИТ. Дипл. инж.-мех. по авиац. двиг. (УАИ, 1959). Канд. техн. наук по тепл. двиг. (МАИ, 1978). Иссл. в обл. газогидр. течений и систем упр. энерг. установок.
Соловьев Алексей Александрович, асп. каф. авиац.
теплотехн. и теплоэнергет. Дипл. магистр техн. и технол. по энергомашиностроению (УГАТУ, 2005). Готовит дис. в обл. газодин. и термодин. закруч. потоков.
Гурин Сергей Владимирович, асп. каф. прикл. гидромех. Дипл. инж.-мех. по тепл. электрич. станциям (УГАТУ, 2004). Готовит дис. в обл. энергосберег. технол.
Пархимович Александр Юрьевич, асп. каф. прикл. гидромех. Дипл. инж.-мех. по гидравл. машинам (УГАТУ, 2004). Готовит дис. в обл. вихревых течений.
