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

Владикавказский математический журнал 2024, Том 26, Выпуск 3, С. 5-23

УДК 519.64

DOI 10.46698/l0699-2536-6844-a

НАЧАЛЬНО-КРАЕВЫЕ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ВЛАГОПЕРНОСА С ДРОБНЫМИ ПРОИЗВОДНЫМИ РАЗНЫХ ПОРЯДКОВ И НЕЛОКАЛЬНЫМ

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

М. Х. Бештоков1

1 Институт прикладной математики и автоматизации КБНЦ РАН, Россия, 360000, Нальчик, ул. Шортанова, 89 а E-mail: beshtokov-murat@yandex.ru

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

Ключевые слова: уравнение влагопереноса, краевые задачи, дробная производная Герасимова — Капуто, нелокальный источник, разностные схемы, устойчивость и сходимость разностных схем.

AMS Subject Classification: 65N06, 65N12, 65R20.

Образец цитирования: Бештоков М. Х. Начально-краевые задачи для уравнения влагоперноса с дробными производными разных порядков и нелокальным линейным источником // Владикавк. мат. журн.—2024.—Т. 26, вып. 3.—C. 5-23. DOI: 10.46698/l0699-2536-6844-a.

1. Введение

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

Wt =

здесь W — влажность в долях единицы, х — глубина, Ь — время, D(W) — коэффициент диффузивности почвенной влаги.

© 2024 Бештоков М. Х.

Диффузионная модель, предполагающая, что если в начальный момент задана неравномерная влажность, то должен возникнуть поток влаги из более влажных в менее влажные слои, часто не оправдывается. Объяснить, когда и при каких условиях происходит движение влаги в прямом и обратном направлении, можно с помощью модифицированного уравнения диффузии (или уравнения влагопереноса, уравнения Аллера [2]):

Wt = )ШХ + ЛШхг)х, (*)

где А — положительная постоянная.

Вопросы моделирования переноса влаги в почвогрунтах [1, 2], течения жидкости в трещиновато-пористых средах [3, 4], движения подземных вод со свободной поверхностью в многослойных средах [5, 6], теплопроводности в двухтемпературных системах [7] и течения некоторых неньютоновских жидкостей [8] связано с необходимостью исследования краевых задач для уравнения вида (*).

Дифференциальные уравнения дробного порядка служат хорошим математическим аппаратом для более точного описания систем и процессов природы, для которых необходимо учитывать предысторию (память, наследственные свойства) процесса. Характеристиками, позволяющими учитывать память процесса в таких уравнениях, являются функции памяти, которые представляют собой ядра интегралов, определяющих операторы дробного интегро-дифференцирования. К примеру, такой функцией памяти может быть степенная функция. Показатель степенной функции памяти определяет порядок производной и связан с фрактальной размерностью среды, в которой протекает исследуемый процесс [9-12].

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

Так, в работе [13] предложены и исследованы математические модели водного режима в почвогрунтах с фрактальной структурой. В основе этих моделей лежат дифференциальные уравнения влагопереноса с дробной по времени производной. Математически обосновывается переход от уравнения с целочисленной производной к уравнению с дробной производной в смысле Герасимова — Капуто. Для задачи Коши выписывается единственное представление решения.

Численным методам решения различных краевых задач для дифференциальных уравнений дробного порядка посвящены работы [14-19].

В настоящей работе исследуются начально-краевые задачи для уравнения влагопере-носа дробного порядка с нелокальным линейным источником и переменными коэффициентами. Основным методом исследования является методом энергетических неравенств. При предположении существования регулярного решения для каждой из рассмотренных первой и третьей начально-краевых задач получена оценка в дифференциальной форме, откуда следуют единственность и непрерывная зависимость решения от входных данных исходной задачи. На равномерной сетке в соответствие каждой дифференциальной задаче ставится соответствующая разностная схема порядка аппроксимации 0(Н2 + т2) при а = в и 0(Н2 + т2-тах{«,в>) при а = в, где 0 < а, в < 1 — порядки дробных производных Герасимова — Капуто. Для каждой разностной задачи в предположении существования решения получена оценка в разностной форме, из чего следуют единственность и устойчивость решения разностной задачи по правой части и начальным данным. В силу линейности рассматриваемых начально-краевых задач полученные оценки в разностной

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

Настоящая работа является продолжением серии работ автора в этом направлении [19-24].

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

В замкнутой области С}т = : 0 ^ ж ^ 0 ^ £ ^ Т} рассмотрим следующую

задачу для обобщенного уравнения влагопереноса:

/ . ,ди\ а д ( . . ди\ . ,ди I

-д(х, Ь)и(х, Ь) + J р(£, Ь)и(£, Ь) ^ + /(х, Ь), 0 < х < I, 0 < Ь < Т,

(1)

о

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

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

где

0 < Со < к(х,Ь), п(х) < С1, |д(х,Ь)|, |г(х,Ь)|, |гх(х,Ь)|, |кх(х,Ь)|, |р(х,Ь)| < С2, (4)

= г(1-7) /о "¿-т'ут ~~ дробная производная в смысле Герасимова — Капуто порядка 7, 0 < 7 < 1.

Отметим, что приведенная выше конструкция дробной производной введена [25] итальянским механиком М. Капуто в 1967 г. Поэтому за рубежом ее называют дробной производной Капуто. Однако правильнее называть ее дробной производной Герасимова — Капуто, так как еще в 1948 г. советский механик А. Н. Герасимов рассматривал [26] подобные выражения:

t

1 [ ит(х,т)

Г(7) У (Ь - т)т

^т.

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

Также отметим, что рассматриваемое нелокальное уравнение влагопереноса, ввиду присутствия интеграла по пространственной переменной и производной дробного порядка в смысле Герасимова — Капуто, относится к классу нагруженных уравнений [27].

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

Обозначим через Мг (г = 1, 2,...) положительные постоянные числа, зависящие только от входных данных исходной задачи.

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

(д^и, и) = ((ких)х , и) + (п(х)их)х , и) + (ти,х(х, г), и)

-ым«,) , (5)

где ( и,у) = /0 uvdx, (и, и) = ||и(-,г)\\0 = ||и|\0, где и, V — заданные на [0,1] функции.

Лемма 1. Для всякой функции и(х), обращающейся в нуль при х = 0, справедлива оценка

/2

н н2 ^ н н2 IIIIо ^ "^"Н^жНо-

< Действительно, в силу условия и(0, г) = 0 справедливо неравенство

х

и(х,г) = J их(х,г) dx. 0

Возведем в квадрат обе стороны, тогда на основании неравенства Коши — Буняковского получим

2

х

00 Проинтегрируем обе части по х от 0 до /, тогда по формуле Дирихле получим:

I I / х

■2(х г) йх ^ I I х I и1(е, г) I йх

00

Jи2(х,г) йх ^ у |х у и2((,г)

I I I

/2 1 /2

— / и2(х, ¿) (1х — — / (хих(х,Ь)) и2(х,1)г1х.

0 0 0

Следовательно, \\и\\$ ^ >

Лемма 2. Для всякой функции у(х), заданной на равномерной сетке

йь = {хг = гН, г = 0,..., N Н = //М, х0 = 0, хN = /}

и обращающейся в нуль при х = 0, справедлива оценка

\|у|\2 < /2\Ы12.

< Действительно, в силу условия у(0,г) =0 справедливо неравенство

Уг = Уо + Ух,зН = ^ Ух'зИ.

',в — / J ух,в' 8=1 8=1

х

Возведем в квадрат обе стороны, тогда на основании неравенства Коши — Шварца получим

(г \ 2 i г г

8=1 / 8=1 8=1 8=1 Умножим теперь обе части на Н и просуммируем по г от 1 до N

N

N

N N

2

г=1

г=1 V 8=1

г=1 8=1

I2 I

N

"%Н < г2||ух ]|Х

8=1

для любого Н € [0,1]. Следовательно,

11у||0 < 12Пух]|Х- >

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

4 4

«-х)х ••) = / » (ких)х = и^Ю -I *

(6) (7)

4 4

дм (пих)х , и) ^ У идд/: (пих)х Йх = пид0в/:их |0 - J п(х)ихд0itих ^х

< щд^ихФоАих)2 <1х

(8)

I

(гих,и) = У гиих(1х ^ ^ J и2 йх + е J их с1х ^ М1(е)\\и\\1 + е\\и: 0

22 ю + е||их ||о,

00 — (^(х, Ь)и, и) ^ С2|и|Х, \ 1 1

р(£,Ь)и(£,Ь) ^У и У р(£,Ь)и(£,Ь) Й^х

0 0

( / 4 1

2'

2

р(е,Ь)и(е,Ь) ^

/

4 4

+ J <1х < М3||и||§,

00 I

/1 1

2 11 — 11о 1- 2 П./ по-

(9) (10)

(11)

(12)

Учитывая преобразования (6)—(12), из (5) c учетом (2) находим

I

^dotlHIo + J v(ux)2 dx + c0\\ux\\l ^ e\\ux\\l + M3(e)\\u\\l +

(13)

Выбирая e = Щ-, из (13) получаем

да ЦиУо + düt f n(nx )2 dx + \\ux ||2 < M4\\u\\0 + Mb\\f Ю-

(14)

I случай. Пусть а > в, тогда, применяя к обеим частям неравенства (14) оператор дробного интегрирования ^й", получаем

!\u\\o + D-(a-ß)\\ux\\0 < M6D0-°

\u\2 + M7 (D-ta\\f \\0 + \\uo(x)\2)

(15)

где D0tYu =

1 ft udr

_ r t г_

Г(7) J0 (i-r)!-7

yot

0 <Y < 1.

На основании [17, леммы 2] из (15) находим оценку

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

u\\? < M8(D0ta\\f \\0 + \\uo(x)ПО),

(16)

2

где \ u\ 1

\u\\2 + D0t

— (а-в) 11 i|2 \ \ ux \ \ 0'

II случай. Пусть а = в, тогда, применяя к обеим частям неравенства (14) оператор

(17)

дробного интегрирования получаем

\|и|\2 + \\их |\2 < 1|и|\2 + М10 \\0 + 1|ио(х)\2 + ||и0(х)\0

На основании [17, леммы 2] из (17) находим оценку

u\\W2i(0,o < M11 f D0ta\f\\2 + \Mx)\\W2i(0,

(18)

II и 2 11 112 II II 2

где 11u\\(0,1) = \\u\0 + \\ux\\0-

III случай. Пусть а < ß, тогда, применяя к обеим частям неравенства (14) оператор

дробного интегрирования D0f, получаем

u

x\ 0

0 + D¿t(в-в)ПuП0 < M12D~ßHu^ + M13 (D~tß\\f \0 + \\u0(x)\

В силу леммы 1 из (19) получаем

\\ux\\2 + D0t(ß ^МЮ < M14D-f \\UxП0 + M15 (D-f \\f\\2 + \\u0(x)\2 На основании [17, леммы 2] из (20) находим оценку

\\u\\2 < M16 (D~ß\\f \\2 + \\u0(x)\2

Ii ii2 11 Ii2 1 (ß~а) и n2

где \\u\2 = WuxE + D0t 'Wu^.

(19)

(20)

(21)

Справедлива следующая теорема:

Теорема 1. Если к(х^) € С1>°{<^т), ф) € С1[0,1], г{х, д(х, р(х, ¡{х, €

с{Ят), и{х,г) е с2-0 (дг) п с1,0 (дг), д^и{х,г) е с(дг), д^ихх{х,г) е с(дт) и выполнены условия (4), то для решения задачи (1)-(3) справедливы оценки: (16) в случае, когда а > в; (18) в случае, когда а = в; (21) в случае, когда а <

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

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

Для приближенного решения задачи (1)—(3) применим метод конечных разностей. В замкнутой области введем равномерную сетку шьт х шт, где шь = {хг = г/г,

г = 0,..., /V, Л, = 1/Щ, шт = {I, = эт, э = 0,1,.. .,эо, т = Т/^}.

На равномерной сетке шит дифференциальной задаче (1)-(3) поставим в соответствие разностную схему порядка аппроксимации О (Л2 + т2) при а = в и 0(й2 + т2-тах(«,в}) при а = в:

А&,+„У = Х? (а?у^)+ (7.У,).,. + 6-а?+ 6+?а^уЦ

N

- ^У(СТ) + Е ^У^ + , (х, Г) € шнт, (22)

У0СТ) = У^ = 0, (23)

у(ж, 0)= зд(ж), (24)

где А = г(2-7) Уг ~ дискретный аналог дробной производной Гераси-

мова — Капуто порядка 7, 0 < 7 < 1, обеспечивающий порядок точности 0(т3-7) при а = 1 - и 0(т2-7) при о- = 0.5 [18],

а07'ст) = а1-7, а(7'ст) =(1 + а)1-7 - (1 - 1 + а)1-7 , 1 ^ 1,

5(7,^ =

1

г 2-7

при j = 0

при j > 0

(/ + а)2"7 -(1-1+ а)2"7 - - (/ + а)1"" + (/ - 1 + а)1"7

= а(7>^). с0 = а0 ;

а07'ст) + б17'ст), 5 = 0;

с^ = I а!7,ст) + б17+^) - б!7'ст), 1 < 5 < j - 1;

1 ^ 1,

-8 ^8 ~ "«+1 "8 (7 = 1 — ^, при а = [3 и и = 0.5 при а ф /3

а(7'ст) - 6?7'ст), 5 = j,

2 I 5 = 0, Ж,

= к(Хг-0.5,^+а), а(+1) = 7г = Г]{Хг-0.5), = Г , ,

, Г? +<г )

$ = f (хг,гз +а), й1 = ц^г^з+а ), р\ = р(х^ tj+a),

г(х,гз +а) = Г+(х,3+а) + Г-(х,3+а), 1г(х,гз +а)| = Г+(х^з+а) - г-(х,^+а), Г+(х,гз +а) = 0.5(г(х,гз+а) + 1г(х,гз +а)|) ^ 0, Г-(х,гз +а) = 0.5(г(х,гз+а) - 1г(х,гз +а)|) ^ 0,

= Щх^з+г) = ~ Разностное число Рейнольдса.

Оценку решения задачи (22)-(24) найдем методом энергетических неравенств, тогда введем скалярные произведения и норму в виде:

N-1

г=1

N

(и, V) = ^ иг(и, V] = ^ и^гН, (и, и) = (1,и2)

г=1

= 1, и =

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

Д-,+„ у,у(а^ = (х(аух^) х ,у(т)) + (д^ (7ух)х ,у(т)) + (ъ-ау<х'\у + (ь+а(+1)уха),у(а)) - (йу(а),у(а)) + ^£рву(а)Я,у(а^ + ($,у(а)) .

(25)

С учетом [18, леммы 1] преобразуем суммы, входящие в тождество (25), тогда получим

1

(а> -А:

2 1

2

Ю)

(26)

(а)) у,(а)\ = хау(а) у(а)

)х,у^ ) = хаух 'у ах(-1), №)2

N

,(а)

< - (ах«^у(а)

аух % (ху(а))« 1

(а) (а)

ах«, ух у( )

(1 + НМ1)

^ ^)2

(27)

<,+„ (7ух)х ,у(а)) = 7у(а)<,+„ух^ - Гт,у(<а)<,+„(ух)

<

^ о 0-

■,7+ст

х] |0;

(28)

а),у(а^ £

С2

0

N

N

^у^ п,у(а)) <

+ Ы

<

чв=0 N

N

2'^

в=0 в=0

чв=0

N N

в=0 в=0 21

о + 2

Из (25) с учетом (26)-(30) получим

(29)

(30)

д0%+ст I

0-

+

3+0

у(а ,у(а)

■хло < - (ахх,у«а)у(а)

+ (ь-ау«а) ,у(а))

о ^ 2

0

2

1

1

2

2

0

0

2

2

0

Оценим первое, второе и третье слагаемые в правой части (31), тогда получаем

, Л+Л+П»^) »-112 . лж чИ ,(ст) М2

(ст) (ст)

ох^Ух У( )

+ (Га,уХст)у(ст)) + (W+^yH) < е||уХст)] |0 + M5(£)||y(CT) (32)

С учетом (32) при е = Щ^- из (31) находим

||у||§ + Ы1о + ЦУЙо < м6Цу(-) II2 + М7Ы1

Перепишем последнее в другой форме:

Aj IMl0 + <,+„ 11Ух]|§ < M8 |yj+1|0 + MCT ||yj ||0 + М7ЦИЮ. (33)

I случай. Пусть a > ß, тогда на основании [22, леммы 7] из (33) получаем неравен-

Ъ -L .DW

|yj+1< Ms (1И10+j 1И' 110) • (34)

II случай. Пусть a = ß, тогда в силу [22, лемма 7] из (33) получаем неравенство

Ну5'+1|W2i(o,i) < Ms( |y0|W2i(o,i) + 0тах '||2), (35)

^e m2wim = \\y% + \rn\i

III случай. Пусть a < ß, тогда в силу (23) и [22, лемма 2, лемма 7] из (33) получаем неравенство

НУ!+1]|0 < Мэ(||УХ||0 + j ^'НО • (36)

Теорема 2. Пусть выполнены условия (4), тогда существует такое малое то = T0(c0,c1, c2,a, о), что если т ^ т0, то для решения разностной задачи (22)-(24) справедливы оценки: (34) в случае, когда a > ß; (35) в случае, когда a = ß; (36) в случае, когда a < ß.

Исследуем вопрос о разрешимости задачи (22)-(24). Для этого рассмотрим однородную задачу (^ = 0, u0(x) = 0), решение которой заведомо существует

AjУ = X (ajУ^)х>. + Aj (Yiyx)x,i

+ ajyg + aj+iyX;? - djУ^ + E Р,У<ст)h, (x, t) € ыЛ>т, (37)

N

ПуИ + П ^ - ^У^ + ^ Р У(0

«=0

Уо^ = У* =0, у(х, 0)=0. (38)

Пусть у(х^) — одно из решений однородной задачи (37)-(38). Из неравенства (34) следует, что ||у^+1М2 < 0, из (35) — ||у^+1М2 < 0, а из (36) — ЦуХ+1у0 < 0, но ||у^+1М2 = 0 и 11Ух+11|0 = 0 лишь при у = 0, х € гу^,. Поэтому из оценок (34) при а > в, (35) при а = в, (36) при а < в следуют, что единственным решением однородной задачи (37)-(38) является у = 0. Тем самым решение задачи (22)-(24) при любых ио(х) существует и единственно.

Таким образом, из оценок (34) при а > в, (35) при а = в, (36) при а < в следуют единственность и устойчивость решения задачи (22)-(24) по начальным данным и правой части.

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

Д ^г = х] (*>¿г)^ + (Ъг,)хЛ + Ъ-!]гЦ + Ъ+!а^

N

— 4 ^ + ^ Р**^ П + Ф! > (х>*) е , (39)

zi+1\\2w!(п/) < M max. \Ф'II2- (43)

4a) = 4°° = 0, (40)

z(x, 0) = 0, xeoJh, (41)

где Ф = Oh2 + т2) при a = ß и Ф = ö(h2 + т2-max{a>ß}) при а = ß — погрешности аппроксимации дифференциальной задачи (1)-(3) разностной схемой (22)-(24) в классе решении u = u(x,t) задачи (1)-(3).

В силу того, что задача (39)-(41) линейна, тогда

I случай. Пусть a > ß, тогда применяя оценку (34) к решению задачи (39)-(41), получаем неравенство

\\zi+1\\2 < M max \Ф'\\2- (42)

11 110 11 110

II случай. Пусть a = ß, тогда применяя оценку (35) к решению задачи (39)-(41), получаем неравенство

\Ui+i|

\wl(o,i)^ ^ ojj 11" "о

III случай. Пусть a < ß, тогда применяя оценку (36) к решению задачи (39)-(41), получаем неравенство

||4+11|2 < М max 1|фЛ|2 (44)

x o^j ' i1 110

Из полученных оценок (42)-(44) следует сходимость решения разностной задачи (22)-(24) к решению дифференциальной задачи (1)-(3) так, что существует такое т0, что при т ^ т0 справедливы оценки:

1) в случае, когда a > ß: \\yi+1 - ui+1\\0 < M (h2 + т2-maxKß}) ;

2) в случае, когда a = ß: \\yi+1 — ui+1\\Wi(0 i) ^ M (h2 + т2);

3) в случае, когда a < ß: \\yl+1 - 4+1]|o ^ M + , где M — const > 0, не зависящая от h и т.

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

Второе краевое условие в (2) заменим условием третьего рода, тогда вместо условия (2) имеем

/<•№0=°. (45)

[- n(1,t)= ß(t)u(l,t) - ß(t),

где

\ß\ < C2, П(х, t) = k(x, t) ux + nd0atux- (46)

Из (5) с учетом (6)-(12) находим

I

1 С 1

dotMÎ + ft J r](ux)2dx + colluxllo < uU(x,t)fQ + e\\ux\\2 + Mi(e)||u||§ + - ||/||§. (47)

Оценим первое слагаемое в правой части (47) с учетом (45)

иП(х, t) |0 = n(l, t)u(l, t) = u(l, t) (p(t) - e(t)u(l, t)) = -e(t)u2(l, t) + ^(t)u(l, t)

1 1 (48)

< M2u2(l,t) + -n2{t) < M3(e)\\u\\2 + e\\ux\\2 + -/x2(i).

Из (47) с учетом (48) при е = у находим

I

ayullg + d0e У n(ux)2dx + ||ux 112 < M4^uy§ + M5 (У/112 + At)) . (49)

о

I случай. Пусть a > в, тогда, применяя к обеим частям неравенства (49) оператор дробного интегрирования D0ta, получаем

l|u||0 + D0t(a-e)||ux||2 < MzDoaHufo + M7{potT (||/1|0 + At)) + ||uo(x)||2)• (50)

На основании [17, леммы 2] из (50) находим оценку

||u||2 < M8(D-ta (||/1|2 + Ц2(t)) + ||uo(x)||2)• (51)

II случай. Пусть a = в, тогда применяя к обеим частям неравенства (49) оператор дробного интегрирования D0ta, получаем

||u||2 + ||ux||2 < M9D0ta|u12 + Mio(Do"ta (||/1|0 + ^(t)) + ||uo(x)|2 + ||u0(x)||0)• (52)

На основании леммы 2 [17] из (52) находим оценку

||u||Wi(o,i) < Mn(D0ta (||/1|0 + At)) + ||uo(x)||Wi(0Л) • (53)

III случай. Пусть a < в, тогда применяя к обеим частям неравенства (49) оператор дробного интегрирования D0te, получаем

||ux||0 + D0t(e-a)||u|2 < Mi2D0te||u||0 + Mi3 (Dof (||/||0 + At)) + ||u0(x)||0) • (54) В силу леммы 1 из (54) получаем

||ux||0 + D-t(e-a)||u||o < Mi4D-f ||ux|0 + Mi5 [D-f (||/1|2 + At)) + ||u0(x)||0) • (55)

На основании леммы 2 [17] из (55) находим оценку

||u||2 < Mi6 (yDôf (||/1|2 + At)) + ||u0(x)|2) • (56)

Теорема 3. Если k(x,t) € C1,0(Ôr), ф) € С^ОД q(x,t),f(x,t) € C(QT), u(x,t) € C2>° (Qt) n c1'0 (QT), d&u(x,t) € C(QT), d%tuxx(x,t) € C(QT) и выполнены условия (4), (46), то для решения задачи (1), (3), (45) справедливы оценки: (51) в случае, когда a > в; (53) в случае, когда a = в; (56) в случае, когда a < в.

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

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

На равномерной сетке дифференциальной задаче (1), (3), (45) поставим в соответствие разностную схему порядка аппроксимации 0(Л2 + т2) при а = в и 0(Л2 + т2-тах{«,в>) при а = в:

Л&,+„у = Х- (4у^) ж>г + Л^, (т)^ + Г а-уй + а^УЙ

N

- <у(ст) + ЕР^Я + , (х, *) € , (57)

«=о

у(0,£) =0, £ € Шт, х = 0,

(58)

XNЗДуЦ^ + Л0^+(Т Ыу*^^ = + УN -Е- х = Ъ

(59)

у(х, 0) = ио(х), х £ шь, (60)

где в7 = в7' + 0.5Л—, Д7 = Д7 + 0.5Л—•

Найдем оценку методом энергетических неравенств, для этого умножим уравнение (57) скалярно на у. Тогда, учитывая (25)-(31), после нетрудных преобразований получаем

+ М2 ||у(ст) ]|0 + (Г а- у£> ,у(ст)) + (&+' а^ ,у(ст))

(61)

где (и,и] = «^¿Я, Я = ] ^^ 2 ^' (и, и] = (1, и2] = 11и]|§, (и,и) = Е^=11 «^¿Л.

I Я, 2 = 0, N,

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

х^ + Д&

у

И

N

= [ХмаУ{х% + (7»Ух,лг)) Ум

3 + &

р + 0.5Л— - вуЙ° - 0.5Му^ - 0.5Л ( Л«^УN - ЕЯ

N

«=0

И

у*

<

-0.5Лу£> А^ум + Мз(в) + вЦу^]!2

2 N

+ М^р2 + 0.5Лу^— - 0.5М (у^) + 0.5Лу^ Е

«=0

Учитывая (62) при е = Щ1, из (61) находим

N

- (&+7' а7+1у^ ,у(ст)

+ Е Я,у(ст)

ч«=0

(62)

0

(66)

Из (63) после несложных преобразований с учетом неравенства Коши с е получим

Л j\\yfo + Ш20 + Ua)]\l < М6\\у^]\1 + М7 (|М|0 + /X2) • (64)

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

Д j\\у}\20 + Д j< Mi\\v3+% + М9СТ||У]|2 + М7 (Ml + /х2) . (65)

I случай. Пусть a > ß, тогда на основании [22, лемма 7] из (65) получаем

ЦУ+Ч12 < MW(||y0]|2 + max (|М|2 + а) .

0 V 0 O^j'^j J

II случай. Пусть a = ß, тогда на основании [22, лемма 7] из (65) получаем

lk'+1]lW(o,0 < Mo (l|y0]|W2i(0,i) + 0jj M + M2)) , (67)

где \\y°}\w^(o,i) = Wy°fo + UfoIII случай. Пусть a < ß, тогда из (65) с учетом леммы 2 получаем неравенство

A^J\yf0 + < Мп||^+1]1о + МШТо + М15 (М§ + А • (68)

На основании [22, лемма 7] из (67) находим

ы ]0 < Mi^ ||ух^]|2+j (|kj' ]|2+M2) ). (69)

Теорема 4. Пусть выполнены условия (4), (46), тогда существует такое малое т0 = т0(с0,c1, c2,a,a), что если т ^ т0, то для решения разностной задачи (57)-(60) справедливы оценки: (66) в случае, когда a > ß; (67) в случае, когда a = ß; (69) в случае, когда a < ß.

Из оценок (66) при a > ß, (67) при a = ß, (69) при a < ß следуют единственность и устойчивость решения задачи (57)—(60) по начальным данным и правой части.

Пусть u(x,t) — решение задачи (1), (3), (45) y(xi,tj) = yj — решение разностной задачи (57)-(60). Для оценки точности разностной схемы (57)-(60) рассмотрим разность zj = yj — uj, где uj = u(xi,tj). Тогда, подставляя y = z + u в соотношения (57)-(60), получаем задачу для функции z :

A j z = yj (ajzf) ^ + Aj (7iZ,)Xtt + b-jajzg + b+j j^}

n (70)

— djz(a) + Y, Psz(sa)h + , (x, t) € uhtT,

s=0

z(0,t)=0, tewT, x = 0, (71)

— (xn aN z%N + Aj (YN z,,n )) = ßzN +0.5^ Aj zn —S^ PSz'si) ^ —v, x = l, (72)

z(x, 0) = 0, xeuh, (73)

где Ф = O(h2 + т2), V = O(h2 + т2) при а = в и Ф = O(h2 + т2-max{«>e}), v = 0(h2 + r2-max{a,e^ при а = в — погрешности аппроксимации дифференциальной задачи (1), (3), (45) разностной схемой (57)-(60) в классе решении u = u(x,t) задачи (1), (3), (45).

В силу того, что задача (70)-(73) линейна:

I случай. Пусть а > в, тогда применяя оценку (66) к решению задачи (70)-(73), получаем неравенство

+1]|0 < M^ax j']|2 + v2) . (74)

II случай. Пусть а = в, тогда применяя оценку (67) к решению задачи (70)-(73), получаем неравенство

1И+1]|W21(0,1) < M j (11Ф j' ]12+V 2). (75)

III случай. Пусть а < в, тогда применяя оценку (68) к решению задачи (70)-(73), получаем неравенство

||^]|2 < ^тах (ЦФ7!2 + г/2) • (76)

Из полученных оценок (74)-(76) следует сходимость решения разностной задачи (57)-(60) к решению дифференциальной задачи (1), (3), (45) так, что существует такое т0, что при т ^ то справедливы оценки:

1) в случае, когда а > в: ||yj+1 - uj+1]|2 < M (h2 + т2-max{«>e}( ;

2) в случае, когда а = в: ||yj+1 — uj+1]|^(01) ^ M (h2 + т2) ;

3) в случае, когда а < ß: \\yl+1 - 4+1]|о ^ М № + т2~т^а'^) , где M — const > 0, не зависящая от h и т.

Замечание. Для проведения численных экспериментов для задачи (1)-(3) необходимо привести схему (22)-(24) к расчетному виду. Как видно, слагаемое Y1Я в (22) приводит к нарушению трехдиагональной структуры матрицы коэффициентов разностной схемы (22)-(24), решение данной системы возможно с помощью метода Гаусса или какого-либо итерационного метода. Однако, если брать значения искомой сеточной функции при вычислении указанной суммы с нижнего слоя, тогда схема (22)-(24) приводится к системе разностных уравнений с трехдиагональной матрицей коэффициентов, решение данной системы возможно уже с помощью метода прогонки [28, с. 40].

6. Результаты численного эксперимента для задачи (1)—(3)

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

Ниже в таблице при различных значениях параметров а = 0.01; 0.5; 0.99, в = 0.01; 0.5; 0.99 и уменьшении размера сетки приведены максимальное значение погрешности (г = у - и) и порядок сходимости (ПС) в норме |[-]|0, когда Л = т. Погрешность уменьшается в соответствии с порядком аппроксимации. Порядок сходимости бу-

где ¿1 и ¿2 — погрешности, соответствующие

дем определять по формуле ПС = log2 шагам 0, 5h, h.

Таблица 1

а 13 К тах \\г3] о 0 <]<т ПС в |[.]0

0.01 0.01 1/20 0.000714141

1/40 0.000178643 1.9991

1/80 0.000044667 1.9998

1/160 0.000011167 2.0000

0.5 1/20 0.000835511

1/40 0.000209953 1.9926

1/80 0.000052735 1.9932

1/160 0.000013260 1.9917

0.99 1/20 0.000821499

1/40 0.000206122 1.9948

1/80 0.000051690 1.9955

1/160 0.000012983 1.9933

0.01 0.5 1/20 0.001031277

1/40 0.000287001 1.8453

1/80 0.000082060 1.8063

1/160 0.000024206 1.7613

0.5 1/20 0.000746853

1/40 0.000186587 2.0010

1/80 0.000046589 2.0018

1/160 0.000011634 2.0016

0.99 1/20 0.001016075

1/40 0.000283307 1.8426

1/80 0.000081198 1.8028

1/160 0.000024024 1.7570

0.01 0.99 1/20 0.000810871

1/40 0.000211283 1.9403

1/80 0.000057013 1.8898

1/160 0.000016374 1.7999

0.5 1/20 0.000809828

1/40 0.000211435 1.9374

1/80 0.000057191 1.8864

1/160 0.000016464 1.7964

0.99 1/20 0.000766004

1/40 0.000191709 1.9984

1/80 0.000047931 1.9999

1/160 0.000011980 2.0003

7. Заключение

В работе рассмотрены начально-краевые задачи для одномерного по пространству уравнения влагопереноса с краевыми условиями первого и третьего родов с двумя операторами дробного дифференцирования разных порядков а и в и нелокальным линейным источником. Исследование начально-краевых задач проводится с помощью метода энергетических неравенств. При предположении существования регулярного решения для каждой из рассмотренных первой и третьей начально-краевых задач получены априорные оценки в дифференциальной форме при различных соотношениях а и в, откуда следуют единственность и непрерывная зависимость решения от входных данных исходной задачи. На равномерной сетке в соответствие каждой дифференциальной задаче ставится соответствующая разностная схема порядка аппроксимации О (Л2 + т2) при а = в и О(Л2 + т2-тах{а>в}) при а = в, где 0 < а, в < 1 — порядки дробных производных Герасимова — Капуто. Для каждой разностной задачи в предположении существования решения получена оценка в разностной форме, из чего следуют единственность и

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

Литература

1. Чудновский А. Ф. Теплофизика почв.—М.: Наука, 1976.—353 с.

2. Hallaire M. Le potentiel efficace de l'eau dans le sol en regime de dessechement // L'Eau et la Production Vegetale. Paris: Institut National de la Recherche Agronomique.—1964.—Vol. 9.—P. 27-62.

3. Баренблатт Г. И. Движение жидкостей и газов в природных пластах.—М.: Недра, 1984.—447 c.

4. Шхануков М. Х. О некоторых краевых задачах третьего порядка, возникающих при моделировании фильтрации жидкости в пористых средах // Дифференц. уравния.—1982.—Т. 18, № 4.— С. 689-699.

5. Colton D. L. On the analytic theory of pseudoparabolic equations // Quart. J. Math.—1972.—Vol. 23.— P. 179-192.

6. Дзекцер Е. С. Уравнения движения подземных вод со свободной поверхностью в многослойных средах // Докл. АН СССР.—1975.—Т. 220, № 3.—С. 540-543.

7. Chen P. J., Gurtin M. E. On a theory of heat conduction involving two temperatures // Journal of Applied Mathematics and Physics (ZAMP).—1968.—Vol. 19, № 4.—P. 614-627. DOI: 10.1007/BF01594969.

8. Ting T. W. Certain non-steady flows of second-order fluids // Arch. Ration. Mech. Anal.—1963.— Vol. 14.—P. 1-26.

9. Нахушев А. М. Дробное исчисление и его применение.—М.: Физматлит, 2003.—272 с.

10. Учайкин В. В. Метод дробных производных.—Ульяновск: Изд-во «Артишок», 2008.—512 с.

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

12. Podlubny I. Fractional Differential Equations.—San Diego: Acad. Press, 1999.—368 p.

13. Беданокова С. Ю. Уравнение движения почвенной влаги и математическая модель влагосодержа-ния почвенного слоя, основанная на уравнении Аллера // Вестн. Адыг. гос. ун-та. Сер. 4: Естеств.-мат. и техн. науки.—2007.—Т. 4.—С. 68-71.

14. Головизнин В. М., Киселев В. П., Короткий И. А. Численные методы решения уравнения дробной диффузии с дробной производной по времени в одномерном случае.—М.: ИБРАЭ РАН, 2003.

15. Лафишева М. М. Разностные методы решения краевых задач для дифференциальных уравнений дробного порядка // Журн. вычисл. матем. и матем. физ.—2008.—Т. 48, № 10.—С. 1878-1887.

16. Diethelm K., Walz G. Numerical solution of fractional order differential equations by extrapolation // Numerical Algorithms.—1997.—Vol. 16.—P. 231-253. DOI: 10.1023/a:1019147432240.

17. Алиханов А. А. Априорные оценки решений краевых задач для уравнений дробного порядка // Дифференц. уравнения.—2010.—Т. 46, № 5.—C. 658-664.

18. Alikhanov A. A. A new difference scheme for the time fractional diffusion equation // J. Comput. Phys.—2015.—Vol. 280.—P. 424-438. DOI: 10.1016/j.jcp.2014.09.031.

19. Бештоков М. Х. Краевые задачи для уравнения соболевского типа дробного порядка c эффектом памяти // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки.—2022.—Т. 26, № 4.—С. 607-629.

20. Бештоков М. Х. Краевые задачи для вырождающихся и невырождающихся уравнений Соболевского типа с нелокальным источником в дифференциальной и разностной трактовках // Диффе-ренц. уравнения.—2018.—Т. 54, № 2.—С. 249-266. DOI: 10.1134/S0374064118020115.

21. Beshtokov M. Kh. The third boundary value problem for loaded differential Sobolev type equation and grid methods of their numerical implementation // IOP Conf. Ser.: Mater. Sci. Eng.—2016.—Vol. 158, № 1.—P. 1-6. DOI: 10.1088/1757-899X/158/1 /012019.

22. Бештоков М. Х., Водахова В. А. Нелокальные краевые задачи для уравнения конвекции-диффузии дробного порядка // Вестник Удмуртского университета. Математика. Механика. Компьютерные наукифиз.—2019.—Т. 29, № 4.—С. 459-482. DOI: 10.20537/vm190401.

23. Бештоков М. Х. Численное исследование начально-краевых задач для уравнения соболевжого типа с дробной по времени производной // Журн. вычисл. матем. и матем. физ.—2019.—Т. 59, № 2.—С. 185-202. DOI: 10.1134/S0044466919020054.

24. Бештоков М. Х. Краевые задачи для нагруженного модифицированного уравнения влагопереноса дробного по времени порядка с оператором Бесселя и разностные методы их решения // Вестн. Удмурт. ун-та. Матем. Мех. Компьют. науки.—2020.—Т. 30, № 2.—С. 158-175.

25. Caputo M. Linear models of dissipation whose Q is almost frequency independent-II // Geophysical Journal International.—1967.—Vol. 13, № 5.—P. 529-539. DOI: 10.1111/j.1365-246X.1967.tb02303.x.

26. Герасимов А. Н. Обобщение линейных законов деформации и их приложение к задачам внутреннего трения // Прикл. матем. и механика.—1948.—Т. 12, № 3.—С. 251-260.

27. Нахушев А. М. Нагруженные уравнения и их применение.—М.: Наука, 2012.—231 с.

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

Статья поступила 28 марта 2023 г. бештоков мурат хамидбиевич

Институт прикладной математики и автоматизации КБНЦ РАН, ведущий научный сотрудник отдела вычислительных методов, РОССИЯ, 360004, Нальчик, ул. Шортанова, 89 a E-mail: beshtokov-murat@yandex.ru https://orcid.org/0000-0003-2968-9211

Vladikavkaz Mathematical Journal 2024, Volume 26, Issue 3, P. 5-23

INITIAL-BOUNDARY PROBLEMS FOR THE MOISTURE TRANSFER EQUATION WITH FRACTIONAL DERIVATIVES OF DIFFERENT ORDERS AND A NON-LOCAL LINEAR SOURCE

Beshtokov, M. Kh.1

1 Institute of Applied Mathematics and Automation KBSC RAS, 89 a Shortanova St., Nalchik 360004, Russia E-mail: beshtokov-murat@yandex.ru

Abstract. The work is devoted to initial-boundary value problems for the fractional order moisture transfer equation with a non-local linear source and variable coefficients. Under the assumption of the existence of a regular solution for each of the considered first and third boundary value problems, an a priori estimate in differential form is obtained, which implies the uniqueness and continuous dependence of the solution on the input data of the original problem. With each differential problem associate a difference scheme on a uniform grid. Assuming the existence of a solution for each difference problem, an a priori estimate in a difference form is obtained, which implies the uniqueness and stability of the solution to the difference problem with respect to the right-hand side and initial data. Due to the linearity of the initial-boundary value problems under consideration, the obtained estimates in difference form allow us to state the convergence of the solution of each difference problem to the solution of the original differential problem (assuming the existence of the latter in the class of sufficiently smooth functions) at a rate equal to the order of the approximation error. Numerical calculations are carried out, illustrating the theoretical results obtained in the work.

Keywords: moisture transfer equation, boundary value problems, Gerasimov-Caputo fractional derivative, non-local source, difference schemes, stability and convergence of difference schemes.

AMS Subject Classification: 65N06, 65N12, 65R20.

For citation: Beshtokov, M. Kh. Initial-Boundary Problems for the Moisture Transfer Equation with Fractional Derivatives of Different Orders and a Non-Local Linear Source, Vladikavkaz Math. J., 2024, vol. 26, no. 3, pp. 5-23 (in Russian). DOI: 10.46698/l0699-2536-6844-a.

References

1. Chudnovskiy, A. F. Teplofizika pochv [Thermal Physics of Soils], Moscow, Nauka, 1976, 353 p. (in Russian).

2. Hallaire, M. Le Potentiel Efficace de L'Eau Dans le Sol en Regime de Dessechement, L'Eau et la Production Vegetale. Paris: Institut National de la Recherche Agronomique, 1964, vol. 9, pp. 27-62.

3. Barenblatt, G. I. Dvizheniye zhidkostey i gazov v prirodnykh plastakh [Movement of Liquids and Gases in Natural Reservoirs], Moscow, Nedra, 1984, 447 p. (in Russian).

4. Shkhanukov, M. Kh. Some Boundary Value Problems for a Third-Order Equation That Arise in the Modeling of the Filtration if a Fluid in Porous Media, Differentsial'nye uravneniya [Differential Equations], 1982, vol. 18, no 4, pp. 689-700 (in Russian).

5. Colton, D. L On the Analytic Theory of Pseudoparabolic Equations, The Quarterly Journal of Mathematics, 1972, vol. 23, pp. 179-192.

6. Dzektser, E. S. Equations of Motion of Groundwater with a Free Surface in Multilayer Media, Doklady Akademii Nauk SSSR, 1975, vol. 220, no. 3, pp. 540-543. (in Russian).

7. Chen, P. J. and Gurtin, M. E. On a Theory of Heat Conduction Involving Two Temperatures, Journal of Applied Mathematics and Physics (ZAMP), 1968, vol. 19, no. 4, pp. 614-627. DOI: 10.1007/BF01594969.

8. Ting, T. W. Certain Non-Steady Flows of Second-Order Fluids, Archive for Rational Mechanics and Analysis, 1963, vol. 14, pp. 1-26.

9. Nakhushev, A. M. Drobnoe ischislenie i ego primenenie [Fractional Calculus and its Application], Moscow, Fizmatlit, 2003, 272 p. (in Russian).

10. Uchaykin, V. V. Metod drobnykh proizvodnykh [Method of Fractional Derivatives], Ul'yanovsk, Izdatel'stvo Artishok, 2008, 512 p. (in Russian).

11. Samko, S. G., Kilbas, A. A., and Marichev, O. I. Integraly i proizvodnye drobnogo poryadka i nekotorye ikh prilozheniya [Fractional Integrals and Derivatives and Some of Their Applications], Minsk, Nauka i Tekhnika, 1987, 688 p. (in Russian).

12. Podlubny, I. Fractional Differential Equations, San Diego, Academic Press, 1999, 368 p.

13. Bedanokova, S. Yu. The Equation of Soil Moisture Movement and the Mathematical Model of Moisture Content of the Soil Layer, Based on the Hallaire's Equation, Vestnik Adygeyskogo Gosudarstvennogo Universiteta. Seriya 4: Yestestvenno-Matematicheskiye i Tekhnicheskiye Nauki, 2007, vol. 4, pp. 68-71. (in Russian).

14. Goloviznin, V. M., Kiselev, V. P. and Korotkij, I. A. Chislennye metody resheniya uravneniya drobnoj diffuzii s drobnoj proizvodnoj po vremeni v odnomernom sluchae [Computational Methods for One-Dimensional Fractional Diffusion Equations], Moscow, Nuclear Safety Institute of the Russian Academy of Sciences, 2003, 35 p. (in Russian).

15. Lafisheva, M. M. Difference Methods for Solving Boundary Value Problems for Differential Equations of Fractional Order, Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki, 2008, vol. 48, no. 10, pp. 1878-1887 (in Russian).

16. Diethelm, K. and Walz, G. Numerical Solution of Fractional Order Differential Equations by Extrapolation, Numerical Algorithms, 1997, vol. 16, pp. 231-253. DOI: 10.1023/a:1019147432240.

17. Alikhanov, A. A. A Priori Estimates for Solutions of Boundary Value Problems for Fractional-Order Equations, Differential Equations, 2010, vol. 46, no. 5, pp. 660-666. DOI: 10.1134/S00166110050058.

18. Alikhanov, A. A. A New Difference Scheme for the Time Fractional Diffusion Equation, Journal of Computational Physics, 2015, vol. 280, pp. 424-438. DOI: 10.1016/j.jcp.2014.09.031.

19. Beshtokov, M. Kh. Boundary Value Problems for a Sobolev Type Equation of Fractional Order with Memory Effect, Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences, 2022, vol. 26, no. 4, pp. 607-629. DOI: 10.14498/vsgtu1942.

20. Beshtokov, M. Kh. Boundary Value Problems for Degenerating and Non-Degenerating Sobolev-Type Equations with a Nonlocal Source in Differential and Difference Forms, Differential Equations, 2018, vol. 54, no. 2, pp. 250-267. DOI: 10.1134/S0012266118020118.

21. Beshtokov, M. Kh. The Third Boundary Value Problem for Loaded Differential Sobolev Type Equation and Grid Methods of Their Numerical Implementation, IOP Conference Series: Materials Science and Engineering, 2016, vol. 58, no. 1, pp. 1-6. DOI: 10.1088/1757-899X/158/1/012019.

22. Beshtokov, M. Kh. and Vogahova, V. A. Nonlocal Boundary Value Problems for the Fractional Order Convection-Diffusion Equation, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2019, vol. 29, no. 4, pp. 459-482. DOI: 10.20537/vm190401.

23. Beshtokov, M. Kh. Numerical Analysis of Initial-Boundary Value Problem for a Sobolev-Type Equation with a Fractional-Order Time Derivative, Computational Mathematics and Mathematical Physics, 2019, vol. 59, no. 2, pp. 175-192. DOI:10.1134/S0965542519020052.

24. Beshtokov, M. Kh. Boundary Value Problems for a Loaded Modified Time-Fractional Moisture Transfer Equation with a Bessel Operator and Difference Methods for Their Solution, Vestnik Udmurtskogo universiteta. Matematika. Mekhanika. Komp'yuternyye nauki, 2020, vol. 30, no. 2, pp. 158-175 (in Russian).

25. Caputo, M. Linear Models of Dissipation whose Q is almost Frequency Independent-II, Geophysical Journal International, 1967, vol. 13, no. 5, pp. 529-539. DOI: 10.1111/j.1365-246X.1967.tb02303.x.

26. Gerasimov A. N. Generalization of Linear Deformation Laws and their Application to Problems of Internal Friction, Prikladnaya Matematika i Mekhanika, 1948, vol. 12, pp. 251-260 (in Russian).

27. Nakhushev, A. M. Nagruzhennyye uravneniya i ikh primeneniye [Loaded Equations and their Application], Moscow, Nauka, 2012, 231 p. (in Russian).

28. Samarskii, A. A. Teoriya raznostnykh skhem [Theory of Difference Schemes], Moscow, Nauka, 1983, 616 p. (in Russian).

Received March 28, 2023 Murat Kh. Beshtokov

Institute of Applied Mathematics and Automation KBSC RAS, 89 a Shortanova Str., Nalchik 360000, Russia, Leading Researcher Computational Methods Department E-mail: beshtokov-murat@yandex.ru https://orcid.org/0000-0003-2968-9211