Метод контрольного объема

Метод контрольного объема можно рассматривать как частный случай метода взвешенных невязок. Основная идея метода контрольного объема легко понятна и поддается прямой физической интерпретации. Расчетную область разбивают на некоторое число непересекающихся контрольных объемов таким образом, что каждая узловая точка содержится в одном контрольном объеме. Дифференциальное уравнение интегрируют по каждому контрольному объему. Для вычисления интегралов используют кусочные профили, которые описывают изменение Ф между узловыми точками. В результате находят дискретный аналог дифференциального уравнения, в который входят значения Ф в нескольких узловых точках. Полученный подобным образом дискретный аналог выражает закон сохранения Ф для конечного контрольного объема точно так же, как дифференциальное уравнение выражает закон сохранения для бесконечно малого контрольного объема. Одним из важных свойств метода контрольного объема является то, что в нем заложено точное интегральное сохранение таких величин, как масса, количество движения и энергия на любой группе контрольных объемов и, следовательно, на всей расчетной области. Это свойство проявляется при любом числе узловых точек, а не только в предельном случае очень большого их числа. Таким образом, даже решение на грубой сетке удовлетворяет точным интегральным балансам.

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

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

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

Хотя при учете конвекции добавляется только один новый член, вводимый в настоящей главе, его аппроксимация оказывается достаточно сложной. Конвективный член непосредственно связан с диффузионным, и их целесообразно рассматривать как целое. Именно поэтому настоящая глава названа «Конвекция и диффузия». Необходимо помнить, что слово «диффузия» используется здесь в обобщенном значении. Оно включает не только диффузию химических компонент, вызванную градиентами концентрации. Диффузионный поток, вызванный градиентом обобщенной переменной Ф, определяется как — Г д Ф /дхj для конкретной величины Ф он может представлять собой диффузионный поток химических компонент, тепловой поток, вязкое напряжение и т. д. Обобщенное дифференциальное уравнение содержит член (д/дхj)( Г д Ф /дхj), который определяется как диффузионный член. Фактически это выражение представляет собой сумму трех составляющих по трем координатным направлениям, однако удобно рассматривать их совместно как единый диффузионный член. Аналогичные рассуждения справедливы для конвективного члена (д/дхj)(риj Ф ). Отметим одно из свойств конвективно-диффузионной задачи. Поскольку заданное поле течения должно удовлетворять уравнению неразрывности

(1.21)

 

то общее дифференциальное уравнение

(1.22)

можно записать в виде:

(1.23)

 

На данном этапе имеются все составляющие, необходимые для получения дискретного аналога, соответствующего общему дифференциальному уравнению (1.22).

Сначала найдем аппроксимацию двухмерного уравнения, но ту же процедуру применим и к трехмерному случаю.

Рассмотрим контрольный объем на рис. 3. Используя опыт, приобретенный при анализе одномерной задачи для получения суммарного теплового потока Jе, и предположив, что найденное выражение применимо ко всей грани контрольного объема площадью Δу×1,сможем сразу записать дискретный аналог для двухмерной задачи.

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

Рис 3. Контрольный объем (заштрихованная область) для двухмерного случая

Уравнение (3.2) в двухмерной форме можно представить в виде

(3.1.1)

где Jx и Jу — суммарные (конвекция плюс диффузия) потоки, определенныеследующим образом:

(3.1.2)

где и, v — компоненты скорости в направлениях осей х и у. Интегрирование уравнения (3.1.1) по контрольному объему, показанному на рис. 3.1.1, дает

(3.1.3)

 

где источниковый член линеаризован стандартным образом, а для нестационарного члена ρP и ФР полагаются преобладающими по всему контрольному объему. Старые значения (на предыдущем шаге по времени) обозначаются рр° и ФР°.

В соответствии с полностью неявным способом аппроксимации все другие величины (не имеющие верхних индексов) рассматриваются как новые значения. Величины Jе, Jw, Jn и Js представляют собой проинтегрированные по граням контрольного объема суммарные потоки, т. е. Jе= ∫ Jdyx, для грани е и т. д.

Аналогичным образом можно проинтегрировать уравнение неразрывности (3.1) по контрольному объему

(3.1.4)

где Fе, Fw, Fп и Fs — массовые расходы потока через грани контрольного объема.

Если ρи в точке е считается преобладающей по всей грани е контрольного объема, то

(3.1.5)

Аналогично

(3.1.6)

Если умножить уравнение (3.1.4) на ФР и вычесть его из уравнения (3.1.3), получим

(3.1.7)

Этот способ с использованием (3.1.3) и (3.1.4) позволяет получить (3.1.7), являющееся дискретным аналогом комбинации уравнений (3.1) и (3.2), на основе которых получено (3.3).

Другим способом можно получить дискретный аналог из уравнения (3.3), но этот способ не так удобен.

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

(3.1.8)

где

(3.1.9)

 

Здесь D и Dw подобно Fе и Fw, содержат площадь Δy×1 граней е и w.

Используя выражения для Jn - FnФР и J s- FsФP, можно записать окончательный вид дискретного аналога. Уравнения (3.1.8) свидетельствуют, что правило относительно суммы соседних коэффициентов легко выполняется.

Когда заданные поля скорости и плотности удовлетворяют дискретному аналогу уравнения неразрывности, приведенные выше вывод и вывод, основанный на (3.1.4), будут давать идентичные дискретные аналоги. Однако когда заданные поля не удовлетворяют уравнению неразрывности, эти два вывода дают разные уравнения и приводят к разным решениям.

Окончательный вид дискретного аналога. Двухмерный дискретный аналог можно записать в следующем виде:

(3.1.10)

где

(3.1.11)

Здесь и обозначают известные значения для времени t, а все другие величины ( , , , , и т. д.) представляют собой неизвестные величины для времени t+Δt.

Массовые расходы Fе, Fw, Fп и Fs определены уравнениями (3.1.5)-(3.1.6). Соответствующие проводимости представим в виде

(3.1.12)

 

А числа Пекле

(3.1.13)

Уравнения количества движения

Еще раз напомним, что при заданном поле давления решение уравнений количества движения можно получить с помощью дискретного аналога уравнения для обобщенной переменной Ф. Для уравнения количества движения Ф обозначает одну из составляющих скорости и коэффициенту Г и свободному члену S следует придать соответствующий смысл. При использовании шахматной сетки дискретные аналог и уравнений количества движения несколько отличаются от дискретных аналогов уравнений для других Ф, рассчитываемых в узлах основной сетки. Однако это отличие относится к несущественным деталям. Оно связано с использованием для аппроксимации уравнений количества движения контрольных объемов на шахматной сетке.

Рис. 3.2.1 Контрольный объем для и.

Рис. 3.2.2 Контрольный объем для v.

Контрольный объем для уравнения количества движения в направлении оси х показан на рис. 3.2.1. Если иметь в виду точки для нахождения только составляющей u, в этом контрольном объеме нет ничего необычного. Его грани лежат между точкой е и соответствующими соседними точками для u. Однако он смещен по отношению к обычному контрольному объему, расположенному вокруг узловой точки P. Смещение объема произошло только в направлении оси х таким образом, что перпендикулярные этому направлению грани проходят через основные узловые точки Р и Е. Отсюда видно одно из главных достоинств шахматной сетки: разность р P- рЕ можно использовать для расчета силы давления, действующей на контрольный объем для скорости и.

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

(3.2.1)

Здесь число соседних членов зависит от размерности задачи. Для двухмерной задачи на рис.3.2.1 показаны четыре точки вне контрольного объема, значения скорости в которых входят в соседние члены (3.2.1); в трехмерном случае войдут шесть соседних значений и. Значения коэффициентов аnb связаны с влиянием совместных конвективных и диффузионных процессов на гранях контрольного объема. Член b определяется так же, как в (3.1.11),но градиент давления не включен в составляющие источникового члена SС или SР. С градиентом давления связан последний член в (3.2.1). Так как требуется определить поле давления, было бы нецелесообразно включать давление в источниковый член уравнения количества движения. Член (р W- рЕ) Ае представляет собой силу давления, действующую на контрольный объем для и, а Ае — площадь поверхности, на которую действует этот перепад давления. В двухмерном случае Ае =Δy×1, в трехмерном Ае =Δy× Δz

Уравнения количества движения в других направлениях аппроксимируются таким же образом. На рис. 3.2.2 показан контрольный объем для уравнения количества движения в направлении оси у; он смещен вдоль оси у. Дискретный аналог будет иметь вид

(3.2.2)

где (р P- рN) Ап — соответствующая сила давления. В трехмерном случае аналогичное уравнение можно записать для составляющей скорости w.

Уравнения количества движения можно решить только в том случае, если поле давления задано или каким-то образом найдено. Если при решении использовать неверное поле давления, найденное поле скорости не будет удовлетворять уравнению неразрывности.

Выразим такое поле скорости, полученное с использованием приближенного поля давления р*, через u *, v*, w*. Это поле скорости находится в результате решения следующих уравнений:

(3.2.3)

В этих уравнениях составляющим скорости и давлению приписан верхний индекс *.

Отметим, что точка t лежит на сеточной линии, направленной вдоль оси z и проходящей через узловые точки Р и Т.

Поправки скорости и давления

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

p = p * + p ' (3.3.1)

где р' назовем поправкой давления. Надо выяснить, как будут изменяться составляющиескорости в соответствии с таким изменением давления. Аналогичным образом введемсоответствующие поправки составляющих скорости:

u = u * + u '; v = v * + v '; w = w * + w '(3.3.2)

Вычитая (3.2.3) верхнее уравнение из (3.2.2), получаем

au ' =Σ au ' + (p ' − p ') A. (3.3.3)

Теперь выбросим член Σ аnbu’nb из уравнения.

В результате получим

a u ' = (p ' − p ') A (3.3.4)

или

u = d p − p (3.3.5)

где

(3.3.6)

Уравнение (3.3.5) назовем поправочной формулой для скорости. Его можно переписать в виде

(3.3.7)

Отсюда видно, какой должна быть поправка к значению скорости ие*, определяемаяпоправками давления, чтобы получилось значение ие.

Аналогичным образом запишем поправочные формулы для других составляющихскорости:

(3.3.8)

 

Уравнение для поправки давления

Преобразуем уравнение неразрывности в уравнение для поправки давления.

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

(3.4.1)

Проинтегрируем это уравнение по показанному на рис. 3.4.1 заштрихованному контрольному объему (на рисунке показано только двухмерное сечение). Напомним, что такой контрольный объем использовался.при выводе дискретного аналога уравнения для обобщенной переменной Ф.

При интегрировании члена дρ/дt предположим, что значение плотности во всем контрольном объеме равно ρP. Также будем считать, что значение массовой скорости на всей грани контрольного объема определяется значением составляющей скорости ие в точке е. В соответствии с полностью неявной аппроксимацией предположим, что новые значения скорости и плотности (в момент времени t +Δ t) преобладают на всем шаге по времени; старое значение плотности ρP0 (в момент времени t) будет входить только из-за наличия члена дρ!дt.

В этих предположениях интегрирование уравнения (3.4.1) дает

(3.4.2)

Рис. 3.4.1 Контрольный объем для уравнения неразрывности (заштрихованная область — контрольный объем)

Если теперь вместо всех составляющих скорости подставить их выражения из поправочных формул для скорости [таких, как (3.3.7)-(3.3.8)], то после группировки соответствующих членов получим следующее уравнение для сеточных значений р':

(3.4.3)

Где

(3.4.4)

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

Из (3.4.4) можно видеть, что член b уравнения для поправки давления по существу равен (со знаком минус) левой части дискретного аналога уравнения неразрывности (3.4.2),записанного через значения составляющих скорости с индексом *. Равенство b =0 означает, что эти составляющие вместе с имеющимся значением (ρP0 - ρP) удовлетворяют уравнению неразрывности и не требуется никакой коррекции давления. Таким образом, член b представляет собой «источник массы», который должен быть скомпенсирован поправками давления (через соответствующие поправки скорости).

Итак, мы получили все уравнения, необходимые для определения составляющихскорости и давления. Теперь можно рассмотреть весь алгоритм решения в целом.

Алгоритм SIMPLE

Процедура, разработанная для расчета поля течения, получила название SIMPLE(Semi-Implicit Method for Pressure-Linked Equations), что означает полунеявный метод для связывающих давление уравнений. Обсудим смысл этого названия несколько позднее.

Последовательность операций. Ниже перечислены основные операции в порядке их выполнения.

1. Задание поля давления р*.

2. Решение уравнений движения, таких, как уравнения (3.2.3), для получения и*,v*, w*.

3. Решение уравнения для р'.

4. Расчет р из уравнения (3.3.1) путем добавления р' к р*.

5. Расчет и, v, w с учетом соответствующих значений со звездочкой и с помощью формул для поправки скорости (3.3.7)-(3.3.8).

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

7. Представление скорректированного давления р как нового р*, возвращение к пункту 2 и повторение всей процедуры до тех пор, пока не будет получено сходящееся решение.