Эта публикация цитируется в
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