2025-11-24T06:52:18.097833

Optimal $L^2$ error estimation for the unfitted interface finite element method based on the non-symmetric Nitsche's methods

Chen, Liu, Zhang
This paper establishes optimal error estimates in the $L^2$ for the non-symmetric Nitsche method in an unfitted interface finite element setting. Extending our earlier work, we give a complete analysis for the Poisson interface model and, by formulating a tailored dual problem that restores adjoint consistency, derive the desired bounds.
academic

Оптимальная оценка ошибки L2L^2 для метода неподогнанных интерфейсных конечных элементов на основе несимметричных методов Нитше

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

  • ID статьи: 2510.12151
  • Название: Optimal L2L^2 error estimation for the unfitted interface finite element method based on the non-symmetric Nitsche's methods
  • Авторы: Gang Chen (Сычуаньский университет), Chaoran Liu (Сычуаньский университет), Yangwen Zhang (Университет Луизианы в Лафайете)
  • Классификация: math.NA cs.NA
  • Дата публикации: 15 октября 2025 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2510.12151

Аннотация

В данной работе установлены оптимальные оценки ошибки L2L^2 для несимметричного метода Нитше в контексте неподогнанного интерфейсного метода конечных элементов. Путём расширения более ранних работ авторов проведён полный анализ модели интерфейса Пуассона и получены желаемые оценки посредством конструирования специализированной двойственной задачи, восстанавливающей согласованность сопряжённого оператора.

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

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

  1. Значимость интерфейсных задач: Интерфейсные задачи имеют широкое применение в инженерии и науке, включая многофизические связи, многофазные потоки и механику композитных материалов. Дифференциальные уравнения в частных производных для этих задач связаны через интерфейс, который разделяет область на подобласти с различными материальными свойствами.
  2. Сложность численного решения: Традиционные подогнанные методы конечных элементов требуют строгого соответствия сетки геометрии интерфейса при работе со сложными или эволюционирующими интерфейсами, что приводит к значительному увеличению вычислительных затрат.
  3. Преимущества неподогнанных методов: Неподогнанные методы конечных элементов не требуют согласования сетки с интерфейсом, решая задачу путём разрезания элементов сетки, пересекающих интерфейс, и применения специализированных стратегий обработки интерфейсных условий.

Ограничения существующих методов

  1. Симметричный метод Нитше: Требует достаточно большого параметра стабилизации, чувствителен к выбору параметра
  2. Несимметричный метод Нитше: Хотя устраняет строгие требования к параметру штрафа, теоретический анализ даёт только субоптимальные оценки ошибки L2L^2, несмотря на то, что численные результаты указывают на возможность оптимальности

Теоретические препятствия

Введение несимметричных членов дополнительно нарушает согласованность сопряжённого оператора вариационной формы, что затрудняет прямое применение классического двойственного аргумента (техника Обена-Нитше) для оценок в норме L2L^2.

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

  1. Построение улучшенной двойственной схемы аргументации: Применимой к несимметричным вариационным формам с объединением специализированного анализа регулярности вблизи интерфейса для преодоления теоретических препятствий, вызванных несимметричностью
  2. Строгое доказательство оптимальной оценки ошибки в норме L2L^2: Для несимметричного метода Нитше классической интерфейсной задачи
  3. Восстановление согласованности сопряжённого оператора: Путём конструирования специализированной двойственной задачи, решающей ключевую проблему нарушения согласованности сопряжённого оператора несимметричностью
  4. Полный анализ регулярности: Установление комплексной теории регулярности для интерфейсной задачи и двойственной задачи

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

Постановка задачи

Рассмотрим эллиптическую интерфейсную задачу с погруженным интерфейсом Γ\Gamma в выпуклой многоугольной/многогранной области Ω⊂Rd\Omega \subset \mathbb{R}^d (d=2,3d = 2, 3):

{−∇⋅(μ∇u)=fв Ω[[μ∇u⋅n]]=gNна Γ[[u]]=gDна Γu=0на ∂Ω\begin{cases} -\nabla \cdot (\mu\nabla u) = f & \text{в } \Omega \\ [[\mu\nabla u \cdot n]] = g_N & \text{на } \Gamma \\ [[u]] = g_D & \text{на } \Gamma \\ u = 0 & \text{на } \partial\Omega \end{cases}

где μ=μi\mu = \mu_i является кусочно-постоянной на Ωi\Omega_i (i=1,2i = 1, 2), [[u]]=u1−u2[[u]] = u_1 - u_2 обозначает скачок, [[μ∇u⋅n]]=μ1∇u1⋅n−μ2∇u2⋅n[[\mu\nabla u \cdot n]] = \mu_1\nabla u_1 \cdot n - \mu_2\nabla u_2 \cdot n обозначает скачок потока.

Несимметричный метод Нитше

Определим билинейную форму:

a(u,v)=∑i=12μi(∇ui,∇vi)Ωi−⟨{{μ∇u⋅n}},[[v]]⟩Γ+⟨[[u]],{{μ∇v⋅n}}⟩Γ+⟨c0h−1[[u]],[[v]]⟩Γa(u,v) = \sum_{i=1}^2 \mu_i(\nabla u_i, \nabla v_i)_{\Omega_i} - \langle\{\{\mu\nabla u \cdot n\}\}, [[v]]\rangle_\Gamma + \langle[[u]], \{\{\mu\nabla v \cdot n\}\}\rangle_\Gamma + \langle c_0h^{-1}[[u]], [[v]]\rangle_\Gamma

где взвешенное среднее определяется как:

{{v}}=w1v1+w2v2,w1=μ2μ1+μ2,w2=μ1μ1+μ2\{\{v\}\} = w_1v_1 + w_2v_2, \quad w_1 = \frac{\mu_2}{\mu_1 + \mu_2}, \quad w_2 = \frac{\mu_1}{\mu_1 + \mu_2}

Параметр стабилизации выбирается как c0={{μ}}=w1μ1+w2μ2c_0 = \{\{\mu\}\} = w_1\mu_1 + w_2\mu_2.

Ключевые технические инновации

1. Переконструирование двойственной задачи

Традиционная двойственная задача: a(ψ,v)=(g,v)a(\psi, v) = (g, v)

Двойственная задача, построенная в данной работе: a∗(u,v)=(f,v)a^*(u, v) = (f, v), где

a∗(u,v)=∑i=12μi(∇ui,∇vi)Ωi+⟨{{μ∇u⋅n}},[[v]]⟩Γ−⟨[[u]],{{μ∇v⋅n}}⟩Γ+⟨c0h−1[[u]],[[v]]⟩Γa^*(u,v) = \sum_{i=1}^2 \mu_i(\nabla u_i, \nabla v_i)_{\Omega_i} + \langle\{\{\mu\nabla u \cdot n\}\}, [[v]]\rangle_\Gamma - \langle[[u]], \{\{\mu\nabla v \cdot n\}\}\rangle_\Gamma + \langle c_0h^{-1}[[u]], [[v]]\rangle_\Gamma

удовлетворяющая соотношению: a(v,u)=a∗(u,v)a(v, u) = a^*(u, v), тем самым восстанавливая согласованность сопряжённого оператора.

2. Схема анализа регулярности

Установлена теория регулярности для задач с граничными условиями Робина на интерфейсе:

Теорема: Для задачи с граничными условиями Робина на интерфейсе

{−∇⋅(μ∇u)=fв Ωγ[[u]]+μ1∇u1⋅n=gDна Γ[[μ∇u⋅n]]=gNна Γ\begin{cases} -\nabla \cdot (\mu\nabla u) = f & \text{в } \Omega \\ \gamma[[u]] + \mu_1\nabla u_1 \cdot n = g_D & \text{на } \Gamma \\ [[\mu\nabla u \cdot n]] = g_N & \text{на } \Gamma \end{cases}

когда Ω\Omega является выпуклой областью и γ≥γ0>0\gamma \geq \gamma_0 > 0, справедливо

μ1∥u∥2,Ω1+μ2∥u∥2,Ω2≤C(1+γ0−1)max⁡{μ1−1,μ2−1,1}(1+(μ1/μ2)1/2)(∥f∥0+∥gN∥1/2,Γ+∥gD∥1/2,Γ)\mu_1\|u\|_{2,\Omega_1} + \mu_2\|u\|_{2,\Omega_2} \leq C(1 + \gamma_0^{-1})\max\{\mu_1^{-1}, \mu_2^{-1}, 1\}(1 + (\mu_1/\mu_2)^{1/2})(\|f\|_0 + \|g_N\|_{1/2,\Gamma} + \|g_D\|_{1/2,\Gamma})

Экспериментальная установка

Теоретическая верификация

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

  1. Оценка ошибки в энергетической норме:
∥u−uh∥∗≤Chk(μ11/2∣u1∣k+1,Ω1+μ21/2∣u2∣k+1,Ω2)\|u - u_h\|_{*} \leq Ch^k(\mu_1^{1/2}|u_1|_{k+1,\Omega_1} + \mu_2^{1/2}|u_2|_{k+1,\Omega_2})
  1. Оценка ошибки в норме L2L^2:
∥u−uh∥0≤Chk+1max⁡{μ1−1/2,μ1−3/2,μ2−1/2,μ2−3/2}(μ11/2∣u1∣k+1,Ω1+μ21/2∣u2∣k+1,Ω2)\|u - u_h\|_0 \leq Ch^{k+1}\max\{\mu_1^{-1/2}, \mu_1^{-3/2}, \mu_2^{-1/2}, \mu_2^{-3/2}\}(\mu_1^{1/2}|u_1|_{k+1,\Omega_1} + \mu_2^{1/2}|u_2|_{k+1,\Omega_2})

В частности, когда μ1,μ2≥1\mu_1, \mu_2 \geq 1:

∥u−uh∥0≤Chk+1max⁡{μ1−1/2,μ2−1/2}(μ11/2∣u1∣k+1,Ω1+μ21/2∣u2∣k+1,Ω2)\|u - u_h\|_0 \leq Ch^{k+1}\max\{\mu_1^{-1/2}, \mu_2^{-1/2}\}(\mu_1^{1/2}|u_1|_{k+1,\Omega_1} + \mu_2^{1/2}|u_2|_{k+1,\Omega_2})

Расширение на методы без штрафа

Статья также анализирует случай без штрафа посредством члена фиктивного штрафа (ghost penalty):

sh(uh,vh)=γg∑i=12∑E∈Ehi∑l=1k⟨μih2l−1[[DnEluih]],[[DnElvih]]⟩Es_h(u_h, v_h) = \gamma_g \sum_{i=1}^2 \sum_{E \in \mathcal{E}_h^i} \sum_{l=1}^k \langle\mu_ih^{2l-1}[[D^l_{n_E}u_{ih}]], [[D^l_{n_E}v_{ih}]]\rangle_E

получив оптимальные оценки ошибки того же порядка.

Результаты экспериментов

Основные теоретические результаты

  1. Оптимальный порядок сходимости: Достигнут оптимальный порядок сходимости O(hk+1)O(h^{k+1}) в норме L2L^2, что впервые реализовано для несимметричного метода Нитше
  2. Оценки, независимые от коэффициентов: Зависимость константы ошибки от коэффициента скачка μ1/μ2\mu_1/\mu_2 получена точно
  3. Доказательство стабильности: Установлено условие inf-sup для несимметричной билинейной формы

Технические прорывы

  1. Восстановление согласованности сопряжённого оператора: Путём переконструирования двойственной задачи устранены дополнительные граничные члены, вызванные несимметричными членами
  2. Совершенствование теории регулярности: Установлены оценки регулярности, независимые от параметров, для задач с граничными условиями Робина на интерфейсе

Связанные работы

Линия развития

  1. 2002 год: Hansbo и Hansbo введли двумерный неподогнанный метод конечных элементов на основе симметричной формы Нитше
  2. 2016-2018 годы: Burman развил стабилизированные методы в рамках симметричного метода Нитше, установив оптимальные оценки ошибки в норме L2L^2
  3. 2017 год: Huang и др. введли несимметричную форму Нитше, устраняющую строгие требования к параметру штрафа, но получившую только субоптимальные оценки L2L^2
  4. 2017 год: Boiveau ввёл метод без штрафа с несимметричной формой Нитше

Уникальность вклада данной работы

Данная работа впервые теоретически доказывает оптимальность несимметричного метода Нитше в норме L2L^2, заполняя разрыв между теоретическими результатами и численными результатами.

Выводы и обсуждение

Основные выводы

  1. Успешно установлены оптимальные оценки ошибки L2L^2 для несимметричного метода Нитше
  2. Путём конструирования специальной двойственной задачи восстановлена согласованность сопряжённого оператора
  3. Разработана полная схема теории регулярности

Ограничения

  1. Анализ ограничен интерфейсными задачами типа Пуассона
  2. Требуется, чтобы интерфейс был C2C^2 гладким
  3. Область должна быть выпуклой многоугольной/многогранной

Направления будущих исследований

  1. Расширение на более общие эллиптические операторы
  2. Рассмотрение случаев негладких интерфейсов
  3. Разработка стратегий адаптивного уточнения сетки

Глубокая оценка

Преимущества

  1. Теоретический прорыв: Решение долгосуществующего теоретического разрыва в несимметричном методе Нитше
  2. Техническая инновация: Переконструирование двойственной задачи обладает высокой степенью инновационности
  3. Строгость анализа: Анализ регулярности является всеобъемлющим и глубоким, доказательства строги
  4. Практическая ценность: Обеспечение прочной теоретической основы для неподогнанных методов

Недостатки

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

Влияние

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

Области применения

  1. Численное моделирование многофазных потоков жидкости
  2. Механический анализ композитных материалов
  3. Интерфейсные задачи теплопередачи и массопередачи
  4. Многорегиональное моделирование в биомедицине

Библиография

Статья цитирует 37 связанных работ, охватывающих важные работы в ключевых областях методов конечных элементов, интерфейсных задач и методов Нитше, обеспечивая прочную теоретическую основу для исследования.

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