CC BY
10Проектирование, сооружение и эксплуатация систем трубопроводного транспорта
УДК 519.63+533.6
УЧЕТ ВЛИЯНИЯ ЦЕНТРОБЕЖНОЙ СИЛЫ ПРИ ЧИСЛЕННОМ МОДЕЛИРОВАНИИ ВОСХОДЯЩЕГО ЗАКРУЧЕННОГО ПОТОКА ГАЗА
ACCOUNTING INFLUENCE OF CENTRIFUGAL FORCE IN THE NUMERICAL MODELING OF RISING SWIRLED GAS FLOW
С. П. Бау тин, А. Г. Обухов
S. P. Bautin, A. G. Obukhov
Уральский государственныйуниверситет путей сообщения, г. Екатеринбург
Тюменский государственный нефтегазовый университет, г. Тюмень
Ключевые слова: полная система уравнений Навье — Стокса: начальные и краевые условия: газодинамические характеристики; центробежная сила Key words: complete system ofNa\ier — Stokes equations; initial and boundary conditions; 'gas dynamic characteristics; centrifugal force
При исследовании сложных течений газа используются модели сжимаемой сплошной среды, основанные на численном решении полной системы уравнений Навье — Стокса [1, 2]. Такие модели наиболее адекватно описывают физические процессы в указанных течениях, поскольку последовательно учитывают диссипативные свойства сжимаемой сплошной среды - вязкости и теплопроводности.
Теоретические и численные исследования, проведенные в работах [2-6], подтвердили предложенную в [7] общую схему возникновения и последующего функционирования восходящего закрученного потока. Основная идея предложенной в [7] схемы возникновения восходящего закрученного потока заключается в том, что в результате локального прогрева поверхности суши или водной поверхности появляется восходящий поток воздуха. Замещающее его радиальное течение под действием силы инерции (силы Кориолиса) приобретает осевую закрутку.
В серии работ [8-10] предложенная схема получила экспериментальное подтверждение. В этих работах свободный вихрь инициировался нагревом снизу металлической круглой плиты пламенем газовой горелки. Появление радиального движения воздуха и его последующая осевая закрутка может инициироваться как локальным нагревом подстилающей поверхности, так и холодным вертикальным продувом. Последний способ получения восходящего закрученного потока был успешно реализован в лабораторных условиях [11].
При правильном выборе начальных и граничных условий [12, 13] численное решение полной системы уравнений Навье — Стокса позволило смоделировать возникно-
92
Нефть и газ
№ 4, 2015
вение и развитие сложных трехмерных нестационарных течении сжимаемого вязкого теплопроводного газа в восходящем закрученном потоке [14—18].
Во всех упомянутых выше исследованиях при численном моделировании предполагалась малой и не учитывалась центробежная сила, действующая на частицы газа как результат вращения Земли.
Целью данной работы является последовательный учет центробежной силы при численных расчетах газодинамических характеристик трехмерного нестационарного течения сжимаемого вязкого теплопроводного газа в восходящем закрученном потоке, вызванном вертикальным холодным продувом.
Полная система уравнений Навъе — Стокса без учета влияния центробежного силы.
При моделировании сложных течений упругой сплошной среды, обладающей дис-сипативными свойствами (вязкостью и теплопроводностью), используется полная система уравнений Навье — Стокса, которая будучи записанной в безразмерных переменных с учетом действия сил тяжести и Кориолиса в векторной форме имеет вид [14]
р(+У-Ур + р(ИуУ~ О,
Т „ 1 — V р + -
ур У V
1 -4
Т(+У-УТ + (у-\)ТЛуУ~ ' -АТ +
рАУ-Ц
р
2 р
{> О"
(1)
(н, - V. )" +(■'_, и, )" ] || 32 (а,. I V,. )" I ( II. + И ... }■ I (V. + " , )" _|]| :
где значения безразмерных коэффициентов вязкости и теплопроводности следующие: д, =0,001, к() «1,458333//,,.
Эта система в дифференциальной форме передает законы сохранения массы, импульса и энергии в движущейся сплошной среде.
В системе (1): /— время; х >'. г— декартовы координаты; р— плотность газа;
V = (и, V, — вектор скорости газа с проекциями на соответствующие декартовы оси; Т— температура газа; £ = — вектор ускорения силы тяжести, а
g = const> 0; У = (а\' -Ь\\>aii.hu} — вектор ускорения Кориолиса, где
а = 20зт ц/. Ь = 20соз ц/. О = |о|: О, — вектор угловой скорости вращения Земли; ц/— широта точки О — начала декартовой системы координат х\~2() , вращающейся вместе с Землей.
Полная система уравнений Навье — Стокса с учетом влияния силы инерции. В книгах [7, 14] приведена система уравнений газовой динамики при учете действия сил тяжести и Кориолиса, но без учета действия центробежной силы:
Рц6 =-/>Ох(Охг),
где г— радиус-вектор точки М в неподвижной системе координат с центром в точке С (рис. 1).
Рис. 1. Рассматриваемые системы координат
№ 4, 2015
Нефть и газ
" 93
Для записи этой силы в подвижной системе координат с центром в точке О имеем г = СМ = СО + ОМ = (0,0,Кзем) + (х, у, 2), где Кзем — радиус Земли.
Следовательно, в подвижной системе координат центробежная сила задается соот-
ношением
С учетом свойства векторного произведения
= axb +ахс
(2)
получаем нужную для вычислений формулу
Рф = -/?Qx(Qx(OA^))-/?Qx(Qx(x,v,z)). Каждое из слагаемых в формуле (2) вычисляется отдельно с учетом координат вектора Q в подвижной системе координат [6, 7] (см. рис. 1):
Q = (0, Qcosi//, Qsini//).
Последовательное вычисление векторных произведений для первого слагаемого из формулы (2) дает следующее:
1) Qx(0,0,R J
i j k 0 Qcosi// Qsini// 0 0 R ..
QR eos w i ■.
зем >
2) Q x (QR3eM cos ц/, 0,0):
i J к
0 Qcosi// Qsini// QR cosí// 0 0
зем г
1 .
Л
=: О- Кзем eos i//sin i//у - 0-'R3eM cos" у/к = 0;-sin2i//;-cos- ц/
■. 2 } Последовательное вычисление векторных произведений для второго слагаемого из формулы (2) дает следующее:
У к
1) Qx(x,v,z) = 0 Qcosi// Qsini// =
X у Z
= Q (z cos ц/ - у sin i//) i + Qx sin у/ j — Qx eos ц/ к ; ' j к
2) Q2
0 cosí// sin у/
zcosц/ -у sin ц/ xsini// -xcosi//
==Q2 ^-xcos21//-XSÍÜ2 i//)/ +(zcosi//sini//-vsin21
-(zcos2 y/- vsinЦ/COSI//)кJ = : Q2 (-x; zsin у/cosí// - vsin2 i//;-zcos2 у/ + vsin у/cosí//)
Итак,
F¡<6 = -/)йх(Пх/| = -pQrR3eM ^0;^-sin2i//;-cos2 у/ |
-PQ2(-
x; z sin i// eos i// - v sin21//; - z eos21// + y sin i// eos i//
)=
94
Rummim.......................................................
ефть и газ
!№ 4, 20.15
• ( 1
1
- pDrxi -pQr\ — R3eM sin2i//+ — zsin2i// - vsiir у/ \ j +
2
2
1
+pQr cos" ц/ +zcos" i//-— vsin2i//jA-.
Учет центробежной силы, вызванной вращением Земли, сводится к добавлению во второе уравнение полной системы уравнений Навье — Стокса центробежного ускоре-
f„ 1 .
-Qx (Qxr) = -QrR3eM I 0;—sin2i//;-cos2 ц/ | V 2 j
x; z sin ц/ cos ц/ -v sin2 i//;-zcos2 ц/ + v sin ц/cos ц/
У-
= Qrxi -Q2[ — R3eM sin2i// + — zsin2i// - vsin2 у/ 2 2
+ Q2| Д cos2 i// + zcos21//-— vsin2i// |A-
При этом полная система уравнений Навье — Стокса с учетом действия центробежной силы будет иметь вид
р,+Г-Ур+рсИуГ=0,
vt+(v-v)v+—Vp+-Vr=g-2QxF-Qx(Qxr)+
v ' ур у
AV .
P
(3)
Tl+V-VT+(y-l)TdvV = ^Xr+fJ°r(r ^
-ь(«г -ж)2 + (у„-м>:)-]+![("„ + (».-+^Т+{у--+м'у)~ :■
В качестве начальных условий при описании конкретного течения сжимаемого вязкого теплопроводного газа в случае постоянных значений коэффициентов вязкости и теплопроводности взяты функции, задающие точное решение [12] системы (1):
¿V.
и = 0, v = 0, w = 0, T0(z) = l-kz, к . I ■■= 0,0065 —, х00 = 50 7„, =288"К
p0(z) = (l-kz)v-1-, v=r-g --const > 0
(4)
(5)
Расчетная область представляет собой прямоугольный параллелепипед с длинами сторон х° = 1, V0 = 1 и =0,04 вдоль осей Ох, Оу и Ог соответственно. Для плотности на всех шести гранях параллелепипеда (х = 0, х = х° , у = 0, у = у0,
г = 0, 2 = ставится «условие непрерывности» потока [13]. Краевые условия для компонент вектора скорости газа соответствуют «условиям непротекания» для нормальной составляющей вектора скорости и «условиям симметрии» для двух других компонент вектора скорости течения [13]. Для температуры на всех шести гранях задаются условия теплоизоляции [13]. Продув газа через вертикальную трубу моделируется заданием вертикальной скорости течения газа в зависимости от времени t в виде
и
4, 2015
RI...........—""¿'..........ни.....................................................
есрть и газ
95
мЧ?) =0,003-[1-ехр (-10/)] (6)
через квадратное отверстие размером 0,1x0,1 в центре верхней грани расчетной области. Расчеты проводились при следующих входных параметрах: масштабные размерные значения плотности, скорости, расстояния и времени равны соответственно
Ао =1,2928-^-, г1т =333—, * =50 , =*• /и„„ =0.15с . м с
Разностные шаги по трем пространственным переменным Ах = Лу = 0,005 (размерное значение 0,25 м), А: = 0,004 (размерное значение 0,2 м), а шаг по времени N. =0,001 (размерное значение 0,00015 с).
Сопоставление 'результатов численного решения полной системы уравнений Навъе — Стокса без учета и с учетом действия центробежной силы.
Для сопоставления результатов численного решения полной системы уравнений На-вье — Стокса без учета влияния центробежной силы (1) и с учетом влияния центробежной силы (3) были проведены расчеты течений газа в восходящем закрученном потоке, инициированном вертикальным продувом газа вверх [15, 16] при одинаковых входных данных. Что касается визуального сопоставления результатов расчетов газодинамических параметров численного решения системы (1) и системы (3), то различия между ними не фиксируются. Сопоставление характерных безразмерных численных значений газодинамических характеристик возникающего течения на высоте 1 м в фиксированный момент времени, соответствующий 30000 расчетному шагу, приводится в таблице.
Сопоставление характерных безразмерных численных значений газодинамических характеристик
Газодинамический параметр Расчет без учета центробежной силы Расчет с учетом центробежной силы Отличие, %
Рmin 0,999559274 0,999559299 2,5-10"6
Т . min 0,998668086 0,998668132 4,6 ТО"6
Pmin 0,998227947 0,998228018 7,1 10"6
и max 0,0054516 0,0054337 0,328
и min -0,0052687 -0,0052865 0,338
V max 0,0053851 0,0053799 0,097
V min -0,0051108 -0,0051168 0,117
W max 0,0015431909 0.0015431826 5,4-10"4
W 3,1141 10"7 3,1142-10"7 3,2-10"3
3,0093.10"7 3,0083-10"7 3,3 10"2
Из сопоставления результатов расчета основных газодинамических характеристик без учета и с учетом влияния центробежной силы видно, что значения основных термодинамических параметров отличаются в среднем на 4,7-10~6 %. Центробежная сила практически не влияет на термодинамические параметры. Значения скоростных характеристик отличаются в среднем на 0,22 % для х-компоненты и у-компоненты и на 5,4-10~4% для г-компоненты. То есть влияние центробежной силы на вертикальную составляющую скорости течения частиц газа значительно меньше, чем на две другие компоненты. Отличие же значений энергетических характеристик не превосходит сотой доли процента. Таким образом, можно сделать вывод о том, что влияние центробежной силы на все газодинамические параметры сложных течений весьма незначительно.
Исследования поддержаны Министерством образования и науки РФ (проект № 3023).
9 б"
нефть и газ
№ 4, 2015
Список литературы
1. Баутин С. П. Представление решений системы уравнений Навье — Стокса в окрестности контактной характеристики // Прикладная математика и механика. - 1987. - Т. 51. - Вып. 4. - С. 574-584.
2. Баутин С. П. Характеристическая задача Коши и ее приложения в газовой динамике. - Новосибирск: Наука, 2009. - 368 с.
3. Баутин С. П, Обухов А. Г. Математическое моделирование разрушительных атмосферных вихрей. -Новосибирск: Наука, 2012. - 152 с.
4. Баутин С. П., Обухов А. Г. Математическое моделирование и численный расчет течений в придонной части тропического циклона // Вестник Тюменского государственного университета. Физико-математические науки. Информатика. -2012. -№ 4. - С. 175-183.
5. Обухов А. Г. Математическое моделирование и численные расчеты течений в придонной части торнадо // Вестник Тюменского государственного университета. Физико-математические науки. Информатика. - 2012. - № 4. -С. 183-189.
6. Баутин С. П., Обухов А. Г. Математическое моделирование придонной части восходящего закрученного потока // Теплофизика высоких температур. - 2013. - Т. 51. - № 4. - С. 567-570.
7. Баутин С. П. Торнадо и сила Кор полиса. - Новосибирск: Наука, 2008. - 96 с.
8. Вараксин А. Ю., Ромаш М. Э., Копейцев В. Н, Горбачев М. А. Моделирование свободных тепловых вихрей: генерация, устойчивость, управление // Теплофизика высоких температур. - 2010. - Т. 48. -№ 6. - С. 965-972.
9. Вараксин А. Ю., Ромаш М. Э., Копейцев В. Н, Горбачев М. А. Физическое моделирование воздушных смерчей: некоторые безразмерные параметры // Теплофизика высоких температур. - 2011. - Т. 49. -№ 2. - С. 317-320.
10. Вараксин А. Ю., Ромаш М. Э., Копейцев В. Н. Торнадо. - М.: Физматлит, 2011. - 312 с.
11. Баутин С. П., Баутин К. В., Макаров В. Н. Экспериментальное подтверждение возможности создания потока воздуха, закрученного силой Кориолиса //Вестник УрГУПС. - 2013. - № 2 (18). - С. 27-33.
12. Баутин С. П., Обухов А. Г. Одно точное стационарное решение системы уравнений газовой динамики // Известия вузов. Нефть и газ. - 2013. - № 4. - С. 81-86.
13. Баутин С. П., Обухов А. Г. Об одном виде краевых условий при расчете трехмерных нестационарных течений сжимаемого вязкого теплопроводного газа // Известия вузов. Нефть и газ. - 2013. - № 5. - С. 55-63.
14. Баутин С. П., Крутова И. Ю., Обухов А. Г., Баутин К. В. Разрушительные атмосферные вихри: теоремы, расчеты, эксперименты. - Новосибирск: Наука; Екатеринбург: Изд-во УрГУПС, 2013. - 215 с.
15. Абдубакова JI. В., Обухов А. Г. Численный расчет скоростных характеристик трехмерного восходящего закрученного потока газа //Известия высших учебных заведений. Нефть и газ. - 2014. - № 3. - С. 88-94.
16. Обухов А. Г., Абдубакова JI. В. Численный расчет термодинамических характеристик трехмерного восходящего закрученного потока газа // Вестник Тюменского государственного университета. Физико-математические науки. Информатика. - 2014. - № 7. - С. 157-165.
17. Абдубакова JI. В., Обухов А. Г. Численный расчет термодинамических параметров закрученного потока газа, инициированного холодным вертикальным продувом //Известия вузов. Нефть и газ. - 2014. - № 5 - С. 57-62.
18. Абдубакова JI. В., Обухов А. Г. Расчет плотности, температуры и давления трехмерного восходящего закрученного потока газа при вертикальном продуве // Нефтегазовое дело. - 2014. - Т. 12. - № 3. - С. 116-122.
Сведения об авторах
Баутин Сергей Петрович, д. физ.-мат. н., профессор кафедры «Прикладная математика», Уральский государственный университет путей сообщений, г. Екатеринбург, тел 89122221223, e-mail: SBau-tin@usurt. ru
Обухов Александр Геннадьевич, д. физ.-мат. н., профессор кафедры «Высшая математика», Тюменский государственный нефтегазовый университет, г. Тюмень, тел 89220014998, e-mail: [email protected]
Information about the authors Bautin S. P., Doctor of Physics and Mathematics, professor of the chair «Applied Mathematics», the Ural State University of Communication Lines, phone: 89122221223, e-mail: [email protected]
Obukhov A. G., Doctor of Physics and Mathematics, professor of the chair «Higher Mathematics», Tyumen State Oil and Gas University, phone: 89220014998, e-mail: ao-bukhov@tsogu. ru
