ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2018
• ФИЗИКА •
Вып. 3 (41)
УДК 532.5 PACS 47.27. nf
Расчет охлаждения природного газа
в вихревом
А. В. Блиновa, С. Н. Пещеренкоa, М. П. Пещеренкоb
a Пермский национальный исследовательский политехнический университет 614039, Пермь, Комсомольский проспект, 29а email: [email protected]
b АО Новомет-Пермь, 614056, ул. Шоссе Космонавтов, 395 email: [email protected]
Нефть состоит из смеси углеводородов, находящихся в жидком состоянии в пластовых условиях и в виде газожидкостной смеси на устье скважины. Газовую фазу нефти называют также попутным газом. Существует проблема его доставки потребителям, поскольку строительство газопроводов обычно экономически не оправдано из-за небольших объемов добычи попутного газа на нефтяных скважинах. Перспективным способом решения этой проблемы является сжижение газа в вихревых трубах и его транспортировка потребителям автомобильным транспортом. Вихревые трубы не имеют движущихся деталей и практически не требуют обслуживания. Исследования по разработке данной технологии стали проводиться относительно недавно и пока не привели к серийным изделиям. Основная трудность состоит в низкой точности инженерных расчетов рабочих характеристик вихревых труб, которые проводятся по полуэмпирическим методикам и требуют испытаний опытных образцов изделий. Решением этой проблемы может стать моделирование процессов охлаждения газа в вихревой трубе средствами вычислительной гидродинамики. В представленной работе реализован именно такой подход. В полной постановке задача о моделировании вихревой трубы включает расчет течения как внутри трубы, так и на выходе из холодного и горячего концов, что требует больших вычислительных ресурсов. Нам удалось свести задачу к моделированию течения только внутри вихревой трубы. Для этого мы продлевали отводящие трубопроводы из холодного и горячего концов так, чтобы давление газа на выходе из них было бы близко к атмосферному. Для тестирования методики использовали экспериментальные данные, полученные в работах других авторов, а также на опытном образце вихревой трубы, разработанной авторами. Вычисленные локальные характеристики течения, такие как поля скоростей (осевая и тангенциальная компоненты), давления и температуры имели тот же вид, как в опубликованных экспериментальных работах. Интегральная характеристика вихревой трубы - величина охлаждения, полученная при численном моделировании, отличалась от экспериментальных данных не более чем на 4-7%.
Ключевые слова: природный газ; сжижение; вихревая труба; численное моделирование.
Поступила в редакцию 06.07.2018; принята к опубликованию 21.08.2018
The numerical modeling of natural gas cooling in vortex tube
A. V. Blinova, S. N. Peshcherenkoa, M. P. Peshcherenkob
"Perm National Research Polytechnic University, Perm, email: [email protected] b Novomet-Perm, Perm, 614056, Perm email: [email protected]
© Блинов А. В., Пещеренко С. Н., Пещеренко М. П., 2018
распространяется на условиях лицензии
Creative Commons Attribution 4.0 International (CC BY 4.0).
Crude oil consists of a mixture of hydrocarbons as a liquid at reservoir conditions and as a gas-liquid mixture at the wellhead. The gas phase of crude oil is called associated gas. There is a great problem how to transfer it to consumers, because pipelining makes no economic sense on account of low gas production at oil wells. An advanced way to solve this problem is to liquefy gas in vortex tubes and then transport it to the consumers by vehicles. Vortex tubes have no moving parts, so they will have no maintenance service. Research on this application of vortex tubes have been started recently and have not yet led to serial products. The main difficulty is the low accuracy of engineering calculations of vortex tube performance, which are based on the semi-empirical methods and require prototype testing of the product. The best way to solve this problem is modeling the gas cooling processes in the vortex tube by means of computational fluid dynamics. Such methodology is suggested in this paper. The complete formulation of the vortex tube modeling includes the calculation of the flow inside the tube and at the cold and hot outflows, which requires large computational resources. We could reduce the problem to air flow modelling inside the vortex tube only. At the outflows we extended the pipeline to get a reasonable approximation for the output gas pressure being close to atmospheric one. Some simulations have been done for air cooling with vortex tubes which were described previously in other experimental works. Also we have performed and described our own experiment. The local flow characteristics such as velocity field (axial and tangential components), pressure and temperature, which were obtained in calculations, have been the same as for published experimental works. Integral characteristics of the vortex tube - the magnitude of cooling from numerical simulation have differed from experimental data not above than 4-7%.
Keywords: natural gas; liquefaction; vortex tube; numerical modelling.
Received 06.07.2018; accepted 21.08.2018 doi: 10.17072/1994-3598-2018-3-45-54
1. Введение
Попутный нефтяной газ состоит в основном из метана. Одним из способов его доставки потребителям состоит в его сжижении на скважине (объем уменьшается примерно в 500 раз) и транспортировке автомобильным транспортом. Сжижению предшествует первичная обработка, включающая процессы сушки и сепарации [1], поскольку он содержит воду и механические примеси [2].
В промышленности применяют два способа сжижения природного газа [3]. Первый способ (изобарический процесс): газ охлаждают внешним источником до температуры начала конденсации, затем конденсируют, отбирая тепло с помощью того же внешнего источника. При втором способе охлаждение получают за счет процесса, проходящего в самом газе: изоэнтальпического (дросселирование) или изоэнтропического (расширение в детандере). В вихревых трубах реализуется изо-энтальпический процесс. По сравнению с детандерами вихревые трубы имеют ряд преимуществ -простоту конструкции, отсутствие движущихся деталей, более высокую надежность и долговечность, а также меньшую себестоимость [4]. Схематически вихревая труба представлена на рис. 1.
На входе газ по касательной к цилиндрической поверхности втекает в вихревую камеру, откуда закрученный поток выходит в трубу, дальний конец которой имеет узкое кольцевое отверстие вдоль ее цилиндрической стенки. В трубе газ, движущийся вдоль цилиндрических стенок, сжимается и нагревается, Нагретый газ выходит через коль-
цевое отверстие. Другая часть газа, движущаяся вблизи оси трубы, после взаимодействия со стенкой контрольного клапана (см. рис. 1), расширяется, охлаждается и возвращается в вихревую камеру, а затем покидает ее через центральное отверстие.
Рис. 1. Вихревая труба
Для создания кольцевого зазора на горячем конце вихревой трубы обычно используют клапан в форме полного или усеченного конуса, положение которого регулируется. Величина зазора между клапаном и стенкой вихревой трубы влияет на соотношение массовых расходов холодного и горячего потоков. Основной величиной, характеризующей эффект охлаждения в вихревой трубе, является величина изменения температур:
ЛГ = Т _ Т
Л1 Тт 1с,
где Тпп - температура на входе в вихревую трубу, Тс - на выходе холодного газа. В настоящее время проектирование вихревых труб осуществляется на основании приближенных инженерных расчетов и экспериментов [5, 6].
Известные методики расчета вихревых труб средствами вычислительной гидродинамики преимущественно ориентированы на получение качественной картины, т.е. объяснение эффекта охла-
ждения [7-20]. Основным их недостатком является упрощение полей течения и граничных условий на выходе горячего потока.
В [9, 10], [14], [16-19] было упрощено поле скоростей газа, его считали 2D осесимметричным, не учитывался реально трехмерный характер течения. Средняя относительная погрешность была порядка 30%. Ее удалось уменьшить до 15%, проводя расчеты для 3D сектора вихревой трубы [8], [11, 12], [15].
Применяемый в этих расчетах способ задания граничных условий на выходе горячего потока является, по сути, подгонкой результатов расчета к экспериментальным данным. Например, в [5], [16] значение давления на выходе горячего потока брали из экспериментальных данных, а в [8-13], [1719] давление на выходе горячего потока варьировали до достижения экспериментального значения расхода.
Целью данной работы стала разработка методики расчета вихревых труб, лишенной перечисленных недостатков, и ее тестирование, как на опубликованных экспериментальных данных, так и данных испытания опытного образца вихревой трубы, предназначенной для охлаждения попутного нефтяного газа. Ниже приведены: объяснение выбора модели течения (модели турбулентности) и граничных условий, описание расчетной сетки, полученных результатов расчетов и их сопоставления с экспериментальными данными.
2. Методика численного расчета вихревых труб
2.1. Выбор модели турбулентности
Типичное значение числа Рейнольдса в вихревых трубах [1], [6], [7] составляет порядка 104, что свидетельствует о развитом турбулентном характере течения. Модели турбулентности, применяемые в современной вычислительной гидродинамике, имеют полуэмпирический характер (содержат подгоночные параметры) и оптимизированы под определенные типы течений. Так, однопараметри-ческая модель Спаларта-Алламараса создавалась для расчета обтеканий профилей и течений в каналах переменного сечения, двухпараметрическая k - е - для течений с отрывами, а ее модификации k - т - для расчета пристенных течений, SST - для задач обтекания и тепловых расчетов.
Согласно [21] расчеты по стандартной k - е модели приводят к значительному расхождению с экспериментальными данными для течения жидкости через отверстия (в нашей задаче это области входа газа в вихревую трубу, а также выходы холодного и горячего потоков). В [21] была предложена модификация k - е модели: выведено новое уравнение для е, а коэффициент пропорциональности СИ в выражении для турбулентной вязкости V/;
(который в стандартной к - е модели является константой) стал функцией разности градиентов скоростей. Эта модель, называемая realizable к - е, реализована в пакете ANSYS Fluent, который и был применен нами для расчета вихревых труб. Приведем ее краткое описание модели течения.
В основной части расчетной области, кроме пристеночной, система уравнений имеет вид:
дй д I__\
- + — {uiuj ) =
dt дх..
1 dp д ---£_ + —
р дх,. дх,
(
(v + vt)
(
дй ди —L + —-
дх дх v j ' У
-2
3 j
Л
=0, v=с, kl,
дх И s
дк д
— +-(ки ) =
дt дх, 1
дх.
((
v+-
vv
дк
к У
дх
(
+ v
ди
ди —- + —
дх дх
V j ' У
дuj дх
дs д , _ ч
— +-(sui) =
дt дх, 1
v + -
ук у
дх..
j У
+ CJ2s
( дй дй Л2 - + 1
дх, дх,
V 1 ' У
- с
к+-jvss
дt
^ 77 2 Л
pT - p+^ p р 2
( (
дх,
и
дх
< 1 р
V,PCP
Pr
i
V v д T дх
Cp (T - Tef ) + y
уУ
( дйй. дй ^ + —^
дх дх V j ' У у
где и1 - средняя скорость, р - среднее давление, С - функция разности градиентов скоростей, связывающая турбулентную кинематическую вязкость V; с турбулентной кинетической энергией k и скоростью диссипации е, С - функция средней скорости турбулентной кинетической энергии, С2 -функция средней скорости диссипации, М„ - молярная масса воздуха, Я - универсальная газовая постоянная, ср - удельная теплоемкость воздуха, Те^ = 300 К, х - коэффициент теплопроводности, Рг = 0.85 - турбулентное число Прандтля, которое определено из экспериментальных данных; ск = 1.1, ск = 1.2.
Выбор модели для расчета пристеночной области зависит от размера ячейки первого слоя сетки, который характеризуется безразмерным параметром у+:
s
+ л/К / р У =—-:
V
где у - размер ячейки, ^ - напряжение сдвига у стенки:
дЦ.
ду
у=0
РСТ ^ Ша 1п "
а Utan - средняя скорость вдоль стенки.
Если у+ <1, то пристеночная область определяется условием Яе < 200 [26], где
Яе у =
Течение жидкости в пристеночной области описывается той же системой уравнений, но из которой исключаются уравнения для е и VI. Скорость диссипации энергии находится из выражения
е =
к 3
укС-4 (1 - ехр (-Яе у/ 2кС-))
где к = 0.4187 - константа фон Кармана. Турбулентная вязкость VI вычисляется так:
к2
vt=^cм-+
е
+ (1 + Хе)С-- 1/4ук(1 - ехр(-Яеу/70))у/к,
4 = 1 е 2
(
1 + Ш
'' Яе^ - 200 Л |Яе - 2001
I у \
\
агсШ0.98
На выходе горячего потока задавали геометрию проточных каналов в полном соответствии с экспериментальной установкой. Учитывались все пазы, стыки датчиков и втулок для шлангов. В качестве примера на рис. 2 приведена геометрия расчетной области для предложенного нами опытного образца вихревой трубы.
Если у+ >1, то течение в пристеночной области описывается пристеночными функциями: скорость и температура вычисляются по эмпирическим функциям, зависящим от турбулентной кинетической энергии k.
При выборе численной аппроксимации дифференциальных уравнений модели взяли второй порядок аппроксимации производных в уравнениях Рейнольдса и первый порядок в уравнениях для k и е. Вычисления проводили на 8 ядрах (2 ГГц) и 32 Гб оперативной памяти. Время одного расчета было порядка 8 ч.
2.2. Граничные условия и модификация геометрии расчетной области
В соответствии с типичными условиями эксперимента [21, 22] на входе в вихревую трубу задавали расход или полное давление Pm и температуру Тт.
Для холодного потока, чтобы не моделировать затопленную струю газа на выходе, удлиняли выходной патрубок так, чтобы давление на его выходе было близко к атмосферному.
Рис. 2. Геометрия расчетной области
2.3. Построение расчетной сетки
Объемные элементы имели форму тетраэдров. Погранслой содержал 5 слоев призматических элементов (с треугольником в основании). Высота элементов росла в направлении от стенки в объем с коэффициентом роста 1.2.
Качество сетки характеризовали коэффициентом асимметрии SQ, который показывает, насколько форма ячейки отличается от формы соответствующего ей правильного многоугольника:
SQ = тах
в -в в -в .
шах е . е шт
180°- в
в
где вe - угол между смежными ребрами грани правильной ячейки (60°для граней тетраэдров и основания призм, 90° для боковых сторон призм), 6шп и 6шах - минимальное и максимальное значение угла между смежными ребрами ячеек сетки. SQ может принимать значения от нуля до единицы, чем ближе значение коэффициента к нулю, тем лучше качество сетки. Качество сетки считается низким, если SQ > 0.98. В наших расчетах SQ был меньше 0.8.
Размер элемента объемной сетки во входном патрубке и выходном для холодного потока не превышал 0.8 мм, внутри вихревой трубы и в области выхода горячего потока - 2 мм.
В кольцевом зазоре между контрольным клапаном и вихревой трубой сетка состояла не менее чем из 4 элементов, а размер ячейки не превышал 0.15 мм.
К =
V
V
3. Тестирование методики
Поскольку в настоящее время нет доступных экспериментальных данных по применению вихревых труб для охлаждения природного газа, тестировали предложенную методику на экспериментальных данных, полученных при охлаждении воздуха и взятых из [20], [23].
3.1. Сравнение с экспериментальными данными
Согласно [20] вихревая труба состояла из следующих элементов: входного патрубка (диаметра Di„), вихревой камеры (диаметра Dvc и длины Lvc), выходного патрубка на выходе охлажденного потока (диаметра Dc), контрольного клапана на выходе нагретого потока в виде конуса с углом в и кольцевым зазором шириной z. Значения геометрических параметров приведены в табл. 1.
Испытания в статье [20] проводили при фиксированном массовом расходе на входе Qin = 8.4 г/с. Давление на холодном и горячем выходах вихревой трубы было равно атмосферному.
Расчетная сетка состояла из 1.5 млн. элементов. После 1850 итераций решение стало устойчивым. Сравнение результатов расчета температур T и расходов Q с экспериментальными данными приведено в табл. 2 при Tin=293.7 K, Pin=5.2 атм, Qin=8.4 г/с. Индекс 'in' относится к входному потоку, 'с' - к холодному выходу, 'А' - к горячему выходу. Видно, что расчетные данные согласуются с экспериментом [20] с погрешностью, не превышающей 5%.
По сравнению с [20] конструкция вихревой трубы в [23] имеет отличия: диаметр вихревой камеры равен диаметру трубы, а входной патрубок (диаметра Din) был не один, как в конструкции [20], а два. Давление на входе в вихревую трубу Pin = 10 атм.
Значения геометрических параметров вихревой трубы приведены в табл. 3.
Расчетная сетка состояла из 1.6 млн элементов. Устойчивое решение было после 1600 итераций. Эксперимент показал, что охлаждение воздуха было на 15.0 град, в расчете - на 16.1 град, отличие 7%.
3.2. По данным испытаний опытного образца вихревой трубы
3.2.1. Проведение эксперимента
Был изготовлен и испытан опытный образец вихревой трубы, схема которого приведена на рис. 3. Входной патрубок имел прямоугольное сечение L1xL2. Геометрические размеры указаны в табл. 4.
При проведении испытаний варьировали давление воздуха в подводящем патрубке вихревой трубы в интервале от 2 до 5 атм. Из холодного и горячего концов трубы воздух выходил в атмосферу. В установившемся режиме измеряли значения давления, температуры и массового расхода на входе и обоих выходах из вихревой трубы.
Рис. 3. Конструкция опытного образца вихревой трубы
Din, мм Dvc, мм Dvt, мм Dc, мм
3 30 16 12
Lvt, мм Lvc, мм в, град z, мм
205 4 30 0.4
Табл. 2. Сравнение экспериментальных результатов [20] и расчета
Tc, K Th, K Qc, г/с Qh, г/с
Эксперимент 275.1 303.8 2.71 5.64
Расчет 271.2 300.2 2.65 5.75
Табл. 3. Размеры вихревой трубы [23]
Din, Dvc Dvh Lvh Dc, в,
мм мм мм мм град
12 24 240 10 30
z, мм
0.5
Табл. 4. Размеры вихревой трубы
Lb мм L2, мм Dvc, мм Dvt, мм Dc, мм
6 4,8 42 10 5
Lvb мм Lc, мм в, град z, мм Dk, мм
220 130 30 0.5 5
Испытания проводили в два этапа. На первом этапе при заданном расходе воздуха Qin = 1040 г/сут и температуре входного потока Тп = 18.9°С меняли положение контрольного клапана так, что ширина зазора г варьировалась в интервале 0.1-1 мм, и измеряли величину охлаждения ДТ = ТПп-Тс. Целью испытания был поиск такого г, при котором эффект охлаждения достигал максимума. Результаты приведены на рис. 4, откуда видно, что максимум ДТ был при г ~ 0.5 мм.
Рис. 4. Влияние ширины кольцевого зазора г на величину охлаждения АТ
конца трубы: Ох и О - массовые потоки газа на выходе из холодного конца трубы и на входе в трубу. Полученные результаты приведены на рис. 7. Максимальное охлаждение было достигнуто при ^ ~ 0.2. Такое положение точки максимума функции ДТ(м) типично для вихревых труб [3].
При сравнении экспериментальных и расчетных данных мы сопоставляли только интегральные характеристики течения в вихревых трубах, поскольку в доступных нам экспериментах измерялись только интегральные характеристики. Теперь проанализируем локальную структуру течения в вихревых трубах, используя полученные расчетные данные.
Рис. 5. Зависимость величины охлаждения AT от входного давления Pin при z=0.5 мм, сплошная линия - эксперимент, пунктирная - расчет
3.2.2. Результаты расчетов
Расчетная сетка состояла из 6.2 млн элементов. Решение установилось после 2500 итераций. Результаты расчета AT(Pin) приведены на рис. 5. Видно, что отличие экспериментальных и расчетных значений было не более 4%, а абсолютная ошибка - не более 0.5°C.
Охлаждение газа AT увеличивается с ростом Pin, практическое значение имеет максимальное значение AT для конкретной вихревой трубы. Результаты расчета (рис. 6) показали, что для данного опытного образца AT ~ 60°. Максимум охлаждения достигается при Pin ~ 13-15 атм и при дальнейшем повышении давления AT не меняется.
К основным рабочим характеристикам вихревых труб относится также зависимость AT(m), где ¡г = Gx / G - доля газа выходящего из холодного
Рис. 6. Расчетная зависимость величины охлаждения ДТ от входного давления Рп при г = 0.5 мм
Рис. 7. Зависимость величины охлаждения от доли газа [ вытекающего из холодного конца трубы, сплошная линия - эксперимент, пунктирная - расчет
На рис. 8 и 9 приведены типичные поля давлений и температур в сечении перпендикулярном оси трубы и проходящем через входной патрубок.
Видно, что оба поля приближенно центрально симметричные. Давление и температура принимают максимальное значение на стенках трубы, минимальное - в центре трубы.
Рис. 8. Распределение давления в поперечной плоскости
Рис. 9. Распределение статической температуры в поперечной плоскости
На рис. 10 показано распределение температуры в поперечном сечении вихревой трубы. Видно, что зона охлаждения газа занимает относительно небольшой объем вихревой трубы и расположена в непосредственной близости к вихревой камере.
Для количественного анализа поля скоростей строили сечения перпендикулярные оси трубы на расстоянии 1, 3, 5, 10 и 20 диаметров трубы (обозначены 1Б, 3Б, 5Б, ...) вниз по основному потоку: от вихревой камеры к горячему концу трубы. В каждом сечении компоненты скорости (осевую и2, тангенциальную иТ, радиальную иг осредняли в окружном направлении и получали их зависимости от безразмерного радиуса г.
80
=. 40
н
о о а. о
2 -40
ф о О
-80
\ ч зсКЛ" ч. \ ч Ч ч \ ч
■-........100 - 20Й ч ч^^ ч \ ч * ' V
— — — // я/ V
0 0.2 0.4 0.6 0.8 1 Безразмерный радиус
Рис. 11. Зависимость осевой компоненты скорости от радиуса на расстояниях от вихревой камеры 1Б, 3Б, 5Б, ¡дБ, 20Б
Рис. 10. Распределение статической температуры в продольной плоскости
Рис. 12. Направление осевой скорости (схема)
На рис. 11 приведена зависимость осевой компоненты скорости и2(г). Видно, что в сечениях, находящихся на расстоянии 1Б, 3Б и 5Б от вихревой камеры, и2 меняет знак: на малых радиусах она положительна и направлена в сторону холодного конца трубы, на больших - в противоположную сторону. Схематично изменение направления и2 в зависимости от радиуса и расстояния от вихревой камеры показано на рис. 12. Снова, как и на рис. 10, видим, что область формирования охлажденного потока локализована вблизи вихревой
камеры. В осевом направлении она начинается не более чем на 5Б от вихревой камеры и имеет радиус не более 0.4 от радиуса трубы или 0.2Б.
Тангенциальная компонента скорости ит немонотонно зависит от г и имеет локальный максимум при г ~ 0.7 (см. рис. 13). Как ит, так и иг уменьшаются по мере увеличения расстояния от вихревой камеры. В экспериментах [5] наблюдали зависимости, аналогичные приведенным на рис. 11 и 13.
Безразмерный радиус
Рис. 13. Зависимость тангенциальной компоненты скорости от радиуса на расстояниях от вихревой камеры 1Б, 3Б, 5Б, ¡ОБ, 20Б
4. Выводы
Предложен способ расчета рабочего процесса вихревых труб методами вычислительной гидродинамики. Сравнение результатов расчета с известными экспериментальными данными показало, что отличие не превышало 4-7%.
Разработан опытный образец вихревой трубы. Результаты его расчета и стендовых испытаний отличались не более чем на 4 %.
Показано, что вычисленные поля скоростей и давлений имеют качественно такой же вид, что и в известных из литературы экспериментальных исследованиях [5].
Список литературы
1. Бирюк В. В., Веретенников С. В., Гурьянов А. И., Пиралишвили Ш. А. Вихревой эффект. Технические приложения. Т.2 М.:Научтехлитиздат, 2014. 490 с.
2. Катц Д. Л., Корнелл Д., Кобаяши Р., Поеттманн Ф. Х., Вери Д. А., Еленбаас Д. Р., Уайнауг Ч. Ф. Руководство по добыче, транспорту и переработке природного газа. М.: Недра, 1965. 677 с.
3. Федорова Е. Б. Современное состояние и развитие мировой индустрии сжиженного природного газа: технологии и оборудование. М.: РГУ нефти и газа им. И.М. Губкина, 2011. 159 с.
4. Жидков М. А., Гусев А. П., Рябов А. П., Вшивцев С. П., Кошовец Н. В., Коляко Б. Г., Соболь Е. П., Староконев Г. С. Трехпоточная вихревая труба - эффективное газодинамическое устройство для подготовки природного газа к транспорту // Нефть, Газ & Энергетика. 2007. Вып. 6. С. 18-22.
5. Меркулов А. П. Вихревой эффект и его применение в технике. М.: Машиностроение, 1969. 185 с.
6. Пиралишвили Ш. А. Вихревой эффект. Физическое явление, эксперимент, теоретическое моделирование. Т. 1. М.:Научтехлитиздат, 2013. 339 с.
7. Казанцева О. В., Пиралишвили Ш. А., Фузее-ва А. А. Численное моделирование закрученных течений в вихревых трубах // Теплофизика высоких температур 2005. Т. 43. № 4. С. 606-611.
8. Akhesmeh S., Pourmahmoud N., Sedgi H. Numerical study of the temperature separation in the Ranque-Hilsh vortex tube // American Journal of Engineering and Applied Sciences. 2008. Vol. 1. № 3. P. 181-187. DOI: 10.3844/ajeassp.2008.181.187
9. Aljuwayhel N. F., Nellis G. F., Klein S. A. Parametric and internal study of the vortex tube using a CFD model // International Journal of Refrigeration. 2005. Vol. 28. P. 442-450.
DOI: 10.1016/j.ijrefrig.2004.04.004
10. Azizi S. H., Andalibi M. R., Saleh S. R., Kahrom M. Numerical simulation of different turbulence models aiming at predicting the flow and temperature separation in a Ranque-Hilsh vortex tube // Thermal science. 2014. Vol. 18. N. 4. P. 1159-1171. DOI: 10.2298/TSCI110727201A
11. Behera U., Paul P. J., Kasthurirengan S., Karunanithi R., Ram S. N., Dinesh K. Jacob S. CFD analysis and experimental investigations towards optimizing the parameters of Ranque-Hilsh vortex tube // International Journal of Heat and Mass Transfer. 2005. Vol. 48. P. 1961-1973. DOI: 10.1016/j .ijheatmasstransfer.2004.12.046
12. Bramo A., Pourmahmoud N. A numerical study on the effect of length to diameter ratio and stagnation point on the performance of counter flow Ranque-Hilsh vortex tubes // Australian Journal of Basic and Applied Sciences. 2010. Vol. 4. N. 10. P. 49434957.
13. Deshmukh B. S., Chhapkhane N. K. Computational fluid dynamics analysis and experimental investigations of optimum geometry for the cold end orifice and snail entry of vortex tube. Journal of Mechanical and Civil Engineering. 2014. Vol. 11. N. 4. P. 2431. DOI: 10.9790/1684-11432431
14. Dutta T., Sinhamahapatra K. P., Bandyopadh-yay S. S. Comparison of different turbulence models in predicting the temperature separation in a Ranque-Hilsh vortex tube. International Journal of Refrigeration. 2010. Vol. 33. P. 783-792. DOI: 10.1016/j.ijrefrig.2009.12.014
15. Dutta T., Sinhamahapatra K.P., Bandyopadh-yay S. S. Numerical investigation of gas species and
energy separation in the Ranque-Hilsh vortex tube using real gas model // International Journal of Refrigeration. 2011. Vol. 34. P. 2118-2128. DOI: 10.1016/j.ijrefrig.2011.06.004.
16. Linhart J., Kalal M., Matas R. Numerical study of vortex tube properties // 16th international symposium on transport phenomena, Prague, 2005. p. 7.
17. Marques C. H., Isoldi L. A., dos Santos E. D., Rocha L. A. O. Constructial design of a vortex tube for several inlet stagnation pressures // Thermal Engineering. 2012. Vol. 11. N. 1. P. 85-92.
18. Pouraria H., Park W.-G. Numerical investigation on cooling performance of Ranque-Hilsh vortex tube // Thermal Science. 2014. Vol. 18. N. 4. P. 11731189. DOI: 10.2298/TSCI120610052P.
19. Pourmahmoud N., Akhesmeh S. Numerical investigation of the thermal separation in a vortex tube // World Academy of Science, Engineering and Technology. 2008. Vol. 19. P. 399-405.
20. Gao C. Experimental study on the Ranque-Hilsh vortex tube / PhD Dissertation in Mechanical Engineering. Eindhoven: Technische Universiteit, 2005. 151 p.
21. Shih T.-H, Liou W. W., Shabbir A., Yang Z., Zhu J. A new k-e eddy-viscosity model for high Reynolds number turbulent flows. Model development and validation // Computers and Fluids. 1995. Vol. 24. N. 3. P. 227-238.
22. Camire J. Experimental investigation of vortex tube concepts. / Master thesis in Applied Science. Vancouver: The University of British Columbia, 1995. 139 p.
23. Raut S. S., Dnyaneshwar N. G., Bhimate C. D., RautM. A., Uplarkar S. A., Patunkar P. P. An experimental modeling andinvestigation of change in working parameters on the performance of vortex tube // International Journal of Advanced Mechanical Engineering. 2014. Vol. 4. N. 3. P. 343-348.
References
1. Biryuk V. V., Veretennikov S. V., Guryanov A. I., Piralishvili Sh. A. Vikhrevoy effekt. Tekhnicheskiye prilozheniya (Vortex effect. Technical applications), vol. 2. Moscow: Nauchtekhlitizdat, 2014. 490 p. (In Russian).
2. Katz D. L., Cornell D., Kobayashi R., Poett-mann F. H.,Vary J. A., Elenbaas J. R., Weinaug C. E. Handbook of natural gasengineering. New-York: McGrawHill, 1959. 677 p.
3. Gronner J. Successful experience with vortex tube technology at the EPE cavity storage of RWE energy. 23rd World Gas Conference. Amsterdam. 2006. 11 p.
4. Zhidkov M. A., Gusev A. P., Ryabov A. P., Vshiv-tsev S. P., Koshovets N. V., Kolyako B. G., Sob-ol Ye. P., Starokonev G. S. Trekhpotochnaya vikhrevaya truba - effektivnoye gazodinamich-eskoye ustroystvo dlya podgotovki prirodnogo ga-
za k transportu (Three-flow vortex tube - efficient gas-dynamic tool for natural gas treatment). Neft, Gaz & Energetika (Oil, Gas and Energy), 2007, no. 6, pp. 18-22 (In Russian).
5. Merkulov A. P .Vikhrevoy effect i ego primeneniye v tekhnike (Vortex effect and its application for engineering). Moscow: Mashinostroyeniye, 1969. 185 p. (In Russian).
6. Piralishvili Sh. A. Vikhrevoy effekt. Fizicheskoye yavleniye, eksperiment, teoreticheskoye modeliro-vaniye (Vortex effect. Physical phenomenon, experiment, theoretical modelling), vol. 1. Moscow: Nauchtekhlitizdat, 2013. 339 p. (In Russian).
7. Kazanceva O. V., Piralishvili Sh. A., Fu-zeeva A. A. Chislennoe modelirovanie techeniy v vihrevih trubah (Numerical simulation of swirling flows in vortex tubes) Teplofizika vysokih temperatur (High temperature thermophysics), 2005, vol. 43, no. 4, pp. 606-611. (In Russian).
8. Akhesmeh S., Pourmahmoud N., Sedgi H. Numerical study of the temperature separation in the Ranque-Hilsh vortex tube. American Journal of Engineering and Applied Sciences, 2008, vol. 1, no. 3, pp. 181-187.
DOI: 10.3 844/aj eassp.2008.181.187
9. Aljuwayhel N. F., Nellis G. F., Klein S. A. Parametric and internal study of the vortex tube using a CFD model. International Journal of Refrigeration, 2005, vol. 28, pp. 442-450. DOI: 10.1016/j.ijrefrig.2004.04.004
10. Azizi S. H., Andalibi M. R., Saleh S. R., Kahrom M. Numerical simulation of different turbulence models aiming at predicting the flow and temperature separation in a Ranque-Hilsh vortex tube. Thermal science, 2014, vol. 18, no. 4, pp. 1159-1171. DOI: 10.2298/TSCI110727201A.
11. Behera U., Paul P. J., Kasthurirengan S., Karunanithi R., Ram S. N., Dinesh K. Jacob S. CFD analysis and experimental investigations towards optimizing the parameters of Ranque-Hilsh vortex tube. International Journal of Heat and Mass Transfer, 2005, vol. 48, pp. 1961-1973. DOI: 10.1016/j.ij heatmasstransfer.2004.12.046
12. Bramo A., Pourmahmoud N. A numerical study on the effect of length to diameter ratio and stagnation point on the performance of counter flow Ranque-Hilsh vortex tubes. Australian Journal of Basic and Applied Sciences, 2010, vol. 4, pp. 49434957.
13. Deshmukh B. S., Chhapkhane N. K. Computational fluid dynamics analysis and experimental investigations of optimum geometry for the cold end orifice and snail entry of vortex tube. Journal of mechanical and civil engineering, 2014, vol. 11, no 4, pp. 24-31. DOI: 10.9790/1684-11432431
14. Dutta T., Sinhamahapatra K. P., Bandyopadh-yay S. S. Comparison of different turbulence models in predicting the temperature separation in a Ranque-Hilsh vortex tube. International Journal
of Refrigeration, 2010, vol. 33, pp. 783-792. DOI: 10.1016/j.ijrefrig.2009.12.014
15. Dutta T., Sinhamahapatra K. P., Bandyopadh-yay S. S. Numerical investigation of gas species and energy separation in the Ranque-Hilsh vortex tube using real gas model. International Journal of Refrigeration, 2011, vol. 34, pp. 2118-2128. DOI: 10.1016/j .ijiefrig.2011.06.004
16. Linhart J., Kalal M., Matas R. Numerical study of vortex tube properties. 16th international symposium on transport phenomena, Prague, 2005. 7 p.
17. Marques C. H., Isoldi L. A., dos Santos E. D., Rocha L. A. O. Constructial design of a vortex tube for several inlet stagnation pressures. Thermal Engineering, 2012, vol. 11, no 1, pp. 85-92.
18. Pouraria H., Park W.-G. Numerical investigation on cooling performance of Ranque-Hilsh vortex tube. Thermal Science, 2014, vol. 18, no. 4, pp. 1173-1189. DOI: 10.2298/TSCI120610052P
19. Pourmahmoud N., Akhesmeh S. Numerical investigation of the thermal separation in a vortex tube.
World Academy of Science, Engineering and Technology, 2008, vol. 19, pp. 399-405.
20. Gao C. Experimental study on the Ranque-Hilsh vortex tube. PhD Dissertation in Mechanical Engineering. Eindhoven: Technische Universiteit, 2005, 151 p.
21. Shih T.-H., Liou W. W., Shabbir A., Yang Z., Zhu J. A new k-e eddy-viscosity model for high Reynolds number turbulent flows. Model development and validation. Computers and Fluids, 1995, vol. 24, no. 3, pp. 227-238.
22. Camire J. Experimental investigation of vortex tube concepts. Master thesis in Applied Science. Vancouver: the University of British Columbia, 1995, 139 p.
23. Raut S. S., Dnyaneshwar N. G., Bhimate C. D., Raut M. A., Uplarkar S. A., Patunkar P. P. An experimental modeling andinvestigation of change in working parameters on the performance of vortex tube. International Journal of Advanced Mechanical Engineering, 2014, vol. 4, no. 3, pp. 343-348.
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Блинов А. В., Пещеренко С. Н., Пещеренко М. П. Расчет охлаждения природного газа в вихревой трубе // Вестник Пермского университета. Физика. 2018. № 3 (41). С. 45-54. doi: 10.17072/1994-3598-2018-3-4554
Please cite this article in English as:
Blinov A. V., Peshcherenko S. N., Peshcherenko M. P. The numerical modeling of natural gas cooling in vortex tube. Bulletin of Perm University. Physics, 2018, no. 3 (41), pp. 45-54. doi: 10.17072/1994-3598-2018-3-45-54