Гибридная разностная схема - Википедия - Hybrid difference scheme

В гибридная разностная схема[1][2] - метод, используемый в численном решении для конвекция – диффузия проблемы. Впервые он был представлен Spalding (1970). Это комбинация центральная разностная схема и схема разности встречных волн поскольку он использует благоприятные свойства обеих этих схем.[3][4]

Вступление[5]

Гибридная разностная схема - это метод, используемый при численном решении задач конвекции-диффузии. Эти проблемы играют важную роль в вычислительная гидродинамика. Его можно описать общим частным уравнением следующим образом:[6]

(1)

Где, является плотность, - вектор скорости, это коэффициент диффузии и это исходный термин. В этом свойстве уравнения возможно температура, внутренняя энергия или компонент вектора скорости в направлениях x, y и z.

Для одномерного анализа задачи конвекции-диффузии в установившемся режиме и без источника уравнение сводится к

(2)

С граничными условиями, и , где L - длина, и - данные значения.

Генерация сетки

Интегрирующее уравнение 2 над контрольный объем содержащий узел N, и используя Теорема Гаусса т.е.

(3)

Дает следующий результат,

= (4)

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

= 0 (5)

Теперь давайте определим переменные F и D для представления конвекционный поток массы и диффузионная проводимость на гранях клеток,

и (6)

Следовательно, уравнения (4) и (5) преобразуются в следующие уравнения:

(7)
(8)

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

(9)

При низком числе Пекле (| Pe | <2) в потоке преобладает диффузия. При большом числе Пекле в потоке преобладает конвекция.

Схема разницы по центру и против ветра[3][7]

Рис. 1. Сетка, используемая для дискретизации в схеме центральных разностей.

В приведенных выше уравнениях (7) и (8), мы видим, что требуемые значения находятся на гранях, а не на узлах. Следовательно, для этого требуются приближения.

В схеме центральной разности мы заменяем значение на лицевой стороне средним значением в соседних узлах,

и (10)
Рис. 2: Сетка, используемая для дискретизации в схеме разницы против ветра для положительного числа Пекле (Pe> 0)
Рис. 3: Сетка, используемая для дискретизации в схеме разницы против ветра для отрицательного числа Пекле (Pe <0)

Помещая эти значения в уравнение (7) и переставляя, получаем следующий результат:

(11)

куда,

В схеме Upwind мы заменяем значение на лицевой стороне значением на соседнем восходящем узле. Например, для потока справа (Pe> 0), как показано на диаграмме, мы заменим значения следующим образом;

и (12)

А для Pe <0 ставим значения, как показано на рисунке 3,

и (13)

Помещая эти значения в уравнение (7) и переставляя, мы получаем то же уравнение, что и уравнение (11), со следующими значениями коэффициентов:

Гибридная разностная схема[3][7]

Рис. 4: Диаграмма, показывающая изменение любого свойства (ϕ) по длине (L) при разных числах Пекле (Pe)

Гибридная разностная схема Сполдинга (1970) представляет собой комбинацию центральной разностной схемы и разностной схемы против ветра. В нем используется центральная разностная схема второго порядка точности для малых чисел Пекле (| Pe | <2). Для больших чисел Пекле (| Pe |> 2) используется разностная схема Upwind, которая первого порядка точности, но учитывает конвекцию жидкости.

Как видно на рисунке 4, для Pe = 0 это линейное распределение, а для высоких Pe принимает значение выше по потоку в зависимости от направления потока. Например, значение на левой грани при различных обстоятельствах будет таким:

за (14)
за (15)
за (16)

Подставляя эти значения в уравнение (7) получаем такое же уравнение (11) со значениями коэффициентов:

Преимущества и недостатки

Он использует благоприятные свойства схемы центрального перепада и против ветра. Он переключается на разностную схему против ветра, когда центральная разностная схема дает неточные результаты для больших чисел Пекле. Он дает физически реалистичное решение и доказал свою полезность при прогнозировании практических потоков. Единственный недостаток, связанный с гибридной разностной схемой, заключается в том, что точность с точки зрения Серия Тейлор ошибка усечения это только первый заказ.

Смотрите также

Рекомендации

  1. ^ Патанкар, Сухас В. (1980). Числовой теплообмен и поток жидкости (14. полиграф. Ред.). Бристоль, Пенсильвания: Тейлор и Фрэнсис. ISBN  9780891165224.
  2. ^ Versteeg, H.K .; Малаласекера, В. (2007). Введение в вычислительную гидродинамику: метод конечных объемов (2-е изд.). Харлоу: Прентис Холл. ISBN  9780131274983.
  3. ^ а б c Скарборо, Дж. Б. (1958) Численный математический анализ, 4-е изд., Издательство Университета Джона Хопкинса, Балтимор, Мэриленд.
  4. ^ Сполдинг, Д. (1972). Новая конечно-разностная формулировка дифференциального выражения, включающего как первую, так и вторую производные, Int. J. Numer. Методы Eng., Vol. 4.
  5. ^ Поллард А. и Сиу А. Л. У. (1982). Расчет некоторых ламинарных течений с использованием различных схем дискретизации, Ж. вычисл. Методы Прил. Мех. Eng., Vol. 35.
  6. ^ Боррис, Дж. П., Брук, Д. Л. (1976). Решение уравнения неразрывности методом переноса с поправкой на поток, J. Comput. Phys., Vol. 16.
  7. ^ а б Роуч, П.Дж. (1976) Вычислительная гидродинамика, Эрмоса, Альбукерке, Нью-Мексико.

внешняя ссылка