CC BY
88
УДК 629.1
Э. И. Зангиров, М. В. Мухаметов, А. В. Свистунов, А. А. Ниндина
ИДЕНТИФИКАЦИЯ ТЕМПЕРАТУРНЫХ ХАРАКТЕРИСТИК КВАЗИИЗОТЕРМИНЕСКОГО ВИХРЕВОГО
РЕГУЛЯТОРА ДАВЛЕНИЯ ГАЗОВ
Рассматривается актуальная проблема получения рабочей характеристики квазиизотермического вихревого регулятора давления газов. Приведены результаты численного и физического моделирования. Представлены результаты их анализа и вывод линериализованного рабочего уравнения регулятора. Природный газ; вихревая труба; дросселирование; математическая модель; экспериментальные исследования; численное моделирование
ВВЕДЕНИЕ
Повышение энергетической эффективности трубопроводных систем распределения и транспортировки природного газа является приоритетным направлением развития энергетики и энергетической эффективности промышленности Российской Федерации и входит в перечень критических технологий, определяющих процесс создания энергосберегающих систем транспортировки, распределения и потребления энергетических ресурсов.
При дросселировании газа происходит снижение его температуры, что приводит к отложению твердых кристаллогидратов в проточной части дросселирующей аппаратуры (регуляторов давления газа), вследствие чего она перестает работать, что может привести к созданию аварийных ситуаций и полной остановке газораспределительной станции (ГРС).
Для проектирования дросселирующих устройств с применением современных технологий вычислительной техники необходима разработка математической модели, описывающей газодинамические и термодинамические процессы, происходящие в газовых потоках, а также построение рабочей характеристики, позволяющей прогнозировать работу устройства на всех режимах.
Процесс дросселирования газа в вихревой трубе, исследуемый с 1937 г., обладает рядом особенностей. Рассмотрены и проанализированы научные труды, публикации, результаты
Контактная информация: 8 (347) 273-09-44 Исследование выполнено при поддержке Министерства образования и науки Российской Федерации, соглашение 14.B37.21.0337 «Научное обоснование, создание и исследование энергосберегающих вихревых технологий фазоразделения, подогрева и редуцирования при транспортировке попутного и природного газа».
исследований, посвященные проблемам изучения и практического применения устройств, реализующих вихревой эффект таких ученых, как В. С. Мартыновский, В. П. Алексеев,
A. П. Меркулов, А. И. Леонтьев, Ш. А. Пира-лишвили, Ю. Д. Райский, Л. Е. Тункель,
B. М. Бродянский, А. Д. Суслов, Р. Я. Амиров, Р. Х. Мухутдинов, М. А. Жидков, А. М. Русак, Ю. М. Ахметов, В. А. Целищев, С. В. Гурин, А. А. Соловьев, А. Ю. Пархимович, из зарубежных современных исследователей вихревого эффекта можно выделить Шепера (Scheper G. V.), Ван Димтера (Van Deemter J. J.), Такахама (Takahama H.) и др.
В работах указанных авторов рассматривались процессы в вихревых устройствах (рис. 1), предусматривающих температурное разделение (температурную стратификацию) потока газа [1].
Prop, Тгор, Grop Рхоп, Тхощ взят
Рвх, Твх, Gbx Рис. 1. Противоточная вихревая труба
В вихревую трубу поток газа подается с параметрами Р Т В результате снижения давления газ перейдет в состояние с параметрами Р'
Т Дросселирование исходного потока в вихревой трубе позволяет получить на выходе поток газа с параметрами нагретого потока Ргор, Тгор и массовым расходом Огор и охлажденного потока РХОЛ, ТХОЛ с массовым расходом ОХОЛ. Полученные в результате прохождения вихревой трубы параметры газа значительно отличаются от параметров газа, полученных после прохождения им дросселирующих устройств (снижение и превышение температуры соответственно). При
РЕАЛИЗАЦИЯ ПРОЦЕССА
квазиизотермического
ДРОССЕЛИРОВАНИЯ ГАЗОВ
Исходя из равенства полных удельных энтальпий входящего в вихревую трубу потока газа и суммы истекающих (нагретого и охлажденного) из нее потоков газа, характеристикой квазиизотермического процесса является температура смешанного потока, а также перепад температуры по всей системе смешения потоков после стратификации в вихревой трубе при cp = const:
С = +(1 -цЖ . (1)
Анализ экспериментальных исследова-
ний противоточной вихревой трубы согласно А. П. Меркулову, представленный в виде графика в Р-Т координатах процесса дросселирования и перераспределения температуры потоков газа после вихревой трубы, изображен на рис. 3 [3].
Сравнительный анализ испытаний вихревых труб [4] позволил выявить зоны, в которых смесевая температура равна входной (зона изотермического дросселирования) либо превышает ее (зона квазиизотермического дросселирования). Качественная зависимость температуры газа по тракту вихревой трубы от доли охлажденного потока газа цх представлена на рис. 4.
Область изотермического дросселирования (Твых=Твх)
Рис. 3. Принципиальное сравнительное представление в Р-Т координатах процесса дросселирования давления газа при протекании газа через гидросопротивление (дроссель эффект Джоуля-Томсона) и процесса дросселирования в вихревой трубе с учетом температурной стратификации
смешении истекающих из трубы потоков газа логично получение потока газа с параметрами РдТ , ТДТ , однако, как было показано в работах Ю. М. Ахметова и С. В. Гурина, существует область, в которой возможна реализация квазиизо-термического процесса дросселирования.
В противоточной вихревой трубе с внешним смешением потоков (рис. 2) массовый расход газа на выходе GВЬIХ = GтОp + GХОЛ определится как сумма массовых расходов нагретого и охлажденного потоков. Интенсивность энергетического разделения газов в вихревой трубе (вихревой эффект), как известно, оценивают по зависимости избыточных величин температур
газа AT
: T* - T и AT
охлажденного потока н
гор
G,.
G,
T* - т* от доли
гор вх
относительно мас-
сового расхода на входе G Доля подогретого
потока при этом составит
( (( - ((x ) [2].
Рем, Тем, Gcm
.....
I
Рхол, Тхол, Gxoj
Ргор, Тгор, Сгор
Рвх, Твх, Овх Рис. 2. Противоточная вихревая труба с внешним смешением потоков
Как видно, смесевая температура газа Тсм на всем участке смешения превышает температуру, полученную при снижении давления газа путем его дросселирования (Дроссель-эффект).
Исследование температурных характеристик, получаемых в результате дросселирования газа в вихревых трубах, подтверждают возможность получения квазиизотермического процесса, для осуществления которого необходимо создание условий интенсивной температурной стратификации входного потока и последующего смешения газа в определенных массовых пропорциях по формуле (1) с максимальным исключением теплообменных процессов между стратифицированным потоком и окружающей средой.
Анализ механизма процесса температурной стратификации в вихревых трубах позволил выявить несоответствие в процессе дросселирования дифференциального дроссель-эффекта (эффекта Джоуля-Томсона) в двух областях: в области потока с низкой температурой и в области потока с высокой температурой. Область, в которой температура потока превышает температуру входа, можно использовать в качестве источника предварительного подогрева для повышения входной температуры. Тогда подогрев входного газа, за счет горячего газа, полученного после стратификации, приводит к тому, что температура газа перед стратификацией непре-
рывно повышается, реализуя положительную обратную связь [5].
ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ КВРД
На основании проведенного анализа вихревых труб сотрудниками УГАТУ и ОАО «Научноисследовательский институт технологии» была разработана конструкция и изготовлен экспериментальный образец квазиизотермического регулятора давления газов (КВРД) (твердотельная модель представлена на рис. 5), позволяющий осуществлять дросселирование газа без снижения его температуры [6].
Рис. 5. Твердотельная модель экспериментального образца КВРД:
1 - корпус; 2 - винтовая головка; 3 - камера энергетического разделения; 4 - закручивающее устройство (улитка); 5 - регулирующие клинья;
6 - тормозное устройство; 7 - устройство возврата потока; 8 - диафрагма; 9 - диффузор
-------------------------- д(7г-7вх)
Зона
I
квазиизотермического процесса (км > №х)
Область изотермического процесса (Км = йх)
Область Дроссель-эффекта (эффект Джоуля-Т омсона)
Д (7х-7вх)
г 1
п
- Г 1
|<- >1
С 1 1 — и - л “ » 1 О
чч и к - е—£ ы 8 < о 4. н-8 і
і ^ / » ■- іщ
4 і о і Я 1
- - І к и > 1 1 і і 4
£ Д ^ их
О
0.2
0,4
0.6
0.8
Рис. 4. Результаты испытаний различных вихревых труб:
□ ОД - согласно экспериментам А. П. Меркулова;
О О О О О в - согласно экспериментам Р. X. Мухутдинова (разное входное давление)
Для прогнозирования работы подобных устройств и их проектирования на уровне современных технологий необходимо получение рабочей характеристики на всем диапазоне их работы (обобщенная рабочая характеристика).
Получение обобщенной рабочей характеристики КВРД осуществляется согласно следующим этапам:
• верификация результатов численного моделирования согласно результатам имеющихся экспериментальных исследований, введение при необходимости поправочных коэффициентов;
• проведение численного моделирования КВРД на всем диапазоне входных и выходных давлений;
• получение обобщенной рабочей характеристики КВРД.
Экспериментальные исследования образца выполнялись в лаборатории «Газодинамика высоких давлений» УГАТУ в 2011 г. В качестве рабочего тела использовался сжатый воздух. Абсолютное входное давление в экспериментах составляло 0,7 МПа. Выходное давление составляло 0,1 МПа. Температура воздуха на входе в регулятор в процессе экспериментов изменялась в интервале 289...278 К.
Схема измерений лабораторного стенда представлена на рис. 6.
г, к
ЗОЯ
Р, 0,1 МПа
Рис. 6. Базовая схема измерений режимных параметров КВРД:
Рвх, Твх - давление и температура потока газа на входе в КВРД; Рхол, Тхол - давление и температура потока газа перед закручивающим устройством КВРД; Тгор - температура газа перед устройством возврата потока; Р , Т - давление и
♦'Г Г 7 вых вых
температура потока газа на выходе из КВРД
Экспериментальное исследование, при котором давление входного потока газа составляло 0,7 МПа; давление потока газа на выходе
0,1 МПа; температура входного потока газа изменяется в диапазоне 292,3...277,8 К представлено на рис. 7.
35
30
25 Ли
-0-7хол
20 Ггор
1? -о-Гвых
-*-Ръх
10 -6-Р\0П
5 V Рпых
(? сек
Рис. 7. Экспериментальное исследование КВРД (Р = 0,7 МПа, Р = 0,1 МПа)
4 вх 7 7 вых 7 '
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ КВРД
Для проведения численного моделирования использовалась математическая модель, опубликованная ранее [7], с учетом влияния теплообмена в газе и стенках КВРД.
Теплопередача в газе описывается посредством следующей общей формы уравнений сохранения энтальпии:
5 . 7 . дР дР дР (2)
---(э Ии) = — + и — + т,.,. — > у ’
дх^ дг 1 дЛ - дЛ
где И - статическая энтальпия, определяемая уравнением:
И = СрТ - СРТо, (3)
где Т - температура; Ср - средняя удельная теплоемкость при постоянном давлении для температуры Т; с0р - удельная теплоемкость при постоянном давлении для температуры Т0.
Основное уравнение для полной энтальпии (Н) можно получить сложением уравнения сохранения механической энергии и уравнения статической энтальпии:
др# д . „ ч дРи- (4)
— +-----(р и ,Н + мл,,) =---(4)
дЛ дх} 1 111 дх,
где Н =1 • и и + И.
2 г г
Применение уравнения переноса полной энтальпии и тем более статической энтальпии рекомендуется осуществлять для потоков с числами Экерта (Ес) порядка 1,0 и более. Величина этого коэффициента определяется как:
и2
Ес =
с Т
р
(5)
где и - величина характеристической скорости, а Т - величина характеристической (абсолютной) температуры.
Наряду с моделированием диффузии тепла в текучей среде (с помощью уравнений Редлиха-
Квонга и Сазерленда), моделируется также теплопередача в твердых телах с помощью уравнения
^=—Ги—V в„. (6)
д,г дх.. { дх1.) н
где е = сТ - удельная теплоемкость, Т - температура, х . - теплопроводность, вн - удельное (в единице объема) тепловыделение источника тепла.
Конвективный теплообмен между поверхностью твердых тел и текучей средой моделируется при расчете пограничного слоя потока текучей среды. В результате в соответствии с законом Стефана-Больцмана излучаемое с единицы поверхности тепло определяется как
вн = е • а0 • (7)
где е - степень черноты поверхности, с0 - постоянная Стефана-Больцмана, Т - температура поверхности.
Учет влияния теплообмена каждой поверхности позволяет определить (с учетом рассчитываемого фактора видимости) разность между приходящими и уходящими (т.е. отраженными и излучаемыми) тепловыми потоками.
Твердотельная модель расчетной области КВРД для численного моделирования была создана в пакете трехмерного твердотельного моделирования БоНёЦЬ^.
Расчетная сетка имеет следующие параметры:
• количество ячеек равно 6515776;
• размер минимальной ячейки равен 0,1 мм3;
• размер максимальной ячейки равен 4,8 мм3.
Дополнительно в расчетной сетке, в областях с максимальными градиентами скорости и давления, применялось локальное измельчение сетки. Измельчение использовалось в области за срезом вводного сопла, в области за диафрагмой (рис. 8).
Случайная погрешность нивелируется количеством измерений и имеет значение на порядок ниже приборной, и по этой причине не учитывается.
Результат численного моделирования позволил визуализировать течение потока газа в наружном и внутреннем контурах КВРД, определить структуру и характер взаимодействия вихревых потоков (рис. 9). Это позволило определить картину теплового взаимодействия внутреннего и внешнего контуров, области подогрева и охлаждении стенок внутреннего контура и определить характер изменения параметров потока газа по тракту регулятора.
Рис. 9. Результат численного моделирования
Визуализация теплового взаимодействия контуров позволила определить направления теплообмена на различных участках корпуса внутреннего контура. Были четко определены характерные места подогрева и охлаждения внутреннего контура, что позволяет выработать конкретные рекомендации по организации теплообмена по контуру положительной обратной связи для получения квазиизотермических процессов дросселирования давления при проектировании регуляторов.
Результаты численного и физического моделирования представлены в таблице.
Результаты физического и численного моделирования
Результаты физического моделирования
Твх, Т 1 хол? Тгор, Твых, Р™ Р ± хол? р 1 ВЫХ?
к к К К МПа МПа МПа
292, 2 296,1 301,8 293,3 0,7 0,72 0,1
Результаты численного моделирования
Рис. 8. Расчетная сетка КВРД с локальным измельчением
Гвх, К Тхо^ К Ттоо, к твых, К
292,2 293 298,6 291,5
Сравнительный анализ полученных данных позволил выявить, что сходимость результатов численного и физического эксперимента по составленной математической модели не превышает погрешности 2%.
ЛИНЕАРИЗОВАННАЯ РАБОЧАЯ ХАРАКТЕРИСТИКА КВРД
В результате проведенной идентификации экспериментальных данных по испытаниям КВРД может быть получено математическое описание зависимости температуры газа на выходе регулятора от температуры газа на входе (т.е. Твых =ДТвх)). Для значений барометрических параметров потока газа Р = 0,7 МПа, Р = 0,1 МПа
А вх 7 7 вых 7
уравнение рабочего процесса КВРД будет иметь вид:
Т = (-0,0043 • (Т )2 + 0,0767 • (Т ) + 0,5223) X
вых вх вх
X (Твх) + (0,0363 • (Твх)2 -0,5657- (Тю) +
+ 9,5545) (3)
С учетом изменения барометрических параметров потока газа на входе и выходе КВРД обобщенное уравнение рабочего процесса КВРД Т = /(Т , Р , Р ) может быть представлено в
вых вх вх вых
виде:
Т = [к (Т ) • АТ + Ь (Т )] +
вых 14 вху вх 14 вху л
+ [ АР) • АР + ВД] - к3 • АР (4)
Предлагаемое уравнение учитывает:
1) влияние величины температуры на входе в регулятор Твх;
2) влияние величины давления на входе в регулятор Рвх;
3) влияние величины давления на выходе из регулятора Рвых;
4) влияние величины перепада АР = Р - Р
вх вых
давления на регуляторе;
ВЫВОДЫ
Разработанная математическая модель ква-зиизотермического вихревого регулятора давления позволяет с высокой точностью описывать газодинамические и термодинамические процессы, происходящие в газовых потоках.
В результате проведенной идентификации и верификации результатов численного моделирования получено эмпирическое уравнение
температурных параметров регулятора при неизменных барометрических параметрах.
Предложено (в общем виде) нелинейное уравнение температурных параметров регулятора с учетом изменения входных и выходных параметров потока газа.
СПИСОК ЛИТЕРАТУРЫ
1. Меркулов А. П. Вихревой эффект и его применение в технике. М: Машиностроение, 1969. 183 с.
2. Пиралишвили Ш. А. Вихревой эффект. Эксперимент, теория, технические решения / под ред. А. И. Леонтьева. М.: УНПЦ «Энергомаш», 2000. 412 с.
3. Соловьев А. А. Исследование возможности получения изотермического процесса при дросселировании в вихревом регуляторе давления газа // Вестник УГАТУ, Т. 8, № 1. Уфа, 2006. С. 3-6.
4. Свистунов А. В. и др. Исследование влияния тормозного устройства на структуру потока и параметры изотермического вихревого регулятора давления // Вестник УГАТУ. 2011. Т. 14, №4 (44). С. 149-153.
5. Пат. 2237918 РФ, МПК G05D16/00. Регулятор давления газа с положительной обратной связью (варианты) / Асадуллин М. З. и др. 2003114801/28; заявлено 19.05.2003; опубл. 10.10.2004. Бюл.28.(НИИТ)
6. Ахметов Ю. М. и др. Опыт реализации ква-зиизотермического редуцирования в вихревых регу-ляторх давления энергетических систем // Вестник УГАТУ 2007. Т. 9. №6 (24). С. 3-6.
7. Ахметов Ю. М., Пархимович А. Ю., Свистунов А. В. Численное моделирование процессов стратификации в изотермическом вихревом регуляторе с внутренним смешением // Вестник УГАТУ Т. 10, №14 (37), 2010.
ОБ АВТОРАХ
Зангиров Эрнест Ирекович, магистр той же каф. (УГАТУ, 2012). Дипл. бакалавра (УГАТУ, 2011). Иссл. в обл. газогидр. течений и систем упр. энерг. установок.
Мухаметов Мансур Вазирович, магистр той же каф. (УГАТУ, 2012). Дипл. бакалавра (УГАТУ, 2011). Иссл. в обл. многофазных течений.
Свистунов Антон Вячеславович, ассист. той же каф., маг. каф. прикл. гидромех. (УГАТУ, 2009). Дипл. бакалавра (УГАТУ, 2007). Канд. техн. наук по гидравл. машинам (УГАТУ, 2012). Иссл. в обл. га-зогидр. течений и систем упр. энерг. установок.
Чиндина Александра Александровна, магистр той же каф. (УГАТУ, 2011). Дипл. бакалавра (УГАТУ 2010). Иссл. в обл. тормозных устройств КВРД.
