Материал из РУВИКИ — свободной энциклопедии

Разрывный метод Галёркина

Разрывный метод Галёркина (англ.  discontinuous Galerkin method, сокращенно DGM) — метод решения операторных уравнений, в основном дифференциальных уравнений. Является развитием классического метода конечных элементов (МКЭ), основанного на вариационной постановке Галёркина.

История метода[править | править код]

Разрывный метод Галёркина был впервые предложен и в начале 70-ы годов XX века, как метод решения дифференциальных уравнений в частных производных, в 1973 году Ридом (англ.) и Хиллом был предложен вариант метода для решения гиперболического уравнения переноса нейтронов. Первая формулировка метода для решения эллиптических задач не может быть определена единственной публикацией, однако на развитие метода оказали сильное влияние Иво Бабушка (англ.) и Жак-Луи Лионс (англ.). Для уравнений 4 порядка вариант метода представил Бэйкер в 1977 году. Так же своим развитием метода обязан публикациям Арнольди, Бреззи, Кокбёрна и Марини.

Основные определения[править | править код]

Конечным элементом называется тройка пространств , где:

  •  — геометрический конечный элемент — множество точек пространства (часто, говоря «конечный элемент», подразумевают именно это множество);
  •  — множество координатных функций — базис пространства, которое будет аппроксимировать решение;
  •  — множество степеней свободы — базис сопряженного пространства. .

Идея метода и отличие от МКЭ[править | править код]

Рассмотрим идею метода для решения ДУ второго порядка в области . В отличие от метода Галёркина, где выполняется постановка в слабой форме, в DGM выполняется постановка в слабой слабой форме (ультра слабой форме, англ. ultra weak variational formulation).
Представим исходное уравнение в виде двух уравнений первого порядка. В зависимости от природы уравнений это можно сделать несколькими способами, которые приведут к различным вариационным постановкам. Далее на расчётной области строим сетку , выполняем вариационную постановку Галёркина для каждой подобласти при этом будет использоваться четыре пространства: два пространства(координатное и проекционное) для самой функции и два для её производной. После этого уравнения суммируются по всей области и из получившихся системы из двух уравнений каким-нибудь способом исключается одно.
Данное описание является весьма общим и неоднозначным, поскольку метод всегда подстраивается под конкретные задачи и получение ультра слабой вариационной постановки зависит от природы процесса и цели решения уравнения.

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

Примеры для конкретных видов уравнений[править | править код]

Уравнение теплопроводности[править | править код]

Рассмотрим простейший случай стационарного уравнения теплопроводности:



 — коэффициент теплопроводности,  — правая часть уравнения.
Выполним замену и тем самым сведём уравнение второго порядка к двум уравнениям первого порядка:


На расчётной области введём пространство Лебега с соответствующим ему скалярным произведением: . И соответствующие ему конечноэлементные пространства:
 — пространство скалярных функций, для аппроксимации решения
 — пространство векторных функций для аппроксимации градиента решения
Введённые пространства являются пространствами Соболева (скалярным и векторным) с соответствующей нормой. Из этих пространств выберем тестовые функции и для каждого уравнения выполним постановку Галёркина на отдельном элементе, получим систему уравнений в слабой форме[1]:


Функции  — это численные потоки, которые могут быть определены по-разному (что ведёт к различным методам) и должны удовлетворять следующим условиям:

  • Численные потоки однозначны на границе
  • Численные потоки представимы в виде: , где  — гладкая функция удовлетворяющая заданным первым краевым условиям.

Для упрощения записи вводят оператор среднего и оператор скачка, определяющие поведение функций на границе элементов:

Теперь просуммируем все полученные уравнения для каждой подобласти и получим два уравнения для всей области:


Воспользуемся свойством[3]: и в результате получим ультра слабую вариационную постановку для исходного уравнения:

Осталось определить численные потоки. Определение численных потоков связанно с задачей и требованиями к решению и ведёт к различным методам, например:

Функция и область IP-метод[4] Стабилизированный IP-метод NIPG[5]
на
на
на и на


Уравнения Максвелла в гармоническом режиме[править | править код]

Подход к построению ультра слабой вариационной постановки для уравнений Максвелла может быть различным: систему уравнений первого порядка можно получить напрямую из самих уравнений Максвелла или сведя эти уравнения к уравнению Гельмгольца, а затем сделав замену, аналогичную замене для уравнения теплопроводности, получить систему первого порядка. В данном случае мы воспользуемся первым способом. Система уравнений Максвелла в гармоническом режиме с частотой , в одном из простейших случаев выглядит как:


Оба уравнения выполняются в расчётной области . Краевые условия: . Домножим скалярно оба уравнения на пробные функции , определённые на соответствующем элементе . Функции из этого же пространства будут использоваться как базисные. Для их определения используем присоединённую систему уравнений Максвелла[6]:


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


Учитывая систему уравнений для пробных функций данное выражение упрощается до:


Введём обозначения:

Вектора Матрицы




Теперь задача ставится как нахождение векторов для всех элементов, удовлетворяющих следующим уравнениям [6]:


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

Литература[править | править код]

  • Arnold D.N., Brezzi F., Cocburn B., Mariri D. «Unified analysis of DCM for elliptic problems», 2002

Примечания[править | править код]

  1. Ю. И. Шокин, Э. П. Шурина, Н. Б. Инткина. Современные многосеточные методы. — НГТУ, 2012. — 98 с. — ISBN 978-5-7782-2119-2.
  2. Значения функций берутся как предел к границе изнутри области, то есть , где - единичная внешняя нормаль.
  3. Arnold D.N., Brezzi F., Cocburn B., Mariri D. Unified analysis of DCM for elliptic problems. — 2002.
  4. англ. Internal penalty, метод внутренних штрафов, метод Бауманна-Одена
  5. Не симметричный IP-метод
  6. 1 2 T. Huttunen, M. Malinen, P. Monk. Solving Maxwell’s Equations using Ultra Weak Variational Formulation (англ.). — 2006.