ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2017 Математика и механика № 49
УДК 532.5.032; 536.22 Б01 10.17223/19988621/49/11
А.В. Шваб, Н.И. Гичёва
МОДЕЛИРОВАНИЕ СВОБОДНОЙ И ВЫНУЖДЕННОЙ КОНВЕКЦИИ В ВИХРЕВОЙ КАМЕРЕ ХИМИЧЕСКОГО РЕАКТОРА
Исследуется влияние естественной конвекции на характер течения и теплообмен в вихревой камере химического реактора. С использованием двух подходов реализовано численное моделирование течения, получены поля линий тока, скоростей и температуры с учётом и без учёта свободной конвекции. Достоверность расчётов проверена сопоставлением результатов различных методов и сравнением с аналитическими результатами.
Ключевые слова: гидродинамика, теплообмен, вихревая камера, свободная конвекция, химический реактор.
Развитие промышленности требует повышения качеств материалов. Это обуславливает интерес науки к тугоплавким металлам, имеющим очень высокую износостойкость. Одним из самых востребованных в различных областях промышленности тугоплавких металлов является вольфрам. Однако получение вольфрама и придание изделиям из него нужной формы весьма затруднительно. Поэтому существует необходимость разработки новых способов выделения чистого вольфрама из его соединений. Перспективным является метод получения вольфрама путем осаждения из парогазовой фазы с помощью восстановления его фторидов или хлоридов. Процесс осаждения вольфрама из смеси гексафторида вольфрама и водорода протекает по реакции
WF6 + 3Н2 = W + 6ИБ, (1)
которая начинает идти при температурах порядка 300 °С и позволяет получать осадки любой толщины.
В данной работе приведен численный расчёт процессов гидродинамики и теп-лопереноса с учётом естественной конвекции в химическом реакторе, имеющем цилиндрическую форму (рис. 1, а). Впуск газов, имеющих температуру Т0, в реактор осуществляется через цилиндрическую трубу в центре со скоростью и0, а выпуск - через кольцевой канал. Верхняя и нижняя стенки сосуда приведены во вращательное движение с угловой скоростью ю. Нижняя стенка имеет температуру Т1.
Течение и теплоперенос газа в вихревой камере могут быть описаны системой уравнений, включающей уравнения Навье - Стокса, уравнение неразрывности, и уравнение переноса температуры. В безразмерной форме в цилиндрической системе координат для осесимметричного случая система имеет вид
2 ТИ,, ,, \ 1 ( Я2*, 1 А., Я2*, ,, \ и 2 ,, 2
диг + д(иг) + д(игиг) __др + _1_ д/ дг дг дг Яе
д иг 1 диг д и„ и„ - +
у дг2 г дг дг2
+
иф иг
ди2 д(игиг) + д(иг)2 __др + 1 д- я2- Л
д/ дг дг дг Яе
д иг 1 диг д иг
.„'2 V Йг* Д-2
У дг г дг дг )
диф + д(игиф) + д(и!иф) __ 2игиф + X д/ дг дг г Яе
2
д иф + 1 диф + дг 2 г дг
) г
игиг Ог
г + Яе2
д 2иф дг 2 иф]
dur du7 ur —- + —- + — = 0 ; dr dz r
de + d (ur 8) + (duz 8) = 1 dt dr dz Pr- Re
2 e +1 ae + a^e ^
dr2 r dr dz2
В качестве масштабов обезразмеривания для получения системы (2) взяты радиус вихревой камеры, входная линейная скорость, а также максимальный перепад температур. Характер гидродинамических и тепловых процессов в данной модели зависит от значений критериев Рейнольдса Яе = рЦК/ц, Прандтля Рг = срц/^ и Грасгофа вг = £р(Г1-Г0)рЯ3/ц2.
0.7
0.8
0.6
0.4
0.2
177777777X77777777777777^
0 0.3 0.8 1
Л 6. _А 5.
1, 4
0
Рис. 1. Вихревая камера химического реактора: расчётная область и граничные условия Fig. 1. Vortex chamber of the chemical reactor: computational domain and boundary conditions
0.2 0.4 0.6 0.8
Численное моделирование процессов в вихревой камере проводится двумя способами: в переменных «вихрь - функция тока» и в физических переменных «скорость - давление».
Для решения в переменных «вихрь - функция тока» система (2) тождественно преобразуется к следующему виду:
dQ . d(urП) d(uzП) 2u9 duф
dt
dr
dz
\ ^ f d2 n 1 ал а^п
dz Re I dr r dr dz
П
Gr ae
Re2 dr
dy = d2 y + a2 y 1 ay п dt dr2 dz2 r dr
(3)
dum d(uu) d(uu )
dt
™rU ф
dr
dz
+ J_ r Re
v dr2 r dr dz2 r2 y
ae + a (ur e) + (duz e) = 1
dt dr dz Pr- Re
a2 e +1 de + d2
dr2 r dr dz2
1
1
Рассматриваемая задача стационарна, однако её решение проводится эволюционным методом до установления по времени. Поэтому уравнения приводятся в нестационарном виде. Время в уравнениях фиктивно и, по сути, заменяет значение итерационного параметра. Для решения системы (3) с вышеприведёнными граничными условиями используется неявная обобщённая схема переменных направлений в Д-форме [1]. Она абсолютно устойчива и имеет второй порядок точности по времени. Данный двумерный метод основан на введении промежуточного временного слоя с целью построения многомерной неявной схемы, в которой требуется обращение только трёхдиагональной матрицы. Рассмотрим применение обобщённой схемы переменных направлений на примере уравнения для переноса тепла. Представим значение температуры на новом временном слое как суперпозицию старого значения и малого изменения Д9:
е"+1 _е" +де. (4)
Уравнение переноса для Д9 запишем в следующем виде:
де 1 д(иг де) 1 д(иг де) _ 1 д2 де 1 д2де _ + + 2 _
д/ 2 дг 2 дг 2Яе-Рг дг2 2Яе- Рг дг
Риге + _ 1
Яе-Рг
Гд2е" дV 1 де"Л —г+—т"+—^
У дг2 дг2 г дг
(5)
Проведём всё возмущение вначале в направлении г, а затем - в направлении г:
де* +1 д(игде*) 1 д2де* _
2
д/ 2 дг 2Яе-Рг дг
д(иге") д(иге") 1 Г д2е" д2е" 1 де"
V дг2 дг2 г дг
+
дг дг Яе-Рг,
** ** 2 ** *
** ** 2 е** *
д/ 2 дг 2Яе-Рг дг2 д/
(6)
де** 1 д(игде**) 1 1 д2де** де*
" + ----Г" . (7)
Таким образом, получаем
е"+1 _е"+де**. (8)
Полученная система линейных алгебраических уравнений решается методом прогонки.
Метод расщепления [2] применяется к уравнениям, записанным в физических переменных (2). В нём вместо уравнения неразрывности решается уравнение Пуассона для давления. Это делается для того, чтобы выделить влияние давления в одно уравнение. Уравнения для скорости представляются в следующем виде:
V "+1 _ V " "
-+ К _ Б __Ур" _У(8р). (9)
д/
Здесь К и Б - слагаемые, определяющие соответственно конвекцию и диффузию, 5р - разность между значениями давления на "-м и ("+1)-м временных слоях. Прибавим и отнимем в числителе промежуточное сеточное значение скорости V :
V "+1 _ V" + V _ V „
-+ К _ Б __Ур" _У(5р). (10)
д/
Разделим слагаемые в выражении (10) на два уравнения, содержащие соответст-
венно значения на п-м слое и промежуточные значения. Будем иметь
V — V п
+ С — Г = —Vрп ; (11)
Д/
и — V
= —^(8р). (12)
Vn+1 — V
Д/
Из уравнения (11) можно найти промежуточное сеточное значение скорости. Подействуем оператором градиент на выражение (12):
УУ п+1 — УУп 2
---= -У2 (8р). (13)
Д/
В силу выполнения уравнения неразрывности выражение сведётся к виду
2 УУ
У2(5р) = —. (14)
Д/
Получаем уравнение Пуассона для поправки к давлению. Метод его решения аналогичен используемому для уравнений переноса. Для поправки к давлению 5р при численном решении на всех границах задаются граничные условия второго рода (Неймана). Исходя из представленных рассуждений, выражения для значений скорости и давления на новом временном слое будут определяться уравнениями
Урп+1 =Урп +У(5р); (15)
УУ п+1 = V — Д/ -У(8р). (16)
Для аппроксимации конвективных и диффузионных слагаемых применяется экспоненциальная схема [2]. Решение реализуется на разнесённых сетках [3]. Для единственности решения задаются следующие граничные условия. На входе:
г 2
у(г) = О = 0, иг = 0, иг =—1, иф = 0, 0 = 0; на оси симметрии:
ди Я0
У = 0, О = 0, иг = 0, —^ = 0, иф = 0, — = 0 ;
дг дг
на нижней стенке камеры:
У = 0, о= 2(у — у-), иг = 0, и2 = 0, иф = Яю-г, 0 = 1, г Дп
где Яю = Я-ю/и - обратное число Россби. На правой стенке:
У = 0, О= 2(У —2 У -), иг = 0, иг = 0, иф= 0, д0 = 0;
г Дп дг
на верхней стенке камеры:
У = г2, О= 2(У— У -), иг = 0, иг = 0, иф= Яю - г, д0 = 0; 2 гДп2 г г ф дг
на внешней стенке входного и внутренней стенке выходного каналов:
у = О = 2(У^-), иг = 0, иг = 0, иф= 0, д0 = 0; 2 гДп дг
на выходе для всех функций:
0.
дг
Достоверность результатов была проверена тремя способами: сравнением результатов разных методов решений, сопоставлением с аналитическими результатами и исследованием на сеточную сходимость.
Ниже сопоставлены графики распределений характеристик потока в различных сечениях, полученные двумя методами решения (рис. 2). Цифры, обозначающие кривые, соответствуют сечениям, в которых построены распределения (см. рис. 1, б). Сплошной линией обозначены результаты расчётов в переменных «скорость - давление», точками - распределения полученные вычислениями в переменных «вихрь - функция тока». Рис. 2, а содержит сечения аксиальной составляющей скорости, 2, б - окружной составляющей, 2, в - температуры. Параметры потока, при которых проводилось данное сопоставление, следующие: Яе = 3, Рг = 1, Яю = 0.5, вг = 0. Рисунок 2 демонстрирует хорошее совпадение результатов, полученных различными методами.
Рис. 2. Сопоставление результатов расчётов, полученных разными методами Fig. 2. Comparison of the calculation results obtained by different methods
На выходе из реактора при малых числах Рейнольдса можно провести тестовое исследование на сравнение результатов с аналитической зависимостью для установившегося течения жидкости в кольцевом канале [4]. Течение описывается формулой
_ 2 (1 -г2)1п, _(1 _г12)1пг , (17)
г т ^ 2^ л '
(1 + г12)1п г1 + (1 _ г12)
где ит - среднее значение скорости. В данной задаче кольцевой канал имеет место на выходе из вихревой камеры (сечение 7 на рис. 1, б). Сравнение графиков распределения и2 дано на рис. 3, а. Сплошной линией обозначено распределение, полученное с использованием разработанной модели, точками - точные значения. Совпадение графиков также говорит о правильности расчётов: полученные результаты совпадают с аналитической зависимостью.
Ещё одним способом проверки достоверности модели является решение классической задачи о течении вблизи бесконечного вращающегося диска [5]. Рассматривается течение вблизи плоского диска, равномерно вращающегося с угловой скоростью вокруг оси, перпендикулярной к плоскости диска. Жидкость вдали от диска принимается покоящейся. Вследствие трения слой жидкости, непосредственно прилегающий к диску, увлекается последним и под действием центробежной силы отбрасывается наружу от диска. Взамен отброшенной жидкости к диску притекает в осевом направлении новая жидкость, которая также увлекается диском и опять отбрасывается наружу. С использованием ранее описанных методов были получены распределения скоростей вблизи диска. На рис. 3, б представлено сравнение численных результатов разработанной модели (сплошные линии) с аналитическим решением данной задачи, полученным В. Г. Кохрэном (обозначено точками). График демонстрирует хорошее совпадение решений, что также указывает на достоверность разработанной модели.
Рис. 3. Сравнение расчётов с аналитическими результатами Fig. 3. Comparison of the calculation results with analytical solutions
Далее представлены результаты численного решения. На рис. 4 изображены линии тока в вихревой камере при различных значениях критерия Рейнольдса при отсутствии (верхний ряд, вг = 0) и при наличии (нижний ряд, вг = 105) свободной конвекции. Параметры потока для всех графиков: Рг = 1, Яю = 0.5, Яе = 3, 20, 100
(соответственно по столбцам). Графики показывают, что свободная конвекция существенно влияет на гидродинамику закрученного течения, - изменяется положение циркуляционной зоны линий тока газа, усиливается вихревое движение.
z 0.8 0.6 0.4 0.20-
z
0.8 0.6 0.4
0 0.2 0.4 0.6 0.8 r
z
0.8
0.6 0.4 0.2 0
I
z
0.8 0.6 0.4 0.2 0
0 0.2 0.4 0.6 0.8 r
z
0.8
0.6 0.4 0.2 0
z
0.8 0.6 0.4 0.2 0
0 0.2 0.4 0.6 0.8 r
0 0.2 0.4 0.6 0.8 r
Рис. 4. Влияние числа Рейнольдса на распределение линий тока при отсутствии и при наличии свободной конвекции Fig. 4. Effect of the Reynolds number on the distribution of streamlines in the absence and presence of free convection
0
Рис. 5 показывает изолинии окружной составляющей скорости иф при различных значениях параметра закрутки и в зависимости от наличия естественной конвекции. Верхний ряд на рис. 5 иллюстрирует поля скорости, полученные вычислениями без учёта свободной конвекции (вг = 0), нижний ряд - с её учётом (вг = 105). Критериальные параметры потока, использованные при визуализации: Яе = 50, Рг = 1, Яю = 0.1, 0.5, 5. При учёте свободной конвекции поле окружной скорости значительно изменяется, область больших значений окружной составляющей скорости перемещается в направлении к оси симметрии.
На следующем рисунке (рис. 6) представлено температурное поле в камере реактора при различных величинах критерия Прандтля и в зависимости от наличия естественной конвекции. Верхний ряд графиков иллюстрирует распределение температуры в камере без учёта свободной конвекции (вг = 0), нижний - с её учётом (вг = 105). Параметры потока: Яе = 50, Яю = 0.5, Рг = 0.2, 1, Рг = 20. Рисунок 6 позволяет судить о том, что при наличии свободной конвекции температура в большей части вихревой камеры возрастает, что свидетельствует об интенсификации теплоотдачи на нижней стенке вихревой камеры.
Графики на рис. 7 являются изотермами в сечениях в вихревой камере, полученными без учёта свободной конвекции (кривые 1, 2, 3), и с её учётом (кривые 4, 5, 6). Номера кривых соответствуют номерам сечений, в которых они построены, приведённых на рис. 1, б.
0 0.2 0.4 0.6 0.8 r
0 0.2 0.4 0.6 0.8 Г
0 0.2 0.4 0.6 0.8 r
Рис. 5. Влияние числа Россби на распределение изолиний окружной составляющей скорости при отсутствии и при наличии свободной конвекции Fig. 5. Effect of the Rossby number on the distribution of peripheral velocity isolines in the absence and presence of free convection
1 0.8' 0.8
1 0.80.8
0 0.2 0.4 0.6 0.8 r
0 0.2 0.4 0.6 0.8 r
z 0.8
0.6
0.4
0.2
0
z 0.8
0.6
0.4
0.2
0
1
0.8 0.6 0.4 0.2 0
0 0.2 0.4 0.6 0.8 r
Рис. 6. Влияние числа Прандтля на распределение температуры
при отсутствии и при наличии свободной конвекции Fig. 6. Effect of the Prandtl number on the temperature distribution in the absence and presence of free convection
Изотермы иллюстрируют более быстрое и равномерное прогревание в большей части вихревой камеры при наличии естественной конвекции, а также свидетельствуют о необходимости её учёта при расчёте подобных химических реакторов.
z
0
0
0
0
0
Рис. 7. Влияние свободной конвекции на изотермы в указанных сечениях в вихревой камере Fig. 7. Effect of free convection on the isotherms in the sections indicated in the vortex chamber
Исследования, проведённые на основе разработанной модели течения в вихревой камере химического реактора, показывают необходимость учёта свободной конвекции в задачах расчёта аппаратов такого рода. Расчёты показали, что наличие естественной конвекции существенно влияет на гидродинамику потока и интенсифицирует теплообмен в вихревой камере.
ЛИТЕРАТУРА
1. Роуч П. Вычислительная гидромеханика. М.: Мир, 1977. 618 с.
2. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости: пер. с англ. под ред. Виленского. М.: Энергоатомиздат, 1984. 152 с.
3. Chung T.J. Computational fluid dynamics. Cambridge University Press, 2002. 1012 p.
4. Петухов Б.С. Теплообмен и сопротивление при ламинарном течении жидкости в трубах. М.: Энергия, 1967. 411 с.
5. ШлихтингГ. Теория пограничного слоя: пер. с нем. М.: Наука, 1974. 712 с.
Статья поступила 17.05.2017 г.
Shvab A.V., Gicheva N.I. (2017) MODELING OF THE NATURAL AND FORCED CONVECTION IN A VORTEX CHAMBER OF THE CHEMICAL REACTOR. Tomsk State University Journal of Mathematics and Mechanics. 49. pp. 114-123
DOI 10.17223/19988621/49/11
The subject of this research is a chemical reactor for producing tungsten. A physical and mathematical model of fluid flow and heat transfer in the chemical reactor under forced and free convection has been described and developed by two methods. The numerical simulation was carried out using "vortex-stream function" and "velocity-pressure" variables. The velocity and temperature distribution in reactor were obtained. The parametric studies of the effect of the Reynolds, Prandtl, and Rossby criteria on the flow characteristics were also performed. The graphs presented show that the natural convection has a significant impact on the flow hydrodynamics and intensifies the heat transfer.
The reliability of the calculations was verified by comparing the results obtained by the methods mentioned above. Furthermore, the created model was applied for numerical solving of a classical test problem of velocity distribution in an annular channel and that of the rotating infinite disk in a stationary liquid. The research results are in a good agreement with the analytical solutions.
Keywords: hydrodynamics, heat transfer, rotating flow, vortex chamber, natural convection, chemical reactor.
SHVAB Alexander Veniaminovich (Doctor of Physics and Mathematics, National Research Tomsk State University, Tomsk, Russian Federation) E-mail avshvab@inbox.ru
GICHEVA Natal'ya Ivanovna (National Research Tomsk State University, Tomsk, Russian Federation) E-mail tashka-zenit@mail.ru
REFERENCES
1. Roache P.J. (1976) Computational Fluid Dynamics. Hermosa: Albuquerque.
2. Patankar S. (1980) Numerical Heat Transfer and Fluid Flow. New York: Hemisphere.
3. Chung T.J. (2002) Computational Fluid Dynamics. Cambridge University Press.
4. Petukhov B.S. (1967) Teploobmen i soprotivlenie pri laminarnom techenii zhidkosti v trubakh [Heat transfer and resistance in a laminar flow in pipes]. Moscow: Energiya.
5. Schlichting H. (1951) Grenzschicht-Theorie [Boundary layer theory]. Karlsruhe: G. Braun.