Вестник КРАУНЦ. Физ.-мат. науки. 2020. Т. 32. № 3. C. 15-28. ISSN 2079-6641
МАТЕМАТИКА
УДК 519.622.2 Научная статья
Численное решение краевой задачи для нагруженного параболического уравнения с нелокальными граничными
условиями
В. М. Абдуллаев1'2
1 Азербайджанский государственный университет нефти и промышленности, AZ1010, г. Баку, пр-т Азадлыг, 20
2 Институт систем управления НАН Азербайджана, AZ114, г. Баку, ул. Вахабзаде, 9
E-mail: vaqif_ab@rambler.ru
В работе с использованием метода прямых исследуется численное решение краевой задачи относительно нагруженного параболического уравнения с нелокальными краевыми условиями. Получены расчетные формулы и приводится алгоритм для решения задачи. Приведены результаты численного решения двух тестовых задач, иллюстрирующие эффективность предложенного подхода.
Ключевые слова: нагруженные дифференциальные уравнения, нелокальные условия, метод прямых, перенос краевых условий.
DOI: 10.26117/2079-6641-2020-32-3-15-28
Поступила в редакцию: 05.08.2020 В окончательном варианте: 28.09.2020
Для цитирования. Паровик Р. И. Численное решение краевой задачи для нагруженного параболического уравнения с нелокальными граничными условиями // Вестник КРАУНЦ. Физ.-мат. науки. 2020. Т. 32. № 3. C. 15-28. DOI: 10.26117/2079-6641-2020-32-3-15-28
Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru)
© Абдуллаев В.М., 2020
Введение
Дифференциальные уравнения, в которых участвуют значения искомой функции в некоторых заданных точках, поверхностях (точечные, поверхностные нагружения) называют нагруженными [1]-[4]. Нагруженные дифференциальные уравнения описывают процессы, в которых их состояние в какой-либо точке и/или в какой-либо момент может оказывает влияние на весь процесс в целом. Такими уравнениями описываются процессы в биологии, экологии, подземной гидрогазодинамики. Нелокальность условий обусловлена практической невозможностью производить замеры или воздействовать на процесс в отдельно взятой точке распределенного в пространстве объекта или мгновенно во времени [5]-[7].
Финансирование. Исследование выполнялось без финансирования
К нагруженным уравнениям приводят также математические постановки задач оптимального управления с обратной связью [8], в которых обратная связь происходит дискретно в отдельные моменты времени или для объектов с распределенными в пространстве параметрами с отдельными его точками [9, 10].
В работах [11, 12] исследованы вопросы существования, единственности решения различных классов задач, описываемых нагруженными дифференциальными уравнениями. В работах [7, 12] для нагруженных систем, описываемых в достаточно общем виде уравнениями с обыкновенными и частными производными параболического типа, исследованы разностные аппроксимирующие задачи, получены априорные оценки для их решений и доказана их сходимость к решению непрерывных задач. В работах [7, 12] для решения разностных уравнений при нагружении по пространственной переменной предлагается использовать методы прогонки [13], но каких-либо конкретных схем и результатов численного решения не приводится.
В работе исследовано применение различных схем конечноразностной аппроксимации краевой задачи относительно точечно нагруженного параболического уравнения [14]. В работах [7, 15] исследовалось численное решения нагруженных параболических уравнений при нагружении по пространственной переменной и была предложена схема решения соответствующих конечноразностных аппроксимирующих задач.
В работах [11, 16, 17, 18] исследовались проблемы, связанные с численным решением обратных задач и задач оптимального управления процессами, описываемых нагруженными дифференциальными уравнениями.
В данной работе на примере нагруженного параболического уравнения исследуется подход к численному решению соответствующих краевых задач, основанный на применении схем метода прямых.
Постановка задачи
Рассматривается следующая краевая задача, описываемая нагруженным уравнением параболического типа:
^^ - ^"uXr1 - tt) = f (x, t), (x, t) e П; (1)
при начальном условии
u(x, 0) = ф(x) , 0 < x < 1, (2)
и нелокальных граничных условиях:
an (t)u(0, t) + ai2(t)^^Of1 + 711 (t)u(1, t) + Yi2(t)= ^'), 0 < t < T, (3)
«21 (t)u(0,t) + a22(t)^^^ + Y2i(t)u(1,t) + Y22(t)= &(t), 0 < t < T, (4)
Здесь u(x, t)-искомая функция, определенная в области П = {0 < x < 1, 0 < t < T}; t и x являются соответственно временной и фазовой координатами; T и 1 — заданные положительные числа; xs e (0, 1),s = 1,2,..., k — заданные точки нагруже-ния; f (x,t), Bs(x,t), s = 1,2,...k, «j(t), в^ (t), Yj(t), г, j = 1,2 — заданные непрерывно-дифференцируемые по своим аргументам функции на отрезке [0, T], а ф(x) — заданная непрерывно-дифференцируемая функция на отрезке [0,1].
Требуется определить непрерывно-дифференцируемую по г е [г0, Т] и дважды непрерывно-дифференцируемую по х е (0,1) функцию и(х,г), удовлетворяющую условиям
(1) - (4).
В работе [12] исследованы вопросы существование, единственности решения для нагруженного параболического уравнения для различных видов краевых условий. Приведены схемы конечно-разностной аппроксимации для нагруженных дифференциальных уравнений параболического типа, показана их сходимости к точному решению. Подобные схемы исследования могут быть использованы и для конкретно рассматриваемой задачи.
Метод решения
На отрезке [0, T] возьмем точки tj = jT, т = T/N и проведем прямые t = tj,
j = 0,1,...,N. Обозначая через Uj(x) = w(x,tj)и заменяя ЩхЛ | разностным отношением, получим следующие нагруженные обыкновенные дифференциальные уравнения второго порядка:
d2UXi(Xl = 1 Uj(x) - ЕВ(x)Uj(xs)-
dx T s=1
-1 Uj-1(x) -fj(x), j = 1,2,...,N, U0(x) = ф(х). (5)
В силу (3)-(4) для этой системы уравнений имеются нелокальные граничные условия:
<Uj (0) + а/2 j + 7л Uj (1) + j ^U^"11 = #, j = 1,2,..., N, (6)
<Uj(0) + aj+ ¿Uj(1) + Y222^ = , j = 1,2,...,N. (7)
*—) + у и-(1) + У ^)
¿х + ^ (1) + Г22 ¿х
Для численного решения задачи (2)-(4) предлагается использовать следующий подход, использующий операцию сдвига неразделенных условий [19, 20], который обобщает известную процедуру переноса краевых условий [21].
Определение. Будем говорить, что функции «/(х) , о>(х), 'х), ?2(х), в/ (х), я = 1,2,...,к, ^'(х),' = 1,2,...,N осуществляют сдвиг краевых условий (3), (4), если для решения нагруженной системы (2) и'(х) в любой точке х отрезка [0,1 ] имеет место:
а/ (х)и'(х) + а{ (х) ^^ + £ 0/'(х) '/)+
"х /=1
+7/ (х)и-(1) + (х) ^ = <§'(х), -= 1,2,..., N. (8)
Здесь: а{(х), а2(х),у[(х), (х),в/(х), / = 1,2,...,к, £''(х),./ = 1,2,...,N — произвольные вектор функции, удовлетворящие условиям:
Пы- а{\ (0)
а1 (0) =
а21(0)
= aj,
а2 (0) =
а12 J
(0) а22(0)
а12 7
а;
22
= а2,
(9)
М Й - ■ = 1' М Й)М ■ I-■ (10)
ß/(0)4 ßs2i 10) I - ( 0 ) ' s - 1'2'...'k Ij(0)- ( ^(0) ) -( I1 I - I■. (Щ
Несложно проверить, что соотношения (8) при х = 0 в силу условий (9)---(7)
совпадают с краевыми условиями (3), (4).
Функции а/(х), «2(х),?1 (х), (х),в/(х), £'(х), 5 = 1,2,...,к,у' = 1,2,...,N, позволяющие осуществить сдвиг краевых условий, определяются не единственным образом. Один из возможных вариантов получения функций, осуществляющих такой перенос, предлагается в нижеприводимой теореме.
Теорема. Пусть функции а/(х), а2(х), 'х), )2(х),в/(х), £у'(х), / = 1,2,...,к, у 1,2,...,N определяются из решения следующих задач Коши и соотношений:
= ^ (х) «1 (х) - - «2 (х), «1 (0) = «1, (12)
ах 0 1 т 2 1 1
^a2^ - S0(x)a2 (x) - a/ (x) «2 (0) - «2, (13)
ил
d ß/ (x)
Ux
- S0(x)ßi(x) - «2(x)BS(x), ß/(0) - 0's - 1'2...'k' (14)
j- - S0(x)|i(x) + «2 (x)(-1 Ui-1(x) - jx))' |j(0) - (15)
dm(x) - S0 (x)mj(x)' mj(0) - I' (16)
Ux
y/(x) - mi(x)y1' y2(x) - m7'(x)Y2'
где
S0 (x) -
1 ■ ■ * k
- a2 (x) a/* (x) + a/ (x) a2* (x) + a/ (x) £ B/ (x)ß/ (x)
T s-1
-Ii(x)a2'(x)(-1 j (x) -jx))] [R(x)]-1'
L
k
R(x) - a/(x)a/*(x) + a2(x)«2*(x) + £ß/(x)ß/*(x) + jx)|j*(x)'
s-1
S0(x) e E2, j = 1,2,..,N, I— 2-мерный вектор, составленный из единиц. Тогда для функции Uj ^—решения задачи (2)-(4) в каждой точке x e [0,1] имеет место (8), причем R(x) = const = R(x0), x e [0,1 ].
Доказательство. Пусть для Uj(x), j = 1,2,...,n, являющихся решением задачи (2)-(4), имеет место соотношение (8), где a1(x), a^(x), y1(x), Y2(x),e/(x), £j(x),
/ = 1,2,...,к, ' = 1,2,...,N, пока произвольные дифференцируемые функции, удовлетворяющие лишь условиям (9) - (7). Поддиференцируем (8)
da( (x)
dx
Uj (x) + a{ (x)
dUj (x) da/ (x) dUj (x)
dx
dx dx
+ ajx)
d 2Uj (x)
dx2
+ £ dß/ (x) Uj (x) + Uj (1) + d(x) dUj (1) _ d<§ j (x)
s=1
dx
dx
dx dx
dx
Учитывая (2), произведем группировку, в результате получим:
da1 (x) 1 т
dx
+ ^ а/ (x)
Uj (x) + £
s=1
+ a2 (x)BS (x)
Uj (xs)+
+
d a2(x) dx
+ a{ (x)
dUj (x) dx
j - a2 (x)( -1 Uj-1(x) - fj (x)
+djx) Uj (1)+diM dUj (i) = 0.
+
их их их
Предполагая, что и'(х) ^0, и пользуясь произвольностью функций а1(х),а2(х), (х),)^(х),в/(х),^'(х), / = 1,2,...,к, потребуем равенства нулю выражений в квадратных скобках:
da{ (x) 1 j da2 (x)
= --a2 (x), —= -aj (x), dx
dx
d ßS (x) dx
= -a^(x)B/ (x), s = 1,2,...,k,
^ = -a2(x)(-1 Uj-1(x) - fj(x)) ,
dY1(x)= 0 dr2 (x)
= 0.
их их
Умножая обе части (8) на произвольную функцию т'(х) такую, что
будем иметь:
mj (0) = I, dUj (x)
mj (x)a!(x)U-/ (x) + mj (x)a^(x)
+ mj (x) £ ß/ (x) Uj (xs)+
S=1
+т' (х) у( (х) и' (1) + т' (х) ^ (х) = т' (х) ' (х), ' = 1, 2,..., N.
ил
Введем обозначения:
(17)
(18)
(19)
(20)
(21)
(22)
g (x) = mj(x) a{ (x), (x) = mj (x) a^ (x), q/ (x) = mj(x)ß/ (x), s = 1, 2,..., k,
(23)
f
г{ (х) = ту (х) (х), г^(х) = ту (х) ^ (х), к' (х) = ту (х) %' (х), (24)
после чего (21) примет вид:
^(х)^ (х) + ¿(х) + Е (х) и' (х/)+ +г-(x)U•/' (1) + г2(х) = ку' (х), у' = 1,2,..., N. (25)
ил
Причем из (9) - (7), (15) - (19) следует
£1(0) = «1, ¿(0) = «2, (0) = 0, / = 1,2...,к,
ку' (0) = % ', г|(0) = ^, г2(0) = , У = 1,2,..., N. Потребуем от ту(х)выполнения условия
k
а/ (x) а/* (x) + aj (x) «2* (x) + £ в/(x)в,7'* (x) + % (x)%j* (x) ) = con/t, x g [0, l]. (26)
s=1
Продифференцировав (26) и (5),(6) и учтя (15) - (19), получим: ^(х). £ (х) + £1(х) и£1:(х)+ ' £>(х) + 'х) и«2"(х)
dx 1 1 dx dx 2 2 dx
k
+ £
s=1
dq/(x) ■ dq/ *(x) q/ (x)
dx q/ (x)+ dx q/(x)
+ ' ' (x) + k'(x) = 0, (27)
dx dx
^ = '(x) + m'(x)^ = ' (m'(x))-1 (x) -m'(x)1 а2(x) = dx dx dx dx т
= ' (m'(x))-1 *iM - 1«2(x), (28)
= / «j Ci)+m Cl) =
dx dx dx
= ^mF (m/'(x))-1 gi(x) -'X(x) = ^md^ (mj(x))-1 gj(x) -g1(x), (29) ^ = 'в/(x) + m'(x)' = ' (m'(x))-1 q,(x) -m/(x)a2(x)B,(x) =
= ^mdr^ (mj(x))-1 q/'(x) -gj(x)B'(x), (30)
dkj (x) dmj (x) % / / d%j (x)
dx dx dx
= ^т^ ^(х))-1 к'(х) + ^(х)(-1 и'-1 -/'(х)) . (31)
Подставив (24) - (31) в (23), после несложных преобразований и группировок
получим:
dm^ (mj(x))-M g1(x)gf(x) + (x) j(x) + £q/(x)qf (x) + kj(xjx) ) +
S=1
+1 g2 (x)gi* (x) + (x)gf (x) +
L
k /1 +g1(x) £ B/ (x)q/ (x) - kj (x)g2(xW - - Uj-1(x) - /j (x)
s=1 ^ T
+
(т'«Г1 (х)^Г(х)+ ^(х)'(х)+ £«/'(х)^'* (х) + к'(х)к>(х^ + +1 (х) + ^ (х)^* (х)+
С
+^1(х) £в/(х)в/(х) -'х)'х) (-1 и'-1(х) -/(х)) * = 0
Пользуясь произвольностью функции т (х), принимем значение выражения в квадратных скобках равным нулю, тогда
^^ Их))-1 = 50(х), (32)
где
1k S0(x) = ~zg2(x)g1+(x) + g1(x)gJ2*(x) + g (x) £ BS (x)q/ (x) -. s=1
-kj(x)g2(x^-1 Uj-1(x) -/j(x)^J [Rj(x)]-1, Rj (x) = g (x)jx) + (x) jx) + £q/ (x)qf(x) + kj (xjx).
s= 1
j
Подставляя (33) в (24)-(31) и учитывая (5), (6), переобозначив функции ^(х), (х), г1(х), г2(х), q/(х), к'(х), / = 1,2,...,к, соответственно через о^(х), а>(х), (х), т2(х) в/(х), ^'(х), / = 1,2,...,к, получим утверждение теоремы.
Из этой теоремы следует, что решение задачи (2) - (4) можно привести к решению задач Коши относительно системы (8) - (12), состоящей из 2(к + 4) дифференциальных уравнений. В результате их решения, запомнив значения функций а1 (х), а2(х), (х), (х), в/(х), ^'(х), / = 1,2,...,к, в точках нагружения
х1,х2...,хк и х/ получаем линейную систему алгебраических уравнений порядка 2(к+
1) относительно Ц'(*), Ц'Ы,..,^(хЛ),, ^^г21,...,^Цах^ Ц(/), ии&х(/1.
Таким образом, предлагаемый подход к решению задачи (1)-(4) заключается в следующем. Применяя метод прямых исходная задача приводится к решению на каждом временном слое г/ = 'т нагруженного дифференциального уравнения второго порядка (2) с нелокальными условиями (3), (4). Для их решения сначала каким-либо численным методом решаются задачи Коши (8)-(12) и определяются значения функций «'(х), «>(х), 'х), т2(х), в/(х), %'(х), 5 = 1,2,...,к, в точках нагружения х1,х2,...,хк и х = /. Осуществляя сдвиг условий (3), (4) посредством сооотношения (8) в точки нагружения х1, х2, ... , хк и х = /, получим систему алгебраических уравне-
ни порядка 2(к +1) относительно 2(к + 1) неизвестных и'(х1), и-/'(х2),...,и-/'(хк),^х^ ^,
ииихх2),...,иииххк),и-/'(/), иЦх(/). Для окончательного определения функции и'(х), необходимо каким-либо численным методом решить задачу Коши (1)-(4) в обратном во времени направлении от х = / до х = 0, учитывая найденные значения в точках нагружения и'(х1),...,и'(хк) и и'(/).
Результаты численных экспериментов
Применим предложенный подход к следующей тестовой задаче. Задача 1.
Эф^О = д Ф,г) + ги(0.25,г) + (г- 1)и(0.5,г) + х-0.75?2- 1.81257?-0.75,
(х,г) е О = {(х,г) :0 < х < 1, 0 < г < 1}, при начальном и нелокальных граничных условиях:
и(х, 0)= х2 + 1, 0 < х < 1, ди(0,г) = -2и(0,г) + и(1,г),
д х
ди(1,г) =(2г + 4)и(0, г) - м(1,г), 0 < г < 1. д х
Точками нагружения являются: хх = 0.25, х2 = 0.5. Несложно проверить, что функция и(х,г) = х2 + хг + 1 является точным решением рассмотриваемой задачи.
В табл. 1 приведены полученные приближенные й(х,г) и точные значения и*(х,г). Жирным шрифтом выделены значения функции в точках нагружения. Численные эксперименты проведены с использованием метода прямых. Временной отрезок [0; 1] разделен на 8 частей, т = 1/8. Решение задач Коши проводилось методом Рунге-Кутта четвертого порядка при числе разбиений, равном 200 (йх = 0.005), а для решения алгебраической системы применен метод Гаусса с выбором главного члена.
Таблица 1
Полученное (ñ(x,t)) и точно (w*(x,t)) решения задачи 1.
t x ñ(x, t) и*(х, t) t x ñ(x, t) M*(x, t)
0.125 0.00 1.0003 1.0000 0.625 0.00 0.9992 1.0000
0.25 1.0935 1.0938 0.25 1.2180 1.2188
0.50 1.3122 1.3125 0.50 1.5618 1.5625
1.00 2.1252 2.1250 1.00 2.6244 2.6250
0.25 0.00 0.9998 1.0000 0.75 0.00 0.9990 1.0000
0.25 1.1248 1.1250 0.25 1.2490 1.2500
0.50 1.3748 1.3750 0.50 1.6241 1.6250
1.00 2.2497 2.2500 1.00 2.7491 2.7500
0.375 0.00 0.9997 1.0000 0.875 0.00 0.9987 1.0000
0.25 1.1559 1.1563 0.25 1.2801 1.2813
0.50 1.4371 1.4375 0.50 1.6863 1.6875
1.00 2.3648 2.3750 1.00 2.8837 2.8750
0.50 0.00 0.9994 1.0000 1.00 0.00 0.9985 1.0000
0.25 1.1869 1.1875 0.25 1.3111 1.3125
0.50 1.4993 1.5000 0.50 1.7483 1.7500
1.00 2.4994 2.5000 1.00 2.9974 3.0000
Задача 2. Пусть процесс описывается следующим нагруженным дифференциальным уравнением:
о>ф,0 = d ^t) + (x + i)costu(0.2,t) + (2x + 1) sintu(0.6,t) + 2u(0.8,t)+
+ (x2 - x - 1) cos t - 0.02(x + 1) sin2t - (0.72x + 0.36) sin21 - (2x + 2.28) sint - 2, (x,t) e П = {(x,t) : 0 < x < 1, 0 < t < 1}, при следующих начальном и нелокальных граничных условиях:
u(x, 0) = 1, 0 < x < 1,
2tu(0, t) + u(1, t) + ^^Í^líi + 3= 7 sint + 2t + 1, 0 < t < 1, dx dx
и(0, t) + (t + 1)и(1, t) + 2^^^ +1 Mbil = (3t + 1) sint +1 + 2, 0 < t < 1
dx dx
Несложно проверить, что все необходимые условия на исходные данные, функции, участвующие в постановке задачи, выполнены, а точным решением задачи является в отличие от решения задачи 1 более сложная функция u(x,t) = x2sint + 1. Пространственные нагружения сосредоточены в трёх точках x1 = 0.2, x2 = 0.6, x3 = 0.8. В таблице 2 приведены полученные приближенные ñ(x,t) и точные значения w*(x,t). Жирным шрифтом выделены значения функции в точках нагружения.
Таблица 2
Полученное (й(х,t)) и точное (w*(x,t)) решения задачи 2
t x й(х, t) м*(х, t) t x U(x, t) M*(x,t)
0.1 0.1 1.0010 1.0010 0.6 0.1 1.0056 1.0056
0.2 1.0040 1.0040 0.2 1.0225 1.0226
0.3 1.0089 1.0090 0.3 1.0507 1.0508
0.4 1.0159 1.0160 0.4 1.0903 1.0903
0.5 1.0249 1.0250 0.5 1.1412 1.1412
0.6 1.0359 1.0359 0.6 1.2036 1.2033
0.7 1.0488 1.0489 0.7 1.2774 1.2767
0.8 1.1273 1.0639 0.8 1.3626 1.3614
0.9 1.0808 1.0809 0.9 1.4595 1.4574
0.2 0.1 1.0020 1.0020 0.7 0.1 1.0064 1.0064
0.2 1.0079 1.0079 0.2 1.0257 1.0258
0.3 1.0178 1.0179 0.3 1.0579 1.0580
0.4 1.0317 1.0318 0.4 1.1031 1.1031
0.5 1.0496 1.0497 0.5 1.1612 1.1611
0.6 1.0715 1.0715 0.6 1.2324 1.2319
0.7 1.0974 1.0973 0.7 1.3165 1.3157
0.8 1.1273 1.1271 0.8 1.4139 1.4123
0.9 1.1613 1.1609 0.9 1.5243 1.5218
0.3 0.1 1.0029 1.0030 0.8 0.1 1.0071 1.0072
0.2 1.0117 1.0118 0.2 1.0286 1.0287
0.3 1.0265 1.0266 0.3 1.0645 1.0646
0.4 1.0472 1.0473 0.4 1.1148 1.1148
0.5 1.0738 1.0739 0.5 1.1796 1.1793
0.6 1.1064 1.1064 0.6 1.2588 1.2582
0.7 1.1450 1.1448 0.7 1.3526 1.3515
0.8 1.1895 1.1891 0.8 1.4609 1.4591
0.9 1.2402 1.2394 0.9 1.5840 1.5811
0.4 0.1 1.0038 1.0039 0.9 0.1 1.0078 1.0078
0.2 1.0155 1.0156 0.2 1.0313 1.0313
0.3 1.0349 1.0350 0.3 1.0705 1.0705
0.4 1.0622 1.0623 0.4 1.1254 1.1253
0.5 1.0973 1.0974 0.5 1.1961 1.1958
0.6 1.1403 1.1402 0.6 1.2827 1.2820
0.7 1.1911 1.1908 0.7 1.3851 1.3838
0.8 1.2499 1.2492 0.8 1.5034 1.5013
0.9 1.3167 1.3154 0.9 1.6378 1.6345
0.5 0.1 1.0047 1.0048 1.0 0.1 1.0084 1.0084
0.2 1.0191 1.0192 0.2 1.0336 1.0337
0.3 1.0431 1.0431 0.3 1.0758 1.0757
0.4 1.0766 1.0767 0.4 1.1348 1.1346
0.5 1.1199 1.1199 0.5 1.2107 1.2104
0.6 1.1728 1.1726 0.6 1.3037 1.3029
0.7 1.2354 1.2349 0.7 1.4137 1.4123
0.8 1.3078 1.3068 0.8 1.5409 1.5385
0.9 1.3900 1.3883 0.9 1.6852 1.6816
ISSN 2079-6641
Заключение
В работе предложен численный метод решения краевой задачи относительно нагруженного параболического уравнения с нелокальными краевыми условиями. Получены формулы и приводится алгоритм для решения задачи. Предложенный подход основан на методе прямых для аппрокцимации дифференциального уравнения и операции сдвига условий, развивающем известный метод переноса краевых условий. Подход позволяет свести решение исходной задачи к решению на каждом временном слое (прямой) задач Коши относительно обыкновенных систем дифференциальных уравнений и системы линейных алгебраических уравнений. Приведены результаты численного решения двух тестовых задач, иллюстрирующие эффективность предложенного метода.
Конкурирующие интересы. Конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответсвенность. Автор участвовал в написании статьи и полностью несет ответственность за предоставление окончательной версии статьи в печать.
Список литературы/References
[1] Нахушев А.М., Уравнения математической биологии, Высшая школа, М., 1995,
305 с. [Nahushev A.M., Uravnenija matematicheskoj biologii, Vvsshaja shkola, Moskva, 1995, 305 pp.]
[2] Нахушева В. А., Дифференциальные уравнения математических моделей нелокальных процессов, Наука, М., 2006, 173 с. [Nahusheva V. A., Differencial'nye uravnenija matematicheskih modelej nelokal'nyh processov, Nauka, M., 2006, 173 pp.]
[3] Нахушев А. М., Нагруженные уравнения и их применение, Наука, М., 2012, 232 с. [Nahushev A.M., Nagruzhennve uravnenija i ih primenenie, Nauka, M., 2012, 232 pp.]
[4] Нахушев А. М., Борисов В. Н., "Краевые задачи для нагруженных параболических уравнений и их приложения к прогнозу уровня грунтовых вод", Дифференц. уравнения, 13:1 (1977), 105-110. [Nahushev A.M., Borisov V. N., "Kraevve zadachi dlja nagruzhennvh parabolicheskih uravnenij i ih prilozhenija k prognozu urovnja gruntovvh vod", Differenc. uravnenija, 13:1 (1977), 105-110].
[5] Худалов М.З., "Нелокальная краевая задача для нагруженного уравнения параболического типа", Владикавк. матем. журнал, 4:4 (2002), 59-64. [Hudalov M.Z., "Nelokal'naja kraevaja zadacha dlja nagruzhennogo uravnenija parabolicheskogo tipa", Vladikavk. matem. zhurn., 4:4 (2002), 59-64].
[6] Дикинов Х.Ж., Керефов А. А., Нахушев А.М., "Об одной краевой задаче для нагруженного уравнения теплопроводности", Дифференц. уравнения, 12:1 (1976), 177-179. [Dikinov H.Zh., Kerefov A.A., Nahushev A.M., "Ob odnoj kraevoj zadache dlja nagruzhennogo uravnenija teploprovodnosti", Differenc. uravnenija, 12:1 (1976), 177-179].
[7] Шхануков-Лафишев М.Х., "Локально-одномерная схема для нагруженного уравнения теплопроводности с краевыми условиями III рода", Ж. вычисл. матем. и математической физики, 49:7 (2009), 1223-1231. [Shhanukov-Lafishev М.Н., "Lokal'no-odnomernaja shema dlja nagruzhennogo uravnenija teploprovodnosti s kraevvmi uslovijami III roda", Zh. vvchisl. matem. i matematicheskoj fiziki, 49:7 (2009), 1223-1231].
[8] Егоров А. И., Основы теории управления, Физматлит, М., 2004, 504 с. [Egorov A.I., Osnovv teorii upravlenija, Fizmatlit, M, 2004, 504 pp.]
[9]
Distributed-Parameter Processes", Autom. Remote Control, 73:9 (2012), 1443-1455.
of the heating process control", Autom. Remote Control, 78:9 (2017), 1585-1599.
[11] Кожанов А. И., "Нелинейные нагруженные уравнения и обратные задачи", Ж. вычисл. матем. и математической физики, 44:4 (2004), 694-716. [Kozhanov A.I., "Nelinejnve nagruzhennve uravnenija i obratnve zadachi", Zh. vvchisl. matem. i matem. fiz., 44:4 (2004), 694-716].
[12] Алиханов А. А., Березков А.М., Шхануков-Лафишев М.Х., "Краевые задачи для некоторых классов нагруженных дифференциальных уравнений и разностные методы их численной реализации", Ж. вычисл. матем. и математической физики, 48:9 (2008),
1619-1628. [Alihanov A. A., Berezkov A.M., Shhanukov-Lafishev М.Н., "Kraevve zadachi dlja nekotorvh klassov nagruzhennvh differencial'nyh uravnenij i raznostnve metodv ih chislennoj realizacii", Zh. vvchisl. matem. i matematicheskoj fiziki, 48:9 (2008), 1619-1628].
[13] Самарский А. A., Николаев Е. С., Методы решения сеточных уравнений, Наука, М.,
1978, 532 с. fSamarskij A. A., Nikolaev Е. S., Metodv reshenija setochnvh uravnenij, Nauka, M., 1978, 532 pp.]
[14]
equations", Comput. Math. Math. Phvs., 56:1 (2016), 93-105.
[15] Бондарев Э.А., Воеводин А.Ф., "Разностный метод решения начально-краевых задач для нагруженных дифференциальных и интегро-дифференциальных уравнений", Диф-ференц. уравнения, 36:11 (2000), 1560-1562. [Bondarev Je. A., Voevodin A.F., "Raznost-nvj metod reshenija nachal'no-kraevyh zadach dlja nagruzhennvh differencial'nyh i integro-differencial'nyh uravnenij", Differenc. uravnenija, 36:11 (2000), 1560-1562].
[16] Дженалиев М.Т., "Оптимальное управление линейными нагруженными параболическими уравнениями", Дифференц. уравнения, 25:4 (1989), 641-651. [Dzhenaliev М.Т., "Optimal'noe upravlenie linejnvmi nagruzhennvmi parabolicheskimi uravnenij ami", Differenc. uravnenija, 25:4 (1989), 641-651].
[17]
Loaded Lumped Systems", Comput. Math. Math. Phvs., 46:9 (2006), 1487-1502.
[18]
ential equations", Vestnik KRAUNC. Fiz.-mat. nauki, 22:2 (2018), 33-44.
[19]
boundarv-value problems for ordinarv differential equations", Comput. Math. Math. Phvs., 54:7 (20l4), 1096-1109.
[20]
dinarv Differential Equations", Comput. Math. Math. Phvs., 44:9 (2004), 1585-1595.
[21] Абрамов А. А., "О переносе граничных условий для систем линейных обыкновенных дифференциальных уравнений (вариант метода прогонки)", Ж. вычисл. матем. и матем. физ., 1:3 (1961), 542-545. [Abramov А. А., "О perenose granichnvh uslovij dlja sistem linejnvh obvknovennvh differencial'nyh uravnenij (variant metoda progonki)", Zh. vvchisl. matem. i matem. fiz., 1:3 (1961), 542-545].
Список литературы (ГОСТ)
[1] Нахушев А. М. Уравнения математической биологии. М.: Высшая школа. 1995. 305 с.
[2] Нахушева В. А. Дифференциальные уравнения математических моделей нелокальных процессов. М.: Наука. 2006. 173 с.
[3] Нахушев А. М. Нагруженные уравнения и их применение. М.: Наука. 2012. 232 с.
[4] Нахушев А.М., Борисов В. Н. Краевые задачи для нагруженных параболических уравнений и их приложения к прогнозу уровня грунтовых вод // Дифференц. уравнения. 1977. Т.13. №1. С. 105-110.
[5] Худалов М. З. Нелокальная краевая задача для нагруженного уравнения параболического типа // Владикавк. матем. журнал. 2002. Т.4, №4. С. 59-64.
[6] Дикинов Х. Ж., Керефов А. А., Нахушев А. М. Об одной краевой задаче для нагруженного уравнения теплопроводности // Дифференц. уравнения. 1976. Т. 12, №1. С. 77-79.
[7] Шхануков-Лафишев М. Х. Локально-одномерная схема для нагруженного уравнения теплопроводности с краевыми условиями III рода // Ж. вычисл. матем. и математической физики. 2009. Т. 49, №7. С. 1223-1231.
[8] Егоров А. И. Основы теории управления. М.: Физматлит. 2004. 504 с.
[9] Aida-zade K. R., Abdullaev V. M. On an Approach to Designing Control of the Distributed-Parameter Processes // Autom. Remote Control. 2012. Vol.73, No.9. Pp. 1443-1455.
[10] Aida-zade K. R., Abdullaev V. M. Optimizing placement of the control points at synthesis of the heating process control // Autom. Remote Control. 2017. vol. 78. No.9. pp. 1585-1599.
[11] Кожанов А. И. Нелинейные нагруженные уравнения и обратные задачи // Ж. вычисл. матем. и математической физики. 2004. T.44, №4. С.694-716.
[12] Алиханов А. А., Березков А. М., Шхануков-Лафишев М. Х. Краевые задачи для некоторых классов нагруженных дифференциальных уравнений и разностные методы их численной реализации // Ж. вычисл. матем. и математической физики. 2008. T. 48. № 9. С.1619-1628.
[13] Самарский А. A., Николаев Е. С. Методы решения сеточных уравнений. М.: Наука, 1978.
[14] Abdullayev V. M., Aida-zade K. R. Finite-difference methods for solving loaded parabolic equations // Comput. Math. Math. Phys. 2016. vol.56. No. 1. P. 93-105.
[15] Бондарев Э. А., Воеводин А. Ф. Разностный метод решения начально-краевых задач для нагруженных дифференциальных и интегро-дифференциальных уравнений // Диффе-ренц. уравнения. 2000. T.36. №11. С. 1560-1562.
[16] Дженалиев М. Т. Оптимальное управление линейными нагруженными параболическими уравнениями // Дифференц. уравнения. 1989. T. 25. №4. С.641-651.
[17] Abdullaev V. M., Aida-zade K. R. Numerical Solution of Optimal Control Problems for Loaded Lumped Systems // Comput. Math. Math. Phys. 2006. vol. 46, No.9. P.1487-1502.
[18] Abdullayev V. M. Numerical solution to parametric identification problems partial differential equations // Vestnik KRAUNC. Fiz.-mat. nauki. 2018. vol. 22. No.2. P. 3344. DOI: 10.18454/2079-6641-2018-22-2-33-44
[19] Abdullaev V. M., Aida-zade K. R. Numerical method of solution to loaded nonlocal boundary-value problems for ordinary differential equations // Comput. Math. Math. Phys. 2014. vol. 54, No.7. P.1096-1109.
[20] Abdullaev V. M., Aida-zade K. R. On the Numerical Solution of Loaded Systems of Ordinary Differential Equations // Comput. Math. Math. Phys. 2004. vol. 44, No. 9. P. 1585-1595.
[21] Абрамов А. А. О переносе граничных условий для систем линейных обыкновенных дифференциальных уравнений (вариант метода прогонки) // Ж. вычисл. матем. и матем. физ. 1961. T.1. №3. С. 542-545.
Vestnik KRAUNC. Fiz.-Mat. Nauki. 2020. vol. 32. no. 3. pp. 15-28. ISSN 2079-6641
MATHEMATICS
MSC 65M20 Research Article
Numerical solution of a boundary value problem for a loaded parabolic equation with nonlocal boundary conditions
V. M. Abdullayev1'2
1 Azerbaijan State Oil and Industry University, Azadlig avenue 20, AZ1010, Baku, Azerbaijan
2 Institute of Control Systems, Azerbaijan National Academy of Sciences, 9 B. Vahabzade str. Baku, AZ1141, Azerbaijan
E-mail: vaqif_ab@rambler.ru
In the work, we propose a numerical method of solution to the boundary-value problem with respect to the loaded parabolic equation with nonlocal boundary conditions. We have obtained formulas and derived an algorithm for the solution of the problem. We provide the results of numerical solution to two test problems, which illustrates the efficiency of the approach proposed.
Key words: loaded differential equations, nonlocal conditions, method of lines, transfer of boundary conditions.
DOI: 10.26117/2079-6641-2020-32-3-15-28
Original article submitted: 05.08.2020 Revision submitted: 28.09.2020
For citation. Abdullayev V. M. Numerical solution of a boundary value problem for a loaded parabolic equation with nonlocal boundary conditions. Vestnik KRAUNC. Fiz.-mat. nauki. 2020,32: 3,15-28. DOI: 10.26117/2079-6641-2020-32-3-15-28
Competing interests. The authors declare that there are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. All authors contributed to this article. Authors are solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by all authors.
The content is published under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/deed.ru)
© Abdullayev V. M., 2020
Funding. The study was carried out without funding