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

03-10-2023

Перейти к: навигация, поиск

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

История метода

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

Основные определения

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

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

Идея метода и отличие от МКЭ

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

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

Примеры для конкретных видов уравнений

Уравнение теплопроводности

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



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

\left \{ \begin{array}{lcr} \nabla u = \sigma \\ - \nabla \cdot \sigma = f \\ \end{array} \right.


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

\left \{ \begin{array}{lcr} \int_{K_i} \sigma \cdot \tau dx \ = - \int_{K_i} u \nabla \cdot \tau dx \ + \int_{\partial K_i} \hat u_{K_i} (n_{K_i} \cdot \tau) ds\ ~\forall \tau \in \Sigma_i\\ \int_{K_i} \sigma \cdot \nabla v dx\ = \int_{K_i} fvdx\ + \int_{\partial K_i} (\hat \sigma_{K_i} \cdot n_{K_i}) v ds\ ~\forall v \in V_i\\ \end{array} \right.


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

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

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

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

\left \{ \begin{array}{lcr} \int_{\Omega} \sigma \cdot \tau dx \ = - \int_{\Omega} u \nabla \cdot \tau dx \ + \sum_{K_i \in K}\int_{\partial K_i} \hat u_{K_i} (n_{K_i} \cdot \tau) ds\ ~\forall \tau \in \Sigma_h\\ \int_{\Omega} \sigma \cdot \nabla v dx\ = \int_{\Omega} fvdx\ + \sum_{K_i \in K} \int_{\partial K_i} (\hat \sigma_{K_i} \cdot n_{K_i}) v ds\ ~\forall v \in V_h\\ \end{array} \right.


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

\int_{\Omega} \nabla u \cdot \nabla v dx\ + \sum_{K_i} \int_{\partial K_i}{ [\hat u - u] \cdot \left \langle \nabla v \right \rangle ds} + \sum_{K_i} \int_{\partial K_i}{\left \langle \hat u - u \right \rangle [\nabla v] ds} = \int_{\Omega} {fv dx} + \sum_{K_i} \int_{\partial K_i}{ [v] \left \langle \hat \sigma \right \rangle ds} + \sum_{K_i} \int_{\partial K_i}{ \left \langle v \right \rangle [\sigma] ds},~\forall v \in V_h

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

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


Уравнения Максвелла в гармоническом режиме

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

\left \{ \begin{array}{lcr} -i\omega \epsilon \mathbf{E} - \nabla \times \mathbf{H} = 0\\ -i\omega \mu \mathbf{H} + \nabla \times \mathbf{E} = 0 \\ \end{array} \right.


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

\left \{ \begin{array}{lcr} i\omega \epsilon \mathbf{\xi^{K_i}} + \nabla \times \mathbf{\psi^{K_i}} = 0\\ i\omega \mu \mathbf{\psi^{K_i}} - \nabla \times \mathbf{\xi^{K_i}} = 0 \\ \end{array} \right.


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

\int_{K_i} { \left( \mathbf {E}\cdot \overline{i\omega \epsilon \mathbf{\xi^{K_i}} + \nabla \times \mathbf{\psi^{K_i}}} + \mathbf{H} \cdot \overline{i\omega \mu \mathbf{\psi^{K_i}} - \nabla \times \mathbf{\xi^{K_i}}} \right) dx} = \int_{\partial K_i} { \left( \mathbf{n^{K_i}} \times \mathbf{H} \cdot \mathbf{\xi^{K_i}} - \mathbf{n^{K_i}} \times \mathbf{E} \cdot \mathbf{\psi^{K_i}} \right) ds}


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

\int_{\partial K_i} { \left( \mathbf{n^{K_i}} \times \mathbf{H} \cdot \mathbf{\xi^{K_i}} - \mathbf{n^{K_i}} \times \mathbf{E} \cdot \mathbf{\psi^{K_i}} \right) ds} = 0


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

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




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

\int_{\partial K_i} {\chi^{K_i} \cdot \overline{\gamma^{K_i}} ds} - \sum_{\Gamma_{ij} = \partial K_i \cap \partial K_j \neq \varnothing} \int_{\Gamma_{ij}} { \chi^{K_j} \cdot \overline{F_{K_i}(\gamma^{K_i})}ds} - \sum_{\Gamma_{i} = \partial K_i \cap \partial \Omega \neq \varnothing} \int_{\Gamma_{i}} { \chi^{K_i} \cdot \overline{F_{K_i}(\gamma^{K_i})}ds} = \int_{\partial K_i} {\mathbf{g} \cdot \overline{F_{K_i}(\gamma^{K_i})} ds} ~~ \forall \gamma^{K_i}


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

Литература

  • Ю. И. Шокин, Э. П. шурина, Н. Б. Инткина «Современные многосеточные методы», 2012
  • Arnold D.N., Brezzi F., Cocburn B., Mariri D. «Unified analysis of DCM for elliptic problems», 2002
  • T. Huttunen, M. Malinen, P. Monk «Solving Maxwell’s Equations using Ultra Weak Variational Formulation», 2006

Примечания

  1. ↑ Значения функций берутся как предел к границе изнутри области, то есть , где - единичная внешняя нормаль.
  2. ↑ Arnold D.N., Brezzi F., Cocburn B., Mariri D. «Unified analysis of DCM for elliptic problems», 2002
  3. ↑ англ. Internal penalty, метод внутрених штрафов, метод Бауманна-Одена
  4. ↑ Не симметричный IP-метод

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

© 2011–2023 stamp-i-k.ru, Россия, Барнаул, ул. Анатолия 32, +7 (3852) 15-49-47