Журнал вычислительной математики и математической физики, 2021, T. 61, № 2, стр. 256-267

Угловой пограничный слой в краевых задачах для сингулярно возмущенных параболических уравнений с кубическими нелинейностями

И. В. Денисов *

Тульский государственный педагогический университет им. Л.Н. Толстого
300026 Тула, пр-т Ленина, 125, Россия

* E-mail: den@tspu@mail.ru

Поступила в редакцию 04.06.2000
После доработки 23.07.2000
Принята к публикации 16.09.2020

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

Аннотация

Для сингулярно возмущенного параболического уравнения

${{\epsilon }^{2}}\left( {{{a}^{2}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} - \frac{{\partial u}}{{\partial t}}} \right) = F(u,x,t,\epsilon )$

в прямоугольнике рассматривается задача с краевыми условиями I рода. Предполагается, что в угловых точках прямоугольника функция $F$ относительно переменной $u$ является кубической. Строится полное асимптотическое разложение решения при $\varepsilon \to 0$ и обосновывается его равномерность в замкнутом прямоугольнике. Библ. 10.

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

ВВЕДЕНИЕ

В работе рассматривается начально-краевая задача вида

(0.1)
${{\epsilon }^{2}}\left( {{{a}^{2}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} - \frac{{\partial u}}{{\partial t}}} \right) = F(u,x,t,\epsilon ),\quad (x,t) \in \Omega ,$
(0.2)
$u(x,0,\epsilon ) = \phi (x),\quad 0 \leqslant x \leqslant 1,$
(0.3)
$u(0,t,\epsilon ) = {{\psi }_{1}}(t),\quad u(1,t,\epsilon ) = {{\psi }_{2}}(t),\quad 0 \leqslant t \leqslant T.$

Через $\Omega $ обозначен прямоугольник $\{ (x,t)\,{\text{|}}\,0 < x < 1,0 < t < T\} $.

Ранее в работах [1]– [5] проведены исследования в предположении, что в угловых точках прямоугольника $(0,0)$ и $(1,0)$ функция $F$ является квадратичной относительно переменной $u$ на промежутке от корня вырожденного уравнения до граничного значения. Построено полное асимптотическое приближение решения задачи (0.1)–(0.3) при $\epsilon \to 0$ и обоснована равномерность этого приближения в замкнутом прямоугольнике с точностью любого порядка.

В работе [6] была предпринята попытка распространить полученные результаты на произвольную монотонную нелинейность. Для доказательства существования подходящих решений нелинейных задач использовался метод верхних и нижних решений. Однако выполнение необходимых неравенств в общем случае доказать не удалось и пришлось их выполнение постулировать. Кроме рассмотренных ранее квадратичных нелинейностей этим неравенствам удовлетворяли и некоторые другие нелинейные функции, что, несомненно, являлось продвижением в исследовании общей задачи. Однако класс подходящих нелинейных функций не был полностью изучен. В данной статье рассматриваются кубические нелинейности. Для предлагаемых барьерных функций необходимые неравенства доказываются. Строится полное асимптотическое приближение решения задачи (0.1)–(0.3) при $\epsilon \to 0$ и обосновывается равномерность этого приближения в замкнутом прямоугольнике с точностью любого порядка.

1. ПОСТАНОВКА ЗАДАЧИ

Пусть, как и в предыдущих работах, выполнены следующие условия.

Условие 1. Функции $F(u,x,t,\epsilon )$, $\phi (x)$, ${{\psi }_{1}}(t)$ и ${{\psi }_{2}}(t)$ являются достаточно гладкими и в угловых точках прямоугольника $\Omega $ выполняются условия согласованности начально-краевых значений:

$\phi (0) = {{\psi }_{1}}(0),\quad \phi (1) = {{\psi }_{2}}(0).$

Условие 2. Вырожденное уравнение $F(u,x,t,0) = 0$ в замкнутом прямоугольнике $\bar {\Omega }$ имеет решение, которое обозначается как $u = {{\bar {u}}_{0}}(x,t)$.

Заметим, что в силу нелинейности это уравнение может иметь и другие решения.

Условие 3. Производная $F_{u}^{'}({{\bar {u}}_{0}}(x,t),x,t,0) > 0$ в замкнутом прямоугольнике $\bar {\Omega }$.

Условие 4. Начальная задача

(1.1)
$\frac{{d{{\Pi }_{0}}}}{{d\tau }} = - F({{\bar {u}}_{0}}(x,0) + {{\Pi }_{0}},x,0,0),\quad {{\Pi }_{0}}(x,0) = \phi (x) - {{\bar {u}}_{0}}(x,0),$
имеет решение ${{\Pi }_{0}}(x,\tau )$ при $\tau \geqslant 0$, удовлетворяющее условию ${{\Pi }_{0}}(x,\infty ) = 0$ (здесь параметр $x \in [0,1]$).

Условие 5. Для систем

(1.2)
$\frac{{d{{z}_{1}}}}{{dy}} = {{z}_{2}},\quad {{a}^{2}}\frac{{d{{z}_{2}}}}{{dy}} = F({{\bar {u}}_{0}}(k,t) + {{z}_{1}},k,t,0),$
прямые ${{z}_{1}} = {{\psi }_{{1 + k}}}(t) - {{\bar {u}}_{0}}(k,t)$ пересекают сепаратрисы, входящие в точку покоя $({{z}_{1}},{{z}_{2}}) = (0,0)$ при $y \to \infty $ (здесь $t$ – параметр, $k = 0$ или $1$).

Условий 1–5 недостаточно, чтобы гарантировать существование решения задачи (0.1)–(0.3) для произвольной функции $F(u,x,t,\epsilon )$. Поэтому требуются дополнительные условия, гарантирующие возможность построения асимптотики решения. Эти условия будут заключаться в выборе определенного класса функций. Решение задачи (0.1)–(0.3) строится согласно методу угловых пограничных функций (см. [7]) в виде суммы

(1.3)
$u(x,t,\epsilon ) = \bar {u} + (\Pi + Q + Q{\kern 1pt} *{\kern 1pt} ) + (P + P{\kern 1pt} *{\kern 1pt} ),$
где $\bar {u}$ обозначает функцию, называемую регулярной частью асимптотики. Эта функция представляет решение задачи во внутренней части прямоугольника $\Omega $ без учета граничных условий. Пограничные функции $\Pi $, $Q$ и $Q{\kern 1pt} *$ осуществляют гладкий переход от регулярной части к граничным условиям на сторонах прямоугольника $\Omega $: $t = 0$, $x = 0$ и $x = 1$ соответственно. Угловые пограничные функции $P$ и $P{\kern 1pt} *$ сглаживают невязки, вносимые пограничными функциями вблизи вершин прямоугольника $\Omega $: $(0,0)$ и $(1,0)$ соответственно.

2. РЕГУЛЯРНАЯ И ПОГРАНСЛОЙНАЯ ЧАСТИ АСИМПТОТИКИ

Формальная процедура построения регулярной части асимптотики и погранслойных функций хорошо отработана (см. [7]), однако, чтобы не обращаться к другим источникам, повторим ее схематично. В уравнении (0.1) функция $F$ заменяется выражением, аналогичным (1.3):

(2.1)
$F(u,x,t,\epsilon ) = \bar {F} + (\Pi F + QF + Q{\kern 1pt} *{\kern 1pt} F) + (PF + P{\kern 1pt} *{\kern 1pt} F).$

Выражения (1.3) и (2.1) подставляются в уравнение (0.1), которое разделяется на шесть частей: регулярную, три погранслойных и две угловых. Регулярная часть асимптотики $\bar {u}$ строится в виде ряда по степеням $\epsilon $:

(2.2)
$\bar {u}(x,t,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}{{\bar {u}}_{k}}(x,t).$

Коэффициент ${{\bar {u}}_{0}} = {{\bar {u}}_{0}}(x,t)$ выбирается в соответствии с условиями 2 и 3, а последующие функции ${{\bar {u}}_{k}},k \geqslant 1$, строятся рекуррентно.

Пусть $\partial \Omega $ обозначает границу прямоугольника $\Omega $ без стороны $t = T$. Регулярная часть $\bar {u}(x,t,\epsilon )$ асимптотики дает решение задачи (0.1)–(0.3) внутри прямоугольника $\Omega $, но на границе $\partial \Omega $ функция $\bar {u}(x,t,\epsilon )$, вообще говоря, не совпадает с начальным и граничными значениями. В связи с этим возникает так называемая “невязка”.

Погранслойная часть асимптотики вводится для устранения невязок регулярной части с начальным и граничными условиями. Погранслойные функции $\Pi $, $Q$ и $Q{\kern 1pt} *$ определяются из уравнений, в которых переходят к растянутым переменным

$\xi = \frac{x}{\epsilon },\quad {{\xi }_{*}} = \frac{{1 - x}}{\epsilon },\quad \tau = \frac{t}{{{{\epsilon }^{2}}}}.$

Функции $\Pi $, $Q$ и $Q{\kern 1pt} *$ ищутся в виде рядов

(2.3)
$\Pi (x,\tau ,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}{{\Pi }_{k}}(x,\tau ),$
(2.4)
$Q(\xi ,t,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}{{Q}_{k}}(\xi ,t),$
(2.5)
$Q{\kern 1pt} *({{\xi }_{*}},t,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}Q_{k}^{*}({{\xi }_{*}},t).$

На стороне $t = 0$ невязки в начальном условии (0.2) призвана устранить функция $\Pi = \Pi (x,\tau ,\epsilon )$. При переходе от переменных $(x,t)$ к переменным $(x,\tau )$ прямоугольник $\Omega $ при $\epsilon \to 0$ растягивается до полуполосы $0 \leqslant x \leqslant 1$, $0 \leqslant \tau < \infty $.

Для ${{\Pi }_{0}} = {{\Pi }_{0}}(x,\tau )$ получается задача

(2.6)
$ - \frac{{\partial {{\Pi }_{0}}}}{{\partial \tau }} = F\left( {{{{\bar {u}}}_{0}}(x,0) + {{\Pi }_{0}},x,0,0} \right),\quad {{\Pi }_{0}}(x,0) = \phi (x) - {{\bar {u}}_{0}}(x,0).$

Здесь $x$ играет роль параметра. В силу условия 4 эта задача имеет решение, для которого в силу условия 3 справедлива экспоненциальная оценка убывания вида

(2.7)
$\left| {{{\Pi }_{0}}(x,\tau )} \right| \leqslant Cexp( - \kappa \tau ),$
где $C$ и $\kappa $ – некоторые положительные числа.

Задачи для определения функций ${{\Pi }_{k}} = {{\Pi }_{k}}(x,\tau )$, $k \geqslant 1$, получаются линейными:

(2.8)
$ - \frac{{\partial {{\Pi }_{k}}}}{{\partial \tau }} = F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(x,0) + {{\Pi }_{0}},x,0,0} \right){{\Pi }_{k}} + {{\tilde {\pi }}_{k}},\quad {{\Pi }_{k}}(x,0) = - {{\bar {u}}_{k}}(x,0).$

Функции ${{\tilde {\pi }}_{k}}$ рекуррентно выражаются через функции ${{\Pi }_{j}}$, $j < k$, и их производные. Поэтому, если для функций ${{\Pi }_{j}}$, $j < k$, справедливы оценки вида (2.7), то для функций ${{\tilde {\pi }}_{k}}$ справедливы оценки того же вида.

Если величина $\phi (x) - {{\bar {u}}_{0}}(x,0)$ не равна тождественно нулю, то решения задач (2.8) имеют вид

${{\Pi }_{k}}(x,\tau ) = - U(x,\tau ){{\bar {u}}_{k}}(x,0) - U(x,\tau )\int\limits_0^\tau \,\mathop {\left( {U(x,\sigma )} \right)}\nolimits^{ - 1} {{\tilde {\pi }}_{k}}(x,\sigma )d\sigma ,$
где
$U(x,\tau ) = exp\left( { - \int\limits_0^\tau \,F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(x,0) + {{\Pi }_{0}}(x,\lambda ),x,0,0} \right)d\lambda } \right)$
есть фундаментальное решение $(U(x,0) = 1)$ соответствующего однородного уравнения, и справедлива оценка вида
$\left| {U(x,\tau )\mathop {\left( {U(x,\sigma )} \right)}\nolimits^{ - 1} } \right| \leqslant Cexp( - \kappa (\tau - \sigma )),$
где переменные $0 \leqslant x \leqslant 1,$ $0 \leqslant \sigma \leqslant \tau $, а постоянные $C$ и $\kappa $ – положительные числа. Эта оценка позволяет для функции ${{\Pi }_{k}}(x,\tau )$ получить оценку вида (2.7).

Если величина $\phi (x) - {{\bar {u}}_{0}}(x,0) \equiv 0$, то ${{\Pi }_{0}}(x,\tau ) \equiv 0$. Коэффициенты при ${{\Pi }_{k}}$ в задачах (2.8) оказываются постоянными и положительными, т.е. задачи упрощаются.

Таким образом определяются коэффициенты ряда (2.3), и функция $\Pi (x,\tau ,\epsilon )$ устраняет невязки с начальным условием (0.2) на стороне $t = 0$.

Построенная регулярная часть асимптотики вносит невязки и в граничные условия. На стороне $x = 0$ невязки в граничных условиях призвана устранить функция $Q = Q(\xi ,t,\epsilon )$, где $\xi = x{\text{/}}\epsilon $ – растянутая переменная. При переходе от переменных $(x,t)$ к переменным $(\xi ,t)$ прямоугольник $\Omega $ при $\epsilon \to 0$ растягивается до полуполосы $0 \leqslant \xi < \infty $, $0 \leqslant t \leqslant T$.

Задача для ${{Q}_{0}} = {{Q}_{0}}(\xi ,t)$ имеет вид

(2.9)
${{a}^{2}}\frac{{{{\partial }^{2}}{{Q}_{0}}}}{{\partial {{\xi }^{2}}}} = F\left( {{{{\bar {u}}}_{0}}(0,t) + {{Q}_{0}},0,t,0} \right),\quad {{Q}_{0}}(0,t) = {{\psi }_{1}}(t) - {{\bar {u}}_{0}}(0,t),\quad {{Q}_{0}}(\infty ,t) = 0,$
где $t$ играет роль параметра. Уравнение (2.9) эквивалентно системе (1.2), в которой следует положить ${{z}_{1}} = {{Q}_{0}}(\xi ,t),$ $k = 0,$ $y = \xi $. Условия затухания выделяют решения уравнения (2.9), для которых справедливы экспоненциальные оценки убывания вида
(2.10)
$\left| {{{Q}_{0}}(\xi ,t)} \right| \leqslant Cexp( - \kappa \xi ),$
где $C$ и $\kappa $ – положительные числа. Так как возможен переход с сепаратрисы на сепаратрису, то решение задачи (2.9) не единственно. Однако такие случаи мы исключаем и рассматриваем только монотонные решения.

Задачи для определения функций ${{Q}_{k}}(\xi ,t),$ $k \geqslant 1$, линейны:

(2.11)
${{a}^{2}}\frac{{{{\partial }^{2}}{{Q}_{k}}}}{{\partial {{\xi }^{2}}}} = F_{u}^{'}({{\bar {u}}_{0}}(0,t) + {{Q}_{0}},0,t,0){{Q}_{k}} + {{\tilde {q}}_{k}},\quad {{Q}_{k}}(0,t) = - {{\bar {u}}_{k}}(0,t),\quad {{Q}_{k}}(\infty ,t) = 0.$

Функции ${{\tilde {q}}_{k}}$ рекуррентно выражаются через функции ${{Q}_{j}}$, $j < k$, и их производные. Поэтому, если для функций ${{Q}_{j}}$, $j < k$, справедливы экспоненциальные оценки вида (2.10), то для функций ${{\tilde {q}}_{k}}$ справедливы оценки того же вида.

Если величина ${{\psi }_{1}}(t) - {{\bar {u}}_{0}}(0,t)$ не равна тождественно нулю, то решение задачи (2.11) имеет вид

${{Q}_{k}}(\xi ,t) = - \frac{{\Phi (\xi ,t)}}{{\Phi (0,t)}}{{\bar {u}}_{k}}(0,t) - \frac{{\Phi (\xi ,t)}}{{a(0,t)}}\int\limits_0^\xi \,\frac{{d\lambda }}{{{{\Phi }^{2}}(\lambda ,t)}}\int\limits_\lambda ^\infty \,\Phi (\sigma ,t){{\tilde {q}}_{k}}(\sigma ,t)d\sigma ,\quad \Phi (\xi ,t) = \frac{{\partial {{Q}_{0}}(\xi ,t)}}{{\partial \xi }},$
и для него справедлива оценка вида (2.10).

Если величина ${{\psi }_{1}}(t) - {{\bar {u}}_{0}}(0,t) \equiv 0$, то ${{Q}_{0}}(\xi ,t) \equiv 0$, а коэффициент при ${{Q}_{k}}$ в уравнениях (2.11) оказывается постоянным и положительным, т.е. задача упрощается.

Таким образом определяются коэффициенты ряда (2.4), и функция $Q(\xi ,t,\epsilon )$ устраняет невязки в граничном условии на стороне $x = 0$.

Регулярная часть асимптотики вносит невязки в граничные условия и на стороне $x = 1$. Эти невязки устраняет функция $Q{\kern 1pt} * = Q{\kern 1pt} *({{\xi }_{*}},t,\epsilon )$, ${{\xi }_{*}} = (1 - x){\text{/}}\epsilon $, которая строится в виде ряда (2.5). Коэффициенты этого ряда $Q_{k}^{*}$ определяются аналогично коэффициентам ряда (2.4) и для них справедливы экспоненциальные оценки убывания вида

${\text{|}}Q_{k}^{*}({{\xi }_{*}},t){\kern 1pt} {\text{|}} \leqslant Cexp( - \kappa {{\xi }_{*}}),$
где $C$ и $\kappa $ – положительные числа.

Таким образом, погранслойная часть асимптотики определяется полностью. Однако каждая в отдельности погранслойная функция, устраняя невязки на соответствующей стороне, в свою очередь вносит невязки на примыкающие стороны прямоугольника. Так, погранслойные функции ${{\Pi }_{k}}(x,\tau )$, устраняя невязки в начальном условии на стороне $t = 0$, вносят дополнительные невязки в граничные условия на сторонах $x = 0$ и $x = 1$. Эти невязки существенны только вблизи угловых точек (0,0) и (1,0), а далее, с ростом $t$, они экспоненциально затухают. Аналогичное влияние функции ${{Q}_{k}}(\xi ,t)$ и $Q_{k}^{*}({{\xi }_{*}},t)$ оказывают на начальное условие на стороне $x = 0$.

3. УГЛОВАЯ ЧАСТЬ АСИМПТОТИКИ. ОСНОВНЫЕ ПРОБЛЕМЫ

С целью устранения невязок с начальным и граничными условиями вблизи угловых точек $(0,0)$ и $(1,0)$ прямоугольника $\Omega $ вводятся угловые пограничные функции $P(\xi ,\tau ,\epsilon )$ и $P{\kern 1pt} *{\kern 1pt} ({{\xi }_{*}},\tau ,\epsilon )$. Ввиду того, что не существует универсального метода их нахождения, задачи для определения функций $P(\xi ,\tau ,\epsilon )$ и $P{\kern 1pt} *{\kern 1pt} ({{\xi }_{*}},\tau ,\epsilon )$ доставляют основные трудности. Практически каждая такая задача требует разработки новых методов решения.

Угловые пограничные функции будем искать в виде рядов

$P(\xi ,\tau ,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}{{P}_{k}}(\xi ,\tau ),\quad P{\kern 1pt} *{\kern 1pt} ({{\xi }_{*}},\tau ,\epsilon ) = \sum\limits_{k = 0}^\infty \,{{\epsilon }^{k}}P_{k}^{*}({{\xi }_{*}},\tau ).$

Рассмотрим угловую точку (0,0). В окрестности этой точки невязки в условия (0.2), (0.3) вносят функции $\Pi (x,\tau ,\epsilon )$ и $Q(\xi ,t,\epsilon )$. Для устранения этих невязок служит функция $P(\xi ,\tau ,\epsilon )$, которая определяется из уравнения

${{\left. {{{\epsilon }^{2}}\left( {{{a}^{2}}\frac{{{{\partial }^{2}}P}}{{\partial {{x}^{2}}}} - \frac{{\partial P}}{{\partial t}}} \right)} \right|}_{\begin{subarray}{l} x = \varepsilon \xi \\ t = {{\varepsilon }^{2}}\tau \end{subarray} }} = PF.$

Задача для определения ${{P}_{0}}(\xi ,\tau )$ ставится в первой четверти

$\mathbb{R}_{ + }^{2}: = \{ (\xi ,\tau )\,{\text{|}}\,\xi > 0,\;\tau > 0\} $
плоскости переменных $(\xi ,\tau )$ и имеет вид

(3.1)
$\begin{gathered} {{a}^{2}}\frac{{{{\partial }^{2}}{{P}_{0}}}}{{\partial {{\xi }^{2}}}} - \frac{{\partial {{P}_{0}}}}{{\partial \tau }} = F\left( {{{{\bar {u}}}_{0}}(0,0) + {{\Pi }_{0}}(0,\tau ) + {{Q}_{0}}(\xi ,0) + {{P}_{0}}(\xi ,\tau ),0,0,0} \right) - \\ \, - F\left( {{{{\bar {u}}}_{0}}(0,0) + {{\Pi }_{0}}(0,\tau ),0,0,0} \right) - F\left( {{{{\bar {u}}}_{0}}(0,0) + {{Q}_{0}}(\xi ,0),0,0,0} \right),\quad \xi > 0,\quad \tau > 0, \\ \end{gathered} $
(3.2)
${{P}_{0}}(0,\tau ) = - {{\Pi }_{0}}(0,\tau ),\quad {{P}_{0}}(\xi ,0) = - {{Q}_{0}}(\xi ,0),\quad {{P}_{0}}(\xi ,\tau ) \to 0\quad {\text{при}}\quad \xi + \tau \to \infty .$

Для функций ${{P}_{k}}(\xi ,\tau ),$ $k \geqslant 1$, в области $\mathbb{R}_{ + }^{2}$ получаются линейные задачи

(3.3)
${{a}^{2}}\frac{{{{\partial }^{2}}{{P}_{k}}}}{{\partial {{\xi }^{2}}}} - \frac{{\partial {{P}_{k}}}}{{\partial \tau }} = F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(0,0) + {{\Pi }_{0}}(0,\tau ) + {{Q}_{0}}(\xi ,0) + {{P}_{0}}(\xi ,\tau ),0,0,0} \right){{P}_{k}} + {{h}_{k}},$
(3.4)
${{P}_{k}}(0,\tau ) = - {{\Pi }_{k}}(0,\tau ),\quad {{P}_{k}}(\xi ,0) = - {{Q}_{k}}(\xi ,0),\quad {{P}_{k}}(\xi ,\tau ) \to 0\quad {\text{при}}\quad \xi + \tau \to \infty ,$
где неоднородности ${{h}_{k}}$ удовлетворяют экспоненциальным оценкам убывания вида
(3.5)
$\left| {{{h}_{k}}(\xi ,\tau )} \right| \leqslant C\exp ( - \kappa (\xi + \tau )),$
если такого же вида оценкам удовлетворяют функции ${{P}_{0}}, \ldots ,{{P}_{{k - 1}}}$. Здесь $C$ и $\kappa $ – некоторые положительные числа.

Если число $\phi (0,0) - {{\bar {u}}_{0}}(0,0) = 0$, то решением задачи (2.6) при $x = 0$ будет функция ${{\Pi }_{0}}(0,\tau ) \equiv 0$, решением задачи (2.8) при $t = 0$ будет функция ${{Q}_{0}}(\xi ,0) \equiv 0$. Решением задачи (3.1), (3.2) будет функция ${{P}_{0}}(\xi ,\tau ) \equiv 0$, а коэффициент в задачах (3.3), (3.4) будет постоянным и положительным:

$F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(0,0) + {{\Pi }_{0}}(0,\tau ) + {{Q}_{0}}(\xi ,0) + {{P}_{0}}(\xi ,\tau ),0,0,0} \right) = F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(0,0),0,0,0} \right) > 0.$

В этом случае решения задач (3.3), (3.4) выписываются в явном виде и для них получаются экспоненциальные оценки вида (3.5).

Если число $\varphi (0,0) - {{\bar {u}}_{0}}(0,0) \ne 0$, то, вообще говоря, неизвестно, имеет или нет задача (3.1), (3.2) решение и удовлетворяет ли решение в случае существования экспоненциальной оценке вида (3.5). Кроме этого, в задачах (3.3), (3.4) коэффициент

$F_{u}^{'}: = F_{u}^{'}\left( {{{{\bar {u}}}_{0}}(0,0) + {{\Pi }_{0}}(0,\tau ) + {{Q}_{0}}(\xi ,0) + {{P}_{0}}(\xi ,\tau ),0,0,0} \right)$
может в зависимости от вида функции $F$ и величины $\phi (0,0)$ принимать как положительные, так и отрицательные значения.

Даже в случае разрешимости задач (3.1)–(3.4) доказательство того, что задача (0.1)–(0.3) имеет решение, все равно остается проблемой. Это связано с тем, что, не зная явного вида функции ${{P}_{0}}(\xi ,\tau )$, мы не можем знать явного вида коэффициента $F_{u}^{'}$, который может оказаться как положительным, так и отрицательным. Если не накладывать дополнительных условий, то это обстоятельство, вообще говоря, не позволит обосновать построенную асимптотику решения.

Таким образом, в результате реализации метода угловых пограничных функций для исследования задачи (0.1)–(0.3) мы перешли к задачам (3.1)–(3.4). Дальнейшие исследования предполагают разрешение, по крайней мере, трех основных проблем.

1. Имеет ли задача (3.1), (3.2) решение, удовлетворяющее экспоненциальной оценке убывания вида (3.5)?

2. Если задача (3.1), (3.2) имеет решение, удовлетворяющее экспоненциальной оценке вида (3.5), то имеют ли задачи (3.3), (3.4) решения, удовлетворяющие подобным оценкам?

3. Если задачи (3.1)–(3.4) разрешимы, т.е. если ряд (1.3) может быть построен, то имеет ли задача (0.1)–(0.3) решение, для которого этот ряд будет асимптотическим представлением при $\epsilon \to 0$ в замкнутом прямоугольнике $\bar {\Omega }$?

Задачи для угловых пограничных функций $P_{k}^{*}({{\xi }_{*}},\tau ),$ ставятся аналогично.

4. КУБИЧЕСКАЯ НЕЛИНЕЙНОСТЬ

Будем рассматривать угловую точку $(0,0)$ прямоугольника $\Omega $ и определим оператор $L$:

$L(P): = {{a}^{2}}\frac{{{{\partial }^{2}}P}}{{\partial {{\xi }^{2}}}} - \frac{{\partial P}}{{\partial \tau }} - F({{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} + P) + F({{\bar {u}}_{0}} + {{\Pi }_{0}}) + F({{\bar {u}}_{0}} + {{Q}_{0}}),$
где

$P = P(\xi ,\tau ),\quad F(u) = F(u,0,0,0),\quad {{\bar {u}}_{0}} = {{\bar {u}}_{0}}(0,0),\quad {{\Pi }_{0}} = {{\Pi }_{0}}(0,\tau ),\quad {{Q}_{0}} = {{Q}_{0}}(\xi ,0).$

Задачу (3.1), (3.2) можно переписать в операторной форме:

(4.1)
$L({{P}_{0}}) = 0\quad {\text{в}}\;{\text{области}}\quad \mathbb{R}_{ + }^{2},$
(4.2)
${{P}_{0}}(0,\tau ) = - {{\Pi }_{0}}(0,\tau ),\quad {{P}_{0}}(\xi ,0) = - {{Q}_{0}}(\xi ,0),\quad {{P}_{0}}(\xi ,\tau ) \to 0\quad {\text{при}}\quad \xi + \tau \to \infty .$

Для доказательства существования решения задачи (4.1), (4.2) используется метод верхних и нижних решений (см. [8]– [10]), который заключается в том, что задача

$L(Z) = 0\quad {\text{в}}\;{\text{области}}\quad \Omega ,\quad Z = h\quad {\text{на}}\;{\text{границе}}\;\quad \partial \Omega $
имеет решение $Z$ в промежутке между барьерными функциями ${{Z}_{ - }} \leqslant Z \leqslant {{Z}_{ + }}$, если в области $\Omega $ выполняются неравенства
$L({{Z}_{ + }}) \leqslant 0,\quad L({{Z}_{ - }}) \geqslant 0,\quad {{Z}_{ - }} \leqslant {{Z}_{ + }},$
а на ее границе

${{Z}_{ - }} \leqslant h \leqslant {{Z}_{ + }}.$

Барьерные функции для задачи (4.1), (4.2) желательно строить с расчетом, чтобы коэффициент $F_{u}^{'}$ в задачах (3.3), (3.4) был положительным. Это обеспечит существование решений ${{P}_{k}}(\xi ,\tau )$ с оценками вида (3.5). В качестве возможного верхнего решения задачи (4.1), (4.2) рассмотрим функцию

${{Z}_{ + }}(\xi ,\tau ) = 0.$

Для определенности будем считать, что $\phi (0,0) > {{\bar {u}}_{0}}(0,0)$, т.е. граничное значение в точке $(0,0)$ лежит правее корня вырожденного уравнения. Тогда значения ${{\Pi }_{0}}$ и ${{Q}_{0}}$ принадлежат промежутку $(0,\phi - {{\bar {u}}_{0}}]$, где $\phi - {{\bar {u}}_{0}} > 0$, и на границе области $\mathbb{R}_{ + }^{2}$ выполняются необходимые для верхнего решения неравенства:

${{Z}_{ + }}(0,\tau ) = 0 \geqslant - {{\Pi }_{0}},\quad {{Z}_{ + }}(\xi ,0) = 0 \geqslant - {{Q}_{0}}.$

Внутри области $\mathbb{R}_{ + }^{2}$ нужно доказать неравенство

(4.3)
$L({{Z}_{ + }}) = - F({{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}}) + F({{\bar {u}}_{0}} + {{\Pi }_{0}}) + F({{\bar {u}}_{0}} + {{Q}_{0}}) \leqslant 0.$

Для краткости обозначим ${{\Pi }_{0}} = y$, ${{Q}_{0}} = z$, $L({{Z}_{ + }}) = L$, где $L = L(y,z)$. Тогда неравенство (4.3) примет вид

$L = - F({{\bar {u}}_{0}} + y + z) + F({{\bar {u}}_{0}} + y) + F({{\bar {u}}_{0}} + z) \leqslant 0,\quad y,z \in (0,\phi - {{\bar {u}}_{0}}].$

Необходимые условия экстремума

$\frac{{\partial L}}{{\partial y}} = - F{\kern 1pt} '({{\bar {u}}_{0}} + y + z) + F{\kern 1pt} '({{\bar {u}}_{0}} + y) = 0,$
$\frac{{\partial L}}{{\partial z}} = - F{\kern 1pt} '({{\bar {u}}_{0}} + y + z) + F{\kern 1pt} '({{\bar {u}}_{0}} + z) = 0,$
приводят к равенствам

(4.4)
$F{\kern 1pt} '({{\bar {u}}_{0}} + y + z) = F{\kern 1pt} '({{\bar {u}}_{0}} + y) = F{\kern 1pt} '({{\bar {u}}_{0}} + z).$

На границе области из-за симметрии функции $L(y,z)$ можно рассмотреть только две из четырех сторон квадрата $(0,\phi - {{\bar {u}}_{0}}] \times (0,\phi - {{\bar {u}}_{0}}]$.

1. При $y = 0$ значения

$L(0,z) = - F({{\bar {u}}_{0}} + z) + F({{\bar {u}}_{0}}) + F({{\bar {u}}_{0}} + z) = 0.$

2. При $y = \phi - {{\bar {u}}_{0}}$ значения

$L(\phi - {{\bar {u}}_{0}},z) = - F(\phi + z) + F(\phi ) + F({{\bar {u}}_{0}} + z).$

Необходимое неравенство $L(\phi - {{\bar {u}}_{0}},z) \leqslant 0$ выполняется, если

(4.5)
$F({{\bar {u}}_{0}} + u) + F(\phi ) \leqslant F(\phi + u),\quad u \in (0,\phi - {{\bar {u}}_{0}}].$

Для проверки условия (4.5) нужно рассмотреть два куска графика функции $F(u)$: один – при $u \in ({{\bar {u}}_{0}},\phi ]$, а другой – при $u \in (\phi ,\phi + \phi - {{\bar {u}}_{0}}]$. Первый кусок нужно сдвинуть параллельно вертикальной оси на величину $F(\phi )$, в связи с чем появляется функция

${{F}_{1}}(u) = F({{\bar {u}}_{0}} + u) + F(\phi ),\quad u \in (0,\phi - {{\bar {u}}_{0}}].$

Второй кусок нужно сдвинуть влево на величину $\phi - {{\bar {u}}_{0}}$, в связи с чем появляется функция

${{F}_{2}}(u) = F(\phi + u),\quad u \in (0,\phi - {{\bar {u}}_{0}}].$

Графики обеих функций выходят из точки $(0,F(\phi ))$. Условие (4.5) оказывается эквивалентным простому условию

${{F}_{1}}(u) \leqslant {{F}_{2}}(u),\quad u \in (0,\phi - {{\bar {u}}_{0}}].$

В качестве возможного нижнего решения задачи (4.1), (4.2) рассмотрим функцию

${{Z}_{ - }}(\xi ,\tau ) = - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} .$

Так как значения ${{\Pi }_{0}}$ и ${{Q}_{0}}$ принадлежат промежутку $(0,\phi - {{\bar {u}}_{0}}]$, то необходимые для нижнего решения неравенства выполняются на границе области $\mathbb{R}_{ + }^{2}$:

${{Z}_{ - }}(0,\tau ) = - 2\sqrt {(\phi - {{{\bar {u}}}_{0}}){{\Pi }_{0}}} \leqslant - {{\Pi }_{0}},\quad {{Z}_{ - }}(\xi ,0) = - 2\sqrt {(\phi - {{{\bar {u}}}_{0}}){{Q}_{0}}} \leqslant - {{Q}_{0}}.$

Внутри области $\mathbb{R}_{ + }^{2}$ нужно доказать неравенство $L({{Z}_{ - }}) \geqslant 0$. Имеем

$L({{Z}_{ - }}) = {{a}^{2}}\frac{{{{\partial }^{2}}{{Z}_{ - }}}}{{\partial {{\xi }^{2}}}} - \frac{{\partial {{Z}_{ - }}}}{{\partial \tau }} - F({{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} ) + F({{\bar {u}}_{0}} + {{\Pi }_{0}}) + F({{\bar {u}}_{0}} + {{Q}_{0}}).$

Производная

$\frac{{{{\partial }^{2}}{{Z}_{ - }}}}{{\partial {{\xi }^{2}}}} = - 2\sqrt {{{\Pi }_{0}}} \frac{{{{d}^{2}}\sqrt {{{Q}_{0}}} }}{{d{{\xi }^{2}}}} = - 2\sqrt {{{\Pi }_{0}}} \frac{d}{{d\xi }}\left( {\frac{1}{{2\sqrt {{{Q}_{0}}} }}\frac{{d{{Q}_{0}}}}{{d\xi }}} \right) = - 2\sqrt {{{\Pi }_{0}}} \left( { - \frac{1}{{4{{Q}_{0}}\sqrt {{{Q}_{0}}} }}\mathop {\left( {\frac{{d{{Q}_{0}}}}{{d\xi }}} \right)}\nolimits^2 + \frac{1}{{2\sqrt {{{Q}_{0}}} }}\frac{{{{d}^{2}}{{Q}_{0}}}}{{d{{\xi }^{2}}}}} \right).$

Задача для определения функции ${{Q}_{0}} = {{Q}_{0}}(\xi ,0)$ получается из (2.9) при $t = 0$:

${{a}^{2}}\frac{{{{d}^{2}}{{Q}_{0}}}}{{d{{\xi }^{2}}}} = F({{\bar {u}}_{0}} + {{Q}_{0}}),\quad {{Q}_{0}}(0,0) = \phi - {{\bar {u}}_{0}},\quad {{Q}_{0}}(\infty ,0) = 0.$

Отсюда

$\frac{{{{d}^{2}}{{Q}_{0}}}}{{d{{\xi }^{2}}}} = \frac{1}{{{{a}^{2}}}}F({{\bar {u}}_{0}} + {{Q}_{0}}),\quad \mathop {\left( {\frac{{d{{Q}_{0}}}}{{d\xi }}} \right)}\nolimits^2 = \frac{2}{{{{a}^{2}}}}\int\limits_0^{{{Q}_{0}}} \,F({{\bar {u}}_{0}} + u)du.$

Поэтому

$\begin{gathered} \frac{{{{\partial }^{2}}{{Z}_{ - }}}}{{\partial {{\xi }^{2}}}} = - 2\sqrt {{{\Pi }_{0}}} \left[ { - \frac{1}{{4{{Q}_{0}}\sqrt {{{Q}_{0}}} }}\frac{2}{{{{a}^{2}}}}\int\limits_0^{{{Q}_{0}}} \,F({{{\bar {u}}}_{0}} + u)du + \frac{1}{{2\sqrt {{{Q}_{0}}} }}\frac{1}{{{{a}^{2}}}}F({{{\bar {u}}}_{0}} + {{Q}_{0}})} \right] = \\ \, = \frac{1}{{{{a}^{2}}}}\sqrt {{{\Pi }_{0}}{{Q}_{0}}} \left( {\frac{1}{{Q_{0}^{2}}}\int\limits_0^{{{Q}_{0}}} \,F({{{\bar {u}}}_{0}} + u)du - \frac{1}{{{{Q}_{0}}}}F({{{\bar {u}}}_{0}} + {{Q}_{0}})} \right). \\ \end{gathered} $

Производная

$\frac{{\partial {{Z}_{ - }}}}{{\partial \tau }} = - 2\sqrt {{{Q}_{0}}} \frac{{d\sqrt {{{\Pi }_{0}}} }}{{d\tau }} = - \frac{{\sqrt {{{Q}_{0}}} }}{{\sqrt {{{\Pi }_{0}}} }}\frac{{d{{\Pi }_{0}}}}{{d\tau }}.$

Задача для определения функции ${{\Pi }_{0}} = {{\Pi }_{0}}(0,\tau )$ получается из (2.6) при $y = 0$:

$ - \frac{{d{{\Pi }_{0}}}}{{d\tau }} = F({{\bar {u}}_{0}} + {{\Pi }_{0}}),\quad {{\Pi }_{0}}(0,0) = \phi - {{\bar {u}}_{0}}.$

Поэтому

$\frac{{\partial {{Z}_{ - }}}}{{\partial \tau }} = \frac{{\sqrt {{{Q}_{0}}} }}{{\sqrt {{{\Pi }_{0}}} }}F({{\bar {u}}_{0}} + {{\Pi }_{0}}) = \sqrt {{{\Pi }_{0}}{{Q}_{0}}} \frac{1}{{{{\Pi }_{0}}}}F({{\bar {u}}_{0}} + {{\Pi }_{0}}).$

В результате для $L({{Z}_{ - }})$ получается выражение

(4.6)
$\begin{gathered} L({{Z}_{ - }}) = \sqrt {{{\Pi }_{0}}{{Q}_{0}}} \left( {\frac{1}{{Q_{0}^{2}}}\int\limits_0^{{{Q}_{0}}} \,F({{{\bar {u}}}_{0}} + u)du - \frac{1}{{{{Q}_{0}}}}F({{{\bar {u}}}_{0}} + {{Q}_{0}}) - \frac{1}{{{{\Pi }_{0}}}}F({{{\bar {u}}}_{0}} + {{\Pi }_{0}})} \right) - \\ \, - F({{{\bar {u}}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} ) + F({{{\bar {u}}}_{0}} + {{\Pi }_{0}}) + F({{{\bar {u}}}_{0}} + {{Q}_{0}}). \\ \end{gathered} $

Теорема 4.1. Пусть функция $F(u)$ имеет вид $F(u) = {{u}^{3}} - {{b}^{3}}$, где $b > 0$, и граничное значение $\varphi (0,0) > {{\bar {u}}_{0}}(0,0)$. Тогда задача (4.1), (4.2) имеет решение ${{P}_{0}}(\xi ,\tau )$, удовлетворяющее экспоненциальной оценке убывания вида (3.5).

Доказательство. Считаем, что ${{\bar {u}}_{0}} = b$. С помощью метода верхних и нижних решений покажем, что существует решение задачи (4.1), (4.2), заключенное между барьерными функциями следующего вида:

${{Z}_{ - }}(\xi ,\tau ) = - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} \leqslant {{P}_{0}}(\xi ,\tau ) \leqslant 0 = {{Z}_{ + }}(\xi ,\tau ).$

Для доказательства того, что функция ${{Z}_{ + }}(\xi ,\tau ) = 0$ является верхним решением задачи, нужно проверить выполнение условия (4.4), которое для рассматриваемой функции имеет вид

${{(b + y + z)}^{2}} = {{(b + y)}^{2}} = {{(b + z)}^{2}},$
так как ${{\bar {u}}_{0}}(0,0) = b$. Это условие выполняется в единственной точке $y = z = 0$, которая не является внутренней для области. Поэтому точек экстремума нет. Условие (4.5), очевидно, выполняется. Таким образом, функция ${{Z}_{ + }}(\xi ,\tau ) = 0$ является верхним решением задачи.

Для доказательства того, что функция ${{Z}_{ - }}(\xi ,\tau ) = - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} $ является нижним решением задачи, нужно проверить выполнение неравенства $L({{Z}_{ - }}) \geqslant 0$. С этой целью запишем (4.6) с учетом замены ${{\bar {u}}_{0}} \to b$, ${{\Pi }_{0}} \to y$, ${{Q}_{0}} \to z$, $L({{Z}_{ + }}) \to L$ в виде

(4.7)
$L = \sqrt {yz} \left( {\frac{1}{{{{z}^{2}}}}\int\limits_0^z \,F(b + u)du - \frac{1}{z}F(b + z) - \frac{1}{y}F(b + y)} \right) - F(b + y + z - 2\sqrt {yz} ) + F(b + y) + F(b + z).$

Доказательство неравенства $L \geqslant 0$ значительно упрощается, если отбросить первое слагаемое и рассмотреть более сильное неравенство

(4.8)
$ - \sqrt {yz} \left( {\frac{1}{z}F(b + z) + \frac{1}{y}F(b + y)} \right) - F(b + y + z - 2\sqrt {yz} ) + F(b + y) + F(b + z) \geqslant 0.$

В отличие от $L$ левая часть (4.8) является функцией, симметричной относительно $y$ и $z$. Ее можно привести к зависимости от переменных

(4.9)
$s = \sqrt {yz} \in (0,\phi - b],\quad t = \frac{{y + z}}{2} \in (0,\phi - b].$

Имеем

$F(b + u) = {{(b + u)}^{3}} - {{b}^{3}} = u[{{(b + u)}^{2}} + b(b + u) + {{b}^{2}}] = u({{u}^{2}} + 3bu + 3{{b}^{2}}),$
$\begin{gathered} F(b + y) + F(b + z) = ({{y}^{3}} + 3b{{y}^{2}} + 3{{b}^{2}}y) + ({{z}^{3}} + 3b{{z}^{2}} + 3{{b}^{2}}z) = {{(y + z)}^{3}} - \\ - \;3yz(y + z) + 3b{{(y + z)}^{2}} - 6byz + 3{{b}^{2}}(y + z) = 8{{t}^{3}} - 6{{s}^{2}}t + 12b{{t}^{2}} - 6b{{s}^{2}} + 6{{b}^{2}}t, \\ \end{gathered} $
$\begin{gathered} \frac{1}{z}F(b + z) + \frac{1}{y}F(b + y) = ({{z}^{2}} + 3bz + 3{{b}^{2}}) + ({{y}^{2}} + 3by + 3{{b}^{2}}) = \\ \, = {{(y + z)}^{2}} - 2yz + 3b(y + z) + 6{{b}^{2}} = 4{{t}^{2}} - 2{{s}^{2}} + 6bt + 6{{b}^{2}}, \\ \end{gathered} $
$F(b + y + z - 2\sqrt {yz} ) = F(b + 2t - 2s).$

В результате замены (4.9) неравенство (4.8) примет вид

$ - s(4{{t}^{2}} - 2{{s}^{2}} + 6bt + 6{{b}^{2}}) - F(b + 2t - 2s) + 8{{t}^{3}} - 6{{s}^{2}}t + 12b{{t}^{2}} - 6b{{s}^{2}} + 6{{b}^{2}}t \geqslant 0,$
или

(4.10)
$H: = 8{{t}^{3}} - 6{{s}^{2}}t + 12b{{t}^{2}} - 6b{{s}^{2}} + 6{{b}^{2}}t - 4s{{t}^{2}} + 2{{s}^{3}} - 6bst - 6{{b}^{2}}s - F(b + 2t - 2s) \geqslant 0.$

Найдем образ квадрата $G:(0,\phi - b] \times (0,\phi - b]$ при отображении

$(y,z) \to \left( {\sqrt {yz} ,\frac{{y + z}}{2}} \right).$

Так как

$\frac{{y + z}}{2} \geqslant \sqrt {yz} ,$
то оба треугольника $y < z$ и $z < y$, составляющие квадрат $G$, отображаются в угол $t > s$. Две стороны квадрата $G$:
$y = 0,\quad z \in (0,\phi - b]\quad {\text{и}}\quad z = 0,\quad y \in (0,\phi - b],$
переходят в один и тот же отрезок

(4.11)
$s = 0,\quad t \in \left( {0,\frac{{\phi - b}}{2}} \right].$

Диагональ квадрата $G$: $z = y$, $y \in (0,\phi - b]$, переходит в отрезок биссектрисы

(4.12)
$t = s,\quad s \in (0,\phi - b].$

И, наконец, оставшиеся две стороны квадрата $G$:

$y = \phi - b,\quad z \in (0,\phi - b]\quad {\text{и}}\quad z = \phi - b,\quad y \in (0,\phi - b],$
переходят в один и тот же кусок параболы

(4.13)
$t = \frac{{{{s}^{2}}}}{{2(\phi - b)}} + \frac{{\phi - b}}{2},\quad s \in \left( {0,\frac{{\phi - b}}{2}} \right].$

Итак, образом квадрата $G$ является криволинейный треугольник $\Gamma $, ограниченный линиями (4.11)–(4.13) и накрываемый дважды.

Исследуем функцию $H$, определяемую в (4.10), на экстремумы в треугольнике $\Gamma $. Необходимые условия экстремума

$\frac{{\partial H}}{{\partial s}} = - 12st - 12bs - 4{{t}^{2}} + 6{{s}^{2}} - 6bt - 6{{b}^{2}} + 2F{\kern 1pt} '(b + 2t - 2s) = 0,$
$\frac{{\partial H}}{{\partial t}} = 24{{t}^{2}} - 6{{s}^{2}} + 24bt + 6{{b}^{2}} - 8st - 6bs - 2F{\kern 1pt} '(b + 2t - 2s) = 0$
приводят к соотношениям
$F{\kern 1pt} '(b + 2t - 2s) = 6st + 6bs + 2{{t}^{2}} - 3{{s}^{2}} + 3bt + 3{{b}^{2}} = 12{{t}^{2}} - 3{{s}^{2}} + 12bt + 3{{b}^{2}} - 4st - 3bs,$
которые выполняются только при $t = s$. Поэтому внутри треугольника $\Gamma $ точек экстремума нет. Остается проверить выполнение неравенства (4.10) на границе треугольника $\Gamma $.

1. На стороне

$s = 0,\quad t \in \left( {0,\frac{{\phi - b}}{2}} \right]$
значения

$H(0,t) = 8{{t}^{3}} + 12b{{t}^{2}} + 6{{b}^{2}}t - F(b + 2t) = F(b + 2t) - F(b + 2t) = 0 \geqslant 0.$

2. На стороне

$t = s,\quad s \in (0,\phi - b]$
значения

$H(s,s) = 8{{s}^{3}} - 6{{s}^{3}} + 12b{{s}^{2}} - 6b{{s}^{2}} + 6{{b}^{2}}s - 4{{s}^{3}} + 2{{s}^{3}} - 6b{{s}^{2}} - 6{{b}^{2}}s - F(b) = 0 \geqslant 0.$

3. Для проверки неравенства (4.10) на стороне

$t = \frac{{{{s}^{2}}}}{{2(\phi - b)}} + \frac{{\phi - b}}{2},\quad s \in \left( {0,\frac{{\phi - b}}{2}} \right]$
удобнее вернуться к переменным $(y,z)$ и рассмотреть одну из двух сторон квадрата $G$:

$y = \phi - b,\quad z \in (0,\varphi - b]\quad {\text{или}}\quad z = \phi - b,\quad y \in (0,\phi - b].$

Для определенности будем рассматривать сторону $z = \phi - b,$ $y \in (0,\phi - b]$. Нужно доказать неравенство (4.8) при $z = \phi - b$. Запишем левую часть этого неравенства в виде

$\left( {1 - \sqrt {\frac{y}{z}} } \right)F(b + z) + \left( {1 - \sqrt {\frac{z}{y}} } \right)F(b + y) - F((b + {{(\sqrt z - \sqrt y )}^{2}}) \geqslant 0,$
и сделаем замены

(4.14)
$t = \sqrt {\frac{y}{z}} ,\quad 0 < t < 1,\quad a = z,\quad a > 0.$

В результате получим неравенство

$(1 - t)F(b + a) - \frac{{1 - t}}{t}F(b + a{{t}^{2}}) - F(b + a{{(1 - t)}^{2}}) \geqslant 0.$

Здесь

$F(b + u) = uh(u),\quad {\text{где}}\quad h(u) = {{u}^{2}} + 3bu + 3{{b}^{2}}.$

Соответственно имеем неравенство

$(1 - t)ah(a) - (1 - t)ath(a{{t}^{2}}) - a{{(1 - t)}^{2}}h(a{{(1 - t)}^{2}}) \geqslant 0,$
которое после деления на $(1 - t)a$ принимает вид

(4.15)
$h(a) - th(a{{t}^{2}}) - (1 - t)h(a{{(1 - t)}^{2}}) \geqslant 0.$

Преобразуем левую часть этого неравенства:

$\begin{gathered} h(a) - th(a{{t}^{2}}) - (1 - t)h(a{{(1 - t)}^{2}}) = ({{a}^{2}} + 3ba + 3{{b}^{2}}) - t({{a}^{2}}{{t}^{4}} + 3ba{{t}^{2}} + 3{{b}^{2}}) - \\ \, - [{{a}^{2}}{{(1 - t)}^{4}} + 3ba{{(1 - t)}^{2}} + 3{{b}^{2}}] + t[{{a}^{2}}{{(1 - t)}^{4}} + 3ba{{(1 - t)}^{2}} + 3{{b}^{2}}] = \\ = a[a + 3b - a{{(1 - t)}^{4}} - 3b{{(1 - t)}^{2}} + a(4{{t}^{4}} - 6{{t}^{3}} + 4{{t}^{2}} - t) + 3b(2{{t}^{2}} - t)]. \\ \end{gathered} $

Поэтому неравенство (4.15) эквивалентно неравенству

(4.16)
$g: = a + 3b - a{{(1 - t)}^{4}} - 3b{{(1 - t)}^{2}} + a(4{{t}^{4}} - 6{{t}^{3}} + 4{{t}^{2}} - t) + 3b(2{{t}^{2}} - t) \geqslant 0.$

Изучим зависимость значений $g$ от параметра $a$ на промежутке $a > 0$. Производная

$\frac{{\partial g}}{{\partial a}} = 1 - {{(1 - t)}^{4}} + 4{{t}^{4}} - 6{{t}^{3}} + 4{{t}^{2}} - t = t(3{{t}^{3}} - 2{{t}^{2}} - 2t + 3) > 0,$
поэтому $g$ возрастает на промежутке $a > 0$ и ее величина

$g(a) > g(0) = 3b - 3b{{(1 - t)}^{2}} + 3b(2{{t}^{2}} - t) = 3b({{t}^{2}} + t) > 0.$

Таким образом, (4.16), а вместе с ним и (4.8) при $z = \phi - b$, являются верными неравенствами. Это завершает доказательство неравенства (4.10), и можно утверждать, что функция ${{Z}_{ - }}(\xi ,\tau ) = - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} $ является нижним решением задачи (4.1), (4.2).

Обе барьерные функции ${{Z}_{ - }}(\xi ,\tau )$ и ${{Z}_{ + }}(\xi ,\tau )$ удовлетворяют экспоненциальным оценкам убывания вида (3.5), поэтому решение ${{P}_{0}}(\xi ,\tau )$ задачи (4.1), (4.2) также удовлетворяeт оценке того же вида. Теорема доказана.

Теорема 4.2. Пусть функция $F(u)$ имеет вид $F(u) = {{u}^{3}} - {{b}^{3}}$, где $b > 0$, и граничное значение $\phi (0,0) > {{\bar {u}}_{0}}(0,0)$. Тогда задачи (3.3), (3.4) имеют решения ${{P}_{k}}(\xi ,\tau )$, удовлетворяющие экспоненциальным оценкам убывания вида (3.5).

Доказательство. Исследуем знак производной $F{\kern 1pt} '({{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} + {{P}_{0}})$. Так как

${{P}_{0}} \geqslant - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} ,$
то

${{\Pi }_{0}} + {{Q}_{0}} + {{P}_{0}} \geqslant {{\Pi }_{0}} + {{Q}_{0}} - 2\sqrt {{{\Pi }_{0}}{{Q}_{0}}} = \mathop {(\sqrt {{{\Pi }_{0}}} - \sqrt {{{Q}_{0}}} )}\nolimits^2 \geqslant 0.$

Поэтому величина

${{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} + {{P}_{0}} = b + {{\Pi }_{0}} + {{Q}_{0}} + {{P}_{0}} \geqslant b,$
вследствие чего производная $F{\kern 1pt} '({{\bar {u}}_{0}} + {{\Pi }_{0}} + {{Q}_{0}} + {{P}_{0}}) > 0$ и задачи (3.3), (3.4) имеют решения с оценками вида (3.5). Теорема доказана.

Отметим, что функции $P_{k}^{*}({{\xi }_{*}},\tau ),\;k \geqslant 0$, определяются аналогично.

5. ОЦЕНКА ОСТАТОЧНОГО ЧЛЕНА

Итак, ряд (1.3) оказывается полностью построенным при дополнительном условии 6:

Условие 6. В угловых точках $(k,0)$, $k = 0,1$, прямоугольника $\Omega $ функция $F(u)$ имеет вид $F(u) = {{u}^{3}} - {{b}^{3}}$, где положительное число $b$ не обязательно одно и то же, и граничные значения $\phi (k,0) > {{\bar {u}}_{0}}(k,0)$.

Теорема 5.1. Пусть выполнены условия 1–6. Тогда для достаточно малых $\varepsilon $ задача (0.1)–(0.3) имеет решение $u(x,t,\varepsilon )$, для которого ряд

$\sum\limits_{k = 0}^\infty \,{{\varepsilon }^{k}}({{\bar {u}}_{k}}(x,t) + {{\Pi }_{k}}(x,\tau ) + {{Q}_{k}}(\xi ,t) + Q_{k}^{*}({{\xi }_{*}},t) + {{P}_{k}}(\xi ,\tau ) + P_{k}^{*}({{\xi }_{*}},\tau ))$
является асимптотическим представлением при $\varepsilon \to 0$ в замкнутом прямоугольнике $\bar {\Omega }$.

Доказательство основано на разрешимости задач для пограничных функций ${{\Pi }_{k}}$, ${{Q}_{k}}$, $Q_{k}^{*}$, ${{P}_{k}}$ и $P_{k}^{*}$ при $k \geqslant 1$ и повторяет доказательство соответствующей теоремы из статьи [3].

Выражаю искреннюю благодарность В.Ф. Бутузову за плодотворное обсуждение полученных результатов.

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

  1. Денисов И.В. Первая краевая задача для квазилинейного сингулярно возмущенного параболического уравнения в прямоугольнике // Ж. вычисл. матем. и матем. физ. 1996. Т. 36. № 10. С. 56–72.

  2. Денисов И.В. Оценка остаточного члена в асимптотике решения краевой задачи // Ж. вычисл. матем. и матем. физ. 1996. Т. 36. № 12. С. 64–67.

  3. Денисов И.В. Угловой пограничный слой в краевых задачах для сингулярно возмущенных параболических уравнений с квадратичной нелинейностью // Ж. вычисл. матем. и матем. физ. 2017. Т. 57. № 2. С. 255–274.

  4. Денисов А.И., Денисов И.В. Угловой пограничный слой в краевых задачах для сингулярно возмущенных параболических уравнений с нелинейностями // Ж. вычисл. матем. и матем. физ. 2019. Т. 59. № 1. С. 102–117.

  5. Денисов А.И., Денисов И.В. Угловой пограничный слой в краевых задачах для сингулярно возмущенных параболических уравнений с немонотонными нелинейностями // Ж. вычисл. матем. и матем. физ. 2019. Т. 59. № 9. С. 1581–1590.

  6. Денисов И.В. Угловой пограничный слой в краевых задачах для сингулярно возмущенных параболических уравнений с монотонной нелинейностью // Ж. вычисл. матем. и матем. физ. 2018. Т. 58. № 4. С. 1–11.

  7. Васильева А.Б., Бутузов В.Ф. Асимптотические методы в теории сингулярных возмущений. М.: Высш. школа, 1990.

  8. Amann H. On the existence of positive solutions of nonlinear elliptic boundary value problems // Indiana Univ. Math. J. 1971. V. 21. № 2. P. 125–146.

  9. Sattinger D.H. Monotone methods in nonlinear elliptic and parabolic boundary value problems // Indiana Univ. Math. J. 1972. V. 21. № 11. P. 979–1000.

  10. Amann H. // Nonlinear Analysis: coll. of papers in honor of E.H. Rothe / Ed. by L. Cesari et al. New York etc: Acad press, cop. 1978. XIII. P. 1–29.

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