Способ минимизации схемной диффузии в численной модели аэродинамики Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2011 Математика и механика № 2(14)

УДК 519.63

М.А. Бубенчиков

СПОСОБ МИНИМИЗАЦИИ СХЕМНОЙ ДИФФУЗИИ В ЧИСЛЕННОЙ МОДЕЛИ АЭРОДИНАМИКИ

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

Ключевые слова: конвективно-диффузионная задача, естественные переменные, разности против потока, метод простой итерации, минимизация схемной или искусственной диффузии.

В работе [1] сформулированы семь положений, характеризующих схемную диффузию, появляющуюся при решении конвективно-диффузионных задач.

1. Схемная или искусственная диффузия имеет место, когда поток наклонен по отношению к линиям сетки и существует ненулевой градиент зависимой переменной в направлении по нормали к потоку.

2. Приближенное выражение для коэффициента искусственной диффузии в двухмерном случае дано в [2]:

D _ pV ArAz sin29

схем 4(r sin3 9 + Az cos3 9)

где V - модуль вектора скорости; 0 - угол наклона (от 0 до 90°) вектора скорости к направлению оси 0z. Из этого соотношения видно, что искусственная диффузия не появляется, если результирующий поток направлен вдоль одной из сеточных линий; кроме того, искусственная диффузия является максимальной, когда направление потока составляет угол 45° с линиями сетки.

3. Вклад искусственной диффузии можно уменьшить, используя меньшие шаги Az и А и располагая сетку (если это возможно) так, чтобы сеточные линии более или менее совпадали с направлением потока.

4. Поскольку реальная диффузия имеет место во многих задачах, то достаточно сделать искусственную диффузию малой по сравнению с реальной.

5. Использование центрально-разностной схемы не является средством избавления от искусственной диффузии. Как упоминалось ранее, центрально-разностная схема дает совершенно нереальные решения, если рассматриваются большие числа Пекле.

6. Основной причиной возникновения искусственной диффузии является практика обращения с потоком через каждую грань контрольного объема как с локально-одномерным.

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

Технология минимизации

Для минимизации схемной диффузии нами будет использован подход, аналогичный Рейсби [3]. Вспомним, что конвективная часть в стационарных уравнениях переноса (V• УФ) есть субстанциональная производная, которая является скоростью изменения переносимой субстанции во времени при перемещении жидкой или газообразной частицы вдоль ее траектории. Продолжая работу по изучению аэродинамики циклонной камеры [4] в рамках (у, 4, W, Т)-описания, в дальнейшем будем полагать, что Ф = W, 4, T. Здесь W - окружная компонента скорости; 4 - завихренность; T - температура. Таким образом

(-уф)_

- (дФ дФ^

В плоском случае V _(Vs ,Vn), УФ _ I----, I. Если в качестве координат (s, п)

V дs дп J

взять ортогональные координаты (z, r), то

/ - \ дФ дФ

(V •УФ)_^— + V—.

дz дт

Если же s - координата, отсчитываемая вдоль траектории, а n - по ее нормали, то

'5Ф дФ^ БФ /т- ,,дФ

дs

V _(V,0), УФ_(—, —| и — _(V•УФ) V-

V дs дп J Dt ' ’

Таким образом, в естественных координатах (л, п), связанных с траекторией жидкой частицы, запись конвективной части упрощается как в плоском, так и в пространственном случаях. Локальные значения естественных координат определяют ориентацию в пространстве конвективной ячейки, которая в общем случае не совпадает с разностной ячейкой диффузии.

Стало быть, подход Рейсби предполагает следующее представление конвективных членов в разностных уравнениях конвективного переноса

дФ дФА Ф, , -Ф* , ч

и 1Г+" »),,, = +0 (А5). (

Здесь V = + V2 ; Ф* - значение переносимой субстанции, найденное на гра-

нице разностной ячейки с наветренной стороны; Ал- - расстояние от центрального узла ячейки (г, г,) до найденной точки на границе (г*, г* ).

Как следует из представленного рисунка, Дл = -\]г* + г*2 . При этом начало отсчета системы координат (г*,г* ) берется в узле (г, г,). Угол наклона вектора ско-

( V А

рости к оси г определяется соотношением 0 = атС^ I ц ), а угол наклона диагона-

( Дг А

ли разностной ячейки к оси 0г есть у = аг^ I — I.

чДг )

Рис. 1. Диффузионная разностная ячейка (одинарная линия) и конвективная разностная ячейка (двойная линия)

Значение искомой величины Ф* в точке (г*, г* ) можно определить с использованием приведенной ниже таблицы.

Т аблица 1

Значение искомой величины с наветренной стороны

0 г* г* Ф,

[ -п ,-п+у] Аг Аг tg 0 Ф* _ Фі+:1,/ + А*Т(Ф‘+1"/+1 - Ф‘+1^./')

г П [-п+у, - -] Аг 0 Аг Ф*_ Ф і,]+1 +~А*г (Ф,+1"і+1 - Ф ‘"І+1)

г П [ - "2 ’ Аг ctg 0 Аг Ф*=фу+1 +^ К-, -ф/+1)

И,0] -Аг -Аг tg 0 ф*=фі-1, ,■+АТ (Ф,-и+1 -Ф,-и)

[0, у] -Аг -Аг tg 0 ф*=фі-1, ,+АТ (ф,-1,,-1 -фі-1, ,)

г П [у, 2] -Аг ^ 0 -Аг Ф*=фу-1+|г ((-и-1 -фу-1 )

[ ^ п-у] -Аг ^ 0 -Аг Ф*=Ф,1+^Аг (Ф-+1,,-1-Ф у-1)

[п-у, п ] Аг Аг tg 0 Ф*=Ф<+и +А- (Ф‘+и-1-Ф+1,,)

Тестовый пример

Рассмотрим следующую конвективно-диффузионную задачу:

и " + г “ = с

дг дг

2Ф 5 2Ф 1 5ФЛ

—_ +—Г + —д-

2 г дг

дг дг

Ф|г = Фг, где I = 1,4; {Гг} - прямоугольник.

(2)

(3)

Аппроксимируя конвективные члены в (2) по способу (1), а диффузионные члены с использованием центральных разностей, получим следующий разностный аналог рассматриваемого уравнения переноса:

ф -ф ^Ф,, -2Ф + Ф ,

1 ’ * 1 г+1,1 г,1 г-1,

Ау

І, І

Ф -2Ф + Ф ,

+ г, ,+1 г,] г, ,-1 +_

Аг 2

1 Фг,]+1 - 2Ф

., ]-1

Л

(4)

Аг 2

2Аг

Выражая из последнего соотношения Фг, найдем рекуррентную формулу для

перерасчета значений искомой сеточной функции во всех внутренних углах прямоугольной сетки:

Ф. , =

г, ]

V ІФ*

Ау

- + Б

1 Фг,]+1 - 2Ф

Фг+1,І +Фг-1,І , Фг,]+1 + Фг,]-1 + Аг 2 2Б 2Б

., ]-і

2Аг

Ау Аг2 Аг2

(5)

. .1 (І+1,]'+ 0

►

(Л,)

(г-1,,-1

Рис. 2. Фрагмент разностной сетки с направлением потока под углом 45° и ступенчатым изменением Ф

Как следует из второго пункта сформированных выше положений, схемная диффузия будет максимальной при 0 = 45о. Пусть Б = 0, Аг = Аг и поток будет направлен по диагонали к разностной ячейке. Покажем, что в этом случае ступенчатый характер распределения Ф не изменится при итерационном применении схе-

мы (5). Полагая в (5) Б = 0, получим

Тогда окончательно результатом применения (5) будет

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

На рис. 3 показаны расчетные распределения температуры Ф = Т, полученные по формуле (5) при Б = а = 0 (а - коэффициент температуропроводности среды) в однородном поле скорости, для которого и = 1, V = 1 (м/с). При этом на левой и верхней границах прямоугольника области определения температура была 350 К, а на правой и нижней - 350 К. Как видно из рисунка, профиль температуры остается ступенчатым, причем этот результат получается уже после первой глобальной итерации. Кроме того, если пролонгировать результат применения (5) к узлам верхней и правой границ прямоугольника, то получим полное распределение температуры без задания граничных условий на этих участках.

0,4

Рис. 3. Ступенчатый профиль температуры в однородном потоке, полученный с использованием рекуррентной формулы (5)

Заключение

В рамках (у, 4, W, Т)-описания динамики вязкой среды предложена простейшая технология минимизации численной диффузии, которая хорошо сочетается с методом простой итерации, применяемым для решения систем эллиптических уравнений.

ЛИТЕРАТУРА

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

2. De Vahl Davis G. andMallinson G.D. False Diffusion in Numerical Fluid Mechanics / Univ. of New South Wales, School of Mech. and Ind. Eng. Sydney, 1972.

3. Raithby G.D. Skew upstream differencing schemes for problems involving fluid flow // Computational Methods in Applied Mechanics and Engineering. 1976. V. 9. No. 2. P. 153-164.

4. БубенчиковМ.А., Иванова И.А. Расчет аэродинамики циклонной камеры // Вестник Томского государственного университета. Математика и механика. 2011. № 1 (13). С. 67-73.

Статья поступила 29.03.2011г.

Bubenchikov M. A. METHOD FOR MINIMIZING THE CIRCUIT DIFFUSION IN THE NUMERICAL MODEL OF AERODYNAMICS. Through linearization of the trajectory of a fluid particle on a fragment corresponding to a difference cell, and as a result of natural variables associated with the trajectory of particles, as well as the use of the differences against the flow, the convective terms were approximated with minimization of the circuit or artificial diffusion. The test calculations confirm this result. It is shown that the proposed computational technique is suitable for use in conjunction with the method of simple iteration.

Keywords: convection-diffusion problem, natural variables, differences against the flow, method of simple iteration, minimization of the circuit or artificial diffusion.

BUBENCHIKOV Mikhail Alekseevich (Tomsk State University)

E-mail: michael121@mail.ru