2025-11-23T07:43:17.092257

Spectral analysis of the stiffness matrix sequence in the approximated Stokes equation

Ferri, Giraudo, Loi et al.

In the present paper, we analyze in detail the spectral features of the matrix sequences arising from the Taylor-Hood $\mathbb{P}_2$-$\mathbb{P}_1$ approximation of variable viscosity for $2d$ Stokes problem under weak assumptions on the regularity of the diffusion. Localization and distributional spectral results are provided, accompanied by numerical tests and visualizations. A preliminary study of the impact of our findings on the preconditioning problem is also presented. A final section with concluding remarks and open problems ends the current work.

academic

Спектральный анализ последовательности матрицы жёсткости в приближённом уравнении Стокса

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

ID статьи: 2510.25252
Название: Spectral analysis of the stiffness matrix sequence in the approximated Stokes equation
Авторы: Samuele Ferri, Chiara Giraudo, Valerio Loi, Miroslav Kuchta, Stefano Serra-Capizzano
Учреждения: Университет Инсубрии (Италия), Университет Осло (Норвегия), Лаборатория исследований Simula (Норвегия), Университет Упсалы (Швеция)
Классификация: math.NA (Численный анализ), cs.NA
Дата подачи: 29 октября 2025 г.
Ссылка на статью: https://arxiv.org/abs/2510.25252v1

Аннотация

В данной статье проводится детальный анализ спектральных характеристик последовательности матриц, возникающих при дискретизации двумерной задачи Стокса с переменной вязкостью методом конечных элементов Тейлора-Худа $\mathbb{P}_2$ - $\mathbb{P}_1$ . При наложении только слабых предположений о регулярности коэффициента диффузии статья предоставляет результаты локализации и распределения спектра с сопровождающими численными тестами и визуализацией. Кроме того, в статье проводится предварительное исследование влияния этих результатов на проблему предобусловливания, и в заключение предлагаются замечания и открытые вопросы.

Исследовательский контекст и мотивация

Контекст проблемы

Уравнение Стокса описывает течение вязкой несжимаемой жидкости без инерционных эффектов и широко применяется в:

Геодинамике (моделирование ледяных щитов)
Планетологии
Двухфазных потоках
Исследовании неньютоновских жидкостей
Задачах взаимодействия жидкости и твёрдого тела

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

Основная проблема

Дискретизация уравнения Стокса приводит к симметричной седловой системе: $\mathcal{A}_h \begin{bmatrix} u_h \\ p_h \end{bmatrix} = \begin{bmatrix} g_h \\ 0 \end{bmatrix}, \quad \text{где} \quad \mathcal{A}_h = \begin{bmatrix} A_h & B_h^T \\ B_h & 0 \end{bmatrix}$

Стандартный предобусловливатель имеет вид $\mathcal{B}_h = \text{diag}(A_h, M_h)^{-1}$ , где $M_h$ — матрица масс, взвешенная по $\mu^{-1}$ . Ключевые проблемы:

Производительность предобусловливателя сильно зависит от гладкости вязкости $\mu$
При наличии сильного контраста вязкости (например, $\mu_1/\mu_0 \sim 10^8$ ) число обусловленности резко возрастает
Существующая теория предоставляет только пессимистические границы числа обусловленности, не позволяя точно предсказать скорость сходимости решателя Крылова

Исследовательская мотивация

Данная статья имеет целью использовать теорию спектрального распределения Вейля и теорию обобщённых локальных матриц Тёплица (GLT):

Точно охарактеризовать спектральное распределение последовательностей матриц $\{A_n\}_n$ и $\{B_n\}_n$
Понять влияние переменной вязкости на спектральные характеристики
Обеспечить теоретическую основу для проектирования робастных предобусловливателей

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

Полная теория спектрального анализа: Впервые предоставляется строгий GLT спектральный символ для дискретизации методом Тейлора-Худа $\mathbb{P}_2$ - $\mathbb{P}_1$ задачи Стокса с переменной вязкостью
Технические инновации:
- Преобразование сложной блочной структуры в форму, пригодную для GLT, с использованием матриц перестановок и полуортогональных сжатий
- Решение технической проблемы нерегулярных размеров блоков матриц (через "extradimensional approach")
- Применение техники $L^1$ аппроксимации для разрывных функций вязкости
Явные формулы символов:
- Спектральный символ матрицы жёсткости $A_n$ : $f(x,y,\theta_1,\theta_2) = \mu(x,y)\hat{G}(\theta_1,\theta_2)$ (матричнозначная функция размера $16\times 16$ )
- Символ сингулярных значений матрицы дивергенции $B_n$ : $G(\theta_1,\theta_2)$ (матричнозначная функция размера $8\times 4$ )
Численная верификация:
- Визуализация подтверждает соответствие теоретического символа и фактического распределения собственных значений/сингулярных значений
- Подтверждено отсутствие выбросов собственных значений при экстремальном контрасте вязкости ( $10^8$ )
Приложение к предобусловливанию: Предложена предварительная стратегия предобусловливания на основе GLT, численные эксперименты показывают значительное снижение числа итераций PGMRES

Подробное описание методов

Определение задачи

Вход: Дискретизация методом конечных элементов Тейлора-Худа задачи Стокса с переменной вязкостью $-\nabla \cdot (\mu\nabla u) + \nabla p = g, \quad -\nabla \cdot u = 0, \quad u|_{\partial\Omega} = 0$

Выход: Спектральные символы последовательностей матриц $\{A_n\}_n$ , $\{B_n\}_n$ (функции, описывающие асимптотическое распределение собственных значений/сингулярных значений)

Ограничения: $\mu \in L^1(\Omega)$ , $\mu > 0$ почти всюду

Основная теоретическая база: Теория GLT

Определение GLT-последовательности

Последовательность матриц размера $d$ -слой $s$ -блок $\{A_n\}_n$ является GLT-последовательностью, обозначаемой $\{A_n\}_n \sim_{\text{GLT}} \kappa$ , если:

Существует единственная измеримая функция-символ $\kappa: [0,1]^d \times [-\pi,\pi]^d \to \mathbb{C}^{s\times s}$
Для любой компактно поддерживаемой непрерывной функции $F$ : $\lim_{n\to\infty} \frac{1}{d_n}\sum_{j=1}^{d_n} F(\lambda_j(A_n)) = \frac{1}{\mu_d(D)}\int_D \frac{1}{s}\sum_{i=1}^s F(\lambda_i(\kappa(t)))dt$

Алгебраические свойства GLT (ключевые аксиомы)

GLT 1-4: Последовательности Тёплица, последовательности диагональной выборки, последовательности нулевого распределения являются GLT-последовательностями
GLT 5: GLT-последовательности образуют $*$ -алгебру (замкнуты относительно сложения, умножения, сопряжённого транспонирования)
GLT 6: Аппроксимирующие классы последовательностей (a.c.s.) сохраняют свойство GLT

Техническая схема

Первый шаг: Анализ блочной структуры матрицы жёсткости $A_n$

Начальная структура: $A_n = \begin{bmatrix} A_{x,n} & 0 \\ 0 & A_{y,n} \end{bmatrix}, \quad A_{x,n} = A_{y,n}$

$A_{x,n}$ имеет периодическую блочную структуру размера $(8n^2 - O(n)) \times (8n^2 - O(n))$ : $\hat{A}_{x,n} = \text{tridiag}(A_1, A_0, A_1^T)$ где $A_0$ и $A_1$ — блочные матрицы размера $8n \times 8n$ .

Технические трудности:

Нерегулярные размеры блоков ( $2n \times 2n$ , $(2n-2) \times 2n$ и т.д. в смешанном виде)
Наличие возмущений ранга $O(n)$ на границе

Решение (Extradimensional Approach):

Расширение: Добавление строк/столбцов размера $O(n)$ для унификации всех блоков в размер $2n \times 2n$ , получение $\tilde{A}_{x,n}$ (размер $8n^2 \times 8n^2$ )
Сжатие: Через полуортогональную матрицу $P_n$ такую, что $P_n^* P_n = I$ , получаем $A_n = P_n^* \tilde{A}_n P_n + E_n$ где $\{E_n\}_n$ — последовательность нулевого распределения ( $\text{rank}(E_n) = O(n)$ )

Теоретическое обоснование (Теорема 3.11): Если $\lim_{n\to\infty} n'/n = 1$ , то $\{X_n\}_n \sim_\lambda f \Leftrightarrow \{Y_{n'}\}_{n'} \sim_\lambda f, \quad Y_{n'} = P_n^* X_n P_n$

Второй шаг: Преобразование подобия перестановкой

Применение матрицы перестановки $\Pi_{4,n} = I_n \otimes \Pi_{2n,4,2}$ даёт: $\Pi_{4,n}\hat{A}_{x,n}\Pi_{4,n}^* = \text{tridiag}(T_n(g_1), T_n(g_0), T_n(g_1)^T) \odot M$

где:

$T_n(g_i)$ — матрица Тёплица, порождённая символом $g_i(\theta_2)$ (блок размера $8\times 8$ )
$M$ — матрица выборки вязкости
$\odot$ обозначает произведение Адамара

Ключевые символы (матричнозначные функции размера $8\times 8$ ):