Метод граничных элементов (МГЭ): снижение размерности задачи на единицу

Снижение размерности: от объема к поверхности

Метод конечных элементов (МКЭ) и метод конечных разностей (МКР) безусловно доминируют в современных инженерных расчетах. Однако оба эти метода являются полностью объемными: они требуют построения густой расчетной сетки во всей 3D-области математического моделирования. Если вам нужно рассчитать обтекание подводной лодки в открытом океане или напряжение вокруг микротрещины в огромном стальном блоке, объемная сетка съест колоссальное количество памяти суперкомпьютера. В таких специфических случаях на помощь приходит элегантный математический аппарат — Метод граничных элементов (МГЭ, или Boundary Element Method, BEM).

Основная геометрическая идея МГЭ заключается в снижении размерности задачи ровно на одну физическую единицу. Для массивного трехмерного тела сетка строится только на его внешней двумерной поверхности (оболочке). Для плоского двумерного сечения крыла самолета сетка представляет собой просто одномерную кривую, очерчивающую внешний контур. Это радикально, на несколько порядков, уменьшает количество узлов и неизвестных переменных, позволяя решать задачи высочайшей сложности на обычных ПК. Но возникает парадокс: как можно математически рассчитать то, что происходит внутри объема, зная только числовые данные на его поверхности?

Интегральные уравнения и фундаментальные решения

Математический секрет кроется в глубоких теоремах классического векторного анализа (таких как формулы Грина, теорема Бетти и принцип Сомильяны). Метод граничных элементов напрямую применим только к тем дифференциальным уравнениям (обычно это линейные однородные уравнения), для которых в математике аналитически известно так называемое фундаментальное решение. Фундаментальное решение (функция Грина) описывает, как среда реагирует на бесконечно малый точечный источник возмущения (например, на единичную силу или единичный точечный электрический заряд).

Математически дифференциальное УЧП строго переписывается в виде граничного интегрального уравнения. Интегрирование ведется исключительно по поверхности (границе) физического тела. В каждом узле поверхностной сетки ставится неизвестное значение (либо самой искомой функции, либо ее нормальной производной — физического потока). В результате численной дискретизации интеграла получается полностью заполненная (плотная), несимметричная система линейных алгебраических уравнений (СЛАУ). Решив эту граничную СЛАУ методом Гаусса, мы находим все недостающие неизвестные на поверхности. После этого значение температуры, механического напряжения или электрического потенциала в абсолютно любой внутренней точке объема может быть моментально вычислено прямым аналитическим интегрированием уже известных граничных данных. МГЭ является абсолютным чемпионом в акустике, сейсмологии и задачах с бесконечными областями.