RUS  ENG
Полная версия
ЖУРНАЛЫ // Компьютерные исследования и моделирование // Архив

Компьютерные исследования и моделирование, 2020, том 12, выпуск 6, страницы 1291–1306 (Mi crm849)

Эта публикация цитируется в 5 статьях

ЧИСЛЕННЫЕ МЕТОДЫ И ОСНОВЫ ИХ РЕАЛИЗАЦИИ

Метод численного решения одной стационарной задачи гидродинамики в конвективной форме в $L$-образной области

В. А. Рукавишниковa, А. В. Рукавишниковb

a Вычислительный центр ДВО РАН, Россия, 680000, г. Хабаровск, ул. Ким ЮЧена, д. 65
b Хабаровское отделение Института прикладной математики ДВО РАН, Россия, 680000, г. Хабаровск, ул. Дзержинского, д. 54

Аннотация: Большой класс задач описывает физические процессы, протекающие в невыпуклых областях, содержащих угол больший 180 градусов на границе. Решение в окрестности такого угла сингулярно, а его отыскание, при использовании классических подходов, влечет за собой потерю точности. В представленной работе рассмотрены стационарные, линеаризованные с помощью итераций Пикара несжимаемые уравнения Навье–Стокса течения вязкой жидкости в конвективной форме в $L$-образной области. Определено $R_{\nu}$-обобщенное решение задачи в специальных множествах весовых пространств. Для нахождения приближенного $R_{\nu}$-обобщенного решения построен специальный метод конечных элементов. Во-первых, пространства конечно-элементных функций удовлетворяют закону сохранения массы в сильном смысле, то есть в узлах сетки. Для этой цели используется Скотт–Вогелиус конечно-элементная пара. Выполнение закона сохранения массы ведет к отысканию более точного с физической точки зрения решения.Во-вторых, базисные функции конечномерных пространств дополнены весовыми функциями как множителями, которые совпадают с расстоянием от точки до вершины тупого угла в $\delta$-окрестности точки сингулярности и радиусом $\delta$ вне ее. Степень весовой функции, как и параметр $\nu$ в определении $R_{\nu}$-обобщенного решения, так и радиус $\delta$-окрестности точки сингулярности являются свободными параметрами метода. Специально подобранная их комбинация приводит к увеличению порядка сходимости приближенного решения к точному решению задачи почти в два раза по сравнению с классическими подходами и достигает единицы по шагу сетки в нормах весовых пространств Соболева. Таким образом, установлено, что скорость сходимости не зависит от величины угла.

Ключевые слова: задача гидродинамики с сингулярностью, весовой метод конечных элементов.

УДК: 519.6

Поступила в редакцию: 20.01.2020
Исправленный вариант: 16.02.2020
Принята в печать: 10.11.2020

DOI: 10.20537/2076-7633-2020-12-6-1291-1306



© МИАН, 2024