О ЧИСЛЕННОМ РЕШЕНИИ НАЧАЛЬНО-КРАЕВЫХ ЗАДАЧ ДЛЯ УРАВНЕНИЯ КОНВЕКЦИИ-ДИФФУЗИИ С ДРОБНОЙ ПРОИЗВОДНОЙ КАПУТО И НЕЛОКАЛЬНЫМ ЛИНЕЙНЫМ ИСТОЧНИКОМ Текст научной статьи по специальности «Математика»

www.volsu.ru

МОДЕЛИРОВАНИЕ, ИНФОРМАТИКА И УПРАВЛЕНИЕ

DOI: https://doi.oгg/10.15688/mpcm.j'volsu.2020.4.4

УДК 519.63 Дата поступления статьи: 16.07.2020

ББК 22.16 Дата принятия статьи: 06.11.2020

О ЧИСЛЕННОМ РЕШЕНИИ НАЧАЛЬНО-КРАЕВЫХ ЗАДАЧ ДЛЯ УРАВНЕНИЯ КОНВЕКЦИИ-ДИФФУЗИИ С ДРОБНОЙ ПРОИЗВОДНОЙ КАПУТО

_ _______ о __л

И НЕЛОКАЛЬНЫМ ЛИНЕЙНЫМ ИСТОЧНИКОМ1

03

ГО

из

Аслан Мартинович Апеков

Кандидат физико-математических наук,

о

см

сЗ ведущий научный сотрудник отдела вычислительных методов,

рз Кабардино-Балкарский научный центр РАН

^ aslkbsu@yandex.ru

° https://orcid.org/0000-0002-6269-3717

§ ул. Шортанова, 89а, 360000 г. Нальчик, Российская Федерация

о

э

Мурат Хамидбиевич Бештоков

о Кандидат физико-математических наук, доцент,

§ ведущий научный сотрудник отдела вычислительных методов,

Кабардино-Балкарский научный центр РАН

X beshtokov-murat@yandex.ru

5' https://orcid.org/0000-0003-2968-9211

о ул. Шортанова, 89а, 360000 г. Нальчик, Российская Федерация

о Е-

э

со ЦЗ

Зарьяна Владимировна Бештокова

т Младший научный сотрудник отдела вычислительных методов,

щ Кабардино-Балкарский научный центр РАН

< zarabaeva@yandex.ru

@ ул. Шортанова, 89а, 360000 г. Нальчик, Российская Федерация

Замир Валериевич Шомахов

Кандидат физико-математических наук,

ведущий научный сотрудник отдела вычислительных методов, Кабардино-Балкарский научный центр РАН shozamir@yandex.ru https://orcid.org/0000-0001-5738-2626

ул. Шортанова, 89а, 360000 г. Нальчик, Российская Федерация

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

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

Введение

В последние годы возрос интерес к исследованию дифференциальных уравнений дробного порядка, в которых неизвестная функция содержится под знаком производной дробного порядка. Это обусловлено как развитием самой теории дробного интегро-дифференцирования, так и приложениями таких конструкций в различных областях науки [6; 12; 13]. Интегралы и производные нецелого порядка и дробные интегро-диффе-ренциальные уравнения находят множество применений в современных исследованиях в теоретической физике, механике и прикладной математике. Дробное математическое исчисление является мощным инструментом, который может быть использован для получения динамических моделей, в которых интегро-дифференциальные операторы по времени и координатам описывают степенную долгосрочную память и пространственную нелокальность сложных сред и процессов [6].

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

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

Краевые задачи для дифференциальных уравнений в частных производных с нелокальным (интегральным) источником возникают при изучении диффузии частиц в турбулентной плазме, переноса влаги в почвогрунтах, интегральный член возникает также при описании функции распределения по массам капель за счет микрофизических процессов конденсации, коагуляции, дробления и замерзания капель в конвективных облаках [2; 7; 9; 10].

Численным методам решения краевых задач для уравнения диффузии дробного порядка посвящены работы [1; 3; 4; 8; 11].

1. Постановка первой краевой задачи

В прямоугольнике = {(х,Ь) :0 < х < I, 0 < Ь < Т} рассмотрим первую краевую задачу для уравнения конвекции-диффузии с дробной производной Капуто порядка а и нелокальным линейным источником

д ( ди\ ди [х

д^и = \к(х,1)—)+ т(х,Ь) —--д(х,1)и — p(s,t)u(s,t)ds + /(х,1),

их у их у их Уд

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

и(0,Ь) = и(1,Ь) = 0, 0 < Ь < Т, (2)

и(х, 0) = и0(х), 0 < х < I, (3)

где

0 <Со < k(x,t),rx(x,t) < ci, \r(x,t)\, \q(x,t)\, \p(x,t)\, \kx{x,t)\< c2, (4)

д&и = щ—О) I (¡-X)«1 — дробная производная в смысле Капуто порядка а, 0 < а <

< 1, Ci = const > 0, г = 0,1, 2.

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

В работе будем использовать обозначения Mi = const > 0, i = 1, 2,..., которые зависят только от входных данных рассматриваемой задачи.

2. Априорная оценка в дифференциальной форме

Для получения априорной оценки решения задачи (1)-(3) в дифференциальной форме умножим уравнение (1) скалярно на u:

= ((kUx )x,u) + (rUx,u^ — (gu,u^ -(/ p{s,t)u(s,t)ds,u) + (^f,uj, (5)

J 0

где = /0 abdx, = ||a||0, где a,b — заданные на [0,/] функции.

Преобразуем интегралы, входящие в тождество (5), пользуясь неравенством Коши с е, леммой 1 [1], получим

(datu,u) > 1 (i,öo>2) = 1 д&ыц, (6)

(^(ких)х,и^ = J и(киdx = ukux:\l0 — J ku^dx, (rUx,uj = J ruxudx = 2ru210 — rxU2dx < 1 rw2|0 — 2||w|| — Puds,и) < 1||w||0 + Puds)

'0 < Ml

2 10 2 ./0 x " 2 1 _ 1 / / Гх \2 2""110 ■ 2V ^0

el rx 1

I u2dsdx + 1||w||0 < М2Ци\\2.

(7)

(8)

00

— ^ qu,u^j = — J qu2dx < с2Ци\\0,

{f,u)

2

I fudx < 2INI0 + 21 Учитывая преобразования (6)-(11), из (5) находим

|2 I o„ ll„. 112 i „ ll„.l|2 ^ o„./„„. i „,„.2|Z

2 + 2cb|ux|0 + сьМ0 < 2ukux + ru210 + М3Ци\\0 +

0

Из (12) с учетом (2) получаем

№110 + N10 + К\\2 <м, и 2 + м5\

Применяя к обеим частям (13) оператор дробного интегрирования а, находим

N о + ОДМ1 адо < м6одН12 + м7(од\/\\2 + \kWWo).

На основании леммы 2 [1] из (14) находим априорную оценку

N о + ОДН1 адо < м (од\/\\2 + \М*)\\о),

(9) (10) (11)

(12)

(13)

(14)

(15)

где D^u

udt

—^f —>_

Г( a) J (i-т)1-« 0

— дробный интеграл Римана — Лиувилля порядка a,

0 < a < 1; М = const > 0, зависящая только от входных данных (1)-(4). Теорема 1. Если к(х, t) E Cl'°(QT), r(x, t), q(x, t), p(x, t), f(x, t) E С(QT), u(x, t) E E C2'0(QT) П Сl'0(QT), d0atu(x, t) E С(QT) и выполнены условия (4), тогда для решения задачи (1)-(3) справедлива априорная оценка (15).

Из оценки (15) следуют единственность и непрерывная зависимость решения задачи (1)—(3) от входных данных.

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

На равномерной сетке whT дифференциальной задаче (1)-(3) поставим в соответствие разностную схему порядка аппроксимации 0(h2 + т2) :

А,

0

З+а

У = К fexa)) + bi 3a3iyx,i + b+al+lу{°} — diy(u) —

^ p^h + ф' , (x, t) E WhtT, (16)

s=0

a

И Н А Уо ) = У¥ = 0,

у(х, 0) = щ(х),

(17)

(18)

где У = п-,) £ — дискретный аналог дробной производной Капуто по-

8=0

рядка о, 0 < а < 1, обеспечивающий высокий порядок точности 0(т3-а) [8].

_(а,а) _ _1—а (а,а) ¿п - а , Li]

(l + а)1 а- (l - 1 + а)1 а, I > 1

2 - а

а)2- а - (I - 1 + а)

2- а

(I + а)1-а + (1 - 1 + а)

1

1> 1,

А (а, а) (а, а)

при j = 0, с0 ) = а0 );

■а(а,а) + ^ s = 0,

при J> 0, = ^ а^ + bi+f - Ъ^, 1 < s < j - 1,

а^ — ^ 8 = ^

о> = к(хг-о.5,«+°), V = , Ф = /(^'П а =1 — 2,

1 „ 0.5%|

к = --—, Л = —---разностное число Реинольдса,

1 + л К

г0 = г(0,г) = г(0+и) < 0, гм = г(1,г) = > 0, = р(х

Л а,а)

1 — а

> (s + а)-а > 0, у

(а)

ау-

j+I

+ (1 - а)уj, dl -d(xг,t'+а).

Априорную оценку найдем методом энергетических неравенств, для этого введем скалярные произведения и норму:

N-1

N

(u,v) = h, (u,v] =^2uiVih, {и,и) = \\и\\0.

i=1

i=1

Умножим теперь (16) скалярно на у(а) :

(л0%+ау, У(а)) = (к(ауХа))х, У{а)) + {Ь-аУ^, У{а)) +

+ (ь+а^1)уХа), у(а)) — (а), у(а)) — Р1у[а)к, + (Ф, у(а)). Преобразуем суммы, входящие в тождество (19), с учетом (17) и леммы 1 [8]

(Ао%+.y, У(а)) > ^ А

ot

'j+ а

(19)

(20)

(х(ау^)х, у(а)^ = хау(^а)у(а)

-(ах(-1), (у ¡а))2] <-(ах-х, уФуЫ

, (ку(%

(1 + h М1)

(а) (а) акх, уУх ;у( )

(а) 2

ах, (уХ ))

(21)

1

1

2

i

2

а

1

— (Va), у(а)) < С2 | | y(<j) | | 2; (22)

(£ Piy(sa)h, У(а)) < 2 | | У(&) 11 2 + (2, (£ P^(a)h)1 < M^Hg, (23) =0 2 2 =0

(ф, < 2 | |у(а) | | 2 + 2 11 ф 11 2. (24)

Принимая во внимание преобразования (20)-(24), из (19) находим

(2, А0\.+а(У2)) + (1+hМ)) (ак, (j^)2] < —(axx, У^У^] + (b~ay(^)) +

2 (1 +

+ (b+a(+l) ^, У(о)) +Мз||у(^)||2 + 2| |ф| |2. (25)

Из (25) находим

А0%+а | | У | | 2 + M4 Л0 < — (a*x, yxa)y(a)} + (b~a, y^y^) +

+ (б+а(+1) у^, уМ) +М5\\у(^)\\0 + М6\\ф\\2. (26) Преобразуем первое, второе и третье слагаемые в правой части (26). Тогда получим

-(ах„ уГУ*)] + (Ь~а, у^у^) + (ь+а(+1)у^, у^) < £\\+ М?\\у^\\0. (27)

Учитывая (27), из (26) получаем

\ \У\ \ 0 + М4\\у^а)]|2 < е\\у^а)]|0 + Ме%\\у°\\0 + Мб\\ф\\о. (28)

Выбирая в последнем е = -М, из (28) получаем

А0%+, \ \У\ \ 0 + \ \< М9\\Л \0 + М1о\\ф\\2. (29)

Перепишем (29) в другой форме

А0\.+а | | у | | 0 < ЩЛу^Х + ВДИЮ + Мю| |ф| |2. (30)

На основании леммы 7 [3] из (30) получаем

| | yi+l | | 0 <м(| |у°||0 + max ||ф'| , (31)

\ 0<3'<3 I

где М = const > 0, не зависящая от h и т.

Теорема 2. Пусть выполнены условия (4), тогда существует такое т0, что если т < т0, то для решения разностной задачи (16)-(18) справедлива априорная оценка (31).

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

Пусть и (ж, ¿) — решение задачи (1)-(3), у(х1, ) = у3 — решение разностной задачи (16) —(18). Для оценки точности разностной схемы (16)-(18) рассмотрим разность zj = = у3 — и3, где и3 = и(х1, tj). Тогда, подставляя у = г + и в соотношения (16)-(18), получаем задачу для функции г

= К' (^4а)) + №.. + 45 - - Е Р^ + ^, (32)

^ у x.i —Г

8=0

^ = z(va) = 0, (33)

z(x, 0) = 0, (34)

где Ф = 0(V + т2^ — погрешности аппроксимации дифференциальной задачи (1)-(3)

разностной схемой (16)-(18) в классе решений и = и(х, t) задачи (1)-(3).

В силу линейности (1)-(3) и (32)-(34), применяя априорную оценку (31) к решению задачи (32)-(34), получаем неравенство

| | | | 2 <Мтах.|ф'||g, (35)

0<i '<]

где М = const > 0, не зависящая от h и т.

Из априорной оценки (35) следует сходимость решения разностной задачи (16)-(18) к решению дифференциальной задачи (1)-(3) в смысле нормы | |zi+l | | 0 на каждом слое так, что существует такое т0, что при т < т0 справедлива оценка

^■+1 -^+i||o < M(h2 + т2).

4. Постановка третьей краевой задачи и априорная оценка в дифференциальной форме

Рассмотрим теперь третью краевую задачу для уравнения (1)

{

к(0, t)ux(0, t) = ß1(i)u(0, t) - H1(i), -k(l, t)ux( l, t) = ß2(t)u(l, t) - Ц2СО,

(36)

где

0 < со <k< ci, | ß 1, $2,r,q, p, rx, kxl < C2. (37)

Умножим уравнение (1) скалярно на u:

,и) = {jkUx)x,^j + (rUx,u) qu,uj p(s, t)u{s, t)ds,u) + (f,uj. (38)

о

Преобразуем второе слагаемое в правой части (38)

/ \ [1 с? [1 [1 с?

(rux,u) = ruuxdx < u2dx + е / u2xdx < е||ux||2 + т2|lull0. (39)

v ' Jo 4е J0 Jo 4е

Учитывая преобразования (7)—(11), (39), из (38) находим

\д,Ы1 + со\\их\\2о < ик (х, г)их(х, г)\0 + £\М0 + М{\\и\\1 + М2\\Д1 (40) Оценим первое слагаемое в правой части (40)

иких\0 = k(l, t)ux(l, t)u(l, t) - к(0, t)ux(0, t)u(0, t) = u(l, t)(Ц2СО - $2(t)u(l, t)^ +

+u(0, t)[Mt) - $i(t)u(0, t)j = -$2(t)u2(l, t) + Ц-2(t)u(l, t) - $i(t)u2(0, t) + Mt)u(0, t) <

<M3(u2(0, t)+u2(l, t)) +2 (v-2(t) + ^(t)) < г\К\\20 + Ml\\u\\l + 1 (£(1) + ^2(t)). (41)

Учитывая (41), из (40) при e = получим

№112 + \\ux\\0 < M5\\u||0 + M6(||f\\20 + tf(t) + v-2(t)). (42)

Применяя к обеим частям неравенства (42) оператор дробного интегрирования D—а, получаем

\\ull2 + D-*\\ux\2 < m5d-«\\u||2 + m7(d-«01 f\\2 + v2(t) + vl(t)) + \Mx)\\2). (43) На основании леммы 2 [1] из (43) находим априорную оценку

\\ull 2 + D-*\\ux\2 < M (р-* (|| л\2 + At) + v2(t)) + lMs)\\2), (44)

где M = const > 0, зависящая от входных данных задачи (1), (36), (3). Теорем а 3. Если k(x,t) E С 1,2(QT), f(x, t), q(x, t), р(х, t), f(x, t) E С (QT), u(x, t) E E C2'2(Qt) П С1,2{Qt), datu(x, t) E С(QT) и выполнены условия (4), (37), тогда для решения задачи (1), (36), (3) справедлива априорная оценка (44).

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

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

На равномерной сетке whT дифференциальной задаче (1), (36), (3) поставим в соответствие разностную схему порядка аппроксимации 0(h2 + Т2)

г

.. ,= I » . АКsr Ми, ,J

= К[Ша)) + Ъ-*4у£} + КЧ+1 yxs - - Е PsVs^h + ф^, (45)

\ у x,i —■*

8=2

КХИ ¿а2 = V iv2u) + 0.5hp2y2u) + 0.5hA«.+ a У2 - 0.25h2 р2у2а) - Vi, tE шх, x = 0, (46)

г N

-KNaNy^l = в2УN + Е Psvffh + 0.5hAatj+a Vn - E P&^h - V2, tE x = I, (47)

8=2 8=2

у(х, 0) = щ(х), X G wh,

(48)

где

ß 1^+a) = ß^j+a) + 0.5Ы0, ß 2^+a) = ß2(iJ+a) + 0.5hd?N,

j+a) = Hl(ij+a) + 0.5h^ö, Ц2(^ j+a) = Ц2 (¿j+a) + 0.5h^N.

Перепишем (45)-(48) в операторной форме

!

У = Л( tj+a )У (a) + Ф, у(х, 0) = щ(х), X G Wh,

да

0

(49)

где

Л(* j+a) y(a)

Луг(а) = к(ау£Л + b-ау{°] + Ь+а(+1)y(a) - dy(a) - £ p^(a)h, г = 1, N - 1,

^ 'x s=0

аух

- (a)_ к0а1 ¿q - ß iy0a) - 0.25h2p0y0a)

л-yi;

, = 0,

л+ (a)

0.5h

-KNaNyiXVN - ß2У(м] - 0.5h £ Piy(a)h v =----, i = N,

N

Ф =

Ф = фг, г = 1, N — 1, 2

ф = |ß'1, ^ = 0,

ф+ = rß2, г = N, h

к =

0.5h f 1

к = + 0.5h|r| , 1 + к 1к

К0 = 1 | 0.5h| г о | , ко.5

к N = 1 | 0.5h|rN | км-0.5

, Г0 < 0,

rv - 0,

t* = ¿j+1/2.

Умножим (49) теперь скалярно на y(a)

Д0\.+ ,У, y(a)] = ^(ij+a)y(a), y(a)] + [ф, y(a)

(50)

N

r 1 N % % 10.5h, г = 0,N

где m,i; ^ w^.h, h = < . [и, u\ = [1, u2\ = |[u\|g, (u, v\ = l^UiVih.

i=0 yh,i = 0,N, i=1

Оценим суммы, входящие в (50):

Д01 j+a У ^ У

(a)

1

> -> 2

1, Д0%+а (У2)

(51)

A(tj+a) y(a), y(a)] = (A(ij+a)y(a), y(a)) +0,5hy0a^-y0a) +0,5hy{N^+y = (K(ayXa))x, y(a)) + (Va^, y(a)) + (W^ y(a), y(a))-

-(dy(a),y(a)) - p^a)h,+ K0a1 yia0y0a) - ß 1(y0a))2-

N

-0, 25h2p0?/0a) - KNavyXvNy{N - ß2(у^)2 - 0, 5h£ p^h =

s=0

(.ayx°\ + (б-a, yxW) + (b+a(+1), у(а)у(а)) -

(О ,

ZtfA у(а)

s=2

- V 1(У2и))2 - V2(уN )

(52)

Преобразуем слагаемые в правой части (52)

^М\ ] J_fhJ°) „м

(ау К + (b-a, yxW) +

у^УП =

= -\ак(-1), (у)2

+ 1 +1),^(О) х + hM:

( О ( а)

«К, yx Г )

+ (б-а, yx°]y(a)) +

ка

(У О)2} + e\\ у^)]\2 + МЦ[у (а)]|2,

(53)

(а), у(а)

Е Piy^h, У(а)

s=2

- в 1 (у2а))2 - в2(УN])

(а)\2 N

d, (у(а))2} - вг(у2О))2 - в2(уN)2 -

Е p^h, у(а)

s=2

<

< e\\yf Щ + Mei[y(а)]|2 +

Е Piy^h^j

< e\\у{°]]I2 + М11[у(а)]122. (54)

Учитывая (53), (54), из (52) находим

Mt3+a)y(а),у(а)] < - МА\\+ 2e\\+ ме|[у(о)]|2.

ф у(а)] = (Ф, у(а)) + 0.5hу2а)ф- + 0.5hyN^+ = [ф, У(а)] + rf + ц2yN <

(55)

< e\\+ MeI[y(а)]|2 + M7[|[ф]|2 + Hi2 + Ц*2). Принимая во внимание преобразования (51)-(56), из (50) находим

(56)

Д0%+,|М|2 + \\Л < eM8\\+ м^у(а)]|2 + MW(\M2 + H1 + H-2). (57)

Из (57) при e

1

2 М8

получим

л0%+.|[у]^ + \ \у^]И < Мц\[У^^ + Ми(|[ф]Ц + н2 + и|).

!2 + \\ у^

На основании леммы 7 [3] из (58) находим априорную оценку

(58)

(

| [y3+1] И <М | [у°] !2 + max | [ф'] ^ + ц2 + н2 ,

2<j'<j

где М = const > 0, не зависящая от h и т.

(59)

Теорема 4. Пусть выполнены условия (4), (37), тогда существует такое т, что если т < т0, то для решения разностной задачи (45), (48) справедлива априорная оценка (59).

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

Пусть и(х, ¿) — решение задачи (1), (36), (3), у(х^, Ь]) = у3 — решение разностной задачи (45)-(48). Для оценки точности разностной схемы (45)-(48) рассмотрим разность = у1 — и], где и] = и(хг, Ь]). Тогда, подставляя у = г + и в соотношения (45)-(48), получаем задачу для функции г

г

A0\+^ = К' Ша)) + Ь~3a\zx,i + К3а>+1 z^ - ^ - Z P^(a)h + Ф*', (x, t) G wh>x,

\ ' хЛ —„

8=0

(60)

к0ах= (i4a) + 0.5hAOt.+az0 - 0.25fr2p04a) - Vb t G шт, ж = 0, (61)

N

(o)_ = в r(a) + 0 5h Aa - ^ P> rH

=0

z(x, 0) = 0, (63)

где Ф = 0(h2 + т2), V1 = 0(h2 + т2), V2 = 0(h2 + т2) — погрешности аппроксимации дифференциальной задачи (1), (36), (3) разностной схемой (45)-(48) в классе решений и = u(x, t) задачи (1), (36), (3).

Применяя априорную оценку (59) к решению задачи (60)-(63), получаем неравенство

|Н+1]|0 < М0max. (|ф'Щ + vi'2 + vf) , (64)

где М = const > 0, не зависящая от h и т.

Из априорной оценки (64) следует сходимость решения разностной задачи (45)-(48) к решению дифференциальной задачи (1), (36), (3) в смысле нормы |[zJ+1]|0 на каждом слое так, что существует такое т0, что при т < т0 справедлива оценка

-xnünz^I = ß 2 4 + 0.5hA«.+a zNPlz^h - V2, t G шт, ж = I, (62)

ly+i -uJ+i]|0 < m(V + т2).

6. Алгоритм численного решения

Для численного решения (1), (36), (3) приведем разностную схему (45)-(48) к расчетному виду. Тогда уравнение (45) приводится к следующему виду

ЛуГ! - Сгу3+1 + вгу+ = -F3 , г = 1, N - 1, (65)

где

Аг = тся^а3 - Thab-3 аг, Вг = TaK3ia3i+1 + Thab++3 аг+1,

T1-aia,a)

Сг = Аг + Вг + h2 V °° , + Tah2 £, г г г Г(2 - а) г'

т1-а i-1 г

?3 — Л Л.J nn.J I RB.J I U2_1_ V^ Ia,ah„,s+1

F3 = ААгУ3-1 - ССгу3 + В В гУг+i + h2T^ - h2 Г^у Е 4-7 (^ - tf) - £ Р^ h,

^ ' s=0 s=0

AAi = т(1 - a)Kia^i - Th(1 - а)b- 3 ai} BBi = т(1 - а)к^а{+1 + Th(1 - а)Ь+а+

С С = AAi + BBi - h

Краевое условие (46) принимает вид

Т1-а Г(а,а) ; Т с2

Г(2 - а)

+ т(1 - a)h2d3.

У2 = К1У1 + Ц1,

(66)

где

К1

таж2а1

т ак2а{ + ah т|3 [ + 0.5h2 ^(2-

а)

Ц1

\i1hT - (1 - a)hт|31У32 + т(1 - а)х2Сц(у{ - у2) + 0.5h

Т1-а г(а,а) ; Т с2

Г(2 - а)

У 2 - 0.2ЬЬ2Р2У,0-

-0.5h2

Т1-а 3-1

Т(2 - а)

ЛГЛу 0+1 -

s=2

так2а\ + аhт(3 \ + 0.5h

Т1-а Ла,а) 1 ; Т с2

Г(2 - а) _

Краевое условие (47) принимает вид

yN = К yN -1 + Ц2,

(67)

где

К =

таж3а3

так3а3 + ahт(32 + 0.5h2 -(2-

а)

Ц2

N

H2hT- (1 - а)Ы&2yN - т(1 - а)кхах(lfN - yN-1) - Р*тЛh + 0.5h

s=2

Т1-а г(а,а) i Т c2

Г(2 - а)

N

-05K2^-г V с™ (уN+1 - У%)

Г(2 - а) ^ 3-3

=2

так3а33 + аЫ$32 +

h2

Т1-а г(а,а) 1 I on

2 Г(2 - а) _

Таким образом, если в (45)-(48) суммы вычислять с нижнего слоя, тогда разностная схема (45)-(48) приводится к трехдиагональной системе линейных алгебраических уравнений (65)-(67), решение которой легко находится известным методом прогонки.

7. Результаты численного эксперимента

Коэффициенты уравнения и граничных условий третьей краевой задачи (1), (36), (3) подбираются таким образом, чтобы точным решением задачи была функция u( x, ) = = t3еx.

Ниже в таблице при различных значениях параметров а = 0,01; 0,5; 0,99 и уменьшении размера сетки приведены максимальное значение погрешности ( = - u) и порядок сходимости (ПС) в нормах Ц]^ и \\ ■ \\c(whT), где \\y\\c('u,hT) = max

(xit tj )ewhi

когда

h = т. Погрешность уменьшается в соответствии с порядком аппроксимации 0(h2 + т2).

Порядок сходимости определяется по следующей формуле: ПС= log hi 77—л-, где Zi —

h2 |[ Z2^2

это погрешность, соответствующая h .

Изменение погрешности и порядка сходимости в нормах |[ ]|0 и \\ ■ \\с(гиНт) при уменьшении размера сетки при различных значениях а = 0, 01; 0, 5; 0, 99, когда к = т

а h max |[2J]|о 0<]<m ПС в |[.]|о Wz\\c (whT) ПС в \\ ■ W C(whT)

0,01 1/10 0,011251484 0,014635951

1/20 0,002769085 2,0226 0,003672308 1,9948

1/40 0,000686142 2,0128 0,000919450 1,9978

1/80 0,000170728 2,0068 0,000230015 1,9990

1/160 0,000042578 2,0035 0,000057522 1,9996

0,5 1/10 0,023164160 0,030433575

1/20 0,005793221 1,9995 0,007617871 1,9982

1/40 0,001448032 2,0003 0,001902556 2,0014

1/80 0,000362003 2,0000 0,000474990 2,0020

1/160 0,000090513 2,0000 0,000118602 2,0018

0,99 1/10 0,030896591 0,040483502

1/20 0,007714381 2,0018 0,010121342 1,9999

1/40 0,001926257 2,0018 0,002528356 2,0011

1/80 0,000481216 2,0010 0,000631654 2,0010

1/160 0,000120260 2,0005 0,000157832 2,0007

Замечание. Основной результат работы заключается в доказательстве априорных оценок для решения первой и третьей начально-краевых задач для уравнения конвекции диффузии с дробной производной Капуто и нелокальным линейным источником интегрального вида как в дифференциальном, так и разностном виде. Полученные неравенства означают устойчивость решения относительно входных данных. В силу линейности рассматриваемых задач эти неравенства позволяют утверждать сходимость приближенного решения к точному (в предположении существования последнего в классе достаточно гладких функций). Задачи с нелокальным источником вида (1)-(3) возникают при исследовании формирования микроструктуры конвективных облаков. В частности, при описании микрофизических процессов коагуляции (объединение мелких диспергированных частиц в большие по размеру агрегаты), взаимодействия капель и кристаллов, дробления и замерзания капель.

ПРИМЕЧАНИЕ

1 Исследование выполнено при финансовой поддержке РФФИ и ГФЕН Китая в рамках научного проекта № 20-51-53007.

СПИСОК ЛИТЕРАТУРЫ

1. Алиханов, А. А. Априорные оценки решений краевых задач для уравнений дробного порядка / А. А. Алиханов // Дифференциальные уравнения. — 2010. — № 46 (5). — C. 660-666.

2. Ашабоков, Б. А. Конвективные облака: численные модели и результаты моделирования в естественных условиях и при активном воздействии / Б. А. Ашабоков,

A. В. Шапавалов. — Нальчик : Изд-во КБНЦ РАН, 2008. — 252 с.

3. Бештоков, М. Х. К краевым задачам для вырождающихся псевдопараболических уравнений с дробной производной Герасимова — Капуто / М. Х. Бештоков // Известия вузов. Математика. — 2018. — № 10. — C. 3-16.

4. Бештоков, М. Х. К краевым задачам для интегро-дифференциальных уравнений дробного порядка / М. Х. Бештоков, Ф. А. Эржибова // Математические труды. — 2020. — № 23 (1). — C. 16-36.

5. Самарский, А. А. Теория разностных схем / А. А. Самарский. — М. : Наука, 1983. — 616 с.

6. Тарасов, В. Е. Модели теоретической физики с интегро-дифференцированием дробного порядка / В. Е. Тарасов. — Ижевск : Изд-во Ижев. ин-та компьютер. исследований, 2011. — 568 с.

7. Численное моделирование облаков / Е. Л. Коган, И. П. Мазин, Б. Н. Сергеев,

B. И. Хворостьянов. — М. : Гидрометеоиздат, 1984. — 186 с.

8. Alikhanov, A. A. A new difference scheme for the time fractional diffusion equation / A. A. Alikhanov // Journal of Computational Physics. — 2015. — № 2805. — P. 424-438.

9. Berry, E. X. An analysis of gloud drop growth by collection / E. X. Berry, R. L. Reinhardt // J. Atmos. Sci. — 1974. — № 31 (7). — P. 1825-1831. — DOI: 10.1175/1520-0469(1974)031<1814:AA0CDG>2.0.C0;2.

10. Berry, E. X. Cloud droplets growth by collection / E. X. Berry // J. Atmos. Sci. — 1967. — № 24 (6). — P. 688-701.

11. Beshtokov, M. Kh. Difference Methods for Solving Local and Nonlocal Boundary Value Problems for a Loaded Fractional Order Heat Equation. / M. Kh. Beshtokov, M. Z. Khudalov // Stability, Control and Differential Games. Lecture Notes in Control and Information Sciences — Proceeding. — Cham : Springer Nature, 2020. — P. 187-201. — DOI: 10.1007/978-3-030-42831-0.

12. Miller, K. S. An Introduction to the Fractional Calculus and Fractional Differential Equations / K. S. Miller, B. Ross. — New York : Wiley, Wiley and Sons, 1993. — 376 p.

13. Oldham, K. B. The Fractional Calculus. Theory and Applications of Differentiation and Integration to Arbitrary Order / K. B. Oldham, J. Spanier. — New York : Academic Press, 1974. — 240 p.

REFERENCES

1. Alikhanov A.A. Apriornye otsenki resheniy kraevykh zadach dlya uravneniy drobnogo poryadka [A Priori Estimates of Solutions to Boundary Value Problems for Fractional Order Equations]. Differentsialnye uravneniya, 2010, no. 46 (5), pp. 660-666.

2. Ashabokov B.A., Shapavalov A.V. Konvektivnye oblaka: chislennye modeli i rezultaty modelirovaniya v estestvennykh usloviyakh i pri aktivnom vozdeystvii [Convective Clouds: Numerical Models and Simulation Results Under Natural and Active Conditions]. Nalchik, Izd-vo KBNTs RAN, 2008. 252 p.

3. Beshtokov M.Kh. K kraevym zadacham dlya vyrozhdayushchikhsya psevdoparabolicheskikh uravneniy s drobnoy proizvodnoy Gerasimova — Kaputo [To Boundary-Value Problems for Degenerating Pseudoparabolic Equations with Gerasimov — Caputo Fractional Derivative]. Izvestiya vuzov. Matematika, 2018, no. 10, pp. 3-16.

4. Beshtokov M.Kh., Erzhibova F.A. K kraevym zadacham dlya integro-differentsialnykh uravneniy drobnogo poryadka [To Boundary-Value Problems for Integro-Differential Equations of Fractional Order]. Matematicheskie trudy, 2020, no. 23 (1), pp. 16-36.

5. Samarskiy A.A. Teoriya raznostnykh skhem [Theory of Difference Schemes]. Moscow, Nauka Publ., 1983. 616 p.

6. Tarasov V.E. Modeli teoreticheskoy fiziki s integro-differentsirovaniem drobnogo poryadka [Models of Theoretical Physics with Fractional Order Integro-Differentiation]. Izhevsk, Izd-vo Izhev. in-ta kompyuter. issledovaniy, 2011. 568 p.

7. Kogan E.L., Mazin I.P., Sergeev B.N., Khvorostyanov V.I. Chislennoe modelirovanie oblakov [Numerical Simulation of Clouds]. Moscow, Gidrometeoizdat, 1984. 186 p.

8. Alikhanov A.A. A New Difference Scheme for the Time Fractional Diffusion Equation. Journal of Computational Physics, 2015, no. 2805, pp. 424-438.

9. Berry E.X., Reinhardt R.L. An Analysis of Gloud Drop Growth by Collection. J. Atmos. Sci, 1974, no. 31 (7), pp. 1825-1831. DOI: 10.1175/1520-0469(1974)031<1814:AA0CDG>2.0.C0;2.

10. Berry E.X. Cloud Droplets Growth by Collection. J. Atmos. Sci, 1967, no. 24 (6), pp. 688-701.

11. Beshtokov M.Kh., Khudalov M.Z. Difference Methods for Solving Local and Nonlocal Boundary Value Problems for a Loaded Fractional Order Heat Equation. Stability, Control and Differential Games. Lecture Notes in Control and Information Sciences — Proceeding. Cham, Springer Nature, 2020, pp. 187-201. DOI: 10.1007/978-3-030-42831-0.

12. Miller K.S., Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations. New York, Wiley, Wiley and Sons, 1993. 376 p.

13. Oldham K.B., Spanier J. The Fractional Calculus. Theory and Applications of Differentiation and Integration to Arbitrary Order. New York, Academic Press, 1974. 240 p.

ON THE NUMERICAL SOLUTION OF INITIAL-BOUNDARY VALUE PROBLEMS FOR THE CONVECTION-DIFFUSION EQUATION WITH A FRACTIONAL CAPUTO DERIVATIVE AND A NONLOCAL LINEAR SOURCE

Aslan M. Apekov

Candidate of Physical and Mathematical Sciences, Leading Researcher, Department of Computational Methods, Kabardino-Balkarian Scientific Center of RAS aslkbsu@yandex.ru

https://orcid.org/0000-0002-6269-3717

Shortanova St, 89A, 360000 Nalchik, Russian Federation

Murat Kh. Beshtokov

Candidate of Physical and Mathematical Sciences, Associate Professor,

Leading Researcher, Department of Computational Methods,

Kabardino-Balkarian Scientific Center of RAS

beshtokov-murat@yandex.ru

https://orcid.org/0000-0003-2968-9211

Shortanova St, 89A, 360000 Nalchik, Russian Federation

Zaryana V. Beshtokova

Junior Researcher, Department of Computational Methods, Kabardino-Balkarian Scientific Center of RAS zarabaeva@yandex.ru

Shortanova St, 89A, 360000 Nalchik, Russian Federation

Zamir V. Shomakhov

Candidate of Physical and Mathematical Sciences, Leading Researcher,

Department of Computational Methods,

Kabardino-Balkarian Scientific Center of RAS

shozamir@yandex.ru

https://orcid.org/0000-0001-5738-2626

Shortanova St, 89A, 360000 Nalchik, Russian Federation

Abstract. In a rectangular domain the first and third initial-boundary value problems are studied for the one-dimensional with respect to the spatial variable convection-diffusion equation with a fractional Caputo derivative and a nonlocal linear source of integral form. Using the method of energy inequalities, under the assumption of the existence of a regular solution, a priori estimates are obtained in differential form, which implies the uniqueness and continuous dependence of the solution on the input data of the problem. On a uniform grid, two difference schemes are constructed that approximate the first and third initial-boundary value problems, respectively. For the solution of the difference problems, a priori estimates are obtained in the difference interpretation. The obtained estimates in difference form imply uniqueness and stability, as well as convergence at a rate equal to the order of the approximation error. An algorithm for the approximate solution of the third boundary value problem is constructed, numerical calculations of test examples are carried out, illustrating the theoretical results obtained in this work.

Key words: initial boundary value problems, a priori estimation, convection-diffusion equation, fractional order differential equation, fractional Caputo derivative.