2025-11-10T03:15:01.076525

Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method

Adjerid, Lin, Meghaichi
The Frenet apparatus is a new framework for constructing high order geometry-conforming immersed finite element functions for interface problems. In this report, we present a procedure for constructing the local IFE bases in some detail as well as a new approach for constructing orthonormal bases using the singular value decomposition of the local generalized Vandermonde matrix. A sample implementation in MATLAB is provided to showcase the simplicity and extensionability of the framework.
academic

Построение базисных функций для метода конформной к геометрии погруженного конечного элемента

Основная информация

  • ID статьи: 2510.12018
  • Название: Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method
  • Авторы: Slimane Adjeril (Virginia Tech), Tao Lin (Virginia Tech), Haroun Meghaichi (Ohio State University)
  • Классификация: math.NA cs.NA
  • Дата публикации: 15 октября 2025
  • Ссылка на статью: https://arxiv.org/abs/2510.12018

Аннотация

В данной работе предложена новая схема, основанная на аппарате Френе, для построения высокопорядковых функций конформного к геометрии погруженного конечного элемента (GC-IFE) для решения задач на интерфейсах. Статья подробно описывает процесс построения локальных базисных функций IFE, а также новый метод построения ортонормированных базисов с использованием разложения по сингулярным значениям (SVD) локальной обобщённой матрицы Вандермонда. Приведены примеры реализации на MATLAB, демонстрирующие простоту и масштабируемость предложенной схемы.

Научный контекст и мотивация

Постановка проблемы

Данное исследование решает эллиптическую задачу на интерфейсе: −∇⋅(β∇u)=f, на Ω−∪Ω+-\nabla \cdot (\beta\nabla u) = f, \text{ на } \Omega^- \cup \Omega^+ где расчётная область Ω⊂R2\Omega \subset \mathbb{R}^2 разделена интерфейсом Γ\Gamma на два подмножества Ω−\Omega^- и Ω+\Omega^+, а коэффициент диффузии β\beta является кусочно-постоянной функцией.

Научная мотивация

  1. Ограничения существующих методов IFE: традиционные методы погруженного конечного элемента могут только приближённо удовлетворять условиям скачка на интерфейсе
  2. Требование высокопорядковой точности: необходимость построения функций IFE произвольного полиномиального порядка
  3. Требование конформности к геометрии: необходимость функций IFE, точно удовлетворяющих условиям скачка на интерфейсе
  4. Вычислительная устойчивость: необходимость системы базисных функций с хорошим числом обусловленности

Значимость

Задачи на интерфейсах широко встречаются в многофазных течениях, материаловедении, биомедицине и других областях. Высокоточные численные методы их решения имеют важное значение для инженерных приложений.

Основные вклады

  1. Предложена схема построения GC-IFE на основе аппарата Френе: использование репера Френе из дифференциальной геометрии для упрощения условий скачка на интерфейсе
  2. Разработаны два метода построения базисных функций:
    • метод начального построения на основе специального формата
    • обобщённый метод построения на основе идеи расширения
  3. Предложены два алгоритма ортогонализации:
    • метод на основе SVD матрицы масс (Подход 1)
    • метод на основе SVD обобщённой матрицы Вандермонда (Подход 2)
  4. Предоставлена полная реализация на MATLAB: демонстрирует практичность и масштабируемость метода

Описание методики

Преобразование Френе

Основная идея заключается в использовании преобразования Френе для преобразования физических координат (x,y)(x,y) в локальные координаты (η,ξ)(\eta,\xi): [ηξ]→PΓ[x(η,ξ)y(η,ξ)]=g(ξ)+ηn(ξ)\begin{bmatrix} \eta \\ \xi \end{bmatrix} \xrightarrow{P_\Gamma} \begin{bmatrix} x(\eta,\xi) \\ y(\eta,\xi) \end{bmatrix} = g(\xi) + \eta n(\xi)

где g(ξ)g(\xi) — параметризация кривой интерфейса, n(ξ)n(\xi) — вектор нормали.

Преобразование условий скачка на интерфейсе

В локальной системе координат условия скачка на интерфейсе принимают вид:

  1. ⟦u^⟧Γ^KF=0\llbracket \hat{u} \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  2. ⟦β^u^η⟧Γ^KF=0\llbracket \hat{\beta}\hat{u}_\eta \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  3. ⟦β^∂j∂ηjL(u^)⟧Γ^KF=0,j=0,1,…,m−2\llbracket \hat{\beta}\frac{\partial^j}{\partial\eta^j}\mathcal{L}(\hat{u}) \rrbracket_{\hat{\Gamma}_{K_F}} = 0, j = 0,1,\ldots,m-2

Методы построения базисных функций

Метод 1: построение на основе специального формата

Построение базисных функций вида: ϕ^i,j(η,ξ)=1β^(η,ξ)qj(ηηh)pi(ξ−ξmidξh)\hat{\phi}_{i,j}(\eta,\xi) = \frac{1}{\hat{\beta}(\eta,\xi)}q_j\left(\frac{\eta}{\eta_h}\right)p_i\left(\frac{\xi-\xi_{mid}}{\xi_h}\right)

Для случаев 1≤j≤m,0≤i≤m1 \leq j \leq m, 0 \leq i \leq m эти функции автоматически удовлетворяют условиям скачка.

Метод 2: обобщённое построение с расширением

Использование более общей формы:

\hat{\lambda}_j^-(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^- R_i(\eta,\xi), & \eta < 0 \\ \hat{\lambda}_j^+(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^+ R_i(\eta,\xi), & \eta > 0 \end{cases}$$ Коэффициентные матрицы $C^{\pm}$ определяются путём решения системы линейных уравнений. ### Ортогонализация и реконструкция #### Подход 1: SVD матрицы масс Пусть $M_q = V_1\Lambda V_1^T$, выбираем $Q_1 = V_1\Lambda^{-1/2}$ так, чтобы реконструированная матрица масс была единичной матрицей. #### Подход 2: SVD обобщённой матрицы Вандермонда Выполняем SVD разложение $\tilde{V} = \sqrt{W}V$: $\tilde{V} = U_2\Sigma V_2^T$, выбираем $Q_2 = V_2\Sigma^{-1}$. Теоретический анализ показывает, что Подход 2 более устойчив при численной реализации, поскольку избегает SVD разложения матрицы с квадратом числа обусловленности. ## Экспериментальная установка ### Тестовые задачи Рассматривается область $\Omega = [-1,1]^2$, интерфейс представляет собой $x^2 + y^2 = r_0^2$ ($r_0 = 1/\sqrt{3}$), тестовая функция: $$u(x,y) = \begin{cases} \frac{1}{\beta^+}\cos(2\pi r^2), & r > r_0 \\ \frac{1}{\beta^-}\cos(2\pi r^2) + \cos(2\pi r_0^2)\left(\frac{1}{\beta^+} - \frac{1}{\beta^-}\right), & r < r_0 \end{cases}$$ ### Критерии оценки - порядок сходимости ошибки проекции в норме $L^2$ - число обусловленности матрицы масс - численная устойчивость построения базисных функций ## Результаты экспериментов ### Проверка сходимости Для различных полиномиальных порядков $m = 1,2,3,4$ ошибка проекции в норме $L^2$ демонстрирует оптимальный порядок сходимости: | $m$ | $N = 16$ | $N = 32$ | $N = 64$ | $N = 128$ | $N = 256$ | Порядок сходимости | |-----|----------|----------|----------|-----------|-----------|--------| | 1 | 8.14E-02 | 2.08E-02 | 5.23E-03 | 1.31E-03 | 3.28E-04 | ≈2.0 | | 2 | 9.29E-03 | 1.19E-03 | 1.50E-04 | 1.88E-05 | 2.36E-06 | ≈3.0 | | 3 | 8.80E-04 | 5.74E-05 | 3.62E-06 | 2.27E-07 | 1.42E-08 | ≈4.0 | | 4 | 7.55E-05 | 2.39E-06 | 7.49E-08 | 2.34E-09 | 7.33E-11 | ≈5.0 | ### Эффект улучшения числа обусловленности Сравнение чисел обусловленности матрицы масс до и после реконструкции показывает: - число обусловленности при начальном построении растёт экспоненциально с порядком полинома - Подход 2 стабильно снижает число обусловленности до 1.0 - Подход 1 демонстрирует численную неустойчивость при высоких порядках ### Эффект предобусловливания Для матрицы коэффициентов $A$ и $\tilde{A}$ предобусловливатели Якоби и строчной нормализации значительно улучшают число обусловленности, особенно при высоких порядках и малых размерах сетки. ## Связанные работы ### Развитие методов IFE - традиционные методы IFE: могут только приближённо удовлетворять условиям на интерфейсе - расширенные методы IFE: высокопорядковые методы на основе расширения Коши - GC-IFE в данной работе: методы конформные к геометрии, точно удовлетворяющие условиям на интерфейсе ### Методы ортогонализации - построение ортогональных полиномов в методах виртуальных элементов - ортогональные базисные функции в методах разрывного Галёркина - применение идей ортогонализации к схеме IFE в данной работе ## Заключение и обсуждение ### Основные выводы 1. **Теоретический вклад**: установлена теоретическая схема GC-IFE на основе аппарата Френе 2. **Алгоритмические инновации**: предложены два метода построения и две стратегии ортогонализации 3. **Численная верификация**: доказаны оптимальная сходимость и численная устойчивость метода 4. **Практическая ценность**: предоставлена полная реализация на MATLAB ### Ограничения 1. **Геометрические ограничения**: требуется достаточная гладкость кривой интерфейса ($C^3$) 2. **Вычислительная сложность**: высокие затраты на сборку матриц при высоких порядках 3. **Масштабируемость**: в настоящее время применимо только к двумерным эллиптическим задачам на интерфейсах ### Направления будущих исследований 1. Расширение на трёхмерные задачи и сложные геометрии 2. Применение к зависящим от времени задачам на интерфейсах 3. Расширение на задачи теории упругости, системы Стокса и другие 4. Разработка более эффективных алгоритмов реализации ## Глубокая оценка ### Преимущества 1. **Теоретическая строгость**: аппарат Френе из дифференциальной геометрии обеспечивает прочную теоретическую базу 2. **Методологические инновации**: впервые реализован высокопорядковый метод IFE, точно удовлетворяющий условиям скачка на интерфейсе 3. **Численная устойчивость**: ортогонализация и реконструкция значительно улучшают вычислительную устойчивость 4. **Полнота реализации**: предоставлены подробные описания алгоритмов и код на MATLAB 5. **Достаточная верификация**: теоретический анализ и численные эксперименты взаимно подтверждают друг друга ### Недостатки 1. **Ограниченная область применения**: в настоящее время ограничена двумерными эллиптическими задачами на интерфейсах 2. **Анализ вычислительной сложности**: анализ сложности вычислений при высоких порядках недостаточно глубок 3. **Строгие требования к геометрии**: высокие требования к гладкости кривой интерфейса 4. **Отсутствие примеров практического применения**: недостаточно примеров применения к сложным инженерным задачам ### Влияние 1. **Академическая ценность**: предоставляет новую теоретическую схему и технологический путь для методов IFE 2. **Практическая ценность**: предоставляет эффективный инструмент для высокоточного решения задач на интерфейсах 3. **Воспроизводимость**: полная реализация кода облегчает использование и расширение другими исследователями 4. **Потенциал развития**: создаёт прочную основу для последующих исследований ### Области применения 1. **Многофазные течения**: численное моделирование интерфейсов различных жидкостей 2. **Материаловедение**: высокоточные расчёты задач на интерфейсах композитных материалов 3. **Биомедицина**: задачи теплопередачи и массопередачи на границах тканей 4. **Геофизика**: распространение волн на границах геологических слоёв ## Список литературы Статья ссылается на 22 работы из смежных областей, включая дифференциальную геометрию, численный анализ, методы конечных элементов и другие, обеспечивая прочную теоретическую базу для данного исследования.