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

Владикавказский математический журнал 2023, Том 25, Выпуск 3, С. 15-35

УДК 519.63

DOI 10.46698/c8748-9711-0633-d

ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ НЕЛОКАЛЬНЫХ КРАЕВЫХ ЗАДАЧ ДЛЯ ОБОБЩЕННЫХ НАГРУЖЕННЫХ УРАВНЕНИЙ АЛЛЕРА

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

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

Аннотация. Работа посвящена нелокальным краевым задачам для одномерных по пространству нагруженных уравнений Аллера с переменными коэффициентами и двумя операторами дробного дифференцирования Капуто с порядками а и в. Подобные задачи возникают в практике регулирования солевого режима почв, когда расслоение верхнего слоя достигается сливом слоя воды с поверхности затопленного на некоторое время участка. Для численного решения поставленных задач на равномерной сетке построены разностные схемы. Методом энергетических неравенств при различных соотношениях между порядками дробной производной Капуто а и в получены априорные оценки в дифференциальной и разностной трактовках для решений нелокальных краевых задач. Из полученных априорных оценок следуют единственность и устойчивость решения по правой части и начальным данным, а также сходимость решения разностной задачи к решению соответствующей исходной дифференциальной задачи (в предположении существования решения дифференциальной задачи в классе достаточно гладких функций) со скоростью O(h2 +т2) при а = в и O(h2 +т2-max{a,e>) при а = в. В работе также приводится алгоритм численного решения нелокальной краевой задачи для нагруженного уравнения Аллера с переменными коэффициентами и оператором Бесселя.

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

AMS Subject Classification: 65N06, 65N12.

Образец цитирования: Бештоков М. Х. Численные методы решения нелокальных краевых задач для обобщенных нагруженных уравнений Аллера // Владикавк. мат. журн.—2022.—Т. 25, вып. 3.— C. 15-35. DOI: 10.46698/c8748-9711-0633-d.

1. Введение

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

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

порядка. Одним из приложений этой теории является теория дифференциальных уравнений с дробными производными. Подробное описание применения дробного исчисления к различным областям науки и техники на современном этапе даны в монографиях [1—4].

Вопросы моделирования фильтрации и течения жидкости в трещиновато-пористых средах [5, 6], двухфазного течения в пористых средах с динамическим капиллярным давлением [7], переноса влаги [8, 9], движения подземных вод со свободной поверхностью в многослойных средах [10, 11], теплопроводности в двухтемпературных системах [12] и течения некоторых неньютоновских жидкостей [13] связаны с необходимостью исследования краевых задач для дифференциального уравнения в частных производных третьего порядка псевдопараболического типа (уравнение Аллера) и более общего класса уравнений — уравнение типа Соболева [14].

В работе [15] предложены и исследованы математические модели водного режима в почвогрунтах с фрактальной структурой. В основе этих моделей лежат дифференциальные уравнения Аллера с дробной по времени производной. Для задачи Коши выписывается единственное представление решения.

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

В данной работе рассматриваются нелокальные краевые задачи для нагруженных обобщенных уравнений Аллера с переменными коэффициентами и дробной производной Капуто. Основной результат работы заключается в доказательстве априорных оценок для решений задач как в дифференциальном, так и в разностном виде при различных соотношениях между порядками дробной производной Капуто (а, в). Полученные неравенства означают устойчивость решения относительно входных данных. В силу линейности рассматриваемых задач эти неравенства позволяют утверждать сходимость приближенного решения к точному (в предположении существования последнего в классе достаточно гладких функций). Подобные задачи возникают также в практике регулирования солевого режима почв, когда расслоение верхнего слоя достигается сливом слоя воды с поверхности затопленного на некоторое время участка [21, с. 233]. В работе [22] для описания движения примеси в потоке однородной жидкости используется дифференциальное уравнение дробного порядка. Если на поверхности поля имеется слой воды постоянной толщины Л, то на верхней границе с учетом фрактальности почвенной среды следует задать условие

где c — концентрация соли в почвенном растворе, D — коэффициент диффузивности, A, h = const > 0.

Настоящая работа является непосредственным продолжением серии работ автора, посвященных разностным методам решения нелокальных краевых задач для обобщенных уравнений Аллера [23-27].

В замкнутом прямоугольнике СЦТ = : 0 ^ ж ^ 0 ^ £ ^ Т} рассмотрим

нелокальную краевую задачу для нагруженного обобщенного уравнения Аллера:

D — + Ад ■ = hd&c,

2. Постановка нелокальной краевой задачи для нагруженного уравнения Аллера

t) + Зм ^ (г](х) + r(x, q(x, t)u(x0,t) + f(x, t),

(1)

0 <x <l, 0 < t ^ T,

n(0,t) = eii(t)u(0,t) + ei2daiu(0,i) - ^i(i), 0 < t < T, (2)

-n(1,t)= P2l(t)u(l,t) + P22dgtu(l,t) - /J-2 (t), 0 < t < T, (3)

u(x, 0) = u0(x), 0 ^ x ^ l, (4)

где

0 < c0 ^ k(x,t),n(x) ^ c1, в12,в22 = const ^ 0, |eil(t)|, |^2l(t)|, |r(x,t)|, |q(x, t)|, |kx(x,t)|, |rx(x,t)| < C2, (5)

= r(i-7) Jo U(t-T)-y ~ дробная производная в смысле Капуто порядка 7,

0 < y < 1, a ^ в, x0 — произвольное фиксированное число, n(x,t) = kux + d0t(n(x)ux), Cj = const > 0, i = 0,1, 2.

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

По ходу изложения будем также использовать положительные постоянные числа Mj,

1 = 1, 2,... , зависящие только от входных данных рассматриваемой задачи.

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

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

(<9otu, u) = ((kux)x, u) + (dot(nux)x, u) + (rux, u) - (qu(xo, t), u) + (f, u), (6)

где скалярное произведение и норма имеют вид (a, b) = J0 abdx, (а, а) = ||а||0, а, b — заданные на [0, l] функции.

Пользуясь неравенством Коши с е, леммой 1 из [19], после несложных преобразований [27] с учетом

i i — (qu(xo, t), u) = — J qu(xo, t)udx ^ ^ J {u2(xo, t) + u2) dx 00

i

^ Щи2(х0,1) + Ц ju2 dx ^ Mi\\u\\2q +e1-^ \\ux\\20

из (6) находим

1 1

о о

Оценим первое слагаемое в правой части (7). Тогда получим

+ и(0, *) (^1 (*) - вии(0, *) - в12д0Х0, *)) < М4 (^2 + ^2) (8)

+ М51Н12 + еМ6\\их 11^ - 0,4)

Учитывая (8), из (7) находим I

да НЮ + У п4ы2 dx + ||их 110 + в12до>2(0,£)+ ^бОи2^, *)

о

< М7||и||2 + еМ8||их||2 + Мд (||/1|2 + ^2(*) + , (9)

Выбирая е = ^Щ) из (9) находим

I

до" ||и|2 + /(их)2 dx + ||их||2 + д^и2(0, ¿) + о

< М1о||и||2 + ми (||/1|2 + + ^2(*)), (10)

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

||и||2 + д0-4(а"в)|их|2 + я0Г||их||2 + и2(0,*)+и2(М)

< М12^0"а|и|2 + М^"^||/1|2 + + + |Ь(х)||2). (11)

На основании [19, лемма 2] из (11) находим априорные оценки:

1) в случае, когда а > в

||и||2 + в)|их|2 + Я°Г||их||2 < ||/1|2 + + ^2(*)) + ||ио(х)|^, (12)

2) в случае, когда а = в

где ||u||Wi(o г) = IMIci + ll«x|lo, M = const > 0, зависящее только от входных данных (1)—(4), Dq^u = р^у Jq — дробный интеграл Римана — Лиувилля порядка 7, 0 < y < 1. 7

Теорема 1. Если k(x,t) € Cl>°(QT), ф) € С1 [ОД r(x,t),q(x,t), f(x,t) € C(QT), u(x,t) € C2|1(Qr) П C1'1^), d&ufat) <E C(QT), d^uxx{x,t) € C(QT) и выполнены условия (5), тогда для решения задачи (1)—(4) справедливы априорные оценки (12) при а > в и (13) при а = в.

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

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

Для решения задачи (1)—(4) применим метод конечных разностей. Построим монотонную схему второго порядка точности, содержащую односторонние производные, учитывающие знак r(x, t). Для этого рассмотрим вместо уравнения (1) следующее уравнение с возмущенными коэффициентами [28, с. 170]:

dot« = Kj(fcux)x + d0t(nux)x + rux - q(x, t)«(xo, t) + f (ж, t), (14)

где x = ^рд, -R = _ разностное число Рейнольдса.

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

Л^,У = а*уХ^)х + Л^ (7^Ух)х + 6-

+ 6+(аг+1»х5 - dj (У(0Г)х- + У(0+1Х+) + &, (15)

коа1»х°о + ^^ (71Ух,о) = в11»оа) + (уг(0ст)х" + У(0+1Х+) + /512ла^^.+стУо - Дъ (16)

- амуХ°" + ^^ (7мУх,м)) = в21У(^) + 0.5Мм (у^Чо + У(0Г+1Х+0)

+ в022Ла4.+„Ум - Д2, (17)

у(х, 0) = По(х), X € Ш/г, (18)

где

г(0, ¿) < 0, г(М) ^ 0, (19)

отметим, что условия (19) нужны для обеспечения порядка аппроксимации О (Л2 + т2) разностной схемы (15)-(18),

в12 = в12 + 0.5Л, Д^+а) = Д1(() + 0.5Л^о, в22 = в22 + 0.5Л, ) = Д2 ) + 0.5Л(,

= к(х^0.5,^+а), Ъ = Г](%г-О.б), ^ = , = ДхгЛ'+^)

УИ = ^ + (1 - , = д(х*, ^), ао7'ст) = а1 -7, г(хЛ'+а) = г+ОМ^ + г-(хД), |г(х,-()| = г+(х,(+ст) - г-(х,(+Д

а|7'ст) = (I + а)1-7 - (I - 1 + а)1-7, I > 1,

ь^ = + ^)2"7 - а -1 + ^)2"7] + + (1-1 + , г > 1,

при 3 =0, Со7'ст) = ао7'ст);

'ао7^ + ЬМ 8 = 0,

при 3 > 0, = ^ а!7'ст) + б^ - б!7'ст), 1 < 5 < 3 - 1,

- б(7'ст), 5 = 3,

(т.ст) ^ 1 - I \-т ^ п - х*0+1 - хо + хо - х*о ^ ^ 47' ' > —2~(в + <г) 7 > 0, х^0 = ---, х+ = ---, хзд < х0 < Хг0+Ъ

о- = 1 - при а = fi и а = 0.5, при а ф fi, Л^.+<тУ = Y,l=о tf-s'yi - дискретный

аналог дробной производной Капуто порядка 7, 0 < 7 < 1, [20] с порядком аппроксимации 0(т3-7) при а = 1 — и 0(т2-7) при а = 0.5. Введем скалярные произведения и норму:

. v^ Î0.5h, i = 0,N, , п ^ п г 2П

[u,v]=2_^ UiVift, п = < ^ i = On = Uîvih, [u,u| = [1,u J =

u = 11, u2! = Ы |0.

i=0

i=1

Перепишем (15)-(18) в операторной форме:

у(х,0) = и0(х), XEWh,

(20) (21)

где

Ayf) = к(ayf^ + rayf + ft+a^y^ - d^x" + y^),

А+У<^] = \{~ ЯМаМу(А ~ fayp - °-bhdN (vï oCT4~o + violât)) '

i = 0, i = N,

= Abij+CT (7i Ух)х

ôy=\5 Уо = Л (A0ij+(J (712^,0) - /5i2A^.+(jyo),

i = 1,N - 1,

i = 0,

S+yN = ^ ( - AotJ+CT (4nVx,n) - P22^mi+aVN), i = N,

«j + CT * /

^ = , 2

i = 1,N - 1,

Ф = = ^ (jLtiCij+o-) + 0.5/г,^), г = 0,

= ^ Mij+a) + 0.5Vw), г = N.

к =

1 +

0.5%| k 1

к0

KN

0.5%p| fco.6 1

1+

0.5%дг| fcjV-0.6

Г0 < 0,

rN ^ 0,

ro = r(0, tj +a ) = ^ 0, rN = r(XN ,tj+a ) = Г^ ^ 0.

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

д:

0tj+CT y, y

и

(22)

Пользуясь леммой 1 из [20], первой разностной формулой Грина, неравенством Коши — Буняковского и е-неравенством [28, с. 110], после некоторых преобразований с учетом

+\ id,y(ст)]2 < (ЙЧ)2 + (yi:li<)2 + ^ [i. (уИ)2] < ЩЫШ + ^ЦУГС

из (22) находим

дигу]|2+<+„W2+1[уХст)]|0 < МЗ|ГУ(-)]|2+]|2+м5(|Г^]|2+^2+^2), (23) где |[y]|1 = |[y]|2 + У2 + yN.

1

I случай. Пусть a > /3, тогда выбирая е = jщ, из (23) получаем

Agj J[y]|J + Aj |Ы1о + ШЮ < + M7(|[^]|0 + д? + Д?). (24)

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

AjJMIi < M8-|[yj+1]|; + M9°|[yj]|1 + M7(|M|0 + Д? + Д?). (25)

На основании [25, лемма 4] из (25) получаем

|[yj+1]|0 < M(|[y0]|? + max (|[^']|? + Д? + Д?)), (26)

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

II случай. Пусть a = в, тогда перепишем (23) в следующем виде:

Aotj.jM|? < Mfo|[yj+1]|? + Mfi|[yj]|? + M7(|[^]|0 + + Д?), (27)

где 1 [y] 10 = |[y]|W2i(0,1) + у0 + yN, |[y]|W2i(0>1) = |[y]|0 + Ы0.

На основании леммы 4 [25] из (27) получаем

|[yj+1]|^(0Л < 4|[y0]|W2i(0,i) + j ']|0 + + Д?))' (28)

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

Теорема 2. Пусть выполнены условия (5), (19), тогда существует такое Т0, что если т ^ т0, то для решения разностной задачи (15)—(18) справедливы априорные оценки (26) при a > в и (28) при a = вИз априорных оценок (26), (28) следуют единственность и устойчивость решения задачи (15)—(18) по начальным данным и правой части.

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

Aj y = x^y^) ^+A0tj+CT (7iy2 )xaiyg a^ -jy^ x" +y(0+ix+), (29) K0aiyX°0 + Aj (7iyx,0) = в^0 + 0.5hd^y(Qo) x" + y(0+iX+) + вю Aj y0, (30)

- aNy^i + Aotj+„ (YNyx,Nd = в21 y(V) + 0-5hdN (y(0O)X-~ + уЙ-1х»+) + A?A0tj+„yw, (31)

y(x, 0) = 0, (32)

Пусть y(x,t) — одно из решений однородной задачи (29)-(32). Из неравенства (26) следует, что |[yj+1]|0 < 0, а из (28) — |[yj+1]|^(0Л < 0, но |[yj+1]|0 = 0 и |[yj+1]|Wi(0,0 = 0 лишь при y = 0, x € Поэтому из априорных оценок (26) при a > в, (28) при a = в следуют, что единственным решением однородной задачи (29)-(32) является y = 0. Тем самым решение задачи (15)-(18) при любых /¿1, Д, u0(x) существует и единственно.

Таким образом, из априорных оценок (26) при a > в, (28) при a = в следуют существование, единственность и устойчивость решения задачи (15)-(18) по начальным данным и правой части.

Пусть и(ж,£) — решение задачи (1)-(4), у(ж^,^^) = У — решение разностной задачи (15)-(18). Для оценки точности разностной схемы (15)-(18) рассмотрим разность = У* — и^, где и^ = и(жг,^). Тогда, подставляя у = г + и в соотношения (15)-(18),

получаем задачу для функции z: Aj z = К (x + Aj friz*), + « z$

+ a+izgi - dj (z£>*r + z^x*) + Ф?, (33)

KoaizX°o + Aj(Yizx>o) = в^ + 0.5hdo (z^x" + 4+i<) + ^Ajzo - Pi, (34)

- aw+ Aetj+CT (YNzi,N^ = e2izi? + 0-5hdN (4*4- + 40+iX;+)

+ e022Agt.+ff zw - P2, (35) z(x, 0) = 0, xGUJh, (36)

где Ф = O(h2 + т2), Pi = O(h2 + т2), P2 = O(h2 + т2), при a = в, и Ф = O(h2 + т2-max{«>e}), pi = O(h2+т2-max{«>e}), p2 = o(h2+т2—max{a>e}), при a = в - погрешности аппроксимации дифференциальной задачи (1)-(4) разностной схемой (15)-(18) в классе решении u = u(®,t) задачи (1)-(4).

В случае, когда a > в, применяя априорную оценку (26) к решению задачи (33)-(36), получаем неравенство

|[zj+i]|o < M max. (|[Ф'"]Ю + vj'2 + vf), (37)

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

В случае, когда a = в, применяя априорную оценку (28) к решению задачи (33)-(36), получаем неравенство

|[zj+i]lW(o,i) < Mj (|[ф?"]io + vj'2 + v2^), (38)

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

Из априорных оценок (37), (38) следует сходимость решения разностной задачи (15)-(18) к решению дифференциальной задачи (1)-(4) так, что существует такое to, что при т ^ To справедливы оценки:

1) в случае, когда a > в

|[yj + i - uj + i]|o < M(h2 + т2-тах{а>в}) ;

2) в случае, когда a = в

|[yj+i - Uo,) < M(h2 + т2).

Замечание 1. Полученные в данной работе результаты справедливы и в случае, когда уравнение (1) имеет вид

д ( дди \ д ( дди \ дди Р

\ / \ / 5=1

0 < ж < I, 0 < £ < Т,

если потребовать выполнения условия |д8| ^ С2.

Замечание 2. При построении алгоритма приближенного решения дифференциальной задачи (1)—(4) следует использовать метод параметрической прогонки [29, с. 131]. В виду того, что уравнение (1) содержит слагаемое д(ж,£)«(жо , ¿) нарушается трехдиаго-нальная структура матрицы коэффициентов разностной схемы (15)—(18) и использовать обычный метод прогонки не представляется возможным.

5. Постановка нелокальной краевой задачи для нагруженного уравнения Аллера с оператором Бесселя

В замкнутом прямоугольнике СЦТ = {(х,£) : 0 ^х ^ 1,0 Т} рассмотрим

следующую нелокальную краевую задачу:

90> = — xmk(x, — 9gt (хтц{х)

0t xm dx \ v ' dx J xm 0t dx \ ' dx J

du

+ r(x, q(x> t)u(x0,t) + /(x, t), 0 < x < I, 0 < t ^ T,

(39)

lim xmn(x, t) = 0, 0 < t < T, (40)

-П(1, t) = в1 (t)u(l, t) + ÄÖgXZ, t) - д (t), 0 < t < T, (41)

u(x, 0) = uo(x), 0 < x < l, (42)

где 0 ^ m ^ 2, в2 = const ^ 0.

При x = 0 ставится условие ограниченности решения |u(0, t)| < то, которое эквивалентно условию (40), равносильному в свою очередь тождеству n(0,t) = 0 [28, c. 173], если функции r(0,t), k(0), q(0,t), f (0,t) конечны.

Уравнение Аллера с оператором Бесселя (39) возникает при переходе от трехмерного уравнения Аллера к цилиндрическим (сферическим) координатам в случае, когда решение u = u(r) не зависит ни от z, ни от

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

Получим априорную оценку методом энергетических неравенств, для этого умножим уравнение (39) скалярно на хти :

(до>,Л) = ((Лиж)ж,и) + (д0ижтпи4,и)

+ (гиж,хти) - (ди(жо,*),Л) + (/,хти). (43)

После некоторых несложных преобразований с учетом

1 I

- („«(^(.Л) = -/^«(аМ)«^ = -«(OM)/

oo

/ 1 \ 2

1 9/ 1 I f т \ II — 112 И HL

^ - u (xQ,t) + - I х qudxj ^ Mi||a; 2 u||0 + еМ2||ж 2

из (43) находим

7 7

■ и\Ц + i J r]dQt(x^~ux)2 dx + J xmkv2dx

2 -uhI-^IIo + ^ j Vdot{x™ux)2 dx + j xmku2x oo

, |7 ii m 112 ii m u2 In m II2

< xmuU(x,t)\0 + M3\\x~u\\0 +еМ4\\х~их\\0 + - ||ж~/||д. (44)

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

ЛП(ж, t) |o = lmu(l, t)n(l, t) = lmu(l, t) (^(t) - в^М^, t) - в2да4и(1, t))

= lmu(l,%(t) - ru2^,^(t) - ^u^dotu^t) (45)

im в

<--7^d$tu2(l,t) +M5\\xfu\\l+eM6\\x^ux\\l + M7^(t).

Учитывая (45), из (44) находим

7

/-v 11 m. 112 I 8 / m \ 2 и m. 112 9

9ot|p2u||0+ / r]dQt[x^ux) dx + \\x 2 ux\\Q+/32dfitu (l,t) o

и m 112 ll?ri 112 /и m 112 / %

< Ms\\x^u\\Q + sM9\\x^ux\\0 + Мю(||ж^/||0 +/x2(i))- (46) Выбирая e = щ, из (46) получим

7

И m ||2 / 8 / — \ 2 и m 112 rv9/

9ot|p2u||0+ / rjd^t[x^ux) dx + ||ж 2 иж||0 + c)Qtu [l,t) o

^Mn\\xfu\\20 + M12(\\xff\\20+^2(t)). (47)

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

и т. 112 , т^—(а—в)\\ ш. 112 , T-._Q.ii 22. ||2 , 2п

\\Х 2 и||0 + Doty И;\\Х 2 их\\0 + -D0i F 2 ИИ1о +и \ЧЧ

^ M13D^ta\\xfu\\20 + Mu^ta(\\xf f\\l + /лЩ + Цж^иоЦо). (48)

На основании [19, лемма 2] из (48) получаем априорные оценки:

1) в случае, когда a > в

II — II2 , г-,—(а—/3) || 22. 112 , „-„и 22. 112

||ж 2 и||0 + Doty > >\\х 2 иж||0 + Dota\\x 2 иж||0

< м(^0-"(||жт/||0 +fx2(t)) + ||Жтзд(ж)||0); (49)

2) в случае, когда a = в

11 га 112 _^м ZH. 112 / _г\ (II — 112 О / \ \ II IZ1 112 \

\\х 2 u\\wj(o,i) ^-^ot \\х 2 их\\о ^ МyDot (JIж 2 У11q + + F 2 W) 11^1(0,0J> (5°)

— цО Ш. мО IZ1 мО

где ||ж 2 г) = ||ж 2 u\\q + ||ж2иж||о, М = const > 0, зависящая только от входных

данных задачи (39)-(42).

Теорема 3. Если к{х,г) € С1,0(дг), ф) € С1 [ОД г{х, д(х, /(ж, *) € С(<5Т), •и(жД € С2-1^) п С1-1^), € С(дт), д&ихх(х,$ € С(дт) и выполнены

условия (5), тогда для решения задачи (39)-(42) справедливы априорные оценки (49) при а > в и (50) при а = вИз априорной оценки (49) следуют единственность и устойчивость решения по правой части и начальным данным.

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

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

= ^ С''"...:."^" ) ,. + ^ KW*»*,.), + ^ (хГ-о.5

+ ^Г o.sfli+iy^) " < + 2ЙА+) + ^, (ж, t) € Wfc,.

(51)

xoaiy^ + A^.+CT(7iyä;0) = (AotJ+CT2/o + do(y£40 + - ßi, (52)

а^уЦг-(yn) = Äyj? +0.5hd^y(oCT)x"o + + Ä^t^yw-Д2, (53)

у(ж, 0) = Uo(x), x GcJh, (54)

где

0.5h

/?! = xßi(tj+l7), ß2 = xß 2 + 0.5h, fil = ^ 1 Д2 = Xß(tj+a) + 0.5h,LpJN, Hö =-ггтп—i—) если rn+<J ^ 0, xjv =--—77— если ^ 0,

^ (m+l)fcS+6ff 1 + fcjv-о.Б

r = r+ + r-, |r| = r+ — r-, r+ = 0.5(r + |r|) ^ 0, r- = 0.5(r — |r|) ^ 0,

СЦ = k(Xi-0.5,tj+a), Ъ = фг-О.б), = -

j+ö"

г = l,iV — 1, ? jxj/"'. г = 1, JV — 1, 10.5h, ¿ = 0,

* [qt+a, i = 0,N, ^ [//+<J, г = 0,AT, ~ \h, i = l,N-l,

_ 1 p _ 0.5h|ri|xi

^г — ~—;—TT, -Щ —

1 + Ri ' ki-

i—0.5

mim — 1 )h2 ——- _ 0.5hm

= H----И—, г = 1, iV — 1, я = 1 +

24ж2 ' I 1 - '

Найдем априорную оценку методом энергетических неравенств, для этого перепишем (51)-(54) в операторном виде

к у = + ¿у + Ф, (55)

1

у(ж, 0) = Ио(ж),

(56)

где

= Хг, Ж € Ши, К = <

11, ж = 0,1,

К» = 1 +

ш(ш — 1)Л,2 24ж|

ж

Жт-0.57гуж,г)

¿у = <

5 УО = д^+<7(71Ух,О),

0.5Л 2

(ж, ¿) € ж = 0,

= + X = 1,

(а)

Л(^+а )уа = <

д - (а) т + 1 / (а)

Л Уо = "77777" к^Ух1

л+у(а) = -

0.5Л 2

Л

км у^ + Ау^ + °-5Мм (у(0а)жг0 + у(а+1 ж+))

ж = 0, ж = 1,

Ф

Р = Рг,

(ж, ¿) € ш^т,

- т + 1 _

= ж = 0,

+ —

0.5Л 1

0.5Л

Д2,

ж = 1.

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

(Л(^+а)у(0) ,жту(0)

+

жту(а)

+ (Ф,жту(а)], (57)

где

N N (

«¿^Я, 11 и] |2 = ^ и2Я, Я =< г=1 г=1 I

0.5Л, г = 0,#, Л, г = 0, N.

Пользуясь неравенством Коши с е, после нетрудных преобразований из (57) получаем

к

— Да

2 '

/ т. \2 Ж 2

у)'

+ А13

+ I 2 , 04.

7 +

/_2П \2

Ж 2

"2/я)'

+

ж агк,

1 + ЛМ1

< м2(и + \\TbiX) - (¿(у^^ + У^-К),^) (58)

+ - х{хмаму^м + хУх,м)) ~ ж™5У^ {*о(ИУ^о + Дад^Ст^х.о))

+ (р, жту(а)) + у^ (Д2 - в1у!°а) - 0-5Мм (у^ж-о + у^ж+О) - в2Д^ум).

X

1

Рассмотрим третье, четвертое и шестое слагаемые в правой части (58):

(ж^ - (хмамУ^ + Дад^ Ыш,м)) ~ х^5у{0а) {ж0а1У^]0 + (ъУх,о))

+ (Д2 - - 0.5Мм (у^*" + у^) - МО^ум = (ж^ - х^уМ [д2 - - 0.5М* (уг(0ст)ж" + у^Л - /32Д

УМ

+- /31 уЙ° - 0.5ММ (у(оСТ)®-о + Уй-Ао) - в2А:

уМ

т ( ®0.5у0

(59)

/?2 /_177, тД л а; 2 ^

- у - ^м) - 4(т + 1)

Учитывая преобразования (59), из (58) находим

Жа5ДМ.,-+ст у0-

+

(я™У)2 Н

4(ш + 1) ж™5А^.+стУО2 +

0

X в2

(60)

Преобразуем первое и пятое слагаемые в левой части (60) с учетом ж™_05 ^ ^ж™-

да

(ж™у)2

+

3 в2

~ жм) у I А"

04

7 +

(ум )2

гп Л2\ , 0.5Л, тл„ , 2 , I , /-т,

(хЬГ) + -м + ^ -м) Ы)2 > ^ (1, Аа

(61)

2

+ -4 (УМ)2 >± (1, у)2) (ум)2 > ± Ао£.+а. ||Ж?

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

< МиЦжТуМ]!2 + еМ12||жТ^>]|2 + М13(||жТ+ Д2 + Д2

04

( 22. \2Ч Ж 2

0.5Н

т -| 12 0"

I 22 и 12 , II 22. (<т)-| 12

Р2Уй п+ Р2Уй

'х Л0

(62)

где ||ж™у]|2 =

* Го

|Ж 2 у||„ + ( ж025у0 ) .

I случай. Пусть а > /3, тогда выбирая е = 2м12 > из получаем

жТу]|2 + + ЦжТ

< МиЦх^у^Х + М15(||ж%||0 +Д2 + Д2). (63)

2

2

2

2

2

2

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

Aotj+A\^y]\l < + Mr7||xfy]|i + ^llo + Д2 +/¿2), (64)

На основании [25, лемма 4] из (64) получаем

\x^yj+1]\l ^ M^\\xfy°]\l +0™ах (|ИИ1о + Д2 + ), (65)

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

II случай. Пусть а = в, тогда перепишем (62) в следующем виде:

< H-^fgll^^lla + +Д? + Д1). (66)

где ||жту]|2 = ||жту]|2 + + (ж^5у0) ■

Пользуясь [25, лемма 4], из (66) находим априорную оценку

\\xfyj+1]\l < м[\\х^ +Д2) (67)

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

Теорема 4. Пусть выполнены условия (5), (19), тогда существуют такие то, ho, что если т ^ т0, h ^ ho, то для решения разностной задачи (51)—(54) справедливы априорные оценки (65) при а > в и (67) при а = в.

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

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

_д а К f m А 1 л в / m Л Ь j / m

= ^та (^г-О.б^й J^, + ^да 0ij+CT \xi-0.blizx,i)x + ^

+ (*но.БОн-1*$) - di + illXZ) + > i) e

(68)

«ьа^Ц + <„„(71^,0) = ^ (а&,+„*> + ¿о(4Ч0 + ^^гъ)) " (69) -кмам40* - д& .+„(7м) = (4°ж" + ^ж+О) + & да .+„¿м -¿2 (70)

г(х,0) = 0, хешн, (71)

где Ф =( О(Л2 + т2), г>1 =) О(Л2 + г2), г>2 = О(Л2 + г2) при а = в и Ф = О(Л2 + т2 - тах{«>в}), г>1 = О(Л,2 + т2 - тах{а>в}), р2 = О(Л2 + т2 - тах{а>в}) при а = в — погрешности аппроксимации дифференциальной задачи (39)-(42) разностной схемой (51)-(54) в классе решении и = и(ж,£) задачи (39)-(42).

В силу того, что задача (68)-(71) линейна, применяя априорные оценки (65), (67) к задаче (68)-(71), получаем оценки:

1) в случае, когда а > ß

+ ¿2 + *>22), (72)

+ ¿2 + *>22), (73)

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

Из оценок (72), (73) следуют единственность и устойчивость решения разностной схемы (51)—(54) по правой части и начальным данным, а также сходимость решения разностной задачи (51)-(54) к решению дифференциальной задачи (39)-(42) со скоростью O (h2 + т2) так, что если существуют такие To,ho, то при т ^ то, h ^ ho справедливы априорные оценки:

1) в случае, когда а > ß

||ж? 1 _ UJ+1)]|1 < M\\xf II 0(h2 + r2-max{a;/3}) ^ M (h2 + r2"1""^'^);

2) в случае, когда а = ß

\\xf (У+1 -UJ'+1)]|2 < M\\xT\\0(h2 +T2) < ~M(h2 + r2),

где M = const > 0, не зависящая от ft и т,

Замечание 3. Полученные в данной работе результаты справедливы и в случае, когда уравнение (39) имеет вид

öo> = — #- i) ^ + — S6 #■ (хтф)^)

0t ж™ dx V v 7 dx J ж™ 0t dx V /v 'dxj

+ г(ж, qs(x, t)u(x8,t) + /(ж, i), 0 < ж < I, 0 < t < T,

s=1

если потребовать выполнения условия |qs| ^ С2.

11 — о _i_ 1 п 12 / 11 222. о/ 112

ж 2 L ^ M max ж 2 Ф-ЧI 11 jiI 711 110

где M = const > 0, не зависящая от h и т; 2) в случае, когда а = в

11 — i А-1 "112 /|| 222. о/ 112

ж 2 L ^ М max ж 2 Ф-7

" Jl2 o<7'<7 V" 110

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

Для численного решения задачи (39)—(42) воспользуемся методом параметрической прогонки [29, с. 131] и приведм разностную схему (51)—(54) к расчетному виду. Тогда уравнение (51) приводится к следующему виду:

Агу£1 - СгУ1+1 + Вгу£1 - 1г2тахГ4 [у^1 ж" + у^ж+) = -В*, г = МУ^Т, (74) где

т 1-вс(в>^) _. . А = таз4+ ъхТ-о.ъ - гЪхтхТ-о.ъ^ 3а\,

В = тстК xm

т 1-в. .

г+0.5аг+1 + 7г+1жг+0.5 — /5) ^ T^l<TXi+0-5ai+1>

1-а («,а)

Сг = Аг + Вг + к2щхТ г(2_°а)

т 1-а *-1 , ч

П = ААгУ1г - ССгУ1 + ВВгУ1+1 + - Н*хТ Е с<£? - у|)

ж^ т 1-в * -1

+ Е4-? ((7тУтГ+1 - ЫгУыУ) - Т(1 - а)^хТШох- +

жтт 1-в ч

Е Ш + 7тЬГ+1 - ((7, +

жт т 1-в *_1

+ г"(2-уЗ) Е^Ьй-О^1 " (ЪУг-гУ),

т1 с • •

АМ = т(1 - (т)фТ-О.Ъа\ - ЪхТ-0.5 ° - -

.1-в с(в,о) с0

Г(2 -/3)

т 1-в с(в>а) , ,

ВВг = т(1 - - Ъ+1ХТ+0.Ь + - о)х™+ о.ф^а1+1>

т 1-«с(«,а) 2— т ' с0

СС, = + ББг - /^¿Ж?

Г(2 - а) ' Краевое условие (52) принимает вид

уо = К11у1 + К12уго + К1зуго+1 + Д1, (75)

где

* , т1 в Сп

ГО-ХОЙ! +71 Г(2—/3)

тггъгпП3 А- л--, т со 4- а5/г г с0

г<тхоа1 + 71 г(2_/3) + т+1 г(2_а)

°'5 - жо)

т+1

^12 = -~

тшуГСП1 А- л--, О -I- а5/г т со

г<тхоа1 + 71 г(2_/3) + т+1 г(2_а)

0.5

Йо^Лт (жо - ж»0 )

=

_ _т+1_

„ 9 ,г1-а!,,(а!><Т)

тггтугсп1 А- л--, о 4- а5/г т со

г<тхоа1 + 71 г(2_/3) + т+1 г(2_а)

т 1-вс(в'а) • • 0 5Л2 т 1-а с(а'а) ДхЛг + г(1 - а)к0а!Ы - - Ы " + ^ТТ Т(2^Г

0 5л2 т 1-а , ч т 1-в *-1 ,,

- ^ ет ' - + щгд Е«^' «™>'+1 - <™>')

т

1-в *-1 0 5

Е ((7ИА>)в+1 - (712/оГ) - -— Л2т( 1 - <7)4,®гу<0

" г(2^9) " ПУ0) " 1 ПУ0) >~ ш + 1

0 5

Н2Т{ 1 - (т)сг0Ж+2/ад+1

т + 1

7 т1-вс0в'а) 0 . 5Л2 т1-а с0а'а) т + 1 Г(2 - а)

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

уМ = К21ум-1 + ^22уго + ^23 уго+1 + /2, (76)

где

тахма'м + 7дг г(2_п^ ^21 = -—ГТТГ^-

К22 =

7 Г ~ 7 (~ 7 ) т1-ас(а'ст)

таяма'м + 7дг Г(2°р) + + + 0.5/г.;/г Г(2_°а)

—— Жо)

7 Т1-/Зс(/З><т) _ • , • , г1-а!с(а!>'т)

таяма'м + 7м Г(2°р) + + + 0.5/г.;/г Г(2_°а)

-0.5^маНт(ж0 - ®г0+1)

>¿23 = -

7 т 1-ес(в,ст) ~ 7 (~ 7 ) т '

таяма'м + 7дг Г(2°р) + + + 0.5/г.;/г Г(2_°а)

/2 = /2Нт - (1 - а)Нхтв1уМ - т(1 - а)жмам (уМ - уМ-1)

т 1-в с(в'ст) т 1-ас(а>ст)

+7м г^.0^) (ум " Ум-1) + М + г(2_°а) Ум

т 1-а 7-1 „ , т 1-в 7-1 , ^ /

+ 0.5Л)Л £ " Ум) " г^у Е (ЫлгГ1

«=0 у ' 7 «=0 -1-е 7-1

7-«

-(7мум)5) - 0.5Л,2г(1 - а)(1мх^уг0 + ^ _ ^ ((7мУм-1)5+1 - (7МУМ-1Г)

-0.5Н2 т (1 - ж+ уго+1

Г(2 - в) «=0 т 1-е с(в'ст) . _ . т 1-ас(а,ст) 1

Г(2 - в) 1 ' ^ . г(2 - а)

Условия устойчивости метода прогонки

А = 0, В = 0, С ^ | Аг | + г = 1, - 1, К11 < 1, К21 < 1

выполнены при а = 0.

Таким образом, с учетом (74)-(76), разностная схема (51)-(54) приводится к следующей системе линейных алгебраических уравнений:

' Агу7+1 - Сгу7+1 + - Н2тажГ^ (угО+Ч- + у^^ж+о) = ,

у0 = К11у1 + К12уго + К1зуго+1 + /¿11, ум = К21ум-1 + К22уго + ^23уго+1 + /21, ^ у° = ис(ж).

Решение системы (77) будем искать в виде

(77)

Уг = аг+1Уг+1 + /?г+1Уад + 7г+1Узд+1 + i = 0,N -1. (78)

Найдем теперь аг, вг, 7й ¿г, г = 1,^. Из условия (75) следует, что

«1 = К11, в1 = ^12, 71 = ^13, ¿1 = /11.

Подставляя yi = ai+iyi+i + вг+1Уад + 7,+ 1У,о+1 + ¿i+1, yi-1 = aiyi + вгУго + YiУад+1 + ¿i в (74), получим

В, „ Ai Bi —

io

«г+1 — -"I-, Pi+l — г< л ,

Ci — Aiai Ci — Ai ai (79)

AiYi — Fj + Ai^i

7i+i =--TZ-> dW =

Ci — Ai ai Ci - Aiai

Выразим неизвестные уг, г = 0, Ж, через уг0, у»0+1 следующим образом:

уг = Н уад + Т у^0+1 + Фг . (80)

В (80) найдем Нм, Тм, Фм • Тогда учитывая условие (76), а также ум = Нмуг0 + Тмуг0+1 + Фм, ум-1 = амум + вмуг0 + 7муг0+1 + ¿м, получим

тг К21вм + К22 „ К217м + К21^м + Д21

-пм = —-, ^М = —-, ФМ = —-•

1 - К21ам 1 - К21 ам 1 - К21ам

Найдем теперь Нг, Тг, Фг. Тогда, подставляя (80) в (78), получим

Щ = оц+Ъ = Ф* = «¿+^¿+1+&+1, г = N — 1,0. (81)

Выразим уг0, уг0+1 через Нг, Тг, Фг. Для этого рассмотрим выражения

уг0 = Нг0 уг0 + Тг0 у»0+1 + ФЗД > (82)

уг0+1 = Нг0+1уг0 + Тг0+1у»0+1 + Фг0+1. (83)

Учитывая (82), (83), получим

- я»0+1ф»0 + Фго+1(1 - нгр) _ Тг0 д, Фг0 /„.ч

Уг0+1 — п тт \(Л гр \ гр тт > УЗД — тт Уго+1 л тт '

(1 - Нг0)(1 - Тг0+1) - Т¿0Нг0+1 1 - Нг0 1 - Нг0

Из (80), (84) находим решение уг системы (77).

Литература

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

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

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

4. Podlubny I. Fractional Differential Equations.—San Diego: Academic Press, 1999.—368 p.

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

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

7. Cuesta C., van Duijn C. J., Hulshof J. Infiltration in porous media with dynamic capillary pressure: travelling waves // Eur. J. Appl. Math.—2000.—Vol. 11, № 4.—P. 381-397. DOI: 10.1017/s0956792599004210.

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

9. 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.

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

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

12. Chen P. J., Gurtin M. E. On a theory of heat conduction involving two temperatures // J. Appl. Math. Phys. (ZAMP).—1968.—Vol. 19, № 4.—P. 614-627. DOI: 10.1007/BF01594969.

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

14. Свешников А. Г., Альшин А. Б., Корпусов М. О., Плетнер Ю. Д. Линейные и нелинейные уравнения соболевского типа.—М.: Физматлит, 2007.—736 c.

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

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

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

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

19. 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.

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

21. 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.

22. Нерпин С. В., Чудновский А. Ф. Энерго- и массообмен в системе растение-почва-воздух.—Л.: Гидрометеоиздат, 1975.—358 c.

23. Нигматуллин Р. Р. Особенности релаксации системы с «остаточной» памятью // Физ. твердого тела.—1985.—Т. 27, № 5.—C. 1583-1585.

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

25. 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.

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

27. Бештоков М. Х. Краевые задачи для обобщенного модифицированного уравнения влагопереноса и разностные методы их численной реализации // Прикл. матем. и физ.—2020.—Т. 52, № 2.— С. 128-138. DOI: 10.18413/2687-0959-2020-52-2-128-138.

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

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

30. Воеводин А. Ф., Шугрин С. М. Численные методы расчета одномерных систем.—Новосибирск: Наука, Сиб. отд-е, 1981.—208 с.

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

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

Vladikavkaz Mathematical Journal 2023, Volume 25, Issue 3, P. 15-35

NUMERICAL METHODS FOR SOLVING NONLOCAL BOUNDARY VALUE PROBLEMS FOR GENERALIZED LOADED HALLAIRE EQUATIONS

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 nonlocal boundary value problems for one-dimensional space-loaded Hallaire equations with variable coefficients and two Caputo fractional differentiation operators with orders a and p. Similar problems arise in the practice of regulating the salt regime of soils, when the stratification of the upper layer is achieved by draining the water layer from the surface of a site flooded for some time. Difference schemes are constructed for the numerical solution of the problems posed on a uniform grid. Using the method of energy inequalities for various relations between the orders of the fractional Caputo derivative a and p, we obtain a priori estimates in differential and difference interpretations for solutions of nonlocal boundary value problems. The obtained a priori estimates imply uniqueness and stability of the solution with respect to the right-hand side and the initial data, as well as the convergence of the solution of the difference problem to the solution of the corresponding original differential problem (assuming the existence of a solution to the differential problem in the class of sufficiently smooth functions) at a rate of O(h2 + t2) for a = p and O(h2 + t2-maxia>^}) for a = p. The paper also presents an algorithm for the numerical solution of a nonlocal boundary value problem for a loaded Hallaire equation with variable coefficients and a Bessel operator.

Keywords: nonlocal boundary value problems, a priori estimate, Hallaire equation, loaded equations, generalized moisture transfer equation, Caputo fractional derivative.

AMS Subject Classification: 65N06, 65N12.

For citation: Beshtokov, M. Kh. Numerical Methods for Solving Nonlocal Boundary Value Problems for Generalized Loaded Hallaire Equations, Vladikavkaz Math. J., 2023, vol. 25, no. 3, pp. 15-35 (in Russian). DOI: 10.46698/c8748-9711-0633-d.

References

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

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

3. 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).

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

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

6. 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).

7. Cuesta, C., van Duijn, C. J. and Hulshof, J. Infiltration in Porous Media with Dynamic Capillary Pressure: Travelling Waves, European Journal of Applied Mathematics, 2000, vol. 11, no. 4, pp. 381-397. DOI: 10.1017/s0956792599004210.

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

9. 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.

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

11. 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).

12. 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.

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

14. Sveshnikov, A. G., Alshin, A. B., Korpusov, M. O. and Pletner, Yu. D. Lineynyye i nelineynyye uravneniya sobolevskogo tipa [Linear and Nonlinear Equations of the Sobolev Type], Moscow, Fizmatlit, 2007, 736 p. (in Russian).

15. 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).

16. Goloviznin, V. M., Kiselev, V. P. and Korotkiy, I. A. Chislennye metody resheniya uravneniya drobnoy diffuzii s drobnoy proizvodnoy 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).

17. Lafisheva, M. M., Shhanukov-LaBshev, M. H. Locally One-Dimensional Difference Schemes for the Fractional Order Diffusion Equation, Comput. Math. Math. Phys., 2008, vol. 48, no.10, pp. 1875-1884.

18. 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.

19. 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/S0012266110050058.

20. 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.

21. Nerpin, S. V. and Chudnovskiy, A. F. Energo- i massoobmen v sisteme rasteniye-pochva-vozdukh [Energy and Mass Transfer in the Plant-Soil-Air System], Leningrad, Gidrometeoizdat, 1975, 358 p. (in Russian).

22. Nigmatulin, R. R. Relaxation Features of a System with Residual Memory, Fizika Tverdogo Tela [Physics of the Solid State], 1985, vol. 27, no. 5, pp. 1583-1585 (in Russian).

23. 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.

24. 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.

25. Beshtokov, M. Kh. To Boundary Value Problems for Degenerating Pseudoparabolic Equations with Gerasimov-Caputo Fractional Derivative, Russian Mathematics, 2018, vol. 62, pp. 1-14. DOI: 10.3103/S1066369X18100018.

26. Beshtokov, M. Kh. Boundary Value Problems for the Generalized Modified Moisture Transfer Equation and Difference Methods for Their Numerical Realization, Prikladnaya Matematika i Fizika [Applied Mathematics and Physics], 2020, vol. 52, no. 2, pp. 128-138 (in Russian). DOI: 10.18413/2687-09592020-52-2-128-138.

27. 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).

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

29. Voyevodin, A. F. and Shugrin, S. M. Chislennyye metody rascheta odnomernykh sistem [Numerical Methods for Calculating One-Dimensional Systems], Novosibirsk, Nauka, Siberian Branch, 1981, 208 p. (in Russian).

Received April 18, 2022 Murat Kh. Beshtokov

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