Разностные методы решения краевых задач для уравнений параболического типа с дробной производной в граничных условиях тема автореферата и диссертации по математике, 01.01.02 ВАК РФ

Худалов, Марат Захарович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Владикавказ МЕСТО ЗАЩИТЫ
2003 ГОД ЗАЩИТЫ
   
01.01.02 КОД ВАК РФ
Диссертация по математике на тему «Разностные методы решения краевых задач для уравнений параболического типа с дробной производной в граничных условиях»
 
Автореферат диссертации на тему "Разностные методы решения краевых задач для уравнений параболического типа с дробной производной в граничных условиях"

Худалов Марат Захарович

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

01.01.02 — "Дифференциальные уравнения"

Автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

Нальчик — 2003

Работа выполнена в Северо-Осетинском государственном университете им. К. Л. Хетагурова.

Научный руководитель доктор физико-математических наук,

профессор Ш хану ков Мухамед Хабалович

Официальные оппоненты: доктор физико-математических наук,

профессор Алероев Темирхан Султанович; кандидат физико-математических наук, доцент Нахушева Зарема Адамовна

Ведущая организация: ' Таганрогский радиотехнический университет

Защита состоится "13 " сентября 2003 г. в часов на заседании диссертационного совета КР 002.195.45 в Научно-исследовательском институте прикладной математики и автоматизации Кабардино - Балкарского научного центра Российской академии наук по адресу: 360000, г. Нальчик, ул. Шортанова, 89А.

С диссертацией можно ознакомиться в научной библиотеке НИИ ПМА КБНЦ РАН.

Автореферат разослан августа 2003 г.

Ученый секретарь диссертационного совета, к.ф.-м.н. / ПсхуА.В.

л?

2ooJ-fl 177S5" -з-

Общая характеристика работы

Актуальность темы. Использование дробных производных для описания и изучения физических процессов стохастического переноса стало в последние годы одной из популярных областей физики, многие проблемы фильтрации жидкости в сильно-пористых (фрактальных) средах приводят также к необходимости изучения краевых задач для уравнений в частных производных дробного порядка.

Уравнениям в частных производных дробного порядка посвящены работы Нахушева A.M., Нигматулина P.P., Чукбара К.В., Шханукова М.Х., Гекиевой С.Х., Керефова A.A., Алероева Т.С., Кумыковой С.К., Псху A.B., Климентовой В.Б. и других. "

Настоящая диссертационная работа посвящена построению и исследованию разностных схем для уравнений параболического типа с дробной производной в граничных условиях при любом а 6 (0,1), где а -порядок дробной производной.

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

Методы исследования. Результаты работы получены с использованием метода энергетических неравенств, теорем вложения в дифференциальной и разностной трактовках, метода Ротэ и принципа максимума.

Научная новизна. В диссертационной работе получены следующие результаты:

1. Для решения уравнения параболического типа общего вида с

дробной производной по времени в граничнь .at углошмк подучена ,апрп

РОС. НАЦИОНАЛ

БИБЛИОТЕКА |

С-Петербург !

ОЭ ««*VjiJ ;

орная оценка, откуда следует единственность и непрерывная зависимость решения от правой части и начального условия.

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

3. Для решения нагруженного уравнения параболического типа с нелокальным условием получены априорные оценки в дифференциальных и разностных трактовках, откуда следует устойчивость и сходимость разностной схемы.

4. Для решения нелокальной по времени задачи для обобщенного уравнения переноса дробного порядка получена априорная оценка. Откуда следует единственность и непрерывная зависимость решения от входных данных.

5. На основе принципа максимума получены оценки в равномерной метрике для решения разностной задачи, аппроксимирующей нелокальную задачу для обобщенного уравнения переноса дробного порядка, откуда следует сходимость разностной схемы в равномерной метрике.

Практическая и теоретическая ценность. Работа является теоретической. Результаты работы могут быть применены в теории дифференциальных уравнений дробного порядка, а так же при изучении физических процессов стохастического переноса и проблем фильтрации жидкости в сильно-пористых (фрактальных) средах.

Апробация работы. Результаты работы докладывались на семинаре по современному анализу, информатике и физике НИИ ПМА КБНЦ РАН, на Второй Международной конференции "Нелокальные краевые задачи и родственные проблемы математической биологии, информа-

тики и физики" (Нальчик, 3-7 декабря 2001 г.), на семинарах математических кафедр СОГУ, на Международном Российско-Узбекском симпозиуме " Уравнения смешанного типа и родственные проблемы анализа и информатики".

Публикации. Основные результаты работы опубликованы в работах [1]-[6]. Из них [1],[5],[6] в соавторстве. Работа [1] выполнена в соавторстве с Юртиным И.И., которому принадлежит постановка задачи, [2] выполнена в соавторстве с М.Х. Шхануковым и Э.Е. Рейдером, Э.Е. Рейдеру принадлежит постановка задачи, М.Х. Шханукову общие указания о путях решения, в работе [3] выполнена в соавторстве с М.Х. Шхануковым, М.Х. Шханукову принадлежит постановка задачи.

Структура и объем работы. Работа состоит из введения, трех глав и списка литературы. В первой главе - шесть параграфов, во второй - одинадцать параграфов, в третьей - двенадцать, список литературы содержит 74 наименований. Объем диссертации - 85 страниц, набранных с использованием пакета Т^Х.

Содержание работы

Во введении дан обзор публикаций по теме диссертации и приведены основные результаты диссертации.

Глава 1. В области С?г = {(я, : 0 < х < /, 0 < £ < Т} рассмотрена задача

0 < х < 1,0 < I < Т, (1)

= /?2(*)и(/,<) +/?з(*)£>0" «(I, О -1Л2(1)

(2)

и(х,0) = и0(х),

(3)

где *(*,<) ^ С1 > О, ШШШШ^ с2, 0 < о < 1,. оператор дробного интегро - дифференцирования (в смысле Римана -Лиувилля).

Коэффициенты уравнения (1) и граничных условий (2) удовлетворяют следующим условиям гладкости: к(х,1) 6 С2,1 (От), С2,1 (От), где Ст'п[Ят)- класс функций, имеющих тп производных по х, и п производных по I, непрерывных в От-

Для решения задачи (1)-(3) получена априорная оценка

«

о

где ||и||> = /и>(х,№, ||«|||1<г1 = /Ци(х,т)№т. о о

Из оценки (4) следует единственность решения и непрерывная зависимость решения от входных данных.

Для решения задачи (1)-(3) построена разностная схема

Уг=к№+ф, (5)

у(в)=«гу+(1-<г)1/, у =

Ау = <

л » =-о^л-■ г = 0'

Лу = (ауг)х -¿у, 1 ^ г ^ N - 1, А+ Йл/Уг.л, + 0.5ЛДгулг Дз ла . м

0.5/г

0.5/г

где

1 17 = г(2 _ а) £ _

-

Ф = Ч> = f{xi,tj+0, при 1 ^ ТУ - 1,

, 0) = и0(х,),

а, = = А = ща + ад, д2 = ма +

Для погрешности г = у —и схемы (5) справедлива априорная оценка

из которой следует сходимость решения разностной задачи (5) к решению дифференциальной задачи (1)-(3) со скоростью 0(Л2 + г). Из полученной оценки следует справедливость теоремы.

Теорема!.. Пусть коэффициенты уравнений (1) и граничных условий (2) удовлетворяют следующим условиям: к(х, <) ^ с\ > 0, д ^

0, А- ^ Д» > 0,1 = ],2; %,*) € С3,1Ш, ?,/ € С2'1^). Тогда реше-

ние разностной задачи (1)-(3) сходится к решению дифференциальной задачи (5) в равномерной метрике со скоростью 0(Л2 + г).

Глава 2. Во второй главе рассматриваются нелокальные краевые задачи для уравнений параболического типа.

В области <3т — {(х, <) : 0 < х < /, 0 < I < Т] рассмотрена задача

о

- *(/,= &(*)«(/,*) + 03(1)01^(1,г) -

и(х,0) = и0(х), 0 ^ х ^ /,

(8)

где а е (0,1), ах е (О,/], к(х,1) > С1 > О, \д\, |&|, |/%| $ са.

Для решения задачи (11)-(13) справедлива априорная оценка

Игё + ЫЦд, ( 11/11^,+|К(х)||02 +1 {ц\{т) + ц\(г)) йт

(9)

Из оценки (9) следует единственность решения рассматриваемой задачи и непрерывная зависимость решения от правой части и начальных данных.

В замкнутой области 0т = {(¡М) : 0 ^ х ^ / , 0 ^ £ ^ Т} введена сетка

ш/,т ХШТ, й7Л = {х{ = гй,г = 0,1,... , Л^},

_ I Т

шт -- =зт,з =0,1,... Т7.т= —

N За

и задаче (6)-(8) поставлена в соответствие разностная схема

Уt

= АуЮ + ф,

(10)

Лу=

а 1 т

0.5/Н

1=0

Ау = {ау£)х - ¿у, г€о)Л,

I -- Ащ^уЕ К:?«-']-".-) *>■

а< = = /?х = ¿о + = <*лг + о^л. <*г = т/г, <г- посто-

янное число,

-, г = 0, г?г;

Л, 1 <: г ^ т — 1,

<р =

= + х = 0' уз = /(ж,-х = ик,

у(х{, 0) = ио(г,-)-

С помощью принципа максимума для решения разностной задачи (10) получена априорная оценка

Ну^Нс* < я ^(М<;<)1 + Ы*;<)1) + ||у°НсЛ + £ г||^'||Сл, Р.

из которой следует сходимость разностной схемы со скоростью 0(/г2+г) в равномерной метрике.

В области фт — {(аг, г) : 0 < х < 0 < I < Т} рассмотрена нелокальная краевая задача для нагруженного уравнения

ди _ д дЬ дх

г) \ т

+ £«*«(&,*)+ /(*.*). (П)

дХ { (12) ы(х, 0) = «о(аг), (13)

«^-постоянные числа, ^ сх > 0, 0 < £1 < £2 < • ■ - < £т <

фиксированные точки интервала (0,/).

Задачи типа (11)-(13) возникают при изучении переноса примеси вдоль русла рек.

Для решения задачи (11)—(13) справедлива оценка аналогичная оценке

Ы\1 + \Ы\1я, |||/1Цд( + 1М*)11о + I . (14)

Рассмотрен метод Ротэ, позволяющий свести решение нестационарной задачи (11)—(13) на каждом слое ^ к решению задачи для обыкновенного дифференциального уравнения

I

о

(15)

у(х, 0) = и0(х), к(х, ¿) ^ С1 > 0,

где

У1 = {у ~ у)/т, у — у1, у = у7-1, г = ТЦо- шаг сетки по времени.

Для решения задачи (15) имеет место дискретный аналог априорной оценки (14):

1М1о + £ Н(» + Л*- < м Е (и/иг + ^У) + г + М13 ¿'=1 )

(21)

Для погрешности г — у — и имеем задачу д { <9-2 \ 171

=ш (к{х'1) ш)+£ ак2&'+ 4 ' к=\

I

о

ф,0) = 0,

где ф = 0(т).

Применяя оценку (16) к решению задачи (17), находим

И^Но + + г)г'|1о ^ М ^^ ||^'||ог1

из которой следует сходимость метода Ротэ со скоростью 0(т). Для решения задачи (11)—(13) построена разностная схема

Ш

= Ау^+ф, у^ = егу + (1 — <т)у, у = у>+К (23).

А_ Ус о - Р1У0 , „ 1 ь „

АУ=\Лу=(аШ)х + У,

к=О

д+„ — № + ^ у .

У~ 0.5Л + ' '

к=1

0 =

я.Ч-Н - & & ~

+ глч+1—г—

у> =

л

А«!®

+ /о(<), ® =

0.5Л

у(а:,0) = и0(х), агб^л,

а; - = ¿¿_1/2(<), I - ^+1/2,

U, i = 1,2,... ,N — 1, аппроксимирующую задачу (11)—(13) с точностью до 0(h2 + тт°), где

Для решения разностной задачи (18) при сг = справедлива при малом г априорная оценка

lli^lß + £ ||(у + y)Ulr ^ М ( £ ll^'llgr + |М*)||3 ) • (19) i'=o \i'=o /

Из оценки (19) следует, что решение разностной задачи (18) сходится к решению дифференциальной задачи (11)-(13) со скоростью 0(h2 + 7"). На основании оценки (19) верна следующая теорема.

Теорема 2. Пусть выполнены условился, t) ^ ci > 0, |a/t|, \ßi\, ^ Со, 0 < fc < fc <...<£„> Тогда в классе гладких коэффициентов уравнения и граничных условиях при малом г ^ то и а = | справедлива априорная оценка

IIS + £ ||(у + У)г]\ог ^ М ( £ ||^'||02г + |K(*)|ß j'=о \j'=0

откуда следует сходимость со скоростью 0(h2 + г2) в норме ||у||ь где

IMIi = IIV+1II2+ t \\{у + у)^]\Ь, IMI2 = к и], [u,w]=

j'=0 t-o

Глава 3. В третьей главе работы рассматриваются нелокальные по

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

типа и обобщенного уравнения переноса с дробной по времени производной. Подобные задачи встречаются при изучении физических процессов с учетом эффекта памяти.

В области С?т = {(£)*) :0<г</, 0<*<Т} рассмотрена задача

Ж = ^ ~ я{х'*)и + °<Х<1>°<*<Т> (2°)

t

к(о, I) 1^(0, г) = I и(о, т)<1т, о <; / ^ т,

(21)

А(М)^(М) = /ММ) - А1(1), o<:t<:т,

и(х,0) = «о (г), О^г^/, (22)

где к(х,г) ^ С1 > 0, |9|,|/?| ^ с2.

Для решения задачи (20)-(22) имеет место априорная оценка

М20 + \Ы\1я, ¡ЦГ\\1я,+1 А»а(г)Л"+1М*)Ио| , (23)

означающая единственность решения задачи и непрерывную зависимость ее решения от входных данных.

Задаче (20)-(22) ставится в соответствие схема Ротэ

< X < I,

ду =£у(0^,)т, ]'=1

^ дх

у(х,0) = и0(х), О^.х^.1,

1о I >

для решения которой справедлив дискретный аналог оценки (23)

из которой следует сходимость метода Ротэ со скоростью О(т). Задаче (20)—(22) ставится в соответствие разностная схема

Ш = + Ф, (24)

\у" = ау + {1-а)у, у = у(аУ = У{*„Ь+1),

где

Ау - <

А_ а1Ух,о — 0, 5Мо2/о 1 ^ ,„ , , п

з'=1

0,5Л

Лу = {аух)х — ¿у, х = шн,

I

0,5/1

Ч>~ - /о — 0,

Д = ж?) + о, ьнйы, г =

у{хи 0) = ио(Х{). Для погрешности схемы (24) справедлива априорная оценка

У о

Н^НЗ + £ 11(5 + < М £ {4 + »1 + г. 3-0 }'=0

(25)

Из оценки (25) следует сходимость решения разностной задачи (24) к решению исходной дифференциальной задачи (20)-(22) со скоростью 0(к2 + г) в норме, стоящей в правой части (25).

- 15В области (¿т = {(£,<) : 0 < х < I, 0 < * < Т}, для обобщенного уравнения переноса рассмотрена задача

Й (Л(®.0^)-9(г.0« + /(®.0. О < х < /, 0<*<Т,(26)

(

ди

^ дх

г=0 Г(1

о

1_ [<0,т)с1т

-а) / (< - т)- ' (27)

-.сг-1

хЫ

= 0, (28)

0%~1и = чъ(х), (29)

где к(х,г) ^ с > 0, 0 < а < 1.

£>о{ы—дробная производная Римана - Лиувилля порядка а. Нелокальные условия (27) являются аналогом нелокального условия (21) в случае обобщенного уравнения переноса с дробной по времени производной (26).

Уравнения типа (26) встречаются при изучении задач фильтрации в сильно пористой (фрактальной) среде, деформационно-прочностных свойств полимерных материалов, физических процессов стохастического переноса.

Для решения задачи (26)-(29) справедлива априорная оценка

+ (30)

Из априорной оценки (30) следует единственность решения задачи (26)-(29).

В области С?т — {(я, t) : 0 < х < I, 0 <t <Т}, рассмотрена задача

= А -$(*,*)«+/(*,0. о < ж < /, о < * < т, (31)

ди к — кдх

1=0 ди * дх

V*-!

ад о,о, = 0,

Х = 1

ОоГ и = и0(х),

где к(х, г) О, 0 < а < 1.

Для решения задачи (31)-(34) справедлива априорная оценка

(33)

(34)

П^По < м Ш\\1я, + 1М*)11о + «о(0)) •

Для решения задачи (31)—(34) с регуляризованной дробной производной построена разностная схема

А= ЛуМ + ф,

у(х{, 0) = 0,

у(<7) = <ху+ (1 - (Г)у, у = (х,-, у = + 1),

где оператор Л определяется следующим образом

(35)

Лу =

а1У*,0 ~ 0, 5Ы0у0 ла 4 \ п У=-~0~5Ь--АщУ(°>Ь)> « = 0,

Лу = {ауц)х - ¿у, ж€и>л,

л+ _ + О.Шмн _ ,

Л 0.5А ' '

0 = у = А*«,

Для разностной схемы (35) верна априорная оценка

||^+1ПсЛ < д- +|д2(^")1) + н^нс, + £ г||^"цсЛ.

(36)

з'=о

Из оценки (36) следует устойчивость и сходимость решения разностной задачи (35) к решению дифференциальной задачи (31)—(34) со скоростью 0(Ь2 + г).

Публикации по теме диссертации

1. Худалов М.З., Юртин Й.И. Определение области загрязнения в процессах диффузии с реакцией // Сб. научных трудов "Нелинейные краевые задачи математической физики и их приложения". Институт математики HAH Украины. - Киев, 1998, с.225-227.

2. Рейдер Э.Е., Худалов М.З., Шхануков-Лафишев М.Х. Метод эквивалентной линеаризации в задачах сферически симметричной крио-деструкции биоткани // Сб. научных трудов "Нелинейные краевые задачи математической физики и их приложения". Институт математики HAH Украины. - Киев, 1999, с.203-205.

3. Худалов М.З., Шхануков М.Х. Об одной нелокальной задаче для уравнения теплопроводности.//Вестник СОГУ, 1999, N1, с 60-61.

4. Худалов М.З. О сходимости разностных схем для уравнения параболического типа с нелокальным условием // Тезисы докладов Второй международной конференции "Нелокальные краевые задачи и родственные проблемы биологии, информатики и физики". - Нальчик, 2001, с. 90.

5. Худалов М.З. Нелокальная краевая задача для нагруженного уравнения параболического типа // Владикавказский математический журнал. Т4., Вып.4, 2002г. с. 59-64.

6. Худалов М.З. Нелокальная задача для обощенного уравнения переноса с дробной производной по времени. Априорная оценка // Вестник СОГУ, серия - Естественные науки, 2002, N2.

7. Худалов М.З. Краевая задача для обощенного уравнения переноса с дробной производной по времени. // Материалы международного Российско-Узбекского симпозиума "Уравнения смешанного типа и родственные проблемы анализа и информатики", Нальчик - Эльбрус, 2003.

Сдано в набор 30.07.2003. Подписано в печать 01.08.2003. Гарнитура Computer Modern. Печать офсетная. Формат 60x84 1/16. Бумага офсетная. Усл. п.л. 0,95. Тираж 120 экз. Заказ №6.

Издательство Северо-Осетинского государственного университета имени K.JI. Хетагурова, 362025, г. Владикавказ, ул. Ватутина, 46.

Полиграфический центр Северо-Осетинского государственного университета имени К.Л. Хетагурова, 362025, г. Владикавказ, ул. Ватутина, 46.

P13 3 8 5

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Худалов, Марат Захарович

Введение.

Глава 1. Краевые задачи для уравнения параболического типа с дробной производной в граничных условиях

§1.1. Задача на полубесконечной полосе.

§1.2. Краевая задача для общего уравнения параболического типа.

§1.3. Априорная оценка.

§1.4. Построение разностной схемы.

§1.5. Погрешность аппроксимации.

§1.6. Устойчивость и сходимость разностной схемы.

Глава 2. Краевая задача для уравнения параболического типа с нелокальным условием.

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

§2.2. Априорная оценка.

§2.3. Построение разностной схемы.

§2.4. Погрешность аппроксимации.

§2.5. Устойчивость и сходимость разностной схемы.

§2.6. Нелокальная краевая задача для нагруженного уравнения параболического типа.

§2.7. Априорная оценка.

§2.8. Метод Ротэ.

§2.9. Разностная схема для нагруженного уравнения параболического типа.

§2.10. Порядок аппроксимации.

§2.11. Устойчивость и сходимость разностной схемы.

Глава 3. Краевая задача для уравнения параболического типа с нелокальным условием по времени на границе.

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

§3.2. Априорная оценка.

§3.3. Метод Ротэ.

§3.4. Построение разностной схемы.

§3.5. Погрешность аппроксимации.

§3.6. Устойчивость и сходимость схемы.

§3.7. Нелокальная задача для обобщенного уравнения переноса

§3.8. Априорная оценка.

§3.9. Нелокальная задача для обобщенного уравнения переноса с дробной по времени производной. Априорная оценка

§3.10. Построение разностной схемы.

§3.11. Погрешность аппроксимации.

§3.12. Устойчивость и сходимость разностной схемы.

 
Введение диссертация по математике, на тему "Разностные методы решения краевых задач для уравнений параболического типа с дробной производной в граничных условиях"

Использование дробных производных для описания и изучения физических процессов стохастического переноса стало в последние годы одной из популярных областей физики, многие проблемы фильтрации жидкости в сильно-пористых (фрактальных) средах приводят также к необходимости изучения краевых задач для уравнений в частных производных дробного порядка [34], [48], [50], [7], [18], [19], [20].

Сама идея обобщения понятия дифференцирования ^^ на нецелые р возникла с самого зарождения дифференциального исчисления, первые работы в этом направлении принадлежат Г. Лейбницу, Я. Бернулли, Л. Эйлеру и Ж. Фурье [58],[59],[64].

В 1832-1837гг. появляется серия работ Лиувилля [64]—[ТО], сделавших его по праву создателем теории дробного исчисления. В работах Б.Римана, Х.Хольмгрена, А. Летникова, А. Грюнвальда [62] было продолжено изучение производных любого порядка.

К первым работам по теории дифференциальных уравнений дробного порядка следует отнести работы L. O'Shaughnessy, S. Mandelbrojt [77],[75]. Задачу типа Коши для уравнения D%xy = f(x,y) рассмотрели E.Pitcher, W.E. Sewell в работе [72], в которой они доказали теорему существования и единственности решения рассматриваемой задачи. В дальнейшем эти результаты были значительно обобщены в работах M.A.Bassam, A.Z.A1. Abedeen, A.Z.A1. Abedeen, H.L.Arora [60]—[63], где получен ряд результатов, аналогичных теоремам из теории обыкновенных дифференциальных уравнений.

В работе А.М.Нахушева [31] изучена задача Штурма - Лиувилля для дифференциального уравнения дробного порядка тп

Ly = у"(х) а0(х)у'(х) + ak(x)D%*k (wk(x)y(x)) + am+i(x)y(x) = f(x), k=1

0 < аь < l,ao(x), am+i, ajfc(x), wk(x), к = 1,2,., m; f(x)— непрерывные на [0,1] функции, Dqx— оператор дробного в смысле Римана-Лиувилля дифференцирования порядка а. К задаче Штурма-Лиувилля

Ро2/'(°) + 0оу(О) = г0, Piy\l) + giy(l) = гь для этого уравнения редуцируются многие прямые и обратные задачи для уравнения смешанного гиперболо-параболического типа.

В работах [1, 2] Т.С. Алероев исследовал спектр задачи Дирихле для дифференциального уравнения и"{х) + a(x)Dgxu = f(x), 0 < х < 1, 0 < а < 1.

Им показано, что задача w(0) +/?u'(0) = it(l) = 0, /3 ^ 0, /(х) = 0, а(х) = А не имеет отрицательных собственных значений.

Еще раньше A.M. Нахушевым в [31] показано, что число Л является собственным значением задачи и(0) = 0, и( 1) = 0, а(х) = А для этого уравнения тогда и только тогда, когда А* является нулем функции Миттаг-Леффлера Е2-а,2(—А).

Ряд работ В.К. Вебера [11] -[16], М.И. Иманалиева, В.К. Вебера [22] посвящен исследованию систем дифференциальных уравнений дробного порядка. Работы [7, 8, 9,17, 52, 53, 54] посвящены построению и исследованию разностных схем для обыкновенных и в частных производных дифференциальных уравнений с дробной производной по времени и пространственной координате.

Бабенко Ю.А. в работе [6] для определения тепловых и диффузионных потоков на границе раздела двух сред применил метод расщепления оператора уравнения на два множителя, каждый из которых содержит производную порядка | по t (см. также [28]).

В монографиях [41],[33] дан достаточно полный обзор работ, посвященных дифференциальным уравнениям дробного порядка. Монография [33] посвящена качественно новым свойствам операторов дробного интегро 5 дифференцирования, и их применения к дифференциальным уравнениям дробного порядка.

В работе [7] рассмотрена краевая задача для уравнения параболического типа с дробной порядка \ производной по времени в граничных условиях. Методом сведения к интегральным уравнениям доказано существование и единственность рассмотренной задачи. Здесь же доказана с помощью принципа максимума сходимость разностной схемы в равномерной метрике.

Диссертационная работа посвящена построению и исследованию разностных схем для уравнений параболического типа с дробной производной в граничных условиях при любом а £ (0,1), где а— порядок дробной О производной.

В ней получены следующие результаты:

1. Для решения уравнения параболического типа общего вида с дробной производной по времени в граничных условиях получена априорная оценка, откуда следует единственность решения и непрерывная зависимость от правой части и начального условия.

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

О условиях.

3. Для решения нагруженного уравнения параболического типа с нелокальным условием получены априорные оценки в дифференциальных и разностных трактовках, откуда следует устойчивость и сходимость разностной схемы.

4. Для решения нелокальной по времени задачи для обобщенного уравнения переноса дробного порядка получена априорная оценка. Откуда следует единственность и непрерывная зависимость решения от входных данных.

5. На основе принципа максимума получены оценки в равномерной мет-t рике для решения разностной задачи, аппроксимирующей нелокальную задачу для обобщенного уравнения переноса дробного порядка, откуда следует сходимость разностной схемы в равномерной метрике.

Перейдем к более подробному изложению содержания диссертации.

 
Список источников диссертации и автореферата по математике, кандидата физико-математических наук, Худалов, Марат Захарович, Владикавказ

1. Алероев Т.С. Задачи Штурма - Лиувилля для дифференциального уравнения первого порядка с дробными производными в младших членах. Дифференц. уравнения 1982, т. 18, N2, с.341.

2. Алероев Т.С. Спектральный анализ одного класса несамосопряженных операторов. Дифференц. уравнения. 1984, т.20, N1, с 171-172.

3. Андреев В.Б. об одном методе численного решения третьей краевой задачи для уравнения параболического типа в р— мерном параллелепипеде. Сб. "Вычисл. методы и программирование", вып. 4. М.: изд-во МГУ, 1967, с. 64-75.

4. Андреев В.Б. О сходимости разностных с расщепляющимся оператором, аппроксимирующим третью краевую задачу для параболического уравнения. ЖВМ и МФ. 1969, т. 9, N2, с. 337-349.

5. Анохин Ю.А., Горстко А.В., Дамещек Л.Ю. и др. Математические модели и методы управления крупномасштабным водным объектом. Новосибирск: Наука, Сибирское отделение. 1987. 187с.

6. Бабенко Ю.И. Тепломассообмен. Метод расчета тепловых и диффузионных потоков. JL: "Химия", 1986, 144 с.

7. Березовский А.А., Шхануков М.Х., Керефов А.А. Краевые задачи для уравнения теплопроводности с дробной производной в граничных условиях и разностные методы их численной реализации. Укр. мат. журн. 1993, t.45,N9, с. 1289-1298.

8. Бечелова А.Р. О сходимости разностных схем для уравнения диффузии дробного порядка. Нелинейные краевые задачи математической физики и их приложения. Сб. научных трудов. Киев, 1996, с. 42-43.

9. Бечелова А.Р. О сходимости разностных схем, аппроксимирующих третью краевую задачу для уравнения диффузии дробного порядка. Нели78• нейные краевые задачи математической физики и их приложения. Сб. научных трудов. Киев, 1996, с 40-41.

10. Бэгли P.JL, Торвик П.Дж. Дифференциальное исчисление, основанное на производных дробного порядка новый подход к расчету конструкций с вязкоупругим демпфированием. Аэрокосмическая техника. 1984, т. 2, е2, с. 84-93.

11. Вебер В.К Асимптотическое поведение решений линейной системы дифференциальных уравнений дробного порядка. Исследования по интегро-дифференциальным уравнениям в Киргизии. Фрунзе: Илим,1983, вып. 16, с 119-125.

12. Вебер В.К. К общей теории линейных систем с дробными производными. Исследования по интегро-дифференциальным уравнениям в Киргизии. Фрунзе: Илим, 1985, вып. 18, с 301-305.

13. Вебер В.К. Линейные уравнения с дробными производными и постоянными коэффициентами в пространстве обобщенных функций. Исследования по интегро-дифференциальным уравнениям в Киргизии. Фрунзе: Илим, 1985, вып. 18, с.306-312.

14. Вебер В.К. Об одном дифференциальном уравнении нецелого порядка.Сб. трудов аспирантов и соискателей Кирг. ун-та. Сер. мат. наук. 1973, вып. 10, с. 7-14.

15. Вебер В.К. Пассивность линейных систем дифференциальных уравнений с дробными производными и квазиасимптотика решений. Исследования по интегро-дифференциальным уравнениям в Киргизии. Фрунзе: Илим, 1983, вып. 16, с. 349-356.

16. Вебер В.К. Структура общего решения системы у^ = Ау, 0 < а < 1. Сб. тр. аспирантов и соискателей Кирг. ун-та. Сер. мат. наук. 1976, вып. 11, с.26-32.

17. Гараперина И.Ю., Бечелова А.Р. Об одном итерационном методе решения первой краевой задачи для дифференциального уравнения с дробной производной. Вестник Кабардино-Балкарского государственного университета. Нальчик, 1996, вып. 1, с. 38-40.

18. Джарбашян М.М. Интегральные преобразования функций в комплексной области. М.: Наука, 1966, 671 с.

19. Джарбашян М.М. Краевая задача для дифференциального оператора дробного порядка типа Штурма-Лиувилля. Изв. АН Арм. ССР, мат. 1970, т. 5, N2, с. 71-97.

20. Джарбашян М.М., Нерсесян А.Б. Дробные производные и задачи Ко-ши для дифференциальных уравнений дробного порядка. Изв. Ан Арм. ССР, мат, 1968, т. 3, N1, с. 3-29.

21. Ильин В. А., Моисеев Е.И. Нелокальная задача для оператора Штурма-Лиувилля в дифференциальной и разностной трактовках

22. Иманалиев М.И., Вебер В.К. Об одном обобщении функции типа Миттаг-Леффлера и его применение. Исследование по интегро-дифференциальным уравнениям в Киргизии. Фрунзе: Илим, 1980, вып.13, с. 49-59.

23. Кочубей А.Ю. Диффузия дробного порядка. Дифференциальные уравнения. 1990, т.26 с. 660-670.

24. Кумыкова С.К. Об одной краевой задаче для уравнения sign\y\muxx + ууу = 0. Дифференциальные уравнения. 1976, т. 12, N1, с. 79-88.

25. Курант Р. Уравнения с частными производными. М.: Мир, 1964, 632 с.

26. Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973, 407 с.

27. Нахушев A.M. Обратные задачи для вырождающихся уравнений и интегральные уравнения Вольтера третьего рода. Дифференциальные уравнения. 1974, т. 10, N1, с. 100-111.

28. Нахушев A.M. Нагруженные уравнения, Дифференциальные уравнения. 1983, т. 19 N1 с. 86-94.

29. Нахушев A.M. Об уравнениях состояния непрерывных одномерных систем и их приложения. Нальчик, 1995, 50 с.

30. Нахушев A.M. Уравнения математической биологии. М.: Высшая школа, 1995, 301 с.

31. Нахушев A.M. Элементы дробного исчисления и их приложения., Нальчик, 2000г. -298с.

32. Нигматуллин P.P. Решение обобщенного уравнения переноса в среде с фрактальной геометрией. Phs.stat. Sol. b. 133. 1986.

33. Самарский А.А. Теория разностых схем. М.: 1989, с. 616.

34. Самарский А.А. Однородные разностные схемы на неравномерных сетках для уравнений параболического типа. Ж. вычисл. матем. и матем. физика., 1963, N2 266-298.

35. Самарский А.А. Аддитивные схемы. Доклад на Международном съезде математиков в Москве. Тезисы докладов, секция 14, вычислит, математика, 1966, с 46-47.

36. Самарский А.А., Гулин А.В. Устойчивость разностных схем. М.: Наука, 1973, 415 с.

37. Самарский А.А. Локально-одномерные схемы на неравномерных сетках. ЖВМ и МФ. 1963, т. 3, N3, с. 431-466.

38. Самарский А.А. Об одном экономическом разностном методе решения многомерного параболического уравнения в произвольной области. ЖВМ и МФ. 1962, т. 2, N5, с. 787-811.

39. Самко С.Г., Килбас А.А., Маричев О.И. Интегралы и производные дробного порядка и некоторые их приложения. Минск: "Наука и техника", 1987, 688 с.

40. Тихонов А.Н. Теоремы единственности для уравнения теплопроводности. Матем. сборник. 1935, т. 42, N2, с. 199-216.

41. Худалов М.З. О сходимости разностных схем для уравнения параболического типа с нелокальным условием. Тез. докл. Второй международной конференции "Нелокальные краевые задачи и родственные проблемы биологии, информатики и физики". -Нальчик, 2001, с.90.

42. Худалов М.З. Нелокальная краевая задача для нагруженного уравнения параболического типа. Владикавказский математический журнал. Т4., Вып.4, 2002г. стр.59-64.

43. Худалов М.З. Нелокальная задача для обощенного уравнения переноса с дробной производной по времени. Априорная оценка. Вестник СОГУ, серия Естественные науки, 2002, N2.

44. Худалов М.З., Шхануков М.Х. Об одной нелокальной задаче для уравнения теплопроводности. Вестник СОГУ, 1999, N1, с 60-61.

45. Худалов М.З., Юртин И.И. Определение области загрязнения в процессах диффузии с реакцией. // Сб. научных трудов "Нелинейные краевые задачи математической физики и их приложения". Институт математики НАН Украины, Киев: 1998, с.225-227.

46. Чукбар К.В. Стохастический перенос и дробные производные. ЖЭТФ. 1995, т. 108, вып. 5(11), с. 1875-1884.

47. Шефер Д., Кефер К. Фракталы в физике. Тр. б-го Междунар. симпоз. по фракталам в физике. (МЦТФ. Триес, Италия, 8-12 июня 1985). М., 1988, с. 62-71.

48. Шогенов В.Х., Кумыкова С.К., Шхануков-Лафишев М.Х. Обобщенное уравнение переноса и дробные производные. Дкол. АМАН. 1996, т.2, N 1, с. 43-45.

49. Шогенов В.Х., Шхануков-Лафишев М.Х., Бештов Х.М. Дробные производные: интерпретация и некоторые применения в физике. Препринт. Сообщения объединенного института ядерных исследований. Дубна, 1997, Р4-97-81.

50. Шхануков-Лафишев М.Х., Бечелова А.Р. Замечание к постановке краевых задач для дифференциальных уравнений с дробными производными. Сб. научных трудов. Киев, 1996, с. 286-287.

51. Шхануков-Лафишев М.Х., Бечелова А.Р. Численное решение третьей краевой задачи для обобщенного уравнения диффузии дробного порядка. Нелокальные задачи и родственные проблемы математической биологии, информатики и физики. Нальчик, 1996, с. 103.

52. Шхануков М.Х. О сходимости разностных схем для дифференциальных уравнений с дробной производной. Докл. РАН. 1996, т. 348, N6, с.746-748.

53. Эйлер Л. De progressionibvs transcendentibvs, sev qvarv termini generales algebrace dari neqvevnt /L.Eulero//Comment. Acad. Sci. Imperials petropolitanae. 1738. T. 5. P. 38-57.

54. Al-Abedeen A.Z. Existence theorem on differtial equations of generalized order// Rafidam J.Sci. Mosul. Univ. Iraq. 1976.Vol.l.P.95-104.

55. Al-Bassan M.A. On fractional calculus and its applications to the theory of ordinary differential equations of generalized order // Lect. Notes in Pure and Appl. Math. Dekker. New York, 1982. Vol. 80. P. 305-331.

56. Al-Bassan M.A. Some existence theorems on differential equations ofgeneralized order // Ibid. 1965. Bd 218. S. 70-78.

57. Fourier J. the Analytical Teory of Heat. N.Y. Dover pull., 1955. 466 p. (First publ.: Theorie Analytique de la Chaleur. A Paris: Chez firmin diclot peer et fils, 1822)

58. Grunwald A.K. Uber "bergenzte"Derivationen und deren Anwendung // Я. angew. Math, und Phys. 1867. Bd 12. S. 441-480.

59. Leibniz G.W. Leibnizan de l'Hospital// Oeuvres Mathematiques de Leibniz. Paris: Libr, de A. Franck, ed. 1853.p.l.Vol.2.P.297-302.

60. Liouville J. Memoire sur le calcul des differentielles a indices quelconques //Ibid. P. 71-162.

61. Liouville J. Memoire sur le changement de la variable independante dans le calcul des differentielles a indices quelquens // J^ 1'Ecole Roy. Polytechn. 1835. T. 15, sent. 24. P. 17-54.

62. Liouville J. Memoire sur le theoreme des fonctions complementaires // J. fur reine und angew. Math. 1834. Bd 11. S. 1-19.

63. Liouville J. Memoire sur l'integration des equations differentielles a indices fractionnaires // Ibid. 1837. T. 15, Na 55, P. 58-84.

64. Liouville J. Memoire sur l'usage Ton pent faire de la formule de Fourier, dans le calcul des differentielles a indices quelconques // Ibid. 1835. Bd 13. S. 219-232.

65. Liouville J. Memoire sur une formule d'analyse // Ibed. 1834. Bd 12, Ne 4, S. 273- 287.

66. Liouville J. Memoire sur quelques questions de geometrie et de mecanique, et sur un nouveau genre de calcul pour resoudre ces questions // J. l'Ecole Roy. Polytechn. 1832.T. 13,sect.21.P. 1-69.

67. Mandelbrojt S. Sulla generalizzazione del calcolo delie variazione // Atti Reale Accad. Naz. Lincei. Rend CI. sci. fis. mat. e natur. Ser. 6. 1925. Vol. 1. P. 151-156.

68. Pitchel E., Sewell W. Existence theorems for solutions of differntial equations of non-integral order. Ibid.1938. Vol.44,N 2. P. 100-107.

69. O'Shaughnessy L. Problem 433 // Amer. Math. Month. 1918. Vol. 25. P. 172-173.

70. Steinig John. The Changes of Sign of Fractions Integrals. Math. Z. 116. 1970. P.183-190.