Доклады Российской академии наук. Математика, информатика, процессы управления, 2021, T. 498, № 1, стр. 10-15

ОБ АТТРАКТОРАХ УРАВНЕНИЙ РЕАКЦИИ–ДИФФУЗИИ В ПОРИСТОЙ ОРТОТРОПНОЙ СРЕДЕ

К. А. Бекмаганбетов 12*, В. В. Чепыжов 34**, Г. А. Чечкин 256***

1 Московский государственный университет имени М.В. Ломоносова, Казахстанский филиал
Нур-Султан, Казахстан

2 Институт математики и математического моделирования
Алматы, Казахстан

3 Институт проблем передачи информации им. А.А. Харкевича Российской академии наук
Москва, Россия

4 Национальный исследовательский университет “Высшая школа экономики”
Москва, Россия

5 Московский государственный университет имени М.В. Ломоносова
Москва, Россия

6 Институт математики с компьютерным центром – подразделение Уфимского федерального исследовательского центра Российской академии наук
Уфа, Россия

* E-mail: bekmaganbetov-ka@yandex.kz
** E-mail: chep@iitp.ru
*** E-mail: chechkin@mech.math.msu.su

Поступила в редакцию 18.02.2021
После доработки 18.02.2021
Принята к публикации 09.03.2021

Полный текст (PDF)

Аннотация

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

Ключевые слова: аттракторы, усреднение, уравнение реакции–диффузии, нелинейные уравнения, слабая сходимость, перфорированная область, быстро осциллирующие члены, странный член

ВВЕДЕНИЕ

Интерес к задачам в перфорированных областях возник в связи прикладными задачами биологии, механики и инженерии. Работы, посвященные асимптотическому анализу таких задач, см., например, [15] и библиографию в этих работах. Особенно интересно изучать задачи, в которых может не выполняться теорема единственности. В таких задачах проводится усреднение соответствующих аттракторов (см. рис. 111, например). Аттракторы характеризуют всю динамику рассматриваемой модели (см., например, монографии [68] и ссылки в них). В работе [9] изучалось усреднение аттракторов скалярных эволюционных уравнений с диссипацией в периодически перфорированной области. В настоящей работе мы рассматриваем начально-краевую задачу для системы нелинейных дифференциальных уравнений с быстро осциллирующими членами в перфорированной области с третьим краевым условием на границе полостей.

Рис. 1.

Аттрактор.

1. ПОСТАНОВКА ЗАДАЧИ И ПРЕДВАРИТЕЛЬНЫЕ ЗАМЕЧАНИЯ

Пусть $\Omega $ – ограниченная область в ${{\mathbb{R}}^{n}}$, $n \geqslant 3$, с кусочно-гладкой границей $\partial \Omega $. Пусть ${{G}_{0}}$ – область, принадлежащая $Y = {{\left( { - \tfrac{1}{2},\tfrac{1}{2}} \right)}^{n}}$, такая, что ${{\bar {G}}_{0}}$ является компактом, диффеоморфным шару.

Пусть $\delta > 0$ и M – некоторое множество, введем следующее обозначение: δM = $\{ x{\kern 1pt} :\;{{\delta }^{{ - 1}}}x \in M\} $. Предположим, что $\varepsilon > 0$ достаточно мало, чтобы

${{\varepsilon }^{{n/(n - 2)}}}{{G}_{0}} \subset \varepsilon Y.$
Для $j \in {{\mathbb{Z}}^{n}}$ определим

$P_{\varepsilon }^{j} = \varepsilon j,\quad Y_{\varepsilon }^{j} = P_{\varepsilon }^{j} + \varepsilon Y,\quad G_{\varepsilon }^{j} = P_{\varepsilon }^{j} + {{\varepsilon }^{{n/(n - 2)}}}{{G}_{0}}.$

Определим область ${{\tilde {\Omega }}_{\varepsilon }} = \{ x \in \Omega {\kern 1pt} :\;\rho (x,\partial \Omega ) > \sqrt n \varepsilon \} $ и множество допустимых индексов

${{\Upsilon }_{\varepsilon }} = \{ j \in {{\mathbb{Z}}^{n}}{\kern 1pt} :\;\,\,G_{\varepsilon }^{j} \cap {{\bar {\tilde {\Omega }}}_{\varepsilon }} \ne \emptyset \} .$

Заметим, что $\left| {{{\Upsilon }_{\varepsilon }}} \right| \cong d{{\varepsilon }^{{ - n}}}$, где $d > 0$ – некоторая постоянная. Рассмотрим область

${{\Omega }_{\varepsilon }} = \Omega {{\backslash }}{{\bar {G}}_{\varepsilon }},\quad {\text{где}}\quad {{G}_{\varepsilon }} = \bigcup\limits_{j \in {{\Upsilon }_{\varepsilon }}} \,G_{\varepsilon }^{j}.$

Введем следующие обозначения:

${{Q}_{\varepsilon }} = {{\Omega }_{\varepsilon }} \times (0, + \infty ),\quad Q = \Omega \times (0, + \infty ).$

Мы изучаем асимптотическое поведение траекторных аттракторов начально-краевой задачи

$\frac{{\partial {{u}_{\varepsilon }}}}{{\partial t}} = H\Delta {{u}_{\varepsilon }} - a\left( {x,\frac{x}{\varepsilon }} \right)f({{u}_{\varepsilon }}) + g\left( {x,\frac{x}{\varepsilon }} \right),\quad x \in {{\Omega }_{\varepsilon }},$
(1)
$\begin{gathered} \frac{{\partial {{u}_{\varepsilon }}}}{{\partial \nu }} + {{\varepsilon }^{{n/(2 - n)}}}B_{\varepsilon }^{j}(x){{u}_{\varepsilon }} = 0,\quad x \in \partial G_{\varepsilon }^{j}, \\ j \in {{\Upsilon }_{\varepsilon }},\quad t \in (0, + \infty ), \\ \end{gathered} $
$\begin{gathered} {{u}_{\varepsilon }} = 0,\quad x \in \partial \Omega , \\ {{u}_{\varepsilon }} = U(x),\quad x \in {{\Omega }_{\varepsilon }},\quad t = 0, \\ \end{gathered} $
где $u\, = \,{{({{u}^{1}}, \ldots ,{{u}^{N}})}^{ \top }}$, $f\, = \,{{({{f}^{1}}, \ldots ,{{f}^{N}})}^{ \top }}$ и g = ${{({{g}^{1}}, \ldots ,{{g}^{N}})}^{ \top }}$, ν – вектор единичной внешней нормали к границе. Функция $a(x,y) \in C(\bar {\Omega } \times {{\mathbb{R}}^{n}})$ такая, что 0 < a0 ≤ ≤ $a(x,y) \leqslant {{A}_{0}}$ с некоторыми постоянными ${{a}_{0}}$, ${{A}_{0}}$, а функция ${{a}_{\varepsilon }}(x) = a\left( {x,\tfrac{x}{\varepsilon }} \right)$ имеет среднее $\bar {a}(x)$ при $\varepsilon \to 0 + $ в пространстве ${{L}_{{\infty ,*w}}}(\Omega )$, т.е.
(2)
$\int\limits_\Omega \,a\left( {x,\frac{x}{\varepsilon }} \right)\varphi (x)dx \to \int\limits_\Omega \,\bar {a}(x)\varphi (x)dx\quad (\varepsilon \to 0 + )$
для любой функции $\varphi \in {{L}_{1}}(\Omega )$. Матрица кросс–диффузии H – это квадратная матрица размера N × N с постоянными коэффициентами, имеющая положительную симметричную часть $\frac{1}{2}(H + {{H}^{ \top }})$ ≥ ≥ αI, где $\alpha > 0$, а I – единичная матрица порядка N. Отметим, что матрица H не обязательно симметрична, также отметим, что принцип максимума может не выполнятся для нашей задачи.

Для вектор-функции $g\left( {x,y} \right)$ будем считать, что функции $g_{\varepsilon }^{i}(x)\, = \,{{g}^{i}}\left( {x,\tfrac{x}{\varepsilon }} \right)\, \in \,{{L}_{2}}(\Omega \, \times \,{{\mathbb{R}}^{n}})$ и имеют средние ${{\bar {g}}^{i}}(x)$ в пространстве $V{\kern 1pt} ' = {{H}^{{ - 1}}}(\Omega )$ при $\varepsilon \to 0 + $, т.е.

(3)
$\int\limits_\Omega \,{{g}^{i}}\left( {x,\frac{x}{\varepsilon }} \right)\varphi (x)dx \to \int\limits_\Omega \,{{\bar {g}}^{i}}(x)\varphi (x)dx\quad (\varepsilon \to 0 + )$
для любой функции $\varphi \in V = H_{0}^{1}(\Omega )$ и для всех $i = 1, \ldots ,N$.

В работе [10] приведены примеры функций вида $a\left( {x,\frac{x}{\varepsilon }} \right)$ и $g\left( {x,\frac{x}{\varepsilon }} \right)$, которые удовлетворяют условиям усреднения (2) и (3).

Здесь $B_{\varepsilon }^{j}(x)$ – диагональная матрица с ограниченными элементами вида

${{b}^{{11}}}\left( {x,\frac{{x - P_{\varepsilon }^{j}}}{{{{\varepsilon }^{{n/(n - 2)}}}}}} \right), \ldots ,{{b}^{{NN}}}\left( {x,\frac{{x - P_{\varepsilon }^{j}}}{{{{\varepsilon }^{{n/(n - 2)}}}}}} \right),\quad j \in {{\Upsilon }_{\varepsilon }},$
где ${{b}^{{kk}}}(x,y) \in C(\Omega \times {{\mathbb{R}}^{n}})$ – 1-периодические по $y$ функции такие, что
(4)
$0 < {{b}_{0}} \leqslant {{b}^{{kk}}}(x,y) \leqslant {{B}_{0}}$
с некоторыми постоянными b0, ${{B}_{0}}$ для всех $k = 1$, 2, ..., N.

Обозначим также вектор

$\bar {B}(x,y): = {{({{b}^{{11}}}(x,y), \ldots ,{{b}^{{NN}}}(x,y))}^{ \top }},$
а диагональную матрицу с элементами ${{b}^{{11}}}(x,y)$, ... ..., bNN(x, y) – через $B(x,y)$.

Предположим, что вектор-функция $f(v)\, \in \,C({{\mathbb{R}}^{N}}$; ${{\mathbb{R}}^{N}})$ удовлетворяет следующим неравенствам:

(5)
$\sum\limits_{i = 1}^N \,{\text{|}}{{f}^{i}}({v}){{{\text{|}}}^{{{{p}_{i}}/({{p}_{i}} - 1)}}} \leqslant {{C}_{0}}\left( {\sum\limits_{i = 1}^N \,{\text{|}}{{v}^{i}}{{{\text{|}}}^{{{{p}_{i}}}}} + 1} \right),$
$\sum\limits_{i = 1}^N \,{{\gamma }_{i}}{\text{|}}{{{v}}^{i}}{{{\text{|}}}^{{{{p}_{i}}}}} - C \leqslant \sum\limits_{i = 1}^N \,{{f}^{i}}({v}){{{v}}^{i}},\quad \forall {v} \in {{\mathbb{R}}^{N}},$
где ${{\gamma }_{i}} > 0$ для всех $i = 1, \ldots ,N$. Для определенности будем считать, что ${{p}_{N}} \geqslant {{p}_{{N - 1}}}$ ≥ … ≥ ${{p}_{1}} \geqslant 2$. Заметим, что выполнение условия Липшица для функции $f(v)$ относительно ${v}$ не предполагается.

Введем следующие обозначения для пространств ${\mathbf{H}}: = {{[{{L}_{2}}(\Omega )]}^{N}}$, ${{{\mathbf{H}}}_{\varepsilon }}: = {{[{{L}_{2}}({{\Omega }_{\varepsilon }})]}^{N}}$, V := := ${{[H_{0}^{1}(\Omega )]}^{N}}$, ${{{\mathbf{V}}}_{\varepsilon }}: = {{[{{H}^{1}}({{\Omega }_{\varepsilon }};\partial \Omega )]}^{N}}$ – множество всех вектор-функций из ${{[{{H}^{1}}({{\Omega }_{\varepsilon }})]}^{N}}$ с нулевым следом на $\partial \Omega $. Нормы в этих пространствах определяются, соответственно, следующим образом:

$\begin{gathered} {\text{||}}{v}{\text{|}}{{{\text{|}}}^{2}}: = \int\limits_\Omega \,\sum\limits_{i = 1}^N \,{\text{|}}{{{v}}^{i}}(x){{{\text{|}}}^{2}}dx, \\ {\text{||}}{v}{\text{||}}_{\varepsilon }^{2}: = \int\limits_{{{\Omega }_{\varepsilon }}} \,\sum\limits_{i = 1}^N \,{\text{|}}{{{v}}^{i}}(x){{{\text{|}}}^{2}}dx, \\ \end{gathered} $
$\begin{gathered} {\text{||}}{v}{\text{||}}_{1}^{2}: = \int\limits_\Omega \,\sum\limits_{i = 1}^N \,{\text{|}}\nabla {{{v}}^{i}}(x){{{\text{|}}}^{2}}dx, \\ {\text{||}}{v}{\text{||}}_{{1\varepsilon }}^{2}: = \int\limits_{{{\Omega }_{\varepsilon }}} \,\sum\limits_{i = 1}^N \,{\text{|}}\nabla {{{v}}^{i}}(x){{{\text{|}}}^{2}}dx. \\ \end{gathered} $

Напомним, что ${\mathbf{V}}{\kern 1pt} ': = {{[{{H}^{{ - 1}}}(\Omega )]}^{N}}$ – двойственное пространство к пространству V, кроме того, ${\mathbf{V}}_{\varepsilon }^{'}$ – двойственное пространство для ${{{\mathbf{V}}}_{\varepsilon }}$.

Пусть ${{q}_{i}} = \frac{{{{p}_{i}}}}{{{{p}_{i}} - 1}}$ для всех $i = 1, \ldots ,N$. Будем использовать следующие векторные обозначения ${\mathbf{p}} = ({{p}_{1}}, \ldots ,{{p}_{N}})$ и ${\mathbf{q}} = ({{q}_{1}}, \ldots ,{{q}_{N}})$, а также определим пространства

$\begin{gathered} {{{\mathbf{L}}}_{{\mathbf{p}}}}: = {{L}_{{{{p}_{1}}}}}(\Omega ) \times \ldots \times {{L}_{{{{p}_{N}}}}}(\Omega ), \\ {{{\mathbf{L}}}_{{{\mathbf{p}},\varepsilon }}}: = {{L}_{{{{p}_{1}}}}}({{\Omega }_{\varepsilon }}) \times \ldots \times {{L}_{{{{p}_{N}}}}}({{\Omega }_{\varepsilon }}), \\ \end{gathered} $
${{{\mathbf{L}}}_{{\mathbf{p}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{\mathbf{p}}}}): = {{L}_{{{{p}_{1}}}}}({{\mathbb{R}}_{ + }};{{L}_{{{{p}_{1}}}}}(\Omega )) \times \ldots \times {{L}_{{{{p}_{N}}}}}({{\mathbb{R}}_{ + }};{{L}_{{{{p}_{N}}}}}(\Omega )),$
${{{\mathbf{L}}}_{{\mathbf{p}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{{\mathbf{p}},\varepsilon }}})\,: = \,{{L}_{{{{p}_{1}}}}}({{\mathbb{R}}_{ + }};{{L}_{{{{p}_{1}}}}}({{\Omega }_{\varepsilon }})) \times \ldots \times {{L}_{{{{p}_{N}}}}}({{\mathbb{R}}_{ + }};{{L}_{{{{p}_{N}}}}}({{\Omega }_{\varepsilon }})).$

Как и в [7], будем исследовать слабые решения начально-краевой задачи (1), т.е. функции

${{u}_{\varepsilon }}(x,s) \in {\mathbf{L}}_{\infty }^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{H}}}_{\varepsilon }}) \cap {\mathbf{L}}_{2}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{V}}}_{\varepsilon }}) \cap {\mathbf{L}}_{{\mathbf{p}}}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{{\mathbf{p}},\varepsilon }}}),$
которые удовлетворяют задаче (1) в смысле обобщенных функций, т.е.
$\int\limits_{{{Q}_{\varepsilon }}} \,\frac{{\partial {{u}_{\varepsilon }}}}{{\partial t}} \cdot \psi dxdt + \int\limits_{{{Q}_{\varepsilon }}} \,H\nabla {{u}_{\varepsilon }} \cdot \nabla \psi dxdt + $
(7)
$\begin{gathered} \, + \int\limits_{{{Q}_{\varepsilon }}} \,{{a}_{\varepsilon }}(x)f({{u}_{\varepsilon }}) \cdot \psi dxdt + \\ \, + {{\varepsilon }^{{n/(2 - n)}}}\sum\limits_{j \in {{\Upsilon }_{\varepsilon }}} \,\int\limits_0^{ + \infty } \,\int\limits_{\partial G_{\varepsilon }^{j}} \,B_{\varepsilon }^{j}(x){{u}_{\varepsilon }} \cdot \psi dxdt = \\ \end{gathered} $
$\, = \int\limits_{{{Q}_{\varepsilon }}} \,{{g}_{\varepsilon }}(x) \cdot \psi dxdt$
для любых функций $\psi \in {\mathbf{C}}_{0}^{\infty }({{\mathbb{R}}_{ + }};{{{\mathbf{H}}}_{\varepsilon }})$. Здесь ${{y}_{1}} \cdot {{y}_{2}}$ означает скалярное произведение векторов y1, ${{y}_{2}} \in {{\mathbb{R}}^{N}}$.

Если ${{u}_{\varepsilon }}(x,t) \in {{{\mathbf{L}}}_{{\mathbf{p}}}}(0,M;{{{\mathbf{L}}}_{{{\mathbf{p}},\varepsilon }}})$, тогда из условия (5) следует, что $f(u(x,t)) \in {{{\mathbf{L}}}_{{\mathbf{q}}}}(0,M;{{{\mathbf{L}}}_{{{\mathbf{q}},\varepsilon }}})$. В то же время, если ${{u}_{\varepsilon }}(x,t) \in {{{\mathbf{L}}}_{2}}(0,M;{{{\mathbf{V}}}_{\varepsilon }})$, тогда

$\Delta {{u}_{\varepsilon }}(x,t) + {{g}_{\varepsilon }}\left( x \right) \in {{{\mathbf{L}}}_{2}}(0,M;{\mathbf{V}}_{\varepsilon }^{'})$.

Поэтому для произвольного слабого решения ${{u}_{\varepsilon }}(x,s)$ задачи (1) имеем 

$\frac{{\partial {{u}_{\varepsilon }}(x,t)}}{{\partial t}} \in {{{\mathbf{L}}}_{{\mathbf{q}}}}(0,M;{{{\mathbf{L}}}_{{{\mathbf{q}},\varepsilon }}}) + {{{\mathbf{L}}}_{2}}(0,M;{\mathbf{V}}_{\varepsilon }^{'}).$

Из теоремы вложения Соболева следует, что

${{{\mathbf{L}}}_{{\mathbf{q}}}}(0,M;{{{\mathbf{L}}}_{{{\mathbf{q}},\varepsilon }}}) + {{{\mathbf{L}}}_{2}}(0,M;{\mathbf{V}}_{\varepsilon }^{'}) \subset {{{\mathbf{L}}}_{{\mathbf{q}}}}(0,M;{\mathbf{H}}_{\varepsilon }^{{ - {\mathbf{r}}}}),$
где пространство ${\mathbf{H}}_{\varepsilon }^{{ - {\mathbf{r}}}}\,: = \,H_{\varepsilon }^{{ - {{r}_{1}}}}\, \times \, \ldots \, \times \,H_{\varepsilon }^{{ - {{r}_{N}}}}$, ${\mathbf{r}}\, = \,({{r}_{1}}, \ldots ,{{r}_{N}})$ и ${{r}_{i}} = max\left\{ {1,n\left( {\frac{1}{{{{q}_{i}}}} - \frac{1}{2}} \right)} \right\}$ для всех $i = 1, \ldots ,N$ (соболевское пространство с отрицательным показателем). Следовательно, для любого слабого решения ${{u}_{\varepsilon }}(x,t)$ задачи (1) имеем $\tfrac{{\partial {{u}_{\varepsilon }}(x,t)}}{{\partial t}} \in {{{\mathbf{L}}}_{{\mathbf{q}}}}(0,M;{\mathbf{H}}_{\varepsilon }^{{ - {\mathbf{r}}}})$.

Замечание 1. Существование слабого решения $u(x,s)$ задачи (1) для любой функции $U \in {{{\mathbf{H}}}_{\varepsilon }}$ и фиксированного ε, такого, что $u(x,0)$ = = U(x), может быть доказано стандартным способом (см., например, [6]). Это решение может быть не единственным, поскольку функция $f(v)$ удовлетворяет условиям (5), (6) и для нее не предполагается выполнение условия Липшица относительно $v$.

Для удобства будем опускать индекс ε в обозначениях пространств, там, где это не вызывает непонимания. Положим ${{E}_{1}} = {{{\mathbf{L}}}_{{\mathbf{p}}}} \cap {\mathbf{V}}$, ${{E}_{0}} = {{{\mathbf{H}}}^{{ - {\mathbf{r}}}}}$, $E = {\mathbf{H}}$ и $A(u) = H\Delta u - a( \cdot )f(u) + g( \cdot )$ и определим банаховы пространства для каждого отрезка $[{{t}_{1}},{{t}_{2}}] \in R$

(8)
$\begin{gathered} {{\mathcal{F}}_{{{{t}_{1}},{{t}_{2}}}}}: = {{{\mathbf{L}}}_{{\mathbf{p}}}}({{t}_{1}},{{t}_{2}};{{{\mathbf{L}}}_{p}}) \cap {{{\mathbf{L}}}_{2}}({{t}_{1}},{{t}_{2}};{\mathbf{V}}) \cap \\ \, \cap {{{\mathbf{L}}}_{\infty }}({{t}_{1}},{{t}_{2}};{\mathbf{H}}) \cap \left\{ {v\left| {\frac{{\partial v}}{{\partial t}}} \right. \in {{{\mathbf{L}}}_{{\mathbf{q}}}}({{t}_{1}},{{t}_{2}};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})} \right\} \\ \end{gathered} $
с нормой

(9)
$\begin{gathered} {\text{||}}v{\text{|}}{{{\text{|}}}_{{{{\mathcal{F}}_{{{{t}_{1}},{{t}_{2}}}}}}}}: = {\text{||}}v{\text{|}}{{{\text{|}}}_{{{{{\mathbf{L}}}_{{\mathbf{p}}}}({{t}_{1}},{{t}_{2}};{{{\mathbf{L}}}_{{\mathbf{p}}}})}}}\; + \;{\text{||}}v{\text{|}}{{{\text{|}}}_{{{{{\mathbf{L}}}_{2}}({{t}_{1}},{{t}_{2}};{\mathbf{V}})}}} + \\ \, + \;{\text{||}}v{\text{|}}{{{\text{|}}}_{{{{{\mathbf{L}}}_{\infty }}(0,M;{\mathbf{H}})}}} + \mathop {\left\| {\frac{{\partial v}}{{\partial t}}} \right\|}\nolimits_{{{{\mathbf{L}}}_{{\mathbf{q}}}}({{t}_{1}},{{t}_{2}};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})} . \\ \end{gathered} $

Положив ${{\mathcal{D}}_{{{{t}_{1}},{{t}_{2}}}}} = {{{\mathbf{L}}}_{{\mathbf{q}}}}({{t}_{1}},{{t}_{2}};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})$, получаем, что ${{\mathcal{F}}_{{{{t}_{1}},{{t}_{2}}}}} \subseteq {{\mathcal{D}}_{{{{t}_{1}},{{t}_{2}}}}}$, а если $u(s)\, \in \,{{\mathcal{F}}_{{{{t}_{1}},{{t}_{2}}}}}$, тогда $A(u(s)) \in {{\mathcal{D}}_{{{{t}_{1}},{{t}_{2}}}}}$. Далее обозначим

$\begin{gathered} \mathcal{F}_{ + }^{{{\text{loc}}}} = {\mathbf{L}}_{{\mathbf{p}}}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{\mathbf{p}}}}) \cap {\mathbf{L}}_{2}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{\mathbf{V}}) \cap \\ \, \cap {\mathbf{L}}_{\infty }^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{\mathbf{H}}) \cap \left\{ {{v}\left| {\frac{{\partial {v}}}{{\partial t}}} \right. \in {\mathbf{L}}_{{\mathbf{q}}}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})} \right\}, \\ \end{gathered} $
$\begin{gathered} \mathcal{F}_{{\varepsilon , + }}^{{{\text{loc}}}} = {\mathbf{L}}_{{\mathbf{p}}}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{{\mathbf{p}},\varepsilon }}}) \cap {\mathbf{L}}_{2}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{V}}}_{\varepsilon }}) \cap \\ \, \cap {\mathbf{L}}_{\infty }^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};{{{\mathbf{H}}}_{\varepsilon }}) \cap \left\{ {{v}\left| {\frac{{\partial {v}}}{{\partial t}}} \right. \in {\mathbf{L}}_{{\mathbf{q}}}^{{{\text{loc}}}}({{\mathbb{R}}_{ + }};\mathop {{{{\mathbf{H}}}_{\varepsilon }}}\nolimits^{ - {\mathbf{r}}} )} \right\}, \\ \end{gathered} $
а через $\mathcal{K}_{\varepsilon }^{ + }$ множество всех слабых решений задачи (1). Напомним, что для любой функции $U \in {\mathbf{H}}$ существует хотя бы одна траектория $u( \cdot ) \in \mathcal{K}_{\varepsilon }^{ + }$ такая, что $u(0) = U(x)$. Следовательно, пространство траекторий $\mathcal{K}_{\varepsilon }^{ + }$ задачи (1) не пусто и достаточно велико.

Ясно, что $\mathcal{K}_{\varepsilon }^{ + } \subset \mathcal{F}_{ + }^{{{\text{loc}}}}$ и пространство траекторий $\mathcal{K}_{\varepsilon }^{ + }$ является трансляционно-инвариантным, т.е. если $u(s) \in \mathcal{K}_{\varepsilon }^{ + }$, тогда и $u(h + s) \in \mathcal{K}_{\varepsilon }^{ + }$ для любых $h \geqslant 0$. Далее, используя норму пространства ${{{\mathbf{L}}}_{2}}({{t}_{1}},{{t}_{2}};{\mathbf{H}})$, определим метрики ${{\rho }_{{{{t}_{1}},{{t}_{2}}}}}( \cdot , \cdot )$ в пространствах ${{\mathcal{F}}_{{{{t}_{1}},{{t}_{2}}}}}$ следующим образом:

$\begin{gathered} {{\rho }_{{0,M}}}(u,v) = \mathop {\left( {\int\limits_0^M \,{\text{||}}u(s) - v(s){\text{|}}{{{\text{|}}}^{2}}ds} \right)}\nolimits^{1/2} , \\ \forall u( \cdot ),v( \cdot ) \in {{\mathcal{F}}_{{0,M}}}. \\ \end{gathered} $

Эти метрики порождают топологию $\Theta _{ + }^{{loc}}$ в пространстве $\mathcal{F}_{ + }^{{{\text{loc}}}}$ (соответственно $\Theta _{{\varepsilon , + }}^{{{\text{loc}}}}$ в $\mathcal{F}_{{\varepsilon , + }}^{{{\text{loc}}}}$). Напомним, что последовательность $\{ {{v}_{k}}\} \subset \mathcal{F}_{ + }^{{{\text{loc}}}}$ сходится к функции $v \in \mathcal{F}_{ + }^{{{\text{loc}}}}$ при $k \to \infty $ в $\Theta _{ + }^{{{\text{loc}}}}$, если ${\text{||}}{{v}_{k}}( \cdot ) - v( \cdot ){\text{|}}{{{\text{|}}}_{{{{{\mathbf{L}}}_{2}}(0,M;{\mathbf{H}})}}} \to 0$ $(k \to \infty )$ для любого $M > 0$. Топология $\Theta _{ + }^{{{\text{loc}}}}$ метризуема и соответствующее метрическое пространство является полным. Мы рассматриваем топологию в пространстве траекторий $\mathcal{K}_{\varepsilon }^{ + }$ задачи (1). Полугруппа сдвигов $\{ S(t)\} $, действующая на $\mathcal{K}_{\varepsilon }^{ + }$, непрерывна в рассматриваемой топологии $\Theta _{ + }^{{{\text{loc}}}}$.

Определим ограниченные множества в $\mathcal{K}_{\varepsilon }^{ + }$, используя банаховы пространства

(10)
$\begin{gathered} \mathcal{F}_{ + }^{b} = {\mathbf{L}}_{{\mathbf{p}}}^{b}({{\mathbb{R}}_{ + }};{{{\mathbf{L}}}_{{\mathbf{p}}}}) \cap {\mathbf{L}}_{2}^{b}({{\mathbb{R}}_{ + }};{\mathbf{V}}) \cap {{{\mathbf{L}}}_{\infty }}({{\mathbb{R}}_{ + }};{\mathbf{H}}) \cap \\ \, \cap \left\{ {{v}\left| {\frac{{\partial {v}}}{{\partial t}}} \right. \in {\mathbf{L}}_{{\mathbf{q}}}^{b}({{\mathbb{R}}_{ + }};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})} \right\}, \\ \end{gathered} $
$\mathcal{F}_{ + }^{b}$ – подпространство пространства $\mathcal{F}_{ + }^{{{\text{loc}}}}$.

Рассмотрим полугруппу сдвигов $\{ S(t)\} $ на $\mathcal{K}_{\varepsilon }^{ + },$ $S(t){\kern 1pt} :\;\mathcal{K}_{\varepsilon }^{ + } \to \mathcal{K}_{\varepsilon }^{ + }$, $t \geqslant 0$.

Пусть ${{\mathcal{K}}_{\varepsilon }}$ означает ядро задачи (1), которое состоит из всех слабых решений $u(s),s \in \mathbb{R}$, ограниченных в пространстве

$\begin{gathered} {{\mathcal{F}}^{b}} = {\mathbf{L}}_{{\mathbf{p}}}^{b}(\mathbb{R};{{{\mathbf{L}}}_{{\mathbf{p}}}}) \cap {\mathbf{L}}_{2}^{b}(\mathbb{R};{\mathbf{V}}) \cap {{{\mathbf{L}}}_{\infty }}(\mathbb{R};{\mathbf{H}}) \cap \\ \, \cap \left\{ {{v}\left| {\frac{{\partial {v}}}{{\partial t}}} \right. \in {\mathbf{L}}_{{\mathbf{q}}}^{b}(\mathbb{R};{{{\mathbf{H}}}^{{ - {\mathbf{r}}}}})} \right\}. \\ \end{gathered} $

Имеет место

Лемма 1. При выполнении условий (5), (6) задача (1) имеет траекторные аттракторы ${{\mathfrak{A}}_{\varepsilon }}$ в топологическом пространстве $\Theta _{ + }^{{{\text{loc}}}}$. Множество ${{\mathfrak{A}}_{\varepsilon }}$ равномерно (по $\varepsilon \in (0,1)$) ограничено в $\mathcal{F}_{ + }^{b}$ и компактно в $\Theta _{ + }^{{{\text{loc}}}}$. Более того,

${{\mathfrak{A}}_{\varepsilon }} = {{\Pi }_{ + }}{{K}_{\varepsilon }},$
ядро ${{\mathcal{K}}_{\varepsilon }}$непусто и равномерно (по $\varepsilon \in (0,1)$) ограничено в ${{\mathcal{F}}^{b}}$. Напомним, что пространства $\mathcal{F}_{ + }^{b}$ и $\Theta _{ + }^{{{\text{loc}}}}$ зависят от ε.

Доказательство этого предложения практически полностью совпадает с доказательством, приведенным в [8], для более частного случая. Существование поглощающего множества, ограниченного в $\mathcal{F}_{ + }^{b}$ и компактного в $\Theta _{ + }^{{{\text{loc}}}}$, доказывается так же, как и в [7].

2. УСРЕДНЕНИЕ АТТРАКТОРОВ

Чтобы определить “странный член” (потенциал в предельном уравнении), рассмотрим задачу

$\begin{gathered} - {{\Delta }_{y}}v = 0,\quad y \in {{\mathbb{R}}^{n}}{{\backslash }}{{G}_{0}}, \\ \frac{{\partial v}}{{\partial {{\nu }_{y}}}} + B(x,y)v = \bar {B}(x,y),\quad y \in \partial {{G}_{0}}, \\ v \to 0,\quad {\text{|}}y{\text{|}} \to \infty , \\ \end{gathered} $
где матрица $B(x,y)$ и вектор $\bar {B}(x,y)$ были определены выше. В этой задаче переменная $x$ играет роль медленного параметра. Определим предельный потенциал по следующей формуле:

(11)
${{V}^{{kk}}}(x) = \int\limits_{\partial {{G}_{0}}} \,\frac{\partial }{{\partial {{\nu }_{y}}}}{{v}^{k}}(x,y)d{{\sigma }_{y}},\quad k = 1, \ldots ,N.$

Усредненная (предельная) задача имеет следующий вид:

(12)
$\begin{gathered} \frac{{\partial u}}{{\partial t}} = H\Delta u - \bar {a}(x)f(u) - V(x)u + \bar {g}(x),\quad x \in \Omega , \\ u = 0,\quad x \in \partial \Omega , \\ u = U(x),\quad t = 0, \\ \end{gathered} $
где $V(x)$ – диагональная матрица с элементами ${{V}^{{kk}}}(x)$, $k = 1, \ldots ,N$. Очевидно, что задача (12) также имеет траекторный аттрактор $\bar {\mathfrak{A}}$ в пространстве траекторий ${{\bar {\mathcal{K}}}^{ + }}$, соответствующем задаче (12), и
$\bar {\mathfrak{A}} = {{\Pi }_{ + }}\bar {\mathcal{K}}{\text{,}}$
где $\bar {\mathcal{K}}$ – ядро задачи (12) в ${{\mathcal{F}}^{b}}$.

Имеет место следующее утверждение о сходимости.

Теорема 1. В топологическом пространстве $\Theta _{ + }^{{{\text{loc}}}}$ справедливо предельное соотношение

${{\mathfrak{A}}_{\varepsilon }} \to \bar {\mathfrak{A}}\quad при\quad \varepsilon \to 0 + .$

Более того

${{\mathcal{K}}_{\varepsilon }} \to \bar {\mathcal{K}}\quad при\quad \varepsilon \to 0 + \quad в\quad {{\Theta }^{{{\text{loc}}}}}.$

Замечание 2. Напомним, что пространства в теореме 1 зависят от ε. Все функции могут быть продолжены внутрь отверстий с сохранением соответствующих норм.

Если рассмотреть уравнения реакции–диффузии, для которых имеет место теорема единственности в задаче Коши. Для этого достаточно предположить, что нелинейный член $f(u)$ в системе уравнений (1) удовлетворяет условию

(13)
$\begin{gathered} (f({{v}_{1}}) - f({{v}_{2}}),{{v}_{1}} - {{v}_{2}}) \geqslant - C{\text{|}}{{v}_{1}} - {{v}_{2}}{{{\text{|}}}^{2}} \\ {\text{для}}\;{\text{всех}}\quad {{v}_{1}},{{v}_{2}} \in {{\mathbb{R}}^{N}}. \\ \end{gathered} $
(см. [7]). В [7] было доказано, что если выполнено (13), тогда уравнения (1) и (12) генерируют динамические полугруппы в H, имеющие глобальные аттракторы ${{\mathfrak{A}}_{\varepsilon }}$ и $\bar {\mathfrak{A}}$, ограниченные в пространстве ${\mathbf{V}} = {\mathbf{H}}_{0}^{1}(\Omega )$ (см. также [8]). Пусть

При этом выполняется следующее

Следствие 1. В условиях теоремы 1 имеет место предельное соотношение

${\text{dis}}{{{\text{t}}}_{{{{{\mathbf{H}}}^{{ - \delta }}}}}}({{\mathfrak{A}}_{\varepsilon }},\bar {\mathfrak{A}}) \to 0\quad (\varepsilon \to 0 + ).$

Список литературы

  1. Марченко В.А., Хруслов Е.Я. Краевые задачи в областях с мелкозернистой границей. Киев: Наукова думка, 1974.

  2. Oleinik O.A., Shamaev A.S., Yosifian G.A. Mathematical Problems in Elasticity and Homogenization. Amsterdam: North–Holland; 1992.

  3. Жиков В.В., Козлов С.М., Олейник О.А. Усреднение дифференциальных операторов. М.: Физматлит, 1993.

  4. Cioranescu D., Murat F. Un terme étrange venu d’ailleurs I & II. In Nonlinear Partial Differential Equations and their Applications. Collège de France Seminar, Volume II & III, ed. H. Berzis, J.L. Lions. Research Notes in Mathematics, 60 & 70, London: Pitman, 98–138 & 154–178; 1982.

  5. Беляев А.Г., Пятницкий А.Л., Чечкин Г.А. Усреднение в перфорированной области с осциллирующим третьим краевым условием // Матем. сб. 2001. Т. 192. № 7. С. 3– 20.

  6. Бабин А.В., Вишик М.И. Аттракторы эволюционных уравнений. М.: Наука, 1989.

  7. Chepyzhov V.V., Vishik M.I. Attractors for equations of mathematical physics. Providence (RI): Amer. Math. Soc.; 2002.

  8. Temam R. Infinite-dimensional dynamical systems in mechanics and physics. Applied Mathematics Series. V. 68. New York (NY): Springer-Verlag; 1988.

  9. Bekmaganbetov K.A., Chechkin G.A., Chepyzhov V.V. “Strange Term” in Homogenization of Attractors of Reaction–Diffusion Equation in Perforated Domain // Chaos, Solitons & Fractals. 2020. V. 140.

  10. Вишик М.И., Чепыжов В.В. Усреднение траекторных аттракторов эволюционных уравнений с быстро осциллирующими членами // Математический сборник. 2001. Т. 192. № 1. С. 13–50.

Дополнительные материалы отсутствуют.

Инструменты

Доклады Российской академии наук. Математика, информатика, процессы управления