Modeling of the natural and forced convection in a vortex chamber of the chemical reactor.pdf Развитие промышленности требует повышения качеств материалов. Это обуславливает интерес науки к тугоплавким металлам, имеющим очень высокую износостойкость. Одним из самых востребованных в различных областях промышленности тугоплавких металлов является вольфрам. Однако получение вольфрама и придание изделиям из него нужной формы весьма затруднительно. Поэтому существует необходимость разработки новых способов выделения чистого вольфрама из его соединений. Перспективным является метод получения вольфрама путем осаждения из парогазовой фазы с помощью восстановления его фторидов или хлоридов. Процесс осаждения вольфрама из смеси гексафторида вольфрама и водорода протекает по реакции WF6 + ЗН2 = W + 6HF, (1) которая начинает идти при температурах порядка 300 °С и позволяет получать осадки любой толщины. В данной работе приведен численный расчёт процессов гидродинамики и теп-лопереноса с учётом естественной конвекции в химическом реакторе, имеющем цилиндрическую форму (рис. 1, а). Впуск газов, имеющих температуру Т0, в реактор осуществляется через цилиндрическую трубу в центре со скоростью U0, а выпуск - через кольцевой канал. Верхняя и нижняя стенки сосуда приведены во вращательное движение с угловой скоростью со. Нижняя стенка имеет температуру 7\. Течение и теплоперенос газа в вихревой камере могут быть описаны системой уравнений, включающей уравнения Навье - Стокса, уравнение неразрывности, и уравнение переноса температуры. В безразмерной форме в цилиндрической системе координат для осесимметричного случая система имеет вид 2 Я/"-, ,, \ 1 (1 А., „Л и 2 ,, 2 д и„ 1 диг д и„ и„ - + МФ и„ + v дг2 г дг dz2 д и 1 ди, д и, -f +--- +-;f р, Pip, д 2 v дг г дг dz j J г Uruz Gr г Re2 d\„ и А г duф | д(игиф) ^ 5(мгмф) _ 2мгмф ^ \ dt дг dz г Re v дг2 г дг dz2 г2 у dur | д(иг) | d(uzur) = dp 1 dt дг dz dr Re ди2 | d(uruz) | д(и2)2 _ dp | 1 1 Л dt дг dz dz Re dur du, и,. -- +-- + - = 0 ; dr dz г 5V 5z2 мге r 1 59 +---h 528 58 | 5(иг8) | (ди,6) _ 1 dt дг dz Рг- Re I dr2 r dr В качестве масштабов обезразмеривания для получения системы (2) взяты радиус вихревой камеры, входная линейная скорость, а также максимальный перепад температур. Характер гидродинамических и тепловых процессов в данной модели зависит от значений критериев Рейнольдса Re - р(//?/ц, Прандтля Рг = Срц/Х и Грасгофа Gr = gP(ri-r0)pR7^2. Тп, и() 1 0.8 i 0.6 0.4 0.2 6 7 3, 6 2, 5 L 4 0.2 0.4 0.6 0.8 0 0.7 0 Рис. 1. Вихревая камера химического реактора: расчётная область и граничные условия Fig. 1. Vortex chamber of the chemical reactor: computational domain and boundary conditions Численное моделирование процессов в вихревой камере проводится двумя способами: в переменных «вихрь - функция тока» и в физических переменных «скорость - давление». Для решения в переменных «вихрь - функция тока» система (2) тождественно преобразуется к следующему виду: 5Q d(urQ.) S(k,Q) 2мф 5иф \ ( d2Q 1 5Q 52Q Q Gr 59 - +--+ r dr 2 Re2 dr ' dz2 dt dr dz r dz Rel dr 5ш д w д ш 1 5ш ^ - =-f + -f-----rQ dt dr dz r dr (3) 5иф | d(uruv) ^ d{uzuv) dt dr dz 58 | d(ur8) | (duz9) dt dr dz 2uruф f 1 f Л 2 5 u. 2 r dr dz2 r~ r Re dz' u„ 9 dr 528 1 58 529^ - + - Pr- Re I dr2 r dr dz2 1 Рассматриваемая задача стационарна, однако её решение проводится эволюционным методом до установления по времени. Поэтому уравнения приводятся в нестационарном виде. Время в уравнениях фиктивно и, по сути, заменяет значение итерационного параметра. Для решения системы (3) с вышеприведёнными граничными условиями используется неявная обобщённая схема переменных направлений в Д-форме [1]. Она абсолютно устойчива и имеет второй порядок точности по времени. Данный двумерный метод основан на введении промежуточного временного слоя с целью построения многомерной неявной схемы, в которой требуется обращение только трёхдиагональной матрицы. Рассмотрим применение обобщённой схемы переменных направлений на примере уравнения для переноса тепла. Представим значение температуры на новом временном слое как суперпозицию старого значения и малого изменения ДО: еи+1=е"+де. (4) Уравнение переноса для Д9 запишем в следующем виде: Д9 1 d(ur А9) 1 d(uzAQ) 1 52А9 1 52А9 _ At 2 дг 2 dz 2 Re-Pr dr2 2 Re-Pr dz2 fd2Qn d2S" р мг9+_ 1 (5) r Re-Pr - + - v dr dz2 r dr Проведём всё возмущение вначале в направлении г, а затем - в направлении z: де* 1 д(игде*) 1 з2де* _ 2 At 2 дг 2Re-Pr дг 5(мге") 5(мге") 1 (s2e" s2e" 1 да" (6) дг dz Re - Pr - + - V дг2 dz2 г дг де** 15(мгде**) 1 1 д2де** де* -+--------^ = -. (7) At 2 dz 2 Re- Pr dz2 At Таким образом, получаем en+1 =е"+де**. (8) Полученная система линейных алгебраических уравнений решается методом прогонки. Метод расщепления [2] применяется к уравнениям, записанным в физических переменных (2). В нём вместо уравнения неразрывности решается уравнение Пуассона для давления. Это делается для того, чтобы выделить влияние давления в одно уравнение. Уравнения для скорости представляются в следующем виде: v"+1-v" -+ К -D =-Vp"-Щ8р). (9) At Здесь KviD- слагаемые, определяющие соответственно конвекцию и диффузию, 8р - разность между значениями давления на п-ш и (п+1)-м временных слоях. Прибавим и отнимем в числителе промежуточное сеточное значение скорости v : - + K-D = -Vp"-V(8p). (10) At Разделим слагаемые в выражении (10) на два уравнения, содержащие соответственно значения на п-м слое и промежуточные значения. Будем иметь v -v At = -V (dp) . = "V(Sp). At Из уравнения (11) можно найти промежуточное сеточное значение скорости. Подействуем оператором градиент на выражение (12): Vv"+1 - W" At В силу выполнения уравнения неразрывности выражение сведётся к виду 9 W V2(§р) = -. At Получаем уравнение Пуассона для поправки к давлению. Метод его решения аналогичен используемому для уравнений переноса. Для поправки к давлению 5р при численном решении на всех границах задаются граничные условия второго рода (Неймана). Исходя из представленных рассуждений, вьфажения для значений скорости и давления на новом временном слое будут определяться уравнениями Vpn+l =Vpn +V(8p); (15) Vv"+1 = - At- У (dp) . (16) . + K-D = -Vp" (П) (12) (13) (14) Для аппроксимации конвективных и диффузионных слагаемых применяется экспоненциальная схема [2]. Решение реализуется на разнесённых сетках [3]. Для единственности решения задаются следующие граничные условия. На входе: 2 Q = 0, и =0, и =-1 \\/(г) = - «„=0, 1 = 0 на оси симметрии: ды ов У = о, Q = О, мг = 0, - = 0, и =0, - = 0 : дг дг на нижнеи стенке камеры: ur= 0, uz= 0, u =Ra>-r, 'w+l гМ2 ц/ = 0, Q = - ' = 1, где Rco - R-co/U- обратное число Россби. На правой стенке: 0 Q=2(yw+1-yJ 8в г An на верхней стенке камеры: ur= 0, uz= 0, u =Rco-r, 2 59 ~dz Q = = 0; Ц1=г An2 на внешней стенке входного и внутренней стенке выходного каналов: 59 5z 'w+l rM2 Q = ur =0, uz = 0, мф = 0, = 0 : Ц1=- 2 на выходе для всех функций: dz Достоверность результатов была проверена тремя способами: сравнением результатов разных методов решений, сопоставлением с аналитическими результатами и исследованием на сеточную сходимость. Ниже сопоставлены графики распределений характеристик потока в различных сечениях, полученные двумя методами решения (рис. 2). Цифры, обозначающие кривые, соответствуют сечениям, в которых построены распределения (см. рис. 1,6). Сплошной линией обозначены результаты расчётов в переменных «скорость - давление», точками - распределения полученные вычислениями в переменных «вихрь - функция тока». Рис. 2, а содержит сечения аксиальной составляющей скорости, 2, б - окружной составляющей, 2, е - температуры. Параметры потока, при которых проводилось данное сопоставление, следующие: Re = 3, Pr= 1, Rm - 0.5. Gr - 0. Рисунок 2 демонстрирует хорошее совпадение результатов, полученных различными методами. Рис. 2. Сопоставление результатов расчётов, полученных разными методами Fig. 2. Comparison of the calculation results obtained by different methods На выходе из реактора при малых числах Рейнольдса можно провести тестовое исследование на сравнение результатов с аналитической зависимостью для установившегося течения жидкости в кольцевом канале [4]. Течение описывается формулой 2 т 2 ч 1 /Л 2 Ч ' (1 + rj )lnrj +(l-rj ) где Um - среднее значение скорости. В данной задаче кольцевой канал имеет место на выходе из вихревой камеры (сечение 7 на рис. 1, б). Сравнение графиков распределения иг дано на рис. 3, а. Сплошной линией обозначено распределение, полученное с использованием разработанной модели, точками - точные значения. Совпадение графиков также говорит о правильности расчётов: полученные результаты совпадают с аналитической зависимостью. Ещё одним способом проверки достоверности модели является решение классической задачи о течении вблизи бесконечного вращающегося диска [5]. Рассматривается течение вблизи плоского диска, равномерно вращающегося с угловой скоростью вокруг оси, перпендикулярной к плоскости диска. Жидкость вдали от диска принимается покоящейся. Вследствие трения слой жидкости, непосредственно прилегающий к диску, увлекается последним и под действием центробежной силы отбрасывается наружу от диска. Взамен отброшенной жидкости к диску притекает в осевом направлении новая жидкость, которая также увлекается диском и опять отбрасывается наружу. С использованием ранее описанных методов были получены распределения скоростей вблизи диска. На рис. 3, б представлено сравнение численных результатов разработанной модели (сплошные линии) с аналитическим решением данной задачи, полученным В. Г. Кохрэном (обозначено точками). График демонстрирует хорошее совпадение решений, что также указывает на достоверность разработанной модели. Рис. 3. Сравнение расчётов с аналитическими результатами Fig. 3. Comparison of the calculation results with analytical solutions Далее представлены результаты численного решения. На рис. 4 изображены линии тока в вихревой камере при различных значениях критерия Рейнольдса при отсутствии (верхний ряд, Gr = 0) и при наличии (нижний ряд, Gr - 105) свободной конвекции. Параметры потока для всех графиков: Рг - 1, Rco - 0.5, Re - 3, 20, 100 (соответственно по столбцам). Графики показывают, что свободная конвекция существенно влияет на гидродинамику закрученного течения, - изменяется положение циркуляционной зоны линий тока газа, усиливается вихревое движение. Z 0.8 0.6 0.4 0.20- Z 0.80.60.4- 0 0.2 0.4 0.6 0.8 г 0 0.2 0.4 0.6 0.8 г Рис. 4. Влияние числа Рейнольдса на распределение линий тока при отсутствии и при наличии свободной конвекции Fig. 4. Effect of the Reynolds number on the distribution of streamlines in the absence and presence of free convection О 0.2 0.4 0.6 0.8 Г Z 0.80.6 0.40.20- Z 0.80.60.40.20- О 0.2 0.4 0.6 0.8 Г о Рис. 5 показывает изолинии окружной составляющей скорости иф при различных значениях параметра закрутки и в зависимости от наличия естественной конвекции. Верхний ряд на рис. 5 иллюстрирует поля скорости, полученные вычислениями без учёта свободной конвекции (Gr - 0), нижний ряд - с её учётом (Gr - 105). Критериальные параметры потока, использованные при визуализации: Re = 50, Рг - 1, Rco = 0.1, 0.5, 5. При учёте свободной конвекции поле окружной скорости значительно изменяется, область больших значений окружной составляющей скорости перемещается в направлении к оси симметрии. На следующем рисунке (рис. 6) представлено температурное поле в камере реактора при различных величинах критерия Прандтля и в зависимости от наличия естественной конвекции. Верхний ряд графиков иллюстрирует распределение температуры в камере без учёта свободной конвекции (Gr = 0), нижний - с её учётом (Gr = 105). Параметры потока: Re = 50, Rco = 0.5, Pr = 0.2,1, Pr = 20. Рисунок 6 позволяет судить о том, что при наличии свободной конвекции температура в большей части вихревой камеры возрастает, что свидетельствует об интенсификации теплоотдачи на нижней стенке вихревой камеры. Графики на рис. 7 являются изотермами в сечениях в вихревой камере, полученными без учёта свободной конвекции (кривые 1, 2, 3), и с её учётом (кривые 4, 5, б). Номера кривых соответствуют номерам сечений, в которых они построены, приведённых на рис. 1, б. О 0.2 0.4 0.6 0.S Г О 0.2 0.4 0.6 0.S Г О 0.2 0.4 0.6 0.S Г Рис. 5. Влияние числа Россби на распределение изолиний окружной составляющей скорости при отсутствии и при наличии свободной конвекции Fig. 5. Effect of the Rossby number on the distribution of peripheral velocity isolines in the absence and presence of free convection э г l O.S 0.6 0.4 0.2 0 0 0.2 0.4 0.6 O.S Г 0 0.2 0.4 0.6 O.S Г 0 0.2 0.4 0.6 O.S Г Рис. 6. Влияние числа Прандтля на распределение температуры при отсутствии и при наличии свободной конвекции Fig. 6. Jsffect of the Prandtl number on the temperature distribution in the absence and presence of free convection Изотермы иллюстрируют более быстрое и равномерное прогревание в большей части вихревой камеры при наличии естественной конвекции, а также свидетельствуют о необходимости её учёта при расчёте подобных химических реакторов. Рис. 7. Влияние свободной конвекции на изотермы в указанных сечениях в вихревой камере Fig. 7. Effect of free convection on the isotherms in the sections indicated in the vortex chamber Исследования, проведённые на основе разработанной модели течения в вихревой камере химического реактора, показывают необходимость учёта свободной конвекции в задачах расчёта аппаратов такого рода. Расчёты показали, что наличие естественной конвекции существенно влияет на гидродинамику потока и интенсифицирует теплообмен в вихревой камере.

Скачать электронную версию публикации