УДК 517.955.8
DOI: 10.14529/mmph220103
АСИМПТОТИКА РЕШЕНИЯ ПЕРВОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ СИНГУЛЯРНО ВОЗМУЩЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ ВТОРОГО ПОРЯДКА ПАРАБОЛИЧЕСКОГО ТИПА
К.Г. Кожобеков, А.А. Шооруков, Д.А. Турсунов
Ошский государственный университет, г. Ош, Киргизская Республика E-mail: [email protected]
Аннотация. Строится полное равномерное асимптотическое разложение по малому параметру решения первой краевой задачи. Первая краевая задача ставится для сингулярно возмущенного линейного неоднородного дифференциального уравнения в частных производных второго порядка с двумя независимыми переменными параболического типа. Задача исследуется на прямоугольнике. Особенности задачи - присутствие малого параметра перед оператором теплопроводности, существование угловых пограничных слоев на нижних углах прямоугольника. Требуется построить равномерное асимптотическое разложение решения первой краевой задачи на прямоугольнике, с любой степенью точности, при стремлении малого параметра к нулю. Асимптотическое разложение решения по малому параметру строится методом Вишика-Люстерника. При решении поставленной задачи нами используются: методы интегрирования обыкновенных дифференциальных уравнений, классический метод малого параметра, метод пограничных функций Вишика-Люстерника и принцип максимума. Как обычно, задача решается в двух этапах: в первом этапе строится формальное разложение решения первой краевой задачи, а во втором этапе оценивается остаточный член полученного разложения и этим доказывается, что полученное разложение действительно является асимптотическим на всем прямоугольнике. В первом этапе формальное асимптотическое решение ищется в виде суммы шести функций (решений): внешнее решение, определенное на всем прямоугольнике, погранслойное решение в малой окрестности нижней стороны прямоугольника, два боковых погранслойных решения в малой окрестности боковых сторон прямоугольника и два угловых погранслойных решения в окрестностях нижних вершин прямоугольника. Все эти погранслойные решения экспоненциально убывают вне пограничных слоев.
Ключевые слова: асимптотическое решение; малый параметр; сингулярно возмущенная задача; первая краевая задача; уравнение теплопроводности; погранслойное решение.
Постановка задачи. В прямоугольнике рассмотрим первую краевую задачу [1, 2]:
( 2 Л
dz(t, x) 2 д z(t, x)
dt дх
2
■s2p(t, x) dz(t'x) + q(t, x)z(t, x) = f (t, x), (t, x) eQ, (1)
dx
z(0, x) = p( x), x e [0,1], (2)
z(t,0) = M(t), z(t,1) = m2(t), t e[0,T], (3)
где s - малый положительный параметр, 0< a - const, Q = {(t, x) 10 < t < T, 0 < x < 1} , q, p,f e C" (Q), cpeC" [0,1], eC" [0,T], q(t, x) > 0: (t, x) eQ, p(0) = M(0), p(1) = ^(0) •
Требуется построить равномерное асимптотическое приближение первой краевой задачи (1)-(3), при s —^ 0.
Решение задачи (1)-(3) будем искать в виде рядов по степеням малого параметра [3-5]:
z(t, x) = U (t, x) + V (t, x) + n1 (t, t) + n2 (t, t]2) + W1(t,t1) + W2(t,t2) (4)
где t = tIs, t = xIX, t2 = (1 -x)IX, X = *Js .
Подставляя соотношение (4) в задачу (1)-(3), получаем следующие задачи:
dU(t, x) _ 2 d2 U(t, x)
dt dx'
2
d U(t x)
-s2p(t, x)-+ q(t, x)U (t, x) = f (t, x), (t, x) eQ (5)
dx
dV(r, x) 2 d2V(r, x) 2 (r, x)
-~sa -+ £ p(r, x)-+ q(r,x)V(r, x) = 0, (r, x) eQj; (6)
8r dx dx
7д2и, (щ) 7-П, (щ) . -П, (гт) -ТЩ - Ч(Щ )П; (Щ ) = I2 7 ( Т +Я3р(г,Щ) -п^, (Щ ) е О-,; (7)
дщ дг дщ 7
дЩ (ТЩ) 2- Щ (т,Щ ) , л ,4 , ,дЩ (т,Щ) , , _
—Н-7— а -1 г + ^(Т'Щ)Щ (Т, Щ) = -^4,Щ) —т--, (Т, Щ) е О,, (8)
дт дщ- дщ
где Ц = {(т, х) | 0 <т< /и~1Т, 0 < х < 1} , = {(Щ )|0 < г < Т, 0 <щ <Я_1},
О,7 = {(т,щ) | 0 <т < Те~\ 0 < щ < Х~*}, . = 1,2.
Подставляя соотношение (4) в начальное условие (2), получаем:
ф) = и(0, х) + V(0, х) + П (0, щ ) + П- (0, т- ) + Щ (0, щ ) + Щ (0, щ- ) ^
V(0, х) = р(х) - и(0, х), (9)
Ж.(0,щ) = -П7(0,щ), . = 1,-. (10)
Теперь подставляя соотношение (4) в граничные условия (3), имеем:
м (0 = и (г, 0) + V (ге-,0) + П С, 0) + П2 (г, Л"1) + Щ (ге_1,0) + Щ (1е~х ,Л~1) , учитывая, что П2(г,2Г1) = 0, Щ (ге~1,2Г1) = 0 - условие для погранслойных функций, получаем:
Пх(г,0) = м(г) - и(г,0), (11)
Щ (т, 0) = -Г(т,0), (12)
аналогично
^ (г) = и (г, 1) + V (га-, 1) + П2 (г, 0) + Щ (ге-, 0) =>
П2(г,0) = ^(г) - и-(г,0), (13)
Щ-(т,0) = -V(т,1) . (14)
В результате мы получили шесть задач:
- из уравнения (5) методом малого параметра однозначно определяем и^,х);
- из (6) и (9) определяем У(т,х);
- из (7), (11) и (13) определяем П7 (г,щ ) ,у = 1, 2;
- из (8), (10), (12) и (14) определяем Щ (т,щ ), у = 1, 2. Начнем с уравнения (5).
Лемма 1. Для решения и (г, х) уравнения (5) справедливо формальное асимптотическое разложение
да
и (г, х) = ^екык (г, х), (15)
к=0
где Ык е Сда(О), к = 0,1,-,... - конкретизируются при доказательстве леммы 1.
Доказательство. Формально подставляя (15) в (5) и приравнивая коэффициенты при одинаковых степенях е, имеем:
Ы0(г, х) = , (16)
д(г, х)
, ч 1
uk (t, x) = —
f
q(t, x)
k
dUk-i(t, x) 2 d Uk-i(t, x) v dt dx2 y
p(t, x) dUk—2 (t, x) q(t, x) dx
k e N, U—1 (t, x) = 0, (17)
заметим, что % e Сю(Q), k = 0,1,2,.... Лемма 1 доказана.
Кожобеков К.Г., Шооруков А.А., Асимптотика решения первой краевой задачи
Турсунов Д.А. для сингулярно возмущенного дифференциального уравнения...
Перейдем к задаче (6), (9). Пусть
да
V(r, x) = ^skvk (г, x), (18)
k=0
где vk (г, x) - пока неизвестные функций. Подставляя (18) в (6) и (9), имеем:
да k=0
^X j, (x) + X г % (x)v.j (г, x)
у ОТ J=0 их J=0
= 0,
X skvk (0, x) = p( x) - X shuk (0, x), k=0 k=0
где qj (x) = аq(0, x)/аГ , q0 (x) = q(0, x), p, (x) = аJp(0, x)/аг. Отсюда следует, что
av0(r, x)
+ %0 (x)v0 (г, x) = 0, V0 (0, x) = cp(x) - U0 (0, x); (19)
%0( x)Vk (г, x) =
аг
aVk (г, x)
бТ (20)
= а2 + £¿РМ) ^^Х} + £х), ^(0,х) = _ик(0,х), к е N.
дх2 у=0 7 дх ]=1 7 7
Решение задачи (19) существует, единственно и представимо в виде:
у0(т,х) = (р(х) _и0(0,х))е^(х)т.
Для ^(т, х) имеем:
бл(г,х) 2 д2^(т,х)
—1 + q0 (х)у, (г, х) = а -^—- + тq1 (х)у0 (т, х), V, (0, х) = -и, (0, х).
дт бх2
Справедлива
Лемма 2. Решение задачи ^ + + + (г,х)еЦ, у(0,х) = у\х), хе[0,1]
( 2 тп+1 ^
р0 (х)г + рх (х)— +... + рп (х)--
2 п +1
у У
существует, единственно и представимо в виде у(т,х) = е-<?о(х№(х) + е-<?о(х)г
где д0(х)>0 *е[°Д]> еС°°[0Д].
Доказательство. Уравнение
^^ + %(х)у(г,х) = (^о(х) + А(х)г +... + ^(х)гй ), (Т,х) е Ц,
запишем в виде
/ . . \ г /
-И
(v(r,x>*°Wr)' =(ft(x) + A(x)r + ... + ^(x)r"),
полученное выражение интегрируем по т, учитывая начальное условие:
т
Ъ(т,х)е^х)т-у\х) = ¡(р0(х) + ... + рп{х)зп)сЬ => Кт,х) = е-^(х)ту\х) + Рп+11(т,х)е-<1о(х)т, о
Т2
где Рп Т(т, х) = р0 (х)т + р (х) — +... + рп (х) тп+ /п +1. Лемма2доказана.
С помощью леммы 2 доказывается существование и единственность решений задач (20). Кроме этого, из леммы 2 следует, что эти решения экспоненциально стремятся к нулю при стремлении тк бесконечности, т. е.:
ук(т,х) = О[ё~9°(х)т), т^да, д0(х)>0: хе[0,1], к = 0,1,2,...
Перейдем к задачам (7), (11) и (13). Здесь две задачи относительно функций П1(г,щ1) и П2(г,т2) - аналогичные, поэтому достаточно рассмотреть одну. Мы рассмотрим задачу относительно функций П1(г,щ1) :
-2 -q^^t,,) = ^ + Л4SM)
s,
Пусть
dt дх
П (t, 0) = fr (t) - U(t, 0), t e [0, T].
да
П l(t ,,1) = £ Л>и (t ,,1),
k=0
, (t,,) eQ
21, (21) (22)
(23)
где n1 k (t, ,1) - пока неизвестные функций, причем lim л1 k (t, ,1) = 0, t e [0, T].
Подставляя (23) в (21) и (22), имеем:
£л— k=0
f 2S2Klk(tt) ^ .
a -;Г2--£j (tKk-j= £Л
s,1 j=0
k+2
V
/
k=0
-t+ Л £0j(t) ax
£ ЛЛли(t, 0) = M(t) - £ Л2\(t, 0), t e [0,T],
к=0
к=0
djq(t,0) .. djp(t, 0) .. . m где q} (t) = ЧУ\ ', pj (t) = J, q0(t) = q(t,0) .
Отсюда для 7T10(t,^1) имеем:
2 S2^1,0
S, 12
- ^0 (t>1,0 = 0, (t, ,1) e Q21, >1,0 (t, 0) = A (t) - "0 (t, 0), lim TTlfi (t, ,1) = 0, t e [0, T].
Решение этой задачи можно представить в виде
>1,0 (t, ,1) = (м (t) - "0 (t, 0)) е
Для л1 k (t, ?/1), k e N имеем:
2 dV
s,1
k _ ^ ¿>1k, 2(t,,1) dn1k-2(t,,1) I— ^0(t) >1,k =£ ,/q (t) >1,k -., (t,1)+£ Pj(t) ——--+—k-,&:, )e Q2^
j=i j=o Sx a >12k (t, 0) = -"k(t,0), >12k-1 (t, 0) = 0, lim >1 k(t, ,1) = 0, t e [0,T].
Справедлива
Лемма 3. Решение задачи
2 а2ж .. -
a —- - q0(t)ж = е
s,2
2
(P1(t),1 +...+pn(t), (t,,1)
e Q
21'
>(t,0) = >0(t), lim >(t,,1) = 0, te [0,T],
существует, единственно и представимо в виде
>(t,,1) = е
Jq o(t) Jqo(t)
a >0(t) + е a
где q0(t)>0 t<=[0,T],q0,P,V еС°°[0,Г].
Лемма 3 доказывается прямым интегрированием, как и лемма 2.
(24)
a
a
a
Кожобеков К.Г., Шооруков А.А., Асимптотика решения первой краевой задачи
Турсунов Д.А. для сингулярно возмущенного дифференциального уравнения...
С помощью леммы 3 доказывается существование, единственность решений задач (24). Для
решений задач (24) справедливы оценки:
( -JüvщЛ --щ
e a
я(щ = O
, щ1 ^да, t е[0,Г], 0 < a, 0 < q0(t).
Перейдем к задачам (8), (10), (12) и (14) для определения угловых погранслойных функций WJ (т,^),] = 1, 2.
Здесь тоже достаточно рассмотреть одну из них, либо задачу для W1 (т,7), либо для W2(т,72), второе исследуется аналогично. Рассмотрим задачу для W1 (т,щ) :
ащ(г,щ) ,а2^(г,щ), , ^Л, Зз , .а^(г,щ)
а
■ - a
аЩ
%(г,щ )W (г,щ ) + Л р(г,;щ )-
ащ
= 0, (г,щ ) еа
з ,
Пусть
W (0,щ ) = -Щ (0,щ ), W (г,0) = —V (г, 0).
W (г,Щ ) = X^k^i,k (г,Щ),
(25)
(26)
(27)
k=0
где w1k (г,щ1 ) - пока неизвестные функций. Подставляя (27) в (25) и (26), получаем:
dw1 ■ а2 w1 j
' ■-a -+ q0wi,j = 0, (г,щ)еОз i, j = 0, 1;
аг
ащ2
2
где qn,m =
sw k 2 а 2w1 k
- - a —2,- + q0w1,k = ф (w1,k-1> w1,k - 2^.^ w'l,o,Г,Щl), (г,Щ1) еаз^ k = 2Д... аг ащ/
wu(0,Щ) = -яи(0,щ),w1,2k(г,0) = -vk(г,0), w>1,2k+1 (г,0) = 0, k = 0,1,..., аn+mq(0,0) Л
,0 < q0 = q(0,0), Фк - линейно зависят от предыдущих w1 j (j<k) и их
производных от щ1, полиномиально зависят от г, и щ . Если ввести обозначение w1 к (г, щ1) = e"^yk (г, щ1) , то
аi,k(r,?7i )
5г
= -q0e q0 yk (г,Щ1)+e
-%0г
3k, щ ) . аХ* ^ щ ) = Ö2 Л k Щ)
аг
ащ2
ащ2
и рассматриваемая задача примет вид:
dyk _2 а2y*
--a
аг ащ2
yk(0,щ) = -Ли(0,щ),y2k(г,0) = -eq0\(г,0), y2k+l(г,0) = 0, k = 0,1,... где Ф, = 0, j = 0, 1 .
Решения задач (28), (29) существуют, единственны и представимы в виде [6]:
да г г да
yk (г, щ ) = - j яи (0, £)0(г, £ щ d + j Wk (t)H(щ1, г -1)dt + j j Фk (t, x^1, x, г - t)dxdt
= Фk (y1,k-1, y1,k - 2 Л,0,г,щ1), (г,щ1) е аз 1 , k = 0 l...
(28) (29)
0 0
где
H (Ц\,г) =
щ1
щ
з/2 e 4г, а(г,x,Vl) = 2a V яг
2
(У1 - x)
exp
4a2 г
exp
f (щ + x)2"
4a2г
^2к(т) = _е vk(т,0),^2к-1(т) = 0, фк(т,7) = фк(Дк-1,Дк-2,...,>'1,0,т,^1). Отсюда следует, что функции к (т, ) экспоненциально убывают при т + 7, ^ ж .
0
0
1
Таким образом, нами построены все функции, входящие в правую часть равенства (4).
Обоснование. Оценим остаточный член разложения (4).
Пусть z(t, х) = zn (г, х) + R(г, х) , где ^г, х) - остаточный член разложения,
Zn (г, х) = ип (г, х) + Vn (т, х) + П1 ,2п+1 (г,т) + П- ,2п+1 (гЩ + Щ ,2п+1 ,2п+1
п п 2 п+1
2Ч,к (
Un (t, x) = £^k"k (t, x),Fn (r, x) = £skVk (r, x), П j,2n+1 (t,, j ) = £ Лклик (t,,j ),
к=0 к=0 к=0
2 п+1
ЩЛ2п+1(т,77) = £ (т,щ).
к=0
Тогда для остаточного члена получим следующую задачу:
£
fSR(t, x) _а2 S2R(t, x)^
at Sx
2
£2p(t,x)SR(^x) + q(t,x)R(t,x) = 0(£n+1), (t,x) eQ, (31) dx
R(0, x) = O(sn+1), x e [0,1], z(t, 0) = z(t, 1) = O(£n1 ), t e [0,T], £ — 0. (32)
Применяя принцип максимума [7], получаем:
|R(t,x) < max iq"1(t,x)O(£n+1), O(sn+1)).
1 1 f(rV=ol >
Отсюда имеем:
(t,x)eQ 0<£<<1
R(t,x) = O(£n+1), (t,x) eQ, £ — 0.
Теорема. Для решения задачи (1)-(3) при стремлении малого параметра к нулю в области О справедливо асимптотическое разложение
*(г, х) = ип (г, х) + Vn (т, х) + П1,2п+1(гт) + П2Лп+х(гщ2) + Щ-п+^щ) + Щ2Лп+х(Т + О(еп+1), где функции ип (г, х) = £п=0 еик (г, х\ К (Т, х) = £п=0 еЧ (Т, x), П7,2п+1(г, щ ) = £-=+1 ,к (г, щ),
^7,2п+1(т,77 ) = £- 1^к-7,к (Т,Л] ) определены выше.
Заключение. Нами построено полное равномерное асимптотическое разложение по малому параметру решения первой краевой задачи для сингулярно возмущенного линейного неоднородного дифференциального уравнения в частных производных второго порядка с двумя независимыми переменными параболического типа. Доказано, что полученное разложение действительно является асимптотическим решением поставленной задачи на всем прямоугольнике. Данная работа для нас является началом исследования бисингулярно возмущенных задач параболического типа, в следующих работах, ссылаясь на эту статью, мы будем исследовать только бисингуляр-ные случаи.
Литература
1. Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. -М.: Наука, 1966. - 724 с.
2. Zauderer, E. Partial Differential Equations of Applied Mathematics / E. Zauderer. - New York etc.: John Wiley & Sons, Inc. - 891 p.
3. Алымкулов, К. Об одном методе построения асимптотических разложений бисингулярно возмущенных задач / К. Алымкулов, Д.А. Турсунов // Известия вузов. Математика. - 2016. -№ 12. - С. 3-11.
4. Вишик, М.И. Регулярное вырождение и пограничный слой для линейных дифференциальных уравнений с малым параметром / М.И. Вишик, Л.А. Люстерник // УМН. - Т. 12, вып. 5(77). -С.3-122.
5. Треногин, В.А. Развитие и приложение асимптотического метода Люстерника-Вишика / В.А. Треногин // УМН. - 1970, Т. 25, вып. 4(154). - С. 123-156.
6. Полянин А.Д. Справочник по линейным уравнениям математической физики. - М.: Физмат. лит., 2001. - 575 с.
Koжобеков K.r., Шооруков A.A., Acимптoтика решения первой краевой задачи
Турсуновß.A. для сингулярно возмущенного дифференциального уравнения...
7. Protter, M.H. Maximum Principles in Differential Equations / M.H. Protter, H.F. Weinberger. -
Prentice Hall, Englewood Cliffs, New Jersey, 1967, 261 p.
Поступила в редакцию 27 декабря 2021 г.
Сведения об авторах
Кожобеков Кудайберди Гапаралиевич - доктор физико-математических наук, ректор Ошско-го государственного университета, г. Ош, Кыргызская Республика, ORCID iD: https://orcid.org/0000-0003-0S56-5113
Шооруков Асылбек Абдибахапович - аспирант, Ошский государственный университет, г. Ош, Кыргызская Республика, ORCID iD: https://orcid.org/ 0000-0002-9550-090X
Турсунов Дилмурат Абдиллажанович - доктор физико-математических наук, директор ВШМОП, Ошский государственный университет, г. Ош, Кыргызская Республика, е-mail: [email protected], ORCID iD: https://orcid.org/0000-0002-6990-1742
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2022, vol. 14, no. 1, pp. 27-34
DOI: 10.14529/mmph220103
ASYMPTOTICS OF THE SOLUTION OF THE FIRST BOUNDARY VALUE PROBLEM FOR A SINGULARLY PERTURBED DIFFERENTIAL EQUATION IN PARTIAL DERIVATIVES OF THE SECOND ORDER OF PARABOLIC TYPE
K.G. Kozhobekov, A.A. Shoorukov, D.A. Tursunov
Osh State University, Osh, Kyrgyz Republic E-mail: [email protected]
Abstract. The article constructs a complete uniform asymptotic expansion in a small parameter of the solution of the first boundary value problem. The first boundary value problem is posed for a singularly perturbed linear inhomogeneous second-order partial differential equation with two independent variables of parabolic type. The problem is investigated on a rectangle. The peculiarities of the problem are the presence of a small parameter in front of the heat conduction operator, the existence of corner boundary layers at the lower corners of the rectangle. It is required to construct a uniform asymptotic expansion of the solution of the first boundary value problem on a rectangle, with any degree of accuracy, as the small parameter tends to zero. The asymptotic expansion of the solution in terms of a small parameter is constructed by the Vishik-Lyusternik method. When solving the problem, we use: methods of integration of ordinary differential equations, the classical method of a small parameter, the Vishik-Lyusternik boundary function method, and the maximum principle. As usual, the problem is solved in two stages: in the first stage, a formal expansion of the solution of the first boundary value problem is constructed; and in the second stage, the remainder of the resulting expansion is estimated and this proves that the resulting expansion is indeed asymptotic over the entire rectangle. In the first stage, a formal asymptotic solution is sought in the form of a sum of six functions (solutions): an external solution defined on the entire rectangle; boundary layer solution in a small neighborhood of the lower side of the rectangle; two lateral boundary layer solutions in a small neighborhood of the lateral sides of the rectangle and two corner boundary layer solutions in the neighborhood of the lower vertices of the rectangle. All these boundary layer solutions exponentially decrease outside the boundary layers.
Keywords: asymptotic solution; small parameter; singularly perturbed problem; first boundary value problem; heat equation; boundary layer solution.
References
1. Tikhonov A.N., Samarskiy A.A. Uravneniya matematicheskoy fiziki (Equations of mathematical physics). Moscow, Nauka Publ., 1966, 724 p. (in Russ.).
2. Zauderer E. Partial Differential Equations of Applied Mathematics. New York etc.: John Wiley & Sons, Inc., 891 p.
3. Alymkulov K., Tursunov D.A. On a method of construction of asymptotic decompositions of bisingular perturbed problems. Russian Mathematics, 2016, Vol. 60, no. 12, pp. 1-8. DOI: 10.3103/S1066369X1612001X
4. Vishik M.I., Lyusternik L.A. Regular degeneration and boundary layer for linear differential equations with small parameter. UspekhiMat. Nauk, 1957, Vol. 12, Iss. 5(77), pp. 3-122. (in Russ.).
5. Trenogin V.A. The development and applications of the asymptotic method of Lyusternik and Vishik. Russian Mathematical Surveys, 1970, Vol. 25, Iss. 4, pp. 119-156. DOI: 10.1070/RM1970v025n04ABEH001262
6. Polyanin A.D. Spravochnik po lineynym uravneniyam matematicheskoy fiziki (A Handbook of Linear Equations in Mathematical Physics), Moscow, Fiz.-mat. lit. Publ., 2001, 575 p. (in Russ.).
7. Protter M.H., Weinberger H.F. Maximum Principles in Differential Equations. Prentice Hall, Englewood Cliffs, New Jersey, 1967, 261 p. DOI:10.1007/978-1-4612-5282-5
Received December 27, 2021
Information about the authors
Kozhobekov Kudaiberdi Gaparalievich is Dr. Sc. (Physics and Mathematics), Rector of Osh State University, Osh, Kyrgyz Republic, ORCID iD: https://orcid.org/0000-0003-0856-5113
Shoorukov Asylbek Abdibahapovich is Post-graduate Student, Osh State University, Osh, Kyrgyz Republic, ORCID iD: https://orcid.org/0000-0002-9550-090X
Tursunov Dilmurat Abdillazhanovich is Dr. Sc. (Physics and Mathematics), Director of the HSIEP, Osh State University, Osh, Kyrgyz Republic, e-mail: [email protected], ORCID iD: https://orcid.org/0000-0002-6990-1742