Аннотация:
Пусть $\mathcal{O}\subset\mathbb{R}^d$ – ограниченная область с границей класса $C^{1,1}$. В пространстве $L_2(\mathcal{O};\mathbb{C}^n)$ рассматривается самосопряженный матричный эллиптический дифференциальный оператор $B_{N,\varepsilon}$, $0<\varepsilon\leqslant1$, второго порядка при условии Неймана на границе. Старшая часть оператора задана в факторизованной форме. Оператор включает члены первого и нулевого порядков. Коэффициенты оператора $B_{N,\varepsilon}$ периодичны и зависят от $\mathbf{x}/\varepsilon$. Изучается обобщенная резольвента $(B_{N,\varepsilon}-\zeta Q_0(\,{\cdot}\,/\varepsilon))^{-1}$, где $Q_0$ – периодическая ограниченная и положительно определенная матрица-функция, а $\zeta$ – комплексный параметр. Получены аппроксимации обобщенной резольвенты по операторной норме в $L_2(\mathcal{O};\mathbb{C}^n)$ и по норме операторов, действующих из $L_2(\mathcal{O};\mathbb{C}^n)$ в класс Соболева $H^1(\mathcal{O};\mathbb{C}^n)$, с двухпараметрическими (относительно $\varepsilon$ и $\zeta$) оценками погрешности. Результаты применяются к изучению поведения решений начально-краевой задачи с условием Неймана для параболического уравнения $Q_0(\mathbf{x}/\varepsilon) \, \partial_t \mathbf{u}_\varepsilon(\mathbf{x},t) = -(B_{N,\varepsilon} \mathbf{u}_\varepsilon)(\mathbf{x},t)$ в цилиндре $\mathcal{O} \times (0,T)$, где $0<T \le \infty$.
Библиография: 51 наименование.
Ключевые слова:
периодические дифференциальные операторы, эллиптические системы, параболические системы, усреднение, операторные оценки погрешности.
Тематика работы – теория усреднения (гомогенизации) периодических дифференциальных операторов (ДО), которой посвящена обширная литература. Укажем в первую очередь книги [1]–[4].
0.1. Класс операторов
Пусть $\Gamma$ – решетка в $\mathbb{R}^d$, $\Omega$ – ячейка решетки $\Gamma$. Для $\Gamma$-периодических функций используем обозначения $\psi^\varepsilon (\mathbf{x}):=\psi (\mathbf{x}/\varepsilon)$, $\overline{\psi}:=|\Omega|^{-1} \int_\Omega \psi (\mathbf{x})\,d\mathbf{x}$. Обозначим $\mathbf{D} := - i \nabla$.
Пусть $\mathcal{O}\subset\mathbb{R}^d$ – ограниченная область с границей класса $C^{1,1}$. В пространстве $L_2(\mathcal{O};\mathbb{C}^n)$ изучается самосопряженный матричный сильно эллиптический ДО второго порядка $\mathcal{B}_{N,\varepsilon}$, $0<\varepsilon\leqslant 1$, при условии Неймана на границе. Старшая часть $A_{N,\varepsilon}$ оператора $\mathcal{B}_{N,\varepsilon}$ задается в факторизованной форме $A_{N,\varepsilon}=b(\mathbf{D})^*g^\varepsilon (\mathbf{x})b(\mathbf{D})$, где $b(\mathbf{D})$ – матричный однородный ДО первого порядка, $g(\mathbf{x})$ – $\Gamma$-периодическая матрица-функция в $\mathbb{R}^d$, ограниченная и положительно определенная. Задача усреднения для оператора $A_{N,\varepsilon}$ изучалась в работах [5], [6]. Сейчас мы рассматриваем более общий класс самосопряженных ДО $\mathcal{B}_{N,\varepsilon}$, включающих младшие члены:
Здесь $a_j(\mathbf{x})$, $j=1,\dots,d$, и $Q(\mathbf{x})$ – $\Gamma$-периодические матрицы-функции, вообще говоря, неограниченные. Строгое определение оператора $\mathcal{B}_{N,\varepsilon}$ дается через соответствующую квадратичную форму, определенную на классе Соболева $H^1(\mathcal{O};\mathbb{C}^n)$. Делаются предположения, обеспечивающие сильную эллиптичность оператора $\mathcal{B}_{N,\varepsilon}$.
Коэффициенты оператора (0.1) быстро осциллируют при малом $\varepsilon$. Типичная задача теории усреднения применительно к оператору $\mathcal{B}_{N,\varepsilon}$ состоит в нахождении аппроксимации при $\varepsilon\to 0$ для резольвенты $(\mathcal{B}_{N,\varepsilon}-\zeta I)^{-1}$ либо обобщенной резольвенты $(\mathcal{B}_{N,\varepsilon}-\zeta Q_0^\varepsilon)^{-1}$. Здесь $Q_0(\mathbf{x})$ – $\Gamma$-периодическая матрица-функция, положительно определенная, ограниченная и ограниченно обратимая.
Пример оператора вида (0.1) рассмотрен в § 11 – это оператор Шрёдингера (11.1) с периодическими быстро осциллирующими коэффициентами (метрикой, магнитным потенциалом и электрическим потенциалом, притом содержащим сингулярное слагаемое $\varepsilon^{-1} v^\varepsilon(\mathbf{x})$).
0.2. Обзор результатов по операторным оценкам погрешности
В серии работ [7]–[9] Бирман и Суслина разработали теоретико-операторный (спектральный) подход к задачам теории усреднения. Изучались операторы
действующие в $L_2(\mathbb{R}^d;\mathbb{C}^n)$. Здесь $g(\mathbf{x})$ – $\Gamma$-периодическая $(m\times m)$-матрица-функция в $\mathbb{R}^d$, ограниченная и положительно определенная; матричный ДО $b(\mathbf{D})$ имеет вид $b(\mathbf{D}) = \sum_{j=1}^d b_j D_j$, где $b_j$ – постоянные $(m \times n)$-матрицы. Предполагается, что $m \geqslant n$ и символ $b(\boldsymbol{\xi})$ имеет максимальный ранг. Многие операторы математической физики допускают запись в виде (0.2). Простейший пример – это оператор акустики $A_\varepsilon = - \operatorname{div}g^\varepsilon (\mathbf{x}) \nabla = \mathbf{D}^* g^\varepsilon (\mathbf{x}) \mathbf{D}$; в этом случае $m=d$, $n=1$ и $b(\mathbf{D}) = \mathbf{D}$. Оператор теории упругости также допускает запись в виде (0.2). Другой пример возникает при изучении системы Максвелла – это оператор вида $A_\varepsilon = \operatorname{rot} \eta^\varepsilon(\mathbf{x}) \operatorname{rot} - \operatorname{\nabla} \nu^\varepsilon(\mathbf{x}) \operatorname{div}$ (при $d=3$). Эти и другие примеры подробно рассмотрены в [7].
В [7] было показано, что при $\varepsilon\to 0$ резольвента $(A_\varepsilon +I)^{-1}$ сходится по операторной норме в $L_2(\mathbb{R}^d;\mathbb{C}^n)$ к резольвенте эффективного оператора $A^0=b(\mathbf{D})^*g^0b(\mathbf{D})$. Здесь $g^0$ – постоянная положительная эффективная матрица. Была установлена оценка
В [9] была получена аппроксимация резольвенты $(A_\varepsilon +I)^{-1}$ по норме операторов, действующих из $L_2(\mathbb{R}^d;\mathbb{C}^n)$ в класс Соболева $H^1(\mathbb{R}^d;\mathbb{C}^n)$:
В этой аппроксимации учтен корректор $K(\varepsilon)$. Оператор $K(\varepsilon)$ содержит быстро осциллирующие множители, а потому зависит от $\varepsilon$. При этом $\|\varepsilon K(\varepsilon)\|_{L_2\to H^1}=O(1)$.
Оценки (0.3), (0.4) точны по порядку. Постоянные в оценках контролируются явно в терминах данных задачи. Подобные результаты получили название операторных оценок погрешности в теории усреднения. Метод работ [7]–[9] основан на применении масштабного преобразования, теории Флоке–Блоха и аналитической теории возмущений.
Впоследствии спектральный метод был распространен Суслиной [10], [11] на случай оператора
действующего в $L_2(\mathbb{R}^d;\mathbb{C}^n)$. Удобно фиксировать вещественный параметр $\lambda$ так, чтобы оператор $B_\varepsilon := \mathcal{B}_\varepsilon+\lambda Q_0^\varepsilon$ был положительно определен. В [10] установлены аналоги оценок (0.3), (0.4):
Здесь $B^0$ – соответствующий эффективный оператор и $\mathcal{K}(\varepsilon)$ – корректор.
К параболическим системам спектральный метод применялся в работах Суслиной [12]–[14]. В [12], [13] был найден старший член аппроксимации операторной экспоненты $e^{-A_\varepsilon t}$, $t>0$, по $(L_2 \to L_2)$-норме, а в [14] была установлена аппроксимация экспоненты по $(L_2 \to H^1)$-норме при учете корректора:
Экспонента от оператора $B_\varepsilon$ изучалась в работе Мешковой [15], где установлены аналоги неравенств (0.8) и (0.9).
Другой подход к получению операторных оценок погрешности в теории усреднения был предложен Жиковым и Пастуховой. В работах [16]–[18] были получены оценки вида (0.3), (0.4) для операторов акустики и теории упругости. Метод, названный авторами “модифицированным методом первого приближения” или “методом сдвига”, основан на анализе первого приближения к решению и введении в задачу дополнительного параметра. Помимо задач в $\mathbb{R}^d$ в работах [17], [18] изучались задачи усреднения в ограниченной области $\mathcal{O}\subset \mathbb{R}^d$ при условии Дирихле либо Неймана на границе. К параболическим уравнениям метод сдвига применялся в [19]. Дальнейшие результаты Жикова, Пастуховой и их учеников отражены в обзоре [20].
В присутствии членов младшего порядка задача усреднения для оператора (0.5) в $\mathbb{R}^d$ изучалась в статье Борисова [21]. Было найдено выражение для эффективного оператора $\mathcal{B}^0$ и получены оценки погрешности вида (0.6), (0.7). При этом предполагалось, что коэффициенты оператора зависят не только от быстрой, но и от медленной переменной. Однако в [21] коэффициенты оператора $\mathcal{B}_\varepsilon$ предполагались достаточно гладкими. Упомянем также работу Сеника [22], в которой изучался несамосопряженный сильно эллиптический оператор второго порядка (с включением членов младшего порядка) на бесконечном цилиндре. Коэффициенты периодичны вдоль цилиндра и быстро осциллируют; получены оценки вида (0.6), (0.7).
Операторные оценки погрешности для эллиптических уравнений второго порядка (без младших членов) в ограниченной области при условиях Дирихле или Неймана изучались многими авторами. Первыми, по-видимому, были Москоу и Вогелиус, установившие оценку (см. [23; следствие 2.2]), допускающую запись в операторных терминах:
Здесь оператор $A_{D,\varepsilon}$ в $L_2(\mathcal{O})$, $\mathcal{O}\subset \mathbb{R}^2$, задан выражением $-\operatorname{div} g^\varepsilon (\mathbf{x})\nabla$ при условии Дирихле на $\partial\mathcal{O}$, а матрица-функция $g(\mathbf{x})$ предполагается $C^\infty$-гладкой. В случае условия Неймана аналогичная оценка получена в [24; следствие 1]. В той же статье найдена аппроксимация с корректором для обратного оператора по норме операторов, действующих из $L_2(\mathcal{O})$ в класс Соболева $H^1(\mathcal{O})$, с оценкой погрешности порядка $O(\sqrt{\varepsilon}\,)$. Ухудшение порядка по сравнению с аналогичным результатом в $\mathbb{R}^d$ объясняется влиянием границы области.
В случае произвольной размерности эллиптические задачи в ограниченной области с достаточно гладкой границей изучались в работах [17] и [18]. Гладкость коэффициентов не предполагалась. Для операторов акустики и упругости при условии Дирихле либо Неймана на границе была получена $(L_2\to H^1)$-аппроксимация резольвенты при учете корректора с оценкой погрешности порядка $O(\sqrt{\varepsilon}\,)$. В качестве следствия была установлена оценка вида (0.10), но с погрешностью порядка $O(\sqrt{\varepsilon}\,)$. Близкие результаты для оператора, заданного выражением $-\operatorname{div} g^\varepsilon (\mathbf{x})\nabla$ в ограниченной области $\mathcal{O}\subset\mathbb{R}^d$ при условии Дирихле либо Неймана на $\partial\mathcal{O}$, были установлены в работах Гризо [25], [26] с помощью метода “анфолдинга”. В [26] для того же оператора впервые была получена точная по порядку оценка (0.10). Для эллиптических систем сходные результаты независимо получены в [27]–[29]. Дальнейшие продвижения и подробный обзор можно найти в работах [5], [6].
Для матричного эллиптического оператора второго порядка (при включении младших членов) в ограниченной области при краевых условиях Дирихле либо Неймана операторные оценки погрешности установлены в работах Ху [30]–[32]. Однако в этих статьях на оператор наложено весьма жесткое условие равномерной эллиптичности. Мы сравним наши результаты с результатами [32] ниже в п. 0.3.
Упомянем также монографию Шена [33], статью [34] и цитированную там литературу.
До сих пор речь шла об аппроксимации резольвент в фиксированной регулярной точке. Аппроксимация резольвенты $(A_\varepsilon -\zeta I)^{-1}$ оператора (0.2) в зависимости от $\varepsilon$ и спектрального параметра $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$ найдена Суслиной [6]. В этой работе также получены двухпараметрические (относительно $\varepsilon$ и $\zeta$) оценки погрешности при усреднении резольвент операторов $A_{D,\varepsilon}$ и $A_{N,\varepsilon}$ вида (0.2), действующих в ограниченной области при условии Дирихле либо Неймана на границе.
Стимулом к получению двухпараметрических оценок послужило изучение усреднения параболических систем. Аппроксимация экспоненты от операторов $A_{D,\varepsilon}$ и $A_{N,\varepsilon}$ найдена в работе Мешковой и Суслиной [35]:
где $\gamma \subset \mathbb{C}$ – контур, обходящий спектр оператора $A_{\flat,\varepsilon}$ в положительном направлении. Это тождество позволяет вывести аппроксимации операторной экспоненты $e^{-A_{\flat,\varepsilon}t}$ из соответствующих аппроксимаций резольвенты $(A_{\flat,\varepsilon}-\zeta I)^{-1}$ с двухпараметрическими (относительно $\varepsilon $ и $\zeta $) оценками погрешности.
Оператор с коэффициентами, периодическими по пространственным переменным и по времени, изучался Генгом и Шеном [36]. В [36] получены операторные оценки погрешности для уравнения
Здесь $\phi =\operatorname{arg}\zeta\in (0,2\pi)$, $|\zeta|\geqslant 1$. Зависимость констант в оценках от угла $\phi$ прослежена. Оценки (0.11), (0.12) равномерны по углу $\phi$ в любой области вида
при сколь угодно малом $\phi_0 >0$. В [37] получены также и оценки при $|\zeta|<1$, $\phi \in (0,2\pi)$.
Отметим, что в последние годы операторные оценки погрешности в различных задачах гомогенизации для дифференциальных операторов привлекают внимание все большего числа исследователей; получено много содержательных результатов. Достаточно подробный обзор современного состояния этой области можно найти в [38; введение].
0.3. Основные результаты
В настоящей работе изучается усреднение эллиптических и параболических задач в ограниченной области при условии Неймана для сильно эллиптического оператора второго порядка вида (0.1) (с включением младших членов). Как уже отмечалось, ранее такие задачи были исследованы в [5], [6], [35] для оператора $A_{N,\varepsilon}$ (без младших членов).
Прежде чем формулировать результаты, удобно перейти к положительно определенному оператору $B_{N,\varepsilon}=\mathcal{B}_{N,\varepsilon}+\lambda Q_0^\varepsilon$, выбирая подходящую постоянную $\lambda$. Пусть $B_N^0$ – соответствующий эффективный оператор.
справедливые при $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и достаточно малом $\varepsilon$. Величины $C(\phi)$ контролируются явно в терминах данных задачи и угла $\phi$. Оценки (0.14), (0.15) равномерны по $\phi$ в любой области вида (0.13) при сколь угодно малом $\phi_0 >0$.
При фиксированном $\zeta$ оценка (0.14) имеет точный порядок $O(\varepsilon)$. Порядок оценки (0.15) хуже, чем в $\mathbb{R}^d$ (ср. (0.7)), из-за влияния границы области.
Корректор в (0.15) в общем случае содержит сглаживающий оператор. Мы выделяем случаи, когда можно использовать более простой корректор. Помимо оценок для обобщенной резольвенты мы находим также аппроксимации по $(L_2\to L_2)$-норме для операторов $g^\varepsilon b(\mathbf{D})(B_{N,\varepsilon}-\zeta Q_0^\varepsilon )^{-1}$, отвечающих потокам. Показано, что в строго внутренней подобласти $\mathcal{O}'$ области $\mathcal{O}$ можно получить аппроксимацию оператора $(B_{N,\varepsilon}-\zeta Q_0^\varepsilon )^{-1}$ по $(L_2(\mathcal{O}) \to H^1(\mathcal{O}'))$-норме с оценкой погрешности точного порядка $O(\varepsilon)$.
Мы находим также аппроксимации оператора $(B_{N,\varepsilon}-\zeta Q_0^\varepsilon)^{-1}$, справедливые в более широкой области изменения спектрального параметра $\zeta$, с оценками погрешности, имеющими другое поведение относительно $\zeta$ (подробнее см. § 9).
Двухпараметрические оценки (0.14), (0.15) применяются к изучению поведения решения начально-краевой задачи для параболического уравнения
при естественном условии (условии Неймана) на $\partial \mathcal{O} \times \mathbb{R}_+$. Здесь $\boldsymbol{\varphi}\in L_2(\mathcal{O};\mathbb{C}^n)$. Решение понимается в обобщенном смысле как решение из “энергетического” класса (ср. [39; гл. 3]).
Показано, что при $\varepsilon\to 0$ решение $\mathbf{u}_\varepsilon (\,{\cdot}\,,t)$ сходится в $L_2(\mathcal{O};\mathbb{C}^n)$ к решению $\mathbf{u}_0(\,{\cdot}\,,t)$ эффективной задачи с постоянными коэффициентами:
справедливую при достаточно малом $\varepsilon$. При фиксированном значении времени $t>0$ эта оценка имеет точный порядок $O(\varepsilon)$. Второй результат для задачи (0.16) – аппроксимация решения $\mathbf{u}_\varepsilon (\,{\cdot}\,,t)$ по энергетической норме:
Здесь $\mathbf{v}_\varepsilon (\,{\cdot}\,,t)=\mathbf{u}_0(\,{\cdot}\,,t)+\varepsilon \mathcal{K}_N(t;\varepsilon )\boldsymbol{\varphi}(\,{\cdot}\,)$ – первое приближение к решению $\mathbf{u}_\varepsilon (\,{\cdot}\,,t)$, оператор $\mathcal{K}_N(t;\varepsilon)$ – корректор. При фиксированном $t$ оценка (0.18) имеет порядок $O(\varepsilon^{1/2})$ из-за влияния границы.
В общем случае корректор содержит сглаживающий оператор. Мы выделяем условия, при которых можно использовать более простой корректор без сглаживающего оператора. Помимо оценки (0.18) получаем аппроксимацию потока $g^\varepsilon b(\mathbf{D})\mathbf{u}_\varepsilon (\,{\cdot}\,,t)$ по $L_2$-норме. В строго внутренней подобласти $\mathcal{O}'\subset \mathcal{O}$ найдена аппроксимация решения $\mathbf{u}_\varepsilon(\,{\cdot}\,,t)$ по норме в $H^1(\mathcal{O}'; \mathbb{C}^n)$ с погрешностью точного порядка $O(\varepsilon)$.
Оценки (0.17) и (0.18) можно переписать в равномерной операторной топологии. В более простом случае, когда $Q_0(\mathbf{x})=\mathbf{1}_n$, формулировки таковы:
Аналогичные результаты для оператора $B_{D,\varepsilon}$ с условием Дирихле получены в работах Мешковой и Суслиной: статья [40] (см. также [41]) посвящена эллиптической задаче Дирихле, а работа [42] посвящена первой начально-краевой задаче для параболического уравнения.
Отметим, что двухпараметрические оценки погрешности для эллиптических задач уже нашли применение к получению операторных оценок погрешности не только в параболических, но и в гиперболических задачах (см. [43]).
Сопоставим наши результаты для эллиптических задач с результатами близкой статьи [32]. Перечислим наши преимущества. Во-первых, мы изучаем сильно эллиптический оператор вида (0.1), а в [32] (как и в [27], [30], [31]) на оператор накладывается весьма ограничительное условие равномерной эллиптичности. Во-вторых, мы включаем в рассмотрение младшие члены с неограниченными коэффициентами (из подходящих $L_p(\Omega)$-классов), а в [32] коэффициенты младших членов предполагаются ограниченными. В-третьих, мы получаем двухпараметрические оценки погрешности (относительно $\varepsilon$ и $\zeta$), а в [32] оценки однопараметрические (по параметру $\varepsilon$). С другой стороны, в [32] некоторые результаты получены в областях с липшицевой границей; допускаются несамосопряженные операторы (с самосопряженной старшей частью).
0.4. Метод доказательства
Для доказательства оценок (0.14), (0.15) используется метод из работ [5], [6]. Он основан на рассмотрении ассоциированной задачи в $\mathbb{R}^d$, использовании оценок (0.11), (0.12) (полученных в [37]), введении поправки типа пограничного слоя и ее тщательном анализе. Существенную техническую роль играет использование сглаживания по Стеклову (заимствованное из работы [18]) и оценок в $\varepsilon$-окрестности границы. Зависимость от спектрального параметра в оценках тщательно отслеживается. Присутствие младших членов с неограниченными коэффициентами вносит дополнительные технические трудности (по сравнению с [6]). Сначала рассматривается случай $\operatorname{Re}\zeta \leqslant 0$, $|\zeta|\geqslant 1$. Мы доказываем оценку (0.15), а затем оценку (0.14), опираясь на неравенство (0.15) и соображения двойственности. Затем переносим уже доказанные оценки из точки $\zeta$ в левой полуплоскости в симметричную точку правой полуплоскости, используя подходящие тождества для резольвент.
Аппроксимации, справедливые в более широкой области изменения параметра $\zeta$, выводятся из уже доказанных оценок в точке $\zeta =-1$ и подходящих резольвентных тождеств.
0.5. Структура работы
Работа состоит из одиннадцати параграфов. В § 1 вводится класс операторов $B_\varepsilon$, действующих в $L_2(\mathbb{R}^d;\mathbb{C}^n)$, и формулируются результаты усреднения обобщенной резольвенты $(B_\varepsilon -\zeta Q_0^\varepsilon)^{-1}$, полученные в [37]. В § 2 описывается класс операторов $B_{N,\varepsilon}$ и определяется эффективный оператор $B_N^0$. В § 3 формулируются основные результаты работы, вводится поправка типа пограничного слоя и устанавливается теорема об аппроксимации решения $\mathbf{u}_\varepsilon = (B_{N,\varepsilon} -\zeta Q_0^\varepsilon)^{-1} \mathbf{F}$ при учете этой поправки. В § 4 содержится вспомогательный материал. В § 5 установлена оценка поправки по $H^1$-норме и найдена аппроксимация (0.15) для обобщенной резольвенты по $(L_2 \to H^1)$-норме в случае $\operatorname{Re}\zeta \leqslant 0$. В § 6 получена оценка поправки по $L_2$-норме и для обобщенной резольвенты получена аппроксимация по $(L_2\to L_2)$-норме с оценкой (0.14) в случае $\operatorname{Re}\zeta \leqslant 0$. В § 7 результаты перенесены в точку $\zeta$ в правой полуплоскости; завершено доказательство основных результатов для обобщенной резольвенты. В § 8 выделены условия, когда сглаживающий оператор в корректоре можно устранить; рассмотрены специальные случаи; получены оценки в строго внутренней подобласти. Оценки, справедливые в более широкой области изменения спектрального параметра, установлены в § 9. Параграф 10 посвящен усреднению решений второй начально-краевой задачи для параболического уравнения. Пример применения общих результатов рассмотрен в § 11.
0.6. Обозначения
Пусть $\mathfrak{H}$, $\mathfrak{H}_*$ – комплексные сепарабельные гильбертовы пространства. Символы $(\,{\cdot}\,,{\cdot}\,)_\mathfrak{H}$ и $\|\,{\cdot}\,\|_\mathfrak{H}$ означают скалярное произведение и норму в $\mathfrak{H}$; символ $\|\,{\cdot}\,\|_{\mathfrak{H}\to\mathfrak{H}_*}$ означает норму линейного непрерывного оператора из $\mathfrak{H}$ в $\mathfrak{H}_*$.
Символы $\langle \,{\cdot}\,,{\cdot}\,\rangle$ и $|\,{\cdot}\,|$ означают соответственно скалярное произведение и норму в $\mathbb{C}^n$, $\mathbf{1}_n$ – единичная $(n\times n)$-матрица. Если $a$ – $(m\times n)$-матрица, то символ $|a|$ означает норму матрицы $a$ как оператора из $\mathbb{C}^n$ в $\mathbb{C}^m$. Для $z\in\mathbb{C}$ через $z^*$ обозначается комплексно-сопряженное число. (Мы используем такое нестандартное обозначение, так как верхняя черта означает среднее значение периодической функции по ячейке периодов.) Используем обозначения $\mathbf{x}=(x_1,\dots, x_d)\in\mathbb{R}^d$, $iD_j=\partial_j =\partial /\partial x_j$, $j=1,\dots,d$, $\mathbf{D}=-i\nabla=(D_1,\dots,D_d)$. Классы $L_p$ функций в области $\mathcal{O}\subset\mathbb{R}^d$ со значениями в $\mathbb{C}^n$ обозначаем через $L_p(\mathcal{O};\mathbb{C}^n)$, $1\leqslant p\leqslant \infty$. Классы Соболева $\mathbb{C}^n$-значных функций в области $\mathcal{O}\subset\mathbb{R}^d$ обозначаются через $H^s(\mathcal{O};\mathbb{C}^n)$. Через $H^1_0(\mathcal{O};\mathbb{C}^n)$ обозначается замыкание класса $C_0^\infty(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathcal{O};\mathbb{C}^n)$. При $n=1$ пишем просто $L_p(\mathcal{O})$, $H^s(\mathcal{O})$ и так далее, но часто мы применяем такие упрощенные обозначения и для пространств вектор-функций или матричнозначных функций.
Используем обозначение $\mathbb{R}_+=[0,\infty)$. Различные оценочные постоянные обозначаются символами $c$, $\mathfrak{c}$, $C$, $\mathcal{C}$, $\mathfrak{C}$, $\mathrm C$, $\beta$, $\gamma$ (возможно, с индексами и значками).
§ 1. Задача усреднения для эллиптического оператора, действующего в $L_2(\mathbb{R}^d;\mathbb{C}^n)$
В этом параграфе формулируются результаты усреднения для эллиптических систем в $\mathbb{R}^d$, полученные в [37].
Через $\widetilde{H}^1(\Omega)$ обозначается подпространство тех функций из $H^1(\Omega)$, $\Gamma$-периодическое продолжение которых на $\mathbb{R}^d$ принадлежит $H^1_{\mathrm{loc}}(\mathbb{R}^d)$. Если $f (\mathbf{x})$ – $\Gamma$-периодическая матрица-функция в $\mathbb{R}^d$, положим $f^\varepsilon (\mathbf{x}):=f (\mathbf{x}/\varepsilon)$, $\varepsilon >0$;
Здесь при определении $\overline{f}$ предполагается, что $f\,{\in}\, L_{1,\mathrm{loc}}(\mathbb{R}^d)$, а при определении $\underline{f}$ считается, что матрица $f$ квадратная и неособая, причем $f^{-1}\in L_{1,\mathrm{loc}}(\mathbb{R}^d)$. Через $[f^\varepsilon ]$ обозначается оператор умножения на матрицу-функцию $f^\varepsilon (\mathbf{x})$.
1.2. Сглаживание по Стеклову
Рассмотрим оператор сглаживания по Стеклову $S_\varepsilon^{(k)}$, действующий в $L_2(\mathbb{R}^d;\mathbb{C}^k)$ (где $k\in\mathbb{N}$) по правилу
Значок $(k)$ будем опускать и писать просто $S_\varepsilon$. Очевидно, $S_\varepsilon \mathbf{D}^\alpha \mathbf{u}=\mathbf{D}^\alpha S_\varepsilon \mathbf{u}$ при $\mathbf{u}\in H^s(\mathbb{R}^d;\mathbb{C}^k)$ для любого мультииндекса $\alpha$ такого, что $|\alpha| \leqslant s$. Отметим неравенство
Предложение 1.2. Пусть $f$ – $\Gamma$-периодическая функция в $\mathbb{R}^d$ такая, что $f\in L_2(\Omega)$. Тогда оператор $[f^\varepsilon ]S_\varepsilon $ непрерывен в $L_2(\mathbb{R}^d)$ и справедлива оценка
В $L_2(\mathbb{R}^d;\mathbb{C}^n)$ рассмотрим оператор $A_\varepsilon$, формально заданный дифференциальным выражением $A_\varepsilon {=}\,b(\mathbf{D})^*g^\varepsilon (\mathbf{x})b(\mathbf{D})$. Здесь $g(\mathbf{x})$ – $\Gamma$-периодическая эрмитова $(m\times m)$-матрица-функция (вообще говоря, с комплексными элементами). Считаем, что $g(\mathbf{x})>0$ и $g,g^{-1}\in L_\infty (\mathbb{R}^d)$. Дифференциальный оператор $b(\mathbf{D})$ имеет вид
где $b_j$, $j=1,\dots,d$, – постоянные матрицы размера $m\times n$ (вообще говоря, с комплексными элементами). Считаем, что $m\geqslant n$ и символ $b(\boldsymbol{\xi})=\sum_{j=1}^d b_j\xi_j$ оператора $b(\mathbf{D})$ имеет максимальный ранг:
Мы изучаем самосопряженный оператор $B_\varepsilon$, старшая часть которого совпадает с $A_\varepsilon$. Чтобы определить младшие члены оператора, введем $\Gamma$-периодические $(n\times n)$-матрицы-функции (вообще говоря, с комплексными элементами) $a_j$, $j=1,\dots,d$, такие, что
Проверим замкнутость формы $\mathfrak{b}_\varepsilon$. Применяя неравенство Гёльдера и теорему вложения Соболева, можно показать (см. [10; (5.11)–(5.14)]), что для любого $\nu>0$ найдутся такие постоянные $C_j(\nu)>0$, что
Если $\nu$ фиксировано, то $C(\nu)$ зависит лишь от $d$, $\rho$, $\alpha_0$, от норм $\|g^{-1}\|_{L_\infty}$, $\|a_j\|_{L_\rho (\Omega)}$, $j=1,\dots,d$, и от параметров решетки $\Gamma$.
Таким образом, форма $\mathfrak{b}_\varepsilon$ замкнута и неотрицательна. Отвечающий ей самосопряженный в $L_2(\mathbb{R}^d;\mathbb{C}^n)$ оператор обозначим через $B_\varepsilon$. Формально можно написать
Эффективный оператор для $A_\varepsilon$ задается выражением $A^0=b(\mathbf{D})^*g^0b(\mathbf{D})$. Здесь $g^0$ – постоянная эффективная матрица размера $m\times m$. Матрица $g^0$ выражается через решение вспомогательной задачи на ячейке. Пусть $\Gamma$-периодическая $(n\times m)$-матрица-функция $\Lambda (\mathbf{x})$ – (слабое) решение задачи
где $\mathbf{l}_k(\mathbf{x})$, $k=1,\dots,m$, – столбцы матрицы $g(\mathbf{x})^{-1}$.
1.6. Эффективный оператор
Чтобы описать усреднение младших членов оператора $B_\varepsilon$, рассмотрим $\Gamma$-периодическую $(n\times n)$-матрицу-функцию $\widetilde{\Lambda}(\mathbf{x})$, являющуюся решением задачи
Здесь $c_*$ – постоянная (1.18). Постоянная $C_L$ зависит только от исходных данных (1.12). Отсюда вытекают оценки для квадратичной формы $\mathfrak{b}^0$ оператора (1.37):
который ограничен как оператор, действующий из $L_2(\mathbb{R}^d;\mathbb{C}^n)$ в $H^1(\mathbb{R}^d;\mathbb{C}^n)$. Это нетрудно установить с помощью предложения 1.2 и включений $\Lambda,\widetilde{\Lambda}\in \widetilde{H}^1(\Omega)$. Отметим, что $\|\varepsilon K(\varepsilon;\zeta)\|_{L_2 \to H^1}=O(1)$ при малом $\varepsilon$ и фиксированном $\zeta$.
В силу предложения 1.2 и включений $\widetilde{g}, g (b(\mathbf{D})\widetilde{\Lambda}) \in L_2(\Omega)$ оператор (1.42) ограничен из $L_2(\mathbb{R}^d;\mathbb{C}^n)$ в $L_2(\mathbb{R}^d;\mathbb{C}^m)$, причем $\| G(\varepsilon;\zeta)\|_{L_2\to L_2}=O(1)$.
Теорема 1.7 (см. [37]). Пусть выполнены условия теоремы 1.6. Пусть операторы $K(\varepsilon;\zeta)$ и $G(\varepsilon;\zeta)$ определены в (1.41) и (1.42) соответственно. Тогда при $0<\varepsilon \leqslant 1$ и $|\zeta| \geqslant 1$ справедливы оценки
Отметим, что условие (2.1) является более ограничительным, чем условие (1.4). Следующее утверждение установлено в книге [44] (см. теорему 7.8 из § 3.7; в этом утверждении границу $\partial \mathcal{O}$ достаточно считать липшицевой).
Предложение 2.2 (см. [44]). Пусть $\mathcal{O}\subset \mathbb{R}^d$ – ограниченная область с липшицевой границей. Условие 2.1 необходимо и достаточно для существования постоянных $k_1>0$, $k_2 \geqslant 0$ таких, что выполнено неравенство Гординга
Замечание 2.3. Постоянные $k_1$, $k_2$ зависят от матрицы $b(\boldsymbol{\xi})$ и от области $\mathcal{O}$, но в общем случае их трудно контролировать явно. Однако часто для конкретных операторов их можно найти. Поэтому в дальнейшем мы будем ссылаться на зависимость других постоянных от $k_1, k_2$.
2.2. Оператор $A_{N,\varepsilon}$
Далее предполагается, что $\mathcal{O}\,{\subset}\, \mathbb{R}^d$ – ограниченная область с границей класса $C^{1,1}$. В пространстве $L_2(\mathcal{O};\mathbb{C}^n)$ рассмотрим оператор $A_{N,\varepsilon}$, формально заданный дифференциальным выражением $b(\mathbf{D})^* g^\varepsilon(\mathbf{x}) b(\mathbf{D})$ при условии Неймана на границе. Зададим строгое определение: $A_{N,\varepsilon}$ есть самосопряженный оператор в $L_2(\mathcal{O};\mathbb{C}^n)$, порожденный квадратичной формой
Из (2.3)–(2.5) видно, что форма $\mathfrak{a}_{N,\varepsilon }$ замкнута и неотрицательна.
2.3. Оператор $B_{N,\varepsilon}$
Пусть $0< \varepsilon \leqslant 1$. Рассмотрим теперь более общий оператор $B_{N,\varepsilon}$, добавляя к $A_{N,\varepsilon}$ члены младшего порядка. Формально оператор $B_{N,\varepsilon}$ задан дифференциальным выражением
Здесь $\rho$ – показатель из условия (1.9), $q=\infty$ при $d=1$, $q=2\rho/(\rho -2)$ при $d\geqslant 2$. Покроем область $\mathcal{O}$ объединением ячеек решетки $\varepsilon\Gamma$, имеющих непустое пересечение с $\mathcal{O}$. Через $N_\varepsilon$ обозначим количество ячеек в этом покрытии. Ясно, что это объединение ячеек содержится в области $\widetilde{\mathcal{O}}$, представляющей собой $2r_1$-окрестность области $\mathcal{O}$, где $2r_1=\operatorname{diam}\Omega$. Поэтому количество ячеек $N_\varepsilon$ можно оценить сверху: $N_\varepsilon\leqslant \mathfrak{c}_1\varepsilon^{-d}$, где величина $\mathfrak{c}_1$ зависит только от области $\mathcal{O}$ и от параметров решетки $\Gamma$. Имеем
В силу компактности вложения $H^1(\mathcal{O};\mathbb{C}^n)\hookrightarrow L_q(\mathcal{O};\mathbb{C}^n)$ для любого $\mu >0$ существует постоянная $\check{\mathcal{C}}(\mu)>0$ такая, что
Если $\mu$ фиксировано, то $\check{\mathcal{C}}(\mu)$ зависит лишь от $d$, $\rho$ и от области $\mathcal{O}$. Из (2.7), (2.10) и (2.11) следует, что для любого $\nu>0$ найдется постоянная $\mathcal{C}(\nu)$ такая, что
Если $\nu$ фиксировано, то $\mathcal{C}(\nu)$ зависит лишь от $d$, $\rho$, $\|a_j\|_{L_\rho(\Omega)}$, $j=1,\dots,d$, от области $\mathcal{O}$ и параметров решетки $\Gamma$.
Аналогично, используя условие (1.10) и неравенство Гёльдера, имеем
где $\check{q} = \infty$ при $d=1$, $\check{q}= 2s/(s-1)$ при $d\geqslant 2$. Отсюда и из компактности вложения $H^1(\mathcal{O};\mathbb{C}^n)\hookrightarrow L_{\check{q}} (\mathcal{O};\mathbb{C}^n)$ следует, что для любого $\nu>0$ существует постоянная $\mathcal{C}_Q(\nu)>0$ такая, что
Если $\nu$ фиксировано, то $\mathcal{C}_Q(\nu)$ зависит лишь от $d$, $s$, $\|Q\|_{L_s(\Omega)}$, от области $\mathcal{O}$ и параметров решетки $\Gamma$.
Неравенства (2.17) и (2.18) показывают, что форма (2.6) замкнута и положительно определена. Самосопряженный оператор в $L_2(\mathcal{O};\mathbb{C}^n)$, порожденный этой формой, обозначаем через $B_{N,\varepsilon}$.
Нам понадобится вспомогательный оператор $\widetilde{B}_{N,\varepsilon}$. Факторизуем матрицу $Q_0$:
2.4. Оценки обобщенной резольвенты $(B_{N,\varepsilon}-\zeta Q_0^\varepsilon)^{-1}$
Наша цель – найти аппроксимацию обобщенной резольвенты $(B_{N,\varepsilon}-\zeta Q_0^\varepsilon )^{-1}$ с двухпараметрическими (по $\varepsilon$ и $\zeta$) оценками погрешности. Считаем, что $\zeta \in \mathbb{C}\setminus\mathbb{R}_+$. Иначе говоря, нас интересует поведение при малом $\varepsilon$ решения
Здесь $c(\phi)$ – величина (1.40). Постоянная $\mathcal{C}_1$ определяется равенством $\mathcal{C}_1^2 = c_4^{-1} \|Q_0^{-1}\|_{L_\infty}(1+\|Q_0\|_{L_\infty} \|Q_0^{-1}\|_{L_\infty})$.
Доказательство. В силу (2.19), (2.21) и неравенства $\widetilde{B}_{N,\varepsilon}> 0$ имеем
Чтобы проверить (2.24), подставим $\boldsymbol{\eta}= \mathbf{u}_\varepsilon$ в тождество (2.22) и воспользуемся оценкой (2.17) и уже доказанным неравенством (2.23). Получаем
где $c_6 = \max \{ 2 d \alpha_1 \|g\|_{L_\infty} +1,\, C_6 \}$.
Таким образом, в силу (2.32) и (2.33) форма (2.25) замкнута и положительно определена. Отвечающий ей самосопряженный оператор в $L_2(\mathcal{O};\mathbb{C}^n)$ обозначим через $B_N^0$.
Теперь мы окончательно фиксируем значение параметра $\lambda$:
где $\lambda_*$ определено в (1.16), $\lambda_0$ – в (2.16), и $\widehat{\lambda}$ – в (2.31). Это обеспечивает неотрицательность оператора $B_\varepsilon$, действующего в $L_2(\mathbb{R}^d; \mathbb{C}^n)$, и положительную определенность операторов $B_{N,\varepsilon}$ и $B_N^0$, действующих в $L_2(\mathcal{O};\mathbb{C}^n)$.
Для удобства дальнейших ссылок назовем “исходными данными” следующие величины:
Ясно, что постоянная (2.34), а также постоянные $c_4$, $c_5$, $c_6$ контролируются через исходные данные (2.35).
В силу условия $\partial\mathcal{O}\in C^{1,1}$ оператор $(B_N^0)^{-1}$ является непрерывным отображением из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^2(\mathcal{O};\mathbb{C}^n)$. При этом выполнена оценка
Здесь постоянная $\widehat{c}$ зависит лишь от исходных данных (2.35). Для оправдания этого факта сошлемся на теоремы о регулярности решений сильно эллиптических систем (см. [45; гл. 4]).
Замечание 2.5. Вместо условия $\partial\mathcal{O}\in C^{1,1}$ можно было бы наложить неявное требование: ограниченная область $\mathcal{O}\subset \mathbb{R}^d$ с липшицевой границей такова, что справедлива оценка (2.36). Для такой области результаты работы остаются в силе. В случае скалярных эллиптических операторов широкие достаточные условия на $\partial \mathcal{O}$, обеспечивающие справедливость оценки (2.36), можно найти в [46] и [47; гл. 7] (в частности, достаточно, чтобы $\partial\mathcal{O}\in C^\alpha$, $\alpha >3/2$).
Факторизуем $\overline{Q_0}=f_0^{-2}$. Отметим, что согласно (2.19)
Положим $(\partial\mathcal{O})_{\varepsilon} \,{:=}\,\{ \mathbf{x}\,{\in}\, \mathbb{R}^d \colon \!\operatorname{dist}\{ \mathbf{x};\partial\mathcal{O}\} \,{<}\,\varepsilon \}$. Выберем числа $\varepsilon_0,\varepsilon_1\in (0,1]$ согласно следующему условию.
Условие 3.1. Пусть число $\varepsilon_0\,{\in}\, (0,1]$ таково, что полоску $(\partial\mathcal{O})_{\varepsilon_0}$ можно покрыть конечным набором окрестностей, допускающих диффеоморфизмы класса $C^{0,1}$, распрямляющие границу $\partial\mathcal{O}$. Обозначим $\varepsilon_1=\varepsilon_0 (1+r_1)^{-1}$, где $2r_1=\operatorname{diam}\Omega$.
Ясно, что $\varepsilon_1$ зависит только от области $\mathcal{O}$ и от параметров решетки $\Gamma$. Отметим, что для справедливости условия 3.1 достаточно, чтобы $\partial\mathcal{O}$ была липшицевой. Мы наложили более ограничительное условие $\partial\mathcal{O}\in C^{1,1}$, чтобы гарантировать справедливость оценки (2.36).
Сформулируем основные результаты работы.
Теорема 3.2. Пусть $\mathcal{O}\subset \mathbb{R}^d$ – ограниченная область с границей класса $C^{1,1}$. Пусть $\zeta =| \zeta | e^{i\phi}\in \mathbb{C}\setminus \mathbb{R}_+$, $|\zeta| \geqslant 1$. Пусть $\mathbf{u}_\varepsilon=(B_{N,\varepsilon} -\zeta Q_0^\varepsilon)^{-1}\mathbf{F}$ и $\mathbf{u}_0=(B_N^0-\zeta \overline{Q_0}\,)^{-1}\mathbf{F}$, где $\mathbf{F} \in L_2(\mathcal{O};\mathbb{C}^n)$. Пусть число $\varepsilon_1$ подчинено условию 3.1. Тогда при $0<\varepsilon\leqslant \varepsilon_1$ справедлива оценка
Чтобы аппроксимировать решение в классе Соболева $H^1(\mathcal{O};\mathbb{C}^n)$, введем корректор. Для этого фиксируем линейный непрерывный оператор продолжения
где постоянная $C_{\mathcal O}^{(l)}$ зависит лишь от $l$ и от области ${\mathcal O}$. Через $R_\mathcal{O}$ обозначим оператор сужения функций в $\mathbb{R}^d$ на область $\mathcal{O}$. Положим
Лемма 3.3. Пусть операторы $K_N(\varepsilon;\zeta)$ и $G_N(\varepsilon;\zeta)$ определены в (3.5) и (3.6) соответственно. Оператор $K_N(\varepsilon;\zeta)$ непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathcal{O};\mathbb{C}^n)$, а оператор $G_N(\varepsilon;\zeta)$ непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $L_2(\mathcal{O};\mathbb{C}^m)$. При $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0< \varepsilon \leqslant 1$ выполнены оценки
Вместе с (2.40) и (2.41) с учетом ограничения $|\zeta|\geqslant 1$ это приводит к оценке (3.7) с постоянной $C_K'= M_1 \alpha_1^{1/2} C_\mathcal{O}^{(1)} \mathcal{C}_1 + \widetilde{M}_1 C_\mathcal{O}^{(0)} \|Q_0^{-1}\|_{L_\infty}$.
где $C_G' = 2 \|g\|_{L_\infty} \alpha_1^{1/2} C_\mathcal{O}^{(1)}$ и $C_G'' = |\Omega|^{-1/2} \alpha_0^{-1/2} C_a \|g\|^{1/2}_{L_\infty} \|g^{-1}\|^{1/2}_{L_\infty} C_\mathcal{O}^{(0)}$. Вместе с (2.40) и (2.41) это дает оценку (3.9) с постоянной $C_G = C_G' \mathcal{C}_1 + C_G'' \|Q_0^{-1}\|_{L_\infty}$. Лемма доказана.
Положим $\widetilde{\mathbf{u}}_0=P_\mathcal{O}\mathbf{u}_0$. Через $\mathbf{v}_\varepsilon$ обозначим первое приближение к решению $\mathbf{u}_\varepsilon$:
Иначе говоря, $\mathbf{v}_\varepsilon = (B_N^0-\zeta \overline{Q_0}\,)^{-1}\mathbf{F} +\varepsilon K_N(\varepsilon;\zeta )\mathbf{F}$, где $K_N(\varepsilon;\zeta)$ – оператор (3.5).
Теорема 3.4. Пусть выполнены условия теоремы 3.2. Пусть матрицы-функции $\Lambda (\mathbf{x})$ и $\widetilde{\Lambda}(\mathbf{x})$ – $\Gamma$’=периодические решения задач (1.20) и (1.30) соответственно, $S_\varepsilon $ – оператор сглаживания по Стеклову (1.1), и пусть $P_\mathcal{O}$ – оператор продолжения (3.3). Положим $\widetilde{\mathbf{u}}_0=P_\mathcal{O}\mathbf{u}_0$. Пусть функция $\mathbf{v}_\varepsilon$ определена в (3.11), (3.12). Тогда при $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0<\varepsilon \leqslant \varepsilon_1$ выполнена оценка
где $K_N(\varepsilon;\zeta)$ – оператор (3.5). Пусть матрица-функция $\widetilde{g}(\mathbf{x})$ определена в (1.22). Для потока $\mathbf{p}_\varepsilon:=g^\varepsilon b(\mathbf{D})\mathbf{u}_\varepsilon$ при $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0<\varepsilon \leqslant \varepsilon_1$ справедлива оценка
Здесь оператор $G_N(\varepsilon;\zeta)$ определен в (3.6). Постоянные $\mathcal{C}_4$, $\mathcal{C}_5$, $\widetilde{\mathcal{C}}_4$ и $\widetilde{\mathcal{C}}_5$ зависят только от исходных данных (2.35).
Следствие 3.5. В условиях теоремы 3.4 при $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0<\varepsilon \leqslant \varepsilon_1$ справедлива оценка
где $B^0$ – оператор (1.37). Тогда $\widetilde{\mathbf{F}}\in L_2(\mathbb{R}^d;\mathbb{C}^n)$ и $\widetilde{\mathbf{F}}|_{\mathcal{O}}=\mathbf{F}$. В силу верхней оценки (1.39), (3.18) и (3.20) имеем
т. е. $\widetilde{\mathbf{u}}_\varepsilon =(B_\varepsilon -\zeta Q_0^\varepsilon)^{-1}\widetilde{\mathbf{F}}$. Объединяя (3.21)–(3.23) и применяя теоремы 1.6 и 1.7, находим, что при $0<\varepsilon \leqslant 1$ и $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $| \zeta | \geqslant 1$, справедливы оценки
Легко видеть, что функционал (3.30) является антилинейным и непрерывным в $H^1(\mathcal{O};\mathbb{C}^n)$. Для первого члена непрерывность следует из предложения 1.2 и включений $\widetilde{g}, g(b(\mathbf{D}) \widetilde{\Lambda}) \in L_2(\Omega)$. Далее, в силу (2.10) и непрерывности вложения $H^1(\mathcal{O};\mathbb{C}^n)\hookrightarrow L_q(\mathcal{O};\mathbb{C}^n)$ (с константой вложения $C(q;\mathcal{O})$) выполнена оценка
где $C_{11}= \mathfrak{c}_1^{1 / \rho} C(q;\mathcal{O}) \bigl( \sum_{j=1}^d \|a_j\|^2_{L_\rho(\Omega)} \bigr)^{1/2}$. Отсюда следует, что второй член в (3.30) непрерывен в $H^1(\mathcal{O};\mathbb{C}^n)$. Аналогично из (2.13) и непрерывности вложения $H^1(\mathcal{O};\mathbb{C}^n)\hookrightarrow L_{\check{q}} (\mathcal{O};\mathbb{C}^n)$ (с константой вложения $C(\check{q};\mathcal{O})$) вытекает оценка
где $C_{12}= \mathfrak{c}_1^{1/2s} C(\check{q};\mathcal{O}) \| Q\|^{1/2}_{L_s(\Omega)}$. Поэтому третье слагаемое в (3.30) непрерывно в $H^1(\mathcal{O};\mathbb{C}^n)$. Для последних двух членов непрерывность очевидна.
Стандартным образом проверяется, что решение $\mathbf{w}_\varepsilon$ задачи (3.29) существует и единственно. Поправку $\mathbf{w}_\varepsilon$ такого типа часто называют “корректором типа пограничного слоя”.
Лемма 3.6. Пусть $\zeta \in \mathbb{C} \setminus \mathbb{R}_+$ и $|\zeta| \geqslant 1$. Пусть $\mathbf{u}_\varepsilon=(B_{N,\varepsilon} -\zeta Q_0^\varepsilon)^{-1}\mathbf{F}$, где $\mathbf{F} \in L_2(\mathcal{O};\mathbb{C}^n)$. Пусть $\widetilde{\mathbf{u}}_\varepsilon$ определено в п. 3.2. Пусть $\mathbf{w}_\varepsilon \in H^1(\mathcal{O};\mathbb{C}^n)$ удовлетворяет тождеству (3.29). Тогда справедливы оценки
Если $\operatorname{Re}\zeta <0$, то в равенстве (3.35) при $\boldsymbol{\eta}=\mathbf{V}_\varepsilon$ возьмем вещественную часть. При рассматриваемых значениях $\zeta$ имеем $c(\phi)=1$ и с учетом (3.36) получаем
С учетом (2.17) и (3.40) отсюда вытекает искомая оценка (3.33) с постоянной $\mathcal{C}_6^2= 81 C_{13}^2 c_4^{-2} + 18 C_{13}^2 c_4^{-1}\| Q_0^{-1}\|_{L_\infty}$. Соотношения (3.33) и (3.40) приводят к неравенству (3.34) с постоянной $\mathcal{C}_7^2= 4 C_{13} \mathcal{C}_6 \| Q_0^{-1}\|_{L_\infty}+ 4 C_{13}^2 \| Q_0^{-1}\|^2_{L_\infty}$. Лемма доказана.
Из леммы 3.6 и оценки (3.28) вытекает следующая теорема, которая показывает, что при учете поправки $\mathbf{w}_\varepsilon$ решение $\mathbf{u}_\varepsilon$ приближается функцией $\mathbf{v}_\varepsilon - \mathbf{w}_\varepsilon$ по норме в $H^1(\mathcal{O};\mathbb{C}^n)$ с погрешностью точного порядка $O(\varepsilon)$.
Теорема 3.7. Пусть выполнены условия теоремы 3.4. Пусть $\mathbf{w}_\varepsilon{\kern1pt}{\in}{\kern1pt}H^1(\mathcal{O};\mathbb{C}^n)$ удовлетворяет тождеству (3.29). Тогда при $0<\varepsilon\leqslant 1$ и $\zeta \in \mathbb{C}\setminus \mathbb{R}_+$, $| \zeta | \geqslant 1$, справедлива оценка
Постоянная $\mathcal{C}_8$ зависит только от исходных данных (2.35).
Дальнейший план доказательства теорем 3.2 и 3.4 таков. Мы сначала докажем оценку (3.14) при $\operatorname{Re} \zeta \leqslant 0$. Затем установим (3.2) также при $\operatorname{Re} \zeta \leqslant 0$, используя уже доказанную оценку (3.14) и соображения двойственности. После этого завершим доказательства теорем, опираясь на подходящие тождества для резольвент, позволяющие переносить уже доказанные оценки из точки $\zeta$ в левой полуплоскости в симметричную точку в правой полуплоскости. (Такой прием ранее применялся в [40; § 10].)
Поэтому для доказательства оценки (3.14) (при $\operatorname{Re} \zeta \leqslant 0$) нужно оценить надлежащим образом норму функции $\mathbf{w}_\varepsilon$ в $H^1(\mathcal{O};\mathbb{C}^n)$.
где $\mathcal{C}_9= \mathcal{C}_7+ C_1 C_{10}$. Поэтому доказательство теоремы 3.2 (при $\operatorname{Re} \zeta \leqslant 0$) сводится к оцениванию функции $\mathbf{w}_\varepsilon$ в $L_2(\mathcal{O};\mathbb{C}^n)$.
§ 4. Вспомогательные утверждения
4.1. Оценки в окрестности границы
В настоящем пункте приводятся вспомогательные утверждения, связанные с оценками интегралов по узкой окрестности $\partial \mathcal{O}$.
Лемма 4.1. Пусть справедливо условие 3.1. Обозначим $\Upsilon_\varepsilon =(\partial\mathcal{O})_\varepsilon\cap \mathcal{O}$. Тогда для любой функции $u\in H^1(\mathcal{O})$ справедлива оценка
Постоянная $\beta$ зависит только от области $\mathcal{O}$.
Лемма 4.2. Пусть выполнено условие 3.1. Пусть $f(\mathbf{x})$ – $\Gamma$-периодическая функция в $\mathbb{R}^d$ такая, что $f\in L_2(\Omega)$. Пусть $S_\varepsilon$ – оператор (1.1). Обозначим $\beta_*=\beta (1+r_1)$, где $2r_1=\operatorname{diam}\Omega$. Тогда при $0<\varepsilon\leqslant\varepsilon_1$ для любой функции $\mathbf{u}\in H^1(\mathbb{R}^d;\mathbb{C}^k)$ справедлива оценка
Лемма 4.2 аналогична лемме 2.6 из [18]. Леммы 4.1 и 4.2 проверены в [28; § 5] при условии $\partial \mathcal{O} \in C^1$, но доказательства переносятся и на случай условия 3.1.
4.2. Традиционная лемма теории усреднения
Нам понадобится следующий вариант традиционной леммы теории усреднения (см., например, [4; гл. 1, § 1]); доказательство этого варианта леммы можно найти в [5; лемма 3.1].
Лемма 4.3. Пусть $f_l(\mathbf{x})$, $l=1,\dots,d,$ – $\Gamma$-периодические $(n \times m)$-матрицы-функции в $\mathbb{R}^d$, причем
где последнее уравнение понимается в смысле теории распределений. Тогда существуют $\Gamma$-периодические $(n \times m)$-матрицы-функции $M_{lj}(\mathbf{x})$ в $\mathbb{R}^d$, $l,j=1,\dots,d$, такие, что
Здесь $C''$ зависит лишь от области $\mathcal{O}$ и решетки $\Gamma$, а $C'$ зависит от тех же параметров и от $\rho$.
$2^\circ$. При дополнительном предположении $a(\mathbf{x})=1$ (т. е. $h=p \in L_2(\Omega)$) при $0< \varepsilon \leqslant \varepsilon_1$ справедлива оценка
4.4. Свойства матриц-функций $\Lambda$ и $\widetilde{\Lambda}$
Следующий результат установлен в [28; следствие 2.4].
Лемма 4.5. Пусть $\Gamma$-периодическое решение $\Lambda (\mathbf{x})$ задачи (1.20) ограничено: $\Lambda\in L_\infty$. Тогда для любой функции $u\in H^1(\mathbb{R}^d)$ при $\varepsilon >0$ выполнена оценка
где $C(\widehat{q};\Omega)$ – норма оператора вложения $H^1(\Omega)\hookrightarrow L_{\widehat{q}}(\Omega)$. Здесь $\widehat{q}=\infty$ при $d=1$ и $\widehat{q}=2p(p-2)^{-1}$ при $d\geqslant 2$.
Следующий результат получен в [37; следствие 3.6].
Лемма 4.7. Пусть $\Gamma$-периодическое решение $\widetilde{\Lambda}(\mathbf{x})$ задачи (1.30) удовлетворяет условию (4.21). Тогда при любом $u\in H^2(\mathbb{R}^d)$ и $0<\varepsilon\leqslant 1$ справедлива оценка
Постоянные $\widetilde{\beta}_1$ и $\widetilde{\beta}_2$ зависят только от $n$, $d$, $\alpha_0$, $\alpha_1$, $\rho$, $\| g\|_{L_\infty}$, $\| g^{-1}\|_{L_\infty}$, от норм $\|a_j\|_{L_\rho (\Omega)}$, $j=1,\dots,d$, а также от параметров решетки $\Gamma$.
§ 5. Оценка поправки $\mathbf{w}_\varepsilon$ в $H^1(\mathcal{O})$ при $\operatorname{Re} \zeta \leqslant 0$
5.1. Случай $\operatorname{Re} \zeta \leqslant 0$. Оценка функционала $\mathcal{I}_\varepsilon[\boldsymbol{\eta}]$
где $\mathfrak{C}:= 2 \alpha_1^{1/2} \| g\|_{L_\infty}$. В силу (1.20)–(1.22) функции $f_l$, $l=1,\dots,d,$ удовлетворяют условиям леммы 4.3, а потому существуют $\Gamma$-периодические матрицы-функции $M_{lj}(\mathbf{x})$ в $\mathbb{R}^d$, $l,j=1,\dots,d$, подчиненные условиям (4.1)–(4.3). Из (4.3) и (5.13) следуют оценки
В силу (4.2) выполнено $f_l^\varepsilon(\mathbf{x}) = \varepsilon \sum_{j=1}^d \partial_j M_{lj}^\varepsilon(\mathbf{x})$, $l=1,\dots,d$. Отсюда и из (5.12) вытекает представление
что проверяется интегрированием по частям при учете равенства $M_{lj} = - M_{jl}$ (при проверке сначала считаем, что $\boldsymbol{\eta} \in H^2(\mathcal{O};\mathbb{C}^n)$, а затем замыкаем результат по непрерывности). Следовательно,
Обозначим последовательные слагаемые справа через $\boldsymbol{\psi}_l^{(1)}(\varepsilon)$, $\boldsymbol{\psi}_l^{(2)}(\varepsilon)$, $\boldsymbol{\psi}_l^{(3)}(\varepsilon)$. Тогда с учетом (4.4) имеем
где $C_{15} = \kappa (d \beta_* \alpha_1)^{1/2} r_0^{-1} \mathfrak{C} (C_8 C_9)^{1/2}$. Аналогично с помощью леммы 4.2 и (1.5), (3.19), (3.20), (4.4) и (5.13) приходим к оценке
где $\gamma_1^{(3)} = C_{14} d^{1/2}$ и $\gamma_2 = (C_{15}+ C_{16}) d^{1/2}$. Вместе с (5.8), (5.11), (5.15) и (5.18) это влечет искомую оценку (5.7) с постоянной $\gamma_1 = \gamma_1^{(1)} + \gamma_1^{(2)} + \gamma_1^{(3)}$. Лемма доказана.
Лемма 5.2. Пусть $\operatorname{Re} \zeta \leqslant 0$ и $|\zeta| \geqslant 1$. Пусть число $\varepsilon_1$ выбрано согласно условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ функционал (5.3) подчинен оценке
Начнем с оценки члена (5.25). В силу неравенства Гёльдера и непрерывности вложения $H^1(\Omega)\hookrightarrow L_q(\Omega)$ (с константой вложения $C(q;\Omega)$) имеем
Очевидно, $\check{f}_l(\mathbf{x})$ – $\Gamma$-периодические $(n \times n)$-матрицы-функции, причем $\check{f}_l \in L_2(\Omega)$ и в силу уравнения (1.30) выполнено $\sum_{l=1}^d D_l \check{f}_l(\mathbf{x})=0$. Заметим, что из (1.20) и (1.34) вытекает представление $V = - \overline{g (b(\mathbf{D}) \widetilde{\Lambda})}$. Поэтому $\int_\Omega \check{f}_l(\mathbf{x})\,d\mathbf{x} =0$, $l=1,\dots,d$. Таким образом, выполнены условия леммы 4.3, в силу которой существуют $\Gamma$-периодические $(n \times n)$-матрицы-функции $\check{M}_{lj}(\mathbf{x})$, $l,j=1,\dots,d$, удовлетворяющие соотношениям $\check{M}_{lj} \in \widetilde{H}^1(\Omega)$, $\int_\Omega \check{M}_{lj}(\mathbf{x})\,d\mathbf{x}=0$,
Следовательно, $\check{f}_l^\varepsilon(\mathbf{x}) = \varepsilon \sum_{j=1}^d \partial_j \check{M}^\varepsilon_{lj}(\mathbf{x})$, $l=1,\dots,d$. Отсюда и из (5.34) вытекает представление
где $C_{19,l} = \alpha_0^{-1/2} \alpha_1^{1/2} \|g\|^{1/2}_{L_\infty} \|g^{-1}\|^{1/2}_{L_\infty} C_a + \|a_l\|_{L_2(\Omega)}$. Отсюда и из (5.36) видно, что
что проверяется интегрированием по частям при учете равенства $\check{M}_{lj} = - \check{M}_{jl}$ (при проверке сначала считаем, что $\boldsymbol{\eta} \in H^2(\mathcal{O};\mathbb{C}^n)$, а затем замыкаем результат по непрерывности). Следовательно,
Обозначим последовательные слагаемые справа через $\check{\boldsymbol{\psi}}_l^{(1)}(\varepsilon)$, $\check{\boldsymbol{\psi}}_l^{(2)}(\varepsilon)$, $\check{\boldsymbol{\psi}}_l^{(3)}(\varepsilon)$. Тогда с учетом (4.4)
где $C_{20,l} = C_8 (2r_0)^{-1} |\Omega|^{-1/2} \bigl( \sum_{j=1}^d (C_{19,l} + C_{19,j})^2 \bigr)^{1/2}$. В силу (4.4) и леммы 4.2 при $0< \varepsilon \leqslant \varepsilon_1$ имеем
где $C_{21,l} = \kappa \beta_*^{1/2} (2 r_0)^{-1} |\Omega|^{-1/2} (C_7 C_8)^{1/2}\bigl( \sum_{j=1}^d (C_{19,l} + C_{19,j})^2 \bigr)^{1/2}$. Аналогично, применяя (4.4) и лемму 4.2, а затем (3.18), (3.19) и (5.40), приходим к оценке
Для оценки функционала (5.52) применим лемму 4.4($2^\circ$) при $a(\mathbf{x})=1$, $p(\mathbf{x})= a_l(\mathbf{x})$. Напомним обозначение $C_a = \bigl( \sum_{l=1}^d \|a_l\|^2_{L_2(\Omega)}\bigr)^{1/2}$. Получаем
Перейдем к рассмотрению члена (5.53). Из (1.30) и (1.34) вытекает представление $V^* = - \sum_{l=1}^d \overline{a_l (D_l \Lambda)}$. Поэтому член (5.53) можно записать в виде
где $\gamma_5^{(4)}= 2 C' \widetilde{C}_a |\Omega|^{1/2} M_2 \alpha_1^{1/2} C_9$, $\gamma_6^{(1)}= C'' |\Omega|^{1/2} M_2 \alpha_1^{1/2} (C_8C_9)^{1/2}$.
Остается рассмотреть член (5.54). Из (1.30) и (1.35) вытекает представление $W = - \sum_{l=1}^d \overline{a_l (D_l \widetilde{\Lambda})}$. Поэтому член (5.54) можно записать в виде
где $\gamma_7^{(2)}= (C_8 C_{23} \alpha_1^{1/2} + C_7 C_{24}) C_{12}$.
Наконец, для оценки члена (5.70) применим лемму 4.4($1^\circ$) при $a(\mathbf{x}) = |Q(\mathbf{x})|^{1/2}$, $p(\mathbf{x}) = |Q(\mathbf{x})|^{1/2}V(\mathbf{x})$ (подразумевается, что $Q(\mathbf{x})= |Q(\mathbf{x})| V(\mathbf{x})$ – полярное разложение матрицы $Q(\mathbf{x})$) и $\rho =2s$. Тогда при $0< \varepsilon \leqslant \varepsilon_1$ выполнена оценка
при $0< \varepsilon \leqslant \varepsilon_1$, где $\gamma_7^{(3)} = 2 {C}' C_8 \|Q\|^{1/2}_{L_s(\Omega)} \|Q\|^{1/2}_{L_1(\Omega)}$, $\gamma_8= {C}'' (C_7C_8)^{1/2} \|Q\|^{1/2}_{L_1(\Omega)}$.
В итоге, сопоставляя (5.67), (5.72), (5.75) и (5.76), приходим к искомой оценке (5.66) с постоянной $\gamma_7 = \gamma_7^{(1)}+ \gamma_7^{(2)} + \gamma_7^{(3)}$. Лемма доказана.
Лемма 5.5. Пусть $\operatorname{Re} \zeta \leqslant 0$ и $|\zeta| \geqslant 1$. Пусть число $\varepsilon_1$ выбрано согласно условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ функционал (5.6) подчинен оценке
с постоянными $\gamma_0' = \gamma_1 + \gamma_3 + \gamma_5 +\gamma_7 + \gamma_9$, $\gamma_0'' = \max\{\gamma_2 + \gamma_4, \gamma_6,\gamma_8\}$. Оценка (5.78) понадобится нам в следующем параграфе для доказательства теоремы 3.2. А для доказательства оценки (3.14) достаточно более грубой оценки
Постоянная $\mathcal{C}_{10}$ зависит лишь от данных задачи (2.35).
Доказательство. В тождество (3.29) подставим $\boldsymbol{\eta} = \mathbf{w}_\varepsilon$ и возьмем мнимую часть от полученного неравенства. С учетом (5.79) при $0< \varepsilon \leqslant \varepsilon_1$ получаем
Теперь возьмем вещественную часть от полученного равенства и, учитывая ограничение $\operatorname{Re} \zeta \leqslant 0$, при $0< \varepsilon \leqslant \varepsilon_1$ получаем
Постоянная $\mathcal{C}_{11}$ зависит лишь от данных задачи (2.35).
Доказательство. Рассмотрим тождество (3.29) и в качестве пробной функции $\boldsymbol{\eta}$ подставим $\boldsymbol{\eta}_\varepsilon = (B_{N,\varepsilon} -\zeta^* Q_0^\varepsilon)^{-1} \boldsymbol{\Phi}$, где $\boldsymbol{\Phi} \in L_2(\mathcal{O};\mathbb{C}^n)$. Тогда левая часть (3.29) совпадает с $(\mathbf{w}_\varepsilon, \boldsymbol{\Phi})_{L_2(\mathcal{O})}$, и тождество приобретает вид
где $\mathcal{C}_{11}^{(2)} \,{=}\, 2 C_9 r_1 \gamma_0 \,{+}\, 2 C_7 \gamma_{10}$. Для оценки первого слагаемого в (6.7) при $0\,{<}\, \varepsilon \,{\leqslant}\, \varepsilon_1$ применим оценку (5.78):
где $C_{25} = (\beta C_8 C_9)^{1/2}$. Далее, в силу леммы 4.2 и оценок (3.18), (3.19) (для $\widetilde{\boldsymbol{\eta}}_0$) с учетом ограничения $|\zeta| \geqslant 1$ при $0< \varepsilon \leqslant \varepsilon_1$ справедлива оценка
Отсюда, используя предложение 1.2 и лемму 4.2, а также (1.5), (1.24), (1.25), (1.32), (1.33) и оценки (3.18)–(3.20) (для $\widetilde{\boldsymbol{\eta}}_0$), с учетом $|\zeta| \geqslant 1$ при $0< \varepsilon \leqslant \varepsilon_1$ получаем
Применяя предложение 1.2, соотношения (1.5), (5.28), (5.29) и оценки (3.18), (3.19) (для $\widetilde{\boldsymbol{\eta}}_0$), с учетом $|\zeta| \geqslant 1$ получаем
где $C_{30}= \bigl(\sum_{l=1}^d C_{17,l}^2\bigr)^{1/2} \alpha_1^{1/2} C_8 +\bigl(\sum_{l=1}^d C_{18,l}^2\bigr)^{1/2} C_7$. Аналогично с учетом (5.73) и (5.74)
В результате приходим к оценке вида (3.1) с постоянной $\mathcal{C}'_3 = \max \{2 (\mathcal{C}_9 + \mathcal{C}_{11}); 2 \| Q_0^{-1}\|_{L_\infty} \}$. Итак, в операторной форме установлена оценка
где $\mathcal{C}_3^{(1)} = 4 \mathcal{C}_3' \|f\|^2_{L_\infty} \|f^{-1}\|^2_{L_\infty}$.
Рассмотрим теперь второе слагаемое в правой части (7.2). Пусть $\boldsymbol{\Phi}_1, \boldsymbol{\Phi}_2 \in L_2(\mathcal{O};\mathbb{C}^n)$. В силу леммы 4.4$(3^\circ)$ при $0< \varepsilon \leqslant \varepsilon_1$ справедлива оценка
где $\mathcal{C}_3^{(2)} = 4 \widetilde{C}''' \mathcal{C}_1 \|Q_0\|_{L_\infty} \|Q_0^{-1}\|_{L_\infty}$.
В результате из (7.2), (7.6) и (7.7) вытекает искомая оценка (3.2) с постоянной $\mathcal{C}_3 = \mathcal{C}_3^{(1)} + \mathcal{C}_3^{(2)}$ в случае $\zeta \in \mathbb{C} \setminus \mathbb{R}_+$, $\operatorname{Re} \zeta > 0$, $|\zeta| \geqslant 1$. С учетом уже установленной оценки (6.24) это завершает доказательство теоремы 3.2.
7.2. Случай $\operatorname{Re} \zeta > 0$. Доказательство оценки (3.14)
Запишем уже доказанную оценку (5.85) в точке $\widehat{\zeta}$:
Обозначим последовательные слагаемые в правой части (7.9) через $\mathfrak{L}_1(\varepsilon;\zeta)$, $\mathfrak{L}_2(\varepsilon;\zeta)$ и $\mathfrak{L}_3(\varepsilon;\zeta)$. Заметим, что $\mathfrak{L}_3(\varepsilon;\zeta) = \mathfrak{T}_2(\varepsilon;\zeta)$.
где $\mathcal{C}_5^{(3)} = 2 c_4^{-1/2} \widetilde{C}''' \|Q_0\|_{L_\infty}(\mathcal{C}_1 \|f\|_{L_\infty}+ c_4^{-1/2} \|Q_0^{-1}\|_{L_\infty})$.
В итоге из (7.9)–(7.11) и (7.16) вытекает оценка (3.14) с постоянной $\mathcal{C}_5 = \mathcal{C}_5^{(1)}+ \mathcal{C}_5^{(2)} + \mathcal{C}_5^{(3)}$ в случае $\zeta \in \mathbb{C} \setminus \mathbb{R}_+$, $\operatorname{Re} \zeta > 0$, $|\zeta| \geqslant 1$. С учетом (5.85) оценка (3.14) полностью доказана.
где $C_{34} = \|g\|_{L_\infty} (d \alpha_1)^{1/2} \mathcal{C}_1 + C_G$. Отсюда и из (3.16) вытекает, что при $\zeta \in \mathbb{C}\setminus \mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0 < \varepsilon \leqslant \varepsilon_1$ выполнено
Это влечет искомую оценку (3.17) с постоянной $\widetilde{\mathcal{C}}'_4 = \widetilde{\mathcal{C}}_4 + (C_{34} \widetilde{\mathcal{C}}_5)^{1/2}$.
§ 8. Устранение сглаживающего оператора. Специальные случаи. Оценки в строго внутренней подобласти
8.1. Устранение оператора $S_\varepsilon$ в корректоре
Покажем, что при некоторых дополнительных предположениях о решениях вспомогательных задач сглаживающий оператор $S_\varepsilon$ в корректоре может быть устранен (заменен тождественным оператором с сохранением порядка погрешности).
Условие 8.1. Предположим, что $\Gamma$-периодическое решение $\Lambda (\mathbf{x})$ задачи (1.20) ограничено, т. е. $\Lambda\in L_\infty (\mathbb{R}^d)$.
Условие 8.2. Предположим, что $\Gamma$-периодическое решение $\widetilde{\Lambda}(\mathbf{x})$ задачи (1.30) таково, что
Некоторые случаи, когда условия 8.1 и 8.2 выполнены автоматически, были выделены в [9; лемма 8.7] и [10; предложение 8.11] соответственно.
Предложение 8.3 (см. [9]). Условие 8.1 заведомо выполнено, если справедливо хотя бы одно из следующих предположений:
$1^\circ)$ $d\leqslant 2$;
$2^\circ)$ размерность $d$ произвольна, а оператор $A_\varepsilon$ имеет вид $A_\varepsilon =\mathbf{D}^* g^\varepsilon (\mathbf{x})\mathbf{D}$, где $g(\mathbf{x})$ – симметричная матрица с вещественными элементами;
$3^\circ)$ размерность $d$ произвольна, и $g^0=\underline{g}$, т. е. справедливы соотношения (1.29).
Предложение 8.4 (см. [10]). Условие 8.2 заведомо выполнено, если справедливо хотя бы одно из следующих предположений:
$1^\circ)$ $d\leqslant 4$;
$2^\circ)$ размерность $d$ произвольна, а оператор $A_\varepsilon$ имеет вид $A_\varepsilon =\mathbf{D}^* g^\varepsilon (\mathbf{x})\mathbf{D}$, где $g(\mathbf{x})$ – симметричная матрица с вещественными элементами.
Замечание 8.5. Если $A_\varepsilon =\mathbf{D}^* g^\varepsilon (\mathbf{x})\mathbf{D}$, где $g(\mathbf{x})$ – симметричная матрица с вещественными элементами, то из [50; гл. III, теорема 13.1] следует, что $\Lambda \in L_\infty$ и $\widetilde{\Lambda} \in L_\infty$. Тогда условия 8.1 и 8.2 выполнены автоматически. При этом норма $\|\Lambda\|_{L_\infty}$ не превосходит величины, зависящей от $d$, $\| g\|_{L_\infty}$, $\| g^{-1}\|_{L_\infty}$ и $\Omega$, а норма $\| \widetilde{\Lambda}\|_{L_\infty}$ оценивается в терминах $d$, $\rho$, $\| g\|_{L_\infty}$, $\| g^{-1}\|_{L_\infty}$, $\| a_j\|_{L_\rho (\Omega)}$, $j=1,\dots,d$, и $\Omega$.
Если выполнены условия 8.1 и 8.2, то оператор (8.1) непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathcal{O};\mathbb{C}^n)$, а оператор (8.2) непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $L_2(\mathcal{O};\mathbb{C}^m)$. Это легко проверить, используя леммы 4.5–4.7.
Наша цель в этом пункте – доказать следующую теорему.
Теорема 8.6. Пусть выполнены условия теоремы 3.4. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть операторы $K_N^0(\varepsilon;\zeta )$ и $G_N^0(\varepsilon;\zeta )$ определены в (8.1) и (8.2). Тогда при $\zeta \in \mathbb{C}\setminus \mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0<\varepsilon\leqslant\varepsilon_1$ справедливы оценки
Постоянные $\mathcal{C}_4$ и $\widetilde{\mathcal{C}}_4$ – те же, что в теореме 3.4. Постоянные $\mathcal{C}_{12}$ и $\widetilde{\mathcal{C}}_{12}$ зависят от исходных данных (2.35), а также от норм $\| \Lambda\|_{L_\infty}$ и $\|\widetilde{\Lambda}\|_{L_p(\Omega)}$.
Для доказательства теоремы 8.6 нам потребуются следующие утверждения (см. [40; леммы 7.7 и 7.8]).
Лемма 8.7. Пусть выполнено условие 8.1. Пусть $S_\varepsilon$ – оператор сглаживания по Стеклову (1.1). Тогда при $0<\varepsilon\leqslant 1$ выполнено
Постоянная $\mathfrak{C}_\Lambda$ зависит только от $m$, $d$, $\alpha_0$, $\alpha_1$, $\| g\|_{L_\infty}$, $\| g^{-1}\|_{L_\infty}$, от параметров решетки $\Gamma$ и нормы $\| \Lambda\|_{L_\infty}$.
Лемма 8.8. Пусть выполнено условие 8.2. Пусть $S_\varepsilon$ – оператор сглаживания по Стеклову (1.1). Тогда при $0<\varepsilon\leqslant 1$ справедлива оценка
Постоянная $\mathfrak{C}_{\widetilde{\Lambda}}$ зависит только от $n$, $d$, $\alpha_0$, $\alpha_1$, $\rho$, $\| g\|_{L_\infty}$, $\| g^{-1}\|_{L_\infty}$, от норм $\| a_j\|_{L_\rho (\Omega)}$, $j=1,\dots,d$, от $p$, $\| \widetilde{\Lambda}\|_{L_p(\Omega)}$ и от параметров решетки $\Gamma$.
Лемму 8.7 легко проверить, используя лемму 4.5, а для проверки леммы 8.8 нужны леммы 4.6 и 4.7.
где $C_{36} = (d\alpha_1)^{1/2}\| g\|_{L_\infty}\| \widetilde{\Lambda}\|_{L_p(\Omega)}C(\widehat{q};\Omega)C_\mathcal{O}^{(1)}\mathcal{C}_2$.
В итоге соотношения (8.7)–(8.10) вместе с (1.22) и (8.2) приводят к искомой оценке (8.4) с постоянной $\widetilde{\mathcal{C}}_{12}= (d\alpha_1)^{1/2}\| g\|_{L_\infty}\mathcal{C}_{12} + C_{35}+C_{36}$. Теорема доказана.
Замечание 8.9. Если выполнено только условие 8.1 (соответственно условие 8.2), то сглаживающий оператор $S_\varepsilon$ можно устранить только в члене корректора, содержащем $\Lambda^\varepsilon$ (соответственно $\widetilde{\Lambda}^\varepsilon$).
8.3. Случай, когда корректор обращается в нуль
Предположим, что $g^0=\overline{g}$, т. е. справедливы соотношения (1.28). Тогда $\Gamma$-периодическое решение задачи (1.20) обращается в нуль: $\Lambda (\mathbf{x})=0$. Пусть кроме этого выполнено равенство
Тогда $\Gamma$-периодическое решение задачи (1.30) также равно нулю: $\widetilde{\Lambda}(\mathbf{x})=0$. Поэтому в рассматриваемом случае оператор (3.5) обращается в нуль, формула (3.14) упрощается и из теоремы 3.4 вытекает следующий результат.
Предложение 8.10. Пусть справедливы равенства (1.28) и (8.11). Тогда при $0<\varepsilon\leqslant \varepsilon_1$ и $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, верна оценка
Предположим теперь, что $g^0=\underline{g}$, т. е. справедливы представления (1.29). Тогда в силу предложения 8.3($3^\circ$) выполнено условие 8.1. При этом согласно [8; замечание 3.5] матрица-функция (1.22) постоянна и совпадает с $g^0$, т. е. $\widetilde{g}(\mathbf{x})=g^0=\underline{g}$. Предположим дополнительно, что справедливо равенство (8.11). Тогда $\widetilde{\Lambda}(\mathbf{x})=0$ и из теоремы 8.6 вытекает следующий результат.
Предложение 8.11. Пусть имеют место соотношения (1.29) и (8.11). Тогда при $0<\varepsilon\leqslant\varepsilon_1$ и $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, верна оценка
Пусть $\mathcal{O}'$ – строго внутренняя подобласть области $\mathcal{O}$. Используя теорему 3.2 и результаты для задачи усреднения в $\mathbb{R}^d$, нетрудно получить аппроксимацию обобщенной резольвенты $(B_{N,\varepsilon}-\zeta Q_0^\varepsilon )^{-1}$ по $(L_2(\mathcal{O})\to H^1(\mathcal{O}'))$-норме с оценкой погрешности точного порядка $O(\varepsilon)$.
Теорема 8.12. Пусть выполнены условия теоремы 3.4. Пусть $\mathcal{O}'$ – строго внутренняя подобласть области $\mathcal{O}$. Введем обозначение $\delta :=\operatorname{dist}\{ \mathcal{O}';\partial \mathcal{O}\}$. Тогда при $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta|\geqslant 1$, и $0<\varepsilon\leqslant\varepsilon_1$ справедливы оценки
Постоянные $\mathcal{C}_{13}$, $\mathcal{C}_{14}$, $\widetilde{\mathcal{C}}_{13}$ и $\widetilde{\mathcal{C}}_{14}$ зависят только от исходных данных (2.35).
Доказательство. Начнем со случая, когда $\operatorname{Re} \zeta \leqslant 0$, $|\zeta| \geqslant 1$.
Фиксируем гладкую срезку $\chi(\mathbf{x})$ такую, что
Далее, функцию $\chi (\mathbf{u}_\varepsilon - \widetilde{\mathbf{u}}_\varepsilon)$ продолжим нулем на $\mathbb{R}^d \setminus \mathcal{O}$ и заметим, что $\mathfrak{U}(\varepsilon) =\mathfrak{b}_{\varepsilon} [\chi(\mathbf{u}_\varepsilon - \widetilde{\mathbf{u}}_\varepsilon), \chi(\mathbf{u}_\varepsilon - \widetilde{\mathbf{u}}_\varepsilon)]$. Тогда из (1.17) и (8.18) вытекает оценка
где $ \mathcal{C}_{14}' = \mathcal{C}_{14}^{(1)}+ \mathcal{C}_{14}^{(2)}$.
Рассмотрим теперь значения $\zeta \in \mathbb{C} \setminus \mathbb{R}_+$, $\operatorname{Re}\zeta > 0$, $|\zeta| \geqslant 1$. Пусть $\widehat{\zeta} = - \operatorname{Re}\zeta + i \operatorname{Im}\zeta$. Воспользуемся тождеством (7.9). Для второго и третьего слагаемых в правой части (7.9) применим прежние оценки (7.11) и (7.16). Оценим первое слагаемое в правой части (7.9):
При дополнительных предположениях о решениях вспомогательных задач сглаживатель $S_\varepsilon$ в корректоре удается устранить.
Теорема 8.13. Пусть выполнены условия теоремы 8.12. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть операторы $K_N^0(\varepsilon;\zeta)$ и $G_N^0(\varepsilon;\zeta)$ определены в (8.1), (8.2). Тогда при $0<\varepsilon\leqslant\varepsilon_1$ и $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$, $|\zeta| \geqslant 1$, справедливы оценки
Постоянные $\mathcal{C}_{13}$ и $\widetilde{\mathcal{C}}_{13}$ – те же, что в теореме 8.12. Постоянные $\mathcal{C}_{15}$ и $\widetilde{\mathcal{C}}_{15}$ зависят от исходных данных (2.35), а также от норм $\| \Lambda\|_{L_\infty}$, $\| \widetilde{\Lambda}\|_{L_p(\Omega)}$.
В теоремах из § 3 и § 8 предполагалось, что $\zeta\in\mathbb{C}\setminus\mathbb{R}_+$ и $|\zeta| \geqslant 1$. В настоящем параграфе устанавливаются результаты, справедливые в более широкой области изменения спектрального параметра. Материал этого параграфа аналогичен § 9 из [40].
9.1. Общий случай
Теорема 9.1. Пусть $\mathcal{O}\subset\mathbb{R}^d$ – ограниченная область с границей класса $C^{1,1}$. Пусть $c_\flat\geqslant 0$ – общая нижняя грань операторов $\widetilde{B}_{N,\varepsilon}=(f^\varepsilon )^*B_{N,\varepsilon}f^\varepsilon $ и $\widetilde{B}_N^0=f_0B_N^0f_0$. Пусть $\zeta\in\mathbb{C}\setminus [c_\flat,\infty)$. Положим $\psi =\operatorname{arg}(\zeta -c_\flat)$, $0<\psi <2\pi$. Введем обозначение
Пусть $\mathbf{v}_\varepsilon$ определено в (3.11), (3.12). Пусть $\mathbf{p}_\varepsilon = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon$. Пусть операторы $K_N(\varepsilon;\zeta)$ и $G_N(\varepsilon;\zeta)$ определены в (3.5) и (3.6). Пусть число $\varepsilon_1$ подчинено условию 3.1. Тогда при $0<\varepsilon \leqslant \varepsilon_1$ справедливы оценки
Постоянные $\mathfrak{C}_1$, $\mathfrak{C}_2$, $\mathfrak{C}_3$, $\widetilde{\mathfrak{C}}_2$ и $\widetilde{\mathfrak{C}}_3$ зависят только от исходных данных (2.35).
Замечание 9.2. 1) Величина $c(\psi)^2| \zeta - c_\flat |^{-2}$ в (9.1) обратна к квадрату расстояния от точки $\zeta$ до $[c_\flat,\infty)$.
2) В силу (2.17), (2.19), (2.32) и (2.37) в качестве $c_\flat$ можно выбрать постоянную $c_4 \|f^{-1}\|_{L_\infty}^{-2} = \frac{1}{2} k_1 \| g^{-1}\|^{-1}_{L_\infty} \| Q_0\|_{L_\infty}^{-1}$.
3) Легко оценить $c_\flat$ сверху. Из вариационных соображений и оценок (2.18) и (2.33) следует, что $c_\flat \leqslant \min\{c_5,c_6\} \|Q_0^{-1}\|_{L_\infty}$. Поэтому число $c_\flat$ контролируется через данные задачи (2.35).
Замечание 9.3. Оценки (9.2)–(9.7) выгодно применять при ограниченных значениях $|\zeta|$ и малом $\varepsilon \rho_\flat(\zeta)$. В этом случае величина $\varepsilon^{1/2} \rho_\flat(\zeta)^{1/2}+ \varepsilon | 1+ \zeta |^{1/2} \rho_ \flat (\zeta )$ контролируется через $\varepsilon^{1/2} \rho_\flat(\zeta)^{1/2}$. При больших значениях $|\zeta|$ оценки из теорем 3.2 и 3.4 предпочтительнее.
Начнем со следующих двух лемм, аналогичных леммам 9.4 и 9.5 из [40].
Лемма 9.4. В предположениях теоремы 9.1 при $0< \varepsilon \leqslant 1$ и $\zeta \in \mathbb{C} \setminus [c_\flat,\infty)$ справедливы оценки
Постоянные $\mathfrak{C}_4$ и $\mathfrak{C}_5$ зависят лишь от исходных данных (2.35).
Доказательство. При наших предположениях спектр оператора $\widetilde{B}_{N,\varepsilon}$ содержится в $[c_\flat,\infty)$, а потому $\|(\widetilde{B}_{N,\varepsilon} - \zeta I)^{-1}\|_{L_2(\mathcal{O}) \to L_2(\mathcal{O})} \leqslant c(\psi) |\zeta - c_\flat|^{-1}$. Вместе с (2.21) это влечет (9.8).
где $\mathfrak{C}_7' = \bigl( M_2 \alpha_1^{1/2} + \widetilde{M}_2 + \widetilde{M}_1\bigr) C_\mathcal{O}^{(1)} \mathfrak{C}_4$ и $\mathfrak{C}_7'' = M_1 \alpha_1^{1/2} C_\mathcal{O}^{(2)} \mathfrak{C}_5$.
В итоге из оценок (9.15) и (9.17) вытекает искомое неравенство (9.16) с постоянной $\mathfrak{C}_7 = \max \{ \mathfrak{C}_7'; \mathfrak{C}_6 + \mathfrak{C}_7''\}$. Лемма доказана.
Обозначим слагаемые в правой части (9.19) через $\mathfrak{J}_1(\varepsilon;\zeta)$ и $\mathfrak{J}_2(\varepsilon;\zeta)$. Аналогично (7.3) и (7.4) имеем
где $\mathfrak{C}_8 = \mathcal{C}'_3 \| f\|^2_{L_\infty}\| f^{-1}\|^2_{L_\infty}(c_\flat +2)^2$.
Оценим теперь норму оператора $\mathfrak{J}_2(\varepsilon;\zeta)$. Пусть $\boldsymbol{\Phi}_1, \boldsymbol{\Phi}_2 \in L_2(\mathcal{O};\mathbb{C}^n)$. В силу леммы 4.4$(3^\circ)$ при $0< \varepsilon \leqslant \varepsilon_1$ выполнено
Обозначим последовательные слагаемые в правой части (9.25) через $\mathcal{L}_1(\varepsilon;\zeta)$, $\mathcal{L}_2(\varepsilon;\zeta)$ и $\mathcal{L}_3(\varepsilon;\zeta)$. Заметим, что $\mathcal{L}_3(\varepsilon;\zeta) = {\mathfrak J}_2(\varepsilon;\zeta)$.
где $\mathfrak{C}_{13}$ зависит только от данных задачи (2.35). Из (9.31) и (9.32) вытекает оценка (9.4). Теорема доказана.
Следствие 9.6. В условиях теоремы 9.1 при $\zeta \in \mathbb{C} \setminus [c_\flat,\infty)$ и $0< \varepsilon \leqslant \varepsilon_1$ справедливы оценки
где $\mathfrak{C}_{15} = 2 \mathfrak{C}_4 + \mathfrak{C}_7$. При $|1+\zeta|^{1/2} \rho_\flat(\zeta)^{1/4} \leqslant \varepsilon^{-1/2}$ используем (9.6) и заметим, что $\varepsilon |1+\zeta|^{1/2} \rho_\flat(\zeta) \leqslant \varepsilon^{1/2} \rho_\flat(\zeta)^{3/4}$. При $|1+\zeta|^{1/2} \rho_\flat(\zeta)^{1/4} > \varepsilon^{-1/2}$ применим (9.35) и заметим, что $(1+ |\zeta|)^{-1/2} \rho_\flat(\zeta)^{1/2} < \varepsilon^{1/2} \rho_\flat(\zeta)^{3/4}$. В результате получаем оценку (9.33) с постоянной $\mathfrak{C}_{14}=\max \{ \mathfrak{C}_2+ \mathfrak{C}_3; 2 \mathfrak{C}_{15}\}$.
Соотношения (9.32) и (9.33) с учетом (1.3) и (1.6) приводят к оценке (9.34) с постоянной $\widetilde{\mathfrak{C}}_{14} = (d \alpha_1)^{1/2} \|g\|_{L_\infty} \mathfrak{C}_{14} + \mathfrak{C}_{13}$. Следствие доказано.
9.3. Устранение сглаживающего оператора
Теорема 9.7. Пусть выполнены условия теоремы 9.1. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть операторы $K_N^0(\varepsilon;\zeta)$ и $G_N^0(\varepsilon;\zeta)$ определены в (8.1), (8.2). Тогда при $0<\varepsilon\leqslant \varepsilon_1$ и $\zeta\in\mathbb{C}\setminus[c_\flat,\infty)$ справедливы оценки
Здесь постоянные $\mathfrak{C}_2$ и $\widetilde{\mathfrak{C}}_2$ – те же, что в теореме 9.1. Постоянные $\mathfrak{C}_{16}$ и $\widetilde{\mathfrak{C}}_{16}$ зависят только от исходных данных (2.35) и от норм $\| \Lambda\|_{L_\infty}$, $\|\widetilde{\Lambda}\|_{L_p(\Omega)}$.
Следующие утверждения проверяются аналогично предложениям 8.10 и 8.11.
Предложение 9.8. Пусть выполнены условия теоремы 9.1. Кроме того, пусть выполнены соотношения (1.28) и (8.11). Тогда при $\zeta\in\mathbb{C}\setminus [c_\flat,\infty)$ и $0<\varepsilon\leqslant \varepsilon_1$ справедлива оценка
Предложение 9.9. Пусть выполнены условия теоремы 9.1. Кроме того, пусть выполнены соотношения (1.29) и (8.11). Тогда при $\zeta\in\mathbb{C}\setminus [c_\flat,\infty)$ и $0<\varepsilon\leqslant \varepsilon_1$ справедлива оценка
Постоянные $\mathfrak{C}_{17}$, $\mathfrak{C}_{18}$, $\widetilde{\mathfrak{C}}_{17}$, $\widetilde{\mathfrak{C}}_{18}$ зависят только от исходных данных (2.35).
Доказательство. Запишем неравенство (8.12) в точке $\zeta=-1$:
где $\mathfrak{C}_{17} = \mathcal{C}_{13} \|f\|_{L_\infty} \|f^{-1}\|_{L_\infty} (c_\flat +2)$, $\mathfrak{C}_{18}' = \mathcal{C}_{14} \|f\|_{L_\infty} \|f^{-1}\|_{L_\infty} (c_\flat +2)$. Для второго и третьего членов в правой части (9.25) используем неравенства (9.27) и (9.30). В результате приходим к оценке (9.38) с постоянной $\mathfrak{C}_{18} = \mathfrak{C}_3 + \mathfrak{C}_{18}'$.
При дополнительных предположениях о решениях вспомогательных задач справедливо следующее утверждение.
Теорема 9.11. Пусть выполнены условия теоремы 9.10. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть операторы $K_N^0(\varepsilon;\zeta)$ и $G_N^0(\varepsilon;\zeta)$ определены в (8.1), (8.2). Тогда при $\zeta \in \mathbb{C} \setminus [c_\flat,\infty)$ и $0 < \varepsilon \leqslant \varepsilon_1$ справедливы оценки
Постоянные $\mathfrak{C}_{17}$ и $\widetilde{\mathfrak{C}}_{17}$ – те же, что в теореме 9.10. Постоянные $\mathfrak{C}_{19}$ и $\widetilde{\mathfrak{C}}_{19}$ зависят от исходных данных (2.35), а также от норм $\| \Lambda\|_{L_\infty}$, $\| \widetilde{\Lambda}\|_{L_p(\Omega)}$.
§ 10. Усреднение решений второй начально-краевой задачи для параболических уравнений
В этом параграфе мы изучаем усреднение решений второй начально-краевой задачи для параболического уравнения (т. е. задачи с условием Неймана). Результаты выводятся из аппроксимаций обобщенной резольвенты с двухпараметрическими оценками погрешности методом интегрирования резольвенты по контуру. Доказательства основных результатов этого параграфа полностью аналогичны доказательствам из статьи [42], где рассматривалась первая начально-краевая задача (с условием Дирихле). Поэтому мы ограничимся формулировками и краткими комментариями, опуская детали.
10.1. Постановка задачи
Рассмотрим (слабое) решение $\mathbf{u}_\varepsilon(\mathbf{x},t)$ задачи
при естественном условии (условии Неймана) на $\partial \mathcal{O} \times \mathbb{R}_+$. Здесь $B_\varepsilon$ – дифференциальное выражение (1.19), коэффициенты которого подчинены предположениям § 1. Предполагается, что $\boldsymbol \varphi \in L_2(\mathcal{O};\mathbb{C}^n)$. Используя (2.19), несложно убедиться, что справедливо представление
Нас интересует поведение решения $\mathbf{u}_\varepsilon(\mathbf{x},t)$ при малом $\varepsilon$, т. е. поведение окаймленной операторной экспоненты $f^\varepsilon e^{-\widetilde{B}_{N,\varepsilon} t} (f^\varepsilon)^*$.
Как в § 9, пусть $c_\flat > 0$ – общая нижняя грань операторов $\widetilde{B}_{N,\varepsilon}=(f^\varepsilon )^*B_{N,\varepsilon}f^\varepsilon $ и $\widetilde{B}_N^0=f_0B_N^0f_0$. Согласно замечанию 9.2 можно фиксировать выбор $c_\flat$ следующим образом:
Здесь постоянные $c_4$, $\widehat{c}$ и $c_\flat$ – те же, что в (2.17), (2.36) и (10.3) соответственно.
10.2. Аппроксимация решения в $L_2(\mathcal{O};\mathbb{C}^n)$
Следующий результат выводится из теорем 3.2 и 9.1 по аналогии с доказательством теоремы 2.2 из [42].
Теорема 10.2. Пусть $\mathbf{u}_\varepsilon(\mathbf{x},t)$ – решение задачи (10.1) с условием Неймана на $\partial \mathcal{O} \times \mathbb{R}_+$. Пусть $\mathbf{u}_0(\mathbf{x},t)$ – решение задачи (10.2) с условием Неймана на $\partial \mathcal{O} \times \mathbb{R}_+$. Пусть число $\varepsilon_1$ подчинено условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ справедлива оценка
Здесь $\gamma\subset\mathbb{C}$ – контур, обходящий спектр оператора $\widetilde{B}_{N,\varepsilon}$ в положительном направлении. Для экспоненты от оператора $\widetilde{B}_N^0$ справедливо аналогичное представление. Так как постоянная (10.3) – общая нижняя грань операторов $\widetilde{B}_{N,\varepsilon}$ и $\widetilde{B}_N^0$, в качестве контура интегрирования удобно выбрать
При $t>0$ оператор (10.12) непрерывно переводит $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathcal{O};\mathbb{C}^n)$. Действительно, согласно (10.6) при $t>0$ оператор $f_0e^{-\widetilde{B}_N^0t}f_0$ непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^2(\mathcal{O};\mathbb{C}^n)$. Следовательно, оператор $b(\mathbf{D})P_\mathcal{O}f_0e^{-\widetilde{B}_N^0 t}f_0$ непрерывно переводит $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathbb{R}^d;\mathbb{C}^m)$, а оператор $P_\mathcal{O}f_0 e^{-\widetilde{B}_N^0 t}f_0$ заведомо непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^1(\mathbb{R}^d;\mathbb{C}^n)$. Остается учесть непрерывность операторов $[\Lambda^\varepsilon ]S_\varepsilon\colon H^1(\mathbb{R}^d;\mathbb{C}^m)\to H^1(\mathbb{R}^d;\mathbb{C}^n)$ и $[\widetilde{\Lambda}^\varepsilon ]S_\varepsilon \colon H^1(\mathbb{R}^d;\mathbb{C}^n) \to H^1(\mathbb{R}^d;\mathbb{C}^n)$, вытекающую из предложения 1.2 и включений $\Lambda,\widetilde{\Lambda}\in \widetilde{H}^1(\Omega)$. Аналогично проверяется непрерывность оператора (10.13) из $L_2(\mathcal{O};\mathbb{C}^n)$ в $L_2(\mathcal{O};\mathbb{C}^m)$.
Положим $\widetilde{\mathbf{u}}_0(\,{\cdot}\,,t):=P_\mathcal{O}\mathbf{u}_0(\,{\cdot}\,,t)$. Через $\mathbf{v}_\varepsilon(\mathbf{x},t)$ обозначим первое приближение к решению $\mathbf{u}_\varepsilon(\mathbf{x},t)$ задачи (10.1):
т. е. $\mathbf{v}_\varepsilon (\,{\cdot}\,,t) = f_0 e^{-\widetilde{B}_N^0t}f_0\boldsymbol{\varphi}(\,{\cdot}\,)+\varepsilon \mathcal{K}_N(t;\varepsilon )\boldsymbol{\varphi}(\,{\cdot}\,)$.
Теорема 10.3. Пусть выполнены условия теоремы 10.2. Пусть функция $\mathbf{v}_\varepsilon(\mathbf{x},t)$ определена в (10.14). Положим $\mathbf{p}_\varepsilon(\mathbf{x},t):= g^\varepsilon(\mathbf{x}) b(\mathbf{D}) \mathbf{u}_\varepsilon(\mathbf{x},t)$. Пусть операторы $\mathcal{K}_N(t;\varepsilon)$ и $\mathcal{G}_N(t;\varepsilon)$ определены в (10.12) и (10.13) соответственно. Тогда при $0<\varepsilon\leqslant\varepsilon_1$ и $t>0$ выполнено
Здесь $K_N(\varepsilon;\zeta )$ и $G_N(\varepsilon;\zeta )$ – операторы (3.5) и (3.6) соответственно.
С помощью тождества (10.17), применяя теорему 9.1 при $\zeta \in \gamma$, $|\zeta| \leqslant \check{c}$, и теорему 3.4 при $\zeta \in \gamma$, $|\zeta| > \check{c}$, легко вывести оценку (10.15). Аналогично на основании тождества (10.18), используя (9.7) при $\zeta \in \gamma$, $|\zeta| \leqslant \check{c}$, и (3.17) при $\zeta \in \gamma$, $|\zeta| > \check{c}$, приходим к (10.16). Теорема доказана.
Замечание 10.4. Пусть $\lambda_1^0$ – первое собственное значение оператора $B_N^0$, и пусть $\sigma>0$ – сколь угодно малое число. Очевидно, число $\lambda_1^0 \| Q_0\|^{-1}_{L_\infty}$ является нижней гранью оператора $\widetilde{B}_N^0$. В силу резольвентной сходимости $B_{N,\varepsilon}$ к $B_N^0$ при достаточно малом $\varepsilon_\circ$ число $\lambda_1^0 \| Q_0\|^{-1}_{L_\infty} - \sigma /2$ является общей нижней гранью операторов $\widetilde{B}_{N,\varepsilon}$ при всех $0< \varepsilon \leqslant \varepsilon_\circ$. Поэтому можно сдвинуть контур интегрирования $\gamma$ так, чтобы он пересекал вещественную ось в точке $\mathfrak{c}_\circ :=\lambda_1^0 \| Q_0\|^{-1}_{L_\infty} - \sigma$ вместо $c_\flat/2$. На этом пути получаются оценки вида (10.7), (10.15), (10.16) с заменой $e^{-c_\flat t/2}$ на $e^{-\mathfrak{c}_\circ t}$ в правых частях. При этом постоянные в оценках станут зависеть от $\sigma$.
10.4. Устранение сглаживателя $S_\varepsilon$ в корректоре
Сглаживающий оператор в корректоре удается устранить, если решения вспомогательных задач подчинены дополнительным условиям. Следующий результат проверяется аналогично теореме 10.3 с помощью теорем 8.6 и 9.7.
Теорема 10.5. Пусть выполнены условия теоремы 10.3. Кроме того, пусть выполнены условия 8.1 и 8.2. Положим
Постоянные $\mathrm{C}_3$ и $\widetilde{\mathrm{C}}_3$ зависят от исходных данных (2.35), от $p$ и от норм $\| \Lambda\|_{L_\infty}$ и $\|\widetilde{\Lambda}\|_{L_p(\Omega)}$.
Замечание 10.6. Если выполнено только условие 8.1 (соответственно условие 8.2), то сглаживающий оператор $S_\varepsilon$ может быть устранен в члене корректора, содержащем $\Lambda^\varepsilon$ (соответственно $\widetilde{\Lambda}^\varepsilon$).
Устранить сглаживатель $S_\varepsilon$ в корректоре возможно также за счет усиления предположения о гладкости границы. Рассмотрим случай $d\geqslant 3$, поскольку при $d\leqslant 2$ применима теорема 10.5 (см. предложения 8.3 и 8.4).
Теорема 10.7. Пусть выполнены условия теоремы 10.2, причем $d\geqslant 3$. Пусть область $\mathcal{O}\subset \mathbb{R}^d$ – ограниченная область с границей класса $C^{d/2,1}$ если $d$ – четное, и класса $C^{(d+1)/2,1}$, если $d$ – нечетное. Пусть операторы $\mathcal{K}_N^0(t;\varepsilon)$ и $\mathcal{G}_N^0(t;\varepsilon)$ определены в (10.19) и (10.20) соответственно. Тогда при $t>0$ и $0<\varepsilon\leqslant\varepsilon_1$ справедливы оценки
Постоянные $\mathrm{C}_4(d)$ и $\widetilde{\mathrm{C}}_4(d)$ зависят только от исходных данных (2.35).
Теорема 10.7 опирается на следующую лемму, аналогичную лемме 2.8 из [42].
Лемма 10.8. Пусть $k \geqslant 2$ – целое число. Пусть $\mathcal{O}\subset\mathbb{R}^d$ – ограниченная область с границей $\partial\mathcal{O}$ класса $C^{k-1,1}$. Тогда при $t>0$ оператор $e^{-\widetilde{B}_N^0t}$ непрерывно переводит $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^q(\mathcal{O};\mathbb{C}^n)$, $0\leqslant q\leqslant k$, и выполнена оценка
Постоянная $\widehat{\mathrm{C}}_q$ зависит только от $q$ и исходных данных (2.35).
Доказательство леммы 10.8 и теоремы 10.7 полностью аналогично случаю первой начально-краевой задачи, рассмотренному в [42; п. 2.7 и § 7]. Ясно, что теорему 10.7 удобно применять при $t$ отделенном от нуля. При малых значениях $t$ порядок множителя $(\varepsilon^{1/2}t^{-3/4}+\varepsilon t^{-d/4-1/2})$ растет с ростом размерности. Это компенсация за устранение сглаживателя.
10.5. Специальные случаи
Выделим специальные случаи. Предположим, что $g^0=\overline{g}$, т. е. выполнены соотношения (1.28). Кроме того, пусть справедливо условие (8.11). Тогда $\Gamma$-периодические решения задач (1.20) и (1.30) равны нулю: $\Lambda (\mathbf{x})=0$ и $\widetilde{\Lambda}(\mathbf{x})=0$. Теорема 10.3 приводит к следующему результату.
Предложение 10.9. Пусть справедливы соотношения (1.28) и (8.11). Тогда в условиях теоремы 10.2 при $t>0$ и $0<\varepsilon\leqslant \varepsilon_1$ справедлива оценка
Предположим теперь, что $g^0=\underline{g}$, т. е. справедливы представления (1.29). Тогда в силу предложения 8.3, $3^\circ$ выполнено условие 8.1. При этом $\widetilde{g}(\mathbf{x})=g^0=\underline{g}$. Предположим дополнительно, что справедливо равенство (8.11). Тогда $\widetilde{\Lambda}(\mathbf{x})=0$ и из теоремы 10.3 с помощью предложения 1.1 несложно вывести следующий результат (ср. доказательство предложения 2.13 из [42]).
Предложение 10.10. Пусть имеют место соотношения (1.29) и (8.11). Тогда в условиях теоремы 10.2 при $0<\varepsilon\leqslant\varepsilon_1$ и $t>0$ верна оценка
Постоянная $\widetilde{\mathrm{C}}'_3$ зависит только от исходных данных (2.35).
10.6. Оценки в строго внутренней подобласти
Следующий результат легко вывести, применяя теоремы 8.12, 9.10 и тождества (10.17), (10.18).
Теорема 10.11. Пусть выполнены условия теоремы 10.3. Пусть $\mathcal{O}'$ – строго внутренняя подобласть области $\mathcal{O}$ и $\delta = \operatorname{dist} \{\mathcal{O}'; \partial \mathcal{O} \}$. Тогда при $0<\varepsilon\leqslant\varepsilon_1$ и $t>0$ выполнены оценки
Постоянные $\mathrm{C}_5$, $\mathrm{C}_6$, $\widetilde{\mathrm{C}}_5$ и $\widetilde{\mathrm{C}}_6$ зависят только от исходных данных (2.35).
Следующий результат проверяется на основании теорем 8.13, 9.11 и тождеств (10.17), (10.18).
Теорема 10.12. Пусть выполнены условия теоремы 10.11. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть операторы $\mathcal{K}_N^0(t;\varepsilon)$ и $\mathcal{G}_N^0(t;\varepsilon)$ определены в (10.19) и (10.20) соответственно. Тогда при $t>0$ и $0<\varepsilon\leqslant \varepsilon_1$ справедливы оценки
Постоянные $\mathrm{C}_5$ и $\widetilde{\mathrm{C}}_5$ – те же, что в теореме 10.11. Постоянные $\mathrm{C}_7$ и $\widetilde{\mathrm{C}}_7$ зависят от исходных данных (2.35), от $p$ и от норм $\| \Lambda\|_{L_\infty}$, $\|\widetilde{\Lambda}\|_{L_p(\Omega)}$.
Отметим, что устранить сглаживатель $S_\varepsilon$ в корректоре в оценках из теоремы 10.11 возможно и без дополнительных условий на матрицы-функции $\Lambda$ и $\widetilde{\Lambda}$. Рассмотрим случай $d\geqslant 3$ (иначе в силу предложений 8.3 и 8.4 применима теорема 10.12). При $t>0$ оператор $e^{-\widetilde{B}_N^0t}$ непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^2(\mathcal{O};\mathbb{C}^n)$ и выполнена оценка (10.6). Кроме того, справедливо свойство “повышения гладкости” внутри области: при $t>0$ оператор $e^{-\widetilde{B}_N^0t}$ непрерывен из $L_2(\mathcal{O};\mathbb{C}^n)$ в $H^l (\mathcal{O}';\mathbb{C}^n)$ при любом целом $l \geqslant 3$. При этом имеют место неравенства
Постоянная $\mathrm{C}'_l$ зависит от $l$ и исходных данных (2.35) (ср. оценку (2.45) из [42] и комментарии к ней).
Используя (10.21), а также свойства матриц-функций $\Lambda$, $\widetilde{\Lambda}$ и оператора $S_\varepsilon$, из теоремы 10.11 можно вывести следующий результат. Доказательство полностью аналогично случаю первой начально-краевой задачи (см. [42; п. 2.10 и § 8]).
Теорема 10.13. Пусть выполнены условия теоремы 10.11, причем $d\geqslant 3$. Пусть операторы $\mathcal{K}_D^0(t;\varepsilon)$ и $\mathcal{G}_D^0(t;\varepsilon)$ определены в (10.19) и (10.20) соответственно. Обозначим
при условии Неймана на $\partial \mathcal{O} \times (0,T)$. Предполагается, что $0< T \leqslant \infty$, $\mathbf{F} \in \mathfrak{H}_r(T):= L_r((0,T); L_2(\mathcal{O};\mathbb{C}^n))$ при некотором $1\leqslant r \leqslant \infty$; $\boldsymbol{\varphi} \in L_2(\mathcal{O};\mathbb{C}^n)$. Тогда
Оценивая интегральный член при $1< r \leqslant \infty$, получаем следующий результат (ср. доказательство теоремы 5.1 из [35]).
Теорема 10.14. Предположим, что $\boldsymbol{\varphi} \in L_2(\mathcal{O};\mathbb{C}^n)$ и $\mathbf{F} \in \mathfrak{H}_r(T)$ при некотором $1< r \,{\leqslant}\,\infty$, где $0< T \leqslant \infty$. Пусть $\mathbf{u}_\varepsilon(\mathbf{x},t)$ – решение задачи (10.22) с условием Неймана на $\partial \mathcal{O} \times (0,T)$. Пусть $\mathbf{u}_0(\mathbf{x},t)$ – решение задачи (10.24) с условием Неймана на $\partial \mathcal{O} \times (0,T)$. Пусть число $\varepsilon_1$ подчинено условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ и $0< t <T$ справедлива оценка
$c_\flat$ – постоянная (10.3). Постоянная $\mathrm{C}_1$ зависит только от исходных данных (2.35), а $c_r$ зависит от $r$ и данных (2.35).
По аналогии с доказательством теоремы 5.2 из [35], несложно вывести из теоремы 10.2 оценку нормы разности $\mathbf{u}_\varepsilon - \mathbf{u}_0$ в классе $\mathfrak{H}_r(T)$.
Теорема 10.15. Предположим, что $\boldsymbol{\varphi} \in L_2(\mathcal{O};\mathbb{C}^n)$ и $\mathbf{F} \in \mathfrak{H}_r(T)$ при некотором $1 \leqslant r \,{<}\, \infty$, где $0< T \leqslant \infty$. Пусть $\mathbf{u}_\varepsilon(\mathbf{x},t)$ – решение задачи (10.22) с условием Неймана на $\partial \mathcal{O} \times (0,T)$. Пусть $\mathbf{u}_0(\mathbf{x},t)$ – решение задачи (10.24) с условием Неймана на $\partial \mathcal{O} \times (0,T)$. Пусть число $\varepsilon_1$ подчинено условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ справедлива оценка
Здесь $\vartheta(\varepsilon,\cdot)$ определено в (10.26), $r^{-1}+ (r')^{-1}=1$. Постоянная $\mathrm{C}_9$ зависит только от исходных данных (2.35), а $c_{r'}$ зависит от $r$ и данных (2.35).
Замечание 10.16. В случае, когда $\boldsymbol{\varphi}=0$ и $\mathbf{F} \in \mathfrak{H}_\infty(T)$, из теоремы 10.14 вытекает оценка
Получим теперь аппроксимацию решения задачи (10.22) в $H^1(\mathcal{O};\mathbb{C}^n)$ с помощью теоремы 10.3. Трудности возникают при рассмотрении интегрального члена в (10.23) из-за сингулярности правой части оценки (10.15) при малом $t$. Считая, что $t \geqslant \varepsilon^2$, мы разбиваем промежуток интегрирования в (10.23) на две части: $(0, t-\varepsilon^2)$ и $(t-\varepsilon^2,t)$. На интервале $(0, t-\varepsilon^2)$ применяем (10.15), а на $(t-\varepsilon^2,t)$ используем оценку
Тогда $\mathbf{u}_{0,\varepsilon}(\,{\cdot}\,,t) = f_0 e^{- \widetilde{B}_{N}^0 \varepsilon^2} f_0^{-1} \mathbf{u}_0(\,{\cdot}\,, t-\varepsilon^2)$. Положим $\widetilde{\mathbf{u}}_{0,\varepsilon}(\,{\cdot}\,,t) = P_\mathcal{O} \mathbf{u}_{0,\varepsilon}(\,{\cdot}\,,t)$, где $P_\mathcal{O}$ – оператор продолжения (3.3). В качестве первого приближения к решению задачи (10.22) возьмем
Теорема 10.17. Предположим, что $\boldsymbol{\varphi} \in L_2(\mathcal{O};\mathbb{C}^n)$ и $\mathbf{F} \in \mathfrak{H}_r(T)$ при некотором $2< r \leqslant \infty$, где $0< T \leqslant \infty$. Пусть $\mathbf{u}_\varepsilon(\mathbf{x},t)$ – решение задачи (10.22) с условием Неймана на $\partial \mathcal{O} \times (0,T)$, $\mathbf{u}_0(\mathbf{x},t)$ – решение задачи (10.24) с условием Неймана на $\partial \mathcal{O} \times (0,T)$ и $\mathbf{p}_\varepsilon(\,{\cdot}\,,t) = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon(\,{\cdot}\,,t)$. Пусть функции $\mathbf{v}_{\varepsilon}$ и $\mathbf{q}_\varepsilon$ определены в (10.27), (10.28) соответственно. Пусть число $\varepsilon_1$ подчинено условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ и $\varepsilon^2 \leqslant t < T$ справедливы оценки
$c_\flat$ – постоянная (10.3). Постоянные $\mathrm{C}_2$, $\widetilde{\mathrm{C}}_2$ зависят только от исходных данных (2.35), а $\check{c}_r$ и $\widetilde{c}_r$ зависят от $r$ и данных (2.35).
Поскольку правая часть в оценке (10.16) имеет меньшую сингулярность при малом $t$, нежели правая часть в (10.15), при $r>4$ удается аппроксимировать поток $\mathbf{p}_\varepsilon(\,{\cdot}\,,t)$ функцией
Следующее утверждение проверяется аналогично доказательству предложения 3.5 из [42].
Предложение 10.18. Пусть $\boldsymbol{\varphi} \in L_2(\mathcal{O};\mathbb{C}^n)$ и $\mathbf{F} \in \mathfrak{H}_r(T)$ при некотором $4< r \leqslant \infty$, где $0< T \leqslant \infty$. Пусть $\mathbf{u}_\varepsilon(\mathbf{x},t)$ – решение задачи (10.22) с условием Неймана на $\partial \mathcal{O} \times (0,T)$ и $\mathbf{p}_\varepsilon(\,{\cdot}\,,t) = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon(\,{\cdot}\,,t)$. Пусть функция $\mathbf{q}^0_\varepsilon$ определена в (10.29). Пусть число $\varepsilon_1$ подчинено условию 3.1. При $0< \varepsilon \leqslant \varepsilon_1$ и $0 < t < T$ справедлива оценка
Здесь $c_\flat$ – постоянная (10.3). Постоянная $\widetilde{\mathrm{C}}_2$ зависит только от исходных данных (2.35), а $\widehat{c}_r$ зависит от $r$ и данных (2.35).
При дополнительных предположениях о решениях вспомогательных задач из теоремы 10.5 несложно вывести следующий результат.
Теорема 10.19. Пусть выполнены условия теоремы 10.17. Кроме того, пусть выполнены условия 8.1 и 8.2. Обозначим
Постоянные $\mathrm{C}_3$, $\widetilde{\mathrm{C}}_3$ зависят только от исходных данных (2.35), от $p$ и от норм $\|\Lambda\|_{L_\infty}$, $\| \widetilde{\Lambda} \|_{L_p(\Omega)}$, а $c'_r$ и $c''_r$ зависят от тех же параметров и от $r$.
Из теоремы 10.11 легко вывести аппроксимации решения и потока в строго внутренней подобласти.
Теорема 10.20. Пусть выполнены условия теоремы 10.17. Пусть $\mathcal{O}'$ – строго внутренняя подобласть области $\mathcal{O}$. Пусть $\delta = \operatorname{dist} \{\mathcal{O}'; \partial \mathcal{O} \}$. Тогда при $0< \varepsilon \leqslant \varepsilon_1$ и $\varepsilon^2 \leqslant t < T$ справедливы оценки
Постоянные $\mathrm{C}_5$, $\mathrm{C}_6$, $\widetilde{\mathrm{C}}_5$, $\widetilde{\mathrm{C}}_6$ зависят только от исходных данных (2.35), а $\mathfrak{c}_r$ и $\widetilde{\mathfrak{c}}_r$ зависят от тех же параметров и от $r$.
Наконец, при дополнительных предположениях о решениях вспомогательных задач из теоремы 10.12 нетрудно извлечь следующий результат.
Теорема 10.21. Пусть выполнены условия теоремы 10.19. Кроме того, пусть выполнены условия 8.1 и 8.2. Пусть функции $\check{\mathbf{v}}_{\varepsilon}$ и $\check{\mathbf{q}}_{\varepsilon}$ определены в (10.30) и (10.31) соответственно. Тогда при $0< \varepsilon \leqslant \varepsilon_1$ и $\varepsilon^2 \leqslant t < T$ справедливы оценки
Постоянные $\mathrm{C}_5$, $\mathrm{C}_7$, $\widetilde{\mathrm{C}}_5$ и $\widetilde{\mathrm{C}}_7$ – те же, что в теореме 10.12. Постоянные $\mathfrak{c}_{r}'$ и $\mathfrak{c}_{r}'' $ зависят от исходных данных (2.35), от $p$ и от норм $\| \Lambda\|_{L_\infty}$, $\| \widetilde{\Lambda}\|_{L_p(\Omega)}$.
§ 11. Пример применения общих результатов
Для эллиптических систем во всем пространстве $\mathbb{R}^d$ рассматриваемый пример изучался в [10] и [37]. Для эллиптических и параболических систем в ограниченной области при условии Дирихле пример рассматривался в [40] и [42] соответственно.
11.1. Скалярный эллиптический оператор с сингулярным потенциалом порядка $\varepsilon^{-1}$
Пусть $n=1$, $m=d$, $b(\mathbf{D})=\mathbf{D}$, а $g(\mathbf{x})$ – $\Gamma$-периодическая симметричная $(d\times d)$-матрица-функция с вещественными элементами, причем $g(\mathbf{x})>0$ и $g,g^{-1}\in L_\infty$. Условие (1.5) справедливо при $\alpha_0=\alpha_1 =1$. Очевидно, условие 2.1 выполнено, причем $k_1=1$, $k_2=0$. Имеем $b(\mathbf{D})^*g^\varepsilon (\mathbf{x}) b(\mathbf{D})=-\operatorname{div}g^\varepsilon (\mathbf{x})\nabla$.
Далее, пусть $\mathbf{A}(\mathbf{x})=\operatorname{col}\{ A_1(\mathbf{x}),\dots,A_d(\mathbf{x})\}$, где $A_j(\mathbf{x})$, $j=1,\dots,d$, – $\Gamma$-периодические вещественные функции, причем
при естественном условии (условии Неймана) на $\partial\mathcal{O}$. Этот оператор можно интерпретировать как оператор Шрёдингера с быстро осциллирующими метрикой $g^\varepsilon$, магнитным потенциалом $\mathbf{A}^\varepsilon$ и электрическим потенциалом $\varepsilon^{-1}v^\varepsilon +\mathcal{V}^\varepsilon$, содержащим сингулярное слагаемое $\varepsilon^{-1}v^\varepsilon$. Точное определение оператора $\mathfrak{B}_{N,\varepsilon}$ дается через квадратичную форму
Комплексные функции $a_j(\mathbf{x})$ заданы выражениями
$$
\begin{equation}
a_j(\mathbf{x})=-\eta_j(\mathbf{x})+i\xi_j(\mathbf{x}),\qquad j=1,\dots, d.
\end{equation}
\tag{11.3}
$$
Здесь $\eta_j(\mathbf{x})$ – компоненты вектор-функции $\boldsymbol{\eta}(\mathbf{x})=g(\mathbf{x})\mathbf{A}(\mathbf{x})$, а функции $\xi_j(\mathbf{x})$ имеют вид $\xi_j (\mathbf{x})=-\partial_j \Phi (\mathbf{x})$, где $\Phi (\mathbf{x})$ – $\Gamma$-периодическое решение задачи $\Delta \Phi(\mathbf{x})=v(\mathbf{x})$, $\int_\Omega \Phi(\mathbf{x})\,d\mathbf{x}=0$. При этом выполнено
Легко убедиться, что функции (11.3) удовлетворяют условию (1.9) с подходящим показателем $\rho '$, зависящим от $\rho$ и $s$, причем нормы $\|a_j\|_{L_{\rho '}(\Omega)}$ контролируются через $\| g\|_{L_\infty}$, $\|\mathbf{A}\|_{L_\rho (\Omega)}$, $\| v\|_{L_s(\Omega)}$ и параметры решетки $\Gamma$. Функция (11.2) удовлетворяет условию (1.10) с подходящим показателем $s'=\min \{ s;\rho/2\}$.
Пусть $Q_0(\mathbf{x})$ – положительно определенная и ограниченная $\Gamma$-периодическая функция. Рассмотрим положительно определенный оператор $\mathcal{B}_{N,\varepsilon}:=\mathfrak{B}_{N,\varepsilon}+\lambda Q_0^\varepsilon$. Здесь постоянная $\lambda$ выбрана из условия (2.34) для оператора, коэффициенты $g$, $a_j$, $j=1,\dots,d$, $Q$ и $Q_0$ которого определены выше. Оператор $\mathcal{B}_{N,\varepsilon}$ можно записать в виде
Опишем эффективный оператор. В рассматриваемом случае $\Gamma$-периодическое решение задачи (1.20) является матрицей-строкой: $\Lambda (\mathbf{x})=i\Psi (\mathbf{x})$, $\Psi (\mathbf{x})=(\psi_1(\mathbf{x}),\dots,\psi_d(\mathbf{x}))$, где $\psi_j \in\widetilde{H}^1(\Omega)$ – решение задачи
Здесь $\mathbf{e}_j$, $j=1,\dots,d$, – стандартные орты в $\mathbb{R}^d$. Ясно, что функции $\psi_j(\mathbf{x})$ вещественнозначные, а элементы матрицы $\Lambda(\mathbf{x})$ чисто мнимые. Согласно (1.22) столбцы $(d\times d)$-матрицы-функции $\widetilde{g}(\mathbf{x})$ – это $g(\mathbf{x})(\nabla \psi_j (\mathbf{x})+\mathbf{e}_j)$, $j=1,\dots,d$. Эффективная матрица определена в соответствии с (1.21): $g^0=|\Omega|^{-1}\int_\Omega\widetilde{g}(\mathbf{x})\,d\mathbf{x}$. Ясно, что $\widetilde{g}(\mathbf{x})$ и $g^0$ имеют вещественные элементы.
Согласно (11.3) и (11.4) периодическое решение задачи (1.30) представляется в виде $\widetilde{\Lambda}(\mathbf{x}) =\widetilde{\Lambda}_1(\mathbf{x})+i\widetilde{\Lambda}_2(\mathbf{x})$, где вещественные $\Gamma$-периодические функции $\widetilde{\Lambda}_1(\mathbf{x})$ и $\widetilde{\Lambda}_2(\mathbf{x})$ являются решениями задач
где $\mathbf{A}^0=(g^0)^{-1}(V_1+\overline{g\mathbf{A}})$ и $\mathcal{V}^0=\overline{\mathcal{V}}+\overline{\langle g\mathbf{A},\mathbf{A}\rangle}-\langle g^0\mathbf{A}^0,\mathbf{A}^0\rangle -W$.
11.2. Эллиптические результаты
Согласно замечанию 8.5 в рассматриваемом случае справедливы условия 8.1 и 8.2, причем нормы $\| \Lambda\|_{L_\infty}$ и $\|\widetilde{\Lambda}\|_{L_\infty}$ оцениваются в терминах данных задачи (11.5). Поэтому можно использовать корректор (8.1), не содержащий сглаживающего оператора:
В соответствии с теоремами 3.2 и 8.6 справедлив следующий результат.
Предложение 11.1. Пусть выполнены предположения п. 11.1. Пусть $\zeta\in\mathbb{C}\setminus \mathbb{R}_+$, $\zeta =| \zeta | e^{i\phi}$, $0<\phi<2\pi$, $|\zeta| \geqslant 1$. Пусть число $\varepsilon_1$ подчинено условию 3.1. Тогда при $0<\varepsilon\leqslant\varepsilon_1$ справедливы оценки
Здесь $c(\phi)$ – величина (1.40). Постоянные $\mathcal{C}_3$, $\mathcal{C}_4$, $\mathcal{C}_{12}$, $\widetilde{\mathcal{C}}_4$ и $\widetilde{\mathcal{C}}_{12}$ зависят только от исходных данных (11.5).
“Другая” аппроксимация оператора $(\mathcal{B}_{N,\varepsilon}-\zeta Q_0^\varepsilon )^{-1}$ получается на основании теорем 9.1 и 9.7.
Предложение 11.2. Пусть выполнены предположения п. 11.1. Пусть $\zeta \in \mathbb{C}\setminus [c_\flat,\infty)$, где $c_\flat = \frac{1}{2}\|g^{-1}\|^{-1}_{L_\infty} \| Q_0\|^{-1}_{L_\infty}$. Пусть $\varrho_\flat (\zeta )$ – величина (9.1). Пусть число $\varepsilon_1$ подчинено условию 3.1. Тогда при $0<\varepsilon\leqslant \varepsilon_1$ имеют место оценки
Постоянные $\mathfrak{C}_1$, $\mathfrak{C}_2$, $\mathfrak{C}_{16}$, $\widetilde{\mathfrak{C}}_2$ и $\widetilde{\mathfrak{C}}_{16}$ зависят только от исходных данных (11.5).
Здесь применимы также теоремы 8.13 и 9.11 об аппроксимациях оператора $(\mathcal{B}_{N,\varepsilon }-\zeta Q_0^\varepsilon )^{-1}$ в строго внутренней подобласти.
11.3. Параболические результаты
Обсудим кратко параболическую задачу. Пусть $u_\varepsilon(\mathbf{x},t)$ – решение начально-краевой задачи
Применяя теоремы 10.2 и 10.5, получаем следующее утверждение.
Предложение 11.3. Пусть выполнены предположения п. 11.1. Пусть $u_\varepsilon(\mathbf{x},t)$ – решение задачи (11.6) и $u_0(\mathbf{x},t)$ – решение усредненной задачи (11.7). Пусть число $\varepsilon_1$ подчинено условию 3.1. Тогда при $0<\varepsilon\leqslant\varepsilon_1$ и $t>0$ справедливы оценки
Здесь $c_\flat = \frac{1}{2}\|g^{-1}\|^{-1}_{L_\infty} \| Q_0\|^{-1}_{L_\infty}$. Постоянные $\mathrm{C}_1$, $\mathrm{C}_3$, $\widetilde{\mathrm{C}}_3$ зависят только от исходных данных (11.5).
Применима также теорема 10.12, дающая аппроксимацию решения $u_\varepsilon(\,{\cdot}\,,t)$ в $H^1(\mathcal{O}')$. Можно рассмотреть и начально-краевую задачу для неоднородного уравнения (аналог (11.6) с дополнительным слагаемым $F(\mathbf{x},t)$ в правой части уравнения) и применить теоремы 10.14, 10.15, 10.19 и 10.21.
Список литературы
1.
A. Bensoussan, J.-L. Lions, G. Papanicolaou, Asymptotic analysis for periodic structures, Stud. Math. Appl., 5, North-Holland Publishing Co., Amsterdam–New York, 1978, xxiv+700 pp.
2.
Н. С. Бахвалов, Г. П. Панасенко, Осреднение процессов в периодических средах. Математические задачи механики композиционных материалов, Наука, М., 1984, 352 с. ; англ. пер.: N. Bakhvalov, G. Panasenko, Homogenisation: averaging processes in periodic media. Mathematical problems in the mechanics of composite materials, Math. Appl. (Soviet Ser.), 36, Kluwer Acad. Publ., Dordrecht, 1989, xxxvi+366 с.
3.
О. А. Олейник, Г. A. Иосифьян, A. С. Шaмaев, Мaтемaтические зaдaчи теоpии сильно неодноpодных упpугих сpед, Моск. гос. ун-т, М., 1990, 312 с. ; англ. пер.: O. A. Oleinik, A. S. Shamaev, G. A. Yosifian, Mathematical problems in elasticity and homogenization, Stud. Math. Appl., 26, North-Holland Publishing Co., Amsterdam, 1992, xiv+398 с.
4.
В. В. Жиков, С. М. Козлов, О. А. Олейник, Усреднение дифференциальных операторов, Физматлит, М., 1993, 464 с. ; англ. пер.: V. V. Jikov, S. M. Kozlov, O. A. Oleinik, Homogenization of differential operators and integral functionals, Springer-Verlag, Berlin, 1994, xii+570 с.
5.
T. A. Suslina, “Homogenization of the Neumann problem for elliptic systems with periodic coefficients”, SIAM J. Math. Anal., 45:6 (2013), 3453–3493
6.
Т. А. Суслина, “Усреднение эллиптических операторов с периодическими коэффициентами в зависимости от спектрального параметра”, Алгебра и анализ, 27:4 (2015), 87–166; англ. пер.: T. A. Suslina, “Homogenization of elliptic operators with periodic coefficients in dependence of the spectral parameter”, St. Petersburg Math. J., 27:4 (2016), 651–708
7.
М. Ш. Бирман, Т. А. Суслина, “Периодические дифференциальные операторы второго порядка. Пороговые свойства и усреднения”, Алгебра и анализ, 15:5 (2003), 1–108; англ. пер.: M. Sh. Birman, T. A. Suslina, “Second order periodic differential operators. Threshold properties and homogenization”, St. Petersburg Math. J., 15:5 (2004), 639–714
8.
М. Ш. Бирман, Т. А. Суслина, “Усреднение периодических эллиптических дифференциальных операторов с учетом корректора”, Алгебра и анализ, 17:6 (2005), 1–104; англ. пер.: M. Sh. Birman, T. A. Suslina, “Homogenization with corrector term for periodic elliptic differential operators”, St. Petersburg Math. J., 17:6 (2006), 897–973
9.
M. Ш. Бирман, Т. А. Суслина, “Усреднение периодических дифференциальных операторов с учетом корректора. Приближение решений в классе Соболева $H^1(\mathbb{R}^d)$”, Алгебра и анализ, 18:6 (2006), 1–130; англ. пер.: M. Sh. Birman, T. A. Suslina, “Homogenization with corrector for periodic differential operators. Approximation of solutions in the Sobolev class $H^1(\mathbb R^d)$”, St. Petersburg Math. J., 18:6 (2007), 857–955
10.
Т. А. Суслина, “Усреднение в классе Соболева $H^1(\mathbb{R}^d)$ для периодических эллиптических дифференциальных операторов второго порядка при включении членов первого порядка”, Алгебра и анализ, 22:1 (2010), 108–222; англ. пер.: T. A. Suslina, “Homogenization in the Sobolev class $H^1(\mathbb{R}^d)$ for second order periodic elliptic operators with the inclusion of first order terms”, St. Petersburg Math. J., 22:1 (2011), 81–162
11.
Т. А. Суслина, “Усреднение эллиптических систем с периодическими коэффициентами: операторные оценки погрешности в $L_2(\mathbb{R}^d)$ с учетом корректора”, Алгебра и анализ, 26:4 (2014), 195–263; англ. пер.: T. A. Suslina, “Homogenization of elliptic systems with periodic coefficients: operator error estimates in $L_2(\mathbb{R}^d)$ with corrector taken into account”, St. Petersburg Math. J., 26:4 (2015), 643–693
12.
Т. А. Суслина, “Об усреднении периодических параболических систем”, Функц. анализ и его прил., 38:4 (2004), 86–90; англ. пер.: T. A. Suslina, “On homogenization of periodic parabolic systems”, Funct. Anal. Appl., 38:4 (2004), 309–312
13.
T. A. Suslina, “Homogenization of a periodic parabolic Cauchy problem”, Nonlinear equations and spectral theory, Amer. Math. Soc. Transl. Ser. 2, 220, Adv. Math. Sci., 59, Amer. Math. Soc., Providence, RI, 2007, 201–233
14.
T. Suslina, “Homogenization of a periodic parabolic Cauchy problem in the Sobolev space $H^1(\mathbb{R}^d)$”, Math. Model. Nat. Phenom., 5:4 (2010), 390–447
15.
Ю. М. Мешкова, “Усреднение задачи Коши для параболических систем с периодическими коэффициентами”, Алгебра и анализ, 25:6 (2013), 125–177; англ. пер.: Yu. M. Meshkova, “Homogenization of the Cauchy problem for parabolic systems with periodic coefficients”, St. Petersburg Math. J., 25:6 (2014), 981–1019
16.
В. В. Жиков, “Об операторных оценках в теории усреднения”, Докл. РАН, 403:3 (2005), 305–308; англ. пер.: V. V. Zhikov, “On operator estimates in homogenization theory”, Dokl. Math., 72:1 (2005), 534–538
17.
В. В. Жиков, “О некоторых оценках из теории усреднения”, Докл. РАН, 406:5 (2006), 597–601; англ. пер.: V. V. Zhikov, “Some estimates from homogenization theory”, Dokl. Math., 73:1 (2006), 96–99
18.
V. V. Zhikov, S. E. Pastukhova, “On operator estimates for some problems in homogenization theory”, Russ. J. Math. Phys., 12:4 (2005), 515–524
19.
V. V. Zhikov, S. E. Pastukhova, “Estimates of homogenization for a parabolic equation with periodic coefficients”, Russ. J. Math. Phys., 13:2 (2006), 224–237
20.
В. В. Жиков, С. Е. Пастухова, “Об операторных оценках в теории усреднения”, УМН, 71:3(429) (2016), 27–122; англ. пер.: V. V. Zhikov, S. E. Pastukhova, “Operator estimates in homogenization theory”, Russian Math. Surveys, 71:3 (2016), 417–511
21.
Д. И. Борисов, “Асимптотики решений эллиптических систем с быстро осциллирующими коэффициентами”, Алгебра и анализ, 20:2 (2008), 19–42; англ. пер.: D. Borisov, “Asymptotics for the solutions of elliptic systems with rapidly oscillating coefficients”, St. Petersburg Math. J., 20:2 (2009), 175–191
22.
N. N. Senik, “Homogenization for non-self-adjoint periodic elliptic operators on an infinite cylinder”, SIAM J. Math. Anal., 49:2 (2017), 874–898
23.
Sh. Moskow, M. Vogelius, “First-order corrections to the homogenised eigenvalues of a periodic composite medium. A convergence proof”, Proc. Roy. Soc. Edinburgh Sect. A, 127:6 (1997), 1263–1299
24.
S. Moskow, M. Vogelius, First order corrections to the homogenized eigenvalues of a periodic composite medium. The case of Neumann boundary conditions, preprint, Rutgers Univ., 1997
25.
G. Griso, “Error estimate and unfolding for periodic homogenization”, Asymptot. Anal., 40:3-4 (2004), 269–286
26.
G. Griso, “Interior error estimate for periodic homogenization”, Anal. Appl. (Singap.), 4:1 (2006), 61–79
27.
C. E. Kenig, Fanghua Lin, Zhongwei Shen, “Convergence rates in $L^2$ for elliptic homogenization problems”, Arch. Ration. Mech. Anal., 203:3 (2012), 1009–1036
28.
М. А. Пахнин, Т. А. Суслина, “Операторные оценки погрешности при усреднении эллиптической задачи Дирихле в ограниченной области”, Алгебра и анализ, 24:6 (2012), 139–177; англ. пер.: M. A. Pakhnin, T. A. Suslina, “Operator error estimates for homogenization of the elliptic Dirichlet problem in a bounded domain”, St. Petersburg Math. J., 24:6 (2013), 949–976
29.
T. A. Suslina, “Homogenization of the Dirichlet problem for elliptic systems: $L_2$-operator error estimates”, Mathematika, 59:2 (2013), 463–476
30.
Qiang Xu, “Uniform regularity estimates in homogenization theory of elliptic system with lower order terms”, J. Math. Anal. Appl., 438:2 (2016), 1066–1107
31.
Qiang Xu, “Uniform regularity estimates in homogenization theory of elliptic systems with lower order terms on the Neumann boundary problem”, J. Differential Equations, 261:8 (2016), 4368–4423
32.
Qiang Xu, “Convergence rates for general elliptic homogenization problems in Lipschitz domains”, SIAM J. Math. Anal., 48:6 (2016), 3742–3788
33.
Zhongwei Shen, Periodic homogenization of elliptic systems, Oper. Theory Adv. Appl., 269, Adv. Partial Differ. Equ. (Basel), Birkhäuser/Springer, Cham, 2018, ix+291 pp.
34.
Zhongwei Shen, Jinping Zhuge, “Convergence rates in periodic homogenization of systems of elasticity”, Proc. Amer. Math. Soc., 145:3 (2017), 1187–1202
35.
Yu. M. Meshkova, T. A. Suslina, “Homogenization of initial boundary value problems for parabolic systems with periodic coefficients”, Appl. Anal., 95:8 (2016), 1736–1775
36.
Jun Geng, Zhongwei Shen, “Convergence rates in parabolic homogenization with time-dependent periodic coefficients”, J. Funct. Anal., 272:5 (2017), 2092–2113
37.
Yu. M. Meshkova, T. A. Suslina, “Two-parametric error estimates in homogenization of second order elliptic systems in $\mathbb{R}^d$”, Appl. Anal., 95:7 (2016), 1413–1448
38.
Т. А. Суслина, “Теоретико-операторный подход к усреднению уравнений типа Шрёдингера с периодическими коэффициентами”, УМН, 78:6(474) (2023), 47–178; англ. пер.: T. A. Suslina, “Operator-theoretic approach to the homogenization of Schrödinger-type equations with periodic coefficients”, Russian Math. Surveys, 78:6 (2023), 1023–1154
39.
О. А. Ладыженская, Краевые задачи математической физики, Наука, М., 1973, 407 с. ; англ. пер.: O. A. Ladyzhenskaya, The boundary value problems of mathematical physics, Appl. Math. Sci., 49, Springer-Verlag, New York, 1985, xxx+322 с.
40.
Yu. M. Meshkova, T. A. Suslina, Homogenization of the Dirichlet problem for elliptic systems: two-parametric error estimates, 2017, arXiv: 1702.00550v4
41.
Ю. М. Мешкова, Т. А. Суслина, “Усреднение задачи Дирихле для эллиптических и параболических систем с периодическими коэффициентами”, Функц. анализ и его прил., 51:3 (2017), 87–93; англ. пер.: Yu. M. Meshkova, T. A. Suslina, “Homogenization of the Dirichlet problem for elliptic and parabolic systems with periodic coefficients”, Funct. Anal. Appl., 51:3 (2017), 230–235
42.
Ю. М. Мешкова, Т. А. Суслина, “Усреднение первой начально-краевой задачи для параболических систем: операторные оценки погрешности”, Алгебра и анализ, 29:6 (2017), 99–158; англ. пер.: Yu. M. Meshkova, T. A. Suslina, “Homogenization of the first initial boundary-value problem for parabolic systems: operator error estimates”, St. Petersburg Math. J., 29:6 (2018), 935–978
43.
Yu. M. Meshkova, “On homogenization of the first initial-boundary value problem for periodic hyperbolic systems”, Appl. Anal., 99:9 (2020), 1528–1563
44.
J. Nečas, Direct methods in the theory of elliptic equations, Transl. from the French, Springer Monogr. Math., Springer, Heidelberg, 2012, xvi+372 pp.
45.
W. McLean, Strongly elliptic systems and boundary integral equations, Cambridge Univ. Press, Cambridge, 2000, xiv+357 pp.
46.
В. А. Кондратьев, С. Д. Эйдельман, “Об условиях на граничную поверхность в теории эллиптических граничных задач”, Докл. АН СССР, 246:4 (1979), 812–815; англ. пер.: V. A. Kondrat'ev, S. D. Èĭdel'man, “On conditions on the boundary surface in the theory of elliptic boundary value problems”, Soviet Math. Dokl., 20 (1979), 561–563
47.
В. Г. Мазья, Т. О. Шапошникова, Мультипликаторы в пространствах дифференцируемых функций, Изд-во Ленингр. ун-та, Л., 1986, 404 с. ; англ. пер.: V. G. Maz'ya, T. O. Shaposhnikova, Theory of multipliers in spaces of differentiable functions, Monogr. Stud. Math., 23, Pitman, Boston, MA, 1985, xiii+344 с.
48.
И. Стейн, Сингулярные интегралы и дифференциальные свойства функций, Мир, М., 1973, 342 с. ; пер. с англ.: E. M. Stein, Singular integrals and differentiability properties of functions, Princeton Math. Ser., 30, Princeton Univ. Press, Princeton, NJ, 1970, xiv+290 с.
49.
V. S. Rychkov, “On restrictions and extensions of the Besov and Triebel–Lizorkin spaces with respect to Lipschitz domains”, J. London Math. Soc. (2), 60:1 (1999), 237–257
50.
О. А. Ладыженская, Н. Н. Уральцева, Линейные и квазилинейные уравнения эллиптического типа, Наука, М., 1964, 538 с. ; англ. пер.: O. A. Ladyzhenskaya, N. N. Ural'tseva, Linear and quasilinear elliptic equations, Academic Press, New York–London, 1968, xviii+495 с.
51.
Т. Като, Теория возмущений линейных операторов, Мир, М., 1972, 740 с. ; пер. с англ.: T. Kato, Perturbation theory for linear operators, Grundlehren Math. Wiss., 132, Springer-Verlag New York, Inc., New York, 1966, xix+592 с.
Образец цитирования:
Т. А. Суслина, “Усреднение эллиптических и параболических уравнений с периодическими коэффициентами в ограниченной области при условии Неймана”, Изв. РАН. Сер. матем., 88:4 (2024), 84–167; Izv. Math., 88:4 (2024), 678–759