Опорные барьерные функции для нелинейных параболических задач Текст научной статьи по специальности «Математика»

ЧЕБЫШЕВСКИИ СБОРНИК

Том 25. Выпуск 2.

УДК 517.9

DOI 10.22405/2226-8383-2024-25-2-235-242

Опорные барьерные функции для нелинейных параболических

задач

А. И. Денисов, И. В. Денисов

Денисов Алексей Игоревич — аспирант, Тульский государственный педагогический университет им. Л. Н. Толстого (г. Тула). e-mail: den_tspu@mail.ru

Денисов Игорь Васильевич — доктор физико-математических наук, профессор, Тульский государственный педагогический университет им. Л. Н. Толстого (г. Тула). e-mail: den_tspu@mail.ru

В рамках нелинейного метода угловых пограничных функций существование решений нелинейных краевых задач доказывается через построение барьерных функций. Барьерные функции конструируются через выделенные специальным образом опорные барьеры. Сами опорные барьеры также могут выступать в роли барьерных функций. Получаемые неравенства в свою очередь представляют самостоятельный функциональный интерес.

Ключевые слова: нелинейные краевые задачи, барьерные функции, функциональные неравенства.

Библиография: 7 названий. Для цитирования:

А. И. Денисов, И. В. Денисов. Опорные барьерные функции для нелинейных параболических задач // Чебышевский сборник, 2024, т. 25, вып. 2, с. 235-242.

The support barrier functions for nonlinear parabolic problems

A. I. Denisov, I. V. Denisov

Denisov Alexey Igorevich — postgraduate student, Tula State Lev Tolstoy Pedagogical University (Tula).

e-mail: den_tspu@mail.ru

Denisov Igor Vasil'evich — doctor of physical and mathematical sciences, professor, Tula State Lev Tolstoy Pedagogical University (Tula). e-mail: den_tspu@mail.ru

Аннотация

CHEBYSHEVSKII SBORNIK Vol. 25. No. 2.

UDC 517.9

DOI 10.22405/2226-8383-2024-25-2-235-242

Abstract

Within the framework of the nonlinear method of angular boundary functions, the existence of solutions to nonlinear boundary value problems is proven through the construction of barrier functions. Barrier functions are constructed through specially designated support barriers. The support barriers themselves can also act as barrier functions. The resulting inequalities, in turn, are of independent functional interest.

Keywords: nonlinear boundary value problems, barrier functions, functional inequalities.

Bibliography: 7 titles.

For citation:

A. I. Denisov, I. V. Denisov, 2024. "The support barrier functions for nonlinear parabolic problems" ,

Chebyshevskii sbornik, vol. 25, no. 2, pp. 235-242.

1. Постановка задачи

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

(д2 и ди\

- = р(и,х,1,£), (х,1) е О,

и(х, 0, е) = ^(х), 0 ^ х ^ 1,

и(о,г,е)= ^(ъ), и(1,г,е)= ^(ъ), о ^ г ^ т,

где О = {(х,Ь) | 0 < х < 1, 0 < Ь < Т} - прямоугольник. Предположим, что выполнены следующие условия.

Условие 1. Функции Р(и,х,Ь,е), ф(х), ф\(Ь) и ф2(Ъ) являются достаточно гладкими и в угловых точках прямоугольника, О выполняются условия согласованност,и на,ч,а,л,ьно-кра,евы,х значений

<р(0) = ф1(0), <р(1)= М0).

Условие 2. Вырожденное уравнение Р(и, х,1, 0) = 0 в замкнутом прямоугольнике О имеет решение, которое обозначается как и = и0(х,Ь).

Условие 3. Производная, ^(и0 (х,Ь),х,Ь, 0) > 0 в замкнутом прямоугольнике О. Условие 4. Начальная, задача,

НП

—0 = -Р(щ(х, 0) + По,х, 0,0), По(х, 0) = у(х) - щ(х, 0), ат

имеет решение П0(х,т) при т ^ 0 удовлетворяющее условию По (ж, те) = 0 (здесь параметр х е [0,1]).

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

2 йх2 ,, . , .

—— = г2, а —— = Ь (и0(кА) + г\,кЛ, 0), ау ау

прямыег\ = ф\+к(^ — и0(к, Ь) пересекают, сепаратрисы, входящие в точку покоя (г\, г2) = (0, 0) при у ^ те (здесь £ - параметр, к = 0 или 1).

Применим нелинейный метод угловых пограничных функций, в рамках которого, как и в линейном случае, решение задачи ищется в виде асимптотического ряда по параметру е ^ 0:

и(х,г,е)= й + (П + + <2*) + (Р + Р *).

Здесь й- регулярная часть асимптотики, играющая роль внутри прямоугольника О, П, и погранслойные функции, играющие роль вблизи сторон прямоугольника О соответственно Ь = 0 х = 0и х = 1, Ри Р *- угловые пограничные функции, играющие роль вблизи вершин прямоугольника О соответственно (0, 0) и (1, 0).

Формальная процедура построения регулярной части асимптотики и погранслойных функций хорошо отработана (см. [1]). Однако для понимания используемых обозначений необходимы минимальные пояснения алгоритма, более подробно см. [2]. Каждая часть асимптотики строится в виде ряда по степеням е. Регулярная часть:

те

к ,

и(х,Ь,е) = ^ ек ии (х,Ь), к=0

погранслойные функции:

тете те

П(х,т,е) = ^ е* Щ (х,т), Я(и,е) = ^ ек Qfc ({,$, <2*(£*,1,е) = ^ ек Я*к (£*,£),

к=0 к=0 к=0 угловые пограничные функции:

р а, т,е) = ^ ^ р* а, г), р * а*, г, е) = ^ ^ ть, г), к=0 к=0

где

^ 1 Ж t

= -, ?* =-, т =

е е е2

- растянутые переменные.

Наиболее сложной в алгоритме является нелинейная задача для нахождения главного члена угловой части асимптотики. Для Р0(^,т) задача ставится в первой четверти плоскости растянутых переменных (£,т):

2

и имеет вид

П := {(М | £> 0, т> 0},

Ь(Р0) = 0 в области М+, (1)

Р0(0,т) = -П0(0,г), Р0(£, 0) = -^(С, 0), (2)

Р0((,т) ^ 0 при £ + т ^ ж, (3)

где

Р(2) := а2- — - Р (ад + П0 + Я0 + 2) + Р (ад + П0) + Р (ад + ^0) (4)

>д2г д2_

~д.ё2 - ~дг

и Р(и) = Р(и, 0, 0, 0). Для определенности будем считать, что в угловой точке (0, 0) граничное значение р больше корня вырожденного уравнения ад. Кроме этого, будем рассматривать случай, когда функция Р(и) положительна на промежутке [ад,^>].

Для доказательства существования решения задачи (1)-(4) используется метод верхних и нижних решений (см. [3] - [5]), который заключается в том, что задача

Р(2) = 0 в области Р,

2 = к на границе дР

имеет решение 2 в границах

2- < 2 < 2+,

если в области И выполняются неравенства

ь(г+) ^ 0, ь(г-) ^ 0, г- ^ г+,

а на ее границе

г- ^ и ^ г+.

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

Определение. Барьерную функцию Z(£, т) назовем опорной, если ее использование в качестве барьера задачи (1)-(4) приводит к неравенству, линейному относительно функции ^ (и).

Опорные барьерные функции впервые появились в работах [6], [7], однако отдельно не изучались. К настоящему моменту удалось выделить три опорных барьерных функции. Рассмотрим эти функции, из которых, при необходимости, конструируются более сложные барьеры.

2. Первая опорная барьерная функция

Первая опорная барьерная функция - нулевая:

= 0.

Для этой функции на границе области М+ выполняются неравенства

гг(0,т) = 0 > -По(0,т), г^, 0) = 0 > -Яо(С, 0).

Поэтому может претендовать только на роль верхнего барьера. Внутри области М+ для верхнего барьера должно выполняться неравенство Ь^1) < 0. Имеем

ВД) = (йо + По + Яо) + ^ (щ + По) + ^ (щ + Яо). (5)

Для краткости обозначим

По(0,т) = з, <2о(£, 0)= I. (6)

В результате выражение (5) примет вид

ь(г1) = (ио + в + г) + ^ (зд + 8) + ^ (йо + г),

п можно сформулировать следующее утверждение. Лемма 1. Функция

= 0

является верхним, барьером для решения Ро = Ро(£,т) задачи (1)-(4), если при любых значениях в и £ из промежутка (0,<р - ио] величина,

р(0) = -р (ио + 8 + г) + р (ио + 8) + р (ио + г) < 0.

(7)

3. Вторая опорная барьерная функция

Вторая опорная барьерная функция имеет вид

И0^0

¿2 =--Г-.

Р - ад

Такой вид продиктован соотношением

По + Qo - ^ = (--So) (l - (l - )(

- -So V V - -So/ v

io + Qo--— = (- - Son 1 - 11--— 111--) ) , (8)

- - So V V - -So/V - -So J J

из которого следует, что величина

So + no + Qo + Z2 e (So; -]. На границе области R+ имеем

Z2(0, т) = -no(0, т), Z2(C, 0) = -Qo(£, 0),

Z2(C, т) — 0 при £ + т —У ж>.

Поэтому Z2 может претендовать на роль как верхнего, так и нижнего барьера. Внутри области R+ имеем

Т(У ) 2 no d2Qo + Qo dno F {_ +n +Q noQo ^ +

L(Z2) = -a2-z——^ +--z—---F So + no + Qo--— +

- - so d£2 - -sodr \ - - so /

+F (So + no) + F (So + Qo) = no -F(So + Qo)--Q^f(So + no) - F (So + П + Qo - +

р -ад Р - S0 V р -110,

(ад + П0) + ^ (ад + ^0) =

= (1 - Р(ад + ^0) + (1 - Р(ггс + П0) - ^ (ад + П0 + ^ - =

V <р - ад / - «0 / V V so)

= (1---— (ад + *) + ( 1---—("^0 + 8) -Р(ад + « + £--— )

- ЗД / \ Р - ЗД / V V - «0 /

при обозначениях (6). Таким образом, верно следующее утверждение. Лемма 2. Функция

7 И0^0

¿2 =--—

Р -ад

является барьером для, решения Р0 = Р0({, г) задачи (1)-(4), если при любых значениях в иЛ из промежутка (0, р - ад] выражение

ВД) = (1--— ) ^(ад + *) + (1--— ) Р(ад + «) - ^ ( ад + 8 + 4--— ) (9)

\ р -щ/ V <Р - и0 / V ¥ - и0 /

сохраняет свой знак.

4. Третья опорная барьерная функция

Третья опорная барьерная функция определяется с помощью погранслойных функций По = По(0, г) и Яо = Яо(С, 0) как

^э = -2^ПоЯо. Выбор такого вида продиктован соотношением

По + Яо - 2^ПоЯо = (^П - л/Оо)2 , (10)

которое показывает, что величина

и0 + По + Яо + е [и0; р),

так как По, фо е (0; <р - йо].

Задача для определения функции По = По(0, т) имеет вид

аПо = ^ (щ+По), По(0,0) = ¡р - и0

йт

(см. [2]). Поэтому

СШо = (щ + По),

(1т

и решение По(0,т) представляет собой функцию, убывающую от По(0, 0) = р - йо > 0 до По(0, те) = 0.

Задача для определения функции фо(С, 0) имеет вид

2 с12Яо

2

(также см. [2]). Поэтому

й^2 а2

и, понижая порядок уравнения, получаем

'ЩоУ_ 2 [

= ^ (и0 + Яо) , Яо(0, 0) = <р - и0, Яо(те, 0)=0, (РЯо = \ Р (щ + Яо),

(Що\ 2

Ы) = р+и)аи.

<!£ ) а- Jо Выбираем случай отрицательной производной:

йЯо 2

2 [Чо

=--2 ^ (ио + и)йи,

а Уо

а

и решение Яо(С, 0^, которое является функцией, убывающей от фо(0,0) = р - щ > 0 до Яо(те, 0) = 0.

На границе области М+ имеем

гэ(0, т) = -2^-йо)П~о < -По,

гэ(£, 0) = -2^(<р - ио)Яо < -Яо, т) ^ 0 при £ + т ^ те,

поэтому Zэ может претендовать только на роль нижнего барьера.

Внутри области М+ имеем

Ь (-2^и0о0) = -2а2Ти0 - 2Уо0 ^-

-Р (ад + П0 + ^0 - 2^Л000) + ^ (ад + П0) + Р (ад + ^0). (11)

Здесь

-2о2^Й0 ^ = -2а2| () = | (f) =

_ _2„2 П (_ 1 (^ V + ^ \ =

= 2а уП0 ^ V # ) +2У00 ^2 )

_ /12 11 N

= -2 а2(-4д0Ж-22 1 р(Й0 + * + 27Ш^(й0 + д0)) =

_ уп -

Р (и0 + и)йи--■= Р (и0 + Qo),

.)0

-2,/0~о= -—1^ = (ад + П0).

йт л/Щ йт у/Пё

Собираем все вместе и с учетом замены (6) получаем

ВД) = — [ Р(й0 + и)йи - —Р({¿0 + ¿) - —\р ({¿0 + в) -70 Vt V 8

-Р (ад + в + г - 2—^) + Р("¿¿0 + з) + Р("¿¿0 + *) =

= - —^ ^ ("¿¿0 + *) + - —0 Р ("¿¿0 + 8) - Р (ад + в + г - 2—^) + — ^ р("¿¿0 +

Таким образом, верно следующее утверждение. Лемма 3. Функция

^э = -2^Й000

является нижним барьером для решения Р0 = Р0({, г) задачи (1)-(4), если при любых значениях в и £ из промежутка (0, р - ад] величина,

Ь Ш = (1 -—(ад + *) + (1 -—("¿¿0 + «) --Р (ад + 8 + г - + — ! ^("¿¿0 + и)йи > 0. (12)

5. Сравнение опорных барьерных функций

При любых значениях 8 и £ из промежутка (0, р - ад] справедливо неравенство

2—^--—— = —* (2 ) > 0.

Ф - «0 \ ^ - "¿¿0 /

Поэтому

-2—^ < - **

^ - «0

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

Лемма 4. Опорные барьерные функции находятся, между собой в следующем, порядке:

< <

а

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

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

2. Денисов А.И., Денисов И.В. Нелинейный метод угловых пограничных функций для сингулярно возмущенных параболических задач с кубическими нелинейностями // Чебышев-ский сборник, 2024, т. 25, вып. 1, с. 25-40.

3. Amann Н. On the Existence of Positive Solutions of Nonlinear Elliptic Boundary Value Problems // Indiana Univ. Math. J. 1971. Vol.21, № 2. P. 125 - 146.

4. Sattinger D.H. Monotone Methods in Nonlinear Elliptic and Parabolic Boundary Value Problems // Indiana Univ. Math. J. 1972. V. 21. № 11. P. 979 - 1000.

5. 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.

6. Денисов И.В. Об асимптотическом разложении решения сингулярно возмущенного эллиптического уравнения в прямоугольнике // Асимптотические методы теории сингулярно возмущенных уравнений и некорректно поставленных задач: Сб. научн. тр. - Бишкек: II.iii.M. 1991. - С. 37.

7. Денисов И.В. Квазилинейные сингулярно возмущенные эллиптические уравнения в прямоугольнике // Журнал вычисл. матем. и матем. физ. - Т.35. № 11. 1995. - С. 1666-1678. (English transl.: Denisov I.V. Quasilinear Singularly Perturbed Elliptic Equations in a Rectangle // Computational Mathematics and Mathematical Physics. 1995. Vol. 35. № 11. pp. 1341-1350.)

REFERENCES

1. Vasilveva, А. В., Butuzov, V. F. 1990. "Asymptotic methods in the theory of singular perturbations", M.: Higher school.

2. Denisov A.I., Denisov I.V., 2024. "Nonlinear method of angular boundary functions for singularly perturbed parabolic problems with cubic nonlinearities" , Chebyshevskii Sbornik, Vol. 25, № . 1, pp. 25-40.

3. Amann H., 1971. "On the Existence of Positive Solutions of Nonlinear Elliptic Boundary Value Problems", Indiana Univ. Math. J., Vol.21, № 2. pp. 125-146.

4. Sattinger D.H., 1972. "Monotone Methods in Nonlinear Elliptic and Parabolic Boundary Value Problems", Indiana Univ. Math. J., Vol.21. № 11. pp. 979-1000.

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

6. Denisov I.V., 1991. "On the asymptotic expansion of the solution of a singularly perturbed elliptic equation in a rectangle", Asymptotic m,et,hods of the theory of singularly perturbed equations and ill-posed problems: Collection of articles, scientific, tr. - Bishkek: Him, p. 37.

7. Denisov I.V., 1995. "Quasilinear Singularly Perturbed Elliptic Equations in a Rectangle", Computational Mathematics and Mathematical Physics, Vol. 35. № 11. pp. 1341-1350.

Получено: 06.03.2024 Принято в печать: 28.06.2024