Челябинский физико-математический журнал. 2023. Т. 8, вып. 1. С. 18-28.
УДК 517.956 Б01: 10.47475/2500-0101-2023-18102
ЗАДАЧА БЕЗ НАЧАЛЬНЫХ УСЛОВИЙ
ДЛЯ УРАВНЕНИЯ С ДРОБНЫМИ ПРОИЗВОДНЫМИ
И ПРОМЕЖУТОЧНЫЕ АСИМПТОТИКИ
В. А. Костин1", Д. В. Костин1'2'6, Х. Алкади1
1 Воронежский государственный университет, Воронеж, Россия 2Воронежский государственный педагогический университет, Воронеж, Россия "[email protected], [email protected]
Показана корректная разрешимость задач без начальных условий для дробно-степенных операторных сумм. Решения задач без начальных условий Я. Б. Зельдович и Г. И. Баренблатт трактуют как промежуточные асимптотики для задач с начальными условиями. На важность таких задач эти авторы указывают в связи с расширением понятия «строгого детерминизма» в статистической физике и квантовой механике и ставят вопрос об изучении свойств явлений, не зависящих от деталей в начальных условиях, проявляющихся при истечении достаточного времени. В данной работе также приводится пример промежуточной асимптотики для уравнения с дробной производной.
Ключевые слова: промежуточная асимптотика, корректная задача, задача Коши, уравнение без начальных условий, сильно непрерывная полугруппа, дробная степень оператора.
Введение
Задачи без начальных условий относятся к классу задач, описывающих установившиеся периодические или переходные процессы, начавшиеся так давно, что начальные данные практически не оказывают влияния на поведение решения в момент наблюдения. В частности, такую задачу решают А. Н. Тихонов, А. А. Самарский в [1] при изучении распространения тепловых волн в почве. На важность таких задач указывают Г. И. Баренблатт и Я.Б.Зельдович в [2], утверждая, что внимание к их теории связано с расширением понятия «строгого детерминизма», обязанным успехам статистической физики и квантовой механики, где этот вопрос рассматривается иначе, нежели в классической физике и механике. Это привело к постановке вопроса о тех свойствах явления, которые не зависят от деталей в начальных условиях. Но такая независимость возможна лишь по истечении достаточного времени. То есть в таком случае можно говорить об асимптотике решений по времени, стремящемуся к бесконечности. Поэтому в [1] такие асимптотики называются «промежуточными», а соответствующие задачи без начальных условий — «вырожденными».
Однако для того, чтобы решение «вырожденной задачи» являлось промежуточной асимптотикой, необходимо, чтобы оно было устойчивым относительно возмущений, то есть такая задача должна быть корректной в некотором классе функций.
Работа выполнена при финансовой поддержке Министерства науки и высшего образования РФ в рамках выполнения государственного задания в сфере науки (номер темы FZGF-2020-0009).
Установление этого факта является ключевым, и здесь используются методы функционального анализа [3-5].
Например, для £ Е (-га, го) единственным ограниченным решением уравнения
+ Аи(£) = 1, Л > 0, (1)
является и(£) = 1/А, а для £ Е [0, го) единственное решение уравнения (1) с условием и(0) = 0 есть функция и0(£) = (1 — е-Л4)/А. Таким образом, и(£) = 1/А — промежуточная асимптотика функции и0 (£): и0(£) ^ и(£) при £ ^ го.
Следующим важным направлением применения задач без начальных условий стало исследование методами дробного интегро-дифференцирования процессов в пористых и фрактальных средах [6-8] (см. также фундаментальный обзор В.В.Учайкина [9]).
Исследованию задач без начальных условий для уравнений с дробной производной посвящена заметка авторов [10], где основное внимание они уделяют применению сильно непрерывных полугрупп и косинус-функций к исследованию задач для уравнений с переменными коэффициентами, не касаясь такого важного понятия, как промежуточная асимптотика. Рассмотрению вопросов, связанных с этим понятием, посвящается настоящая работа.
1. Постановка задачи
Пусть Е — банахово пространство с нормой || ■ ||е = || ■ || и А — линейный замкнутый оператор в Е с областью определения Д(А), плотной в Е.
Обозначим через См,(Ед) банахово пространство вектор-функций /(£), £ Е (—го, го) = К, со значениями в Е и нормой
||/|и = «пр ||е-^/(£)||, М 2 0.
¡ев
Введём в рассмотрение дробно-дифференциальное выражение
п л -ат
Ь„(а,й)и(£) = ати(£),
т= 1
где а = (а1, а2 ..., ап), ат 2 0, <5 = (а1, а2,..., ап), ат Е (0,1), производная Маршо имеет вид
-ат 1 гте
- и(£) = —-, в-(1+а) [и(£ — в) — и(в)]-в,
-£ат Г( — «тЫ 0
Г(г) — гамма-функция Эйлера.
Для / Е См,(Ед) рассмотрим уравнение
Ьп(а,й)и(£) = Аи(£) + / (£). (2)
Решением уравнения (2) будем называть вектор-функцию
и Е с^ед^ (3)
для которой и(£) Е ^(А) при всех £ Е К, на К определены дробные производные
<1Ьат
Маршо ^¡ати и выполняется равенство (2).
Будем говорить, что задача (2), (3) поставлена корректно, если она имеет единственное решение u(t) и для него выполняется оценка корректности
IMIc. ^ M||/||Cm, (4)
где константа M ^ 0 не зависит от /.
Функцию u+ (t) будем называть решением задачи (2), (3) с начальным условием
u(0) = 0, (5)
если она удовлетворяет уравнению (2), включению (3) и условию (5).
Решение u(t) задачи (2), (3) будем называть промежуточной асимптотикой решения u+(t) задачи (2)-(5), если выполняется соотношение
lim ||u(t) - u+(t)Hß = 0. (6)
В работе доказывается корректность задачи (2), (3), откуда следует, что её решение является промежуточной асимптотикой решения задачи (2)-(5). Приводится пример такой асимптотики.
2. Связь дробных степеней операторов и дробных производных
Связь дробных степеней операторов в банаховом пространстве с дробными производными описана в [7, с. 102] следующим образом. Если оператор —А является генератором полугруппы и (г, — А) класса С0 с оценкой
II и (г, —А)р\\ ^ Ме-^|М1, т> 0, г ^ 0, (7)
то для А определены его дробные степени формулой Балакришнана:
1 [ ^
А> = г-1-а[и(г, —А) — а е (0,1). (8)
г(—а) Л
Например, если А задаётся выражением , х е го), и областью опре-
деления О (А) = е Ьр: (р' е Ьр,ш}, где Ьр,ад — весовое лебегово банахово пространство с нормой
p,w
e-wx\,(x)\ pdx
i/p
w > 0, p ^ 1,
то —A является генератором С0-полугруппы U(t, —A),(x) = ,(x — t), t ^ 0, и для , G D(A) равенство
A., = D++, = ( dX )a *x> = ^ / " ^^^^ dt =
1 fx ,(t) — ,(x) dt
Г(—а)1 (х — г)!+а даёт дробную производную Маршо (см. [3, с. 368]).
Замечание 1. Согласно [3, с. 368] формула (8) справедлива и для т = 0 в (7), так как в этом случае и (г, —А) является равностепенно непрерывной полугруппой класса С0 (см. определение на с. 324 и пример на с. 325 в [3]).
зо
3. Обращение суммы генераторов
Напомним, что сильно непрерывная полугруппа линейных ограниченных операторов и (в) (в 2 0) по определению удовлетворяет условиям
1) и(0)р = р, р е Е;
2) и(* + в)р = и(*)и(в)р = и(в)и(*)р, р е Е;
3) Иш ||и(*)р - р||^ = 0, р е Е.
Потребуем также, чтобы для полугруппы выполнялась оценка ||и(з)р||^ ^ Ме-ад5, эд > 0, в 2 0. В таком случае если оператор А является производящим оператором полугруппы, т. е.
в
—и(в)р = Аи(в)р, р е Д(А), то справедливо представление обратного оператора
ргх
А-1р = - и (в,А)рвв.
0
Далее нам понадобится результат Ж. Да Прато и Д. Грисварда об обращении операторной суммы [11, ТеоремаЗ.З, с. 320].
Пусть X — банахово пространство с нормой || ■ ||х; А и В — два замкнутых в X линейных оператора с областями определения Да и О в, причём О в плотно в X. Рассмотрим оператор Ьр = Ар + Вр, р е Да П Дв = Дь.
Теорема 1. (Да Прато — Грисвард). Пусть
(А - А)-1(В - Л)-1р = (В - Л)-1(А - Л)-1р, Л > 0, р е Дь, и существуют такие Ма, Мв > 0, что для всех к е N
||(А - А)-к|| ^ Ма, И(В - А)-к|| ^ МВ, А > 0. (9)
Тогда Ь = А + В допускает замыкание Ь, при этом рь С (0, го),
||(Ь - А)-к|| ^ , А > 0, к е N (10)
(Ь - А)-1 (Да + Дв) С Дь, А > 0. (11)
Примером применения этой теоремы является доказательство корректной разрешимости задачи Коши для уравнения (см. [11, с. 332])
-^ + ВЦ*) - Аи(*) = /(*), * е (0,Т),
ш:
где / (*) — векторная функция со значениями в Е, принадлежащая пространству X = Ьр[(0,Т), Е] с нормой
гТ
1еш
0
Из теоремы получаем следствия.
1/р
, Р 2 1.
Следствие 1. Оператор Ь = А + В является генератором С0-полугруппы операторов и (г, А + В) в силу оценки (10) и соотношения (11).
Следствие 2. Если А и В — генераторы С0-полугрупп и (г, А) и и (г, В) с оценками II и (г, А) II ^ Ы1в-Ш11, II и (г, В )\\ ^ М2в-ад2*, шьш2 > 0, при этом для всех 5,г ^ 0
и (г, А)и (5, В = и (5, В )и (г, А)р, (12)
то оператор Ь = А + В имеет обратный Ь-1 и справедливо равенство
/><х
(А + В)-1 р = — и (5, А) и (5, В (13)
0
Доказательство. Для доказательства заметим, что в силу представления резольвенты [3, с. 45]
(А — А)-1р = Я(А,А)р = е-Л4и (г,А)р«г,
0
из (12) следует коммутируемость резольвент
Я(А,А)Я(А,В)р = е-Л4и (г,АП е-Л*и (з,В)рМ = 00
00
"ОО />оо
00
и соотношение
е-л(4+.)и (¿, А) и (5, В = е-л(4+.)и (5, в )и (г, = д(А, в)я(а, А)р (14)
¿•те
— (А + ВИ и (5, А)и (5, В = 0
/ Аи(5, А)и(в,В)^ — ви(5,В)и(в,А)^ = 00
-Г^и (5, В — -У-- и (з,А)р«8 =
/о ./о
I —[и (5, А) и (5, В = —и (5, А)и (5, В)р|~ = р.
./о
□
Этот результат позволяет обобщить равенство (13) на случай суммы нескольких операторов:
\ — 1
П \ 1 „оо п
р = — Ци (з,Аг)р^,
чг=1 / -70 ,= 1
где и(з,А,) — сильно непрерывные полугруппы с производящими операторами А,, такими, что для любых г,ш е (1,2, ...,п} и (з,,А,)и (зт,Ат)р = и(5т, Ат)и(5,, А,)р.
4. Доказательство корректности задачи без начальных условий
Для доказательства неравества (4) сначала рассмотрим случай дробных степеней операторов А и В, которые удовлетворяют теореме 1. Так как в силу оценок (9) эти операторы являются генераторами равностепенно непрерывных С0-полугрупп (см. [3, с. 324,325, пример 1]), то в соответствии с [3, с. 358] определены дробные степени (-А)а1 и (-В)а2, а е (0,1),г = 1, 2, при этом операторы Аа1 = -(-А)а1 и Ва2 = -(-В)а2 являются генераторами равностепенно непрерывных С0-полугрупп
и (в,А«1 )р = Н8,а1 (£ )и (£,А)рв£, (15)
'0
поо
и (з,Вод )р = к^ а )и (е,В)рве, (16)
0
где к8а1 (£),г = 1, 2 — функции Иосиды, являющиеся обратным преобразованием Лапласа функции е-8р"г, а е (0,1), г = 1, 2.
Применяя (15), (16), по аналогии с (14) можно получить коммутируемость полугрупп и (в, Аа1) и и (*,Ва2), а вместе с тем и представление обратного оператора
ргх
(А«1 + В«2)-1р = и(в,А«1 )и(5,В«2)рв5, р е Д(А«1) п Д^). 0
Далее заметим, что оператор, заданный выражением Вр = -^ и областью определения
Шр
Щ
является генератором сильно непрерывной полугруппы и (в, - )р(*) = р(* - в), с оценкой ||и(в, -)р||м ^ б-М8||р||м. Следовательно, согласно утверждению К. Иосиды [3, с. 358], оператор В имеет дробную степень
D(B) = <| р : р е dp е , (17)
Bар = — Р, а е (0,1), Р е
ва
При этом оператор - Ва является производящим оператором сильно непрерывной полугруппы
и (в, -В а)р(:) = ^~ ^,«(£)и (е, - р(*)ше = !о ~ К,а (£)р(: - £)в£,
где к8,а (£) — функция Иосиды, при этом справедлива оценка
- р
^ / эй*-е)iue ^ ^(еd£ми = 00
м-
18)
Записывая уравнение (2) с учётом (3) в операторной форме
n
^ Amu - Au = f, (19)
'm1
m=1
где Ат = , / = /(*), и учитывая коммутируемость операторов Ат и А на мно-
жестве (17), заключаем, что к уравнению (19) применима теорема 1 и следствие 2
о непрерывной обратимости суммы операторов, влекущее корректность задачи (2) (3) и представление её решения в виде
/>те>
u = —
П U(s,Am)U(s,A)/ds
/0 m=1
hs,ai (Ы — / hs,«m (Cm)U(s,A)/ (t - 6-----£m)d£m . . . d6ds
0 0 0
Отсюда, пользуясь (18), получаем оценку корректности
lull„ ^ M / e-s«iMai-s«2-------wsds
IM ^ / ^ IIM n
0 E am - w
m=1
5. Пример промежуточной асимптотики
Сравним решения задачи без начальных условий и задачи с начальными условиями для уравнения
да d2u(x t)
a—u(x,t) =—+ eMt cos ax, x G R, u ^ 0. (20)
dta ^ dx2
1. В случае задачи без начальных условий решение имеет вид
, . eMt cos ax
u(x,t) =-2. (21)
v y + a2 v 7
Это можно проверить непосредственно, воспользовавшись соотношением
тт( д2 \ 1 [~ -t,t 1 Г -,
Uls,—cos тх = о I e 4s cos T£d£ = 0 / e 4s cos т(x — a)da
1 лте
e-^ cos Tada cos tx = e-T s cos tx (22)
/ns J0
(см. [12, с. 207]). Отсюда следует оценка корректности
Hull,, ^
а^ьа + а2
2. В случае задачи с начальным условием и(х, 0) = 0 для уравнения (20) полугруппа
тт ( д \ . . Гр(г — 5), г ^ 5,
и А 5, —тт р(г) = г
V дгу ^; \0, г < 5,
даёт представление
и+ ( 5, — а|а) р(г) = ^те (ои+(е = £ к^аш — (23)
Пользуясь (22), (23) и формулами, приведёнными в [9, с. 445], получаем вид решения и+(х,г) при г ^ 0:
рте гЬ
u+(x,t) = I eMt I ha>as(C)e d£e ff2sds cos ax 00
те лте
те
pt i ra+i<x> r<x>
= в^Ч e-^-—. ep4 cos ax
JO Ja-i<x JO
e^t I - I -cos ax = — e-^Eaa--cos ax.
а Jo ' \ а J
(24)
л 2пг ара + а л
Здесь Еа,в (¿) — двупараметрическая функция Миттаг-Леффлера. Используя (21) и (24), перейдём к пределу
ШП Ц+М) = 11Ш Г' ^Еа (_^ше =
г^ж ц(х,г) а ^те у0 \ а у
= а^ Гж е-*Еаа (-ше = = 1. (25)
а Уо 'V а ^ а(^а + ^) ' '
Здесь использованы формулы из [9, с. 445]. Из (25) следует (6), так как
|u(x,t) - u+(x,t)||M —
+ а2
u+(x,t)
u(x, t)
Таким образом, u(x,t) даёт промежуточную асимптотику для u+(x,t).
Заметим, что аналогичный результат справедлив, если в правой части уравнения (20) стоит функция e^x sin ax. Отсюда следует, что если правая часть в (20) имеет вид e^xf (x), где / (x) — периодическая функция с рядом Фурье
/ (x) = ат + cos nx + sin nx),
n=1
то также имеет место промежуточная асимптотика.
Отметим, что численная реализация решения (21) проще, чем для (24). Авторы выражают глубокую благодарность В. Е. Федорову за полезные обсуждения затронутой в статье тематики.
Список литературы
1. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М. : Наука, 1966.
2. Баренблатт Г. И., Зельдович Я. Б. Промежуточные асимптотики математической физики // Успехи мат. наук. 1971. Т. 26, вып. 2 (158). С. 115-129.
3. ИосидаК. Функциональный анализ. М. : Мир, 1967.
4. Крейн С. Г. Линейные дифференциальные уравнения в банаховом пространстве. М. : Наука, 1967.
5. Костин А. В., Костин В. А. К теории функциональных пространств Степанова. Воронеж : Издат.-полиграф. центр ВГУ, 2007.
6. Зельдович Я. Б., Соколов Д. Д. Фрактали, подобие, промежуточная асимптотика // Успехи физ. наук. 1985. Т. 146, вып. 3. С. 493-505.
7. СамкоС.Г., КилбасА. А., МаричевО.И. Интегралы и производные дробного порядка и некоторые их приложения. Минск : Наука и техника, 1987.
8. Федоров В. Е., Туров М. М. Дефект задачи типа Коши для линейных уравнений с несколькими производными Римана — Лиувилля // Сиб. мат. журн. 2021. Т. 62, № 5. С. 1143-1162.
9. УчайкинВ.В. Методы дробных производных. Ульяновск : Артишок, 2002.
10. Костин В. А., Алкади Х. О разрешимости задачи без начальных условий для обобщённого уравнения с дробно-степенной суммой // Вестн. ВГУ. Сер. : Математика. Физика. 2022. № 3. С. 82-90.
11. DaPratoG., GrisvardP. Sommes d'operateurs lineaires et equations différentielles operationelles // Journal de Mathematiques Pures et Appliquees. 1975. Vol. 54. P. 305387.
12. ДвайтГ. Б. Таблицы интегралов и другие математические формулы. M. : Наука, 1983.
Поступила в 'редакцию 05.11.2022. После переработки 17.01.2023.
Сведения об авторах
Костин Владимир Алексеевич, профессор, доктор физико-математических наук, профессор кафедры математического моделирования, Воронежский государственный университет, Воронеж, Россия; e-mail: [email protected].
Костин Дмитрий Владимирович, доктор физико-математических наук, профессор кафедры математического моделирования, Воронежский государственный университет, Воронеж, Россия; ведущий научный сотрудник научно-исследовательской лаборатории математических методов теории управления и оптимизации, Воронежский государственный педагогический университет, Воронеж, Россия; e-mail: [email protected]. Алкади Хамса, аспирант, кафедра математического моделирования, Воронежский государственный университет, Воронеж; e-mail: [email protected]
Chelyabinsk Physical and Mathematical Journal. 2023. Vol. 8, iss. 1. P. 18-28.
DOI: 10.47475/2500-0101-2023-18102
PROBLEM WITHOUT INITIAL CONDITIONS FOR EQUATION WITH FRACTIONAL DERIVATIVES AND NTERMEDIATE ASYMPTOTICS
V.A. Kostin1'", D.V. Kostin1'2'6, H. Alkadi1
1 Voronezh State University, Voronezh, Russia
2 Voronezh State Pedagogical University, Voronezh, Russia "[email protected], [email protected]
The correct solvability of problems without initial conditions for fractional-power operator sums is shown. Solutions of problems without initial conditions are interpreted by Ya.B. Zeldovich and G.I. Barenblatt as intermediate asymptotics for problems with initial conditions. These authors point out the importance of such problems in connection with the expansion of the concept of "strict determinism"in statistical physics and quantum mechanics and raise the question of studying the properties of phenomena that do not depend on details in the initial conditions that manifest themselves after sufficient time. The present paper also provides an example of intermediate asymptotics for an equation with a fractional derivative.
Keywords: intermediate asymptotics, well-posed problem, Cauchy problem, equation without initial conditions, strongly continuous semigroup, fractional power of operator.
References
1. Tikhonov A.N., Samarskii A.A. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, Nauka, 1966. (In Russ.).
2. Barenblatt G.I., Zeldovich Ya.B. Intermediate asymptotics in mathematical physics. Russian Matematiccal Surveys, 1971, vol. 26, iss. 2, pp. 45-61.
3. YosidaK. Functional Analysis. Berlin, Gottingen, Heidelberg, Springer-Verlag, 1965.
4. KreinS.G. Linear Differential Equations in Banach Space. American Mathematical Society, 1971.
5. KostinA.V., KostinV.A. K teorii funktsional'nykh prostranstv Stepanova [On the theory of Stepanov's functional spaces]. Voronezh: Publishing and Printing Center of Voronezh State University, 2007. (In Russ.).
6. Zeldovich Ya.B., SokolovD.D. Fractals, similarity, intermediate asymptotics. Soviet Physics Uspekhi, 1985, vol. 28, iss. 7, pp. 608-616.
7. SamkoS.G., KilbasA.A., MarichevO.I. Fractional Integrals and Derivatives. Yverdon, Gordon and Breach, 1987.
8. Fedorov V.E., TurovM.M. The defect of a Cauchy type problem for linear equations with several Riemann — Liouville derivatives. Siberian Mathematical Journal, 2021, vol. 62, no. 5, pp. 925-942.
9. Uchaikin V.V. Metody drobnykh proizvodnykh [Methods of fractional derivatives]. Ulyanovsk, Artishok, 2008. (In Russ.).
10. KostinV.A., AlkadiH. O razreshimosti zadachi bez nachal'nykh usloviy dlya obobschyonnogo uravneniya s drobno-stepennoy summoy [On the solvability of the problem without initial conditions for a generalized equation with a fractional-power sum]. Vestnik Voronezhskogo gosudarstvennogo universiteta. Seriya: Matematika. Fizika [Bulletin of Voronezh State University. Series: Mathematics. Physics], 2022, no. 3, pp. 8290. (In Russ.).
The work was carried out with the financial support of the Ministry of Science and Higher Education of the Russian Federation within the framework of the state task in the field of science (topic number FZGF-2020-0009).
11. DaPratoG., GrisvardP. Sommes d'opérateurs linéaires et équations différentielles opérationelles. Journal de Mathématiques Pures et Appliquées, 1975, vol. 54, pp. 305387.
12. Dwight H.B. Tables of Integrals and Other Mathematical Data. New York, The Macmillan Company, 1947.
Article received 05.11.2022. Corrections received 17.01.2023.