CC BY
9
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
ЯО/ - 10.32743/UniTech.2022.103.10.14415
ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ КРОСС ДИФФУЗИИ С НЕЛОКАЛЬНЫМИ ГРАНИЧНЫМИ УСЛОВИЯМИ
Урунбаев Жасур Эркинович
PhD физ-мат. наук, Самаркандский государственный университет,
Республика Узбекистан, г. Самарканд E-mail: _ jasururunbayev@,gmail.com
NUMERICAL SOLUTION TO THE PROBLEM OF CROSS DIFFUSION WITH NONLOCAL BOUNDARY CONDITIONS
Jasur Urunbaev
PhD in Physics and Mathematics Sciences, Samarkand State University,
Republic of Uzbekistan, Samarkand
АННОТАЦИЯ
В работе исследуется асимптотика автомодельных решений нелинейной системы кросс диффузии, связанной с нелокальными граничными условиями. Получен главный член асимптотики автомодельных решений. Для численного исследования рассматриваемой задачи предложен способ выбора подходящего начального приближения для итерационного процесса. Используя асимптотические формулы в качестве начального приближения для итерационного процесса, произведены численные расчеты и анализ результатов.
ABSTRACT
Asymptotic behavior of self-similar solutions of a nonlinear cross-diffusion system coupled in the non-local boundary conditions is studied in the paper. The main term of the asymptotics of self-similar solutions is obtained. For a numerical study of the problem under consideration, a method is proposed for choosing an appropriate initial approximation for the iterative process. Using asymptotic formulas as an initial approximation for the iteration process, numerical calculations and analysis of the results are conducted.
Ключевые слова: кросс диффузия, критические экспоненты, глобального существование, неограниченные решения, автомодельный анализ.
Keywords: cross diffusion, critical exponents, global existence, unbounded solutions, self-similar analysis.
В настоящей работе исследуются качественные свойства решений нелинейной системы кросс диффузии, связанной с нелинейными граничными условиями
dt dx
du d ( щ _ du
dx I9
(1)
du d ( m _ dü I
— = —|um21— I, x e R , t > 0,
dt dx dx
_umi4 — (0, t) = uqi (0, t)
0, t) = uq2 (0, t), t > 0,
(2)
(x,0) = u0(x), ü(x,0) = u0(x), x e R+, (3)
где т1 > 1, д. > 0(г = 1,2), Щ и и0 (х) - неотрицательные непрерывные функции с компактным носителем в .
Процесс кросс диффузии (перекрестной диффузии) означает, что пространственное перемещение одного объекта, описываемого одной из переменных, происходит за счет диффузии другого объекта, описываемого другой переменной [6].
Модели кросс диффузии встречаются в различных областях естествознания. Например, в физических системах (физике плазмы) [1-3], в химических системах (динамике электролитических растворов), в биологических системах (кросс-диффузионном транспорте, динамике популяционных систем), в экологии (динамике возрастной структуры леса), в сейсмологии - модель Бурриджа-Кнопоффа, описывающая взаимодействие тектонических плит [4-7]. При исследовании биологической популяции и движения тектонических плит активно применяются математические модели с кросс-диффузией [4, 5].
—77 m2 _
1J
Библиографическое описание: Урунбаев Ж.Э. ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ КРОСС ДИФФУЗИИ С НЕЛОКАЛЬНЫМИ ГРАНИЧНЫМИ УСЛОВИЯМИ // Universum: технические науки : электрон. научн. журн. 2022. 10(103). URL: https://7universum.com/ru/tech/archive/item/14415
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
Известно, что системы вырождающихся уравнений могут не иметь классическое решение в области, где и, V = 0. В этом случае изучается обобщенное
решение системы (1), имеющей физический смысл в классе
u(х,t), и(х,t)> 0, 1. u^ UС (R+x( 0, +»))
Доказали, что min {y - r, y - Г } = 0, где
U
дх
и удовлетворяющие системе (1) в смысле распределения [1, 3].
Изучению условий глобальной разрешимости и неразрешимости задачи (1)-(3) при различных значениях числовых параметров посвящено большое количество работ [5-15] (подробно см. в библиографии [6]). Авторы работ [8, 9] изучали условия глобальной разрешимости и неразрешимости по времени решения и установили оценку решения вблизи времени взрыва нелокальной задачи диффузии
и = ихх, V =vxx, X>0,0<Т<0, (4) -и(0,г) = иСр, — 0,г) = uqvf}, 0<г<Т, (5) и(х,0) = и0(х), v(х,0) = х), х>0■ (6)
Доказали, что если pq < (1 — с)(1 — , то всякое решение задачи (4)-(6) является глобальным.
В работе [10] исследованы следующие задачи
и =(ип ) , ц=(ок) , х е , г > 0, (7)
—(ип) (0 г) = ^р (0, г),
(8)
(uk) (0,t) = uq (0,t), t > 0,
и (х,0 ) = и0 (х), v( x,0) = v( х), х е , (9)
Показано, что решение задачи (7)-(8) является глобальным, если pq <(п +1)(к +1)/4. Были получены условия для числовых параметров систем (7)-(9), при которых решение задачи взрывается за конечное время.
Также следует отметить работу [11], в которой исследовалась система (7) со следующими краевыми условиями
— (ип) (0, г) = uavp (0, г),
r =
2 p + k +1 - 2ß
4 pq-(k +1 - 2a)( n +1 - 2ß)
2 p + n +1 - 2ß 4 pq-(k +1 - 2a)( n +1 - 2ß)
1 - r (n -1) 1 - r2 (k -1)
У = —-—, у =-—
л 2 ^ 2
является критической экспонентой типа Фуджита. Данная работа посвящена исследованию асимптотики автомодельного решения задачи (1)-(3). Построены различные автомодельные решения задачи (1)-(3) для случая медленной диффузии (щ, щ > 1), являющиеся асимптотикой решений
рассматриваемой задачи. Для численного исследования предлагаются способы выбора подходящего начального приближения для итерационного процесса, сохраняющие качественные свойства задачи (1)-(3). Также сконструирован итерационный процесс и проведены численные расчеты, показывающие быструю сходимость к точному решению.
Система уравнений (1) при щ > 1 (/ = 1,2) описывает процессы с конечной скоростью распространения возмущений. Уравнения (1) при и (х, г),
Vх, г) = 0 являются вырождающимися, поэтому
задача (1)-(3) допускает обобщенное решение, не имеющее в точках вырождения необходимой гладкости.
Система (1) имеет ограниченные автомодельные решения с компактным носителем следующего вида
u ( ^ t ) = (Г +1 )-a f (£), (10)
u(x,t) = (T +1P p(£), £ = x(T +1)-ß
где T > 0
ß
q - m
q?- m
2q - m -1 2q2 - m -1
1
a
1
a? =
2q1 — щ — 1 2 2 ^ — щ — 1 функции , ф ( ^ ) ) являются решением следую-
щей задачи
Г2 =
а
-(uk) (0,t) = uquß(0,t), t > 0.
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
d
d
^Ш1-1 f] + ßgf + "f = о, nn
dg) dg (11)
f-1 dr)\+ß^ d?+^=о,
dg) dg
y»i-1 f (о ) = f' (0), r dg ' W (12)
f"2 -1 f( о ) = ^'2 (о)'
которая получается после подстановки (10) в (1)-(3) и некоторых упрощений. Рассмотрим следующие функции
1 i
ф(д) = (a -g2)m-1,
(13)
где а > 0.
Теорема 1. Пусть щ > 1 и щ > 1, тогда решение с компактным носителем системы уравнений (11) при £ ^ Vа имеет асимптотику
\m = cj{t;)(l + o{X)),
(14)
где
ß(ш2 -1)
ß(ш -1)
1 (ш -1)
V ( Ш2 -1)
Доказательство. Ищем решение системы уравнений (11) в следующем виде
(15)
где т = - 1п(а), щ (г), w2(т) - неотрицательные и ограниченные функции, т ^го при . После подстановки (15) в (11) получим следующую систему
- а (w1, ^)+
1
d_ dz
2(a-e~T\ m2 -1
т i \ ß
L1 (wp w2 ) + -
Ж
L2 ( wp Ж2 ) +
2(a-e ш -1
L2 ( Ж1, Ж2) +
V Ш2 -1
w.
- w\
где
L ( W1, W2 ) = Ж2Ш1 L2 (W2 ) = W1
-1
w.
Ш2 -1
^ - w1
ш -1
w„
ш -1
- w:
Отметим, что изучение решений последней системы уравнений является равносильным изучению тех решений системы уравнений (11), каждое из которых в некоторой промежутке [То, + , удовлетворяет неравенствам:
Ж (z)> о,-^- - ЖХФ о
m -1
"1 e
Т-\ w1 = о,
4 (a - e- 1
(16)
W:
w2 = о,
VV-, / ч
ш-1 1) 4 (a-e-)
w (z)> о,-^ - w;^ о ш -1
Покажем, что решения (щ (т),щ (т)) системы
уравнений (16) имеют конечные пределы при т ^ . Пусть
К (z)= L1 ( w1» w2 ) > h2 (Z)= L2 ( w1. w2 )■
Тогда систему (16) приведем к виду
(17)
c1 =
С2 =
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
h'T
2(а - e"T) m2 -1
h (т)-
ß
ж
1 1 V m 2 - 1
a e
2(а -e т) щ -1
h (т)-|
4 (а -е"т)
ж = 0,
ж
- ж
а2 е _ + —--г ж = 0,
m-1 1 J 4 (а-е-)
Для анализа решений последней системы уравнений рассмотрим следующие вспомогательные функции
G2 (Т,,2 ) =
2(а -е"т) m2 -1
2(а -е т) m -1
1 2 ,-ß
22
w
V m2 - 1
w.
m -1
1 ж
a e
ж
4 (а - e-)
a e
4 (а-e-)
тж, = 0,
ж = 0,
где К, К2 - действительные числа. Видно, что в соответствующей правой части последнего тождества функций О (т,к), 02 (т,к2) сохраняют знак,
т.е. удовлетворяют одно из неравенств
G1 (т,кх)> 0, G2(т,к2)> 0, G(т,)<0, G2(т,к2)<0,
(18)
в некотором промежутке
[V^^[т^ +да), ТТ ^ . Допустим, что для функций О (т,к), О2 (т,к2) пределы при т ^ +да не существуют. Тогда в силу колеблемости функций О (т, к), О2 (т, к2 ), прямая
О. = к (/ = 1,2) пересекает бесконечное число раз их графиками на интервале , +да) X [Г* , +00) ■
Но, на интервале I тк , +да)х| тк , +да) выполняется
одно из неравенств (18) и поэтому пересечение их графиков бесконечное число раз невозможно. Следовательно, графики функций О (т,к), О2 (т,к2)
пересекают прямую О = к. (1 = 1,2) на интервале
тк ,+да) X ,+°°) только один раз. Тогда, для
функций О (т,к), О (т,к) существует предел при т ^ +да . Следуя (17) для О (м, , О (м , имеем
G = жГ1"
G2 =
ж
( 0 \ m\ -1 ж1
V m 2 -1
П | = ( )
■ ж' 1 = 1 I —+ о(1),
Щ -1
ж
V m1 -1
(ж)
m -1 ж
Щ -1
'(1)
Поэтому необходимо, чтобы
Нш
т^+да
Нш
т^+да
2(а - e"T) m2 -1
,1
1 2
ж,
- ж!
V Щ2 - 1 У
«1 e ~т
1 -ж = 0,
4 (а - e-)
2(а - e~T) m1 -1
22
ж
- ж!
«2 e
2 ж = 0,
--ж--т-т
m1 -1 1У 4 (а - e-)
Отсюда, с учетом следующего предельного перехода
Нш
а - e
= 0,
ß(Щ2 -1) 2
V (m-1)
V 2 У
X .4X1/ (Щ2-1)
ß( Щ -1)
при т ^
г
ж =
e
п
ж =
легко убедиться в том, что
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
Теперь построим численную схему, основанную на метод конечных разностей. Для этого уравнения система (1) аппроксимируется со вторым порядком точности по пространственным координатам и с первым порядком по t. Построим итерационный процесс, во внутренних шагов итерации значения узлов вычисляются методом прогонки. Известно, что выбор подходящего начального приближения для итерационного процесса в общем случае является основной трудностью численного решения нелинейной задачи. При решении конкретных задач используются функции, отражающие некоторые свойства искомого решения и полученное на основе качественного анализа нелинейной задачи. Эта трудность, в зависимости от значения числовых параметров уравнений, преодолевается путем удачного выбора начального приближения, в качестве которого при вычислениях брались выше установленные асимптотические формулы. На основе выше приведенных результатов были произведены численые расчеты. Ниже приведем численные схемы и некоторые результаты вычислительных экспериментов.
Рассмотрим систему уравнений (1) с начальными данными (3) и краевыми условиями (2) и
u ( b, t ) = ф (t), u( b, t) = Ф2 (t),
Для удобства перепишем систему (1) следующим образом
du _ д dt dx du _ d dt dx
K (u) B ( u )
du dx
du
dx
(19)
m -1
где К(и) = ищ1-1, В ( и ) = и Теперь построим равномерную сетку Sh по х с шагом к :
= {х/ = /• А, А>О, / = 1,2,..п-Ь = Ь},
и временную сетку
Гт = у- г, г>0, у =1,2,...,«, ^-г^Т7, Г>0|.
Построим разностную схему. Для этого используем метод баланса и неявную разностную схему:
У
j+1
^ = £[a,+1 (g)(У^1 - yj+1) - а, (g)(yj+1 - y^)], i[С+1 (У)(j -g/+1)-С(У)(gj+1 -j)
g'+1 - g
(20)
/ = 2,3,...,я-1, / = 0,1,...,и-1,
ry,° = u ( x, ,о ),
g,° =u( x ,о ),
г = 0,1,...,и.
(21)
- a1 ( g ) = ( Уо)',
- С ( У ) ^^ = ( gi )'2,
/ = 0,1,...,«-1,
(22)
'yj =Ф (tj),
g«7 = Ф2 (tj),
/ = 2,3,...,«.
(23)
где а (^) , с () вычисляются по одной из следующих формул
а)
a (g ) = к
С, ( g ) = B
g, + g,-1 V 2
'у,- + У,--1X
24)
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
б)
g)= *()+ * (g.-i),
(У ) =
5 ( У.) + B (У.-)
Видно, что системы алгебраических уравнений (25) (20) нелинейны относительно у^+1 и ^+1. Для
численного решения таких систем нелинейных уравнений применимы различные итерационные методы. Используем для них метод простой итерации:
s+1 s
J+1 „J
УГ - yJ = 1 h2
т
s+1 s
gJ+1 - gJ _ 1
а
s+1 s+1 s+1 s+1
J+1 J+1
т
h2
7+1
'i+1
s g
V У
s У
V У
У/++1 - У.+1
s+1 s+1 g.+1 - g.
а
- c
g
V У
s У
V У
У. - У-
s+1 s+1 g.j+1 - gj+
(26)
где ^ = 0,1,2,....
Известно, что итерационные методы требуют подходящего начального приближения, обеспечивающего быструю сходимость к точному решению и сохраняющего физический смысл задач. При этом в качестве подходящих начальных приближений выбираются выше полученные асимптотические формулы.
Значения начальной итерации для каждого шага
^+1 я+1
по времени у , g. берутся из предыдущего шага по 0 0
времени: у ^+1 = у ^, gJ+1 = gJ . При счете по итерационной схеме задается точность итерации, при которой процесс продолжается до тех пор, пока не выполняются условия
max
0< i < n
max
0< i < n
s+1 s
У.-У.
s+1 s
g. s.
< £.
Введем обозначения у = у]. +1, ^ = у]. +1. Тогда разностные уравнения (26) можно записать в виде
An У,-1 - Си У + ВхгУг+1 _-Fh, (27)
s+1
s+1
s+1
A2У.-1 - С2. У. + B2У.+1 _- F2i,
где ЛП , А2г, ВП , ВИ , С1/, , Г , Г2/ учитывая формулы (25), определяются следующим образом:
л т
An _ 1Г а
g
V У
2h:
г/ s \m1 gJJ 1
V
V У
/ 5 N
gJJ
gi-1 V У
-Л
, ^_т Ci
У
V У
2h:
/ ^ \m2 У.+1
V У
да, -1
/ с \Щ2
у/-!'
V У
-1>
5 т
B1. _ тт+,
g
V У
2h:
g.+i
\m
g.
J+1
i- Л
vV У V У у
,^=1J '
У V У
2h
/ s -1
У+?
V У
/ ^ \m2 У/+1
V у
-
s т
С'._h
а I g 1 + а+i I g
+1, с,-
Для численного решения системы алгебраических уравнений (27) применяется метод прогонки. Согласно методу прогонки
У. _ai. (ßn + У.+1), g =a?i (ß + f.+1),
(28)
где с, с2/, Ру, р21 - коэффициенты, которые вычисляются по следующим формулам:
I У 1 + c+i I У
s s
+1, F1. = yJ+1, F 2i _ g/+'.
ß _ Alißß. + Fl.
ß1i+1 _ Г Д
a __Bi-_
a.+1 С -a A 9
ci. ai.Ai. _ B2i 1 _ С -a л '
О _ A?гß?г ß_
СИ «Л Ali
где / = 1,2,...л ■ Значения , а20, Д0, /?20 находятся из краевых условий (22).
2
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
октябрь, 2022 г.
Используя вышеизложенные численные схемы, проведен вычислительный эксперимент. Приведем некоторые результаты численных экспериментов. Шаг сетки достаточно мелкий ^=0.05, число узлов #=2500
.-3
и точность итерации задается в = 10 . Счет проводился до t=2 с шагом т = 0.02 . В качестве начального приближения для итерационного процесса брались формулы (10), (14).
Рисунок 1. Численное решение задачи (1)-(3) при д = 4.75, д = 5.5, щ = 1.15, щ = 1.35
Рисунок 2. Численное решение задачи (1)-(3) при д = 1.75, д = 1.95, щ = 1.55, щ = 1.65
Рисунок 3. Численное решение задачи (1)-(3) при д = 1.8, д = 2.3, щ = 1.75, щ = 1.25
№ 10 (103)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
На рис. 1 -3 представлены графики результатов численного решения задачи (1)-(3) при щ > 1 (/ = 1,2), соответствующей случаю медленной диффузии. При щ > 1 (/ = 1,2), как следует из
асимптотических формул (10), (14) и графиков, перемещение объекта происходит с конечной скоростью. Глубина проникновения диффузионной
октябрь, 2022 г.
волны зависит от времени и фронта волны (точка, в которой и (х, £), и(х, t) обращаются в нуль) для
каждой среды, находящейся в конечной точке:
хф =4а(Т + £.
Список литературы:
1. Z.Q. Wu, J.N. Zhao and J.X. Yin, et al., "Nonlinear Diffusion Equations," World Scientific Publishing Co., Inc., River Edge: 2001. - 500 p.
2. Арипов М.М. Методы эталонных уравнений для решения нелинейных краевых задач. Ташкент: Фан, 1988. -189 C.
3. Калашников А.С. Некоторые вопросы качественной теории нелинейных вырождающихся параболических уравнений второго порядка // УМН. - 1987 - т.42. 2 (254) - C. 135-176.
4. Murray J.D. Mathematical Biology, 3rd ed., Berlin: Springer, 2002-2003.
5. Malchow H, Petrovskii SV, Venturino E. Spatiotemporal patterns in ecology and epidemiology: theory, models, and simulations. London: Chapman & Hall/CRC Press; 2008.
6. M.A. Tsyganov, V.N. Biktashev, J. Brindley, A.V. Holden, G.R. Ivanitsky, Waves in cross-diffusion systems -a special class of nonlinear waves// UFN - 2007 - vol. 177, issue 3 - P. 275-300.
7. Levine, H., The role of critical exponents in blowup theorems// SIAM Rev. - 1990 - 32(2) - P. 262-288.
8. Wang S, Xie C H, Wang M X. Note on critical exponents for a system of heat equations coupled in the boundary conditions.// J Math Analysis Applic, - 1998 - 218 - P. 313-324.
9. Wang S, Xie C H, Wang M X, The blow-up rate for a system of heat equations completely coupled in the boundary conditions. // Nonlinear Anal - 1999 - 35 - P. 389-398.
10. Quiros F, Rossi J D. Blow-up set and Fujita-type curves for a degenerate parabolic system with nonlinear conditions. // Indiana Univ Math J - 2001 - 50 - P. 629-654.
11. Zheng S N, Song X F, Jiang Z X. Critical Fujita exponents for degenerate parabolic equations coupled via nonlinear boundary flux.// J Math Anal Appl - 2004 - 298 - P. 308-324.
12. Rakhmonov Z. On the properties of solutions of multidimensional nonlinear filtration problem with variable density and nonlocal boundary condition in the case of fast diffusion // Journal of Siberian Federal University. Mathematics & Physics 2016 - 9(2), 236-245.
13. Рахмонов З. Оценки решений нелинейной системы уравнений теплопроводности с переменной плотностью и с нелокальным граничным условием // Вестник НУУз - 2016 - №1(2) - C.145-154.
14. Aripov M.M., Matyakubov A.S. To the qualitative properties of solution of system equations not in divergence form of polytrophic filtration in variable density. //Nanosystems: Physics, Chemistry, Mathematics - 2017 - 8(3) - P. 317-322.
15. Aripov M.M., Matyakubov A.S. Self-similar solutions of a cross-diffusion parabolic system with variable density: explicit estimates and asymptotic behavior. //Nanosystems: Physics, Chemistry, Mathematics, 2017 - 8(1) - P. 5-12.
