Журнал вычислительной математики и математической физики, 2019, T. 59, № 2, стр. 264-276
Оценки в классах Гёльдера решения неоднородной задачи Дирихле для сингулярно возмущенного однородного уравнения конвекции-диффузии
В. Б. Андреев *, И. Г. Белухина **
МГУ ВМК
119992 Москва, Ленинские Горы, Россия
* E-mail: andreev@cs.msu.su
** E-mail: belukh@cs.msu.su
Поступила в редакцию 25.03.2018
После доработки 03.09.2018
Аннотация
В полуплоскости рассматривается неоднородная первая краевая задача для однородного сингулярно возмущенного уравнения конвекции-диффузии с постоянными коэффициентами и конвекцией, направленной ортогонально границе полуплоскости с направлением от границы. В предположении принадлежности граничной функции пространству ${{C}^{{2,\lambda }}}$, $0 < \lambda < 1$, получена неулучшаемая оценка ограниченного на бесконечности решения в соответствующей норме Гёльдера. Библ. 5.
1. ВВЕДЕНИЕ
Целью данной работы является получение ${{C}^{{2,\lambda }}}$-оценки решения заданного в полуплоскости однородного сингулярно возмущенного уравнения конвекции-диффузии с постоянными коэффициентами и направленной ортогонально от границы конвекцией через соответствующую норму правой части граничного условия Дирихле. Эта работа является продолжением исследований, начатых в [1], где рассматривалась та же задача, но для неоднородного уравнения с однородным граничным условием. Объединение этих результатов приводит к оценке полностью неоднородной задачи, что открывает короткий путь к получению оценок в ${{C}^{{k,\lambda }}}$, где $k > 2$. Более того, точность полученных оценок позволит их использовать для исследования уравнения с переменными коэффициентами, как это было сделано в [2] для одномерного случая (ср. с [3, гл. III]). Следует отметить, что ранее близкие оценки были получены в [4], но там оценивались коэффициенты Гёльдера только по направлению, задаваемому границей полуплоскости. Для получения этих оценок использовался аппарат преобразования Фурье по касательной переменной и соответствующая теорема о мультипликаторе. Эта красивая техника, к сожалению, не позволяет получить оценки коэффициентов Гёльдера по нормальному к границе направлению и тем самым не позволяет получить наш результат. С другой стороны, оценки из [4] не являются достаточными для исследования уравнения с переменными коэффициентами. Мы же для получения оценок используем аппарат функции Грина, который является более громоздким, однако позволяющим конструктивно получить оценки коэффициентов Гёльдера по любому направлению, т.е. с выписыванием при необходимости постоянных в этих оценках. Наша техника получения оценок предполагает раздельное оценивание коэффициентов Гёльдера по переменным $x$ и $y$. В связи с тем, что оценки по касательному направлению (у нас это направление $y$) уже получены в [4], мы в работе приводим доказательство оценок коэффициентов Гёльдера только по направлению $x$, а для оценок по $y$ формулируем результат и ссылаемся на [4].
Дальнейшее содержание работы таково. В разд. 2 дается постановка задачи и формулируется основной результат работы – теорема 1. В разд. 3 оцениваются коэффициенты Гёльдера по направлению $x$ самого решения и его первых производных. В разд. 4 оцениваются вторые производные.
На протяжении всей статьи для обозначения постоянных, которые не зависят от решения и малого параметра $\varepsilon $, будем использовать строчные буквы c с индексом или без. Одной и той же буквой часто будем обозначать различные постоянные.
2. ПОСТАНОВКА ЗАДАЧИ И ОСНОВНОЙ РЕЗУЛЬТАТ
В правой полуплоскости $\mathbb{R}_{ + }^{2}$ плоскости $OXY$ ищется ограниченное решение задачи
(2.1)
$ - \varepsilon \Delta U + 2\alpha \frac{{\partial U}}{{\partial X}} + qU = 0,\quad (X,Y) \in \mathbb{R}_{ + }^{2},$Анализировать задачу (2.1), (2.2) удобнее в растянутых переменных
Пусть $U(X,Y) = U(\varepsilon x,\varepsilon y) = :u(x,y)$, а $\mathcal{G}(Y) = \mathcal{G}(\varepsilon y) = :g(y)$. В новых переменных задача (2.1), (2.2) примет вид(2.4)
$Lu: = - \Delta u + 2\alpha \frac{{\partial u}}{{\partial x}} + \varepsilon qu = 0,\quad (x,y) \in \mathbb{R}_{ + }^{2},$(2.6)
$u(x,y) = - \frac{1}{\pi }{{e}^{{\alpha x}}}\int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}\left( {a\sqrt {{{x}^{2}} + {{{(\eta - y)}}^{2}}} } \right)}}{{\partial x}}g(\eta )d\eta ,$Чтобы сформулировать основной результат работы, введем некоторые обозначения. Будем через ${{C}^{{k,\lambda }}}$ обозначать пространство $k$ раз непрерывно дифференцируемых в $\mathbb{R}_{ + }^{2}$ функций, $k$-е производные которых удовлетворяют условию Гёльдера с показателем $\lambda $. Если какой-либо индекс $k$ или $\lambda $ будет равен нулю, то будем его опускать и писать просто ${{C}^{\lambda }}$, ${{C}^{k}}$ или $C$. Под полунормой $|\, \cdot \,|$ в пространстве ${{C}^{k}}$ будем понимать сумму максимумов модулей ее $k$-х производных, а в ${{C}^{{k,\lambda }}}$ – сумму коэффициентов Гёльдера этих производных, т.е.
Теорема 1. Для решения $U(X,Y)$ задачи (1), (2) справедливы априорные оценки
(2.8)
${{\left\| U \right\|}_{{{{C}^{{k,\lambda }}}}}} \leqslant c{{\left\| G \right\|}_{{{{C}^{{k,\lambda }}}}}},\quad k = 0,1,2,\quad \lambda \in (0,1).$Доказательству этой теоремы посвящены остальные разделы работы.
3. ОЦЕНКИ РЕШЕНИЯ И МЛАДШИХ ПРОИЗВОДНЫХ
Начнем со вспомогательного утверждения.
Лемма 3.1. Пусть положительные числа $a$ и $\alpha $ связаны соотношением (2.7), а $r > 0$, $0 \leqslant \varphi \leqslant \pi {\text{/}}2$. Тогда при $\delta \geqslant 0$, $\beta \geqslant 0$ справедлива оценка
(3.1)
$f(r,\varphi ) = \mathop {\left( {sin\varphi } \right)}\nolimits^\delta {{r}^{\beta }}{{e}^{{ - \left( {a - \alpha cos\varphi } \right)r}}} \leqslant c\left( {1 + {{\varepsilon }^{{\delta /2 - \beta }}}} \right).$Доказательство. Рассмотрим оцениваемую функцию $f(r,\varphi )$ как функцию, зависящую от $r$, считая $\varphi $ параметром. При $\beta = 0$ очевидно максимум достигается при $r = 0$, поэтому
При $\beta > 0$ функция $f(r,\varphi )$ неотрицательна, обращается в нуль при $r = 0$, экспоненциально убывает к нулю при $r \to \infty $ и поэтому имеет максимум по $r$ в некоторой точке $r{\text{*}}(\varphi )$. Этот максимум, зависящий от $\varphi $ как от параметра, естьПолезной во многих случаях в дальнейшем будет следующая формула:
(3.2)
$\left( {\frac{\partial }{{\partial x}} + a} \right)\frac{{\partial {{K}_{0}}(ar)}}{{\partial x}} = \frac{{a\cos 2\varphi }}{r}{{K}_{1}}(ar) - {{a}^{2}}\cos \varphi \left[ {{{K}_{1}}(ar) - {{K}_{0}}(ar)} \right] - {{a}^{2}}\cos \varphi (1 - \cos \varphi ){{K}_{0}}(ar),$Перейдем к оценкам решения и его первых производных.
Лемма 3.2. Для решения $u(x,y)$ задачи (2.4), (2.5) справедлива оценка
(3.3)
${{\left| u \right|}_{{C_{x}^{\lambda }}}} \leqslant \frac{c}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{\lambda }} + {{\varepsilon }^{\lambda }}{{{\left| g \right|}}_{0}}} \right\}.$Доказательство. Сначала несколько преобразуем представление (2.6) решения $u(x,y)$ задачи (2.4), (2.5), вычтя в подынтегральном выражении из функции $g(\eta )$ постоянное значение $g(y)$ и прибавив его снова. Получим
(3.4)
$u(x,y) = - \frac{1}{\pi }\int\limits_{ - \infty }^\infty {{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}}}{{\partial x}}\left[ {g(\eta ) - g(y)} \right]d\eta - \frac{1}{\pi }{{e}^{{\alpha x}}}g(y)\int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}}}{{\partial x}}d\eta .$Последнее слагаемое теперь можно вычислить, учитывая четность по $(\eta - y)$ подынтегральной функции. После замены переменной интегрирования $\eta {\text{'}} = \eta - y$ и использования формулы ($\nu = 1$, $2\mu + 1 = 0$) последнее слагаемое в (3.4) принимает вид
(3.5)
$ - \frac{1}{\pi }{{e}^{{\alpha x}}}g(y)\int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}}}{{\partial x}}d\eta = {{e}^{{ - (a - \alpha )x}}}g(y).$После замены в (3.4) переменной $\eta {\kern 1pt} {\text{'}} = \eta - y$ при $\eta \geqslant y$ и $\eta {\kern 1pt} {\text{''}} = - \eta + y$ при $\eta \leqslant y$, опуская в дальнейшем штрихи у $\eta $, с учетом (3.5) и четности по $\eta $ функции $\partial {{K}_{0}}{\text{/}}\partial x$, приходим к представлению
(3.6)
$u(x,y) = - \frac{1}{\pi }\int\limits_0^\infty {{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}\left( {a\sqrt {{{x}^{2}} + {{\eta }^{2}}} } \right)}}{{\partial x}}\delta _{\eta }^{2}g(y)d\eta + g(y){{e}^{{ - (a - \alpha )x}}},$Очевидно, что для оценки функции или ее производных в ${{C}^{\lambda }}$ достаточно оценить их величины в $C_{x}^{\lambda }$ и $C_{y}^{\lambda }$. Проведем оценки в $C_{x}^{\lambda }$. Пусть $\bar {x}$, $x$ – две различные точки положительной полуоси $Ox$. Для определенности будем полагать, что $\bar {x} < x$, так что $x - \bar {x} = \Delta > 0$. Оценим разность значений решения $u(x,y)$ из (3.6) в точках $x$ и $\bar {x}$
(3.8)
$\begin{gathered} u(x,y) - u(\bar {x},y) = - \frac{1}{\pi }\int\limits_0^\infty \left[ {{{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}\left( {a\sqrt {{{x}^{2}} + {{\eta }^{2}}} } \right)}}{{\partial x}} - {{e}^{{\alpha \bar {x}}}}\frac{{\partial {{K}_{0}}\left( {a\sqrt {\mathop {\bar {x}}\nolimits^2 + {{\eta }^{2}}} } \right)}}{{\partial{ \bar {x}}}}} \right]\delta _{\eta }^{2}g(y)d\eta + \\ + \;g(y)\left[ {{{e}^{{ - (a - \alpha )x}}} - {{e}^{{ - (a - \alpha )\bar {x}}}}} \right] = :{{V}_{0}}(x,\bar {x},y) + {{V}_{{00}}}(x,\bar {x},y). \\ \end{gathered} $(3.9)
${{V}_{{00}}}(x,\bar {x},y) = g(y)\left[ {{{e}^{{ - (a - \alpha )x}}} - {{e}^{{ - (a - \alpha )\bar {x}}}}} \right].$(3.10)
$\left| {{{V}_{{00}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{\lambda }{{\Delta }^{\lambda }}{{(a - \alpha )}^{\lambda }}{{\left| g \right|}_{0}}.$Теперь рассмотрим оставшуюся часть (3.8), ${{V}_{0}}(x,\bar {x},y)$. Для ее оценки, следуя [3, гл. III], разобьем интервал интегрирования $[0,\infty )$ на два интервала: $[0,\Delta ]$ и $[\Delta ,\infty )$. Далее, в интеграле по $[\Delta ,\infty )$ разность выражений в квадратных скобках представим в виде интеграла по отрезку $[\bar {x},x]$ от соответствующей производной. Получим
где(3.12)
${{V}_{{01}}}(x,\bar {x},y) = - \frac{1}{\pi }\int\limits_0^\Delta \left[ {{{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}}}{{\partial x}} - {{e}^{{\alpha \bar {x}}}}\frac{{\partial {{K}_{0}}}}{{\partial{ \bar {x}}}}} \right]\delta _{\eta }^{2}g(y)d\eta ,$(3.14)
$\left| {\delta _{\eta }^{2}g(y)} \right| \leqslant 2{{\eta }^{\lambda }}{{\left| g \right|}_{\lambda }}.$(3.15)
$\left| {{{V}_{{01}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{\lambda }{{\Delta }^{\lambda }}{{\left| g \right|}_{\lambda }}.$(3.16)
$\left| {{{V}_{{022}}}(\bar {x},x,y)} \right| \leqslant c\frac{{{{\Delta }^{\lambda }}}}{{(1 - \lambda )}}{{\left| g \right|}_{\lambda }}.$При оценке ${{V}_{{021}}}(x,\bar {x},y)$, поступая аналогично, получаем
Далее выполним интегрирование по $\eta $ и используем представление $\left( {\tfrac{\partial }{{\partial z}} + a} \right)\tfrac{{\partial {{K}_{0}}}}{{\partial z}}$ в полярных координатах формулой (3.2), а затем оценим величины $\left[ {{{K}_{1}}(r) - {{K}_{0}}(r)} \right]$, ${{K}_{1}}(r)$, ${{K}_{0}}(r)$ при помощи формул из [1, ф. (А.6), (А.9), (А.10)]. После применения к каждому из слагаемых вида ${{(sin\varphi )}^{\delta }}{{r}^{\beta }}{{e}^{{ - (a - \alpha cos\varphi )r}}}$ в полученном выражении леммы 3.1 при ($\delta = 2$, $\beta = 0$), ($\delta = 2,$, $\beta = 1{\text{/}}2$), ($\delta = 4$, $\beta = 3{\text{/}}2$) соответственно будем иметь(3.17)
$\left| {{{V}_{{021}}}(\bar {x},x,y)} \right| \leqslant c{{\Delta }^{\lambda }}\frac{1}{{(1 - \lambda )}}{{\left| g \right|}_{\lambda }}.$Лемма 3.3. Для производной $\partial u{\text{/}}\partial x$ решения задачи (2.4), (2.5) справедлива оценка
(3.18)
${{\left| {\frac{{\partial u}}{{\partial x}}} \right|}_{{C_{x}^{\lambda }}}} \leqslant c\frac{1}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{{1,\lambda }}} + {{\varepsilon }^{{1 + \lambda }}}{{{\left| g \right|}}_{0}}} \right\}.$Доказательство. Используя представление (2.6) решения $u(x,y)$ задачи (2.4), (2.5), запишем выражение $\partial u{\text{/}}\partial x$ в виде
(3.19)
$\frac{{\partial u}}{{\partial x}}(x,y) = - \frac{1}{\pi }{{e}^{{\alpha x}}}\left( {\frac{\partial }{{\partial x}} + \alpha } \right)\int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}}}{{\partial x}}g(\eta )d\eta .$(3.22)
$\begin{gathered} {{I}_{{10}}}(x,y) = - \frac{1}{\pi }{{e}^{{\alpha x}}}\left( {\frac{\partial }{{\partial x}} + \alpha } \right)\int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}}}{{\partial x}}\left[ {S_{g}^{1}(\eta - y)} \right]d\eta , \\ {{I}_{{11}}}(x,y) = - \frac{1}{\pi }\int\limits_{ - \infty }^\infty {{e}^{{\alpha x}}}\left( {\frac{\partial }{{\partial x}} + \alpha } \right)\frac{{\partial {{K}_{0}}}}{{\partial x}}\left[ {g(\eta ) - S_{g}^{1}(\eta - y)} \right]d\eta , \\ \end{gathered} $Сделаем в интегралах (3.22) замену переменной $\eta {\text{'}} = \eta - y$. Заметим, что величина ${{I}_{{10}}}(x,y)$ с учетом (3.20) имеет вид
Теперь преобразуем ${{I}_{{11}}}(x,y)$. Для этого воспользуемся уравнением для ${{K}_{0}}\left( {a\sqrt {{{x}^{2}} + {{\eta }^{2}}} } \right)$: $ - \Delta {{K}_{0}} + {{a}^{2}}{{K}_{0}} = 0$. Тогда ${{I}_{{11}}}(x,y)$ примет вид
(3.24)
${{I}_{{11}}}(x,y) = - \frac{1}{\pi }\int\limits_{ - \infty }^\infty {{e}^{{\alpha x}}}\left[ { - \frac{{{{\partial }^{2}}{{K}_{0}}}}{{\partial {{\eta }^{2}}}} + {{a}^{2}}{{K}_{0}} + \alpha \frac{{\partial {{K}_{0}}}}{{\partial x}}} \right]\left[ {g(\eta + y) - S_{g}^{1}(\eta )} \right]d\eta .$(3.25)
${{I}_{{11}}}(x,\bar {x},y) = - \frac{1}{\pi }\int\limits_0^\infty \left\{ {{{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}}}{{\partial \eta }}\left[ {g{\text{'}}(y + \eta ) - g{\text{'}}(y - \eta )} \right] + \left( {\alpha \frac{{\partial {{K}_{0}}}}{{\partial x}} + {{a}^{2}}{{K}_{0}}} \right)\delta _{\eta }^{2}g(y)} \right\}d\eta .$Теперь можно переходить к оценке разности $\partial u{\text{/}}\partial x - \partial u{\text{/}}\partial{ \bar {x}}$. С учетом (3.21) имеем
(3.26)
$\left| {{{V}_{{10}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{\lambda }{{\Delta }^{\lambda }}{{(a - \alpha )}^{{1 + \lambda }}}{{\left| g \right|}_{0}}.$Теперь рассмотрим ${{V}_{{11}}}(x,\bar {x},y)$ (см. (3.25)). Снова, следуя [3, гл. III], представим интеграл по $\eta $ в виде суммы интегралов ${{V}_{{111}}}(x,\bar {x},y)$ и ${{V}_{{112}}}(x,\bar {x},y)$ по интервалам соответственно $[0,\Delta ]$ и $[\Delta ,\infty )$, затем в интеграле по интервалу $[\Delta ,\infty )$, как и при рассмотрении ${{\left| u \right|}_{{C_{x}^{\lambda }}}}$, представим разность входящих туда выражений, зависящих от $x$ и $\bar {x}$, в виде интеграла по отрезку $[\bar {x},x]$ от производной по соответствующей переменной. Получим
(3.27)
$\begin{gathered} {{V}_{{111}}}(x,\bar {x},y) = - \frac{1}{\pi }\int\limits_0^\Delta \left\{ {\left[ {{{e}^{{\alpha x}}}\frac{{\partial {{K}_{0}}}}{{\partial \eta }} - {{e}^{{\alpha \bar {x}}}}\frac{{\partial {{K}_{0}}}}{{\partial \eta }}} \right]\left[ {g{\text{'}}(y + \eta ) - g{\text{'}}(y - \eta )} \right] + } \right. \\ + \;\left. {\left[ {{{e}^{{\alpha x}}}\left( {\alpha \frac{{\partial {{K}_{0}}}}{{\partial x}} + {{a}^{2}}{{K}_{0}}} \right) - {{e}^{{\alpha \bar {x}}}}\left( {\alpha \frac{{\partial {{K}_{0}}}}{{\partial{ \bar {x}}}} + {{a}^{2}}{{K}_{0}}} \right)} \right]\delta _{\eta }^{2}g(y)} \right\}d\eta {\kern 1pt} {\kern 1pt} , \\ \end{gathered} $(3.28)
$\begin{gathered} {{V}_{{112}}}(x,\bar {x},y) = - \frac{1}{\pi }\int\limits_\Delta ^\infty d\eta \int\limits_{\bar {x}}^x {{e}^{{\alpha z}}}\left\{ {\left[ {\left( {\frac{\partial }{{\partial z}} + \alpha } \right)\frac{{\partial {{K}_{0}}}}{{\partial \eta }}} \right]\left[ {g{\text{'}}(y + \eta ) - g{\text{'}}(y - \eta )} \right]} \right. + \\ \left. { + \;\left( {\frac{\partial }{{\partial z}} + \alpha } \right)\left( {\alpha \frac{{\partial {{K}_{0}}}}{{\partial z}} + {{a}^{2}}{{K}_{0}}} \right)\delta _{\eta }^{2}g(y)} \right\}dz. \\ \end{gathered} $При рассмотрении же интеграла в (3.28) (по $[\Delta ,\infty )$) заменим внутренний интеграл по теореме о среднем. Получим для ${{V}_{{112}}}(x,\bar {x},y)$ оценку
Дальнейший процесс доказательства схематично выглядит следующим образом. Выполним интегрирование по $\eta $ в (3.29), затем найдем максимумы входящих величин. Для этого оценим максимум модуля суммы суммой максимумов модулей отдельных выражений, предварительно преобразованных к удобному для исследования виду, перейдем в них от декартовых координат к полярным, а затем воспользуемся полученными в [1] представлениями и оценками этих выражений и применим лемму 3.1 к каждому слагаемому вида ${{(sin\varphi )}^{\delta }}{{r}^{\beta }}{{e}^{{ - (a - \alpha cos\varphi )r}}}$.Более подробно. Воспользовавшись формулами из [1, ф. (А.3), (А.9)] и леммой 3.1 при ($\delta = 2$, $\beta = 0$), ($\delta = 2$, $\beta = 1{\text{/}}2$) в первом слагаемом из (3.29), а затем, после преобразования второго слагаемого из той же формулы к виду
(3.31)
$\left| {{{V}_{{111}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{\lambda }{{\Delta }^{\lambda }}{{\left| g \right|}_{{1,\lambda }}}.$Аналогично поступаем и при оценке $\left| {{{V}_{{112}}}(x,\bar {x})} \right|$ из (3.30), т.е. очевидным образом преобразуем оцениваемое выражение к виду
(3.32)
$\left| {{{V}_{{112}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{{(1 - \lambda )}}{{\Delta }^{\lambda }}{{\left| g \right|}_{{1,\lambda }}}.$(3.33)
$\left| {{{V}_{{11}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{{\lambda (1 - \lambda )}}{{\Delta }^{\lambda }}{{\left| g \right|}_{{1,\lambda }}}.$Лемма 3.4. Для производной $\partial u{\text{/}}\partial y$ решения задачи (2.4), (2.5) справедлива оценка
(3.34)
${{\left| {\frac{{\partial u}}{{\partial y}}} \right|}_{{C_{x}^{\lambda }}}} \leqslant c\frac{1}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{{1,\lambda }}} + {{\varepsilon }^{\lambda }}{{{\left| g \right|}}_{1}}} \right\}.$Доказательство следует из оценки (3.3), если в последней заменить $g(y)$ на $g{\text{'}}(y)$, так как производная по $y$ совпадает с решением той же задачи (2.4), (2.5), но с правой частью $g{\text{'}}(y)$ граничного условия (2.5) вместо $g(y)$.
4. ОЦЕНКИ ВТОРЫХ ПРОИЗВОДНЫХ
Продолжим изучение задачи (2.4), (2.5). В этом разделе будут получены необходимые оценки вторых производных в $C_{x}^{\lambda }$.
Лемма 4.1. Для производной ${{\partial }^{2}}u{\text{/}}\partial {{x}^{2}}$ решения задачи (2.4), (2.5) справедлива оценка
(4.1)
${{\left| {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right|}_{{C_{x}^{\lambda }}}} \leqslant c\frac{1}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{{2,\lambda }}} + {{\varepsilon }^{\lambda }}{{{\left| g \right|}}_{2}} + {{\varepsilon }^{{2 + \lambda }}}{{{\left| g \right|}}_{0}}} \right\}.$Доказательство. Для того, чтобы выразить ${{\partial }^{2}}u{\text{/}}\partial {{x}^{2}}$, снова используем представление (2.6) решения $u(x,y)$ задачи (2.4), (2.5). Имеем
(4.3)
$\begin{gathered} \frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}(x,y) = {{I}_{{20}}}(x,y) + {{I}_{{21}}}(x,y),\quad {\text{г д е }} \\ {{I}_{{20}}}(x,y) = - \frac{1}{\pi }{{e}^{{\alpha x}}}\mathop {\left( {\frac{\partial }{{\partial x}} + \alpha } \right)}\nolimits^2 \int\limits_{ - \infty }^\infty \frac{{\partial {{K}_{0}}}}{{\partial x}}S_{g}^{2}(\eta - y)d\eta , \\ {{I}_{{21}}}(x,y) = - \frac{1}{\pi }\int\limits_{ - \infty }^\infty {{e}^{{\alpha x}}}\mathop {\left( {\frac{\partial }{{\partial x}} + \alpha } \right)}\nolimits^2 \frac{{\partial {{K}_{0}}}}{{\partial x}}\left[ {g(\eta ) - S_{g}^{2}(\eta - y)} \right]d\eta , \\ \end{gathered} $Сделаем в (4.3) замену переменной интегрирования $\eta {\text{'}} = \eta - y$. Заметим, что с учетом (4.2) выражение ${{I}_{{20}}}(x,y)$ примет вид
(4.4)
$\frac{1}{\pi }\int\limits_{ - \infty }^\infty \frac{{{{\eta }^{2}}}}{2}\frac{{\partial {{K}_{0}}\left( {a\sqrt {{{x}^{2}} + {{\eta }^{2}}} } \right)}}{{\partial x}}d\eta = - \frac{1}{{2a}}x{{e}^{{ - ax}}}.$(4.5)
${{I}_{{20}}}(x,y) = {{(a - \alpha )}^{2}}{{e}^{{ - (a - \alpha )x}}}g(y) + \frac{1}{{2a}}\left[ {{{{(a - \alpha )}}^{2}}x - 2(a - \alpha )} \right]{{e}^{{ - (a - \alpha )x}}}g{\text{''}}(y).$(4.6)
$\begin{gathered} {{I}_{{21}}}(x,y) = \frac{1}{\pi }\int\limits_{ - \infty }^\infty {{e}^{{\alpha x}}}\left\{ {\left[ {\frac{{\partial {{K}_{0}}}}{{\partial x}} + 2\alpha {{K}_{0}}} \right]\left[ {g(\eta + y) - S_{g}^{2}(\eta )} \right]{\text{''}} - } \right. \\ - \;\left. {\left[ {{{{(a - \alpha )}}^{2}}\frac{{\partial {{K}_{0}}}}{{\partial x}} + 2a\alpha \left( {\frac{\partial }{{\partial x}} + a} \right){{K}_{0}}} \right]\left[ {g(\eta + y) - S_{g}^{2}(\eta )]} \right]} \right\}d\eta . \\ \end{gathered} $(4.7)
${{I}_{{21}}}(x,y) = \frac{1}{\pi }\int\limits_0^\infty {{e}^{{\alpha x}}}\left\{ {{{\mathfrak{F}}_{1}}(x,\eta )\delta _{\eta }^{2}g{\text{''}}(y) + {{\mathfrak{F}}_{2}}(x,\eta )\left[ {\delta _{\eta }^{2}g(y) - {{\eta }^{2}}g{\text{''}}(y)} \right]} \right\}d\eta ,$Теперь можно переходить к оценке разности
(4.9)
$\left| {{{V}_{{20}}}(x,\bar {x},y)} \right| \leqslant \frac{c}{\lambda }{{\Delta }^{\lambda }}\left\{ {{{{(a - \alpha )}}^{{2 + \lambda }}}{{{\left| g \right|}}_{0}} + {{{(a - \alpha )}}^{{1 + \lambda }}}{{{\left| g \right|}}_{2}}} \right\}.$Рассмотрим теперь ${{V}_{{21}}}(x,\bar {x},y)$ (см. (4.7)). Снова, следуя [3, гл. III], представим интеграл по $\eta $ в виде суммы интегралов по интервалам $[0,\Delta ]$ и $[\Delta ,\infty )$, затем в интеграле по интервалу $[\Delta ,\infty )$, представим разность входящих туда выражений, зависящих от $x$ и $\bar {x}$, в виде интеграла по отрезку $[\bar {x},x]$ от производной по соответствующей переменной. Будем иметь
(4.10)
$\begin{gathered} {{V}_{{21}}}(x,\bar {x},y) = {{V}_{{211}}}(x,\bar {x},y) + {{V}_{{212}}}(x,\bar {x},y), \\ {{V}_{{211}}}(x,\bar {x},y) = \frac{1}{\pi }\int\limits_0^\Delta \left\{ {\left[ {{{e}^{{\alpha x}}}{{\mathfrak{F}}_{1}}(x,\eta ) - {{e}^{{\alpha \bar {x}}}}{{\mathfrak{F}}_{1}}(\bar {x},\eta )} \right]\delta _{\eta }^{2}g{\text{''}}(y) + } \right. \\ + \;\left. {\left[ {{{e}^{{\alpha x}}}{{\mathfrak{F}}_{2}}(x,\eta ) - {{e}^{{\alpha \bar {x}}}}{{\mathfrak{F}}_{2}}(\bar {x},\eta )} \right]\left[ {\delta _{\eta }^{2}g(y) - {{\eta }^{2}}g{\text{''}}(y)} \right]} \right\}d\eta , \\ {{V}_{{212}}}(x,\bar {x},y) = \frac{1}{\pi }\int\limits_\Delta ^\infty d\eta \int\limits_{\bar {x}}^x {{e}^{{\alpha z}}}\left( {\frac{\partial }{{\partial z}} + \alpha } \right)\left\{ {{{\mathfrak{F}}_{1}}(z,\eta )\delta _{\eta }^{2}g{\text{''}}(y) + {{\mathfrak{F}}_{2}}(z,\eta )\left[ {\delta _{\eta }^{2}g(y) - {{\eta }^{2}}g{\text{''}}(y)} \right]} \right\}dz. \\ \end{gathered} $Точнее, при рассмотрении величин (1), (2) и (3) используем формулы из [1, ф. (А.9), (А.10)], содержащие оценки ${{K}_{0}}$ и ее производной в полярных координатах, что приводит к необходимости оценить выражения
Рассматривая величину (4), используем формулу из [1, ф. (В.5)], выражающую оцениваемую функцию в полярных координатах, а затем применяем оценки из [1, ф. (А.6), (А.10)]. Это приводит нас к необходимости оценить величину
Для оценки величины (5), после очевидного преобразования $\alpha = a - (a - \alpha )$, используем формулу (3.2), выражающую оцениваемую величину в полярных координатах, а затем применяем оценки из [1, ф. (А.3), (А.9), (А.6), (А.10)]. В результате приходим к необходимости оценить выражение
Аналогично поступаем при рассмотрении величины (6). Здесь мы также после очевидного преобразования $\alpha = a - (a - \alpha )$, использования формулы из [1, ф. (В.5)] для представления оцениваемой величины в полярных координатах и оценок из [1, ф. (А.6), (А.10)] сводим первоначальную оценку к оценке величины выражения
Исследование величины (7) аналогично уже проведенному исследованию величины (5), что приводит к необходимости оценить выражение
И, наконец, поступая аналогично при оценке выражения (8), т.е. в очередной раз заменяем $\alpha = a - (a - \alpha )$, используем формулы из [1, ф. (В.6), (В.5)], оценки из [1, ф. (А.7), (А.9), (А.10)] и применяем лемму 3.1 при соответствующих каждому слагаемому вида $si{{n}^{\delta }}\varphi {{r}^{\beta }}{{e}^{{ - (a - \alpha cos\varphi )r}}}$ значениях $\delta $, $\beta $, приходим к оценке выражения (8) величиной $\left[ {{{c}_{{11}}} + {{c}_{{12}}}(a - \alpha ){{\varepsilon }^{{ - 1/2}}}} \right]$.
Таким образом, после выполнения в (4.11) интегрирования по $\eta $, с учетом (2.7), получим
(4.12)
$\left| {{{V}_{{21}}}(x,\bar {x},y)} \right| \leqslant c\frac{1}{{\lambda (1 - \lambda )}}{{\Delta }^{\lambda }}{{\left| g \right|}_{{2,\lambda }}}.$Лемма 4.2. Для производных ${{\partial }^{2}}u{\text{/}}\partial {{y}^{2}}$ и ${{\partial }^{2}}u{\text{/}}\partial x\partial y$ решения задачи (2.4), (2.5) справедливы оценки
(4.13)
${{\left| {\frac{{{{\partial }^{2}}u}}{{\partial x\partial y}}} \right|}_{{C_{x}^{\lambda }}}} \leqslant c\frac{1}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{{2,\lambda }}} + {{\varepsilon }^{{1 + \lambda }}}{{{\left| g \right|}}_{1}}} \right\}.$(4.14)
${{\left| {\frac{{{{\partial }^{2}}u}}{{\partial {{y}^{2}}}}} \right|}_{{C_{x}^{\lambda }}}} \leqslant c\frac{1}{\lambda }\left\{ {\frac{1}{{(1 - \lambda )}}{{{\left| g \right|}}_{{2,\lambda }}} + {{\varepsilon }^{\lambda }}{{{\left| g \right|}}_{2}}} \right\}.$Доказательство первой из этих оценок следует из оценки (3.18), если в последней заменить $g(y)$ на $g{\text{'}}(y)$, так как смешанная производная совпадает с производной по $x$ решения той же задачи (2.4), (2.5), но с правой частью $g{\text{'}}(y)$ граничного условия (2.5) вместо $g(y)$. Аналогично вторая оценка следует из (3.3) при замене в ней $g(y)$ на $g{\text{''}}(y)$.
5. ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 1
Суммируем оценки, полученные в разд. 1–4 (см. леммы 3.2–4.2).
Лемма 5.1. Для решения $U(X,Y)$ задачи (2.1), (2.2) справедливы оценки
(5.1)
${{\left\| U \right\|}_{{C_{x}^{{k,\lambda }}}}} \leqslant c{{\left\| \mathcal{G} \right\|}_{{{{C}^{{k,\lambda }}}}}},\quad k = 0,\;1,\;2.$Доказательство. Перепишем требующие доказательства неравенства (5.1) более подробно в нерастянутых переменных.
Имеем: при $k = 1$
(5.2)
${{\left| {\frac{{\partial u}}{{\partial x}}} \right|}_{C}} \leqslant \frac{{{{t}^{\lambda }}}}{{(1 + \lambda )}}{{\left| {\frac{{\partial u}}{{\partial x}}} \right|}_{{C_{x}^{\lambda }}}} + 2{{t}^{{ - 1}}}{{\left\| u \right\|}_{C}},\quad t \in (0,\infty )\; - \;{\text{л ю б о е ,}}$(5.3)
${{\left| u \right|}_{{C_{x}^{\lambda }}}} \leqslant {{2}^{{1 - \lambda }}}\left[ {\frac{\lambda }{{(1 + \lambda )}}{{\varepsilon }^{{ - 1}}}{{{\left| {\frac{{\partial u}}{{\partial x}}} \right|}}_{{C_{x}^{\lambda }}}} + (1 - \lambda )(1 + \lambda ){{\varepsilon }^{\lambda }}{{{\left\| u \right\|}}_{C}}} \right].$(5.4)
${{\left| g \right|}_{{{{C}^{\lambda }}}}} \leqslant {{2}^{{1 - \lambda }}}\left[ {\frac{\lambda }{{1 + \lambda }}{{\varepsilon }^{{ - 1}}}{{{\left| g \right|}}_{{{{C}^{{1,\lambda }}}}}} + (1 - \lambda )(1 + \lambda ){{\varepsilon }^{\lambda }}{{{\left| g \right|}}_{C}}} \right],\quad {{\left| g \right|}_{{{{C}^{2}}}}} \leqslant \frac{{{{\varepsilon }^{{ - \lambda }}}}}{{1 + \lambda }}{{\left| g \right|}_{{{{C}^{{2,\lambda }}}}}} + 2\varepsilon {{\left\| g \right\|}_{{{{C}^{1}}}}}.$Замечание 1. Оценки решения $U(X,Y)$ задачи (2.1), (2.2) в $C_{y}^{{k,\lambda }}$ ($k = 0,\;1,\;2$) получаются аналогично оценкам в $C_{x}^{{k,\lambda }}$, но мы их использовать не будем.
Справедливость утверждений (2.8) теоремы 1 следует теперь из (5.1) и полученных в [4] оценок решения в $C_{y}^{{k,\lambda }}$ при $k = 0,\;1,\;2$. Тем самым теорема 1 полностью доказана.
Список литературы
Андреев В.Б. Оценки в классах Гёльдера регулярной составляющей решения сингулярно возмущенного уравнения конвекции-диффузии // Ж. вычисл. матем. и матем. физ. 2017. Т. 57. № 12. С. 1983–2020.
Андреев В.Б. К оценке гладкости регулярной составляющей решения одномерного сингулярно возмущенного уравнения конвекции-диффузии // Ж. вычисл. матем. и матем. физ. 2015. Т. 55. № 1. С. 22–33.
Ладыженская О.А., Уральцева Н.Н. Линейные и квазилинейные уравнения эллиптического типа. М.: Наука, 1971.
Naughton A., Kellogg R.B., Stynes M. Regularity and derivative bounds for a convection-diffusion Problem with a Neumann outflow condition // J. Differential Equations V. 247. 2009. P. 2495–2516.
Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1971.
Дополнительные материалы отсутствуют.
Инструменты
Журнал вычислительной математики и математической физики