2025-11-22T09:58:16.048981
Superconvergent and Divergence-Free Finite Element Methods for Stokes Equation
Chen, Huang, Zhang et al.
Superconvergent and divergence-free finite element methods for the Stokes equation are developed. The velocity and pressure are discretized using $H(\mathrm{div})$-conforming vector elements and discontinuous piecewise polynomials. The discrete formulation employs a weak deviatoric gradient operator built with tangential-normal continuous finite elements for traceless tensors, requiring no stabilization. Optimal and superconvergent error estimates are established. The method connects to nonconforming virtual element and pseudostress-velocity-pressure mixed formulations. Numerical experiments verify the theory.
academicСуперсходящиеся и дивергентно-свободные методы конечных элементов для уравнения Стокса
Основная информация
- ID статьи: 2510.14192
- Название: Superconvergent and Divergence-Free Finite Element Methods for Stokes Equation
- Авторы: Long Chen, Xuehai Huang, Chao Zhang, Xinyue Zhao
- Классификация: math.NA cs.NA
- Дата подачи: 16 октября 2025 г. на arXiv
- Ссылка на статью: https://arxiv.org/abs/2510.14192
Аннотация
В данной работе разработаны суперсходящиеся и дивергентно-свободные методы конечных элементов для решения уравнения Стокса. Поле скорости дискретизируется с использованием H(div)-совместимых векторных элементов, а поле давления — с использованием разрывных кусочно-полиномиальных функций. Дискретная схема использует оператор слабого девиаторного градиента, построенный на основе бесследовых тензоров с касательно-нормальной непрерывностью конечных элементов, без необходимости стабилизации. Установлены оптимальные и суперсходящиеся оценки погрешности. Метод связан с несогласованными виртуальными элементами и смешанной схемой псевдонапряжение-скорость-давление. Численные эксперименты подтверждают теоретические результаты.
Исследовательский контекст и мотивация
Важность проблемы
Уравнение Стокса является фундаментальным уравнением, описывающим движение несжимаемой жидкости, и занимает центральное место в вычислительной гидродинамике. Численное решение этого уравнения сталкивается с двумя основными вызовами:
- Ограничение дивергентной свободы: поле скорости должно удовлетворять div u = 0, что отражает сохранение массы
- Робастность давления: численный метод должен быть робастным при обработке членов давления
Ограничения существующих методов
Классические стабилизированные пары конечных элементов (такие как элементы Тейлора-Худа, элементы MINI, несогласованные элементы P1-P0) имеют следующие проблемы:
- Отсутствие поточечной дивергентной свободы: не могут строго удовлетворять ограничению дивергентной свободы, что ослабляет сохранение массы
- Недостаточная робастность давления: в оценках погрешности отсутствует робастность давления
- Ограничение порядка сходимости: некоторые методы страдают от потери порядка сходимости
Исследовательская мотивация
Целью данной работы является построение метода конечных элементов, который одновременно сохраняет строгое свойство дивергентной свободы и достигает суперсходимости, избегая при этом стабилизации, требуемой традиционными методами.
Основные вклады
- Предложена новая смешанная схема конечных элементов: использует H(div)-совместимые элементы для дискретизации скорости и разрывные кусочно-полиномиальные функции для давления
- Построен оператор слабого девиаторного градиента: на основе бесследовых тензоров с касательно-нормальной непрерывностью конечных элементов, без стабилизации
- Установлены оценки суперсходимости: доказана суперсходимость порядка
- Предоставлены методы постобработки: построены суперсходящиеся приближения скорости
- Установлены эквивалентные соотношения: связь с несогласованными виртуальными элементами и схемой псевдонапряжений
Описание метода
Постановка задачи
Рассмотрим уравнение Стокса на ограниченной области :
-\Delta u - \nabla p = f & \text{в } \Omega \\ \text{div } u = 0 & \text{в } \Omega \\ u = 0 & \text{на } \partial\Omega \end{cases}$$ где $u$ — поле скорости, $p$ — поле давления, $f$ — внешняя сила. ### Построение пространств конечных элементов #### Пространство скорости Используются H(div)-совместимые элементы Равиара-Томаса (RT) или Бреззи-Дугласа-Марини (BDM): $$\stackrel{\circ}{V}^{\text{div}}_{k,\ell} := \{v_h \in H_0(\text{div},\Omega) : v_h|_T \in P_k(T;\mathbb{R}^d) + H_\ell(T)x\}$$ где $\ell = k$ (элементы RT) или $\ell = k-1$ (элементы BDM). #### Пространство напряжений Введено пространство бесследовых тензоров $\Sigma^{tn}_k$ для дискретизации девиаторного градиента: $$\Sigma^{tn}_k := \{\tau_h \in \Sigma^{-1}_k(\mathcal{T}) : [\Pi_F \tau n]_F = 0 \text{ для всех } F \in \mathcal{F}_h\}$$ #### Пространство множителей Лагранжа Для ослабления касательно-нормальной непрерывности введено: $$\Lambda_k = P_k(\mathcal{F}_h;\mathbb{R}^{d-1}), \quad \stackrel{\circ}{\Lambda}_k = P_k(\stackrel{\circ}{\mathcal{F}}_h;\mathbb{R}^{d-1})$$ ### Оператор слабого девиаторного градиента Определен оператор слабого девиаторного градиента $\text{dev grad}_w : \stackrel{\circ}{V}^{\text{div}}_{k,\ell} \times \stackrel{\circ}{\Lambda}_k \to \Sigma^{-1}_k(\mathcal{T})$: Для $(v,\mu) \in H^1(\mathcal{T}_h;\mathbb{R}^d) \times L^2(\mathcal{F}_h;\mathbb{R}^{d-1})$, определяется поэлементно: $$(\text{dev grad}_w(v,\mu), \tau)_T = -(v, \text{div } \tau)_T + (n \cdot v, n^\top \tau n)_{\partial T} + (\mu, \Pi_F \tau n)_{\partial T}$$ ### Смешанная схема конечных элементов Найти $u_h \in \stackrel{\circ}{V}^{\text{div}}_{k,\ell}$, $\lambda_h \in \stackrel{\circ}{\Lambda}_k$, $p_h \in P_\ell(\mathcal{T}_h)/\mathbb{R}$ такие, что: $$(\text{dev grad}_w(u_h,\lambda_h), \text{dev grad}_w(v_h,\mu_h)) + (\text{div } v_h, p_h) = (f, v_h)$$ $$(\text{div } u_h, q_h) = 0$$ для всех $v_h \in \stackrel{\circ}{V}^{\text{div}}_{k,\ell}$, $\mu_h \in \stackrel{\circ}{\Lambda}_k$, $q_h \in P_\ell(\mathcal{T}_h)/\mathbb{R}$. ### Технические инновации 1. **Конструкция без стабилизации**: тщательно построенный слабый оператор избегает штрафных членов, требуемых традиционными методами DG 2. **Суперсходимость**: использование свойств коммутативности достигает сходимости на один порядок выше стандартных методов 3. **Робастность давления**: метод естественным образом обладает робастностью давления 4. **Доступность низких порядков**: поддерживает низкие порядки $(k,\ell) = (0,-1), (0,0), (1,0)$ и т.д. ## Экспериментальная установка ### Численные примеры **Пример 5.1 (2D)**: - Точное решение: $u = \text{curl } \psi_2$, $p = -x^5 - y^5 + \frac{1}{3}$ - где $\psi_2 = x^2(x-1)^2y^2(y-1)^2$ **Пример 5.2 (3D)**: - Точное решение: $u = \text{curl}(\psi_3, \psi_3, \psi_3)^T$, $p = -x^5 - y^5 - z^5 + \frac{1}{2}$ - где $\psi_3 = x^2(x-1)^2y^2(y-1)^2z^2(z-1)^2$ ### Вычислительная область и сетка - Вычислительная область: $\Omega = (0,1)^d$, $d = 2,3$ - Сетка: равномерное разбиение на симплексы - Реализация: на основе пакета MATLAB iFEM ### Показатели оценки - Погрешность скорости: $\|u - u_h\|$ - Погрешность напряжения: $\|\sigma - \sigma_h\|$ - Погрешность давления: $\|p - p_h\|$ - Погрешность постобработки: $\|u - u_h^*\|$, $\|\text{grad}_h(u - u_h^*)\|$ ## Результаты экспериментов ### Основные результаты **Таблица 1: Результаты погрешности в 2D** | h | (k,ℓ) | $\|u-u_h\|$ | порядок | $\|\sigma-\sigma_h\|$ | порядок | $\|p-p_h\|$ | порядок | |---|-------|-------------|---------|-------------------|---------|-------------|---------| | 2⁻³ | (0,0) | 2.988e-03 | - | 3.103e-02 | - | 7.810e-02 | - | | 2⁻⁴ | (0,0) | 1.284e-03 | 1.22 | 1.677e-02 | 0.89 | 3.914e-02 | 1.00 | | 2⁻⁵ | (0,0) | 5.988e-04 | 1.10 | 8.700e-03 | 0.95 | 1.963e-02 | 1.00 | | 2⁻³ | (1,0) | 3.296e-04 | - | 2.447e-03 | - | 7.453e-02 | - | | 2⁻⁴ | (1,0) | 8.382e-05 | 1.98 | 6.305e-04 | 1.96 | 3.760e-02 | 0.99 | | 2⁻⁵ | (1,0) | 2.104e-05 | 1.99 | 1.597e-04 | 1.98 | 1.880e-02 | 1.00 | ### Проверка порядка сходимости Экспериментальные результаты полностью подтверждают теоретические предсказания: - **Скорость и напряжение**: $\|u - u_h\| = \|\sigma - \sigma_h\| = O(h^{k+1})$ (суперсходимость) - **Давление**: $\|p - p_h\| = O(h)$ - **Скорость после постобработки**: $\|u - u_h^*\| = O(h^{k+2})$, $\|\text{grad}_h(u - u_h^*)\| = O(h^{k+1})$ ### Результаты в 3D Эксперименты в 3D также подтверждают эффективность метода, порядки сходимости соответствуют теории. ## Теоретический анализ ### Анализ устойчивости Установлено условие слабой дивергентной устойчивости: $$\inf_{v_h \in \stackrel{\circ}{V}^{\text{div}}_{k,k-1}} \sup_{\tau_h \in \Sigma^{tn}_k} \frac{(\text{div}_w \tau_h, v_h)_{0,h}}{\|\tau_h\|_{\text{div}_w} \|v_h\|} = \alpha > 0$$ ### Оценки погрешности **Теорема**: Предположим, что $u \in H^{k+2}(\Omega;\mathbb{R}^d)$, тогда: $$\|\sigma - \sigma_h\|_{0,h} + \|\text{dev grad}_w(I^{\text{div}}_{k,k}u - u_h, Q_{k,\mathcal{F}_h}\lambda - \lambda_h)\| + \|Q_\ell p - p_h\| \lesssim h^{k+1}|u|_{k+2}$$ ### Свойство коммутативности Ключевое свойство коммутативности: $$Q^{tn}_k \text{dev grad} = \text{dev grad}_w I^{\text{div}}_{k,k}$$ Это является основой для достижения суперсходимости. ## Связанные работы ### Методы дивергентно-свободных конечных элементов - **Элементы Скотта-Фогелиуса**: требуют специальных условий на сетку - **Гладкие пары конечных элементов**: требуют сверхгладких степеней свободы - **Согласованные пары на разделенных сетках**: сложная реализация ### Смешанные методы - **Метод MCS**: требует $k = \ell \geq 1$, данная работа охватывает больше низкопорядковых случаев - **Метод HDG**: требует стабилизации, данная работа без стабилизации - **Виртуальные элементы**: обычно требуют стабилизации, более низкие порядки сходимости ## Заключение и обсуждение ### Основные выводы 1. Успешно построен метод суперсходящихся дивергентно-свободных конечных элементов без стабилизации 2. Достигнута суперсходимость порядка $h^{k+1}$, не зависящая от симметрии сетки 3. Метод обладает робастностью давления и поддерживает низкие порядки 4. Установлена эквивалентность с виртуальными элементами и схемой псевдонапряжений ### Ограничения 1. Теоретический анализ в основном ориентирован на многогранные области 2. Предположение о регулярности $H^2$ ограничивает область применения 3. Сложность реализации относительно высока ### Будущие направления 1. Расширение на уравнения Навье-Стокса 2. Адаптивное уточнение сетки 3. Разработка параллельных алгоритмов 4. Верификация на инженерных приложениях ## Глубокая оценка ### Преимущества 1. **Теоретическая строгость**: полный анализ устойчивости и сходимости 2. **Методологическая инновация**: искусное построение оператора слабого девиаторного градиента 3. **Практическая ценность**: поддержка низких порядков, достижение суперсходимости 4. **Полная экспериментальная верификация**: численные проверки в 2D и 3D ### Недостатки 1. **Сложность реализации**: вычисление слабого оператора относительно сложно 2. **Теоретические ограничения**: требуются сильные предположения о регулярности 3. **Отсутствие верификации на приложениях**: нет проверки на реальных инженерных задачах ### Влияние Данная работа имеет важное теоретическое значение в области численных методов для уравнения Стокса, предоставляя новые идеи для построения высокоточных дивергентно-свободных методов. Свойства суперсходимости и отсутствие необходимости в стабилизации делают метод потенциально полезным для вычислительной гидродинамики. ### Области применения - Расчеты жидкостей, требующие точного сохранения массы - Приложения, предъявляющие высокие требования к точности вычислений - Академические исследования и верификация методов ## Библиография Статья цитирует 56 связанных работ, охватывающих методы конечных элементов, смешанные методы, виртуальные элементы и другие важные направления, обеспечивая полный обзор литературы.