Fletcher-2-rus (1185919), страница 68
Текст из файла (страница 68)
Глава 17 Несжимаемые вязкие течения В этой главе не будет делаться никаких предположений об относительной величине компонент скорости. Следовательно, укороченные уравнения Навье — Стокса (гл. 16) становятся неприменимы и необходимо рассматривать полные уравнения Навье — Стокса; при этом предполагается, что течение несжимаемое. Численные методы, рассматриваемые в настоящей главе, будут применяться к задачам, в которых нет выделенного направления течения, например задача о вентиляции помещений.
Кроме того, во многих случаях будут возникать области возвратного течения. Если эти области велики или их появление связано с нестационарностью течения (например, течение за уступом), маршевые методы (итерационные), рассмотренные в гл. 16, неприменимы для решения подобных задач. Предположение о несжимаемости течения приводит к дополнительным вычислительным трудностям.
В уравнение неразрывности (! 1.13) входят лишь компоненты скорости. Поэтому в данном случае нет прямой связи с давлением, которая в случае сжимаемых течений осуществляется через плотность. Для расче~а несжимаемых течений возможны два общих подхода. В первом из них используются исходные переменные и, и, р в двумерном случае, а для решения уравнения неразрывности вводятся специальные процедуры. При обобщении этого подхода на случай трехмерных течений не возникает дополнительных трудностей. Соответствующие методы расчета нестационарных течений, основанные на исходных переменных, описаны в ~ 17.1.
Методы в исходных переменных, предназначенные для расчета стационарных течений, описаны в ~ 17.2. В двумерном случае явного использования уравнения неразрывности можно избежать, если ввести в рассмотрение функцию тока. Кроме того, введение уравнения переноса для завихренности позволяет получить описание течения в переменных завихренность — функция тока. Такой подход описан в 17.3. Обобщение этого подхода на трехмерный случай не столь очевидно, поскольку в случае трех пространственных переменных функция тока не существует. Для описания течения Гл.
!7. Несжимаемые вязкие течения в этом случае используются завихренность и векторный потенциал (п. 17.4.1). Большинство практически интересных течений являются турбулентными, если только число Рейнольдса не слишком мало. Для учета турбулентных эффектов при расчете несжимаемых течений обычно используется либо алгебраическая модель турбулентной вязкости (п. 11.4.2), либо (й — е)-модель (п. 11.5.2). С вычислительной точки зрения использование алгебраической модели требует лишь незначительного изменения алгоритмов, предназначенных для расчета ламинарных течений.
Структура дифференциальных уравнений для й и е (11.95) и (11.96) аналогична структуре уравнений импульса, и дискретизация этих уравнений обычно проводится одинаково. Таким образом, алгоритмы расчета ламинарных вязких течений столь же эффективны (с небольшими изменениями) и для расчета турбулентных течений. Поэтому в настоящей главе явное внимание расчетам турбулентных течений уделяться не будет. В 17.1.
Исходные переменные: нестационарные течения Уравнения, описывающие двумерные нестационарные несжимаемые ламинарные течения, имеют вид — + — =О, ди до дх ду (17.1) ди диз д др ! ! даи дзо Ъ вЂ” + — + — (ио) + — = — ! — + — 1, (17.2) д! дх -ду дх Ке 1 дхз дуя)' до д доз др 1 Г дзо дзо Ъ вЂ” + — (ио) + — + — = — ~ —, + —, ~ . (17.3) д! дх ду ду Ке 1 дхз дуз 1' Уравнения (!7.!) — (17.3) записаны в безразмерном виде; плотность включена в число Рейнольдса Ке. При помощи уравнения (17.1) левые части (17.2) и (!7.3), как в уравнении (11.116), записаны в консервативном виде. Для нестационарных течений требуется определить начальные условия и = ив(х, у) и о = ое(х, у), удовлетворяющие уравнению (17.1).
Граничные условия на твердой поверхности заключаются в отсутствии относительного движения жидкости и твердого тела, что определяет компоненты скорости. Граничных условий для давления на твердой поверхности задавать не надо. Если компоненты скорости определены на всей границе области расчета, как, например, в задаче о движущейся полости, необходимо обеспечить выполнение глобального условия ~ у. пс(з=О, (17.4) с $17.1. Исходные переменные: иестапиоиариые течения 391 где с — граница области расчета. Уравнение (17.4) является глобальным выражением (!7.1), что может быть получено из сравнения уравнений (11.7) и (11.10) при постоянном значении плотности р. Если область расчета содержит открытые границы, как в задаче о течении вблизи уступа (п.
17.3.3), число граничных условий иа открытых границах может быть получено из табл. 11.5. На входной границе необходимо поставить два граничных условия; обычно задаются одна компонента скорости и давление. На выходной границе можно положить равными нулю производные от скорости, граничное условие для давления ставить не надо. Поскольку в уравнения входят лишь производные от давления, его величина может быть определена в одной точке, относительно которой будет осуществляться отсчет давления. Следует подчеркнуть, что описанная выше постановка граничных условий проведена так, что уравнения и граничные условия образуют корректно поставленную задачу, имеющую регулярное решение. Однако для получения регулярного численного решения дискретных уравнений может понадобиться введение дополнительных граничных условий. Описанные в этом параграфе методы основаны главным образом на конечио-разностной дискретизации и решении уравнения Пуассона для давления (п.
17.1.2). Чтобы получить достаточно точное решение без чрезмерного измельчения сетки, может понадобиться более точное дискретное представление конвективных членов (п. 17.1.5). Многие из описанных в данном разделе способов решения, использующие конечно-разностную дискретизацию, могут применяться и при ином способе дискретизации, например в случае спектральных методов (п. 17.1.6). 17.1.1. Разнесенная сетка Численное решение системы (!7.1) — (17.3) часто проводится на разнесенной сетке. Это означает, что различные зависимые переменные определяются в разных точках сетки. В работе 1Реуге(, Тау!ог, 1983) проведено сравнение различных разнесенных сеток применительно к определению давления. В работе [Най!пег, %!1!!агпз, 1980) рассмотрена возможность представления различных мод Фурье (см.
п. 9.2.1) на разнесенных сетках различной конфигурации для уравнений мелкой воды, аналогичных уравнениям Эйлера (гл. 14), Лучшая конфигурация разнесенной сетки представлена иа рис. !7.1. Видно, что давление определяется в центре ячейки, а компоненты скорости — на границах. Такая процедура делает сетку удобной для проведения дискретизации по методу конечных Гл. 17. Несжимаемые вязкие течения 392 объемов ($ 5.2); эта связь будет использована в п.
17.2.3. В результате дискретизации уравнения (17.1) на разнесенной сетке, изображенной на рис. 17.1, получается выражение + ' ~'в ' и — О, (!7.5) Ьх аи которое можно представить в виде изап, а Ау+ одачпзбх — и! Оз аЛу — оь а ьоЛх =О. (!7 6) Уравнение (17.6) является дискретным представлением уравнения (17.4), т. е.
соотношения (17.5) и (17.6) сохраняют массу Рис. 17.1. Разнесенная сетка. на минимальном сеточном размере. Кроме того, из разложения (17.5) в ряд Тейлора в окрестности центра ячейки следует, что порядок аппроксимации (17.5) равен 0(Лхз, Луи) несмотря на использование лишь четырех точек. Дискретизация уравнений (17.2) проводится с помощью конечно-разностных выражений, центрированных относительно точки сетки (1 + 1/2, й). Это позволяет представить др/дх в виде выражения (р;ч.ь а — рс е)/Лх, которое в точке (1+!/2,и) имеет второй порядок точности.
Аналогично (!7.3) дискретизируется центральными разностями относительно точки (1', й+ 1/2) и др/ду представляется в виде (рь ечч †Рь «)/Ау. Использование разнесенной сетки дает возможность связать значения и, и и р в соседних точках. Это также позволяет избежать появления осцилляций в решении, в частности для р, которые могут возникнуть, если центральные разности используются для аппроксимации всех производных на неразнесенной сетке. Осциллирующее решение появляется из-за двух несвязанных на различных точках сетки решений для давления, которое возможно, если центральные разности используются на неразнесенной сетке.
Осцилляции, как правило, возрастают при увеличении числа Рейнольдса, поскольку диссипативные члены $ !7.1. Исходные переменные: нестапнонарные течения 393 посредством которых осуществляется связь значений и и и в соседних точках, в этом случае малы. Из уравнений (17.1)— (17.3) следует, что членов, приводящих к диссипации р, не имеется. Дополнительное преимущество использования разнесенной сетки состоит в том, что уравнение Пуассона (!7.13) для давления автоматически удовлетворяет дискретному представлению интегрального граничного условия (17.4). В этом случае не требуется, как в уравнении (16.98), проводить дополнительную коррекцию правой части уравнения Пуассона. Применение разнесенных сеток имеет также некоторые недостатки.
Компьютерные программы в этом случае труднее интерпретировать. При программировании необходимо связать полный набор независимых переменных с индексами массивов, в которых они хранятся. При использовании разнесенных сеток элементам массивов, в которых хранятся и, и и р с номером (1, ц), могут соответствовать значения и,ччтв е, по е+ия и рь а (рис. 17.1). В случае разнесенной сетки обычно труднее осуществить постановку граничных условий, поскольку по крайней мере одна из зависимых переменных и или и не будет определена на границе. Если используются обобщенные координаты (гл. 12) и сетка непрямоугольная, применение разнесенных сеток усложняется еще больше. Разнесенная сетка, приведенная на рис. 17.1, используется в методе МАС (п.