ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Публикационная активность

(сведения по итогам 2017 г.)
2-летний импакт-фактор РИНЦ: 0,500
2-летний импакт-фактор РИНЦ без самоцитирования: 0,405
Двухлетний импакт-фактор РИНЦ с учетом цитирования из всех
источников: 0,817
5-летний импакт-фактор РИНЦ: 0,319
5-летний импакт-фактор РИНЦ без самоцитирования: 0,264
Суммарное число цитирований журнала в РИНЦ: 6012
Пятилетний индекс Херфиндаля по цитирующим журналам: 404
Индекс Херфиндаля по организациям авторов: 338
Десятилетний индекс Хирша: 17
Место в общем рейтинге SCIENCE INDEX за 2017 год: 527
Место в рейтинге SCIENCE INDEX за 2017 год по тематике "Автоматика. Вычислительная техника": 16

Больше данных по публикационной активности нашего журнале за 2008-2017 гг. на сайте РИНЦ

Вход


Забыли пароль? / Регистрация

Добавить в закладки

Следующий номер на сайте

4
Ожидается:
16 Декабря 2018

Разработка численной схемы стекловаренной печи

Development numerical scheme of glass furnace
Статья опубликована в выпуске журнала № 3 за 2010 год.[ 07.09.2010 ]
Аннотация:В статье приведен пример построения дискретного аналога трехмерной математической модели гидродинамики и теплообмена ванны стекловаренной печи. Выполнены необходимые преобразования уравнений неразрывности и балансов импульса, проведен переход к уравнениям дискретного аналога на разнесенной конечноразностной сетке. Также даны рекомендации по решению уравнений дискретного аналога с использованием итерационного метода и по нахождению граничных условий.
Abstract:The article gives an example of constructing a discrete analogue of three-dimensional mathematical model for bath hydrodynamics and heat exchange in glass furnace. Authors have done all the necessary transformations of the continuity equation and impulse balances equation, as well as transition to the equations of discrete analogue on straddling finite-difference grid. The article gives recommendations on the method of further solution to the equations of discrete analogue using an iterative method, and recommendations on finding the boundary conditions.
Авторы: Швыдкий В.С. (vshvit@isnet.ru) - Уральский государственный технический университет, г. Екатеринбург, , , доктор технических наук, Собянин С.Е. (ssergey@pochtamt.ru) - Уральский государственный технический университет, г. Екатеринбург, ,
Ключевые слова: аппроксимация, дискретный аналог, численная схема, стекловаренная печь, математическая модель
Keywords: approximation, discrete analog, numerical scheme, glass furnace, mathematical model
Количество просмотров: 6734
Версия для печати
Выпуск в формате PDF (5.84Мб)
Скачать обложку в формате PDF (1.43Мб)

Размер шрифта:       Шрифт:

При разработке численной схемы стекловаренной печи необходимо обязательно учесть особенности стекломассы: большую плотность и вязкость, медленное движение (практически ползущее). Учитывая их, можно уверенно отнести стекломассу к несжимаемой жидкости.

Конечно, плотность стекломассы зависит от температуры, а температура в бассейне изменяется, но, во-первых, эта зависимость слабая и, во-вторых, температура стекломассы в ванне изменяется от 1000 до 1580 оС, что приводит к относительно небольшим изменениям плотности. Иными словами, с достаточной точностью можно считать для стекломассы div v=0 и записать уравнения движения в виде:

·     в проекции на ось x

  (1)

·     в проекции на ось y

                (2)

·     в проекции на ось z

                         (3)

Вязкость стекломассы существенно зависит от температуры. Однако уравнения (1)–(3) – это уравнения баланса импульса для бесконечно малого объема. При составлении этого баланса можно считать, что температура стекломассы, следовательно, и ее вязкость, в пределах элементарного объема одинаковы. Тогда m можно вынести из-под знака производных, что приводит к упрощению уравнений.

Например, в уравнении (1) выделим слагаемые

то же сделаем и для уравнений (2) и (3). Заметим, что такой подход не исключает учета зависимости m от температуры, поэтому окончательный вид уравнений движения стекломассы можно представить следующим образом:

Уравнение неразрывности можно упростить. Хотя формально необходимо учитывать наличие внутренних источников (стоков) массы, фактически сделать это нельзя из-за неопределенности кинетики физико-химических превращений и отсутствия математического описания. Поэтому используем уравнение неразрывности в форме

                    (4)

Здесь возникло некоторое противоречие со сделанным утверждением о несжимаемости стекломассы, однако, во-первых, введение плотности под знаки производных не слишком усложняет уравнение, а, во-вторых, это важно для дальнейшего учета теплообмена.

Построение дискретного аналога всегда нужно выполнять для безразмерных уравнений, поскольку в этом случае легче оценить порядок аппроксимации и устойчивость численной схемы. Введем соответствующие безразмерные параметры. Предварительно, однако, приведем уравнения движения к дивергентному виду, объединив с уравнением неразрывности:

Отсчитаем компоненты скорости в долях среднерасходной скорости стекломассы в протоке V0=P/(86,4∙r0∙z2∙b2) м/с, где r0 – масштабное значение плотности стекломассы (кг/м3); P – производительность печи, т/сутки; z2 – высота прото- ка, м; b2 – ширина протока, м. Индекс 0 у теплофизических параметров характеризует значение при масштабной температуре T0=1000 оС. В качестве характерной длины выберем длину ванны L. Тогда имеем

Здесь h – глубина слоя стекломассы в ван- не, м; b – ширина ванны, м; x1 – длина загрузочной части ванны, м. Подстановка этих соотношений в консервативную форму записи уравнения движения в проекции на ось x приводит к выражению

или

                              (5)

Подпись:  
Рис. 1. Контрольный объем для аппроксимации 
уравнения неразрывности                            (6)

                    (7)

где P=p/(r0V02) – безразмерное давление, или число Эйлера; Re=r0V0L/m0 – число Рейнольдса; Fr=V02/(gL) – число Фруда.

Введем обозначения:

F1= G1=

H1= 

F2= G2= 

H2=

F3= G3= 

H3=

Тогда уравнения движения можно переписать в компактном виде:

                  (8)

Для аппроксимации уравнения неразрывности и уравнений системы (8) воспользуемся шахматной, или разнесенной, конечноразностной сеткой. Осуществим переход к разностным соотношениям на различных контрольных объемах (КО). Нужно отметить, что сетка на самом деле едина (в геометрическом плане) и индексация в различных КО изменяется лишь с целью более четкого выделения искомых параметров. В ряде случаев такая условность затрудняет построение дискретных аналогов (возникает неопределенность индексации), поэтому в данной работе следует связать индексацию с геометрическим центром КО, на котором аппроксимируется уравнение неразрывности. В этом центре определим все параметры состояния среды: давление, температуру и концентрацию.

Используя контрольный объем, соответствующий данным рисунка 1, запишем уравнение неразрывности в виде

Подпись:  
Рис. 2. Контрольный объем для аппроксимации 
уравнения импульса по оси x   (9)

Но по заданному условию теплофизические свойства берутся в центре КО и далее остаются неизменными во всем объеме параллелепипеда. Тогда во всех слагаемых плотность одинакова и, следовательно, исключается из (9). В результате имеем

                 (9а)

Можно видеть, что данная форма аппроксимации сохраняет физический смысл уравнения неразрывности как разности между массовыми расходами втекающей в КО и вытекающей из него стекломассы. По сравнению с первоначальной трактовкой изменена индексация компонентов скорости. Согласно физическому смыслу баланса массы определяемыми величинами являются компоненты скорости в геометрическом центре граней, то есть Uj-1/2,k,l, Uj+1/2,k,l, Vj,k-1/2,l, Vj,k+1/2,l, Wj,k,l-1/2 и Wj,k,l+1/2.

Перейдем к уравнениям баланса импульсов. В проекции на ось Х имеем

В новой индексации не все суммарные индексы попадают в целочисленные узлы. В соответствии с рисунком 2 координаты КО для оси Х изменяются в пределах: по Х от Xj до Xj+1; по Y от Yk-1/2 до Yk+1/2; по Z от Zl-1/2 до Zl+1/2. Применяя метод контрольного объема к уравнению х-компо­ненты импульса (8) с КО, показанным на рисун- ке 2, можно получить дискретное уравнение

        (10)

Центральным параметром данного КО является компонент Uj+1/2,k,l, поэтому все аппроксимации выполняются относительно этой величины. После небольших математических преобразований получим следующее уравнение, записанное в каноническом виде:

                        (11)

где  означает все конвективные и диффузионные вклады из соседних узлов.

Можно заметить, что уравнение (11) нелинейно и, следовательно, может быть решено только методом итераций. Теплофизические параметры на гранях КО можно определить либо как среднегеометрическое из соседних узловых значений, либо как среднеарифметическое.

Аналогично рассматриваются и другие координатные направления. КО для оси y лежит в диапазоне: по Х – от Xj -1/2 до Xj+1/2; по Y – от Yk до Yk+1; по Z – от Zl–½ до Zl+1/2 (рис. 3). При помощи КО, изображенного на рисунке 3, дискретная форма уравнения y-компоненты импульса (второе уравнение системы (8)) может быть записана в виде

         (12)

Здесь основной (искомой) переменной является Vj,k+1/2,l. После небольших математических преобразований получим уравнение для оси Y, записанное в канонической форме:

                         (13)

Наконец, в соответствии с КО, представленным на рисунке 4, для z-компоненты импульса (третье уравнение системы (8)) можно записать

        (14)

В данном случае определяющей величиной является компонента Wj,k,l+1/2, поэтому получаем следующее уравнение для оси Z, записанное в канонической форме:

                    (15)

Подведем итоги. Гидродинамическая часть численной модели на данном этапе включает четыре уравнения: одно уравнение неразрывнос- ти (9а):

Подпись:  
Рис. 3. Контрольный объем для аппроксимации 
уравнения импульса по оси y
 
Рис. 4. Контрольный объем уравнения импульсов 
по оси z 

и три уравнения баланса импульсов:

Уравнения баланса импульсов дадут верные значения компонентов скорости, превращающие уравнение неразрывности в тождество, лишь в том случае, когда при их решении будет использоваться правильное поле давления. Но даже тогда для решения системы понадобятся граничные условия. Для решения уравнений дискретного аналога можно воспользоваться методом SIMPLE С. Патанкара. В плане граничных условий основную сложность представляет задание компонентов скорости на свободной поверхности. Для неизотермической задачи здесь можно попытаться воспользоваться эффектом Марангони. Однако для изотермического потока этот выход отсутствует, и приходится изыскивать какие-то другие способы.

Литература

1.   Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. 152 с.

2.   Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 1978. 592 с.

3.   Дзюзер В.Я., Брылина Т.С., Собянин С.Е. Численное моделирование внешнего теплообмена в газоэлектрической стекловаренной печи с подковообразным направлением пламени // Стекло и керамика. 2008. № 5. С. 8–16.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=2574
Версия для печати
Выпуск в формате PDF (5.84Мб)
Скачать обложку в формате PDF (1.43Мб)
Статья опубликована в выпуске журнала № 3 за 2010 год.

Возможно, Вас заинтересуют следующие статьи схожих тематик: