Журнал вычислительной математики и математической физики, 2019, T. 59, № 4, стр. 621-648

Разрушение решений неклассических нелокальных нелинейных модельных уравнений

М. О. Корпусов 12*

1 МГУ физ. ф-т
119992 Москва, Ленинские горы, Россия

2 РУДН
117198 Москва, ул. Миклухо-Маклая, 6, Россия

* E-mail: korpusov@gmail.com

Поступила в редакцию 01.12.2016
После доработки 16.06.2017
Принята к публикации 14.11.2018

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

Аннотация

Рассматривается абстрактная задача Коши для некоторого нелинейного нелокального дифференциально-операторного уравнения первого порядка, причем эта задача Коши является обобщением ряда модельных физических примеров. Для абстрактной задачи Коши получены результаты о существовании непродолжаемого во времени классического решения, при некоторых достаточных условиях получены результаты о разрушении за конечное время, для которого получены двусторонние оценки. Наконец, при некоторых условиях доказана глобальная корректность задачи вне зависимости от величины начальной функции. Библ. 29.

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

1. ВВЕДЕНИЕ

В своих двух классических работах [1] и [2] в 1973, 1974 г. предложен новый энергетический метод исследования возникновения разрушения в двух задачах Коши для абстрактных формально параболического и формально гиперболического уравнений:

$P{{u}_{t}} = - Au + F(u),\quad u(0) = {{u}_{0}},$
$P{{u}_{{tt}}} = - Au + F(u),\quad u(0) = {{u}_{0}},\quad u{\kern 1pt} '(0) = {{u}_{1}},$
где операторы $P$ и $A$ линейные, а оператор $F(u)$ нелинейный. Результаты о разрушении были получены как для классических решений, так и для слабых решений. Отметим также работу [3].

В 1977 г. опубликована работа [4], в которой энергетический метод был использован для решения такой пары абстрактных задач Коши:

$P{{u}_{t}} = - Au + B(u) + F(t,u),\quad u(0) = {{u}_{0}},$
$P{{u}_{{tt}}} = - Au + B(u) - aP{{u}_{t}} + F(t,u),\quad u(0) = {{u}_{0}},\quad u{\kern 1pt} '(0) = {{u}_{1}},$
где оператор $B(u)$ может быть нелинейным. В 1997 г. вышла работа [5], в которой уже рассматривалось следующее уравнение:
${{(P({{u}_{t}}))}_{t}} + A(u) + Q(t,{{u}_{t}}) = F(u),\quad u(0) = {{u}_{0}},\quad u{\kern 1pt} '(0) = {{u}_{1}}.$
В 1998 г. в работе [6] было рассмотрено следующее общее уравнение формально параболического типа:
$Q(t,u,{{u}_{t}}) + A(t,u) = F(t,u),\quad u(0) = {{u}_{0}}.$
В это же время начинает развиваться тематика доказательства разрушения решений формально гиперболических уравнений с положительной энергией. В этой связи отметим следующие работы P. Pucci и J. Serrin [7] и [8]. Отметим также недавние результаты в этом направлении [9]–[12].

Отметим, что важным развитием энергетического метода явилась работа [13], в которой рассматривалась первая начально-краевая задача для уравнения

${{u}_{{tt}}} - \Delta u + a{{u}_{t}}{\text{|}}{{u}_{t}}{{{\text{|}}}^{{m - 1}}} = bu{\text{|}}u{{{\text{|}}}^{{p - 1}}}$
в цилиндре $[0,T] \times \Omega $. Результат был обобщен в работах [14] и [15] на случай следующего нелинейного нелокального уравнения:
${{u}_{{tt}}} - \Delta u + \int\limits_0^t \,g(t - s)\Delta u(x,s)ds\; + \;{\text{|}}{{u}_{t}}{{{\text{|}}}^{{m - 2}}}{{u}_{t}} = {\text{|}}u{{{\text{|}}}^{{p - 2}}}u.$
Наша работа продолжает исследования, начатые в работах [16]–[18]. В работе [18] рассматривалась абстрактная задача Коши следующего вида:
(1.1)
$A\frac{{{{d}^{2}}u}}{{d{{t}^{2}}}} + \frac{d}{{dt}}\left( {{{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u)} \right) + H_{f}^{'}(u) = F_{f}^{'}(u),\quad u(0) = {{u}_{0}},\quad u{\kern 1pt} '(0) = {{u}_{1}},$
где $H_{f}^{'}(u),$ $F_{f}^{'}(u)$ – это производные Фреше нелинейных операторов. В этой работе мы рассмотрели классическое и слабое решения задачи Коши (1.1) и доказали результаты о существовании непродолжаемых решений, а также получили достаточные условия разрушения решений за конечное время. Для доказательства разрушения за конечное время мы пользовали нашу модификацию энергетического метода, изложенную в работе [19].

В настоящей работе мы рассмотрим абстрактную задачу Коши следующего вида:

(1.2)
$\frac{d}{{dt}}\left( {{{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u)} \right) + {{L}_{1}}u + \int\limits_0^t \,dsh(t - s){{L}_{2}}u(s) + DP(u) = F(u),\quad u(0) = {{u}_{0}},$
где операторы ${{A}_{0}}$ и ${{L}_{1}},$ ${{L}_{2}},$ $D$ линейные, а операторы ${{A}_{j}}(u)$ и $F(u),$ $P(u)$ нелинейные. Для доказательства существования сильного решения этой задачи Коши мы применим метод монотонных операторов Браудера–Минти [20], а для доказательства разрушения за конечное время применим метод энергетических оценк, развитый в работе [19].

2. СИСТЕМЫ УРАВНЕНИЙ КВАЗИСТАЦИОНАРНОГО ПОЛЯ В КРИСТАЛЛИЧЕСКИХ ПОЛУПРОВОДНИКАХ

В этой работе мы рассмотрим модельные задачи о возникновении “пробоя” в кристаллических полупроводниках при наличии сильной временной дисперсии.

Начнем с рассмотрения системы уравнений квазистационарного электрического поля. Общая система уравнений в этом случае имеет следующий вид (см., например, [21]):

(2.1)
$\operatorname{div} {\mathbf{D}} = - 4\pi n,\quad \operatorname{rot} {\mathbf{E}} = 0,$
(2.2)
$\frac{{\partial n}}{{\partial t}} = \operatorname{div} {\mathbf{J}} + Q,$
где ${\mathbf{D}}$ – вектор индукции электрического поля, ${\mathbf{E}}$ – вектор напряженности электрического поля, ${\mathbf{J}}$ – вектор плотности тока свободных зарядов, $n$ – концентрация свободных зарядов, наконец, слагаемое $Q$ описывает распределение источников или стоков свободных зарядов.

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

В полупроводниках распределение связанных зарядов определяется вектором поляризуемости среды ${\mathbf{P}}$, который связан с вектором индукции электрического поля ${\mathbf{D}}$ равенством

(2.3)
${\mathbf{D}} = {\mathbf{E}} + 4\pi {\mathbf{P}}.$
При этом распределение плотности связанных зарядов в полупроводнике определяется величиной
(2.4)
$\rho = 4\pi \operatorname{div} {\mathbf{P}}.$
В дальнейшем мы будем рассматривать различные модели кристаллических полупроводников. В одних моделях вектор поляризуемости ${\mathbf{P}}$ связан с вектором напряженности электрического поля ${\mathbf{E}}$, например, таким соотношением
(2.5)
${\mathbf{P}} = {{\varkappa }_{0}}\mathop {\left| {\mathbf{E}} \right|}\nolimits^2 {\mathbf{E}},\quad {{\varkappa }_{0}} > 0,$
а в других моделях мы будем рассматривать непосредственно распределение связанных зарядов $\rho $ в виде некоторой функции от потенциала $\phi $ электрического поля ${\mathbf{E}}$ (который, заметим, существует в предположении поверхностной односвязанности области $\Omega \subset {{\mathbb{R}}^{3}},$ где рассматривается система уравнений (2.1)), например, такого вида
$\rho = \operatorname{div} {\mathbf{P}} = {\text{|}}\phi {{{\text{|}}}^{q}}\phi ,\quad q \geqslant 0.$
Теперь рассмотрим вектор плотности тока свободных зарядов ${\mathbf{J}}$. Во многих моделях его можно рассматривать зависящим от напряженности электрического поля следующим образом:
${\mathbf{J}} = \sigma ({\text{|}}{\mathbf{E}}{\text{|}}){\mathbf{E}},$
где величина $\sigma (r)$ по своему физическому смыслу есть коэффициент проводимости среды, который зависит, вообще говоря, от модуля напряженности электрического поля ${\mathbf{E}}$. Кроме того, в данной работе мы будем рассматривать случай сильной временной дисперсии, т.е. когда зависимость плотности тока свободных зарядов от напряженности поля является нелокальной. Типичный пример такой зависимости следующий:
(2.6)
${\mathbf{J}}(x,t) = \int\limits_0^t \,ds\sigma (t - s){\mathbf{E}}(x,s),\quad \sigma (t) \in \mathbb{C}[0, + \infty ).$
Наконец, обсудим слагаемое $Q$ в правой части уравнения (2.2). Как мы уже говорили, это слагаемое описывает распределение источников или стоков свободных зарядов в полупроводнике и представляет собой функцию от потенциала $\phi $ самосогласованного электрического поля ${\mathbf{E}}$. В настоящей монографии мы рассматриваем тот важный случай, когда функция $Q$ описывает распределение источников свободных зарядов, которое мы будем моделировать следующим образом:
(2.7)
$Q(\phi ) = {\text{|}}\phi {{{\text{|}}}^{q}}\phi ,\quad q > 0.$
Заметим, что мы будем рассматривать два механизма возникновения пробоя в кристаллических полупроводниках – это пробой, вызванный источниками свободных зарядов, распределение которых описывается функцией (2.7), или пробой, вызванный наличием отрицательной дифференциальной проводимости среды (см., например, [22]). Отрицательность дифференциальной проводимости среды учитывается тем, что вектор плотности тока электрического поля ${\mathbf{J}}$ связан с напряженностью поля ${\mathbf{E}}$ соотношением вида
${\mathbf{J}} = - \sigma ({\text{|}}{\mathbf{E}}{\text{|}}){\mathbf{E}},\quad \sigma (r) > 0\quad {\text{п р и }}\quad r > 0.$
Итак, мы будем рассматривать эти два механизма возникновения пробоя в кристаллических полупроводниках.

Теперь мы обсудим формулу связи индукции электрического поля ${\mathbf{D}}$ от напряженности электрического поля ${\mathbf{E}}$. Для этого с учетом формулы (2.3) надо обсудить различные механизмы зависимости вектора ${\mathbf{P}}$ от ${\mathbf{E}}$. Действительно, помимо формулы (2.5), которая учитывает керровскую зависимость вектора ${\mathbf{P}}$ от вектора ${\mathbf{E}}$ можно учесть еще так называемую пространственную дисперсию среды. Давайте обсудим, как учитывается пространственная дисперсия среды. Как известно (см., например, [23]), в этом случае связь ${\mathbf{P}}$ и ${\mathbf{E}}$ является операторной:

(2.8)
${\mathbf{P}} = \hat {\varkappa }{\mathbf{E}},$
где оператор $\hat {\varkappa }$ поляризуемости среды имеет, вообще говоря, следующий нелокальный характер:
(2.9)
$\hat {\varkappa }f(x) = \int\limits_{{{R}^{3}}} dy\varkappa (x - y)f(y).$
Заметим, что формула (2.9) записана для всего трехмерного эвклидова пространства ${{\mathbb{R}}^{3}}$. Предположим, что формула (2.8) рассматривается так же во всем трехмерном пространстве ${{\mathbb{R}}^{3}}$. Тогда формально применим преобразование Фурье к обеим частям равенства (2.8) и получим следующее алгебраическое равенство:
(2.10)
${\mathbf{\hat {P}}}(k) = \varkappa (k){\mathbf{\hat {E}}}(k),$
где $k = ({{k}_{1}},{{k}_{2}},{{k}_{3}}) \in {{\mathbb{R}}^{3}}$ – это волновой вектор, соответствующий рассматриваемому преобразованию Фурье. Теперь предположим, что
$\varkappa (k) \in {{\mathbb{C}}^{\infty }}({{R}^{3}}).$
Тогда эту функцию можно разложить в ряд Тейлора в окрестности точки $(0,0,0) \in {{R}^{3}}$ с остаточным членом в форме Пеано:
(2.11)
$\varkappa (k) = \sum\limits_{{{\alpha }_{1}},{{\alpha }_{2}},{{\alpha }_{3}} = 1,1,1}^{{{n}_{1}},{{n}_{2}},{{n}_{3}}} {{c}_{{{{\alpha }_{1}}{{\alpha }_{2}}{{\alpha }_{3}}}}}k_{1}^{{{{\alpha }_{1}}}}k_{2}^{{{{\alpha }_{2}}}}k_{3}^{{{{\alpha }_{3}}}} + \bar {o}(k_{1}^{{{{n}_{1}}}}k_{2}^{{{{n}_{2}}}}k_{3}^{{{{n}_{3}}}}).$
В этом месте обычно многие физики-теоретики рассматривают вместо формулы (2.11) ее “укороченный” вариант. А именно следующую формулу:
(2.12)
${{\varkappa }_{0}}(k) = \sum\limits_{{{\alpha }_{1}},{{\alpha }_{2}},{{\alpha }_{3}} = 1,1,1}^{{{n}_{1}},{{n}_{2}},{{n}_{3}}} {{c}_{{{{\alpha }_{1}}{{\alpha }_{2}}{{\alpha }_{3}}}}}k_{1}^{{{{\alpha }_{1}}}}k_{2}^{{{{\alpha }_{2}}}}k_{3}^{{{{\alpha }_{3}}}}.$
При этом формула (2.10) примет следующий вид:
(2.13)
${\mathbf{\hat {P}}}(k) = {{\varkappa }_{0}}(k){\mathbf{\hat {E}}}(k).$
Если теперь формально применить обратное преобразование Фурье к равенству (2.13), то с учетом формулы (2.12) мы получим выражение
(2.14)
${\mathbf{P}}(x) = \sum\limits_{{{\alpha }_{1}},{{\alpha }_{2}},{{\alpha }_{3}} = 1,1,1}^{{{n}_{1}},{{n}_{2}},{{n}_{3}}} \,{{c}_{{{{\alpha }_{1}}{{\alpha }_{2}}{{\alpha }_{3}}}}}\mathop {\left( {i\frac{\partial }{{\partial {{x}_{1}}}}} \right)}\nolimits^{{{\alpha }_{1}}} \mathop {\left( {i\frac{\partial }{{\partial {{x}_{2}}}}} \right)}\nolimits^{{{\alpha }_{2}}} \mathop {\left( {i\frac{\partial }{{\partial {{x}_{3}}}}} \right)}\nolimits^{{{\alpha }_{3}}} {\mathbf{E}}(x).$
Таким образом, мы пришли к искомой модельной связи векторов ${\mathbf{P}}$ и ${\mathbf{E}}$. У многих физиков-теоретиков укороченная функция ${{\varkappa }_{0}}(k)$ имеет следующий простой вид:
${{\varkappa }_{0}}(k) = {{\varkappa }_{0}}{\text{|}}k{{{\text{|}}}^{2}},\quad {{\varkappa }_{0}} > 0,\quad {\text{|}}k{{{\text{|}}}^{2}} = k_{1}^{2} + k_{2}^{2} + k_{3}^{2}.$
И тогда выражение (2.14) примет следующий вид:
(2.15)
${\mathbf{P}}(x) = - {{\varkappa }_{0}}{{\Delta }_{x}}{\mathbf{E}}(x),\quad {{\Delta }_{x}} \equiv \frac{{{{\partial }^{2}}}}{{\partial x_{1}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{2}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{3}^{2}}}.$
Заметим, что эти рассуждения, приведшие нас к формуле (2.14), имеют “физический” характер. Кроме того, формула (2.14) выведена только для всего пространства $x = ({{x}_{1}},{{x}_{2}},{{x}_{3}}) \in {{\mathbb{R}}^{3}}$. Однако формально формула (2.14) корректна и для произвольной области $\Omega \subset {{\mathbb{R}}^{3}}$. Поэтому мы будем говорить о пространственной дисперсии и в случае произвольной области изменения переменных $x = ({{x}_{1}},{{x}_{2}},{{x}_{3}}) \in \Omega \subset {{\mathbb{R}}^{3}}$, имея в виду соотношение вида (2.14).

Теперь обсудим еще один фактор, который влияет на общую ситуацию в кристаллических полупроводниках – это наличие внешнего электрического поля. Пусть ${{{\mathbf{E}}}_{0}}$ – это внешнее постоянное электрическое поле. Тогда согласно общей макроскопической теории электричества (см., например, [24]) наличие этого внешнего поля приводит к возникновению тока свободных зарядов ${{{\mathbf{J}}}_{0}},$ направление которого противоположно направлению внешнего поля ${{{\mathbf{E}}}_{0}}{\text{:}}$

${{{\mathbf{J}}}_{0}} = e{{n}_{0}}(\phi ){{{\mathbf{E}}}_{0}},$
где ${{n}_{0}}(\varphi )$ – есть “квазистационарное” распределение плотности свободных зарядов. Мы будем моделировать эту функцию следующим образом:
(2.16)
${{n}_{0}} = {{n}_{0}}(\phi ) = 1 + {{a}_{1}}\phi + {{a}_{2}}{{\phi }^{2}},\quad {{a}_{1}},{{a}_{2}} \in {{\mathbb{R}}^{1}}.$
Наличие этого постоянного внешнего электрического поля ${{{\mathbf{E}}}_{0}}$ приводит к появлению дополнительного слагаемого в уравнении (2.2). Именно, уравнение (2.2) примет следующий вид:
$\frac{{\partial n}}{{\partial t}} = {\text{div}}({\mathbf{J}} + {{{\mathbf{J}}}_{0}}) + Q.$
Таким образом, мы учли самые разнообразные факторы, имеющие место в кристаллических полупроводниках. Теперь давайте обсудим, каким образом мы будем учитывать все эти факторы. Поскольку все эти факторы имеют отношения к векторам ${\mathbf{J}}$ и ${\mathbf{P}},$ то и учитывать все эти факторы мы будем в виде линейной суперпозиции:
(2.17)
${\mathbf{J}} = {{{\mathbf{J}}}_{1}} + \ldots + {{{\mathbf{J}}}_{n}},\quad {\mathbf{P}} = {{{\mathbf{P}}}_{1}} + \ldots + {{{\mathbf{P}}}_{n}}.$
Действительно, каждое слагаемое в (2.17) отвечает за свой эффект, имеющий место в кристаллическом полупроводнике.

Приступим к рассмотрению системы уравнений квазистационарного магнитного поля в кристаллическом полупроводнике. Общая система уравнений в отсутствие тока свободных зарядов имеет следующий вид (см., например, [23]):

$\operatorname{div} {\mathbf{B}} = 0,\quad \operatorname{rot} {\mathbf{H}} = 0,\quad {\mathbf{B}} = {\mathbf{H}} + 4\pi {\mathbf{M}},$
где ${\mathbf{B}},$ ${\mathbf{H}}$ и ${\mathbf{M}}$ – это векторы индукции, напряженности и намагниченности магнитного поля. Данную систему уравнений необходимо дополнить уравнением, связывающим векторы ${\mathbf{M}}$ и ${\mathbf{H}}$. Действительно, в качестве такового уравнения возьмем уравнение Ландау–Лифшица (см. [24]):
(2.18)
$\frac{{\partial {\mathbf{M}}}}{{\partial t}} = \gamma \left[ {{\mathbf{H}},{\mathbf{M}}} \right] + {\mathbf{R}},$
где вектор ${\mathbf{R}}$ характеризует, в частности, распределение “источников” магнитных доменов в кристаллическом полупроводнике–магнетике. Довольно часто уравнение (2.18) можно существенно упростить, воспользовавшись тем, что вектор намагниченности представим в виде следующей суммы:
(2.19)
${\mathbf{M}} = \mathfrak{M} + {\mathbf{m}},\quad {\text{|}}{\mathbf{m}}{\text{|}} \ll \left| \mathfrak{M} \right|,$
т.е. в виде “квазистационарной” намагниченности $\mathfrak{M}$ и малой, быстропеременной добавки ${\mathbf{m}}$. С учетом (2.19) из (2.18) получим уравнение
$\frac{{\partial {\mathbf{M}}}}{{\partial t}} = \gamma \left[ {{\mathbf{H}},\mathfrak{M}} \right] + {\mathbf{R}}.$
Теперь мы учтем некоторые факторы, имеющие место в кристаллических полупроводниках магнетиках.

Прежде всего заметим, что можно учесть, что квазистационарная намагниченность состоит из нескольких слагаемых, каждое из которых отвечает за учет некоторого фактора:

$\mathfrak{M} = {{\mathfrak{M}}_{1}} + \cdot \cdot \cdot + {{\mathfrak{M}}_{n}}.$
Во-первых, как правило, $\mathfrak{M}$ содержит линейную по полю часть вида
${{\mathfrak{M}}_{{1j}}} = {{a}_{j}}{{H}_{j}},\quad {{a}_{j}} \in {{R}^{1}},\quad j = 1,2,3.$
Во-вторых, точно также как и ранее, можно учесть сильную пространственную дисперсию кристаллического полупроводника–магнетика, т.е. следующую связь векторов ${{\mathfrak{M}}_{2}}$ и H:
(2.20)
${{\mathfrak{M}}_{2}} = - {{\chi }_{0}}\Delta {\mathbf{H}},\quad \Delta \equiv \frac{{{{\partial }^{2}}}}{{\partial x_{1}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{2}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{3}^{2}}}.$
Вывод уравнения (2.20) точно такой же, как и вывод уравнения (2.15). Вектор ${\mathbf{R}}$ также может быть записан в виде следующей суперпозиции:
(2.21)
${\mathbf{R}} = {{{\mathbf{R}}}_{1}} + \ldots + {{{\mathbf{R}}}_{n}}.$
Можно учесть наличие “источников” доменов, описываемых формулой вида
${{{\mathbf{R}}}_{1}} = - {{a}_{3}}{\text{|}}{\mathbf{H}}{{{\text{|}}}^{{p - 2}}}{\mathbf{H}},\quad {{a}_{3}} > 0.$
Наконец, учтем наличие релаксационных механизмов или, иначе говоря, сильную временную дисперсию следующего вида:

(2.22)
${{{\mathbf{R}}}_{2}}(t) = {{a}_{4}}\int\limits_0^t \,dsh(t - s){\mathbf{H}}(s),\quad {{a}_{4}} \in {{\mathbb{R}}^{1}},\quad h(t) \in \mathbb{C}[0, + \infty ).$

Теперь мы обсудим граничные условия для рассматриваемых систем векторных уравнений. Конечно, их вид и количество для корректной постановки рассматриваемых далее задач зависят от того факта, учитываем ли мы сильную пространственную дисперсию или нет. В основном мы будем рассматривать тот случай, когда граница материальной среды представляет собой заземленную и идеально проводящую среду. Тогда на границе $\partial \Omega $ области $\Omega $ имеют место следующие равенства:

(2.23)
${{\left. \phi \right|}_{{\partial \Omega }}} = {{\left. {\frac{{\partial \phi }}{{\partial {{n}_{x}}}}} \right|}_{{\partial \Omega }}} = 0,$
где $\phi $ – это потенциал электрического поля. При рассмотрении кристаллического полупроводника–магнетика мы будем рассматривать аналогичные граничные условия вида
(2.24)
${{\left. \psi \right|}_{{\partial \Omega }}} = {{\left. {\frac{{\partial \psi }}{{\partial {{n}_{x}}}}} \right|}_{{\partial \Omega }}} = 0,$
где $\psi $ – это потенциал магнитного поля.

Давайте отдельно обсудим возникновение нелинейных граничных условий для потенциала электрического поля. Действительно, предположим, что на границе кристаллического полупроводника $\partial \Omega $, занимающего область $\Omega \subset {{\mathbb{R}}^{3}}$, имеются либо стоки, либо источники поверхностных зарядов, тогда согласно граничным условиям общего вида получим граничное условие следующего вида:

(2.25)
${{\left. {({\mathbf{E}},{{n}_{x}})} \right|}_{{\partial \Omega }}} = {{\left. {\eta (x)} \right|}_{{\partial \Omega }}},$
где ${{n}_{x}}$ – это внешняя нормаль в точке $x \in \partial \Omega $. Предположим, что распределение поверхностных зарядов моделируется нелинейной степенной функцией вида:
(2.26)
$\eta (\phi (x)) = {{\eta }_{0}}{\text{|}}\phi (x){{{\text{|}}}^{q}}\phi (x),\quad q > 0,$
где ${{\eta }_{0}} > 0$ для стоков, а ${{\eta }_{0}} < 0$ для источников свободных зарядов.

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

2.1. Уравнение нелинейных нелокальных волн ББМБ с источником

Рассмотрим кристаллический полупроводник, занимающий ограниченную поверхностно односвязанную область $\Omega \subseteq {{\mathbb{R}}^{3}}$ с гладкой границей $\partial \Omega \in {{\mathbb{C}}^{{2,\delta }}}$ при $\delta \in (0,1]$ в случае $\Omega \ne {{\mathbb{R}}^{3}}$. Предположим, что распределение связанных зарядов в этом полупроводнике описывается формулой (2.4), т.е.

$\rho = \operatorname{div} {\mathbf{P}} = {{P}_{0}}{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ,\quad {{q}_{1}} \geqslant 0,\quad {{P}_{0}} > 0.$
Предположим, что этот полупроводник находится во внешнем постоянном и однородном электрическом поле ${{{\mathbf{E}}}_{0}}$, причем в эвклидовом пространстве ${{\mathbb{R}}^{3}}$ выбран ортонормированный репер $({{{\mathbf{e}}}_{1}},{{{\mathbf{e}}}_{2}},{{{\mathbf{e}}}_{3}})$ таким образом, чтобы
${{{\mathbf{E}}}_{0}} = {{E}_{0}}{{{\mathbf{e}}}_{1}},\quad {{E}_{0}} \in {{\mathbb{R}}^{1}}.$
Помимо всего прочего предположим, что в полупроводнике имеет место временная дисперсия, задаваемая формулой (2.6). Пусть, кроме того, в полупроводнике имеются распределенные источники свободных зарядов, задаваемые функцией $Q(\phi )$ следующего вида:
$Q(\phi ) = {{Q}_{0}}{\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi ,\quad {{q}_{2}} > 0,\quad {{Q}_{0}} > 0.$
Таким образом, с учетом результатов предыдущего раздела приходим к следующей системе уравнений:
(2.27)
$\operatorname{div} {\mathbf{D}} = - 4\pi n,\quad {\mathbf{D}} = {\mathbf{E}} + 4\pi {\mathbf{P}},\quad {\mathbf{E}} = - \nabla \phi ,$
$\operatorname{div} {\mathbf{P}} = {{P}_{0}}{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ,$
$\frac{{\partial n}}{{\partial t}} = {\text{div}}({{{\mathbf{J}}}_{0}} + {{{\mathbf{J}}}_{1}} + {{{\mathbf{J}}}_{2}}) + {{Q}_{0}}{\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi ,$
(2.28)
${{{\mathbf{J}}}_{0}} = {{\sigma }_{0}}{{n}_{0}}(\phi ){{{\mathbf{E}}}_{0}},\quad {{{\mathbf{J}}}_{1}} = {{\sigma }_{1}}\int\limits_0^t \,dsh(t - s){\mathbf{E}}(s),\quad {{{\mathbf{J}}}_{2}} = {{\sigma }_{2}}{\mathbf{E}},$
где “квазистационарное” распределение свободных зарядов ${{n}_{0}}(\phi )$ моделируется формулой (2.16), т.е. следующим образом:
${{n}_{0}}(\phi ) = 1 + {{a}_{1}}\phi + {{a}_{2}}{{\phi }^{2}},\quad {{a}_{1}},{{a}_{2}} \in {{\mathbb{R}}^{1}}.$
В “единичных” коэффициентах, т.е. когда не имеющие принципиального характера коэффициенты полагаются равными единице, из системы уравнений (2.27), (2.28) получим следующее уравнение:
(2.29)
$\frac{\partial }{{\partial t}}(\Delta \phi \; - \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ) + \Delta \phi + \frac{{\partial \phi }}{{\partial {{x}_{1}}}} + \phi \frac{{\partial \phi }}{{\partial {{x}_{1}}}} + \int\limits_0^t \,dsh(t - s)\Delta \phi (s)\; + \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi = 0,\quad {{q}_{1}} \geqslant 0,\quad {{q}_{2}} > 0,$
где $h(t) \in \mathbb{C}[0, + \infty )$, $x = ({{x}_{1}},{{x}_{2}},{{x}_{3}}) \in \Omega \subseteq {{\mathbb{R}}^{3}}.$ Уравнение (2.29) в случае $\Omega \ne {{\mathbb{R}}^{3}}$ дополним однородным граничным условием
(2.30)
${{\left. \phi \right|}_{{\partial \Omega }}} = 0,$
а также начальным условием

(2.31)
$\phi (x,0) = {{\phi }_{0}}(x).$

2.2. Уравнение нелинейных нелокальных волн ББМБ с нелокальным источником

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

(2.32)
$Q(\phi ) = - {{Q}_{0}}\mathop {\left[ {\int\limits_\Omega \,dx{\text{|}}\nabla \phi {{{\text{|}}}^{2}}} \right]}\nolimits^{{{q}_{2}}} \Delta \phi ,\quad {{q}_{2}} > 0.$
Давайте обсудим физические основы появления такого рода нелинейной и нелокальной зависимости. Действительно, рассмотрим кристаллический полупроводник при наличии отрицательной проводимости, т.е. когда связь векторов тока проводимости ${\mathbf{J}}$ и напряженности электрического поля ${\mathbf{E}}$ следующая:
(2.33)
${\mathbf{J}} = - {{\sigma }_{0}}{\mathbf{E}},\quad {{\sigma }_{0}} > 0.$
Как известно, величина ${{\sigma }_{0}}$ имеет физический смысл коэффициента проводимости среды, который, вообще говоря, зависит от средней по области температуры. Заметим, что эта усредненная температура имеет следующий вид:
$T = {{T}_{0}} + \frac{1}{{{\text{|}}\Omega {\text{|}}}}\int\limits_\Omega {\text{|}}\nabla \phi {{{\text{|}}}^{2}}dx,$
где ${{T}_{0}}$ – это температура фононов. Положив температуру фононов ${{T}_{0}} = 0$ и взяв модельную зависимость ${{\sigma }_{0}} = {{\sigma }_{0}}(T)$ вида
(2.34)
${{\sigma }_{0}} = {{a}_{0}}{{T}^{{{{q}_{2}}}}},\quad {{q}_{2}} > 0,$
мы получим зависимость вида (2.32). Действительно, из (2.33), (2.34) вытекает следующая цепочка равенств:
$Q(\phi ) = - \operatorname{div} {\mathbf{J}} = {{a}_{0}}{{\sigma }_{0}}{{T}^{{{{q}_{2}}}}}\operatorname{div} {\mathbf{E}} = - {{a}_{0}}{{\sigma }_{0}}\mathop {\left[ {\int\limits_\Omega dx{\text{|}}\nabla \phi {{{\text{|}}}^{2}}} \right]}\nolimits^{{{q}_{2}}} \Delta \phi ,$
где ${\mathbf{E}} = - \nabla \phi .$ Таким образом, мы пришли к следующему уравнению:
$\begin{gathered} \frac{\partial }{{\partial t}}(\Delta \phi \; - \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ) + \Delta \phi + \frac{{\partial \phi }}{{\partial {{x}_{1}}}} + \phi \frac{{\partial \phi }}{{\partial {{x}_{1}}}} + \int\limits_0^t \,dsh(t - s)\Delta \phi (s) - \mathop {\left[ {\int\limits_\Omega \,dx{\text{|}}\nabla \phi {{{\text{|}}}^{2}}} \right]}\nolimits^{{{q}_{2}}} \Delta \phi = 0, \\ {{q}_{1}} \geqslant 0,\quad {{q}_{2}} > 0, \\ \end{gathered} $
где $h(t) \in \mathbb{C}[0, + \infty )$, $x = ({{x}_{1}},{{x}_{2}},{{x}_{3}}) \in \Omega \subset {{\mathbb{R}}^{3}},$ которое нужно дополнить граничным и начальным условиями (2.30) и (2.31).

2.3. Нелокальное диссипативное уравнение типа Розенау–Бюргерса с источником

Теперь перейдем к рассмотрению кристаллического полупроводника при учете пространственной дисперсии. Итак, прежде всего учтем влияние двух факторов в векторе поляризуемости P:

(2.35)
${\mathbf{P}} = {{{\mathbf{P}}}_{1}} + {{{\mathbf{P}}}_{2}}.$
Во-первых, учтем сильную пространственную дисперсию
${{{\mathbf{P}}}_{1}} = - {{a}_{0}}\Delta {\mathbf{E}},\quad \Delta \equiv \frac{{{{\partial }^{2}}}}{{\partial x_{1}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{2}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{3}^{2}}},$
а во-вторых, учтем керровскую нелинейность вида
${{{\mathbf{P}}}_{2}} = {{a}_{1}}{\text{|}}{\mathbf{E}}{{{\text{|}}}^{{{{p}_{1}} - 2}}}{\mathbf{E}},\quad {{p}_{1}} \geqslant 2,\quad {{a}_{1}} > 0.$
Предположим, что кристаллический полупроводник, занимающий область $\Omega \subseteq {{\mathbb{R}}^{3}},$ имеет отрицательную дифференциальную проводимость, т.е. ток проводимости ${\mathbf{J}}$ можно представить в следующем виде:
${\mathbf{J}} = {{{\mathbf{J}}}_{1}} + {{{\mathbf{J}}}_{2}},$
где
${{{\mathbf{J}}}_{1}} = {{\sigma }_{0}}[1\; - \;{\text{|}}{\mathbf{E}}{{{\text{|}}}^{{{{p}_{2}} - 2}}}]{\mathbf{E}},\quad {{\sigma }_{0}} > 0,$
а вектор ${{{\mathbf{J}}}_{2}}$ учитывает временную дисперсию
${{{\mathbf{J}}}_{2}}(t) = \int\limits_0^t \,dsh(t - s){\mathbf{E}}(s).$
Таким образом, мы пришли к следующей системе уравнений:
(2.36)
$\operatorname{div} {\mathbf{D}} = - 4\pi n,\quad {\mathbf{D}} = {\mathbf{E}} + 4\pi {\mathbf{P}},\quad {\mathbf{E}} = - \nabla \phi ,$
(2.37)
$\frac{{\partial n}}{{\partial t}} = \operatorname{div} {\mathbf{J}},$
где
${\mathbf{P}} = - {{a}_{0}}\vartriangle {\mathbf{E}} + {{a}_{1}}{\text{|}}{\mathbf{E}}{{{\text{|}}}^{{{{p}_{1}} - 2}}}{\mathbf{E}},\quad {\mathbf{J}} = {{\sigma }_{0}}[1\; - \;{\text{|}}{\mathbf{E}}{{{\text{|}}}^{{{{p}_{2}} - 2}}}]{\mathbf{E}} + \int\limits_0^t \,dsh(t - s){\mathbf{E}}(s).$
Из системы уравнений (2.36), (2.37) приходим к следующему модельному уравнению пятого порядка с единичными коэффициентами:
$\begin{gathered} \frac{\partial }{{\partial t}}( - {{\Delta }^{2}}\phi + \Delta \phi + {{\Delta }_{{{{p}_{1}}}}}\phi ) + \Delta \phi + \int\limits_0^t \,dsh(t - s)\Delta \phi (s) - {{\Delta }_{{{{p}_{2}}}}}\phi = 0, \\ {{\Delta }_{p}}\phi \equiv {\text{div}}({\text{|}}\nabla \phi {{{\text{|}}}^{{p - 2}}}\nabla \phi ),\quad p \geqslant 2,\quad {{p}_{1}} \geqslant 2,\quad {{p}_{2}} > 2. \\ \end{gathered} $
Данное уравнение в случае $\Omega \ne {{\mathbb{R}}^{3}}$ дополним граничными условиями (2.23) и начальным условием

$\phi (x,0) = {{\phi }_{0}}(x).$

2.4. Нелинейное-нелокальное уравнение спиновых волн

Теперь мы приступим к выводу модельного нелокального-нелинейного уравнения спиновых волн в кристаллическом полупроводнике–магнетике (см. уравнения (2.18)(2.22)). Рассмотрим следующие факторы, влияющие на вид квазистационарной намагниченности $\mathfrak{M}:$

(2.38)
$\mathfrak{M} = {{\mathfrak{M}}_{1}} + {{\mathfrak{M}}_{2}} + {{\mathfrak{M}}_{3}},$
где слагаемое ${{\mathfrak{M}}_{1}}$ линейно зависит от напряженности магнитного поля:
${{\mathfrak{M}}_{{1j}}} = {{a}_{j}}{{H}_{j}},\quad {{a}_{j}} \in {{R}^{1}},\quad j = 1,2,3,$
слагаемое ${{\mathfrak{M}}_{2}}$ учитывает сильную пространственную дисперсию вида
${{\mathfrak{M}}_{2}} = - {{a}_{0}}\Delta {\mathbf{H}},\quad {{a}_{0}} > 0,\quad \Delta \equiv \frac{{{{\partial }^{2}}}}{{\partial x_{1}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{2}^{2}}} + \frac{{{{\partial }^{2}}}}{{\partial x_{3}^{2}}},$
наконец, слагаемое ${{\mathfrak{M}}_{3}}$ учитывает “керровскую” нелинейность вида
${{\mathfrak{M}}_{3}} = {{a}_{1}}{\text{|}}{\mathbf{H}}{{{\text{|}}}^{{{{p}_{1}} - 2}}}{\mathbf{H}},\quad {{a}_{1}} > 0,\quad {{p}_{1}} \geqslant 2.$
Заметим, что
$\left[ {{{\mathfrak{M}}_{2}},{\mathbf{H}}} \right] = \left[ {{{\mathfrak{M}}_{3}},{\mathbf{H}}} \right] = 0.$
Величина ${\mathbf{R}}$ из формулы (2.21) учитывает два фактора. Первый – это наличие временной дисперсии
${{{\mathbf{R}}}_{1}}(t) = \int\limits_0^t \,dsh(t - s){\mathbf{H}}(s),$
второй – это наличие “источников” магнитных доменов
${{{\mathbf{R}}}_{2}} = - {{a}_{2}}{\text{|}}{\mathbf{H}}{{{\text{|}}}^{{{{p}_{2}} - 2}}}{\mathbf{H}},\quad {{p}_{2}} > 2.$
Таким образом,
${\mathbf{R}} = {{{\mathbf{R}}}_{1}} + {{{\mathbf{R}}}_{2}}.$
Следовательно, если кристаллический полупроводник–магнетик занимает ограниченную поверхностно-односвязанную область $\Omega \subset {{\mathbb{R}}^{3}},$ то при учете указанных факторов приходим к следующей системе уравнений
$\operatorname{div} {\mathbf{B}} = 0,\quad \operatorname{rot} {\mathbf{H}} = 0,\quad {\mathbf{B}} = {\mathbf{H}} + 4\pi {\mathbf{M}},\quad {\mathbf{H}} = - \nabla \psi ,$
$\frac{{\partial {\mathbf{M}}}}{{\partial t}} = \gamma \left[ {{\mathbf{H}},\mathfrak{M}} \right] + {\mathbf{R}},\quad {\mathbf{M}} = \mathfrak{M} + {\mathbf{m}}.$
Займемся арифметикой:
(2.39)
$\operatorname{div} [{\mathbf{H}},{{\mathfrak{M}}_{1}}] = - ({\mathbf{H}},\operatorname{rot} {{\mathfrak{M}}_{1}}) = - {{\beta }_{1}}\frac{\partial }{{\partial {{x}_{1}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{2}}}}\frac{{\partial \psi }}{{\partial {{x}_{3}}}}} \right) - {{\beta }_{2}}\frac{\partial }{{\partial {{x}_{2}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{3}}}}\frac{{\partial \psi }}{{\partial {{x}_{1}}}}} \right) - {{\beta }_{3}}\frac{\partial }{{\partial {{x}_{3}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{1}}}}\frac{{\partial \psi }}{{\partial {{x}_{2}}}}} \right),$
где
${{\beta }_{1}} = {{a}_{2}} - {{a}_{3}},\quad {{\beta }_{2}} = {{a}_{3}} - {{a}_{1}},\quad {{\beta }_{3}} = {{a}_{1}} - {{a}_{2}},\quad {\text{|}}{{\beta }_{1}}{\text{|}} + \;{\text{|}}{{\beta }_{2}}{\text{|}} + \;{\text{|}}{{\beta }_{3}}{\text{|}} > 0.$
Заметим, что
${{\beta }_{1}} + {{\beta }_{2}} + {{\beta }_{3}} = 0.$
С учетом (2.18)–(2.22) и (2.38), (2.39) приходим к следующему уравнению пятого порядка с единичными коэффициентами:
(2.40)
$\begin{gathered} \frac{\partial }{{\partial t}}( - {{\Delta }^{2}}\psi + \Delta \psi + {{\Delta }_{{{{p}_{1}}}}}\psi ) + \int\limits_0^t \,dsh(t - s)\Delta \psi (s) + {{\alpha }_{1}}\frac{\partial }{{\partial {{x}_{1}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{2}}}}\frac{{\partial \psi }}{{\partial {{x}_{3}}}}} \right) + {{\alpha }_{2}}\frac{\partial }{{\partial {{x}_{2}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{3}}}}\frac{{\partial \psi }}{{\partial {{x}_{1}}}}} \right) + \\ + \;{{\alpha }_{3}}\frac{\partial }{{\partial {{x}_{3}}}}\left( {\frac{{\partial \psi }}{{\partial {{x}_{1}}}}\frac{{\partial \psi }}{{\partial {{x}_{2}}}}} \right) - {{\Delta }_{{{{p}_{2}}}}}\psi = 0,\quad {{\Delta }_{p}}\psi \equiv {\text{div}}({\text{|}}\nabla \psi {{{\text{|}}}^{{p - 2}}}\nabla \psi ),\quad p \geqslant 2,\quad {{p}_{1}} \geqslant 2,\quad {{p}_{2}} > 2, \\ \end{gathered} $
где ${{\alpha }_{1}} + {{\alpha }_{2}} + {{\alpha }_{3}} = 0,$ но ${\text{|}}{{\alpha }_{1}}{\text{|}} + \;{\text{|}}{{\alpha }_{2}}{\text{|}} + \;{\text{|}}{{\alpha }_{3}}{\text{|}} > 0.$ Уравнение (2.40) в случае $\Omega \ne {{\mathbb{R}}^{3}}$ дополним граничными условиями (2.24), а также начальным условием

$\psi (x,0) = {{\psi }_{0}}(x).$

2.5. Одна нелинейная-нелокальная система уравнений

Теперь рассмотрим систему уравнений (2.35)–(2.37) в смысле комплекснозначных функций, т.е. когда

${\mathbf{D}},\;{\mathbf{E}},\;{\mathbf{P}},\;{\mathbf{J}} \in {{\mathbb{C}}^{3}} \equiv {{\mathbb{R}}^{3}} + i{{\mathbb{R}}^{3}},$
$\phi = {{\phi }_{1}} + i\phi ,\quad n = {{n}_{1}} + i{{n}_{2}},$
однако, функцию $h(t)$ будем считать вещественной. Тогда относительно вещественных потенциалов ${{\phi }_{1}}$ и ${{\phi }_{2}}$ получим следующую систему уравнений:
(2.41)
$\begin{gathered} \frac{\partial }{{\partial t}}( - {{\Delta }^{2}}{{\phi }_{1}} + \Delta {{\phi }_{1}} + {\text{div}}({{[{\text{|}}\nabla {{\phi }_{1}}{{{\text{|}}}^{2}} + \;{\text{|}}\nabla {{\phi }_{2}}{{{\text{|}}}^{2}}]}^{{({{p}_{1}} - 2)/2}}}\nabla {{\phi }_{1}})) + \Delta {{\phi }_{1}} + \\ + \;\int\limits_0^t \,dsh(t - s)\Delta {{\phi }_{1}}(s) - {\text{div}}({{[{\text{|}}\nabla {{\phi }_{1}}{{{\text{|}}}^{2}} + \;{\text{|}}\nabla {{\phi }_{2}}{{{\text{|}}}^{2}}]}^{{({{p}_{2}} - 2)/2}}}\nabla {{\phi }_{1}}) = 0, \\ \end{gathered} $
(2.42)
$\begin{gathered} \frac{\partial }{{\partial t}}( - {{\Delta }^{2}}{{\phi }_{2}} + \Delta {{\phi }_{2}} + {\text{div}}({{[{\text{|}}\nabla {{\phi }_{1}}{{|}^{2}} + \;{\text{|}}\nabla {{\phi }_{2}}{{{\text{|}}}^{2}}]}^{{({{p}_{1}} - 2)/2}}}\nabla {{\phi }_{2}})) + \Delta {{\phi }_{2}} + \\ + \;\int\limits_0^t \,dsh(t - s)\Delta {{\phi }_{2}}(s) - {\text{div}}({{[{\text{|}}\nabla {{\phi }_{1}}{{{\text{|}}}^{2}} + \;{\text{|}}\nabla {{\phi }_{2}}{{{\text{|}}}^{2}}]}^{{({{p}_{2}} - 2)/2}}}\nabla {{\phi }_{2}}) = 0, \\ \end{gathered} $
где ${{p}_{1}} \geqslant 2,\;{{p}_{2}} > 2.$ Систему уравнений (2.41), (2.42) в случае $\Omega \ne {{\mathbb{R}}^{3}}$ необходимо дополнить граничными условиями
${{\left. {{{\phi }_{1}}} \right|}_{{\partial \Omega }}} = {{\left. {\frac{{\partial {{\phi }_{1}}}}{{\partial {{n}_{x}}}}} \right|}_{{\partial \Omega }}} = {{\left. {{{\phi }_{2}}} \right|}_{{\partial \Omega }}} = {{\left. {\frac{{\partial {{\phi }_{2}}}}{{\partial {{n}_{x}}}}} \right|}_{{\partial \Omega }}} = 0,$
а также начальными условиями

${{\phi }_{1}}(x,0) = {{\phi }_{{10}}}(x),\quad {{\phi }_{2}}(x,0) = {{\phi }_{{20}}}(x).$

2.6. Нелинейная-нелокальная система уравнений А.П. Осколкова с источником

Приведем без вывода (см., например, [25]) одну важную нелинейную и нелокальную систему уравнений А.П. Осколкова с источником следующего вида:

(2.43)
$\frac{\partial }{{\partial t}}\left( {\Delta {\mathbf{u}} - {\mathbf{u}}} \right) + \Delta {\mathbf{u}} + \int\limits_0^t \,h(t - s)\Delta {\mathbf{u}}(s)ds + ({\mathbf{u}},\nabla ){\mathbf{u}} + \mathop {\left| {\mathbf{u}} \right|}\nolimits^2 {\mathbf{u}} = \nabla p,$
где $p$ – давление в жидкости, ${\mathbf{u}} = ({{u}_{1}},{{u}_{2}},{{u}_{3}})$ – вектор скорости в жидкости. Систему уравнений (2.43) необходимо дополнить условием несжимаемости жидкости:
$\operatorname{div} {\mathbf{u}} = 0,$
а также в случае $\Omega \ne {{\mathbb{R}}^{3}}$ граничным условием
${{\left. {\mathbf{u}} \right|}_{{\partial \Omega }}} = 0$
и, наконец, начальным условием
${\mathbf{u}}(x,0) = {{{\mathbf{u}}}_{0}}(x),$
где $x = ({{x}_{1}},{{x}_{2}},{{x}_{3}}) \in \Omega \subset {{\mathbb{R}}^{3}},$ а область $\Omega $ является ограниченной с гладкой границей $\partial \Omega \in {{\mathbb{C}}^{{2,\delta }}}$ при $\delta \in (0,1].$

2.7. Нелокальное-нелинейное уравнение с нелинейным граничным условием

Теперь рассмотрим кристаллический полупроводник в ограниченной области $\Omega \subset {{\mathbb{R}}^{3}}$ при учете распределения связанных зарядов вида

(2.44)
$\rho = \operatorname{div} {\mathbf{P}} = \phi ,$
источников свободных зарядов, распределение которых описывается функцией $Q(\phi )$ следующего вида:
$Q(\phi ) = {\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi ,\quad {{q}_{2}} > 0,$
и, наконец, при учете тока проводимости вида
(2.45)
${\mathbf{J}} = {{\sigma }_{0}}{\mathbf{E}} + \int\limits_0^t \,dsh(t - s){\mathbf{E}}.$
Предположим, что на границе $\partial \Omega $ области $\Omega $ имеются “стоки” свободных зарядов. В этом случае согласно (2.25) и (2.26) получим граничное условие
${{\left. {\left( {\frac{{\partial \phi }}{{\partial {{n}_{x}}}}\; + \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi } \right)} \right|}_{{\partial \Omega }}} = 0,\quad {{q}_{1}} \geqslant 0,$
а для потенциала $\phi $ в области $\Omega $ из условий (2.44), (2.45) вытекает следующее уравнение:
$\frac{\partial }{{\partial t}}(\Delta \phi - \phi ) + \Delta \phi + \int\limits_0^t \,dsh(t - s)\Delta \phi (s)\; + \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi = 0,$
которое осталось дополнить начальным условием

$\phi (x,0) = {{\phi }_{0}}(x).$

2.8. Нелинейное уравнение с нелинейным эволюционным нелокальным граничным условием

Теперь рассмотрим кристаллический полупроводник в ограниченной области $\Omega \subset {{\mathbb{R}}^{3}}$ при учете распределения связанных зарядов вида

(2.46)
$\rho = \operatorname{div} {\mathbf{P}} = \phi \; + \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ,$
источников свободных зарядов, распределение которых описывается функцией $Q(\phi )$ следующего вида:
(2.47)
$Q(\phi ) = {\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi ,\quad {{q}_{2}} > 0,$
и при учете тока проводимости вида
${\mathbf{J}} = {{\sigma }_{0}}{\mathbf{E}}.$
Наконец, учтем наличие распределения свободных зарядов на границе $\partial \Omega $ области $\Omega $ так, что ток проводимости на границе связан с распределением свободных зарядов вида
(2.48)
$\frac{\partial }{{\partial t}}({\mathbf{J}},{{n}_{x}}) + ({\mathbf{J}},{{n}_{x}}) = \int\limits_0^t \,dsh(t - s)f(x,\phi )(s).$
Из (2.47), (2.48) вытекает следующая модельная задача:
$\frac{\partial }{{\partial t}}(\Delta \phi - \phi \; - \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{1}}}}}\phi ) + \Delta \phi + \;{\text{|}}\phi {{{\text{|}}}^{{{{q}_{2}}}}}\phi = 0,$
${{\left. {\left( {\frac{\partial }{{\partial t}}\frac{{\partial \phi }}{{\partial {{n}_{x}}}} + \frac{{\partial \phi }}{{\partial {{n}_{x}}}} + \int\limits_0^t dsh(t - s)f(x,\phi )(s)} \right)} \right|}_{{\partial \Omega }}} = 0,$
которую дополним начальным условием

$\phi (x,0) = {{\phi }_{0}}(x).$

3. МОДЕЛЬНЫЕ УРАВНЕНИЯ ДЛЯ ФУНКЦИИ h(t)

Во всех полученных в предыдущем разделе уравнениях содержится слагаемое вида

$\int\limits_0^t \,dsh(t - s)\Delta \phi (s),$
описывающее переходные процессы в кристаллических полупроводниках.

В этом разделе мы рассмотрим некоторые модельные выражения для этой функции $h(t)$. Прежде всего заметим, что в качестве этой функции можно взять следующее выражение:

(3.1)
$h(t) = {{h}_{0}}{{e}^{{ - at}}},\quad {{h}_{0}} > 0,\quad a > 0.$
Действительно, формула (3.1) описывает самый распространенный механизм релаксации. Кроме того, можно рассмотреть следующее уравнение:
$h(t) = {{h}_{0}}{{t}^{n}}{{e}^{{ - at}}},\quad {{h}_{0}} > 0,\quad a > 0,\quad n \in N.$
Наконец, можно учесть колебательный режим релаксации следующего вида:
$h(t) = {{h}_{0}}sin({{\omega }_{0}}t){{e}^{{ - at}}},\quad {{h}_{0}} > 0,\quad a > 0,\quad {{\omega }_{0}} > 0,$
или
$h(t) = {{h}_{0}}cos({{\omega }_{0}}t){{e}^{{ - at}}},\quad {{h}_{0}} > 0,\quad a > 0,\quad {{\omega }_{0}} > 0.$
С другой стороны, в физике обычно релаксационные процессы протекают лишь конечное время, т.е. функция $h(t)$ финитна:
$h(t) = \left\{ \begin{gathered} \omega (t)\quad {\text{п р и }}\quad 0 \leqslant t \leqslant {{t}_{0}}; \hfill \\ 0\quad {\text{п р и }}\quad {{t}_{0}} \leqslant t \leqslant + \infty , \hfill \\ \end{gathered} \right.$
где функция $\omega (t) \in {\mathbf{C}}_{0}^{\infty }(R_{ + }^{1})$. В следующих разделах мы будем изучать полученные в этом разделе уравнения при обобщении следующего важного условия, что найдутся постоянные a > 0 и h0 > 0 такие, что

$h(t) \in {{\mathbb{C}}^{{(1)}}}(R_{ + }^{1}),\quad \max \{ \left| {h(t)} \right|,\;\left| {h{\kern 1pt} '(t)} \right|\} \leqslant {{h}_{0}}{{e}^{{ - at}}}\quad {\text{п р и }}\quad t \geqslant 0.$

4. ПРЕДВАРИТЕЛЬНЫЕ УСЛОВИЯ

Рассмотрим банаховы пространства ${{V}_{0}},$ ${{V}_{j}},$ ${{W}_{i}}$ при $j = \overline {1,n} $ и при $i = \overline {1,4} $ относительно соответствующих норм

${{\left\| {\, \cdot \,} \right\|}_{0}},\quad {{\left\| {\, \cdot \,} \right\|}_{j}},\quad {{\left| {\, \cdot \,} \right|}_{i}}$
и с сопряженными банаховыми пространствами $V_{0}^{ * },$ $V_{j}^{ * },$ $W_{i}^{ * }$ относительно соответствующих скобок двойственностей
${{\langle \cdot , \cdot \rangle }_{0}},\quad {{\langle \cdot , \cdot \rangle }_{j}},\quad {{( \cdot , \cdot )}_{i}}$
и соответствующих норм
$\left\| {\, \cdot \,} \right\|_{0}^{ * },\quad \left\| {\, \cdot \,} \right\|_{j}^{ * },\quad {\text{|}}\, \cdot \,{\text{|}}_{i}^{ * }.$
Предположим, что банаховы пространства ${{V}_{0}},$ ${{V}_{j}},$ ${{W}_{i}}$ при $j = \overline {1,n} $ и при $i = \overline {1,4} $ являются рефлексивными и сепарабельными. Предположим также, что
${{A}_{0}}:{{V}_{0}} \to V_{0}^{*},\quad {{A}_{j}}:{{V}_{j}} \to V_{j}^{*},\quad {{L}_{1}}:{{W}_{1}} \to W_{1}^{*},\quad {{L}_{2}}:{{W}_{2}} \to W_{2}^{*},$
$F:{{W}_{3}} \to W_{3}^{*},\quad P:{{V}_{0}} \to {{W}_{4}},\quad D:{{W}_{4}} \to V_{0}^{*}.$
Теперь последовательно укажем условия на все операторные коэффициенты задачи (1.2). Итак,

Условия ${{A}_{0}}{\text{:}}$

(i) оператор ${{A}_{0}}:{{V}_{0}} \to V_{0}^{*}$ является линейным, непрерывным и симметричным, причем имеет место неравенство

$\left\| {{{A}_{0}}u} \right\|_{0}^{*} \leqslant {{M}_{0}}{{\left\| u \right\|}_{0}}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{V}_{0}};$

(ii) оператор ${{A}_{0}}$ является коэрцитивным, причем имеет место неравенство

${{\langle {{A}_{0}}u,u\rangle }_{0}} \geqslant {{m}_{0}}\left\| u \right\|_{0}^{2}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{V}_{0}};$

(iii) величина $\langle {{A}_{0}}u,u\rangle _{0}^{{1/2}}$ является нормой на ${{V}_{0}},$ порождающей рассматриваемую топологию банахова пространства ${{V}_{0}}.$

Условия $A{\text{:}}$

(i) оператор ${{A}_{j}}:{{V}_{j}} \to V_{j}^{*}$ является монотонным и непрерывным оператором;

(ii) оператор ${{A}_{j}}$ дифференцируем по Фреше, причем его производная Фреше

$A_{{if}}^{'}(u):{{V}_{j}} \to \mathcal{L}({{V}_{j}},V_{j}^{*})$
является непрерывным, симметричным, монотонным и неотрицательно-определенным оператором при фиксированном $u \in {{V}_{j}};$

(iii) оператор ${{A}_{j}}$ является положительно-однородным

${{A}_{j}}(ru) = {{r}^{{{{p}_{j}} - 1}}}{{A}_{j}}(u)\quad {\text{п р и }}\quad {{p}_{j}} > 2,\quad r \geqslant 0,\quad u \in {{V}_{j}};$

(iv) справедливы следующие неравенства сверху и снизу:

$\left\| {{{A}_{j}}(u)} \right\|_{j}^{*} \leqslant {{M}_{j}}\left\| u \right\|_{j}^{{{{p}_{j}} - 1}},\quad {{\langle {{A}_{j}}(u),u\rangle }_{j}} \geqslant {{m}_{j}}\left\| u \right\|_{j}^{{{{p}_{j}}}},\quad {{M}_{j}},{{m}_{j}} > 0;$

(v) величина $\langle {{A}_{j}}(u),u\rangle _{j}^{{1/{{p}_{j}}}}$ является нормой на банаховом пространстве ${{V}_{j}},$ порождающей рассматриваемую топологию банахова пространства ${{V}_{j}}.$

Условия $F$:

(i) оператор $F:{{W}_{3}} \to W_{3}^{*}$ является ограниченно липшиц-непрерывным, т.е. имеет место неравенство

$\mathop {\left| {F({{u}_{1}}) - F({{u}_{2}})} \right|}\nolimits_3 \leqslant \mu (R){\text{|}}{{u}_{1}} - {{u}_{2}}{{{\text{|}}}_{3}}\quad {\text{д л я }}\;{\text{в с е х }}\quad {{u}_{1}},{{u}_{2}} \in {{W}_{3}},$
где $\mu = \mu (R)$ есть ограниченная но всяком компакте неубывающая функция своего аргумента;

(ii) оператор $F$ является положительно однородным, т.е.

$F(ru) = {{r}^{{1 + {{q}_{2}}}}}F(u)\quad {\text{п р и }}\quad {{q}_{2}} > 0,\quad r \geqslant 0,\quad u \in {{W}_{3}};$

(iii) оператор $F$ имеет симметричную производную Фреше

$F_{f}^{'}({\kern 1pt} \cdot {\kern 1pt} ):{{W}_{3}} \to \mathcal{L}({{W}_{3}},W_{3}^{*});$

(iv) оператор $F$ удовлетворяет неравенству сверху

${\text{|}}F(u){\text{|}}_{3}^{*} \leqslant M{\text{|}}u{\text{|}}_{3}^{{{{q}_{2}} + 1}}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{W}_{3}}.$

Условия $L{\text{:}}$

(i) оператор ${{L}_{i}}:{{W}_{i}} \to W_{i}^{*}$ при $i = 1,2$ является линейным, непрерывным и симметричным, причем

$\left| {{{L}_{i}}u} \right|_{i}^{*} \leqslant {{D}_{i}}{\text{|}}u{{{\text{|}}}_{i}}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{W}_{i}};$

(ii) оператор ${{L}_{i}}$ является коэрцитивным, причем

${{({{L}_{i}}u,u)}_{i}} \geqslant {{d}_{i}}{\text{|}}u{\text{|}}_{i}^{2}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{W}_{i}};$

(iii) величина $({{L}_{i}}u,u)_{i}^{{1/2}}$ является нормой на ${{W}_{i}},$ порождающей рассматриваемую топологию банахова пространства ${{W}_{i}}.$

Условия $DP{\text{:}}$

(i) оператор $D:{{W}_{4}} \to V_{0}^{*}$ является линейным и непрерывным, причем имеет место неравенство сверху

$\left\| {Du} \right\|_{0}^{*} \leqslant {{D}_{3}}{\text{|}}u{{{\text{|}}}_{4}}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{W}_{4}};$

(ii) оператор $P:{{V}_{0}} \to {{W}_{4}}$ является ограниченно липшиц-непрерывным, т.е.

${\text{|}}P({{u}_{1}}) - P({{u}_{2}}){{{\text{|}}}_{4}} \leqslant {{\mu }_{0}}(R){{\left\| {{{u}_{1}} - {{u}_{2}}} \right\|}_{0}}\quad {\text{д л я }}\;{\text{в с е х }}\quad {{u}_{1}},{{u}_{2}} \in {{V}_{0}},$
где функция ${{\mu }_{0}} = {{\mu }_{0}}(R)$ – ограниченная на всяком компакте неубывающая функция своего аргумента, $R = max\{ {{\left\| {{{u}_{1}}} \right\|}_{0}},{{\left\| {{{u}_{2}}} \right\|}_{0}}\} $;

(iii) справедливо неравенство сверху

${\text{|}}P(u){{{\text{|}}}_{4}} \leqslant {{D}_{4}}\left\| u \right\|_{0}^{{1 + {{q}_{3}}}},\quad {{q}_{3}} \geqslant 0\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{V}_{0}};$

(iv) для всех $u \in {{V}_{0}}$ имеет место неравенство

${{\langle DP(u),u\rangle }_{0}} = 0\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{V}_{0}}.$

Рассмотрим теперь используемые нами банаховы пространства ${{V}_{0}},$ ${{V}_{j}},$ ${{W}_{i}}$ при $j = \overline {1,n} $ и $i = \overline {1,4} $. Пусть $H$ – некоторое сепарабельное гильбертово пространство, отождествленное со своим сопряженным. Предположим, что выполнены следующие условия.

Условия $H{\text{:}}$

(i) имеют место следующие цепочки плотных и непрерывных вложений:

${{V}_{0}}\;\mathop \subset \limits^{ds} \;{{V}_{j}}\;\mathop \subset \limits^{ds} \;H\;\mathop \subset \limits^{ds} \;V_{j}^{*}\;\mathop \subset \limits^{ds} \;V_{0}^{*}\quad {\text{п р и }}\quad j = \overline {1,n} ,$
${{V}_{0}}\;\mathop \subset \limits^{ds} \;{{W}_{i}}\;\mathop \subset \limits^{ds} \;H\;\mathop \subset \limits^{ds} \;W_{i}^{*}\;\mathop \subset \limits^{ds} \;V_{0}^{*}\quad {\text{п р и }}\quad i = 1,2,3;$

(ii) имеют место непрерывные вложения

${{V}_{0}} \subset {{W}_{4}} \subset H \subset W_{4}^{*} \subset V_{0}^{*}.$

Заметим, что в силу пункта (i) из условий $H$ приходим к следующим равенствам скобок двойственностей

${{\langle f,u\rangle }_{0}} = {{\langle f,u\rangle }_{j}}\quad {\text{д л я }}\;{\text{в с е х }}\quad f \in V_{j}^{*},\quad u \in {{V}_{0}},\quad {\text{п р и }}\quad j = \overline {1,n} ,$
${{\langle f,u\rangle }_{0}} = {{(f,u)}_{i}}\quad {\text{д л я }}\;{\text{в с е х }}\quad f \in W_{i}^{*},\quad u \in {{V}_{0}},\quad {\text{п р и }}\quad i = \overline {1,3} .$

Наконец, потребуем выполнения следующих условий.

Условия h(t):

(i) функция $h(t) \in {{\mathbb{C}}^{{(1)}}}[0, + \infty )$;

(ii) функция $h(t)$ удовлетворяет неравенствам

$max\left\{ {{\text{|}}h(t){\text{|}},\left| {h{\text{'}}(t)} \right|} \right\} \leqslant \chi {\text{'}}(t),\quad \chi (t) \leqslant 0,\quad \chi (t) \in {{\mathbb{C}}^{{(1)}}}[0, + \infty );$

(iii) для функции $\chi {\text{'}}(t)$ выполнено следующее свойство:

$\int\limits_0^{ + \infty } \chi {\text{'}}(t)dt < + \infty .$

В дальнейшем мы будем предполагать, что выполнены все указанные условия.

5. ВСПОМОГАТЕЛЬНЫЕ РЕЗУЛЬТАТЫ

Прежде всего рассмотрим следующие функционалы:

(5.1)
${{\psi }_{1}}(u) = {{({{L}_{1}}u,u)}_{1}},\quad {{\psi }_{2}}(u) = {{({{L}_{2}}u,u)}_{2}},\quad {{\psi }_{3}}(u) = {{(F(u),u)}_{3}},$
и, следующий важный функционал
(5.2)
$\Phi (u) = \frac{1}{2}\mathop {\left\langle {{{A}_{0}}u,u} \right\rangle }\nolimits_0 + \sum\limits_{j = 1}^N \frac{{{{p}_{j}} - 1}}{{{{p}_{j}}}}\mathop {\left\langle {{{A}_{j}}(u),u} \right\rangle }\nolimits_j .$
Справедлива (см. [18, лемма 4.1])

Лемма 1. Если оператор $A:X \to X{\text{*}}$ дифференцируем по Фреше и имеет симметричную производную Фреше

$A_{u}^{'}(u):X \to \mathcal{Z}(X,X*)$
и $A(su) = {{s}^{{p - 1}}}A(u)$ при условиях $s \geqslant 0$ и $p \geqslant 2$, где $X$ – это пространство Банаха с сопряженным $X{\text{*}}$ относительно скобок двойственности $\langle \cdot , \cdot \rangle ,$ тогда функционал
$\psi (u) \equiv \left\langle {A(u),u} \right\rangle :X \to {{R}^{1}}$
является непрерывно дифференцируемым по Фреше и его производная Фреше есть

$\psi _{f}^{'}(u) = pA(u)\quad д л я \;в с е х \quad u \in X.$

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

$\psi _{{1f}}^{'}(u) = 2{{L}_{1}}u,\quad \psi _{{2f}}^{'}(u) = 2{{L}_{2}}u,\quad \psi _{{3f}}^{'}(u) = ({{q}_{2}} + 1)F(u),$
$\Phi _{f}^{'}(u) = {{A}_{0}}u + \sum\limits_{j = 1}^n \,({{p}_{j}} - 1){{A}_{j}}(u).$

Имеет место следующая (см. лемма 4.2 работы [18]):

Лемма 2. Пусть выполнены все условия леммы 1 и предположим, что $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];X)$ для некоторого $T > 0,$ тогда функционал

$\psi (u)(t) \equiv \left\langle {A(u),u} \right\rangle \in {{\mathbb{C}}^{{(1)}}}([0,T]).$

Но тогда в силу этой леммы и цепочек плотных и непрерывных вложений из Условия $H$ приходим к выводу о том, что на функциях

$u(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];{{V}_{0}})\quad {\text{п р и }}\quad T > 0$
функционалы

${{\psi }_{1}}(u)(t),\quad {{\psi }_{2}}(u)(t),\quad {{\psi }_{3}}(u)(t),\quad \Phi (u)(t) \in {{\mathbb{C}}^{{(1)}}}([0,T]).$

Справедливо утверждение (лемма 4.3 из [18]):

Лемма 3. Пусть выполнены все условия леммы 1 и предположим, что $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];X)$ для некоторого $T > 0,$ тогда

$\langle (A(u)){\text{'}},u\rangle = \frac{{p - 1}}{p}\frac{d}{{dt}}\langle A(u),u\rangle .$

Поэтому имеют место следующие равенства для функционалов:

(5.3)
$\mathop {\left\langle {\left( {{{A}_{j}}(u)} \right)'u} \right\rangle }\nolimits_j = \frac{{{{p}_{j}} - 1}}{{{{p}_{j}}}}\frac{d}{{dt}}\mathop {\left\langle {{{A}_{j}}(u),u} \right\rangle }\nolimits_j \quad {\text{п р и }}\quad j = \overline {1,n} .$
Справедлива следующая

Лемма 4. Имеет место неравенство

(5.4)
$\left| {{{{\langle DP(u),v\rangle }}_{0}}} \right| \leqslant \frac{\varepsilon }{2}{{\langle {{A}_{0}}v,v\rangle }_{0}} + \frac{{c_{1}^{2}}}{{2\varepsilon }}\langle {{A}_{0}}u,u\rangle _{0}^{{1 + {{q}_{3}}}}\quad д л я \;в с е х \quad u,v \in {{V}_{0}},$
где

${{c}_{1}} = \frac{{{{D}_{3}}{{D}_{4}}}}{{m_{0}^{{1 + {{q}_{3}}/2}}}},\quad \varepsilon > 0.$

Доказательство. Имеет место следующая цепочка неравенств:

$\begin{gathered} \left| {{{{\left\langle {DP(u),v} \right\rangle }}_{0}}} \right| \leqslant \left\| {DP(u)} \right\|_{0}^{ * }\left\| {{{v}_{0}}} \right\| \leqslant {{D}_{3}}{\text{|}}P(u){{{\text{|}}}_{4}}\frac{1}{{m_{0}^{{1/2}}}}\left\langle {{{A}_{0}}v,v} \right\rangle _{0}^{{1/2}} \leqslant \\ \leqslant \;{{D}_{3}}{{D}_{4}}\frac{1}{{m_{0}^{{1/2}}}}\frac{1}{{m_{0}^{{(1 + {{q}_{3}})/2}}}}\left\langle {{{A}_{0}}v,v} \right\rangle _{0}^{{1/2}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{(1 + {{q}_{3}})/2}} = \\ = \;{{c}_{1}}\left\langle {{{A}_{0}}v,v} \right\rangle _{0}^{{1/2}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{(1 + {{q}_{3}})/2}}\frac{\varepsilon }{2}{{\left\langle {{{A}_{0}}v,v} \right\rangle }_{0}} + \frac{{c_{1}^{2}}}{{2\varepsilon }}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1 + {{q}_{3}}}}, \\ \end{gathered} $
здесь мы воспользовались Условием ${{A}_{0}}$ пункта (ii), а также Условием $DP$.

Лемма доказана.

Справедливо следующее утверждение.

Лемма 5. Имеют место неравенства

(5.5)
$\left| {{{{({{L}_{i}}u,v)}}_{i}}} \right| \leqslant {{l}_{i}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}}\left\langle {{{A}_{0}}v,v} \right\rangle _{0}^{{1/2}} \leqslant \frac{\varepsilon }{2}{{\left\langle {{{A}_{0}}v,v} \right\rangle }_{0}} + \frac{{l_{i}^{2}}}{{2\varepsilon }}{{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}}\quad д л я \;в с е х \quad u,v \in {{V}_{0}},\quad \varepsilon > 0.$

Доказательство. В силу свойства (iii) Условия $L$ мы приходим к неравенству Шварца

$\left| {{{{({{L}_{i}}u,v)}}_{i}}} \right| \leqslant ({{L}_{i}}u,u)_{i}^{{1/2}}({{L}_{i}}v,v)_{i}^{{1/2}}.$
С другой стороны, имеет место по Условию $H$ непрерывное вложение $H$. Наконец, в силу Условия ${{A}_{0}}$ п. (iii) величина $\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}}$ есть норма на ${{V}_{0}},$ поэтому имеет место неравенство
$({{L}_{i}}u,u)_{i}^{{1/2}} \leqslant l_{i}^{{1/2}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}}\quad {\text{д л я }}\;{\text{в с е х }}\quad u \in {{V}_{0}}.$
Осталось воспользоваться неравенством Коши–Буняковского с “$\varepsilon $”.

Лемма доказана.

6. ЛОКАЛЬНАЯ РАЗРЕШИМОСТЬ В СИЛЬНОМ ОБОБЩЕННОМ СМЫСЛЕ

Дадим определение сильного обобщенного решения задачи (1.2).

Определение 1. Функцию $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];{{V}_{0}}),$ удовлетворяющую равенству при некотором $T > 0$

(6.1)
$\mathop {\left\langle {R(u),w} \right\rangle }\nolimits_0 = 0\quad {\text{д л я }}\;{\text{в с е х }}\quad t \in [0,T],\quad {\text{и }}\;{\text{в с е х }}\quad w \in {{V}_{0}},$
где
$R(u) \equiv \frac{d}{{dt}}\left( {{{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u)} \right) + {{L}_{1}}u + \int\limits_0^t \,dsh(t - s){{L}_{2}}u(s) + DP(u) = F(u),\quad u(0) = {{u}_{0}} \in {{V}_{0}},$
назовем сильным обобщенным решением задачи Коши (1.2).

Прежде всего нам нужно доказать, что оператор

$A(u) \equiv {{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u):{{V}_{0}} \to V_{0}^{ * }$
является обратимым, причем обратный оператор является липшиц-непрерывным. С этой целью докажем, что оператор $A(u)$ удовлетворяет всем условиям теоремы Браудера–Минти:

(I) оператор $A(u)$ является радиально-непрерывным.

Это утверждение есть следствие непрерывности операторов ${{A}_{0}}$ и ${{A}_{j}}( \cdot )$.

(II) оператор $A(u)$ является сильно монотонным.

Действительно, имеет место следующая цепочка неравенств:

$\begin{gathered} {{\left\langle {A({{u}_{1}}) - A({{u}_{2}}),{{u}_{1}} - {{u}_{2}}} \right\rangle }_{0}} = {{\left\langle {{{A}_{0}}{{u}_{1}} - {{A}_{0}}{{u}_{2}},{{u}_{1}} - {{u}_{2}}} \right\rangle }_{0}} + \sum\limits_{j = 1}^n {{\left\langle {{{A}_{j}}({{u}_{1}}) - {{A}_{j}}({{u}_{2}}),{{u}_{1}} - {{u}_{2}}} \right\rangle }_{j}} \geqslant \\ \geqslant \;{{\left\langle {{{A}_{0}}{{u}_{1}} - {{A}_{0}}{{u}_{2}},{{u}_{1}} - {{u}_{2}}} \right\rangle }_{0}}\; \geqslant \;{{m}_{0}}\left\| {{{u}_{1}} - {{u}_{2}}} \right\|_{0}^{2}. \\ \end{gathered} $

(III) оператор $A(u)$ является коэрцитивным.

Действительно, имеет место следующая цепочка неравенств:

${{\left\langle {A(u),u} \right\rangle }_{0}} = {{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}} + \sum\limits_{j = 1}^n {{\left\langle {{{A}_{j}}(u),u} \right\rangle }_{j}} \geqslant {{m}_{0}}\left\| u \right\|_{0}^{2} + \sum\limits_{j = 1}^n \,{{m}_{j}}\left\| u \right\|_{j}^{{{{p}_{j}}}} \geqslant {{m}_{0}}\left\| u \right\|_{0}^{2}.$
Таким образом, в силу теоремы Браудера–Минти для оператора
$A(u):{{V}_{0}} \to V_{0}^{ * }$
определен обратный оператор
${{A}^{{ - 1}}}(v):V_{0}^{ * } \to {{V}_{0}}.$
Докажем, что оператор ${{A}^{{ - 1}}}({\kern 1pt} \cdot {\kern 1pt} )$ является липшиц-непрерывным. Действительно, в силу сильной монотонности оператора $A(u)$ имеет место следующая цепочка неравенств:
$\begin{gathered} {{m}_{0}}\left\| {{{u}_{1}} - {{u}_{2}}} \right\|_{0}^{2}{{\left\langle {A({{u}_{1}}) - A({{u}_{2}}),{{u}_{1}} - {{u}_{2}}} \right\rangle }_{0}} \leqslant \\ \leqslant \;\left\| {A({{u}_{1}}) - A({{u}_{2}})} \right\|_{0}^{*}{{\left\| {{{u}_{1}} - {{u}_{2}}} \right\|}_{0}} \Rightarrow {{m}_{0}}{{\left\| {{{u}_{1}} - {{u}_{2}}} \right\|}_{0}} \leqslant \left\| {A({{u}_{1}}) - A({{u}_{2}})} \right\|_{0}^{ * }, \\ \end{gathered} $
откуда вытекает искомое неравенство
${{\left\| {{{A}^{{ - 1}}}({{w}_{1}}) - {{A}^{{ - 1}}}({{w}_{2}})} \right\|}_{0}} \leqslant \frac{1}{{{{m}_{0}}}}\left\| {{{w}_{1}} - {{w}_{2}}} \right\|_{0}^{ * }\quad {\text{д л я }}\;{\text{в с е х }}\quad {{w}_{1}},{{w}_{2}} \in V_{0}^{ * }.$
Таким образом, если ввести следующее обозначение:
$A(u) = v,$
то абстрактную задачу Коши (1.2) можно переписать в классе $v(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];V_{0}^{ * })$ в эквивалентном виде:

(6.2)
$\frac{{dv}}{{dt}} = - {{L}_{1}}{{A}^{{ - 1}}}(v) - \int\limits_0^t \,dsh(t - s){{L}_{2}}{{A}^{{ - 1}}}(v)(s) - DP({{A}^{{ - 1}}}(v)) + F({{A}^{{ - 1}}}(v)),$
(6.3)
$v(0) = {{v}_{0}} = A({{u}_{0}}) \in V_{0}^{ * }.$

В проверке корректности формулировки нуждается только начальное условие (6.3). Прежде всего заметим, что в классе $v(t) \in \mathbb{C}([0,T];V_{0}^{ * })$ имеет место следующая цепочка неравенств:

$u({{t}_{1}}) - u{{({{t}_{2}})}_{0}} = {{A}^{{ - 1}}}(v({{t}_{1}})) - {{A}^{{ - 1}}}{{(v({{t}_{2}}))}_{0}} \leqslant v({{t}_{1}}) - v({{t}_{2}})_{0}^{ * } \to 0\quad {\text{п р и }}\quad {{t}_{1}} \to {{t}_{2}}.$
Значит, $u(t) \in \mathbb{C}([0,T];{{V}_{0}})$. В силу непрерывности операторов ${{A}_{0}}$ и ${{A}_{j}}$ приходим к выводу о справедливости следующей цепочки выражений:
${{v}_{0}} = v(0) = \mathop {lim}\limits_{t \downarrow 0} v(t) = \mathop {lim}\limits_{t \downarrow 0} A(u(t)) = A(u(0)) = A({{u}_{0}}).$
Следовательно, начальное условие (6.3) корректно.

Заметим, что в классе $v(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];V_{0}^{ * })$ задача (6.2) эквивалентна интегральному уравнению:

(6.4)
$v(t) = {{v}_{0}} + \int\limits_0^t \,dsG(v)(s),$
где
$G(v)(t) = - {{L}_{1}}{{A}^{{ - 1}}}(v) - \int\limits_0^t \,dsh(t - s){{L}_{2}}{{A}^{{ - 1}}}(v)(s) - DP({{A}^{{ - 1}}}(v)) + F({{A}^{{ - 1}}}(v)).$
Будем искать решение интегрального уравнения (6.4) в классе $v(t) \in C(0,T;V_{0}^{ * })$. С этой целью перепишем интегральное уравнение (6.4) в операторном виде
(6.5)
$v(t) = H(v)(t),$
где
$H(v)(t) = {{v}_{0}} + \int\limits_0^t \,dsG(v)(s).$
Здесь можно воспользоваться результатами работы [26] и доказать следующее утверждение.

Теорема 1. Для любого ${{v}_{0}} \in V_{0}^{ * }$ найдется такое ${{T}_{0}} = {{T}_{0}}({{v}_{0}}) > 0,$ что существует единственное решение $v(t) \in \mathbb{C}([0,{{T}_{0}});V_{0}^{ * })$ уравнения (6.5), причем либо ${{T}_{0}} = + \infty $, либо ${{T}_{0}} < + \infty $ и в последнем случае справедливо следующее предельное свойство:

$\mathop {lim}\limits_{t \uparrow {{T}_{0}}} v_{0}^{ * }(t) = + \infty .$
Отметим, что
$G:\mathbb{C}([0,{{T}_{0}});V_{0}^{ * }) \to \mathbb{C}([0,{{T}_{0}});V_{0}^{ * })$
и поэтому
$\int\limits_0^t \,G(v)(s)ds \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});V_{0}^{ * }).$
Стало быть, решение $v(t)$ уравнения (6.5) принадлежит к классу ${{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});V_{0}^{*})$.

Итак, мы пришли к уравнению с известной правой частью:

(6.6)
${{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u) = v(t) \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});V_{0}^{*})\quad {\text{п р и }}\;{\text{н е к о т о р о м }}\quad {{T}_{0}} > 0.$
Нам понадобится следующая теорема (см. теорему 12.3.3 с. 651 из [27] или теорему 4.2.1 с. 60 из [28]):

Теорема 2. Пусть $P$непрерывно дифференцируемое отображение открытого шара $U = {{B}_{r}}({{x}_{0}})$ в банаховом пространстве $X$ в банахово пространство $Y$. Предположим, что оператор $\Lambda : = P_{f}^{'}({{x}_{0}})$ взаимно однозначно отображает $X$ на все $Y$. Тогда $P$ взаимно однозначно отображает некоторую окрестность $V$ точки ${{x}_{0}}$ на некоторую окрестность $W$ точки $P({{x}_{0}})$, причем отображение $R: = {{P}^{{ - 1}}}:W \to V$ непрерывно дифференцируемо и выполнено равенство

$R_{f}^{'}(y) = \mathop {(P_{f}^{'}({{P}^{{ - 1}}}(y)))}\nolimits^{ - 1} ,\quad y \in W.$
Пусть оператор $P(u): = {{A}_{0}}u + \sum\nolimits_{j = 1}^n {{{A}_{j}}(u)} ,$ $U = X = {{V}_{0}}$, $Y = V_{0}^{*}$. Рассмотрим производную Фреше оператора $P:$
$P_{f}^{'}({{u}_{0}}) = {{A}_{0}} + \sum\limits_{j = 1}^n \,A_{j}^{'}({{u}_{0}}):{{V}_{0}} \to V_{0}^{*},\quad {{u}_{0}} \in {{V}_{0}}.$
Докажем, что для любого ${{u}_{0}} \in {{V}_{0}}$ оператор $P_{f}^{'}({{u}_{0}})$ имеет обратный. С этой целью снова воспользуемся монотонностью рассматриваемых операторов.

(I) оператор $P_{f}^{'}({{u}_{0}})$ является радиально непрерывным.

Это следствие непрерывности оператора ${{A}_{0}}$ и того, что при фиксированном ${{u}_{0}} \in {{V}_{0}}$ оператор $A_{{if}}^{{''}}({{u}_{0}}) \in \mathcal{L}({{V}_{0}};V_{0}^{*})$.

(II) оператор $P_{f}^{'}({{u}_{0}})$ является сильно монотонным.

Действительно, имеет место следующая цепочка неравенств:

$\begin{gathered} \mathop {\left\langle {\left[ {{{A}_{0}} + \sum\limits_{j = 1}^n \,A_{{if}}^{'}({{u}_{0}})} \right]{{v}_{1}} - \left[ {{{A}_{0}} + \sum\limits_{j = 1}^n \,A_{{jf}}^{'}({{u}_{0}})} \right]{{v}_{2}},{{v}_{1}} - {{v}_{2}}} \right\rangle }\nolimits_0 = \mathop {\left\langle {{{A}_{0}}{{v}_{1}} - {{A}_{0}}{{v}_{2}},{{v}_{1}} - {{v}_{2}}} \right\rangle }\nolimits_0 + \\ + \;\sum\limits_{j = 1}^n \mathop {\,\left\langle {A_{{jf}}^{'}({{u}_{0}}){{v}_{1}} - A_{{jf}}^{'}({{u}_{0}}){{v}_{2}},{{v}_{1}} - {{v}_{2}}} \right\rangle }\nolimits_j \geqslant {{m}_{0}}\left\| {{{v}_{1}} - {{v}_{2}}} \right\|_{0}^{2}. \\ \end{gathered} $

(III) Оператор $P_{f}^{'}({{u}_{0}})$ является коэрцитивным.

Но это утверждение есть следствие пункта (II) и линейности этого оператора при фиксированной функции $u(t) \in \mathbb{C}([0,{{{\text{T}}}_{0}});{{V}_{0}}).$ Итак, оператор

$P_{f}^{'}({{u}_{0}}): = {{A}_{0}} + \sum\limits_{j = 1}^n \,A_{{jf}}^{'}({{u}_{0}}):{{V}_{0}} \to V_{0}^{*}$
является обратимым. Следовательно, в силу теоремы 2 из (6.6) имеет место равенство
$u(t) = R(v(t)) \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});{{V}_{0}}).$
Таким образом, из теорем 1, 2 и (6.6) вытекает

Теорема 3. Для любых функций ${{u}_{0}} \in {{V}_{0}}$ найдется такое ${{T}_{0}} = {{T}_{0}}({{u}_{0}}) > 0,$ что для любого $T \in (0,{{T}_{0}})$ существует единственное классическое решение задачи (1.2) класса $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,T];{{V}_{0}}),$ причем либо ${{T}_{0}} = + \infty $, либо ${{T}_{0}} < + \infty $ и в последнем случае имеем

$\mathop {lim}\limits_{t \uparrow {{T}_{0}}} {{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u)_{0}^{*}(t) = + \infty .$

7. РАЗРУШЕНИЕ СИЛЬНОГО ОБОБЩЕННОГО РЕШЕНИЯ ЗАДАЧИ (1.2) ПРИ ${{q}_{2}} + 2 > \bar {p}$

Пусть $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});{{V}_{0}})$ – сильное обобщенное решение задачи (1.2). Прежде всего введем следующие обозначения:

(7.1)
$\Phi (t) = \frac{1}{2}{{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}} + \sum\limits_{j = 1}^n \frac{{{{p}_{j}} - 1}}{{{{p}_{j}}}}{{\left\langle {{{A}_{j}}(u),u} \right\rangle }_{j}},$
(7.2)
$J(t) = {{\left\langle {{{A}_{0}}u{\text{'}},u{\text{'}}} \right\rangle }_{0}} + \sum\limits_{j = 1}^n {{\left\langle {A_{{jf}}^{'}(u)u',u{\text{'}}} \right\rangle }_{j}}.$

Справедлива

Лемма 6. Имеет место следующее неравенство:

(7.3)
${{(\Phi {\text{'}}(t))}^{2}} \leqslant \bar {p}J(t)\Phi (t)\quad п р и \quad \bar {p} = \mathop {max}\limits_{j = \overline {1,n} } {{p}_{j}},\quad t \in [0,{{T}_{0}}).$

Доказательство. В силу условий $A$ и ${{A}_{0}}$ справедливо неравенство Шварца для производных Фреше $A_{{j,u}}^{'}:{{V}_{j}} \to \mathcal{L}({{V}_{j}};V_{j}^{ * })$ операторов ${{A}_{j}}:{{V}_{j}} \to V_{j}^{ * }:$

(7.4)
$\begin{gathered} \left| {{{{\left\langle {\left( {{{A}_{j}}(u)} \right){\kern 1pt} ',u} \right\rangle }}_{j}}} \right| = \left| {{{{\left\langle {A_{{jf}}^{'}(u)u{\kern 1pt} ',u} \right\rangle }}_{j}}} \right| \leqslant \left\langle {A_{{jf}}^{'}(u)u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle _{j}^{{1/2}}\left\langle {A_{{jf}}^{'}(u)u,u} \right\rangle _{j}^{{1/2}} = \\ = \left\langle {\left( {{{A}_{j}}(u)} \right){\kern 1pt} ',u{\kern 1pt} '} \right\rangle _{j}^{{1/2}}{{({{p}_{j}} - 1)}^{{1/2}}}\left\langle {{{A}_{j}}(u),u} \right\rangle _{j}^{{1/2}}, \\ \left| {{{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u} \right\rangle }}_{0}}} \right| \leqslant \left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle _{0}^{{1/2}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}}. \\ \end{gathered} $
Здесь мы воспользовались следующим равенством:
(7.5)
$A_{{jf}}^{'}(v)v = ({{p}_{j}} - 1){{A}_{j}}(v),$
вывод которого приведен в лемме 1. Из (7.4), (7.5) вытекает
$\mathop {\left| {\frac{d}{{dt}}\Phi } \right|}\nolimits^2 \leqslant \mathop {\left| {\left| {{{{\left\langle {{{A}_{0}}u',u} \right\rangle }}_{0}}} \right| + \sum\limits_{j = 1}^N \left| {{{{\left\langle {\left( {{{A}_{j}}(u)} \right){\kern 1pt} ',u} \right\rangle }}_{j}}} \right|} \right|}\nolimits^2 \leqslant \left( {{{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N {{{\left\langle {\left( {{{A}_{j}}(u)} \right){\kern 1pt} ',u{\kern 1pt} '} \right\rangle }}_{j}}} \right)\; \times $
$ \times \;\left( {{{{\left\langle {{{A}_{0}}u,u} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N \,({{p}_{j}} - 1){{{\left\langle {{{A}_{j}}(u),u} \right\rangle }}_{j}}} \right) \leqslant \bar {p}\left( {{{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N {{{\left\langle {\left( {{{A}_{j}}({{u}_{m}})} \right){\kern 1pt} ',u_{m}^{'}} \right\rangle }}_{j}}} \right)\; \times $
$ \times \;\left( {\frac{1}{{\bar {p}}}{{{\left\langle {{{A}_{0}}u,u} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N \frac{{{{p}_{j}} - 1}}{{\bar {p}}}{{{\left\langle {{{A}_{j}}(u),u} \right\rangle }}_{j}}} \right) \leqslant \bar {p}\left( {{{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N {{{\left\langle {\left( {{{A}_{j}}(u)} \right){\kern 1pt} ',u{\kern 1pt} '} \right\rangle }}_{j}}} \right)\; \times $
$ \times \;\left( {\frac{1}{2}{{{\left\langle {{{A}_{0}}u,u} \right\rangle }}_{0}} + \sum\limits_{j = 1}^N \frac{{{{p}_{j}} - 1}}{{{{p}_{j}}}}{{{\left\langle {{{A}_{j}}(u),u} \right\rangle }}_{j}}} \right) = \bar {p}J(t)\Phi (t),$
где $\bar {p} = \mathop {max}\limits_{j = 1,N} {{p}_{j}} > 2$.

Лемма доказана.

Приступим к выводу первого энергетического равенства. Действительно, возьмем в равенстве (6.1) $w = u,$ тогда с учетом обозначения (7.1) получим первое энергетическое равенство:

(7.6)
$\frac{{d\Phi }}{{dt}} + {{({{L}_{1}}u,u)}_{1}} + \int\limits_0^t \,dsh(t - s){{({{L}_{2}}u(s),u(t))}_{2}} = {{(F(u),u)}_{3}},$
здесь мы воспользовались формулой (5.3). Теперь возьмем в равенстве (6.1) $w = u'$ и с учетом обозначения (7.2) получим второе энергетическое равенство:
$J(t) + \frac{1}{2}\frac{d}{{dt}}{{({{L}_{1}}u,u)}_{1}} + \int\limits_0^t \,dsh(t - s){{({{L}_{2}}u(s),u{\kern 1pt} '(t))}_{2}} = - {{\left\langle {DP(u),u{\kern 1pt} '} \right\rangle }_{0}} + \frac{1}{{{{q}_{2}} + 2}}\frac{d}{{dt}}{{(F(u),u)}_{3}},$
где мы воспользовались следующими равенствами:
${{(F(u),u{\kern 1pt} ')}_{3}} = \frac{1}{{{{q}_{2}} + 2}}\frac{d}{{dt}}{{(F(u),u)}_{3}},\quad {{({{L}_{1}}u,u{\kern 1pt} ')}_{1}} = \frac{1}{2}\frac{d}{{dt}}{{({{L}_{1}}u,u)}_{1}}$
для всех $u(t) \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});{{V}_{0}}),$ которые есть следствия леммы 1. Теперь подставим выражение для ${{(F(u),u)}_{3}}$ из первого энергетического равенства во второе энергетическое равенство. Тогда получим следующее равенство:
(7.7)
$\begin{gathered} J(t) = - \frac{{{{q}_{2}}}}{{2({{q}_{2}} + 2)}}\frac{d}{{dt}}{{({{L}_{1}}u,u)}_{1}} - \frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}\int\limits_0^t \,dsh(t - s){{({{L}_{2}}u(s),u'(t))}_{2}} + {{\left\langle {DP(u),u{\kern 1pt} '} \right\rangle }_{0}} + \\ + \;\frac{1}{{{{q}_{2}} + 2}}\Phi {\kern 1pt} {\text{''}} + \frac{{h(0)}}{{{{q}_{2}} + 2}}{{({{L}_{2}}u,u)}_{2}} + \frac{1}{{{{q}_{2}} + 2}}\int\limits_0^t \,dsh{\kern 1pt} '(t - s){{({{L}_{2}}u(s),u(t))}_{2}}. \\ \end{gathered} $
Теперь займемся постепенным оцениванием слагаемых в правой части равенства (7.7). В силу формулы (5.5) при $i = 1$ имеет место первая оценка:
(7.8)
$\begin{gathered} \left| {\frac{{{{q}_{2}}}}{{2({{q}_{2}} + 2)}}\frac{d}{{dt}}{{{({{L}_{1}}u,u)}}_{1}}} \right| \leqslant \frac{{{{q}_{2}}}}{{({{q}_{2}} + 2)}}\left| {{{{({{L}_{1}}u,u{\kern 1pt} ')}}_{1}}} \right| \leqslant \frac{\varepsilon }{2}{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle }_{0}} + \\ + \;\frac{{l_{1}^{2}}}{{2\varepsilon }}\mathop {\left( {\frac{{{{q}_{2}}}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 {{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}} \leqslant \frac{\varepsilon }{2}J(t) + \frac{{l_{1}^{2}}}{\varepsilon }\mathop {\left( {\frac{{{{q}_{2}}}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 \Phi (t). \\ \end{gathered} $
Теперь в силу формулы (5.5) при $i = 2$ имеет место следующая оценка:
$\begin{gathered} \left| {\frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}\int\limits_0^t dsh(t - s){{{({{L}_{2}}u(s),u'(t))}}_{2}}} \right| \leqslant \frac{\varepsilon }{2}{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u'{\kern 1pt} } \right\rangle }_{0}} + \mathop {\left( {\frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 \frac{{l_{2}^{2}a}}{{2\varepsilon }}\int\limits_0^t \,ds\left| {h(t - s)} \right|{{\left\langle {{{A}_{0}}u(s),u(s)} \right\rangle }_{0}} \leqslant \\ \leqslant \;\frac{\varepsilon }{2}J(t) + \mathop {\left( {\frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 \frac{{l_{2}^{2}a}}{\varepsilon }\int\limits_0^t \,ds\left| {h(t - s)} \right|\Phi (s),\quad a: = \int\limits_0^{ + \infty } \left| {h(t)} \right|dt. \\ \end{gathered} $
Теперь в силу формулы (5.4) имеет место следующее неравенство:
$\left| {{{{\left\langle {DP(u),u{\kern 1pt} '} \right\rangle }}_{0}}} \right| \leqslant \frac{\varepsilon }{2}{{\left\langle {{{A}_{0}}u{\kern 1pt} ',u{\kern 1pt} '} \right\rangle }_{0}} + \frac{{c_{1}^{2}}}{{2\varepsilon }}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1 + {{q}_{3}}}} \leqslant \frac{\varepsilon }{2}J(t) + \frac{{c_{1}^{2}{{2}^{{{{q}_{3}}}}}}}{\varepsilon }{{\Phi }^{{1 + {{q}_{3}}}}}(t).$
Теперь в силу формулы (5.5) при $i = 2$ имеет место следующая оценка:
$\left| {\frac{{h(0)}}{{{{q}_{2}} + 2}}{{{({{L}_{2}}u,u)}}_{2}}} \right| \leqslant \frac{{{{l}_{2}}\left| {h(0)} \right|}}{{{{q}_{2}} + 2}}{{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}} \leqslant \frac{{2{{l}_{2}}\left| {h(0)} \right|}}{{{{q}_{2}} + 2}}\Phi (t).$
Опять в силу формулы (5.5) при $i = 2$ имеет место следующая оценка:
(7.9)
$\begin{gathered} \left| {\frac{1}{{{{q}_{2}} + 2}}\int\limits_0^t \,dsh{\kern 1pt} '(t - s){{{({{L}_{2}}u(s),u(t))}}_{2}}} \right| \leqslant \frac{{b{{l}_{2}}}}{{2({{q}_{2}} + 2)}}{{\left\langle {{{A}_{0}}u,u} \right\rangle }_{0}} + \frac{{{{l}_{2}}}}{{2({{q}_{2}} + 2)}}\int\limits_0^t \,ds\left| {h{\kern 1pt} '(t - s)} \right|{{\left\langle {{{A}_{0}}u(s),u(s)} \right\rangle }_{0}} \leqslant \\ \leqslant \;\frac{{b{{l}_{2}}}}{{{{q}_{2}} + 2}}\Phi (t) + \frac{{{{l}_{2}}}}{{{{q}_{2}} + 2}}\int\limits_0^t \,ds\left| {h{\kern 1pt} '(t - s)} \right|\Phi (s),\quad b: = \int\limits_0^{ + \infty } \left| {h{\kern 1pt} '(t)} \right|dt. \\ \end{gathered} $
Таким образом, из оценок (7.8), (7.9) и из равенства (7.7) вытекает следующее неравенство:
$\begin{gathered} \left( {1 - \frac{3}{2}\varepsilon } \right)J(t) \leqslant \frac{1}{{{{q}_{2}} + 2}}\Phi {\kern 1pt} {\text{''}} + \left[ {\frac{{l_{1}^{2}}}{\varepsilon }\mathop {\left( {\frac{{{{q}_{2}}}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 + \frac{{2{{l}_{2}}\left| {h(0)} \right|}}{{{{q}_{2}} + 2}} + \frac{{b{{l}_{2}}}}{{{{q}_{2}} + 2}}} \right]\Phi (t) + \\ + \;\frac{{c_{1}^{2}{{2}^{{{{q}_{3}}}}}}}{\varepsilon }{{\Phi }^{{1 + {{q}_{3}}}}}(t) + \int\limits_0^t \,ds{{h}_{\varepsilon }}(t - s)\Phi (s), \\ \end{gathered} $
где
${{h}_{\varepsilon }}(t) = \mathop {\left( {\frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 \frac{{l_{2}^{2}a}}{\varepsilon }{\text{|}}h(t){\text{|}} + \frac{{{{l}_{2}}}}{{{{q}_{2}} + 2}}{\text{|}}h{\kern 1pt} '(t){\text{|}}.$
Потребуем пока, чтобы
$\varepsilon \in \left( {0,\;\frac{2}{3}} \right),$
тогда из неравенств (7.3) приходим к искомому интегродифференциальному неравенству:
$\Phi \Phi {\text{''}} - \alpha {{(\Phi {\kern 1pt} ')}^{2}} + {{\gamma }_{1}}{{\Phi }^{2}} + {{\gamma }_{3}}{{\Phi }^{{1 + \lambda }}} + {{\gamma }_{2}}\int\limits_0^t \,ds{{h}_{\varepsilon }}(t - s)\Phi (s)\Phi (t) \geqslant 0,$
где
$\alpha = \frac{{{{q}_{2}} + 2}}{{\bar {p}}}\left( {1 - \frac{3}{2}\varepsilon } \right),\quad \lambda = 1 + {{q}_{3}},$
${{\gamma }_{1}} = ({{q}_{2}} + 2)\left[ {\frac{{l_{1}^{2}}}{\varepsilon }\mathop {\left( {\frac{{{{q}_{2}}}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 + \frac{{2{{l}_{2}}\left| {h(0)} \right|}}{{{{q}_{2}} + 2}} + \frac{{b{{l}_{2}}}}{{{{q}_{2}} + 2}}} \right],$
${{\gamma }_{3}} = ({{q}_{2}} + 2)\frac{{c_{1}^{2}{{2}^{{{{q}_{3}}}}}}}{\varepsilon },\quad {{\gamma }_{2}} = {{q}_{2}} + 2.$
Теперь проверим выполнимость условия
$\lambda < 2\alpha - 1.$
Действительно, имеет место следующая цепочка рассуждений:
$\begin{gathered} \lambda < 2\alpha - 1 \Rightarrow 2 + {{q}_{3}} < 2\alpha \Rightarrow \alpha > 1 + {{q}_{3}}{\text{/}}2 \Rightarrow 1 - \frac{3}{2}\varepsilon > \left( {1 + \frac{{{{q}_{3}}}}{2}} \right)\frac{{\bar {p}}}{{{{q}_{2}} + 2}} \Rightarrow \\ \Rightarrow \varepsilon \in \left( {0,\frac{2}{3}\left[ {1 - \left( {1 + \frac{{{{q}_{3}}}}{2}} \right)\frac{{\bar {p}}}{{{{q}_{2}} + 2}}} \right]} \right) \Rightarrow {{q}_{3}} < 2\frac{{{{q}_{2}} + 2 - \bar {p}}}{{\bar {p}}}. \\ \end{gathered} $
Итак, приходим к двум условиям:
(7.10)
${{q}_{2}} + 2 > \bar {p},\quad {{q}_{3}} < 2\frac{{{{q}_{2}} + 2 - \bar {p}}}{{\bar {p}}}.$
Таким образом, в силу результата леммы 8 приходим к следующему утверждению.

Теорема 4. Пусть ${{u}_{0}} \in {{V}_{0}}$ и выполнены условия (7.10), тогда при выполнении следующих начальных условий:

$\Phi {\kern 1pt} '(0) > \mathop {\left( {\frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} + \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{2 + {{q}_{3}}}}}} \right)}\nolimits^{1/2} ,\quad \Phi (0) > 0,$
существует единственное сильное обобщенное решение класса
$u(t) \in {{\mathbb{C}}^{{(1)}}}([0,{{T}_{0}});{{V}_{0}}),$
причем ${{T}_{0}} \in (0,{{T}_{1}}]$ и имеет место предельное выражение
$\mathop {\lim \sup }\limits_{t \uparrow {{T}_{0}}} \Phi (t) = + \infty ,$
где

$\Phi {\kern 1pt} '(0) = {{(F({{u}_{0}}),{{u}_{0}})}_{3}} - {{({{L}_{1}}{{u}_{0}},{{u}_{0}})}_{1}},$
$\Phi (0) = \frac{1}{2}{{\left\langle {{{A}_{0}}{{u}_{0}},{{u}_{0}}} \right\rangle }_{0}} + \sum\limits_{j = 1}^n \frac{{{{p}_{j}} - 1}}{{{{p}_{j}}}}{{\left\langle {{{A}_{j}}({{u}_{0}}),{{u}_{0}}} \right\rangle }_{j}},$
${{T}_{1}} = \frac{{{{\Phi }^{{1 - \alpha }}}(0)}}{A},$
${{A}^{2}} = {{(\alpha - 1)}^{2}}{{\Phi }^{{ - 2\alpha }}}(0)\left[ {{{{(\Phi {\kern 1pt} '(0))}}^{2}} - \frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} - \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{\delta (1 - \alpha ) + 2\alpha }}}} \right],$
${{\gamma }_{4}} = {{\gamma }_{1}} + {{\gamma }_{2}}\left| {\chi (0)} \right|\left[ {\mathop {\left( {\frac{{{{q}_{2}} + 1}}{{{{q}_{2}} + 2}}} \right)}\nolimits^2 \frac{{l_{2}^{2}a}}{\varepsilon } + \frac{{{{l}_{2}}}}{{{{q}_{2}} + 2}}} \right],\quad \delta = 1 + \frac{{1 + {{q}_{3}} - \alpha }}{{1 - \alpha }},$
$\varepsilon \in \left( {0,\;\frac{2}{3}\left[ {1 - \left( {1 + \frac{{{{q}_{3}}}}{2}} \right)\frac{{\bar {p}}}{{{{q}_{2}} + 2}}} \right]} \right).$

Справедлива

Лемма 7. Имеет место двустороннее неравенство

${{M}_{1}}{{\Phi }^{{1/2}}}(t) \leqslant \left\| {A(u)} \right\|_{0}^{*} \leqslant {{M}_{2}}{{\Phi }^{{1/2}}} + \sum\limits_{j = 1}^n \,{{B}_{j}}{{\Phi }^{{({{p}_{j}} - 1)/{{p}_{j}}}}}(t),$
где постоянные ${{M}_{1}},$ ${{M}_{2}}$ и ${{B}_{j}}$ больше нуля и не зависят от $u(t),$

$A(u): = {{A}_{0}}u + \sum\limits_{j = 1}^n \,{{A}_{j}}(u).$

Доказательство. Докажем сначала оценку снизу. Действительно, с одной стороны, имеем

(7.11)
${{\left\langle {A(u),u} \right\rangle }_{0}} = \left\langle {{{A}_{0}}u,u} \right\rangle + \sum\limits_{j = 1}^n {{\left\langle {{{A}_{j}}(u),u} \right\rangle }_{j}} \geqslant \Phi (t).$
С другой стороны, имеем
(7.12)
${{\left\langle {A(u),u} \right\rangle }_{0}} \leqslant \left\| {A(u)} \right\|_{0}^{*}{{\left\| u \right\|}_{0}} \leqslant \left\| {A(u)} \right\|_{0}^{ * }\frac{1}{{m_{0}^{{1/2}}}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}} \leqslant \left\| {A(u)} \right\|_{0}^{ * }\frac{{{{2}^{{1/2}}}}}{{m_{0}^{{1/2}}}}{{\Phi }^{{1/2}}}(t).$
Из неравенств (7.11) и (7.12) вытекает оценка снизу
$\left\| {A(u)} \right\|_{0}^{ * } \geqslant \mathop {\left( {\frac{{{{m}_{0}}}}{2}} \right)}\nolimits^{1/2} {{\Phi }^{{1/2}}}(t) = {{M}_{1}}{{\Phi }^{{1/2}}}(t).$
Докажем теперь оценку сверху. Действительно, согласно определению нормы $\left\| {\, \cdot \,} \right\|_{0}^{ * }$ справедлива следующая цепочка выражений:

$\begin{gathered} \left\| {A(u)} \right\|_{0}^{ * } = \mathop {sup}\limits_{{{{\left\| h \right\|}}_{0}} \leqslant 1} \left| {{{{\left\langle {A(u),h} \right\rangle }}_{0}}} \right| = \mathop {sup}\limits_{{{{\left\| h \right\|}}_{0}} \leqslant 1} \left| {{{{\left\langle {{{A}_{0}}u,h} \right\rangle }}_{0}} + \sum\limits_{j = 1}^n {{{\left\langle {{{A}_{j}}(u),h} \right\rangle }}_{j}}} \right| \leqslant \mathop {sup}\limits_{{{{\left\| h \right\|}}_{0}} \leqslant 1} \left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1/2}}\left\langle {{{A}_{0}}h,h} \right\rangle _{0}^{{1/2}} + \\ + \;\mathop {sup}\limits_{{{{\left\| h \right\|}}_{0}} \leqslant 1} \sum\limits_{j = 1}^n \left\| {{{A}_{j}}(u)} \right\|_{j}^{ * }{{\left\| h \right\|}_{j}} \leqslant {{c}_{3}}{{\Phi }^{{1/2}}}(t) + {{c}_{4}}\sum\limits_{j = 1}^n \left\| u \right\|_{j}^{{{{p}_{j}} - 1}} \leqslant {{c}_{3}}{{\Phi }^{{1/2}}}(t) + {{c}_{5}}\sum\limits_{j = 1}^n {{\left\langle {{{A}_{j}}(u),u} \right\rangle }^{{({{p}_{j}} - 1)/{{p}_{j}}}}} \leqslant \\ \leqslant \;{{c}_{3}}{{\Phi }^{{1/2}}}(t) + {{c}_{6}}\sum\limits_{j = 1}^n \,{{\Phi }^{{({{p}_{j}} - 1)/{{p}_{j}}}}}(t) = {{M}_{2}}{{\Phi }^{{1/2}}} + \sum\limits_{j = 1}^n \,{{B}_{j}}{{\Phi }^{{({{p}_{j}} - 1)/{{p}_{j}}}}}(t). \\ \end{gathered} $

Лемма доказана.

Замечание 1. В силу результата этой леммы имеем ${{T}_{0}} < + \infty $ при выполнении условий теоремы 4.

8. ОЦЕНКА СНИЗУ НА ВРЕМЯ СУЩЕСТВОВАНИЯ СИЛЬНОГО РЕШЕНИЯ ЗАДАЧИ (1.2)

Теперь приступим к получению нижней оценки на время ${{T}_{0}} > 0$ существования решения задачи (1.2).

Напомним, что величина ${{({{L}_{1}}u,u)}_{1}} \geqslant 0$ для всех $u \in {{V}_{0}},$ поэтому из первого энергетического равенства (7.6) вытекает следующее неравенство:

(8.1)
$\frac{{d\Phi }}{{dt}} \leqslant - \int\limits_0^t \,dsh(t - s){{({{L}_{2}}u(s),u(t))}_{2}} + {{(F(u),u)}_{3}}.$
Справедлива следующая цепочка неравенств:
(8.2)
$\left| {{{{(F(u),u)}}_{3}}} \right| \leqslant {{M}_{F}}\left| u \right|_{3}^{{{{q}_{2}} + 2}} \leqslant {{M}_{F}}{{c}_{1}}\left\| u \right\|_{0}^{{{{q}_{2}} + 2}} \leqslant {{M}_{F}}{{c}_{1}}m_{0}^{{1 + {{q}_{2}}/2}}\left\langle {{{A}_{0}}u,u} \right\rangle _{0}^{{1 + {{q}_{2}}/2}} \leqslant {{c}_{2}}{{\Phi }^{{1 + {{q}_{2}}/2}}}(t).$
Как и в случае с оценкой (7.9) можно получить следующее неравенство:
(8.3)
$\left| {\int\limits_0^t \,dsh(t - s){{{({{L}_{2}}u(s),u(t))}}_{2}}} \right| \leqslant {{l}_{2}}\Phi (t) + {{l}_{2}}a\int\limits_0^t \,ds{\text{|}}h(t - s){\text{|}}\Phi (s).$
Таким образом, из (8.2) и (8.3) и неравенства (8.1) вытекает оценка:
(8.4)
$\frac{{d\Phi }}{{dt}} \leqslant {{l}_{2}}\Phi (t) + {{l}_{2}}{{h}_{0}}\int\limits_0^t \,ds\Phi (s) + {{c}_{2}}{{\Phi }^{{1 + {{q}_{2}}/2}}}(t).$
Теперь воспользуемся арифметическим неравенством Гёльдера:
$r\Phi \leqslant {{r}^{{(2 + {{q}_{2}})/{{q}_{2}}}}}\frac{{{{q}_{2}}}}{{{{q}_{2}} + 2}} + \frac{2}{{2 + {{q}_{2}}}}{{\Phi }^{{1 + {{q}_{2}}/2}}}.$
Тогда из (8.4) после интегрирования по времени $t \in [0,T],$ получим неравенство вида
$\Phi (t) \leqslant {{c}_{4}}(T) + {{c}_{5}}(T)\int\limits_0^t \,ds{{\Phi }^{{1 + {{q}_{2}}/2}}}(s),$
из которого в силу леммы Гронуолла–Беллмана–Бихари (см. [29]) получим неравенство
$\Phi (t) \leqslant \frac{{{{c}_{3}}(T)}}{{\mathop {\left[ {1 - \tfrac{{{{q}_{2}}}}{2}c_{4}^{{{{q}_{2}}/2}}(T){{c}_{5}}(T)t} \right]}\nolimits^{2/{{q}_{2}}} }},$
откуда в свою очередь вытекает, что при выполнении условия
$t \in [0,{{T}_{2}}),\quad \frac{{{{q}_{2}}}}{2}c_{4}^{{{{q}_{2}}/2}}({{T}_{2}}){{c}_{5}}({{T}_{2}}){{T}_{2}} = 1$
имеет место оценка
$\Phi (t) \leqslant {{c}_{6}}({{T}_{2}}) < + \infty \quad {\text{д л я }}\;{\text{в с е х }}\quad t \in [0,{{T}_{2}}].$
Следовательно, приходим к неравенству

${{T}_{0}} \geqslant {{T}_{2}}.$

9. ГЛОБАЛЬНАЯ РАЗРЕШИМОСТЬ ЗАДАЧИ (1.2) ПРИ ${{q}_{2}} + 2 \leqslant \bar {p}$

Теперь мы приступим к получению одного достаточного условия глобального во времени существования сильного обобщенного решения задачи (1.2), которое, как мы увидим, близко к необходимому. Действительно, пусть выполнено следующее условие:

$\bar {p} \geqslant {{q}_{2}} + 2,$
причем предположим, что существует такое $j \in \overline {1,n} ,$ что ${{V}_{j}} \subset {{W}_{3}}$ и ${{p}_{j}} = \bar {p}$. Тогда справедлива следующая цепочка неравенств:
(9.1)
$\begin{gathered} \left| {{{{(F(u),u)}}_{3}}} \right| \leqslant {{M}_{F}}\left| u \right|_{3}^{{{{q}_{2}} + 2}} \leqslant {{M}_{F}}{{c}_{j}}\left\| u \right\|_{j}^{{{{q}_{2}} + 2}} = \frac{{{{M}_{F}}{{c}_{j}}}}{{m_{j}^{{({{q}_{2}} + 2)/{{p}_{j}}}}}}\mathop {({{m}_{j}}\left\| u \right\|_{j}^{{{{p}_{j}}}})}\nolimits^{({{q}_{2}} + 2)/{{p}_{j}}} \leqslant \\ \leqslant \;\frac{{{{M}_{F}}{{c}_{j}}}}{{m_{j}^{{({{q}_{2}} + 2)/{{p}_{j}}}}}}\left\langle {{{A}_{j}}(u),u} \right\rangle _{j}^{{({{q}_{2}} + 2)/{{p}_{j}}}} \leqslant m{{\Phi }^{{({{q}_{2}} + 2)/{{p}_{j}}}}}. \\ \end{gathered} $
Теперь в силу оценки (8.3) и неравенства (9.1) приходим к следующему неравенству:
$\frac{{d\Phi }}{{dt}} \leqslant m{{\Phi }^{{({{q}_{2}} + 2)/{{p}_{j}}}}} + {{l}_{2}}\Phi (t) + {{l}_{2}}{{h}_{0}}\int\limits_0^t \,ds\Phi (s).$
Здесь нужно рассмотреть два случая: (I) ${{q}_{2}} + 2 = \bar {p}$ и (II) ${{q}_{2}} + 2 < \bar {p}$. В первом случае мы приходим к следующему неравенству:
$\frac{{d\Phi }}{{dt}} \leqslant ({{l}_{2}} + m)\Phi (t) + {{l}_{2}}{{h}_{0}}\int\limits_0^t \,ds\Phi (s).$
Рассмотрим теперь более сложный второй случай. Здесь можно воспользоваться арифметическим неравенством Гёльдера. Действительно, положим
$\beta = \frac{{\bar {p}}}{{{{q}_{2}} + 2}} > 1,$
поэтому имеет место неравенство
$r{{\Phi }^{{1/\beta }}} \leqslant \frac{{\beta - 1}}{\beta }{{r}^{{\beta /(\beta - 1)}}} + \frac{1}{\beta }\Phi .$
И в результате и во втором случае мы приходим к неравенству вида
(9.2)
$\frac{{d\Phi }}{{dt}} \leqslant {{r}_{1}} + {{r}_{2}}\Phi (t) + {{r}_{3}}\int\limits_0^t \,ds\Phi (s),\quad {{r}_{1}},{{r}_{2}},{{r}_{3}} \geqslant 0.$
Теперь докажем, что решение дифференциального неравенства (9.2) существует глобально во времени. С этой целью проинтегрируем обе части неравенства (9.2) по времени $t \in [0,T]$, тогда справедлива следующая цепочка неравенств:
$\Phi (t) \leqslant {{r}_{1}}t + {{r}_{2}}\int\limits_0^t \,ds\Phi (s) + \int\limits_0^t \,ds(t - s)\Phi (s) \leqslant {{r}_{1}}T + [{{r}_{2}} + T{{r}_{3}}]\int\limits_0^t \,ds\Phi (s).$
Теперь воспользуемся леммой Гронуолла–Беллмана (см. [29]) и получим следующую оценку:
$\Phi (t) \leqslant ({{\Phi }_{0}} + {{r}_{1}}T)exp\{ [{{r}_{2}} + T{{r}_{3}}]t\} \quad {\text{д л я }}\;{\text{в с е х }}\quad t \in [0,T].$
Отсюда, из теоремы 3 и леммы 7 вытекает, что в данном случае ${{T}_{0}} = + \infty .$ Таким образом, нами доказана

Теорема 5. Пусть ${{u}_{0}} \in {{V}_{0}}$ и выполнено неравенство

$\bar {p} \geqslant {{q}_{2}} + 2,\quad \bar {p} = \mathop {max}\limits_{j = \overline {1,n} } {{p}_{j}},$
причем найдется такое $m \in \overline {1,n} ,$ что $\bar {p} = {{p}_{m}}$ и имеет место непрерывное вложение ${{V}_{m}} \subset {{W}_{3}},$ тогда существует единственное сильное обобщенное решение задачи (1.2) класса

$u(t) \in {{\mathbb{C}}^{{(1)}}}([0, + \infty );{{V}_{0}}).$

10. РЕШЕНИЕ ОСНОВНОГО ИНТЕГРОДИФФЕРЕНЦИАЛЬНОГО НЕРАВЕНСТВА

Рассмотрим теперь следующее интегродифференциальное неравенство:

(10.1)
$\Phi \Phi {\kern 1pt} {\text{''}} - \alpha {{(\Phi {\kern 1pt} ')}^{2}} + {{\gamma }_{1}}{{\Phi }^{2}} + {{\gamma }_{3}}{{\Phi }^{{1 + \lambda }}} + {{\gamma }_{2}}\int\limits_0^t \,dsh(t - s)\Phi (s)\Phi (t) \geqslant 0,$
где $\alpha > 1,$ ${{\gamma }_{1}},{{\gamma }_{2}},{{\gamma }_{3}},\lambda \geqslant 0$. Предположим, что найдется такая функция
$\chi (t) \in {{\mathbb{C}}^{{(1)}}}(\mathbb{R}_{ + }^{1}),$
что имеют место следующие неравенства:
(10.2)
$\left| {h(t)} \right| \leqslant \chi {\kern 1pt} '(t)\quad {\text{д л я }}\;{\text{в с е х }}\quad t \in \mathbb{R}_{ + }^{1},$
причем
(10.3)
$\chi (t) \leqslant 0\quad {\text{п р и }}\quad t \geqslant 0.$
С учетом неравенства (10.2) мы получим следующую цепочку соотношений:
$\begin{gathered} \int\limits_0^t \,dsh(t - s)\Phi (s) \leqslant \int\limits_0^t \,ds{\text{|}}h(t - s){\text{|}}\Phi (s) \leqslant \int\limits_0^t \,ds\chi {\kern 1pt} '(t - s)\Phi (s) = - \int\limits_0^t \,d\chi (t - s)\Phi (s) = \\ = - \chi (t - s)\left. {\Phi (s)} \right|_{{s = 0}}^{{s = t}} + \int\limits_0^t \,ds\chi (t - s)\Phi '(s) = - \chi (0)\Phi (t) + \chi (t)\Phi (0) + \int\limits_0^t \,ds\chi (t - s)\Phi {\kern 1pt} {\text{'}}(s). \\ \end{gathered} $
Отсюда получим следующее неравенство:
(10.4)
$\Phi \Phi {\text{''}} - \alpha {{(\Phi {\kern 1pt} ')}^{2}} + {{\gamma }_{4}}{{\Phi }^{2}} + {{\gamma }_{3}}{{\Phi }^{{1 + \lambda }}} \geqslant 0\quad {\text{п р и }}\quad t \in [0,{{t}_{1}}],\quad {{t}_{1}} \in (0,T],$
где
${{\gamma }_{4}} = {{\gamma }_{1}} + {{\gamma }_{2}}{\text{|}}\chi (0){\text{|}},$
причем в предположении $\Phi {\kern 1pt} '(0) > 0$ имеет место неравенство
$\Phi {\kern 1pt} '(t) \geqslant 0\quad {\text{п р и }}\quad t \in [0,{{t}_{1}}].$
Кроме того, потребуем, чтобы
$\Phi (0) > 0.$
Из неравенства (10.4) вытекает следующая цепочка неравенств:
$\begin{gathered} \quad \frac{{\Phi {\kern 1pt} ''}}{{{{\Phi }^{\alpha }}}} - \alpha \frac{{{{{(\Phi {\kern 1pt} ')}}^{2}}}}{{{{\Phi }^{{1 + \alpha }}}}} + {{\gamma }_{4}}{{\Phi }^{{1 - \alpha }}} + {{\gamma }_{3}}{{\Phi }^{{\lambda - \alpha }}} \geqslant 0 \Rightarrow \left[ {\frac{{\Phi {\kern 1pt} '}}{{{{\Phi }^{\alpha }}}}} \right]{\text{'}} + {{\gamma }_{4}}{{\Phi }^{{1 - \alpha }}} + {{\gamma }_{3}}{{\Phi }^{{\lambda - \alpha }}} \geqslant 0 \Rightarrow \\ \quad \Rightarrow \;\frac{1}{{1 - \alpha }}({{\Phi }^{{1 - \alpha }}}){\kern 1pt} {\text{''}} + {{\gamma }_{4}}{{\Phi }^{{1 - \alpha }}} + {{\gamma }_{3}}{{\Phi }^{{\lambda - \alpha }}} \geqslant 0 \Rightarrow Z(t) = {{\Phi }^{{1 - \alpha }}}, \\ \quad - \frac{1}{{\alpha - 1}}Z{\kern 1pt} {\text{''}} + {{\gamma }_{4}}Z + {{\gamma }_{3}}{{Z}^{{(\lambda - \alpha )/(1 - \alpha )}}} \geqslant 0 \Rightarrow Z{\kern 1pt} '' \leqslant (\alpha - 1){{\gamma }_{4}}Z + (\alpha - 1){{\gamma }_{3}}{{Z}^{{(\lambda - \alpha )/(1 - \alpha )}}}. \\ \end{gathered} $(10.5)
Заметим, что
$Z{\kern 1pt} ' = (1 - \alpha ){{\Phi }^{{ - \alpha }}}\Phi {\kern 1pt} ' \leqslant 0\quad {\text{п р и }}\quad t \in [0,{{t}_{1}}].$
Поэтому умножим обе части итогового неравенства из цепочки неравенств (10.5) и получим следующее неравенство:
$Z{\kern 1pt} 'Z{\kern 1pt} '' \geqslant (\alpha - 1){{\gamma }_{4}}ZZ{\kern 1pt} {\text{'}} + (\alpha - 1){{\gamma }_{3}}Z{\kern 1pt} {\text{'}}{{Z}^{{(\lambda - \alpha )/(1 - \alpha )}}}.$
Откуда приходим к выводу, что
(10.6)
$\frac{1}{2}\frac{d}{{dt}}\mathop {(Z{\kern 1pt} ')}\nolimits^2 \geqslant \frac{{(\alpha - 1){{\gamma }_{4}}}}{2}\frac{d}{{dt}}{{Z}^{2}} + \frac{{(\alpha - 1){{\gamma }_{3}}}}{\delta }\frac{d}{{dt}}{{Z}^{\delta }},$
где
$\delta = 1 + \frac{{\lambda - \alpha }}{{1 - \alpha }}.$
Теперь потребуем выполнения условия $\delta > 0.$ Займемся арифметикой:
$\delta > 0 \Rightarrow 1 + \frac{{\lambda - \alpha }}{{1 - \alpha }} > 0 \Rightarrow \frac{{1 + \lambda - 2\alpha }}{{1 - \alpha }} > 0 \Rightarrow \frac{{2\alpha - 1 - \lambda }}{{\alpha - 1}} > 0 \Rightarrow \lambda < 2\alpha - 1.$
Итак, потребуем выполнения следующего условия:
(10.7)
$\lambda < 2\alpha - 1.$
При выполнении условия (10.7) продолжим рассмотрение неравенства (10.6). Проинтегрируем его по $t \in [0,{{t}_{1}}]$ и получим, что
(10.8)
${{(Z{\kern 1pt} ')}^{2}} \geqslant (\alpha - 1){{\gamma }_{4}}{{Z}^{2}} + \frac{{2(\alpha - 1){{\gamma }_{3}}}}{\delta }{{Z}^{\delta }} + {{A}^{2}},$
где
${{A}^{2}} = {{(Z{\kern 1pt} '(0))}^{2}} - (\alpha - 1){{\gamma }_{4}}{{(Z(0))}^{2}} - \frac{{2(\alpha - 1){{\gamma }_{3}}}}{\delta }{{(Z(0))}^{\delta }},\quad \delta = 1 + \frac{{\lambda - \alpha }}{{1 - \alpha }}.$
Потребуем теперь выполнения условия ${{A}^{2}} > 0.$ Займемся арифметикой:
(10.9)
$\begin{gathered} {{A}^{2}} = {{(\alpha - 1)}^{2}}{{\Phi }^{{ - 2\alpha }}}(0){{(\Phi {\kern 1pt} '(0))}^{2}} - (\alpha - 1){{\gamma }_{4}}{{(\Phi (0))}^{{2 - 2\alpha }}} - \frac{{2(\alpha - 1){{\gamma }_{3}}}}{\delta }{{(\Phi (0))}^{{\delta - \delta \alpha }}} = \\ = {{(\alpha - 1)}^{2}}{{\Phi }^{{ - 2\alpha }}}(0)\left[ {{{{(\Phi {\kern 1pt} '(0))}}^{2}} - \frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} - \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{\delta (1 - \alpha ) + 2\alpha }}}} \right]. \\ \end{gathered} $
Заметим, что имеет место следующая цепочка равенств:
$\delta (1 - \alpha ) + 2\alpha = 1 - \alpha + \lambda - \alpha - 2\alpha = 1 + \lambda .$
С учетом этой цепочки приходим из (10.9) к эквивалентному условию ${{A}^{2}} > 0{\text{:}}$
(10.10)
$\Phi {\kern 1pt} '(0) > \mathop {\left( {\frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} + \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{1 + \lambda }}}} \right)}\nolimits^{1/2} .$
При выполнении последнего условия (10.10) из неравенства (10.8) получим, что
(10.11)
$\begin{gathered} {{(Z{\kern 1pt} ')}^{2}} \geqslant {{A}^{2}} \Rightarrow {\text{|}}Z{\kern 1pt} '{\text{|}} \geqslant A \Rightarrow - Z{\kern 1pt} ' \geqslant A \Rightarrow Z{\kern 1pt} ' \leqslant - A \Rightarrow \\ \Rightarrow \;(1 - \alpha ){{\Phi }^{{ - \alpha }}}\Phi {\kern 1pt} ' \leqslant - A \Rightarrow \Phi {\kern 1pt} ' \geqslant \frac{A}{{\alpha - 1}}{{\Phi }^{\alpha }}(t) > 0\quad {\text{п р и }}\quad t \in [0,{{t}_{1}}], \\ \end{gathered} $
поскольку
$\Phi (0) > 0,\quad \Phi {\kern 1pt} '(t) \geqslant 0\quad {\text{п р и }}\quad t \in [0,{{t}_{1}}].$
Таким образом, рассуждая дальше как и ранее, получаем, что
$\Phi {\kern 1pt} '(t) \geqslant 0\quad {\text{д л я }}\;{\text{в с е х }}\quad t \in [0,T],$
но тогда имеет место цепочка неравенств, вытекающая из (10.11):
$Z(t) \leqslant Z(0) - At \Rightarrow {{\Phi }^{{1 - \alpha }}}(t) \leqslant \Phi _{0}^{{1 - \alpha }} - At \Rightarrow {{\Phi }^{{\alpha - 1}}}(t) \geqslant \frac{1}{{\Phi _{0}^{{1 - \alpha }} - At}} \Rightarrow \Phi (t) \geqslant \frac{1}{{\mathop {[\Phi _{0}^{{1 - \alpha }} - At]}\nolimits^{1/(\alpha - 1)} }},$
из которого имеем, что $T < + \infty $. Действительно, найдется такой момент времени
${{T}_{0}} \in [0,{{T}_{1}}],\quad {{T}_{1}} = \frac{{\Phi _{0}^{{1 - \alpha }}}}{A},$
что верно
$\mathop {\lim \sup }\limits_{t \uparrow {{T}_{0}}} \Phi (t) = + \infty .$
Таким образом, справедлива

Лемма 8. Пусть неотрицательная функция $\Phi (t) \in {{\mathbb{C}}^{{(2)}}}([0,T])$ при некотором $T > 0$ удовлетворяет интегродифференциальному неравенству (10.1), при указанных условиях (10.2), (10.3) на функцию $h(t)$ и выполнены начальные условия

$\Phi {\kern 1pt} '(0) > \mathop {\left( {\frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} + \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{1 + \lambda }}}} \right)}\nolimits^{1/2} ,\quad \Phi (0) > 0,$
а также условие
$0 \leqslant \lambda < 2\alpha - 1,$
тогда $T < + \infty $. Более того, найдется такой момент времени
${{T}_{0}} \leqslant {{T}_{1}},\quad {{T}_{1}} = \frac{{{{\Phi }^{{1 - \alpha }}}(0)}}{A},$
что
$\mathop {\lim \sup }\limits_{t \uparrow {{T}_{0}}} \Phi (t) = + \infty ,$
где

$\alpha > 1,\quad {{\gamma }_{3}} \geqslant 0,\quad {{\gamma }_{4}} = {{\gamma }_{1}} + {{\gamma }_{2}}{\text{|}}\chi (0){\text{|}},\quad {{\gamma }_{1}} \geqslant 0,\quad {{\gamma }_{2}} \geqslant 0,\quad \delta = 1 + \frac{{\lambda - \alpha }}{{1 - \alpha }},$
${{A}^{2}} = {{(\alpha - 1)}^{2}}{{\Phi }^{{ - 2\alpha }}}(0)\left[ {{{{(\Phi {\kern 1pt} '(0))}}^{2}} - \frac{{{{\gamma }_{4}}}}{{\alpha - 1}}{{{(\Phi (0))}}^{2}} - \frac{{2{{\gamma }_{3}}}}{{(\alpha - 1)\delta }}{{{(\Phi (0))}}^{{\delta (1 - \alpha ) + 2\alpha }}}} \right] > 0.$

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

  1. Levine H.A. Some nonexistence and instability theorems for solutions of formally parabolic equations of the form $P{{u}_{t}} = - Au + F(u)$ // Arch. Ration. Mech. and Analys. 1974. V. 192. P. 1–21.

  2. Levine H.A. Instability and nonexistence of global solutions to nonlinear wave equations of the form $P{{u}_{{tt}}} = - Au + F(u)$ // Trans. Amer. Math. Soc. 1973. V. 51. P. 371–386.

  3. Straughan B. Further global nonexistence theorems for abstract nonlinear wave equations // Proc. Amer. Math. Soc. 1075. V. 48. № 2. P. 381–390.

  4. Калантаров В.К., Ладыженская О.А. О возникновении коллапсов для квазилинейных уравнений параболического и гиперболического типов // Зап. научн. сем. ЛОМИ. 1977. Т. 69. № 10. С. 77–102.

  5. Levine H.A., Serrin J. Global nonexistence theorems for quasilinear evolution equations with dissipation // Arch. Rational Mech. Anal. 1997. V. 137. № 4. P. 341–361.

  6. Levine H.A., Park S.R., Serrin J. Global existence and nonexistence theorems for quasilinear evolution equations of formally parabolic type // J. Differential Equations. 1998. V. 142. № 1. P. 212–229.

  7. Pucci P., Serrin J. Some new results on global nonexistence for abstract evolution with positive initial energy // Topol. Methods Nonlinear Analys. 1997. V. 10. № 2. P. 241–247.

  8. Pucci P., Serrin J.Global nonexistence for abstract evolution equations with positive initial energy // J. Differential Equations. 1998. V. 150. № 1. P. 203–214.

  9. Levine H.A., Todorova G. Blow up of solutions of the Cauchy problem for a wave equation with nonlinear damping and source terms and positive initial energy in four space-time dimensions // Proc. Amer. Math. Soc. 2001. V. 129. № 3. P. 793–805.

  10. Корпусов М.О. Разрушение решений обобщенного уравнения Клейна–Гордона с сильной диссипацией // Изв. РАН. Сер. матем. 2013. Т. 77. № 2. С. 109–138.

  11. Корпусов М.О. О разрушении решения одной нелинейной системы уравнений с положительной энергией в теории поля // Ж. вычисл. матем. и матем. физ. 2018. Т. 58. № 3. С. 447–458.

  12. Bilgin B.A., Kalantarov V.K. Non-existence of global solutions to nonlinear wave equations with positive initial energy // Commun. Pure Appl. Anal. 2018. V. 17. № 3. P. 987–999.

  13. Georgiev V., Todorova G. Existence of a solution of the wave equation with nonlinear damping and source terms // J. Differential Equations. 1994. V. 109. P. 295–308.

  14. Messaoudi S.A. Blow up and global existence in a nonlinear viscoelastic wave equation // Math. Nachr. 2003. V. 260. P. 58–66.

  15. Messaoudi S.A. Blow-up of positive-initial-energy solutions of a nonlinear viscoelastic hyperbolic equation // J. Math. Anal. Appl. 2006. V. 320. № 2. P. 902–915.

  16. Корпусов М.О. О разрушении решений класса сильно нелинейных уравнений типа Соболева // Известия РАН. Серия математическая. 2004. Т. 68. № 4. С. 151–204.

  17. Корпусов М.О., Свешников А.Г. О разрушении решений класса сильно нелинейных волновых диссипативных уравнений типа Соболева с источниками // Известия РАН. Серия математическая. 2004. Т. 69. № 4. С. 89–128.

  18. Корпусов М.О., Панин А.А. О разрушении решения абстрактной задачи Коши для формально гиперболического уравнения с двойной нелинейностью // Известия РАН. Серия математическая. 2014. Т. 78. № 5. С. 91–142.

  19. Al’shin A.B., Korpusov M.O., Sveshnikov A.G. Blow-up in nonlinear Sobolev type equations. De Gruyter Ser. Nonlinear Anal. Appl., 2011.

  20. Гаевский Х., Грегер К., Захариас К. Нелинейные операторные уравнения и операторные дифференциальные уравнения. М.: Мир, 1978.

  21. Тамм И.Е. Основы теории электричества. М.: Наука, 1989.

  22. Бонч-Бруевич В.Л., Калашников С.Г. Физика полупроводников. М.: Наука, 1990.

  23. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. Теоретическая физика. Том 8. М.: Наука, 1992.

  24. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. Теоретическая физика. Том 10. М.: Физматлит, 2001.

  25. Корпусов М.О., Свешников А.Г. О разрушении решения системы уравнений Осколкова // Матем. сб. 2009. Т. 200. № 4. С. 83–108.

  26. Панин А.А. О локальной разрешимости и разрушении решения абстрактного нелинейного интегрального уравнения Вольтерра // Матем. заметки. 2015. Т. 97. № 6. С. 884–903.

  27. Богачев В.И., Смолянов О.Г. Действительный и функциональный анализ: университетский курс. Ижевск: РХД, 2011.

  28. Картан А. Дифференциальное исчисление. Дифференциальные формы. М.: Мир, 1971.

  29. Демидович Б.П. Лекции по математической теории устойчивости. М.: Наука, 1967.

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