Метод конечных элементов: решение краевых задач на сложных областях

Метод решает краевые задачи двадцать третьей статьи на произвольной геометрии.

Вариационная формулировка и слабое решение

Задача (Delta u = f) сводится к (int abla u cdot abla v dx = int f v dx) для тестовых функций (v). Пространство Соболева (H^1) обеспечивает существование слабого решения. Галёркин: аппроксимация в конечномерном подпространстве (V_h).

Разбиение сложной области на треугольники с линейными базисными функциями на каждом элементе . Схема показывает локальную аппроксимацию.

Сборка глобальной матрицы жесткости

Локальная матрица элемента (K^e_{ij} = int_e abla phi_i cdot abla phi_j dx) собирается в глобальную систему (K mathbf{u} = mathbf{f}). Линейные элементы дают треугольную матрицу, квадратичные — лучшую точность. Решение использует методы пятой статьи.

Теория МКЭ содержится в лекциях MIT https://ocw.mit.edu/courses/18-100b-real-analysis-spring-2025/video_galleries/video-lectures/. Подробности в НГУ http://www.phys.nsu.ru/podvigin/Lectures.pdf.

Адаптивное h-p уточнение сетки

Априорные оценки ошибок (||u - u_h|| leq C h^p) определяют локальное уточнение сетки. hp-метод комбинирует уточнение сетки и повышение степени полинома. Апослериорные индикаторы ошибок управляют адаптацией.

Последовательное уточнение сетки в областях больших градиентов с оценкой ошибки . График показывает эффективность.

Связь с другими методами

При прямоугольной сетке МКЭ совпадает с конечными разностями двадцать восьмой статьи. Спектральные элементы комбинируют МКЭ и спектральные методы двадцать девятой статьи.

Рекомендуемые книги по МКЭ

Полезны материалы по функциональному анализу http://old.math.nsc.ru/~matanalyse/basic2.pdf и векторному анализу https://math.csu.ru/new_files/students/lectures/analit_geom/matveev_posobie_po_vektr_algebra.pdf. Источники содержат приложения.