Острая стабильность спектра подъёма-спуска при сильной сходимости резольвенты

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

• ID статьи: 2511.20971
• Название: Sharp Ascent--Descent Spectral Stability under Strong Resolvent Convergence
• Автор: Marwa Ennaceur
• Классификация: math.NA (Численный анализ), cs.NA, math.FA (Функциональный анализ), math.SP (Спектральная теория)
• Дата подачи: 26 ноября 2025 г. на arXiv
• Ссылка на статью: https://arxiv.org/abs/2511.20971

Аннотация

В статье установлены острые результаты стабильности спектров подъёма и спуска (ascent and descent spectra) несамосопряжённых операторов в рамках сильной сходимости резольвенты (Strong Resolvent Convergence, SRS). Ключевым количественным условием является приведённый минимальный модуль (reduced minimum modulus) $\gamma(T-\lambda)>0$, который гарантирует замкнутость области значений и обеспечивает передачу критерия Каасхука-Тейлора через сходимость графиков операторов. На фундаментальном уровне теория B-Фредгольма расширяет стабильность на степени $(T-\lambda)^m$ при условии, что все промежуточные степени удовлетворяют $\gamma((T-\lambda)^j)>0$ ($1\le j\le m$). Авторы вводят вычислимую диагностическую величину конечных элементов $\gamma_h = \sigma_{\min}(M^{-1/2}(A_h-\lambda M)M^{-1/2})$ как практическую замену $\gamma(T-\lambda)$, сохраняющую последовательную положительность даже в конвективно-доминирующих случаях при использовании стабилизированных схем (таких как SUPG). Численные эксперименты подтверждают, что $\liminf_{h\to0}\gamma_h>0$ является необходимым и достаточным условием спектральной стабильности, тогда как вольтеррины контрпримеры демонстрируют неустранимость условия замкнутости области значений для степеней.

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

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

Статья исследует стабильность тонких спектральных инвариантов (индексов подъёма и спуска) дифференциальных операторов (особенно несамосопряжённых и сингулярно возмущённых) при дискретизации методом конечных элементов. Конкретно:

• Индекс подъёма $\text{asc}(S)$: минимальный порядок, на котором стабилизируется цепь ядер $\ker S \subset \ker S^2 \subset \cdots$
• Индекс спуска $\text{dsc}(S)$: минимальный порядок, на котором стабилизируется цепь областей значений $\text{Ran}(S) \supset \text{Ran}(S^2) \supset \cdots$

Эти инварианты чрезвычайно хрупки в конечномерном случае: нильпотентная матрица $S$ может иметь индекс подъёма 3, но произвольно малое возмущение $S+\varepsilon I$ приводит его к нулю.

2. Важность проблемы

Классический пример (парадокс из введения):

• Непрерывный одномерный оператор переноса $Lu=u'$ ($u(0)=0$) на $L^2(0,1)$ имеет замкнутую область значений и $\text{asc}(L)=1$
• При дискретизации центральной разностью матрица $A_h$ имеет собственные значения, скапливающиеся около нуля (предполагая сходимость спектра), но $\text{asc}(A_h)=\infty$ (катастрофический отказ)
• Схема против потока (upwind) корректно сохраняет $\text{asc}(A_h)=1$ и дискретную замкнутость области значений

Коренная причина: для центральной разности $\gamma_h \sim Ch \to 0$, тогда как для схемы против потока $\gamma_h \ge c > 0$.

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

• Сходимость резольвенты по норме (norm resolvent convergence): отказывает при грубых коэффициентах ($L^\infty$ потенциалы) или сингулярных пределах (например, $\varepsilon \to 0$ в конвективной диффузии)
• Сходимость Мосько: применима только к секториальным операторам, не может обрабатывать чистые конвективные задачи
• Отсутствие вычислимых диагностических величин для предсказания успеха/отказа дискретизации

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

Предоставить единую теоретическую базу:

• Работать в рамках сильной сходимости резольвенты (SRS), которая является естественной и широко применимой
• Связать количественное условие $\gamma>0$ с теорией непрерывных операторов и дискретным численным анализом
• Ввести практическую диагностическую величину $\gamma_h$ для руководства выбором стабилизированных схем

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

1. Острая теорема стабильности (Теорема 3.2): при SRS доказана непрерывность и замкнутость спектров подъёма/спуска с необходимым и достаточным условием $\gamma(S)>0$ или $\limsup_n \gamma(S_n)>0$
2. Механизм передачи через степени (Лемма 3.4): раскрыта необходимость сходимости графиков $G(S_n^m) \xrightarrow{\text{gap}} G(S^m)$ для всех промежуточных степеней $\gamma(S^j)>0$ ($1\le j\le m$)
3. Вычислимая диагностическая величина: $\gamma_h = \sigma_{\min}(M^{-1/2}(A_h-\lambda M)M^{-1/2})$ как дискретный аналог непрерывного условия $\gamma(T-\lambda)>0$, эффективно вычисляемая стандартными инструментами линейной алгебры
4. Три механизма стабильности (Теоремы 3.15-3.16):
   • Секториальный случай: сходимость Мосько $\Rightarrow$ SRS $\Rightarrow$ $\gamma_h \ge \gamma(T-\lambda)>0$
   • SUPG стабилизация: установлена равномерная условие inf-sup $\Rightarrow$ $\gamma_h^{\text{stab}} \ge c>0$ (даже при $\varepsilon \to 0$)
   • Чистый конвективный предел: предложена гипотеза 3.17 с численными доказательствами
5. Контрпримеры и необходимость: оператор Вольтерра (Приложение A) и отказ центральной разности (Предложение 4.12) демонстрируют неустранимость условия замкнутости области значений
6. Алгоритм адаптивного уточнения сетки (Алгоритм 2): стратегия маркировки элементов на основе локальной диагностической величины $\gamma_h(K)$

Детальное описание методов

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

Дан непрерывный оператор $T$ и его последовательность конечноэлементных приближений $T_h$. Требуется определить, когда дискретные спектры подъёма/спуска сходятся к непрерывным: $\sigma_{\text{asc}}(T_h) \to \sigma_{\text{asc}}(T), \quad \sigma_{\text{dsc}}(T_h) \to \sigma_{\text{dsc}}(T)$

Входные данные:

• Слабая форма непрерывного оператора $T$: $a(\cdot,\cdot)$
• Конечноэлементное пространство $V_h \subset H^1_0(\Omega)$
• Спектральный параметр $\lambda \in \mathbb{C}$

Выходные данные:

• Диагностическая величина стабильности $\gamma_h$
• Проверка сохранения индексов подъёма/спуска

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

1. Приведённый минимальный модуль и замкнутость области значений

Для замкнутого оператора $S$ определяется: $\gamma(S) := \inf\{\|Sx\| : x \in D(S), \, \text{dist}(x, \ker S) = 1\}$

Фундаментальная эквивалентность (Като 21): $\gamma(S) > 0 \iff \text{Ran}(S) \text{ замкнута}$

2. Критерий Каасхука-Тейлора

$\text{asc}(S) < \infty \iff \exists m : \text{Ran}(S^m) \cap \ker S = \{0\}$$\text{dsc}(S) < \infty \iff \exists m : \text{Ran}(S) + \ker(S^m) = H$

Эти условия трансверсальности подпространств стабильны при сходимости графиков, но при условии, что $\gamma(S)>0$ гарантирует замкнутость области значений.

3. Сильная сходимость резольвенты и сходимость графиков

Определение SRS: для $\mu \in \rho(T) \cap \bigcap_n \rho(T_n)$, $(T_n - \mu)^{-1}x \to (T - \mu)^{-1}x \quad \forall x \in H$

Ключевая лемма (Лемма 3.3): сходимость графиков $G(S_n) \xrightarrow{\text{gap}} G(S)$ (т.е. $\|P_{G(S_n)} - P_{G(S)}\| \to 0$) сохраняет:

• Сходимость графиков ядер/областей значений
• Непрерывность нетривиального пересечения $E \cap F \ne \{0\}$
• Замкнутость суммы $E+F$ и её коразмерность

4. Препятствия для передачи через степени

Лемма 3.4: если $G(S_n^m) \xrightarrow{\text{gap}} G(S^m)$, то необходимо $\gamma(S^j)>0$ для всех $1\le j \le m$.

Схема доказательства: если $\gamma(S^k)=0$, то отображение графика $\hat{S}^{k-1}: G(S^{k-1}) \to G(S^k), \quad (x, S^{k-1}x) \mapsto (S^{k-1}x, S^kx)$ не имеет нижней границы. Даже если $G(S_n^{k-1}) \to G(S^{k-1})$, образ $\hat{S}^{k-1}(G(S_n^{k-1})) = G(S_n^k)$ не может сходиться к $G(S^k)$.

Конструирование дискретной диагностической величины

Конечноэлементная база

• Матрица масс: $M_{ij} = \int_\Omega \varphi_i \varphi_j$
• Матрица жёсткости/конвекции: $A_{ij} = a(\varphi_j, \varphi_i)$
• M-скалярное произведение: $\langle u, v \rangle_M := u^* M v$

Дискретный приведённый минимальный модуль

$\gamma_h := \inf_{u \ne 0} \frac{\|(A_h - \lambda M)u\|_{M^{-1}}}{\|u\|_M} = \sigma_{\min}(M^{-1/2}(A_h - \lambda M)M^{-1/2})$

Ключевые свойства:

• M-самосопряжённый случай: $\gamma_h = \text{dist}(\lambda, \sigma(A_h, M))$ (спектральный зазор)
• Неправильный случай: $\gamma_h \ge \text{dist}(\lambda, W_M(A_h))$ (расстояние до числовой области)

Гарантия сходимости (Предложение 3.14)

Если $T$ является m-секториальным оператором и $a_h \xrightarrow{\text{Mosco}} a$, то $\liminf_{h \to 0} \gamma_h \ge \gamma(T - \lambda) > 0$

Технические инновации

1. Единая база: объединяет самосопряжённые (Шрёдингер) и несамосопряжённые (конвективная диффузия) операторы в одну теорию через $\gamma$ вместо спектрального зазора
2. Острота количественных условий:
   • Необходимость: контрпримеры Вольтерра и центральная разность
   • Достаточность: конструктивное доказательство Теоремы 3.2
3. Теоретическое объяснение стабилизированных схем: доказано, что равномерное условие inf-sup для SUPG (Предложение B.1) эквивалентно $\gamma_h^{\text{stab}} \ge c>0$
4. Алгоритмизированная диагностика:
   • Алгоритм 1: адаптивный выбор степени $m$
   • Алгоритм 2: уточнение сетки на основе локальной $\gamma_h(K)$
5. Вычислительная стратегия (Замечание 3.8): избегает явного вычисления $(A_h - \lambda M)^m$ (приводит к заполнению и плохой обусловленности), используя вместо этого методы подпространства Крылова для решения обобщённой задачи на собственные значения

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

Наборы данных (модельные операторы)

1. Оператор Шрёдингера (самосопряжённый случай)

$H = -\frac{d^2}{dx^2} + V(x) \text{ на } L^2(0,1), \quad V \in L^\infty(0,1)$

• 1D: равномерная сетка, P1 конечные элементы
• 2D: единичный квадрат $(0,1)^2$ и L-образная область $(-1,1)^2 \setminus [0,1]^2$
• Тестовые параметры: $V \equiv 25$, $\lambda = -1$ или $\lambda = 25$

2. Оператор конвективной диффузии (несамосопряжённый случай)

$Lu = -\varepsilon u'' + \beta u' + cu \text{ на } (0,1), \quad u(0)=u(1)=0$

• Диапазон параметров: $\varepsilon \in [10^{-5}, 1]$, $\beta \in [8, 50]$, $c=0$
• Дискретизация:
  • Стандартный Галёркин (центральная разность)
  • SUPG стабилизация: $\delta_K = \delta h_K / \|\beta\|_{L^\infty(K)}$, $\delta \in (0, 1/2]$

3. Оператор чистого переноса (предельный случай)

$L_0 u = u', \quad u(0)=0$

• Центральная разность vs. одношаговая схема против потока
• Проверка $\gamma_h \to 0$ (центр) vs. $\gamma_h \ge c>0$ (против потока)

Метрики оценки

1. Дискретный приведённый минимальный модуль: $\gamma_h = \sigma_{\min}(M^{-1/2}(A_h - \lambda M)M^{-1/2})$
2. Расстояние до числовой области: $\text{dist}(\lambda, W_M(A_h))$
3. Число обусловленности: $\kappa(A_h - \lambda M)$
4. Индексы подъёма/спуска: $\text{asc}(A_h)$, $\text{dsc}(A_h)$ (через проверку ранга $\text{rank}(A_h^{m+1}) = \text{rank}(A_h^m)$)
5. Скорость сходимости: наклон $|\gamma_\infty - \gamma_h|$ на логарифмическом графике

Детали реализации

• Уточнение сетки: $h = 2^{-k}$, $k=4,\ldots,9$ (1D); $N \times N$, $N=20,40,80$ (2D)
• Вычисление сингулярных значений: обратная итерация ARPACK, допуск $10^{-8}$
• Время CPU (Таблица 5): $N=10^5$ узлов, $m=3$ примерно 31 секунда
• Разложение Холецкого: предобработка матрицы масс $M$
• Адаптивное уточнение: порог маркировки $\theta = \tau \cdot \min_K \gamma_h(K)$

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

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

1. Оператор Лапласа (Таблица 4, Рисунок 2)

• $H = -\partial_x^2$, $\lambda = 25 \in (\pi^2, 4\pi^2)$
• Дискретные собственные значения: $\zeta_1(h) \downarrow \pi^2 \approx 9.87$, $\zeta_2(h) \downarrow 4\pi^2 \approx 39.48$
• $\gamma_h = \min\{|25-\zeta_1(h)|, |\zeta_2(h)-25|\} \to 14.48$ (монотонность Рэлея-Ритца)
• Сходимость: при $h=2^{-4}$ $\gamma_h=14.88$, при $h=2^{-8}$ $\gamma_h=15.13$

2. Конвективная диффузия (Таблица 6, Рисунок 3)

• $\varepsilon=0.02$, $\beta=8$, $c=0$, $\lambda=-1$
• $\gamma_h$ возрастает с 6.8 ($h=2^{-4}$) до 8.7 ($h=2^{-8}$)
• Расстояние до числовой области: $\text{dist}(\lambda, W_M(A_h))$ возрастает с 6.5 до 8.5 (немного ниже $\gamma_h$)

3. Робастность SUPG стабилизации (Таблица 3)

Ключевое открытие: SUPG сохраняет $\gamma_h \approx 8$ при $\varepsilon \to 0$, тогда как без стабилизации $\gamma_h \to 0$.

4. Высокое число конвекции (Рисунок 4b)

• $\beta=50$, $\varepsilon \in [10^{-3}, 1]$, $h=1/200$
• $\min_\varepsilon \gamma_h \approx 79.38 > 0$ (равномерная нижняя граница)

Абляционные эксперименты

Зависимость от степени (Таблица 2, Замечание 3.6)

Схема против потока для $Lu=u'$ ($\lambda=0$):

Интерпретация: несмотря на $\gamma(A_h) \ge c>0$, имеем $\gamma(A_h^m) \sim h^{m-1} \to 0$ (приближение интегрального оператора Вольтерра), подтверждая необходимость Леммы 3.4.

Проверка индекса подъёма (Таблица 11, Пример 4.16)

Смешанные граничные условия $T=-d^2/dx^2$, $u(0)=u'(1)=0$, $\text{asc}(T)=2$:

• Для всех $h \in \{2^{-5}, \ldots, 2^{-9}\}$: $\text{asc}(A_h)=2$ (правильно)
• $\gamma_h \in [2.41, 2.48]$ (равномерно положительно)

Анализ конкретных случаев

Отказ центральной разности (Предложение 4.12, Таблица 10)

$L_0 u = u', \quad \gamma_h \sim 4h \to 0$

• Непрерывный: $\text{asc}(L_0)=1$, $0 \notin \sigma_{\text{asc}}(L_0)$
• Дискретный: $\text{asc}(A_h)=\infty$, $0 \in \sigma_{\text{asc}}(A_h)$ (катастрофическое несоответствие)

Физическая интерпретация: центральная разность лишена численной диссипации, высокочастотные моды (длина волны $\sim h$) загрязняют структуру ядра/области значений.

Стабилизация против потока

$\gamma_h \ge c > 0 \Rightarrow \text{asc}(A_h)=1$ (восстановление правильного поведения)

Двумерное расширение

L-образная область (Таблица 8, Рисунки 5-6)

• Сингулярность повторного входа: регулярность решения $u \in H^{1+\pi/3-\varepsilon}(\Omega)$
• Сходимость $\gamma_h$: 2.31 ($h=2^{-5}$) $\to 2.58 ($h=2^{-9}$)
• Скорость сходимости: $|\gamma_\infty - \gamma_h| \sim O(h)$ (соответствует пониженной регулярности)
• Локальная вариация: $\gamma_h(K) \in [2.31, 2.56]$ (только 11% разницы)

Вывод: даже при геометрической сингулярности $\gamma_h$ остаётся равномерно положительным, подтверждая робастность базы.

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

1. Теория спектра операторов

• Каасхук-Тейлор 20,27: характеризация подпространств подъёма/спуска
• Беркани 10,11: теория B-Фредгольма и существенный спектр
• Като 21: теория возмущений и сходимость графиков

2. Численный анализ

• Бабушка-Озборн 8: методы конечных элементов для задач на собственные значения (сходимость резольвенты по норме)
• Брукс-Хьюз 12: SUPG стабилизация
• Эрн-Геермон 14: обобщённая сходимость Мосько для несогласованных/стабилизированных схем

3. Задачи конвективной диффузии

• Методы числовой области: Густафсон-Рао 16, Трефетен-Эмбри 28
• Условия inf-sup: Хьюстон и др. 17 (методы разрывного Галёркина)

Преимущества данной работы

• Впервые установлены острые необходимые и достаточные условия стабильности спектра подъёма/спуска при SRS
• Единое рассмотрение самосопряжённых и несамосопряжённых операторов
• Предоставлены вычислимые диагностические величины $\gamma_h$ (существующие работы часто ограничиваются абстрактной сходимостью)
• Объяснены механизмы спектральной стабильности стабилизированных схем типа SUPG

Заключение и обсуждение

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

1. Необходимое и достаточное условие: $\liminf_{h \to 0} \gamma_h > 0 \iff$ стабильность спектра подъёма/спуска
2. Передача через степени: требуется проверка всех промежуточных степеней $\gamma((T-\lambda)^j)>0$ ($1\le j \le m$)
3. Роль стабилизации: SUPG гарантирует равномерное условие inf-sup, обеспечивающее $\gamma_h^{\text{stab}} \ge c>0$ даже при $\varepsilon \to 0$
4. Диагностическая величина $\gamma_h$: эффективно вычисляется (методы Крылова), предсказывает успех/отказ дискретизации

Ограничения

1. Чистый конвективный предел ($\varepsilon=0$):
   • Гипотеза 3.17 лишена строгого доказательства (требует выхода за рамки теории Мосько)
   • Численные доказательства убедительны (Таблица 3), но теоретический разрыв существует
2. Предположения о сетке:
   • Требуется квазиравномерность
   • Стратегии локального уточнения для неглад ких областей (Алгоритм 2) нуждаются в дополнительной проверке
3. Вычисления в высоких размерностях:
   • При 3D и $N \sim 10^6$ стоимость разложения Холецкого $O(N^2)$
   • Требуется алгебраический многосеточный (AMG) предобработчик
4. Несекториальные операторы:
   • Пример 4.15 (Гельмгольц) является лишь предварительным исследованием
   • Гарантии SRS для общих несекториальных операторов неясны

Будущие направления

1. Теоретическое совершенствование:
   • Доказательство Гипотезы 3.17 (сходимость графиков при чистой конвекции)
   • Расширение на блочные структурированные операторы (Стокс, Максвелл)
2. Оптимизация алгоритмов:
   • Параллелизация вычисления $\gamma_h^{(m)}$
   • Строгие оценки ошибок для адаптивного уточнения
3. Расширение приложений:
   • Полудискретная стабильность для задач, зависящих от времени
   • Локальная линеаризация нелинейных операторов
4. Программная реализация:
   • Интеграция в модули спектрального анализа FEniCS/deal.II

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

Достоинства

1. Теоретическая инновация

• Острота: необходимое и достаточное условие $\gamma>0$ не может быть ослаблено (контрпример Вольтерра)
• Единство: самосопряжённые/несамосопряжённые, секториальные/несекториальные операторы объединены в одну базу
• Глубина: раскрыта сущность $\gamma$ как "стационарного условия CFL"

2. Практическая ценность

• Операционность: вычисление $\gamma_h$ требует только стандартной линейной алгебры (ARPACK)
• Предсказательность: Таблица 3 чётко различает успешные (SUPG) и неудачные (центральная разность) схемы
• Руководство: Алгоритм 2 предоставляет стратегию адаптивного уточнения сетки

3. Полнота экспериментов

• Многомерная проверка: 1D/2D, гладкие/сингулярные области, $\varepsilon \in [10^{-8}, 1]$
• Количественные сравнения: $\gamma_h$ vs. расстояние до числовой области vs. число обусловленности (Таблица 9)
• Полные контрпримеры: центральная разность, Вольтерра, затухание степеней против потока

4. Ясность изложения

• Логическая цепь (Рисунок 1): численный отказ $\to$ диагностическая величина $\to$ критерий подпространств $\to$ механизм сходимости
• Богатые таблицы: Таблица 1 (соответствие непрерывное-дискретное), Таблица 2 (затухание степеней)
• Подробные приложения: Приложение A (необходимость), B (доказательство SUPG), C (численная стабильность)

Недостатки

1. Теоретические разрывы

• Чистая конвекция: Гипотеза 3.17 не доказана, ограничивает полноту теории
• Несекториальные операторы: гарантии SRS в несекториальном случае зависят от конкретного оператора (Пример 4.15 — лишь частный случай)

2. Вычислительные затраты

• Узкое место в высоких размерностях: вычисление $\gamma_h$ для 3D крупномасштабных задач ($N>10^6$) недостаточно обсуждено
• Вычисление степеней: критерий остановки Алгоритма 1 $|\gamma_h^{(m+1)}/\gamma_h^{(m)} - 1| < \varepsilon_{\text{tol}}$ может отказать для плохо обусловленных задач

3. Ограничения экспериментов

• Ограниченные 2D результаты: только Лапласиан (Таблицы 7-8), отсутствуют подробные данные для 2D конвективной диффузии
• L-образная область: локальная вариация $\gamma_h(K)$ мала (11%), не демонстрирует значительное преимущество адаптивного уточнения

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

• Линейные предположения: нелинейные УЧП требуют локальной линеаризации, передача стабильности неясна
• Зависимость от времени: полудискретная стабильность не рассмотрена

Влияние

Вклад в область

• Численный анализ: первый вычислимый критерий стабильности спектра подъёма/спуска
• Теория операторов: связь теории B-Фредгольма с дискретными приближениями
• Вычисления УЧП: объяснение механизмов спектральной стабильности стабилизированных схем

Практическая ценность

• Инженерные приложения: руководство по выбору схем для конвективно-доминирующих задач (например, CFD)
• Разработка ПО: $\gamma_h$ может служить стандартной диагностической величиной в библиотеках конечных элементов

Воспроизводимость

• Отсутствие кода: открытая реализация не предоставлена (но алгоритмы описаны подробно)
• Воспроизводимость данных: параметры Таблиц 4-11 полные, позволяют независимую проверку

Применимые сценарии

Наиболее подходящие

1. Несамосопряжённые УЧП: конвективная диффузия, Гельмгольц, магнитная гидродинамика
2. Сингулярные возмущения: задачи с пограничными слоями при $\varepsilon \ll 1$
3. Оценка схем: сравнение различных стабилизированных методов (SUPG, DG, против потока)

Требуют осторожности

1. Чистая конвекция ($\varepsilon=0$): теория не полностью охватывает
2. Нелинейные задачи: требуется дополнительная теоретическая поддержка
3. Сверхбольшие масштабы ($N>10^7$): вычислительные затраты могут быть чрезмерны

Неприменимо

1. Задачи с эволюцией по времени: требуется расширение на полудискретизацию
2. Сильная нелинейность: например, турбулентность Навье-Стокса
3. Стохастические УЧП: распространение неопределённости не рассмотрено

Ключевые ссылки

1. Като 21: Perturbation Theory for Linear Operators — классический учебник по теории возмущений и сходимости графиков
2. Каасхук 20, Тейлор 27: характеризация подпространств подъёма/спуска
3. Беркани 10,11: теория B-Фредгольма
4. Брукс-Хьюз 12: оригинальная работа по SUPG стабилизации
5. Эрн-Геермон 14: современный справочник по теории конечных элементов
6. Бабушка-Озборн 8: авторитетный обзор численных методов для задач на собственные значения

Итоговое резюме

Статья устанавливает первое острое необходимое и достаточное условие стабильности спектра подъёма/спуска несамосопряжённых операторов в рамках сильной сходимости резольвенты, связывая абстрактную теорию операторов с практическим численным анализом через вычислимую диагностическую величину $\gamma_h$. Теория строга (необходимость и достаточность, контрпримеры, расширение B-Фредгольма), эксперименты полны (1D/2D, многопараметрические, стабилизированные схемы), изложение ясно. Основные вклады:

1. Раскрытие $\gamma>0$ как количественного порога спектральной стабильности
2. Объяснение успеха методов типа SUPG
3. Предоставление практических алгоритмических инструментов

Ограничения включают теоретический разрыв в случае чистой конвекции и вычислительные затраты в высоких размерностях. В целом, это высококачественная работа на пересечении численного анализа и теории операторов, имеющая важное значение для спектрального анализа дискретизаций УЧП.