Журнал вычислительной математики и математической физики, 2022, T. 62, № 6, стр. 965-976

О нелинейной задаче для системы интегродифференциальных уравнений теории переноса излучения

А. В. Калинин 12*, А. А. Тюхтина 1**

1 ННГУ
603022 Нижний Новгород, пр-т Гагарина, 23, Россия

2 ИПФ РАН
603950 Нижний Новгород, ул. Ульянова, 46, Россия

* E-mail: avk@mm.unn.ru
** E-mail: kalinmm@yandex.ru

Поступила в редакцию 12.12.2021
После доработки 20.01.2022
Принята к публикации 11.02.2022

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

Аннотация

Рассматривается начально-краевая задача для системы нелинейных интегродифференциальных уравнений теории переноса излучения. Доказывается теорема о существовании и единственности решения поставленной задачи. На основании свойств полугрупп изотонных операторов, действующих в условно-полных решетках, устанавливается стабилизация решения задачи при $t \to \infty $. Библ. 17.

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

1. ВВЕДЕНИЕ

Исследование нелинейных процессов переноса излучения приводит к необходимости изучения краевых и начально-краевых задач для нелинейных интегродифференциальных уравнений с частными производными. Вопросы физического, математического и численного моделирования процессов переноса излучения рассматриваются, в частности, в [1]–[7].

В работе [8] изучается нелинейная стационарная система, включающая кинетическое уравнение переноса излучения и уравнения статистического равновесия, возникающие при исследовании модели двухуровневого атома в предположении полного перераспределения излучения по частоте [1]–[3]. Получены строгие результаты о существовании и единственности решения краевой задачи, предложен и обоснован линеаризирующий итерационный алгоритм ее решения.

Соответствующая нестационарная нелинейная система интегродифференциальных уравнений переноса излучения изучалась в работах [9]–[13]. Обсуждались вопросы корректности постановки смешанной задачи для рассматриваемой системы, была установлена стабилизация при $t \to \infty $ решения задачи при произвольных начальных условиях.

В соответствии с физическим смыслом, решения уравнений теории переноса излучения должны быть неотрицательными функциями. В этом случае естественно возникают упорядоченные функциональные пространства, определяемые конусом неотрицательных функций, и при изучении задач возможно применение теории порядковых структур. Как показано в работе [13], вопросы порядковой и метрической стабилизации решений нестационарных задач могут исследоваться на основе общих подходов, связанных со свойствами полугрупп изотонных операторов, действующих в полных и условно полных решетках [14]–[17].

В настоящей работе развиваются результаты работ [9]–[13]. Рассматривается система интегродифференциальных уравнений, содержащая нестационарное уравнение переноса излучения и систему стационарности для модели двухуровнего атома. Изучается корректность постановки начально-краевой задачи для этой системы, исследуется поведение решения задачи при $t \to \infty $ с использованием методов теории упорядоченных пространств.

2. ПОСТАНОВКА ЗАДАЧИ И ОСНОВНЫЕ РЕЗУЛЬТАТЫ

Рассмотрим следующую систему нелинейных интегродифференциальных уравнений теории переноса излучения, соответствующую модели двухуровнего атома [3], [4]:

(2.1)
$\begin{gathered} \frac{1}{c}\frac{{\partial \psi (x,\nu ,\omega ,t)}}{{\partial t}} + (\omega ,\nabla )\psi (x,\nu ,\omega ,t) + \\ + \;h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}\left( {{{B}_{{12}}}{{C}_{1}}(x,t) - {{B}_{{21}}}{{C}_{2}}(x,t)} \right)\psi (x,\nu ,\omega ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}{{A}_{{21}}}{{C}_{2}}(x,t), \\ \end{gathered} $
(2.2)
$\begin{gathered} \left( {{{C}_{{12}}}{{n}_{e}}(x) + {{B}_{{12}}}\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \psi (x,\nu ,\omega ,t)d\omega d\nu } \right){{C}_{1}}(x,t) = \\ = \;{{C}_{2}}(x,t)\left( {{{A}_{{21}}} + {{C}_{{21}}}{{n}_{e}}(x) + {{B}_{{21}}}\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \psi (x,\nu ,\omega ,t)d\omega d\nu } \right), \\ \end{gathered} $
(2.3)
${{C}_{1}}(x,t) + {{C}_{2}}(x,t) = f(x).$
Здесь $x = \{ {{x}_{1}},{{x}_{2}},{{x}_{3}}\} \in G \subset {{\mathbb{R}}^{3}}$, $\omega \in \Omega = \{ \omega = \{ {{\omega }_{1}},{{\omega }_{2}},{{\omega }_{3}}\} \in {{\mathbb{R}}^{3}}:\omega _{1}^{2} + \omega _{2}^{2} + \omega _{3}^{2} = 1\} $; $\nu \in I = [0,{{\nu }_{0}}]$; $t \in [0,T]$, $T > 0$.

Функция $\psi $ – удельная интенсивность излучения, ${{C}_{1}}$ и ${{C}_{2}}$ – пространственные плотности атомов среды, находящихся в основном и в возбужденном состоянии соответственно. Подробная информация о физическом смысле функций и коэффициентов приводится в работах [3], [4].

Пусть $G$ – выпуклое ограниченное множество с гладкой границей $\partial G$, имеющей всюду внешнюю нормаль $n(x)$ ($x \in \partial G$), и диаметром $d > 0$. Система (2.1)–(2.3) дополняется граничным условием

(2.4)
$\psi (x,\nu ,\omega ,t) = 0,\quad x \in \partial G,\;(\omega ,n(x)) < 0,\quad \omega \in \Omega ,\quad \nu \in I,\quad t \in [0,T],$
соответствующим отсутствию внешнего потока частиц, падающего на границу области, и начальным условием

(2.5)
$\psi (x,\nu ,\omega ,0) = {{\psi }^{0}}(x,\nu ,\omega ).$

Предполагается, что $h$, ${{\nu }_{{12}}}$, ${{B}_{{12}}}$, ${{B}_{{21}}}$, ${{C}_{{12}}}$, ${{C}_{{21}}}$, ${{A}_{{21}}}$, ${{\nu }_{0}}$, $c$ – заданные положительные числа, удовлетворяющие условию

(2.6)
${{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}} > 0;$
${{n}_{e}}(x)$, $f(x)$, $x \in G$, $\kappa (\nu )$, $\nu \in I$ – заданные функции, измеримые и неотрицательные почти всюду в своих областях определения, удовлетворяющие условиям

(2.7)
$\operatorname{esssup} {{n}_{e}} = n_{e}^{ * } < \infty ,\quad \operatorname{esssup} f = f* < \infty ,\quad \operatorname{esssup} \kappa (\nu ) = \kappa * < \infty ,\quad \int\limits_I \kappa (\nu )d\nu = 1.$

Определим характеристику $\{ {{l}_{\omega }}\} $ дифференциального оператора $\partial {\text{/}}c\partial t + (\omega ,\nabla )$ системой уравнений

$cdt = \frac{{d{{x}_{1}}}}{{{{\omega }_{1}}}} = \frac{{d{{x}_{2}}}}{{{{\omega }_{2}}}} = \frac{{d{{x}_{3}}}}{{{{\omega }_{3}}}} = \frac{{d{{\omega }_{1}}}}{0} = \frac{{d{{\omega }_{2}}}}{0} = \frac{{d{{\omega }_{3}}}}{0} = \frac{{d\nu }}{0}$
и обозначим через $1{\text{/}}c{{\left( {d{\text{/}}d\tau } \right)}_{\omega }}$ оператор дифференцирования вдоль характеристики $\{ {{l}_{\omega }}\} $:

$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\psi (\xi ,\nu ,\omega ,t) = \frac{1}{c}{{\left. {\frac{{d\psi (x + c\omega (\tau - t),\nu ,\omega ,\tau )}}{{d\tau }}} \right|}_{{\tau = t}}}.$

Пусть $\mathcal{D} = G \times I \times \Omega $. Для произвольного измеримого подмножества $\Pi $ евклидова пространства через ${{K}_{\infty }}(\Pi )$ обозначим конус неотрицательных функций в ${{L}_{\infty }}(\Pi )$; ${{D}_{\infty }}(\mathcal{D} \times [0,T])$ – класс функций $\psi \in {{L}_{\infty }}(\mathcal{D} \times [0,T])$, абсолютно непрерывных вдоль почти каждой характеристики $\{ {{l}_{\omega }}\} $ в ${{D}_{\infty }}(\mathcal{D} \times [0,T])$ и таких, что

$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\psi \in {{L}_{\infty }}(\mathcal{D} \times [0,T]),$
${{\mathfrak{M}}_{T}} = {{D}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])),$
${{\mathcal{K}}_{T}} = {{K}_{\infty }}(\mathcal{D} \times [0,T]) \times {{K}_{\infty }}(G \times [0,T]) \times {{K}_{\infty }}(G \times [0,T])$.

Аналогичные классы впервые были введены в работе В.С. Владимирова [6] для стационарных задач теории переноса и использовались при изучении нестационарных задач в [9], [13].

Предполагается, что начальная функция ${{\psi }_{0}}(x,\nu ,\omega )$ принадлежит классу ${{K}_{\infty }}(\mathcal{D})$.

Решением задачи (2.1)–(2.5) называется функция $\Phi = \{ \psi ,{{C}_{1}},{{C}_{2}}\} \in {{\mathfrak{M}}_{T}}$, удовлетворяющая системе (2.1)–(2.3) и начальным и краевым условиям (2.4), (2.5) почти всюду, дифференциальный оператор $d{\text{/}}cdt + (\omega ,\nabla )$ в (2.1) понимается как оператор $1{\text{/}}c{{\left( {d{\text{/}}d\tau } \right)}_{\omega }}$ дифференцирования вдоль характеристики $\{ {{l}_{\omega }}\} $.

Теорема 1. Пусть коэффициенты системы удовлетворяют условиям (2.6), (2.7). Тогда при любой начальной функции ${{\psi }_{0}} \in {{K}_{\infty }}(\mathcal{D})$ и любом $T > 0$ решение $\Phi \in {{\mathfrak{M}}_{T}}$ задачи (2.1)–(2.5) существует, единственно и непрерывно зависит от начальной функции ${{\psi }_{0}}$. Кроме того, имеет место включение $\Phi \in {{\mathfrak{M}}_{T}} \cap {{\mathcal{K}}_{T}}$.

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

(2.8)
${{\mathcal{R}}_{{12}}}(\psi )(x,t) = {{C}_{{12}}}{{n}_{e}}(x) + {{B}_{{12}}}J(\psi )(x,t),\quad {{\mathcal{R}}_{{21}}}(\psi )(x,t) = {{A}_{{21}}} + {{C}_{{21}}}{{n}_{e}}(x) + {{B}_{{21}}}J(\psi )(x,t),$
(2.9)
$J(\psi )(x,t) = \int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \psi (x,\nu ,\omega ,t){\kern 1pt} d\omega d\nu ,\quad \mathcal{R}(\psi ) = {{\mathcal{R}}_{{12}}}(\psi )(x,t) + {{\mathcal{R}}_{{21}}}(\psi )(x,t),$
(2.10)
$\mathcal{F}(C)(x,\nu ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}({{B}_{{12}}}{{C}_{1}}(x,t) - {{B}_{{21}}}{{C}_{2}}(x,t)),$
(2.11)
$\mathcal{P}(C)(x,\nu ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}{{A}_{{21}}}{{C}_{2}}(x,t).$

В работах [9], [13] подобные обозначения использовались для исследования других классов нестационарных задач теории переноса излучения. Систему (2.1)–(2.3) можно переписать в виде

(2.12)
$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\psi (x,\nu ,\omega ,t) + \mathcal{F}(C)(x,\nu ,t)\psi (x,\nu ,\omega ,t) = \mathcal{P}(C)(x,\nu ,t),$
(2.13)
${{C}_{1}}(x,t) = f(x)\frac{{{{\mathcal{R}}_{{21}}}(\psi )(x,t)}}{{\mathcal{R}(\psi )(x,t)}},\quad {{C}_{2}}(x,t) = f(x)\frac{{{{\mathcal{R}}_{{12}}}(\psi )(x,t)}}{{\mathcal{R}(\psi )(x,t)}}.$

Функция $\Phi = \{ \psi ,{{C}_{1}},{{C}_{2}}\} \in {{\mathfrak{M}}_{T}}$ является решением задачи (2.1)–(2.5) тогда и только тогда, когда $\psi \in {{D}_{\infty }}(\mathcal{D} \times [0,T])$ – решение задачи (2.12), (2.4), (2.5), где ${{C}_{1}}$, ${{C}_{2}}$ в (2.12) определяются по формулам (2.13).

Дифференциальные свойства функционального класса ${{D}_{\infty }}(\mathcal{D} \times [0,T])$ позволяют определить для любого $T > 0$ однопараметрические семейства операторов

${{\gamma }_{t}}:{{D}_{\infty }}(\mathcal{D} \times [0,T]) \to {{L}_{\infty }}(\mathcal{D}),\quad t \in [0,T],$
имеющие смысл “следа” функции $\psi \in {{D}_{\infty }}(\mathcal{D} \times [0,T])$ на множестве $\mathcal{D} \times \{ t\} $.

Если $\psi $ – решение задачи (2.12), (2.4), (2.5) из класса ${{D}_{\infty }}(\mathcal{D} \times [0,T]) \cap {{K}_{\infty }}(\mathcal{D} \times [0,T])$, то значение следа ${{\gamma }_{t}}\psi $ не зависит от $T \geqslant t$ и лежит в ${{K}_{\infty }}(\mathcal{D})$ при всех $t \in [0,T]$.

Определим семейство $\{ {{U}_{t}}\} $, $0 \leqslant t \leqslant \infty $, операторов ${{U}_{t}}:{{K}_{\infty }}(\mathcal{D}) \to {{K}_{\infty }}(\mathcal{D})$ соотношением ${{U}_{t}}{{\psi }_{0}} = {{\gamma }_{t}}\psi $ для всех ${{\psi }_{0}} \in {{K}_{\infty }}(\mathcal{D})$, где $\Phi = \{ \psi ,{{C}_{1}},{{C}_{2}}\} \in {{\mathfrak{M}}_{T}} \cap {{\mathcal{K}}_{T}}$ – решение задачи (2.1), (2.2)–(2.5) с начальным условием ${{\psi }_{0}} \in {{K}_{\infty }}(\mathcal{D})$. Введенное семейство операторов в силу независимости коэффициентов системы (2.1), (2.2), (2.3) от времени обладает полугрупповыми свойствами и называется разрешающей полугруппой задачи (2.1)–(2.5).

Множество ${{K}_{\infty }}(\mathcal{D})$ является условно полной решеткой [14] с отношением частичного порядка $ \succ $, где ${{\psi }_{1}} \succ {{\psi }_{2}}$ тогда и только тогда, когда ${{\psi }_{1}}(x,\nu ,\omega ) \geqslant {{\psi }_{2}}(x,\nu ,\omega )$ почти всюду в $\mathcal{D}$. Формула

(2.14)
$E[\psi ] = \frac{1}{c}\int\limits_I {\int\limits_\Omega {\int\limits_G \psi } } (x,\nu ,\omega )dxd\omega d\nu ,\quad \psi \in {{K}_{\infty }}(\mathcal{D}),$
определяет функционал $E:{{K}_{\infty }}(\mathcal{D}) \to {{\mathbb{R}}^{1}}$, удовлетворяющий условиям
$E[{{\psi }_{1}}] + E[{{\psi }_{2}}] = E[{{\psi }_{1}} \vee {{\psi }_{2}}] + E[{{\psi }_{1}} \wedge {{\psi }_{2}}];$
$E[{{\psi }_{1}}] > E[{{\psi }_{2}}]$ при ${{\Phi }_{1}} \succ {{\Phi }_{2}}$, ${{\Phi }_{1}} \ne {{\Phi }_{2}}$, и является, поэтому, положительной оценкой на ${{K}_{\infty }}(\mathcal{D})$ (см. [14]), превращающей ${{K}_{\infty }}(D)$ в метрическое пространство с функцией расстояния

$\rho ({{\psi }_{1}},{{\psi }_{2}}) = E[{{\psi }_{1}} \vee {{\psi }_{2}}] - E[{{\psi }_{1}} \wedge {{\psi }_{2}}] = \frac{1}{c}{{\left\| {{{\psi }_{1}} - {{\psi }_{2}}} \right\|}_{{{{L}_{1}}(\mathcal{D})}}}.$

Справедлива

Теорема 2. Пусть коэффициенты системы (2.1), (2.2), (2.3) удовлетворяют условиям (2.6), (2.7). Тогда существует единственная стационарная траектория ${{\psi }_{\infty }}(t) = {{\psi }_{\infty }} \in {{K}_{\infty }}(\mathcal{D})$, $0 \leqslant t < \infty $, относительно разрешающей полугруппы операторов $\{ {{U}_{t}}\} $ задачи (2.1), (2.2)–(2.5). Для любых $\psi \in {{K}_{\infty }}(\mathcal{D})$ справедлива оценка

$\rho ({{U}_{t}}\psi ,{{\psi }_{\infty }}) \leqslant \gamma \exp \{ - \mu t\} $
при $t \geqslant {{T}_{0}}$, где положительные постоянные $\gamma $, $\mu $, ${{T}_{0}}$ не зависят от выбора $\psi \in {{K}_{\infty }}(\mathcal{D})$.

3. ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 1

Положим ${{\varphi }_{0}}(x,\nu ,\omega ,t) = {{e}^{{ - \alpha t}}}\psi (x,\nu ,\omega ,t)$, где $\alpha > 0$, $\varphi = \{ {{\varphi }_{0}},{{C}_{1}},{{C}_{2}}\} $.

Система (2.12), (2.13) примет вид

(3.1)
$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}{{\varphi }_{0}}(x,\nu ,\omega ,t) + (\alpha + \mathcal{F}(C)(x,\nu ,t)){{\varphi }_{0}}(x,\nu ,\omega ,t) = {{e}^{{ - \alpha t}}}\mathcal{P}(C)(x,\nu ,t),$
(3.2)
${{C}_{1}}(x,t) = f(x)\frac{{{{\mathcal{R}}_{{21}}}({{e}^{{\alpha t}}}{{\varphi }_{0}})(x,t)}}{{\mathcal{R}({{e}^{{\alpha t}}}{{\varphi }_{0}})(x,t)}},\quad {{C}_{2}}(x,t) = f(x)\frac{{{{\mathcal{R}}_{{12}}}({{e}^{{\alpha t}}}{{\varphi }_{0}})(x,t)}}{{\mathcal{R}({{e}^{{\alpha t}}}{{\varphi }_{0}})(x,t)}}.$
Система (3.1), (3.2) рассматривается при начальных и граничных условиях

(3.3)
${{\varphi }_{0}}(x,\nu ,\omega ,0) = {{\psi }^{0}}(x,\nu ,\omega ),\quad {{\varphi }_{0}}(x,\nu ,\omega ,t) = 0,\quad x \in \partial G,\quad (\omega ,n(x)) < 0.$

Пусть $\{ {{l}_{\omega }}\} = \{ (x + c\omega (\tau - t),\nu ,\omega ,\tau ),\;\tau \in \mathbb{R}\} $ – характеристика, проходящая в момент времени $\tau = t$ через точку $(x,\nu ,\omega )$. Обозначим через $t_{\omega }^{ - }(x,t)$ момент пересечения характеристикой границы области $\mathcal{D} \times [0,T]$ такой, что либо $t_{\omega }^{ - }(x,t) = 0$, либо $t_{\omega }^{ - }(x,t)$ соответствует пересечению той части $\partial G \times I \times \Omega \times [0,T]$, где $(\omega ,n(x)) < 0$. Тогда

${{g}_{0}}(x,\nu ,\omega ,t) = {{\varphi }_{0}}(x + c\omega (t_{\omega }^{ - }(x,t) - t),\;\nu ,\omega ,t_{\omega }^{ - }(x,t)) \equiv \left\{ \begin{gathered} 0,\quad t_{\omega }^{ - }(x,t) > 0, \hfill \\ {{\psi }_{0}}(x - c\omega t,\nu ,\omega ),\quad t_{\omega }^{ - }(x,t) = 0. \hfill \\ \end{gathered} \right.$

Разрешая уравнение (3.1) как линейное относительно ${{\varphi }_{0}}(x,\nu ,\omega ,t)$, получаем

(3.4)
$\begin{gathered} {{\varphi }_{0}}(x,\nu ,\omega ,t) = {{g}_{0}}(x,\nu ,\omega ,t)\exp \left\{ { - \int\limits_{t_{\omega }^{ - }(x,t)}^t {\left( {\alpha + c\mathcal{F}(C)(x + c\omega (s - t),\nu ,s)} \right)ds} } \right\} + \\ + \int\limits_{t_{\omega }^{ - }(x,t)}^t {{e}^{{ - \alpha \tau }}}c\mathcal{P}(C)(x + c\omega (\tau - t),\nu ,\tau )\exp \left\{ {\int\limits_t^\tau {\left( {\alpha + c\mathcal{F}(C)(x + c\omega (s - t),\nu ,s)} \right)} ds} \right\}d\tau . \\ \end{gathered} $

Система уравнений (3.2), (3.4) эквивалентна при заданных условиях на коэффициенты и на рассматриваемом классе функций $\varphi \in {{\mathfrak{M}}_{T}}$ системе с начальными и граничными условиями (3.1)–(3.3).

Пусть ${{\mathcal{L}}_{T}} = {{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])$. Введем операторы

(3.5)
$\begin{gathered} {{A}_{0}}:{{\mathcal{L}}_{T}} \to {{D}_{\infty }}(\mathcal{D} \times [0,T]),\quad {{A}_{i}}:{{\mathcal{L}}_{T}} \to {{\mathcal{L}}_{\infty }}(G \times [0,T]),\quad i = 1,\;2, \\ A:{{\mathcal{L}}_{T}} \to {{\mathfrak{M}}_{T}} \subset {{\mathcal{L}}_{T}},\quad \mathcal{L}\varphi = \{ {{A}_{0}}\varphi ,{{A}_{1}}\varphi ,{{A}_{2}}\varphi \} , \\ \end{gathered} $
операторы ${{A}_{0}}\varphi $, ${{A}_{i}}\varphi $ ($i = 1,\;2$) совпадают с правыми частями уравнений (3.4), (3.2) соответственно.

Положим

$M = \max \left\{ {{{{\left\| {{{\psi }_{0}}} \right\|}}_{{{{L}_{\infty }}(\mathcal{D})}}},{{{\left\| f \right\|}}_{{{{L}_{\infty }}(G)}}}} \right\},\quad \alpha = 2c(b + Ma),$
где $b = h{{\nu }_{{12}}}\kappa {\text{*}}{{A}_{{21}}}{{(4\pi )}^{{ - 1}}}$, $a = h{{\nu }_{{12}}}\kappa {\text{*}}({{B}_{{12}}} + {{B}_{{21}}})(4\pi {{)}^{{ - 1}}}$.

Лемма 1. Пусть ${{A}_{B}}$сужение оператора $A$ на замкнутое множество $B(\mathcal{D} \times [0,T])$,

$B(\mathcal{D} \times [0,T]) = \left\{ {\varphi = \{ \xi ,{{\eta }_{1}},{{\eta }_{2}}\} \in {{\mathcal{K}}_{T}}:{{{\left\| \xi \right\|}}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant M,\;{{\eta }_{1}} + {{\eta }_{2}} = f} \right\},$
$R({{A}_{B}})$ – множество значений оператора ${{A}_{B}}$. Тогда $R({{A}_{B}}) \subset B(\mathcal{D} \times [0,T])$.

Доказательство. Пусть $\varphi = \{ {{\varphi }_{0}},{{C}_{1}},{{C}_{2}}\} \in B(\mathcal{D} \times [0,T])$.

Так как

$\mathcal{F}(C)(x,\nu ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}\frac{{{{A}_{{21}}}{{B}_{{12}}} + ({{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}}){{n}_{e}}(x)}}{{\mathcal{R}({{e}^{{\alpha t}}}{{\varphi }_{0}})(x,t)}} \geqslant 0,$
из (3.2), (3.4) вытекает неотрицательность вектора $A\varphi $. Из (3.2) следует, что ${{A}_{1}}\varphi + {{A}_{2}}\varphi = f$. Далее, ввиду (3.4) и явного вида $\mathcal{F}(C)$, $\mathcal{P}(C)$ имеем
${{A}_{0}}{{\varphi }_{0}}(x,\nu ,\omega ,t) \leqslant \max \left\{ {M,{{{\left\| {\frac{{c{{e}^{{ - \alpha t}}}\mathcal{P}(C)}}{{\alpha + c\mathcal{F}(C)}}} \right\|}}_{{{{L}_{\infty }}(G \times I \times [0,T])}}}} \right\} \leqslant \max \left\{ {M,\frac{{cb}}{\alpha }M} \right\},$
т.е. ${{A}_{0}}{{\varphi }_{0}}(x,\nu ,\omega ,t) \leqslant M$ почти всюду.

Лемма 2. Существует единственное решение задачи

(3.6)
$\varphi (x,\nu ,\omega ,t) = {{A}_{B}}\varphi (x,\nu ,\omega ,t),\quad \varphi \in B(\mathcal{D} \times [0,T]).$

Доказательство. Покажем, что при достаточно больших $k$ отображение $A_{B}^{k} = \underbrace {{{A}_{B}} \cdot \ldots \cdot {{A}_{B}}}_B$ будет сжимающим на $B(\mathcal{D} \times [0,T])$. Пусть

${{\varphi }^{{(i)}}} = \left\{ {\varphi _{0}^{{(i)}},C_{1}^{{(i)}},C_{2}^{{(i)}}} \right\} \in B(\mathcal{D} \times [0,T]),\quad i = 1,\;2,$
и пусть ${{J}_{i}} = J\left( {{{e}^{{\alpha t}}}\varphi _{0}^{{(i)}}} \right)$, ${{\mathcal{R}}_{i}} = \mathcal{R}(\varphi _{0}^{{(i)}})$, $i = 1,\;2$. Используя (3.4) и явный вид ${{\mathcal{R}}_{{12}}}$, ${{\mathcal{R}}_{{21}}}$, получаем
${{A}_{1}}{{\varphi }^{{(1)}}} - {{A}_{1}}{{\varphi }^{{(2)}}} = f\frac{{({{A}_{{21}}}{{B}_{{12}}} + ({{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}}){{n}_{e}})({{J}_{2}} - {{J}_{1}})}}{{{{\mathcal{R}}_{1}}{{\mathcal{R}}_{2}}}},$
$\left| {{{A}_{1}}{{\varphi }^{{(1)}}} - {{A}_{1}}{{\varphi }^{{(2)}}}} \right| \leqslant Mp{{e}^{{\alpha t}}}\left| {\varphi _{0}^{{(1)}} - \varphi _{0}^{{(2)}}} \right|,$
где $p = ({{A}_{{21}}}{{B}_{{12}}} + ({{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}})n_{e}^{ * })A_{{21}}^{{ - 2}}$,

${{A}_{2}}{{\varphi }^{{(1)}}} - {{A}_{2}}{{\varphi }^{{(2)}}} = {{A}_{1}}{{\varphi }^{{(2)}}} - {{A}_{1}}{{\varphi }^{{(1)}}},\quad \left| {{{A}_{2}}{{\varphi }^{{(1)}}} - {{A}_{2}}{{\varphi }^{{(2)}}}} \right| \leqslant Mp{{e}^{{\alpha t}}}\left| {\varphi _{0}^{{(1)}} - \varphi _{0}^{{(2)}}} \right|.$

Пусть, далее, ${{\mathcal{F}}_{i}} = \mathcal{F}({{C}^{{(i)}}})$, ${{\mathcal{P}}_{i}} = \mathcal{P}({{C}^{{(i)}}})$. Из (3.4) вытекает справедливость равенства [15]

$\begin{gathered} {{A}_{0}}{{\varphi }^{{(1)}}}(x,\nu ,\omega ,t) - {{A}_{0}}{{\varphi }^{{(2)}}}(x,\nu ,\omega ,t) = \int\limits_{t_{\omega }^{ - }(x,t)}^t {(c{{e}^{{ - \alpha \tau }}}({{\mathcal{P}}_{1}} - {{\mathcal{P}}_{2}})} (x + c\omega (\tau - t),\nu ,\omega ,\tau ) + \\ + \;c{{A}_{0}}{{\varphi }^{{(2)}}}(x + c\omega (\tau - t),\nu ,\omega ,\tau )({{\mathcal{F}}_{1}} - {{\mathcal{F}}_{2}})(x + c\omega (\tau - t),\nu ,\omega ,\tau )) \times \\ \times \;\exp \left\{ {\int\limits_t^\tau {\left( {\alpha + c{{\mathcal{F}}_{1}}(x + c\omega (s - t),\nu ,s)} \right)} ds} \right\}d\tau . \\ \end{gathered} $
Следовательно,
$\begin{gathered} \left| {{{A}_{0}}{{\varphi }^{{(1)}}}(x,\nu ,\omega ,t) - {{A}_{0}}{{\varphi }^{{(2)}}}(x,\nu ,\omega ,t} \right| \leqslant \\ \leqslant {{\left\| {\frac{{c{{e}^{{ - \alpha t}}}({{\mathcal{P}}_{1}} - {{\mathcal{P}}_{2}}) + c{{A}_{0}}{{\varphi }^{{(2)}}}({{\mathcal{F}}_{1}} - {{\mathcal{F}}_{2}})}}{{\alpha + c{{\mathcal{F}}_{1}}}}} \right\|}_{{{{L}_{\infty }}(G \times I \times [0,T])}}} \leqslant c\frac{{b + Ma}}{\alpha }{{\left\| {C_{2}^{{(1)}} - C_{2}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}}. \\ \end{gathered} $
Таким образом,

${{\left\| {{{A}_{0}}{{\varphi }^{{(1)}}} - {{A}_{0}}{{\varphi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant \frac{1}{2}{{\left\| {C_{2}^{{(1)}} - C_{2}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}} \leqslant \frac{1}{2}{{\left\| {{{\varphi }^{{(1)}}} - {{\varphi }^{{(2)}}}} \right\|}_{{{{\mathcal{L}}_{T}}}}},$
${{\left\| {A_{0}^{k}{{\varphi }^{{(1)}}} - A_{0}^{k}{{\varphi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant \frac{1}{{{{2}^{k}}}}\parallel {{\varphi }^{{(1)}}} - {{\varphi }^{{(2)}}}{{\parallel }_{{{{\mathcal{L}}_{T}}}}},$
${{\left\| {A_{i}^{k}{{\varphi }^{{(1)}}} - A_{i}^{k}{{\varphi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}} \leqslant Mp{{e}^{{\alpha T}}}{{\left\| {A_{0}^{{k - 1}}{{\varphi }^{{(1)}}} - A_{0}^{{k - 1}}{{\varphi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant \frac{{Mp{{e}^{{\alpha T}}}}}{{{{2}^{{k - 1}}}}}{{\left\| {{{\varphi }^{{(1)}}} - {{\varphi }^{{(2)}}}} \right\|}_{{{{\mathcal{L}}_{T}}}}},$
${{\left\| {{{A}^{k}}{{\varphi }^{{(1)}}} - {{A}^{k}}{{\varphi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant \frac{1}{{{{2}^{k}}}}(1 + 4Mp{{e}^{{\alpha T}}}){{\left\| {{{\varphi }^{{(1)}}} - {{\varphi }^{{(2)}}}} \right\|}_{{{{\mathcal{L}}_{T}}}}}.$

Согласно принципу сжимающих отображений, уравнение (3.6) имеет единственное решение $\varphi \in B(\mathcal{D} \times [0,T])$.

Поскольку $A:{{\mathcal{L}}_{T}} \to {{\mathfrak{M}}_{T}}$, решение задачи (3.6) лежит в ${{\mathfrak{M}}_{T}}$ и является решением задачи (3.1)–(3.3). Следовательно, решение задачи (2.1)–(2.5) существует.

Лемма 3. Решение задачи (2.1)–(2.5) единственное и непрерывно зависит от начальных условий.

Доказательство. Пусть ${{\Phi }^{{(i)}}} = \left\{ {{{\psi }^{{(i)}}},C_{1}^{{(i)}},C_{2}^{{(i)}}} \right\} \in {{\mathfrak{M}}_{T}}$, $i = 1,\;2$ – решения задачи (2.1)–(2.5) c начальными условиями

${{\psi }^{{(i)}}}(x,\nu ,\omega ,0) = \psi _{0}^{{(i)}}(x,\nu ,\omega ),\quad i = 1,\;2,$
причем ${{\Phi }^{{(1)}}} \in {{\mathcal{K}}_{T}}$. Положим $\psi = {{\psi }^{{(1)}}} - {{\psi }^{{(2)}}}$, $C = C_{2}^{{(1)}} - C_{2}^{{(2)}} = C_{1}^{{(2)}} - C_{1}^{{(1)}}$, ${{\psi }_{0}} = \psi _{0}^{{(1)}} - \psi _{0}^{{(2)}}$.

Тогда справедливы равенства

(3.7)
$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\psi (x,\nu ,\omega ,t) + \mathcal{F}({{C}^{1}})(x,\nu ,t)\psi (x,\nu ,\omega ,t) = K(x,\nu ,\omega ,t)C(x,t),$
(3.8)
$C(x,t) = H(x,t)\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \psi (x,\nu ,\omega ,t)d\omega d\nu ,$
где

$K(x,\nu ,\omega ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}({{A}_{{21}}} + ({{B}_{{12}}} + {{B}_{{21}}}){{\psi }^{{(2)}}}(x,\nu ,\omega ,t)),$
$H(x,t) = \frac{{{{B}_{{12}}}f(x) - ({{B}_{{12}}} + {{B}_{{21}}})C_{2}^{{(2)}}(x,t)}}{{\mathcal{R}({{C}^{{(1)}}})(x,t)}}.$

Проинтегрировав равенство (3.7), получим

$\begin{gathered} \psi (x,\nu ,\omega ,t) - {{\psi }_{0}}(x,\nu ,\omega ) + c\int\limits_{t_{\omega }^{ - }(x,t)}^t \mathcal{F} ({{C}^{{(1)}}})(x + c\omega (\tau - t),\nu ,\tau )\psi (x + c\omega (\tau - t),\nu ,\omega ,\tau )d\tau = \\ = \;c\int\limits_{t_{\omega }^{ - }(x,t)}^t K (x + c\omega (\tau - t),\tau )C(x + c\omega (\tau - t),\tau )d\tau . \\ \end{gathered} $

Пусть $\Delta (t) = {{\left\| {\psi (x,\nu ,\omega ,t)} \right\|}_{{{{L}_{\infty }}(\mathcal{D})}}}$. Тогда

$\Delta (t) \leqslant \Delta (0) + \beta \int\limits_{t_{\omega }^{ - }(x,t)}^t \Delta (\tau )d\tau ,$
где

$\beta = c\left( {{{{\left\| {\mathcal{F}({{C}^{1}})} \right\|}}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} + {{{\left\| H \right\|}}_{{{{L}_{\infty }}(G \times [0,T])}}}{{{\left\| K \right\|}}_{{{{L}_{\infty }}(G \times [0,T])}}}} \right).$

По лемме Гронуолла заключаем, что

(3.9)
${{\left\| {{{\psi }^{{(1)}}} - {{\psi }^{{(2)}}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant {{e}^{{\beta T}}}{{\left\| {\psi _{0}^{{(1)}} - \psi _{0}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D})}}},$
(3.10)
${{\left\| {C_{2}^{{(1)}} - C_{2}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}} = {{\left\| {C_{2}^{{(1)}} - C_{2}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}} \leqslant {{e}^{{\beta T}}}{{\left\| {\psi _{0}^{{(1)}} - \psi _{0}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D})}}}{{\left\| H \right\|}_{{{{L}_{\infty }}(G \times [0,T])}}}.$

Из (3.9), (3.10) вытекает, в частности, единственность решения задачи (2.1)–(2.5). Кроме того, из (3.7) получаем

${{\left\| {{{{\left( {\frac{d}{{d\tau }}} \right)}}_{\omega }}({{\psi }^{{(1)}}} - {{\psi }^{{(2)}}})} \right\|}_{{{{L}_{\infty }}(\mathcal{D})}}} \leqslant \beta {{e}^{{\beta T}}}{{\left\| {\psi _{0}^{{(1)}} - \psi _{0}^{{(2)}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D})}}}.$
Таким образом, решение задачи (2.1)–(2.5) непрерывно зависит от начальных данных.

4. ОБЩИЕ СВОЙСТВА ПОЛУГРУПП ИЗОТОННЫХ ОПЕРАТОРОВ

Доказательство теоремы о стабилизации решений рассматриваемой смешанной задачи основано на свойствах полугрупп изотонных операторов, действующих в условно полных решетках. Приведем без доказательств необходимые в дальнейшем результаты. Подробное изложение рассматриваемых вопросов содержится в [12], [13].

Пусть $\mathfrak{M}$ – непустое частично упорядоченное множество с отношением порядка $ \prec $, являющееся полной решеткой [14]; $U = \{ {{U}_{t}}\} $, $t \geqslant 0$ – однопараметрическое семейство изотонных операторов, обладающее полугрупповыми свойствами:

${{U}_{0}} \equiv E - {\text{тождественный}}\;{\text{оператор}},$
${{U}_{{t + s}}} = {{U}_{t}}{{U}_{s}},$
${{U}_{t}}u \prec {{U}_{t}}v\quad {\text{при}}\quad u \prec v\quad {\text{для}}\;{\text{всех}}\quad t,s \geqslant 0,\quad u,v \in \mathfrak{M};$
здесь  $U$-траекторией  называется функция $u:[0,\infty ) \to \mathfrak{M}$ такая, что $u(t + s) = {{U}_{t}}u(s)$, $t,s \geqslant 0$; $U$-траектория называется стационарной, если

(4.1)
$u(t) = u(0),\quad t \geqslant 0.$

Элемент $u(0) \in \mathfrak{M}$, удовлетворяющий (4.1) при всех $t \geqslant 0$, называется стационарной точкой полугруппы U.

Теорема 3. Пусть $U = \{ {{U}_{t}}\} $, $t \geqslant 0$, – полугруппа изотонных операторов ${{U}_{t}}:\mathfrak{M} \to \mathfrak{M}$. Тогда множество стационарных траекторий $u:[0,\infty ) \to \mathfrak{M}$ и множество $P = \{ \xi \in \mathfrak{M}:{{U}_{t}}\xi = \xi ,\;t \geqslant 0\} $ соответствующих стационарных точек полугруппы $U$ непусто и $\inf P \in P$, $\sup P \in P$.

Теорема является прямым следствием следующего утверждения [17].

Теорема 4. Пусть $\mathfrak{M}$полная решетка, $\mathfrak{U}$ – семейство изотонных коммутирующих операторов $A:\mathfrak{M} \to \mathfrak{M}$. Тогда множество $P$ общих неподвижных точек семейства операторов $\mathfrak{U}$ непусто и является полной подрешеткой решетки $\mathfrak{M}$.

Изотонный оператор $A:\mathfrak{M} \to \mathfrak{M}$ будем называть оператором класса $C(\mathfrak{M}, \prec )$, если для любых последовательностей

(4.2)
${{\xi }_{1}} \prec {{\xi }_{2}} \prec \; \ldots \; \prec {{\xi }_{n}} \prec \; \ldots \;{{\eta }_{1}} \succ {{\eta }_{2}} \succ \; \ldots \; \succ {{\eta }_{n}} \succ \; \ldots $
элементов из $\mathfrak{M}$ выполняются равенства
$\mathop {\sup }\limits_{n \in N} {\kern 1pt} \{ A{{\xi }_{n}}\} = A{\kern 1pt} \mathop {\sup }\limits_{n \in N} {\kern 1pt} \{ {{\xi }_{n}}\} ,\quad \mathop {\inf }\limits_{n \in N} {\kern 1pt} \{ A{{\eta }_{n}}\} = A\mathop {\inf }\limits_{n \in N} \{ {{\eta }_{n}}\} .$
Полугруппа $U = \{ {{U}_{t}}\} $, $t \geqslant 0$, изотонных операторов ${{U}_{t}}:\mathfrak{M} \to \mathfrak{M}$ называется полугруппой класса $\zeta (\mathfrak{M}, \prec )$, если ${{U}_{t}} \in C(\mathfrak{M}, \prec )$ при всех $t \geqslant 0$.

Теорема 5. Пусть $U = \{ {{U}_{t}}\} $, $t \geqslant 0$, – полугруппа класса $\zeta (\mathfrak{M}, \prec )$; $P \in \mathfrak{M}$ – множество стационарных точек. Тогда порядковый интервал

$\left\langle {\inf P,\sup P} \right\rangle = \{ \xi \in \mathfrak{M}:\inf P \prec \xi \prec \sup P\} $
является притягивающим множеством для всех траекторий, т.е.
$\inf P \prec \mathop {\lim }\limits_{t \to \infty } \inf \{ {{U}_{t}}\xi \} \prec \mathop {\lim }\limits_{t \to \infty } \sup \{ {{U}_{t}}\xi \} \prec \sup P$
при всех $\xi \in \mathfrak{M}$, в частности, если множество $P$ стационарных точек полугруппы $U$ coстоит из одного элемента ${{u}_{\infty }} = \inf P = \sup P$, то при $t \to \infty $ имеет место порядковая стабилизация всех траекторий:

$\mathop {\lim }\limits_{t \to \infty } \{ {{U}_{t}}\xi \} = \mathop {\lim }\limits_{t \to \infty } \sup \{ {{U}_{t}}\xi \} = \mathop {\lim }\limits_{t \to \infty } \inf \{ {{U}_{t}}\xi \} = {{u}_{\infty }}\quad {\text{при}}\;{\text{всех}}\quad \xi \in \mathfrak{M}.$

Положительной оценкой на полной решетке $\mathfrak{M}$ называется функционал $l:\mathfrak{M} \to {{\mathbb{R}}^{1}}$, удовлетворяющий при всех $\xi ,\eta \in \mathfrak{M}$ условиям

$l(\xi ) + l(\eta ) = l(\sup \{ \xi ,\eta \} ) + l(\inf \{ \xi ,\eta \} ),\quad l(\xi ) < l(\eta )\quad {\text{при}}\quad \xi \prec \eta ,\quad \xi \ne \eta .$
Положительный функционал $l:\mathfrak{M} \to {{\mathbb{R}}^{1}}$ превращает решетку $\mathfrak{M}$ в метрическую решетку с функцией расстояния $d(\xi ,\eta ) = l(\sup \{ \xi ,\eta \} ) - l(\inf \{ \xi ,\eta \} )$ [17].

Если положительная оценка $l:\mathfrak{M} \to {{\mathbb{R}}^{1}}$ удовлетворяет следующему условию непрерывности: для любых последовательностей $\{ {{\xi }_{n}}\} _{{n = 1}}^{\infty }$, $\{ {{\eta }_{n}}\} _{{n = 1}}^{\infty }$ элементов из $\mathfrak{M}$, удовлетворяющих (4.2), выполняются равенства

(4.3)
$\mathop {\lim }\limits_{t \to \infty } l({{\xi }_{n}}) = l\left( {\mathop {\sup }\limits_{n \in N} \{ {{\xi }_{n}}\} } \right),\quad \mathop {\lim }\limits_{t \to \infty } l({{\eta }_{n}}) = l\left( {\mathop {\inf }\limits_{n \in N} \{ {{\eta }_{n}}\} } \right),$
то метрическая решетка $\mathfrak{M}$ является метрически полной.

Оператор $A:\mathfrak{M} \to \mathfrak{M}$ является оператором класса $C(\mathfrak{M}, \prec )$ тогда и только тогда, когда $A$ метрически непрерывен.

Теорема 6. Пусть $U = \{ {{U}_{t}}\} $, $t \geqslant 0$, – полугруппа изотонных операторов, $P \subset \mathfrak{M}$ – множество стационарных точек полугруппы, $l:\mathfrak{M} \to {{\mathbb{R}}^{1}}$ – положительная оценка, удовлетворяющая условию непрерывности (4.3), и для каждого $t \geqslant 0$ оператор ${{U}_{t}}$ : $\mathfrak{M} \to \mathfrak{M}$ метрически непрерывен. Тогда порядковый интервал, т.е. $\left\langle {\inf P,\sup P} \right\rangle $ является метрически притягивающим множеством для всех траекторий, т.е.

$\mathop {\lim }\limits_{t \to \infty } d\left( {{{U}_{t}}\xi ,{\kern 1pt} \left\langle {\inf P,\sup P} \right\rangle } \right) = 0$
при всех $\xi \in \mathfrak{M}$, в частности, если множество стационарных точек полугруппы $U$ состоит из одного элемента ${{u}_{\infty }} = \inf P = \sup P$, то при $t \to \infty $ имеет место метрическая стабилизация всех траекторий: $\mathop {\lim }\limits_{t \to \infty } d({{U}_{t}}\xi ,{{u}_{\infty }}) = 0$ при всех $\xi \in \mathfrak{M}$.

Теорема 7. Пусть $\mathfrak{M}$условно полная решетка, ${{\mathfrak{M}}_{0}} \subset \mathfrak{M}$ – полная подрешетка в $\mathfrak{M}$; $l:\mathfrak{M} \to \mathfrak{M}$ – положительная оценка на $\mathfrak{M}$, $U = \{ {{U}_{t}}\} $, $t \geqslant 0$, – полугруппа операторов, удовлетворяющая следующим условиям:

1) существует ${{T}_{1}} > 0$ такое, что ${{U}_{t}}\xi \in {{\mathfrak{M}}_{0}}$ при всех $\xi \in \mathfrak{M}$, $t \geqslant {{T}_{1}}$;

2) при любых $\xi ,\eta \in {{\mathfrak{M}}_{0}}$ таких, что $\xi \prec \eta $ при всех $t \geqslant 0$, выполнено ${{U}_{t}}\xi \prec {{U}_{t}}\eta $;

3) существуют постоянные ${{T}_{2}} > 0$, $\gamma \in (0,1)$ такие, что при любых $\xi ,\eta \in {{\mathfrak{M}}_{0}}$, $\xi \prec \eta $

$l({{U}_{{{{T}_{2}}}}}\eta ) - l({{U}_{{{{T}_{2}}}}}\xi ) \leqslant \gamma (l(\eta ) - l(\xi )).$

Тогда существует единственная стационарная траектория ${{u}_{\infty }}(t) = {{u}_{\infty }} \in {{\mathfrak{M}}_{0}}$, $t \geqslant 0$, а также существуют постоянные $\mu > 0$, ${{T}_{0}} > 0$, $M > 0$ такие, что

$d({{U}_{t}}\xi ,{{u}_{\infty }}) \leqslant M\exp \{ - \mu t\} \quad при\quad t \geqslant {{T}_{0}}.$

5. ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 2

Пусть $\{ {{U}_{t}}\} $, ${{U}_{t}}:{{K}_{\infty }}(\mathcal{D}) \to {{K}_{\infty }}(\mathcal{D})$, ${{U}_{t}}{{\psi }_{0}} = {{\gamma }_{t}}\psi $ – разрешающая полугруппа операторов задачи (2.1)–(2.5). Докажем, что для нее выполнены условия теоремы 7.

Будет использоваться следующее утверждение [13].

Лемма 4. Пусть $E$банахова решетка с нормой $\left\| {\, \cdot \,} \right\|$, монотонной относительно упорядоченности $ \prec $, и конусом ${{E}^{ + }}$ положительных элементов, $\left\| {\left| {\, \cdot \,} \right|} \right\| = \left\| {\, \cdot \,} \right\|$. Пусть $A:E \to E$ – оператор и $\{ {{B}_{\lambda }}\} \in {{\mathfrak{L}}^{ + }}(E)$, $\lambda \in (0,\infty )$, – однопараметрическое семейство операторов со спектральными радиусами $\rho ({{B}_{\lambda }})$, удовлетворяющими условию $\sup \rho ({{B}_{\lambda }}) = \beta < 1$. Тогда если $\left| {A({{y}_{1}}) - A({{y}_{2}})} \right| \prec {{B}_{\lambda }}\left( {\left| {{{y}_{1}} - {{y}_{2}}} \right|} \right)$ для всех ${{y}_{i}} \in E$ таких, что $\left\| {{{y}_{i}}} \right\| \leqslant \lambda $, $i = 1,\;2$, то уравнение

(5.1)
$y = A(y)$
имеет не более одного решения $y \in E$. Если, кроме того, $A:E \to E$ оставляет инвариантным некоторое замкнутое множество $Y \in E$, то решение уравнения (5.1) существует и лежит в $Y$.

Пусть ${{M}_{1}} = h{{\nu }_{{12}}}\kappa {\text{*}}f{\text{*}}{{A}_{{21}}}(4{{\pi }^{{ - 1}}})$, ${{\tau }_{0}} = d{{c}^{{ - 1}}}$. Определим полную подрешетку

${{K}_{0}} = \{ \psi \subset {{K}_{\infty }}(\mathcal{D}):\psi \prec {{M}_{1}}\} .$

Лемма 5. При любых ${{\psi }_{0}} \in {{K}_{\infty }}(\mathcal{D})$ имеет место включение

${{U}_{t}}\psi \in {{K}_{0}},\quad t > {{\tau }_{0}}.$

Доказательство. Пусть $\Phi = \{ \psi ,{{C}_{1}},{{C}_{2}}\} \in {{\mathfrak{M}}_{T}} \cap {{\mathcal{K}}_{T}}$ – решение задачи (2.1)–(2.5). Тогда верно

(5.2)
$\begin{gathered} \psi (x,\nu ,\omega ,t) = {{g}_{0}}(x,\nu ,\omega ,t)\exp \left\{ { - \int\limits_{t_{\omega }^{ - }(x,t)}^t c \mathcal{F}(C)(x + c\omega (s - t),\nu ,s)ds} \right\} + \\ + \;\int\limits_{t_{\omega }^{ - }(x,t)}^t c \mathcal{P}(C)(x + c\omega (\tau - t),\nu ,\tau )\exp \left\{ {\int\limits_t^\tau c \mathcal{F}(C)(x + c\omega (s - t),\nu ,s)ds} \right\}d\tau , \\ \end{gathered} $
при этом имеем

$0 \leqslant \mathcal{F}(C)(x,\nu ,t) \leqslant bpf*,\quad 0 \leqslant \mathcal{P}(C)(x,\nu ,t) \leqslant bf*,\quad x \in G,\quad \nu \in I,\quad t \geqslant 0.$

Так как $\left| {t - t_{\omega }^{ - }(x,t)} \right| \leqslant d{{c}^{{ - 1}}} = {{\tau }_{0}}$, то $t_{\omega }^{ - }(x,t) > 0$ при $t > {{\tau }_{0}}$. Следовательно, справедливо неравенство

$\psi (x,\nu ,\omega ,t) \leqslant {{M}_{1}},\quad t > {{\tau }_{0}},$
что и следовало доказать.

Пусть $\psi _{0}^{1}$, $\psi _{0}^{2} \in {{K}_{0}}$, $\psi _{0}^{1} \succ \psi _{0}^{2}$, ${{\Phi }^{{(i)}}} = \left\{ {{{\psi }^{{(i)}}},C_{1}^{{(i)}},C_{2}^{{(i)}}} \right\} \in {{\mathfrak{M}}_{T}} \cap {{\mathcal{K}}_{T}}$, $i = 1,\;2$ – решения задачи (2.1)–(2.5) c начальными условиями

${{\psi }^{{(i)}}}(x,\nu ,\omega ,0) = \psi _{0}^{i}(x,\nu ,\omega ),\quad i = 1,\;2.$

Положим $\psi = {{\psi }^{{(1)}}} - {{\psi }^{{(2)}}}$, $C = C_{1}^{{(1)}} - C_{1}^{{(2)}}$, ${{\psi }_{0}} = \psi _{0}^{1} - \psi _{0}^{2}$. Тогда справедливы равенства

(5.3)
$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\psi (x,\nu ,\omega ,t) + \mathcal{F}({{C}^{{(1)}}})(x,\nu ,t)\psi (x,\nu ,\omega ,t) = K(x,\nu ,\omega ,t)C(x,t),$
(5.4)
$C(x,t) = H(x,t)\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \psi (x,\nu ,\omega ,t)d\omega d\nu ,$
где

$K(x,\nu ,\omega ,t) = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}\left( {{{A}_{{21}}} + ({{B}_{{12}}} + {{B}_{{21}}}){{\psi }^{{(2)}}}(x,\nu ,\omega ,t)} \right),$
$H(x,t) = f(x)\frac{{{{A}_{{21}}}{{B}_{{21}}} + ({{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}}){{n}_{e}}(x)}}{{\mathcal{R}({{C}^{{(1)}}})(x,t)\mathcal{R}({{C}^{{(2)}}})(x,t)}},$
$0 \leqslant K \leqslant q = h{{\nu }_{{12}}}\frac{{\kappa (\nu )}}{{4\pi }}({{A}_{{21}}} + ({{B}_{{12}}} + {{B}_{{21}}}){{M}_{1}}),\quad 0 \leqslant \mathcal{F}({{C}^{{(1)}}}) \leqslant bp,\quad 0 \leqslant H \leqslant f{\text{*}}p.$

Лемма 6. При любых $\psi _{0}^{1}$, $\psi _{0}^{2} \in {{K}_{0}}$ таких, что $\psi _{0}^{1} \succ \psi _{0}^{2}$, выполнено соотношение

(5.5)
${{U}_{t}}\psi _{0}^{1} \succ {{U}_{t}}\psi _{0}^{2},\quad t \geqslant 0.$

Доказательство. Пусть $\xi = {{e}^{{ - \alpha t}}}\psi $, $\eta = {{e}^{{ - \alpha t}}}C$, где $\alpha = cqpf{\text{*}}$. Тогда

(5.6)
$\frac{1}{c}{{\left( {\frac{d}{{d\tau }}} \right)}_{\omega }}\xi (x,\nu ,\omega ,t) + (\alpha + \mathcal{F}({{C}^{{(1)}}})(x,\nu ,t))\xi (x,\nu ,\omega ,t) = G(x,\nu ,\omega ,t)\eta (x,t),$
(5.7)
$\eta (x,t) = H(x,t)\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \xi (x,\nu ,\omega ,t)d\omega d\nu .$
Уравнение (5.6) эквивалентно уравнению
(5.8)
$\begin{gathered} \xi (x,\nu ,\omega ,t) = {{g}_{0}}(x,\nu ,\omega ,t)\exp \left\{ { - \int\limits_{t_{\omega }^{ - }(x,t)}^t {\left( {\alpha + c\mathcal{F}({{C}^{{(1)}}})(x + c\omega (s - t),\nu ,s)} \right)} ds} \right\} + \\ + \;\int\limits_{t_{\omega }^{ - }(x,t)}^t c K(x + c\omega (\tau - t),\nu ,\omega ,\tau )\eta (x + c\omega (\tau - t),\tau ) \times \\ \times \;\exp \left\{ {\int\limits_t^\tau {\left( {\alpha + c\mathcal{F}({{C}^{{(1)}}})(x + c\omega (s - t),\nu ,s)} \right)} ds} \right\}d\tau . \\ \end{gathered} $
Пусть $\Phi = \{ \xi ,\eta \} $. Систему (5.7), (5.8) можно представить в виде
(5.9)
$A\Phi = \Phi ,$
где оператор $A:{{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T]) \to {{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])$ ставит в соответствие вектору $\Phi $ вектор из правых частей уравнений (5.8) и (5.7).

Пространство ${{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])$ является банаховой решеткой относительно порядка, порожденного конусом ${{K}_{\infty }}(\mathcal{D} \times [0,T]) \times {{K}_{\infty }}(G \times [0,T])$.

Положим

$B(\Phi ) = \left\{ {\int\limits_{t_{\omega }^{ - }(x,t)}^t \eta (x + c\omega {\kern 1pt} (\tau - t),\tau )d\tau ,\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } \xi (x,\nu ,\omega ,t){\kern 1pt} d\omega d\nu } \right\}.$
Оператор $B:{{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T]) \to {{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])$ оставляет инвариантным конус ${{K}_{\infty }}(\mathcal{D} \times [0,T]) \times {{K}_{\infty }}(G \times [0,T])$ и справедливы оценки
$\left\| {{{B}^{{2k - 1}}}} \right\| \leqslant \frac{{\tau _{0}^{{k - 1}}}}{{(k - 1)!}}\max \left\{ {1,\frac{{{{\tau }_{0}}}}{k}} \right\},\quad \left\| {{{B}^{{2k}}}} \right\| \leqslant \frac{{\tau _{0}^{k}}}{{k!}}\quad k = 1,\;2,\; \ldots $
Поэтому спектральный радиус $\rho (B) = \mathop {\lim }\nolimits_{n \to \infty } \sqrt[n]{{\left\| {{{B}^{n}}} \right\|}}$ равен нулю.

Для любых ${{\Phi }_{1}}$, ${{\Phi }_{2}} \in {{L}_{\infty }}(\mathcal{D} \times [0,T]) \times {{L}_{\infty }}(G \times [0,T])$ выполняется неравенство

$\left| {A({{\Phi }_{1}}) - A({{\Phi }_{2}})} \right| \prec \beta B(\left| {{{\Phi }_{1}} - {{\Phi }_{2}}} \right|),\quad \beta = \max \{ q,f{\text{*}}p\} .$
Кроме того, оператор $A$ оставляет инвариантным замкнутое множество

$S = \{ \{ \xi ,\eta \} \in {{K}_{\infty }}(\mathcal{D} \times [0,T]) \times {{K}_{\infty }}(G \times [0,T]):{{\left\| \xi \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}} \leqslant {{\left\| {{{\psi }_{0}}} \right\|}_{{{{L}_{\infty }}(\mathcal{D} \times [0,T])}}}\} .$

Таким образом, для операторов $A$ и ${{B}_{\lambda }} = \beta B$, $\lambda \geqslant 0$ выполнены условия леммы 4 и $\psi \geqslant 0$, $C \geqslant 0$, т.е. ${{\psi }^{{(1)}}} \succ {{\psi }^{{(2)}}}$.

Лемма 7. Пусть $\psi _{0}^{i} \in {{K}_{0}}$, $i = 1,\;2$, $\psi _{0}^{1} \succ \psi _{0}^{2}$, $\psi _{0}^{1} \ne \psi _{0}^{2}$. Тогда

(5.10)
$\rho ({{U}_{{2{{\tau }_{0}}}}}\psi _{0}^{1},{{U}_{{2{{\tau }_{0}}}}}\psi _{0}^{2}) \leqslant \gamma \rho (\psi _{0}^{1},\psi _{0}^{2}),$
где $\gamma \in (0,\;1)$некоторая постоянная, не зависящая от $\psi _{0}^{1}$, $\psi _{0}^{2}$.

Доказательство. Из уравнения (5.3) вытекает справедливость неравенства [13]

(5.11)
$\begin{gathered} \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} {{{{\left( {\frac{d}{{dt}}} \right)}}_{\omega }}} } \psi (x,\nu ,\omega ,t)dxd\nu d\omega dt \geqslant \int\limits_\mathcal{D} {(\psi (x,\nu ,\omega ,2{{\tau }_{0}})} - \psi (x,\nu ,\omega ,0))dxd\nu d\omega + \\ + \;\varepsilon \int\limits_\mathcal{D} {{{\psi }_{0}}} (x,\nu ,\omega )dxd\nu d\omega + c\varepsilon \int\limits_0^{2{{\tau }_{0}}} K (x,\nu ,\omega ,t)C(x,t)dxd\nu d\omega dt, \\ \end{gathered} $
где $\varepsilon = \exp ( - \kappa {\text{*}}f{\text{*}}{{B}_{{12}}}h{{\nu }_{{12}}}d{{(4\pi )}^{{ - 1}}})$. С другой стороны,
$\begin{gathered} \frac{1}{c}\int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} {{{{\left( {\frac{d}{{dt}}} \right)}}_{\omega }}} } \psi (x,\nu ,\omega ,t)dxd\nu d\omega dt = \\ = \;\int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} K } (x,\nu ,\omega ,t)C(x,t)dxd\nu d\omega dt - \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} \mathcal{F} } ({{C}^{{(1)}}})(x,\nu ,t)\psi (x,\nu ,\omega ,t)dxd\nu d\omega dt. \\ \end{gathered} $
Следовательно,
$\begin{gathered} \frac{1}{c}\int\limits_\mathcal{D} {(\psi (x,\nu ,\omega ,2{{\tau }_{0}})} dxd\nu d\omega \leqslant \frac{{1 - \varepsilon }}{c}\int\limits_\mathcal{D} {{{\psi }_{0}}} (x,\nu ,\omega )dxd\nu d\omega + \\ + \;\int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} K } (x,\nu ,\omega ,t)C(x,t)dxd\nu d\omega dt - \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} \,\mathcal{F}} ({{C}^{{(1)}}})(x,\nu ,t)\psi (x,\nu ,\omega ,t)dxd\nu d\omega dt. \\ \end{gathered} $
Используя явный вид $\mathcal{F}({{C}^{{(1)}}})$, $K$ и $C$, получаем
$\begin{gathered} \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} K } (x,\nu ,\omega ,t)C(x,t)dxd\nu d\omega dt = h{{\nu }_{{12}}}\int\limits_0^{2{{\tau }_{0}}} {\int\limits_G {\frac{{g(x)f(x)}}{{{{\mathcal{R}}_{1}}(x,t){{\mathcal{R}}_{2}}(x,t)}}} } J(\psi )(x,t) \times \\ \times \;\int\limits_I {\int\limits_\Omega {\frac{{\kappa (\nu )}}{{4\pi }}} } ({{A}_{{21}}} + ({{B}_{{12}}} + {{B}_{{21}}}){{\psi }^{{(2)}}}(x,\nu ,\omega ,t))d\omega d\nu dxdt = \\ = \;h{{\nu }_{{12}}}\int\limits_0^{2{{\tau }_{0}}} \int\limits_G \frac{{g(x)f(x)}}{{{{\mathcal{R}}_{1}}(x,t){{\mathcal{R}}_{2}}(x,t)}}J(\psi )(x,t)({{A}_{{21}}} + ({{B}_{{12}}} + {{B}_{{21}}})J({{\psi }^{{(2)}}})(x,t))dxdt, \\ \end{gathered} $
$g(x) = {{A}_{{21}}}{{B}_{{12}}} + ({{B}_{{12}}}{{C}_{{21}}} - {{B}_{{21}}}{{C}_{{12}}}){{n}_{e}}(x),$
$\int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} \mathcal{F} } ({{C}^{{(1)}}})(x,\nu ,t)\psi (x,\nu ,\omega ,t)dxd\nu d\omega dt = h{{\nu }_{{12}}}\int\limits_0^{2{{\tau }_{0}}} {\int\limits_G {\frac{{g(x)f(x)}}{{{{\mathcal{R}}_{1}}(x,t)}}} } J(\psi )(x,t)dxdt.$
Таким образом,
$\begin{gathered} \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} K } (x,\nu ,\omega ,t)C(x,t)dxd\nu d\omega dt - \int\limits_0^{2{{\tau }_{0}}} {\int\limits_\mathcal{D} \mathcal{F} } ({{C}^{{(1)}}})(x,\nu ,t)\psi (x,\nu ,\omega ,t)dxd\nu d\omega dt = \\ = \;h{{\nu }_{{12}}}\int\limits_0^{2{{\tau }_{0}}} {\int\limits_G {\frac{{g(x)f(x)J(\psi )(x,t)}}{{{{\mathcal{R}}_{1}}(x,t){{\mathcal{R}}_{2}}(x,t)}}} } ( - ({{B}_{{12}}} + {{B}_{{21}}})J(\psi )(x,t) - ({{C}_{{12}}} + {{C}_{{21}}}){{n}_{e}}(x))dxdt \leqslant 0, \\ \frac{1}{c}\int\limits_\mathcal{D} {\psi (x,\nu ,\omega ,2{{\tau }_{0}})} dxd\nu d\omega \leqslant \frac{{1 - \varepsilon }}{c}\int\limits_\mathcal{D} {{{\psi }_{0}}} (x,\nu ,\omega )dxd\nu d\omega . \\ \end{gathered} $
Последнее неравенство означает, что справедливо неравенство (5.10), где $\gamma = 1 - \varepsilon $.

Из доказанных лемм 5–7 при ${{T}_{0}} = 2{{\tau }_{0}}$ на основании теоремы 7 следует теорема 2.

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

  1. Соболев В.В. Перенос лучистой энергии в атмосферах звезд и планет. М.: Гостехтеоритиздат, 1956.

  2. Белл Д., Глесстон С. Теория ядерных реакторов. М.: Атомиздат, 1974.

  3. Михалас Д. Звездные атмосферы. М.: Мир, 1982.

  4. Иванов В.В. Перенос излучения и спектры небесных тел. М.: Наука, 1969.

  5. Гермогенова Т.А. Локальные свойства решений уравнения переноса. М.: Наука, 1986.

  6. Владимиров В.С. Математические задачи односкоростной теории переноса частиц // Тр. Матем. Ин-та им. В.А. Стеклова АН СССР. 1961. Вып. 61. С. 2–158.

  7. Сушкевич Т.А. Математические модели переноса излучения. М.: БИНОМ, 2006.

  8. Калинин А.В., Морозов С.Ф. Об одной нелинейной краевой задаче теории переноса излучения // Ж. вычисл. матем. и матем. физ. 1990. Т. 30. С. 1071–1080.

  9. Морозов С.Ф., Сумин В.И. Нелинейные интегродифференциальные системы уравнений нестационарного переноса // Сиб. матем. журнал. 1978. Т. 19:4. С. 842–848.

  10. Калинин А.В., Морозов С.Ф. О разрешимости “в целом” нелинейной задачи переноса излучения // Дифференц. ур-ния. 1985. Т. 21. № 3. С. 482–494.

  11. Калинин А.В., Морозов С.Ф. О стабилизации решения нелинейной системы переноса излучения в двухуровневом приближении // Докл. АН СССР. 1990. Т. 311. № 2. С. 343–346.

  12. Калинин А.В., Морозов С.Ф. Задача Коши для одного нелинейного интегродифференциального уравнения переноса // Матем. заметки. 1997. Т. 61. Вып. 5. С. 677–686.

  13. Калинин А.В., Морозов С.Ф. Смешанная задача для нестационарной системы нелинейных интегродифференциальных уравнений // Сиб. матем. журнал. 1999. Т. 40. № 5. С. 1052–1066.

  14. Биркгоф Г. Теория решеток. М.: Мир, 1984.

  15. Канторович Л.В., Акилов Г.П. Функциональный анализ. М.: Наука, 1977.

  16. Красносельский М.А. Положительные решения операторных уравнений. М.: Физматгиз, 1962.

  17. Tarski A.A. A lattice-theoretical fixpoint theorem and its applications // Pacif. J. Math. 1955. V. 5. № 2. P. 285–309.

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