Содержание к диссертации
Введение
Глава 1. Краевые задачи для дифференциальных уравнений дробного порядка и сопутствующим интегральные операторы
1.1. Постановка задачи 25
1.2. Начальные и краевые задачи для дифференциальных уравнений дробного порядка в локальной и нелокальной постановках 28
1.3.О качественных свойствах несамосопряженного интегрального оператора, порожденного уравнениями в частных производных смешанного типа 40
1.4. Об одном операторе, сопутствующем дифференциальному уравнению дробного порядка
1.5. Области в комплексной плоскости, где нет собственных значений оператора дробного интегрирования 58
1.6. Об одном способе приближенного решения уравнений дробного порядка 66
1.7. Применение операционного исчисления к решению краевых задач для дифференциального уравнения дробного порядка 70
Глава 2. Некоторые вопросы полноты систем собственных функций операторов дробного дифференцирования
2.1.О полноте систем собственных функций оператора дробного дифференцирования 75
2.2 . О полноте систем собственных функций оператора второго порядка с дробными производными в младших членах 87
Глава 3. Применение дифференциальных уравнений дробного порядка
3.1. Модель деформационно-прочностных характеристик одного класса полимеров на основе производных дробного порядка... 95
3.2. Анализ решения дифференциального уравнения дробного порядка, моделирующего особенности притока нефти к скважине в трещинном деформируемом пласте 101
Заключение 109
Список литературы 111
- Начальные и краевые задачи для дифференциальных уравнений дробного порядка в локальной и нелокальной постановках
- Области в комплексной плоскости, где нет собственных значений оператора дробного интегрирования
- . О полноте систем собственных функций оператора второго порядка с дробными производными в младших членах
- Анализ решения дифференциального уравнения дробного порядка, моделирующего особенности притока нефти к скважине в трещинном деформируемом пласте
Введение к работе
В связи с важными проблемами уравнений в частных производных смешанного типа и газовой динамики трансзвуковых течений за последние годы существенно возрос интерес к краевым задачам для дифференциальных уравнений дробного порядка. Тому подтверждение монография A.M. Нахушева [75], посвященная основополагающим элементам дробного исчисления и их приложениям к проблемам математического моделирования различных процессов в живых и неживых системах с фрактальной структурой и памятью. Интенсификации фундаментальных и прикладных исследований в области дробного дифференциального исчисления значительно способствовали проблеме описания нелокальных краевых задач со смещением для основных типов уравнений математической физики и биологии [78], выход в свет в 1972 г. книги "Избранные вопросы дифференциальных и интегральных уравнений" коллектива авторов Х.Г. Бжихатлова, И.М. Карасева, И.П. Лесковского, A.M. Нахушева [21], а также весьма содержательного обобщающего труда С.Г. Самко, А.А. Килбаса и О.А. Мари-чева "Интегралы и производные дробного порядка и некоторые их приложения", изданного в г. Минске в 1987 г. [91]. Необходимость развития теории дифференциальных уравнений сделала востребованными монографии Ю.И. Бабеико "Тепло- массообмеп. Метод расчета тепловых и диф-фузионых потоков" [16| и "The Fractional Calculus (Theory and Applications of Differentiation and Integration to Arbitrary Order) "[104] авторов Oldham Keitk В., Spanier Jerome, которая вышла в Лондоне в 1974 г.
Фундаментальное значение в развитии теории дифференциальных уравнений дробного порядка и ее прикладных аспектов безусловно сыграют монографии В.А. Нахушевой [84], А.В.Псху [88], Л.И.Сербиной (93).
Существуют различные определения операции дробного интегрирования и дифференцирования, которые для достаточного гладких функций, как правило, совпадают.
В предлагаемой работе эти операции понимаются в смысле Римана-Лиувилля и следуя А.М.Нахушеву [75] оператор дробного интстродиффе-ренцирования порядка а с началом в точке а Є К и с концом в точке х Є Ш порядка |а|, действующего на функцию ip(t) Є L[A, В] по переменной і, обозначается через D%x . По определению
D>(t) = {
г-а) [ -a) J
а < О,
sign (a; —a) f ip(t)dt
Г(-а) J {x-t)a+v
(p(t), а = О,
signal (Х - a)|^D-M-VM, а > О,
где [а] - целая часть а, которая удовлетворяет неравенству [а] < а < [а]+1; T(z) - гамма-функция Эйлера; а, х Є [А, В].
При 0 < а < 1 класс функций и = и(х), представимых в виде
и(х) = с{х - a)a_1 + D^f[t, u{t)], с = const, (1)
f[x,u(x)] Є L[a,b]
порождает обыкновенное дифференциальное уравнение порядка а:
DaaAt) = f(x,u)1 а<х<Ь. (2)
В частности, интегральное уравнение Абеля [98]
D^f(t) = и(х), 0 < а < 1, а < х < Ъ (3)
в случае, когда дробный интеграл D~xu{t) от правой части и(х) является абсолютно непрерывным на сегменте [а, Ь] порождает следующее простейшее дифференциальное уравнение дробного порядка
Da0xu(t) = f(x). (4)
Через АСа~1[а,Ь] обозначим класс функций и(х), представленный в виде (1). При а = 1 этот класс совпадает с классом АС[а, 6].
В силу (1) для уравнения (2) можно задать видоизмененное начальное условие
]іт{х-аУ~аи{х) = С. (5)
Поскольку D"~l(t — а)а~1 — Г (а) условие (5) можно заменить нелокальным условием
HmZ^W-Ci, (б)
х-а
где Сі, как и С, - заданное число.
Уравнение (4), как формула обращения уравнения Абеля (3), является первоосновой теории дифференциальных уравнений произвольного порядка, хотя предполагается (см. [91, с.615]), что началом развития этой теории следует считать дискуссию о способах решения уравнения
xD^u = и. (7)
Здесь Бг/2и ~ Dqx u(t) - дробная производная порядка 1/2 от функции и = и(х).
Так как
[xDte1/2U(t)] = xD^'2u(t) + D-alJ\(t),
то замена v = D^J u{t) Є АС[0,6] связывает (7) с вырождающимся при х = 0 дифференциальным уравнением неріюго порядка
(ху)' - D^v(t) - v{x), 0
Задаче Коши как о локальной, так и в нелокальной постановках для уравнений вида (2) с произвольным а > 0 посвящены многочисленные работы отечественных и зарубежных исследователей.
При 0 < а < 1 задачу (2), (5) принято называть задачей Коши в локальной постановке, а задачу (2), (6) - задачей Коши в нелокальной постановке.
Задаче Коши в локальной (в том числе классической) и нелокальной постановках для уравнения вида (2) посвящено значительное число работ. Среди работ последних лет следует отметить работы А.А. Килбаса и СМ. Марзана [61], А.А. Ворошилова и А.А. Килбаса [24], В.М. Головизнина и И.А. Короткина [39], В.А. Чадаева [95], [96).
В работах М.М. Држбашяна [40]-[49] исследованы задача типа Коши и задача типа Штурма-Лиувилля для обыкновенных дифференциальных уравнений дробного порядка.
Диссертационная работа выполнена в этом направлении. Здесь исследуются вопросы существования и единственности решений начальных и краевых задач для дифференциальных уравнений дробного порядка в ло-
кальной и нелокальной постановках, а также доказана теорема существования и единственности решения задачи Трикоми с нелокальным условием сопряжения и изучены качественные свойства нссамосопряженпого интегрального оператора, порожденного уравнением смешанного параболо-гиперболического типа. Исследованы вопросы полноты системы собственных функций краевых задач для дифференциальных операторов дробного порядка.
Подобные проблемы находятся в настоящее время в центре внимания многих исследователей, работающих в области дробного интегродиффе-ренцирования и его приложений.
Справедливо отмечено в монографии A.M. Нахушева "Дробное исчисление и его применение"(2003 г.) "... дробное (дифференциальное и интегральное) исчисление в теории фракталов и систем с памятью приобретает такое же важное значение, как и классический анализ в физике (механике) сплошных сред".
Дифференциальным уравнениям дробного порядка посвящены исследования Алеросва Т.С. [1]—[15], Барретта Д.Х. [100], Вольтерра В. [23], Ви-мана А. [106], Глушака А.В. [37]-[38], Джрбашяна М.М. [40]-[47], Зарубина А.Н., Карасева А.А., Килбаса А.А. [91], Летникова А.В., Маричева О.И. [91], Нахушева A.M. [75]-[83], Нахушсвой В.А. [81]-[85], Нерсесяна А.Б. [48]-[49], Псху А.В. [88], Репина О.А., Самко С.Г. [91], Ссрбиной Л.И. [93].
Вопросы распределения пулей функции типа Миттаг-Леффлсра исследованы многими авторами (см. Wiman А. [106], Седлсцкий A.M. [87],
[92], ПоповА.Ю. [86]-[87], Псху А.В. Щ).
Содержательность класса целых функций типа Миттаг-Леффлера и его асимптотические свойства объясняют интерес к вопросу о вещественных пулях, тесно связанному с разрешимостью спектральных задач для уравнений дробного порядка, рассматривался в работах Алсроева Т.С. [3]-[15], Джрбашяна М.М. [41]-[43], Нахушева A.M. [75]-[83]), Нерсесяиа А.Б. [48]-[49], Попова АЛО. [86]-[87], Псху А.В. [89], Ссдлецкого A.M. [87], [92].
Работа состоит из трех глав, посвященных исследованию краевых задач для дифференциальных уравнений дробного порядка и им сопутствующих интегральных операторов.
Сформулируем основные результаты диссертации.
В первой главе рассматриваются краевые задачи для дифференциальных уравнений дробного порядка и им сопутствующие интегральные операторы.
Пусть: 7 — {70)71)72} ~ точка из трехмерного евклидова простран-ства 1R3 с координатами 7j Є]0,1], j = 0,1,2; crj- = ХІ73'-1, / = 0* + 1, к — 0,1,2; р — 1/(72 > 0 j Dqx - оператор дробного дифференцирования порядка \у\ с началом в точке 0, который действует па функцию ip{x) из области его определения D{D$X) С L[0,1], где V(z) - гамма-функция Эйлера; L[0,1] - пространство абсолютно суммируемых функций;
ЩМі) = {
1 [ 4>(t)dt
г(-л)У (я-*)1-*' ^<и'
Введем интегродифференциальные операторы
Д } ~ dx-d-ть)7 - Г(1 - 7о) J (*- О70'
П }~ dx-V-n) <&*П } - Г(1 - 71) J {x- i)71 '
где p - 1 < a < p = 1,2,...;
d~t1-^ d71 d70 D " /W ~ ds-d-ft) dx^dx^f{x)-
Предполагается, что эти операторы имеют смысл, по крайней мере, почти
всюду на [0,1].
Легко видеть, что D^f(x) = D^x"lf{t), fc = 0,1;
(}-(<*) f(r)
-^-IU = DU(t), [а]-К a <„=[<»];
Рассматривается следующая задача: найти решение и(х) уравнения
D{a2)u(x) - [A + q{x)]u(x) = О, (8)
где А - спектральный параметр, q(x) Є С[0,1], из класса С}0,1] П L[0,1], удослегворяюгцее видоизмеиешгому (локальному) начальному условию Копій
V 1-70 f \ X V Ж1_7СИ(^)-^ х /0ч
hm Iі 7ои(ж) = do, lim — = ді, (9)
х-*0 х-*0 ж7'
где 7о > 1 ~ 72 j а ^о ^ <5i - заданные числа.
Решение задачи (8)-(9) ищется в классе С$[0,1] функций, прсдстави-мых в виде v(x)x~s, где v(x) С[0,1], 5 = const < 1.
Доказана
Теорема 1.1. Пусть q(x) удовлетворяет условию Липшица. Тогда в классе Cfi[0,1] задача (8)-(9) имеет и притом единственное решение и = и(х).
Рассматривается задача:
Найти решение и(х) уравнения
d&>u(t) = l\ + q(x)]u{x), (10)
где 8q'x = ^Щх~ - регуляризоваипый оператор дробного дифференцирования порядка 7i с началом в точке 0нс концом в точке х из класса С[0,1] П С*]0, 1], удовлетворяющее краевым условиям
и(0) = *0) и{1) = 62, (И)
где 8q и 5% - заданные числа.
Доказаны
Теорема 1.2. Пусть q(x) Є Cl[Q, 1] u q'{x) > 0, A > -g(0). Тогда задача (10)-(11} безусловно u однозначно разрешима.
Теорема 1.3. Пусть q(x) = const. Тогда условие
Ep(Xq;a2)^0, А > -#) (12)
является необходимым и достаточным условием однозначной разрешимости задачи (Ю)-(И). При выполнении условия (12) её решение и{х)
определяется формулой
и(х) = 50Ep(Xqx^ 1) + x^Ep(Xqx^ a2)[S2 - S0Ep(Xq; l)}/Ep(Xq; а2).
Теорема 1.4. Пусть q(x) = const, Ep(Xq; 1} ф 0. Тогда смешанная задача
lim ' = 6h и(1) = 62,
для уравнения
Dlul(t) = [\ + q(x)]u(x), (13)
имеет единственное решение
и(х) = б^Е^Х^;!) + 51T{a2)Ep(Xqx'72;a2)xai,
So = [fc - hr(~2)Ep(\q;a2)]/Ep(Xq; 1).
Теорема 1.5. Задача Коши и(0) = 6а эквивалентна интегральному уравнению Волыперра второго рода
и(х) = saxuEp{Xqx^;fi0) - СаХ^ЕріКр'^іи) +
-YY{a2)D-^Ep[Xq(x - t)^; a2]q0(t)u(t), (14)
где sa = <$оГ(^о) и ca ~ — d\T{pL\) - произвольные постоянные, Xq = A+g(0), q0(x) = q{x) - q{0),
CO J.
Е,М = Е,[г,Яа] = 1щШІ
fc=0
- функция типа Миттаг-Леффлера по терминологии М.М. Доісрбашяиа, которое в классе С[0,1] всегда разрешимо и притом единственным образом.
При q(x) = const решение и(х) задаётся формулой
и{х) = 6q[1 + \х*Ер(\х*\ 1 + 72)]. (15)
Уравнение (8) при 7о = 71 — 1 , 72 = <* — 1 переходит в уравнение
dgxu{t) = [\ + q(x)]u(x), (16)
где д$х = Dq~ ^2 - оператор, область определения которого принадлежит классу 52[0,1] всех функций и(х), измеримых на [0,1] вместе со своими производными до второго порядка.
Уравнение (8) при q(x) = 0, А = иа известно как дробное осциляци-онное уравнение.
Теорема 1.6. Задача Дирихле (11) для уравнения (16) в классе АСх[0,1] имеет единственное решение и{х), если
q(x)eC%l}1 А>-#), q'(x)>0. (17)
Во втором параграфе доказана теорема существования и единственности решения задачи Три ком и с нелокальным линейным условием сопряжения и исследованы качественные свойства несамосопряженпого интегрального оператора, порожденного уравнением смешанного параболо-гиперболического типа.
На евклидовой плоскости с декартовыми ортогональными координатами х и у рассмотрим модельное уравнение в частных производных смешанного (параболо-гиперболического) типа
1У] дх* ду1+п(~у)> [i*>
*W = (llVy>'
[О, Vy<0 - функция ХевисаЙда.
Уравнение (18) в верхней полуплоскости у > О совпадает с уравнением Фурье
д^ = ду' (19)
а в нижнем - с гиперболическим уравнением
Пусть: П - область, ограниченная отрезками AAq , AqBq , BqB прямых х — 0, у = уо, х ~ г соответственно и характеристиками
. s-, ^ і \ т+2 „ г.»-, > і -\ т+-
АС: х {-у)~ = О, ВС : х + (-у)"5" = г
m + 2^ J> ' m + 2v Jl
уравнения (20); Q.+ = {(х,у) : 0 < х < г, 0 < у < уо} - параболическая
часть смешанной области S1; П_ - часть области "2, лежащая в нижней
полуплоскости у < 0, ограниченная характеристиками АС, 0 < х < г/2 ,
ВС, г/2 < х < г и отрезком АВ; ]0, г[= {(х, 0) : 0 < х < г] ;
и+{х,у), V(x,y)ett+,
и =
(21)
и (ж,у), V(x,y) Є О ; Z?qx - оператор дробного интегродиффереицирования порядка |1| с началом в точке 0 [75, с. 9].
Для уравнения (18) в области Q рассмотрим краевую задачу типа задачи Трикоми с нелокальным условием линейного сопряжения и следующей постановке.
Рассм атр и в ается
Задача 1. Найти регулярное в областях Q+ , Q" решение (21) уравнения (18) со следующими свойствами:
и+ Є C(U+) Л Сг{П+ф, г[), и~ Є С(Q-) П С\а-и]0, г[); (22)
У+(Ж}0) = гГ(х,0) - XD^u~{t,0), 0<х<г;
(23)
<9(и+ — и )
у=0
= 0, 0 < ж < г;
(24)
и+(0,у) = p0fa), и+{г,у) = yv(y), 0 < 2/ < да;
(25)
и ,„ = ^(#), 0 < х < г, (26)
где U -замыкание Q, А - спектральный параметр, ipo(y) и
Нелокальное условие сопряжения задается уравнением (23), а локальные краевые условия (24) и (25) совпадают с условиями Трикоми. При А = 0 задача 1 идентична аналогу задачи Трикоми [20], [51].
Из (23), (26) при х = 0 и (25) при у = 0 следует, что равенство
V>0(0) = V(0) (27)
является необходимым условием согласования граничных данных.
Справедлива
Лемма 1.1. Пусть: существует решение (21) задачи (18); т(х) = = и~(х,0), v(x) = |f ; x~h{x) Є L[0,r]. Тогда
t"{x)-\t{x) = v{x), (28)
r(2P)D^r(t) = [pmv(x) + xtDl-Mt)]m, (29)
t(0) =
Отмстим, что в работах Алеросва Т.С. [3]—[10] подобные утверждения
доказывались исключительно для задач вида
u" + ADjj> = 0, 0<а< 1,
(1.4.20)
u(0) = 0, u(l) = 0.
(1.4.21)
Никакие обобщения для случая 1 < а < 2 полностью не проводились, хотя про возможность обобщения полученных результатов оговаривалось.
Таким образом, теорема 1,10 дает возможность полностью исследовать задачу (1.4.20)-(1.4.21) в случае 1 < а < 2.
И еще одно важное замечание. Оператор
А(,и = rv1) / ^ " *)^1у()^-
~ I x(l-t)i>-lu(t)dt
сходится к оператору
X 1
Aqu = {х- t)u(t)dt - х(1- t)u(t)dt
в обобщенном смысле, т.е. при р —> - \\Ар — Aq\\ -* 0, или, что то же
самое
X X
f^J— j(x - tfrlu{t)dt - ^^ J(x - t)S-lu(t)dt-
x l
- (х- t)u(t)dt + / x(l ~ t)u(t)dt
1 при p -+ - .
Если учесть, что Ad является оператором, обратным по отношению к
обыкновенному дифференциальному оператору второго порядка, то можно
утверждать, что при достаточно малых а можно получить спектр, близкий к обыкновенному дифференциальному оператору 2-го порядка.
1.5. Области в комплексной плоскости, где нет собственных значений оператора дробного интегрирования
Рассмотрим оператор
X і
Аи = j^y J(x - *)Ни(і)«й - jJ-_ Jx(l - tf>-lu(t)dt.
о 0
Положим, что 1 > - , отсюда следует, р > 2.
р 2
Отметим, что при р > 2 квадрат ядра оператора А не является интегрируемым, т.е. уже оператор А не действует из Ьч в hi.
Теорема 1.11. Вне круга с центром в начале координат и радиуса
—, тт оператор А не имеет собственных значений.
Г(1 + р-і)
Доказательство. В теории линейных операторов [62], [94] доказывается существование конечного предела
ро = lim VPi
п-їоо
Число ро называется радиусом ограниченного оператора В. Известно также, что
Ро < \\В\\ и что спектральный радиус ро характеризуется тем, что из неравенства
|А| > Ро
вытекает существование ограниченного оператора
(В -му1.
В частности, если В - вполне непрерывный линейный оператор, то / = |А0[, где Ао - наибольшее по модулю собственное значение оператора В . Оценим спектральный радиус оператора А сверху:
х і
IVі)
Ці + Р-1)'
о о
Из чего и следует доказательство теоремы 1.2.
Более точную оценку можно получить, если воспользоваться следующим известным утверждением:
Пусть N и М - линейные ограниченные операторы, тогда
PQ{N + M)
0{N) + po{M).
Так как наш оператор Л имеет вид
х 1
о о
то за операторы 7V и М примем операторы
^=г^/(*-0'~Ч0^
I Мч = -^~ j x(l-t)T\(t)dt.
Так как оператор N - вольтерров, то p(N) = 0.
Оператор М -одномерный, поэтому для тех р, для которых оператор М - вполне непрерывен, р{М) можно вычислить точно, р(М) = X, где А - собственное значение оператора. Так как К(х,Ь) - ядро оператора М - вырожденное, то характеристическое число оператора М определяется определителем Фредгольма d(X) этого уравнения
d{X) = \l-XKnl
(1+р-1)
гср-1)^!! = - у і(і -1)?-1^ = —ї ^1
о Очевидно, что норма ядра
K{x,t) = x(l-t)>~1
допускает оценку
I*-1
Итак, можно сформулировать следующую теорему.
Теорема 1.12. Вне круга с центром в начале координат и радиуса 1/Г(р-1) оператор А не имеет характеристических чисел.
Теперь мы проанализируем тот факт, что вне круга с центром в начале координат и радиуса 2/Г(р-1) оператор А не имеет собственных значений. Чем больше р, тем меньше 2/Г(р"1), а если учесть, как связаны собственные числа оператора А и нули ЕР(Х,2), то получается интересная картина, как для нулей функции ЕР(Х, 2), так и для собственных функций
и собственных значений оператора А. Значительный интерес представляет их расположение при р -> со , то есть спектр оператора Л при достаточно больших р.
Из приведенных теорем следует, что вне круга с центром в начале координат и радиуса 1/Г(р~1) функция .ЕДА; 2) не имеет нулей.
Если р достаточно мал, то в круге достаточно большого радиуса Ер(Х;2) не имеет нулей. Этот интересный факт можно использовать при изучении нулей ЕР(Х;2).
Вспомним асимптотику нулей функции Ep{z,\i) [41, с. 142]:
log/Л]
со.
7f = е±г'(2тпЦ> jl + o[-
Если р принимает достаточно большое значение, то ясно, что из выписанной асимптотики следует, что некоторые нули Ep(z,2) как угодно близко подходят к соответствующему единичному кругу.
При исследовании распределения нулей Ep{z,p) будет полезно провести анализ, основанный на следах степеней оператора Л, тем более, что даже в том случае, когда оператор Л но является ядерным, то есть и для тех р, для которых ряд
расходится (здесь Хп(р) -нули функции Ер(Х, 2)), следы для степеней Лп
легко вычисляются.
Рассмотрим оператор
X і
Ли = г(^) j(x ~ *)'"luWdi - f^T) J <l - ФМ*№-
о о
Ясно, что если р > 1, то оператор А пс является ядерным, но мы
всегда можем вычислить сумму соответствующих степеней характеристических чисел оператора А.
Теперь покажем, что оператор А при р диссипативиый, т.е. его мнимая компонента
неотрицательная. Здесь і - мнимая единица: і ~ \/~1-
Для этого рассмотрим задачу D :
Г(2-
Vf^+a"w = 0'
u(0) = 0, u(l) = 0.
Оператор
J (і-т)"'1 о
можно преобразовать следующим образом: обозначим через z{r) = и'(т), тогда имеем
/ и" (т) dT=rjirLdT.
J {t-r)0'1 J {t-r)a
о 0
Известно, что [75]
U J~T(l-a)dxJ (x-ty
Г(1 - а)
ДО) , Г f'(t)
+
(х - t)a
отсюда
f №
V{l-a)J [x~t)a о
dt =
і d f f(t) _ /(0)
T(l-a)dxj {x-ty" хаГ(1-а) о
Поэтому
f Ar)
T(l-a)J (t-r о
-dr =
1 d f z(t) , *(0)
-dr —
T(l-a)dtJ (t~r)a taT{l-a)'
Рассмотрим следующую линейную форму
Здесь и(х) - функция, комплексно-сопряженная функции и(х) = Reu(x)-\-+ilmu(x): й(х) = Reu(x) ~ гІггш(а;).
Знак (и, v) - это знак скалярного произведения, т.е.
Ж,Ч*)) = /«(,)№
Очевидно, что если z'(t) = и"{t), то
dx о
j {x-t)a xn [ ' v '
d f z(t)
dx о
fj3r^)-{f^\
d f u'{t)
dx J (x — t) о
-Лй(х)\ - (ШМх]у
Распишем
і / X
-л,адг/ /Si "Wdx-
d f u'{t)
о / о \o
Используя формулу интегрирования по частям, получим, что
1 / X
j^U)*= /Д^іад
о \0 / \0
1 / х \ 1
о \о / о
1 / X
-<0>/
о о \0
1 / X \ 1 / X
о \о / о \0 /
Последний оператор диссипативеп. Введем функцию
д(х), хЄ[0,є], ие(х) = <
и{х), х Є [є, 1], где є > 0, а д(х) такая, что '(ж)| 0 = 0. Очевидно, что
щ [1Лл)^х)г=
r^rdt] ,uf(x) \dx,
(х - 0 Г / "
поэтому оператор, порождающий задачу D , диссипативеп.
Из диссипативности последнего оператора следует, что его собственные значения лежат в правой полуплоскости комплексной плоскости. Из этого факта следует очень важное свойство функции типа Миттаг-Лсффлера.
Следствие 1.1. Все нули функции Ер{\\2) лежат в правой полуплоскости.
1.6. Об одном способе приближенного решения уравнений
дробного порядка
Рассмотрим задачу
и" + АГ(2 - a)Dlu = в(х), (1.6.1)
u(0) = 0, и(1) = 0. (1.6.2)
Здесь Q(x) - известная функция. Задача (1.6.1)-(1.6.2) эквивалентна интегральному уравнению Фредгольма второго рода
и(х) + А \ J х(1 - t)l-au{t)dt - Их - ty-au{t)dt і = F{x), (1.6.3)
і где F(x) — J G(x, t)Q(t)dt, причем
f(l-ar), 0
G{x,t) = <
s(l-t), t
Теорема 1.13. Пусть даны, два интегральных уравнения
ъ
S)ip(s)ds, (1.6.4)
а Ь
Ф{Ь) = Мі) + \[к0{і,з)ф{8)йз1 (1.6.5)
и пусть Ro(t,s}\) -резольвента уравнения (1,6.4) с ядром Ko(tys)
Пусть далее имеют место неравенства ь
\K{t,s)-K0(t,s)\ds
\№-Ш)\<т j\m^x)ds\
Тогда, если выполнено условие
ф|(1 + |Л|М)<1,
то уравнение (1,6.4) имеет единственное решение (p(t) и разность между этим решением и решением ф(і) уравнения (1.6.5) удовлетворяет
оценке
і /ч ~<^ JV|A|(1 + |Л|М)2е 1л т,_ч
|А|(1 + |Л|М)е
где N= sup i/(t)j
Чтобы воспользоваться этой теоремой, в уравнении (1.6.3) ядро
1-а/
-1(1-^-4(1-1-
1, t <х,
G0(x,t) = <
О, t>x заменим на вырожденное ядро ~х(1 — t)[~Q.
Тогда наша задача сведется к исследованию уравнения
(х) + A fx(l - tf-au{t)dt = F{x).
(1.6.6)
Ядро последнего уравнения - вырожденное, его решение задается по
известной формуле
и(х) = F(x) + Л / Я(х, t, X)F(t)dt,
где R(x,t,X) - резольвента Фредгольма:
D(x,t,\)
R(x,t,X) =
d(A) '
00 \т
Здесь D{x,t,X) = Y(~l)m—Bm(x;t), где
m=0
/<"(,f) K(x,a\) ... K(x,ar
Bm{x,t) = j... j
о 0
К(ац,і) K{a2,al) ... K(am+i,am)
dai... da„
причем Kq (x,t) ~ Вычислим є:
K{am,t) К(ат,ац) ... K(am+\,Qm)
x(l - t)l~a, где d(X) = |1 - XKn\, где
f І [-х{1 -1)1-" + {x- t)l~aG{x, t) + x(l - t)l~a] | dt =
Поэтому К*)- Ф)1< 2з(2з-47А)-Итак, доказана теорема. Теорема 1.14. Если 0 < А < (3 — а) (2 — а)/24, то уравнение
и(х) + \1 ( х(\ - t)l-nu(t)di - j(x-t)l-au(t)dt\ =F{x) имеет единственное решение и разность между этим решением и реше-
нием
u=+%y$^fr-*-w
уравнения
(х)
Г(2 - а)
о удовлетворяет оценке
23(23 -47А)
1.7. Применение операционного исчисления к решению краевых задач для дифференциального уравнения дробного
порядка
Для уравнения
и"+ >> = Ли (1.7.1)
рассмотрим задачу
u(0)-/?u'(0) = 0, ы(1) = 0, (1.7.2)
где а ]0,1[, /3 Ж, X - спектральный параметр.
Проинтегрируем уравнение (1.7.1) дважды от 0 до х. Тогда, пользуясь формулами перестановки Дирихле и свойствами гамма-функции Эйлера, получим
XX X
u{t)dt. (1.7.3)
ооо Запишем (1.7.3) в виде
х х I t \ X
Idu'{t) + nb)fd(Iw^rd4 =xIu^dt (L7-4)
о о \о /о
Перепишем (1.7.4) в виде
^^щЬ)(/(Г%^*)1о=А/оМЛ {L7-5)
\о /о
или, что то же самое,
X X
"'W - "'(0) + Г(гЬ) / ^ҐГМ = А / U{t)dt (17-6)
Проинтегрируем (1.7.G) еще раз от нуля до х , получим
X ( t
Ju\t)dt~uWjdt + ~--jj
Г «(О
(t - 0а
d}dt =
о ко
х t
= Л
I ju(t)d№
о о или, что то же самое
(1.7.7)
х ( t
и(х) - и(0) - и'(0)х +
Г(1-а)У W (*-fl о U
dO <й =
= А / fu()ddt.
о о Так как
Г(1-а)/ / (t-ty о 1о
^ )dt =
(1.7.8)
У (* - о,
то (1.7,8) можно переписать в виде
и(х) - и(0) - u'(0)ic + p/g1, ч f(x - ty-au{t)dt -
х t
A f fu()dtdt.
(1.7.9)
о о
Теперь очевидно, что решение уравнения (1.7.9), удовлетворяющее условию
u(0)-/7u'(0) = 0, совпадает с решением уравнения
X t
-A f fu{t)dtdt = u'{Q)(x + P). (1.7.10)
о о
Уравнение (1.7.10) решим с помощью преобразования Лапласа.
Пусть U(s) - изображение функции и(х), т.е. U(s) = и(х) или, что то же самое [50]
U(s)= I' e~3tu(t)dt.
о Естественно, здесь предполагается, что решение уравнения (1.7.10) ищется
в классе функций u(t), где обеспечивается сходимость интегралов Лапласа.
Ясно, что функция
J(x - ty-au(№
является сверткой функций хг~а и функции и{х).
Вспомним, что свертка функций /i(i), /2(і) определяется следующим
образом;
F(t) = j h(r)W~T)dr.
И теперь, чтобы найти изображение функции
j(x-i)l'au(t)dt
остается найти изображение функции xl~Q.
Известно, что изображение функции х'1 при р, > — 1 равняется Г(/і + l)s~1-''. Теперь из теоремы о свертывании оригиналов, т.е. из теоремы, которая гласит, что
t Jh(T)f2(t-T)dr = F1(s)F2(s), о
если fi(r) = Fi(s), /2(т) === F2(s), следует
f{x - t)l-nu(t)dt = U{s)sa-2Y{2 - a). о
Далее найдем изображение функции
X t
jdtju(i)di. о о
Очевидно, что
х і
Jdtju(№? О о
X і
.и (')
л '
Теперь изображающее уравнение имеет вид
1 А
1 , ^
1 +
U(s)
s2-a s2
«ч)(і+
отсюда имеем
w s2 + sa-\
Теперь, пользуясь обратным преобразованием Лапласа, т.е. формулой
(7+ІСО
т = h\ eS'F{s)ds'
(7-гоо
двойственной к формуле
F(s) = J e-stf{t)dt,
получим, что
С+ІО0
и'(0) и(х) =
ds.
esx{l + (3s) s2 + sa -А
2ттг J
с-гоо
Из последнего соотношения следует, что А является собственным значением задачи (1.7.1)-(1.7.2) тогда и только тогда, когда А является нулем функции м(1, А):
С+200
U(^X) = ~KZT I .2 , ..„ \ds-
ц'(0) Г es{l + Ps) 2т J s2 + sa-X'
с—ICO
Таким образом, доказаны теоремы:
Теорема 1.15. Число А является собственным значением задачи (1.7.1)-(1.7.2) тогда и только тогда, когда А является нулем функции
с+гсо
т = /
es(l + (3s)
s1 + sa - А
С-100
ds.
Теорема 1.16. Пусть Xj является нулем функции U(X). Тогда соответствующей собственному значению Xj собственной (функцией задачи (1.7.1)-(1.7.2) является функция
с+гсо
ттп ^ f ^ +fob
ЩХі>х)= J ^T^~xds-
с-гоо
Начальные и краевые задачи для дифференциальных уравнений дробного порядка в локальной и нелокальной постановках
В связи с важными проблемами уравнений в частных производных смешанного типа и газовой динамики трансзвуковых течений за последние годы существенно возрос интерес к краевым задачам для дифференциальных уравнений дробного порядка. Тому подтверждение монография A.M. Нахушева [75], посвященная основополагающим элементам дробного исчисления и их приложениям к проблемам математического моделирования различных процессов в живых и неживых системах с фрактальной структурой и памятью. Интенсификации фундаментальных и прикладных исследований в области дробного дифференциального исчисления значительно способствовали проблеме описания нелокальных краевых задач со смещением для основных типов уравнений математической физики и биологии [78], выход в свет в 1972 г. книги "Избранные вопросы дифференциальных и интегральных уравнений" коллектива авторов Х.Г. Бжихатлова, И.М. Карасева, И.П. Лесковского, A.M. Нахушева [21], а также весьма содержательного обобщающего труда С.Г. Самко, А.А. Килбаса и О.А. Мари-чева "Интегралы и производные дробного порядка и некоторые их приложения", изданного в г. Минске в 1987 г. [91]. Необходимость развития теории дифференциальных уравнений сделала востребованными монографии Ю.И. Бабеико "Тепло- массообмеп. Метод расчета тепловых и диф-фузионых потоков" [16 и "The Fractional Calculus (Theory and Applications of Differentiation and Integration to Arbitrary Order) "[104] авторов Oldham Keitk В., Spanier Jerome, которая вышла в Лондоне в 1974 г. Фундаментальное значение в развитии теории дифференциальных уравнений дробного порядка и ее прикладных аспектов безусловно сыграют монографии В.А. Нахушевой [84], А.В.Псху [88], Л.И.Сербиной (93).
Существуют различные определения операции дробного интегрирования и дифференцирования, которые для достаточного гладких функций, как правило, совпадают. В предлагаемой работе эти операции понимаются в смысле Римана-Лиувилля и следуя А.М.Нахушеву [75] оператор дробного интстродиффе-ренцирования порядка а с началом в точке а Є К и с концом в точке х Є Ш порядка а, действующего на функцию ip(t) Є L[A, В] по переменной і, обозначается через D%x .
Задаче Коши в локальной (в том числе классической) и нелокальной постановках для уравнения вида (2) посвящено значительное число работ. Среди работ последних лет следует отметить работы А.А. Килбаса и СМ. Марзана [61], А.А. Ворошилова и А.А. Килбаса [24], В.М. Головизнина и И.А. Короткина [39], В.А. Чадаева [95], [96).
В работах М.М. Држбашяна [40]-[49] исследованы задача типа Коши и задача типа Штурма-Лиувилля для обыкновенных дифференциальных уравнений дробного порядка.
Диссертационная работа выполнена в этом направлении. Здесь исследуются вопросы существования и единственности решений начальных и краевых задач для дифференциальных уравнений дробного порядка в ло кальной и нелокальной постановках, а также доказана теорема существования и единственности решения задачи Трикоми с нелокальным условием сопряжения и изучены качественные свойства нссамосопряженпого интегрального оператора, порожденного уравнением смешанного параболо-гиперболического типа. Исследованы вопросы полноты системы собственных функций краевых задач для дифференциальных операторов дробного порядка.
Подобные проблемы находятся в настоящее время в центре внимания многих исследователей, работающих в области дробного интегродиффе-ренцирования и его приложений.
Справедливо отмечено в монографии A.M. Нахушева "Дробное исчисление и его применение"(2003 г.) "... дробное (дифференциальное и интегральное) исчисление в теории фракталов и систем с памятью приобретает такое же важное значение, как и классический анализ в физике (механике) сплошных сред".
Области в комплексной плоскости, где нет собственных значений оператора дробного интегрирования
Дифференциальным уравнениям дробного порядка посвящены исследования Алеросва Т.С. [1]—[15], Барретта Д.Х. [100], Вольтерра В. [23], Ви-мана А. [106], Глушака А.В. [37]-[38], Джрбашяна М.М. [40]-[47], Зарубина А.Н., Карасева А.А., Килбаса А.А. [91], Летникова А.В., Маричева О.И. [91], Нахушева A.M. [75]-[83], Нахушсвой В.А. [81]-[85], Нерсесяна А.Б. [48]-[49], Псху А.В. [88], Репина О.А., Самко С.Г. [91], Ссрбиной Л.И. [93].
Вопросы распределения пулей функции типа Миттаг-Леффлсра исследованы многими авторами (см. Wiman А. [106], Седлсцкий A.M. [87], [92], ПоповА.Ю. [86]-[87], Псху А.В. Щ). Содержательность класса целых функций типа Миттаг-Леффлера и его асимптотические свойства объясняют интерес к вопросу о вещественных пулях, тесно связанному с разрешимостью спектральных задач для уравнений дробного порядка, рассматривался в работах Алсроева Т.С. [3]-[15], Джрбашяна М.М. [41]-[43], Нахушева A.M. [75]-[83]), Нерсесяиа А.Б. [48]-[49], Попова АЛО. [86]-[87], Псху А.В. [89], Ссдлецкого A.M. [87], [92]. Работа состоит из трех глав, посвященных исследованию краевых задач для дифференциальных уравнений дробного порядка и им сопутствующих интегральных операторов. Сформулируем основные результаты диссертации. В первой главе рассматриваются краевые задачи для дифференциальных уравнений дробного порядка и им сопутствующие интегральные операторы. Каждое нетривиальное решение этой задачи назовем собственной функцией, а соответствующее значение параметра А - собственным значением. В [44] излагается метод построения биортогоиальных систем собственных и присоединенных функций, порожденных задачей (2.1.3)-(2.1.5). М.М. Джрбашян [44] отмечает, что вопрос полноты систем собственных функций задачи (2.1.3)-(2.1.5) или более тонкий вопрос о том, что составляют ли эти системы базис в 1 (0,1), имеет безусловный интерес, но их решение, по-видимому, сопряжено со значительными аналитическими трудностями.
При q(x) = 0 полнота системы собственных функций была установлена Алероевым Т.С. [3]. Доказательство было получено с помощью теоремы Л и декого [64] [66]. Позже М.П. Маламуд [67]-[68] и его ученики также привели результаты, посвященные вопросам полноты системы собственных и присоединенных функций подобных задач. Их исследования опираются на известный аналог теоремы Наймарка [72]—[74]. Как отмечает М.П. Маламуд, их исследования дают возможность установить полноту только для аналитических потенциалов. Когда q(x) полуограничена, их методика не проходит.
Теперь изучим вопросы, связанные с полнотой системы собственных и присоединенных функций оператора А. Система {tpk} {к = 0,1,...) элементов линейного топологического пространства М называют полной, если замыкание ее линейной оболочки совпадает с М. Или, что то же самое: Систему {ipk} {к = 0,1,...) функций, голоморфных в области G, называют полной, если какова бы ни была голоморфная в G функция f(t), существует последователь!юсть линейных комбинаций рк CnXtpnit) с постоянными коэффициентами, которая равномерно сводится в области G к функции /.
Далее из того, что А-1 диссипативен и ядерный, следует, что система собственных и присоединенных функций оператора А"1 полна в г(0, 1) -Теорема Лидского [G4]-[66J утверждает, что условием полноты системы собственных и присоединенных функций линейного ограниченного оператора является ядерность и диссипативпость.
. О полноте систем собственных функций оператора второго порядка с дробными производными в младших членах
Теперь из теоремы Лидского [64]-[66] следует полнота системы собственных и присоединенных функций оператора А-1, или, что то же самое, полнота системы собственных и присоединенных функций оператора А. Теперь можно выписать систему собственных и присоединенных функций оператора А, т.е. конкретно выписать полную в 1 (0,1) систему функций. Еще раз отметим, что для полуограниченного потенциала q(x) подобный результат получен только в этой работе. 2.2. О полноте систем собственных функций оператора второго порядка с дробными производными в младших членах
В работе [58] были изучены деформационно-прочностные характеристики одного класса ароматических хлоросодержащих полиэфиров, которые обладают высокой прочностью аист 60% и синтезированных в работах Кехарсаевой Э.Р., Микитаева А.К., Алероева Т.С. [55]-[57],
На участке "больших" деформаций строится модель указанных полимеров, основанная на производных дробного порядка. В общем виде стандартная модель линейного вязкоупругого тела записывается в виде [19], [90], где а — cr(t) и є є(і) - напряжение и деформация в момент времени t; b(, EQ, ЕІ - заданные величины.
Линейное дифференциальное уравнение (относительно а при заданном є или наоборот) определяет соотношение между деформацией и напряжением. Легко видеть, что, известная модель Максвелла является частным случаем модели (3.1.1), и, следуя [75], будем называть моделью Бегл и Р. Л. и Торвика П.Дж. [19], [99].
В линейном уравнении (3.1.2) т означает время релаксации, Е - модуль упругости (например, полимера [18]). Впервые дробное исчисление при изучении вязкоупругих свойств полимеров использовано А. Джеман-том в работах [102], [103] и предложил использовать при моделировании многих вязкоупругих материалов производные дробного порядка.
Из модели (3.1.3) не следует уравнение Максвелла, что означает, что модель (3,1.3) лучше описывает деформационно- прочностные характеристики, чем модель Максвелла. Сами авторы модели (3.1.1) утверждают [75], что для моделирования многих вязкоупругих тел достаточно ограничиться моделью а + ЫУ&а = Е0є + ife (3.1.4) которая содержит всего пять параметров: Ъ, ft, EQ, Е\ И а, причем 0 а 1, 0 /? 1. Очевидно, что модели (3.1.1) и (3.1.4) являются частными случаями (3.1.3). Если в модели (3.1.4) положить j3 = Е Е$г (г = EQ Е"1), а затем предположить, что .Е"1 достаточно мало (полимер не имеет мгновенной упругости, а паши полимеры не имеют ее), то сс можно аппроксимировать уравнением, которое является частным случаем модели (3.1.4), а именно a(t) = E0D%te(t). (3.1.5) Уравнение (3.1.5), предложенное Скотт-Блэром [19, с.25], [22, с. 75], называют реологическим уравнением состояния. В [84] показано, что если напряжение a = a(t) удовлетворяет условию 1\та(є) -Iа = 0, то уравнение (3.1,5) можно заменить определяющим уравнением a{t) = Eodfeit), (3.1.6) то 6 - оператор дробного дифференцирования порядка а в смысле M.Caputo[101]: dfeit) = Dfc(t) - т{1Ф)а)іа, 0 а 1. Приведенный анализ в монографии [84] говорит о том, что если па-пряжение и = а (і) и деформация є — є (с) за конечное время не могут обращаться в бесконечность, то механические свойства вязкоупругих материалов, в частности, полимеров, можно описать не уравнением (3.1.4), а реологическим уравнением состояния (3.1.6).
Таким образом, на участке больших деформаций [58], т.е. на участке 0,06 є 0,18 напряжение a(t) из-за достаточной малости Е-1 = 1-Ю3 для наших полимеров, то можно аппроксимировать интегральной кривой следующего уравнения: a(t) = д 65є. На участке вязкого течения а = 88 [58] существенным является тот факт, что значения формы {DQte {i),e{i)) лежат в угле JargAj , если є(0) = = e(t) = 0. Известно [16], что для уравнения (3.1.5) имеет место теорема существования и единственности, т.е. зная одну точку, через которую проходит интегральная кривая, можно восстановить всю интегральную кривую этого уравнения. В работе [75] определяющее уравнение eE = D a (3.1.7) названо обобщенным законом высокоэластичной деформации.
Анализ решения дифференциального уравнения дробного порядка, моделирующего особенности притока нефти к скважине в трещинном деформируемом пласте
Общее решение уравнения (3.1.7) (относительно є = e(t)) при любом а Є [0,1] задается формулой еЕ = СЄ 1 + 07V, (3.1.8) где с - произвольная постоянная величина. Если к дифференциальному уравнению дробного порядка а (3.1.6) присоединить начальное условие limrt- o, (3.1.9) ТО С = QE . Если Го = 0, то модель (3,1.5) с условием (3,1.9) совпадает с моделью (3.1.7). Видим, что определяющее уравнение а + тЩьо = rEDfe, т = ± р = а+1 (3.1.10) является следствием модели (3.1.4). Модель (3.1,10) при а = 0 совпадает с моделью Максвелла (3.1.2) [75, с. 175]. Поэтому на (3,1,10) можно смотреть как на обобщенное уравнение или па его аналог в дробном исчислении.
Теорема. (Мацаева-Паланта [70]). Значения линейной формы (Jau}u) лежат в угле argA —, т.е. значения (Jl u,u) лежат в угле argA — . Это утверждение, безусловно, имеет физическое истолкование, хотя это и выходит за рамки пашей работы. Определенные особенности притока жидкости в проницаемых средах изучаются по данным гидродинамических исследований скважин [17]. При разработке нефтяных залежей используются различные зависимости дебита от перепада давления [97]. Безусловно, в общем случае, эти зависимости не являются линейными. Отклонение рассматриваемых зависимостей от линейной объясняется причинами, основными из которых являются: деформация коллектора, инерционные силы сопротивления, изменение свойств пласта и жидкости, а также действующей толщины пласта.
В призабойной зоне нарушается линейный закон [17]. Существует функциональная зависимость, учитывающая ииерциональ-ные составляющие сопротивления движению жидкости
Здесь Vp - градиент давления, // - динамическая жидкость, v - скорость фильтрации, к - проницаемость среды, )3 - скалярная величина, зависящая от модуля вектора скоростей, /(/?, v) - безразмерная функция, полученная согласно тг-теореме анализа размерности [17]. Если предположить, что функция f(0,v) разлагается в ряд Тейлора и ограничиться двумя первыми членами разложения, получим Vp=_ y- OT. (3.2.1) к к Справедливость полученной двучленной зависимости допускается при описании процесса фильтрации жидкости в трещинном коллекторе [17], [97], с ростом градиента давления изменяется фильтрационная способность коллектора. Это может быть связано с изменением действующей толщины пласта, которая для каждой конкретного коллектора при различных градиентах давления JVp (до критического значения Уркр) различна. При достижении VpKp действующая толщина пласта достигает максимального значения и остается неизменной при дальнейшем увеличении градиента давления (см. рис. 1) [97].
Зависимость действующей толщины пласта от градиента давления Последующее измерение Vp] может привести к изменению некоторых фильтрационных каналов и к уменьшению действующей толщины пласта. Подобные явления часто отмечаются на практике, в частности, на нефтяных залежах Белоруссии и Ставрополья [4], [97].
В диссертации, посвященной исследованию линейных краевых задач для дифференциальных уравнений дробного порядка и их применению, получены следующие основные результаты: 1. Исследованы начальные и краевые задачи для широкого класса линейных дифференциальных уравнений дробного порядка в локальной и нелокальной постановках. 2. Изучены качественные свойства нссамосопряжеиного интегрального оператора, порожденного уравнением в частных производных смешанного (параболо-гиперболического) типа. 3. Для класса интегральных операторов, сопутствующих дифференциальным уравнениям дробного порядка, получены в терминах функций типа Миттаг-Леффлера явные выражения для собственных функций и собственных значений. 4. Доказана теорема, которая позволяет выделить области в комплексной плоскости, где нет собственных значений для класса линейных дифференциальных операторов дробного порядка. 5. Доказана теорема о системе собственных и присоединенных функций класса линейных дифференциальных операторов дробного порядка с полуограниченным потенциалом. 6. Для дифференциального оператора второго порядка с дробными производными в младших членах доказана теорема о полноте системы собственных и присоединенных функций с полуограниченным потенциалом. 7. Для дифференциального оператора второго порядка с дробными производными в младших членах установлена его диссипативностъ в случае полуограниченного потенциала. 8. Проведен конструктивный анализ решения дифференциального уравнения дробного порядка, моделирующего особенности притока нефти к скважинам в трещинном деформируемом пласте.