Vestnik KRAUNC. Fiz.-Mat. nauki. 2023. vol. 43. no. 2. P. 9-19. ISSN 2079-6641
MATHEMATICS
" https://doi.org/10.26117/2079-6641-2023-43-2-9-19
Research Article
Full text in English
MSC 35B44, 35C06, 35K51, 35K61
Global and Blow-Up Solutions for a Nonlinear Diffusion System with a Source and Nonlinear Boundary Conditions
A. A. Alimov1'2*, Z.R. Rakhmonov1
1 National University of Uzbekistan named after Mirzo Ulugbek, 100174, Tashkent, Universitetskaya, 4, Uzbekistan
2 Tashkent branch of the G.V. Plekhanov Russian University of Economics, 100164, Tashkent, Shakhriobod, 3, Uzbekistan
Abstract. In this paper, we study the global solvability and unsolvability of a nonlinear diffusion system with nonlinear boundary conditions in the case of slow diffusion. The conditions for the global existence of the solution in time and the unsolvability of the solution of the diffusion problem in a homogeneous medium are found on the basis of comparison principle and self-similar analysis. We obtain the critical exponent of the Fujita type and the critical global existence exponent, which plays an important role in the study of the qualitative properties of nonlinear models of reaction-diffusion, heat transfer, filtration and other physical, chemical, biological processes. In the global solvability case the principal terms of the asymptotic of solutions are obtained. It is well known that iterative methods require the presence of a suitable initial approximation, resulting in a rapid convergence to the exact solution and preserving qualitative properties of nonlinear processes under study, it is a major challenge for the numerical solution of nonlinear problems. This difficulty, depending on the value of the numerical parameters of the equation is overcome by a successful choice of initial approximations, which are mainly in the calculations suggested taking asymptotic formula.
Key words: blow-up, nonlinear boundary condition, critical exponents, nonlinear diffusion system, asymptotic
Received: 21.05.2023; Revised: 09.06.2023; Accepted: 30.06.2023; First online: 03.07.2023
For citation. Alimov A. A., Rakhmonov Z.R. Global and blow-up solutions for a nonlinear diffusion system with a source and nonlinear boundary conditions. Vestnik KRAUNC. Fiz.-mat. nauki. 2023,43: 2,9-19. EDN: XJQODE. https://doi.org/10.26117/2079-6641-2023-43-2-9-19. Funding. The work was carried out without the support of funds.
Competing interests. There are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. Авторы участвовали в написании статьи и полностью несут ответственность за предоставление окончательной версии статьи в печать.
* Correspondence: A E-mail: akram.alimov18@gmail.com ^
The content is published under the terms of the Creative Commons Attribution 4-0 International License © Alimov A. A., Rakhmonov Z. R., 2023
© Institute of Cosmophysical Research and Radio Wave Propagation, 2023 (original layout, design, compilation)
Вестник КРАУНЦ. Физ.-мат. науки. 2023. Т. 43. №2. C.9-19. ISSN 2079-6641
МАТЕМАТИКА
" https://doi.org/10.26117/2079-6641-2023-43-2-9-19
Научная статья
Полный текст на английском языке
УДК 517.957
Глобальные решения и решения с обострением для
и и 11 о
нелинейной диффузионном системы с источником и нелинейными граничными условиями
А. А. Алимов1'2*, З. Р. Рахмонов1
1 Национальный университет Узбекистана имени Мирзо Улугбека, 100174, г. Ташкент, ул. Университетская, 4, Узбекистан
2 Ташкентский филиал Института им. Плеханова Российский экономический университет, 100164, г. Ташкент, ул. Шахриобод, 3, Узбекистан
Аннотация. В данной работе изучается глобальная разрешимость и неразрешимость одной нелинейной системы диффузии с нелинейными граничными условиями в случае медленной диффузии. Найдены условия глобального существования решения по времени и неразрешимости решения нелинейной задачи диффузии в однородной среде на основе автомодельного анализа и метода сравнения решений. Получены критическая экспонента типа Фуджита, и критическая экспонента глобального существования решения по времени, играющих важную роль при исследованиях качественных свойств нелинейных моделей реакции - диффузии, теплопроводности, фильтрации и других физических, химических, биологических процессов. В случае глобальной разрешимости получен главный член асимптотики автомодельных решений. Известно, что итерационные методы требуют наличия подходящего начального приближения, приводящее быстрой сходимости к точному решению и сохраняющие качественные свойства изучаемых нелинейных процессов, это является основной трудностью для численного решения нелинейных задач. Эта трудность в зависимости от значения числовых параметров нелинейной системы диффузии с нелинейными граничными условиями преодолевается путем удачного выбора начальных приближений, в качестве которых при вычислениях предложено брать полученные асимптотические формулы.
Ключевые слова: обострение, нелинейное краевое условие, критические показатели, нелинейная диффузионная система, асимптотика
Получение: 21.05.2023; Исправление: 09.06.2023; Принятие: 30.06.2023; Публикация онлайн: 03.07.2023
Для цитирования. Alimov A.A., Rakhmonov Z.R. Global and Blow-Up Solutions for a Nonlinear Diffusion System with a Source and Nonlinear Boundary Conditions // Вестник КРАУНЦ. Физ.-мат. науки. 2023. Т. 43. № 2. C. 9-19. EDN: XJQODE. https://doi.org/10.26117/2079-6641-2023-43-2-9-19. Финансирование. Работа выполнялась без поддержки фондов.
Конкурирующие интересы. Конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответственность. Авторы участвовали в написании статьи и полностью несут
ответственность за предоставление окончательной версии статьи в печать.
* Корреспонденция: E-mail: akram.alimov18@gmail.com ф
Контент публикуется на условиях Creative Commons Attribution 4.0 International License © Alimov A. A., Rakhmonov Z. R., 2023
© ИКИР ДВО РАН, 2023 (оригинал-макет, дизайн, составление)
Introduction
In this article, it is dealt with the doubly degenerate parabolic equations with the source
3ut 9 9t dx
9u
dx
m—1
suk, ?
—- I + u?1, x g R+, t>0, ox ' L
i = 1,2
(1)
coupled through nonlinear boundary flux
duk m— 1 9uk
dx dx x=0
= uq—i(0,t), t>0, i = 1,2
(2)
where m > 1, k > 1 and qi, pi > 0 are numerical parameters. The following initial data should be taken into account
uilt=o = uio(x), i = 1,2
(3)
moreover, they are expected to be continuous, non-negative and compact in R+.
Parabolic equations with nonlinearity (1) are found in population dynamics, heat transfer, chemical reactions and so forth. The functions ui (t,x), u2(t,x) serve as two populations' densities in terms of biology while a migration progresses or the thickness of two types of chemical reagents within a chemical process and two different types of materials' temperatures during propagation. Most of all, equations (1) can be used to characterize unsteady flows in a liquid medium with a power-law dependence of shear stress on displacement velocity under polytropic conditions.
It has already established that the local existence of weak solutions for the problem (1)-(3) is determined by the usual integration method and can be easily established as well as the comparison principle (see [9] - [13] and [1] - [5]).
Global existence and blow-up conditions of nonlinear parabolic systems are intensively studied (see [1]- [3], [5,11] and references therein). In particular, there is a great interest in critical Fujita exponents for various non-linear parabolic equations in mathematical physics (see [3,10-12] and references therein).
In [5], investigated following problem
ut = (uk
x > 0, 0 < t < T,
—(uk)x(0,t)= uq(0,t), 0<t<T,
(4)
x>0
u(x,0) = u0(x),
V
and the heat conduction problem
[ut = (|ux|k—1 ux)x, x > 0, 0 < t < T, <— |ux|k—1Vx(0,t) =uq (0,t), 0<t<T, u(x,0) = u0(x), x>0
xx
with k > 1, q > 0 and u0 has compact support. It has been proved that, for problem (4),
1
fl/ i 1 ^ titiiqvq r\C
Fujita type in contrast with (5), the critical exponent of Fujita type qc = 2k and the
2k
critical exponent of global solvability is q0 =
the critical global exponent is q0 = 2(k +1) where qc = k +1 is the critical exponent of
k +1
In [6] Z.R.Rakhmonov and A.I.Tillaev studied the problem
ip(x)ut = (|ux|k-2ux)x + p(x)uß, (x,t) g R+ X (0, +oo) < -|ux|k-2ux(0,t)= um(0,t), t > 0 u(x,0)= u0(x) > 0, x g R+
(6)
with k > 2, ß, m > 0, p(x) = x n, n g R, u0(x) - is a bounded, continuous, nonnegative and nontrivial initial data. They established that:
(2 - n)(k - 1)
-if 0 < ß < 1 and 0 < m <---the problem has global solution;
k- n
(2 - n)(k -1)
-if ß < 1 and m >-the problem has blow-up solution.
k-n
In [9] Zhaoyin Xiang, Chunlai Mu and Yulan Wang studied the next problem
/ 9u 8 ( 9umi pi 2 9umi
9t 9x I 9x 9x
< 9v 9vm2 P2- ~2 9vm2 ^
9t 9^ 9x 9x J
(x,t) g R+ X (0,T)
(7)
dumi pi -2 9umi
9x 9x
9vm2 P2 -2 9vm2
9x dx
= vqi(0,t)
x=0
= uq2(0,t)
x=0
t g (0,T)
(8)
u(x,0) = U0 (x)
, x G R+ (9)
kv(x,0) = V0(x)
where m^ > 1, pi > 2, qi > 0, i = 1,2. They discovered that:
(i) if q1 q2 < ((p1 — 1 )(p2 — 1 )(m1 +1 )(m2 +1 ))/p1p2 then every nonnegative solution of the problem (7)-(9) is global in time;
(ii) if q1 q2 > ((p1 — 1)(p2 — 1)(m1 + 1)(m2 + 1))/p1p2 then the problem (7)-(9) solutions that explode in a finite amount of time exist.
When q1q2 > ((p1 — 1)(p2 — 1)(m1 + 1)(m2 + 1))/p1P2.
(i) if min(a1 + ^1,^2 + ^2} > 0, then solution of the problem (7)-(9) is global in time;
(ii) if max{a/| + ^1,^2 + ^2} < 0, then the solution of the problem (7)-(9) is blow-up.
In [7, 8] nonlinear parabolic equations in approximate derivatives with nonlinear boundary conditions are used to describe mathematical models of nonlinear cross diffusion. There are very few explicit analytical solutions to these nonlinearly coupled partial differential equation systems, so several numerical methods have been used to obtain approximations. Self-similar analysis and the standard equation approach are used to study the qualitative characteristics considering a cross-diffusion system that is nonlinear and has nonlocal boundary conditions. For the slow diffusion scenario, a number of solutions of the cross-diffusion problem have been developed that are self-similar. It is shown that a nonlinear cross-diffusion system of parabolic equations, connected by nonlinear boundary conditions may not have global solutions for some values of numerical parameters. The Fujita-type critical exponent and the global solvability critical exponent are established using self-similar analysis and the comparison principle. Upper limits are achieved for global solutions while lower limits are found for blow-up solutions using the comparison theorem [1].
This article was inspired by the aforementioned works and has two goals. In order to determine the (1)-(3) system's critical global existence curve, firstly one should to build the self-similar super- and subsolutions. On the other hand, with the help of some fresh findings, the critical curve of the Fujita type is hypothesized. Being faced with a system as opposed to a single equation compels us to create some novel strategies.
Degenerate equations don't necessarily need to have classical solutions, as is common knowledge. As a result, its solution is comprehended in a broad way, see [1,2].
Definition 1. The function u(x,t) is seen as the inadequate remedy for the problem
m-19uk
(1) - (3) in a = {(0, +oo) x (0,t)}, if o < ui(x,t) e c(a),
9x
9x
g c (a), i = 1,2
and if it complies with (1)-(3) in the sense of distribution in a, where T> 0 is the longest possible time, see [6].
Main results
This paragraph is intended to build self-similar sub- and super-solutions to (1)-(3), to establish the theorems of Global existence and nonexistence solutions. The first theorem is about the conditions in which the problem (1)-(3) has the global solution. The global existence of solution draws conclusions from the comparison principle.
( — V
Theorem 1. If qiq2 <1 —I (k +1 — pi ) (k +1 — P2), then every nonnegative
solution of the problem (1)-(3) is global in time.
Proof. To prove the theorem self-similar sub-solution has been constructed and showed that it is limited for any t > 0. For this purpose, it has been looked for strict supersolutions of self-similar form
Ui(t,x) = eh2i-1 (M + e—Lixe-h2it ) k (10)
where L, > 0, H2i-i,2i > 0, M = max|||ul+1|; i = 1,2. Substituting (10) into (1)-(2) and using comparison principles it has been obtained:
^ = h2i-i ■ eh2i-i t. (m + e-Lixe-h2it) k + e(h2i-i -H2iH. k. U ■ x ■ H2i x x (M + e-Lixe-h2it) k-1 > H2i-1 eh2i-it (M + e-Lixe-h2it) k > H2i-1 e^i-^Mk
9
9x
9uk
9x
m-1
9u
9x
- I = mLm+1g[h2i-lkm-(m+1)h2i]te-mLixe h2it < mLm+1 e[h2i-i km-(m+1 )h2i]t
upi = epih2i-it ^m + e-Lixe h2i^ k < epih2i-it (M +1)
Pi
duk m- 1 duk
9x 9x x=0
= Lme(h2i-ik-h2i)mt, uqi , | = eqih5-2it (M + 1)"k
The solution u- are considered to be global if inequalities
9U, 9 dt > 9x
9uk
9x
m-1
9uk , p.
cU" I +Upi, i = 1,2
(11)
hold for any x e R+, t > 0. That is why using the computations above in (11) following expressions have been achieved:
H2i-1 eh2i-itMk > mLm+1e[h2i-ikm-(m+1)h2i]t + epih2i-it (m +1)£
5i
Lme(h2i-ik-h2i)mt = eqih5-2it (m + 1) Ik
Li = (M +1)km , qih5-2i = (h2i-1k-%) m, i = 1,2
q1 q2
h2i-1 > h2i-1 km - (m +1 )h2i + piH2i-1, H k - H2 = —H3, H3k - H4 = —H1
mm
h2i >
(km+pi -1) h2i-1 m +1 :
qi (km+p1 -1) Hi qi (km + p2 -1) H3
h1k--H3 >-—-,H3k--H1 >-—-
m m +1 m m +1
computing the last expressions, it can be seen that the last inequality should be always
hold for any m > 1, k > 1 in order to the solution of the problem (1)-(3) was global in
time. Theorem is proved.
2
Remark 1. Theorem demonstrates that q1 q2 = critical global existence of the problem (1)-(3).
m
m +1
(k +1 - p1)(k +1 - P2) is
k
m(p3-i — 1 )(pi + k)
Theorem 2. If 0 < pi < 1 and qi > —--—-—— or pi > 1 and qi <
(Pi - 1)(m +1) then, each of the solutions to (1)-(3) blows up.
m(p3-i - 1)(Pi + k) (Pi - 1)(m +1)
Proof. To prove the theorem the subsolutions of the problem (1)-(3) have been looked for in the next form:
Ui(t,x)= tai fi(£i), Li = xt—ßi,
(12)
1 pi - km where ai = --, pi = --—-—, i = 1,2.
1 -pi (pi -1 )(m +1)
After substitution (12) into (1)-(3) it has been reached the next self-similar inequalities and boundary conditions that should be hold for any ui(t,x) that treated as blow up solutions:
d
dLi
dfk
dLi
m—1
dfk\ df _
i 1 + ßiLiTT" — aifi + fP1 > 0
dLi
dLi
(13)
9uk m— 1 9uk
9x 9x x=0
< u3—i(0,t)
(14)
Let
m+1 m+1 \ mk-1
fi(^i) = AiU m — L m
(15)
Substitution (15) into (13), (14) lead us to the following conditions that show (14) always takes place:
/ k(m + 1)N
I mk — 1
m +1 \ mk m +1
Amk > ßi-m^A^
mk — 1
mk — 1
Ai >
ßi
/ mk — 1
1 k(m +1)
i
mk-1
fP1 = ApM Oi m — ^
m+1 m+1 \ mk-W m+1 m+1 \ mk-1
+
o m — L i ^
-(Pi—1)
<
+
m+1 m+1 \ mk-1
< AP1 oi mk-1 oi m — L m
+
(m+1)(pi-1)
AP1 a mk-1
> «iAi + Am^k(m +1)1
mk — 1
m
m
m
m
m
m
m
m
By taking
(m+IHpi-1)
a mk-1 > aiAi-Pi + Amk—pi I
. /k(m + 1)N
mk — 1
0<pi < 1 and qt >
m(p3—i — 1 )(pi + k)
can be easily checked and ensure that A1 and A2
(Pi - 1)(m +1)
can be taken sufficient to prevent inequalities (13) and (14) are valid. Because of this, if the initial data u1 (x,0), u2(x,0) are large enough that u10(x) > u1 (x,0), u20(x) > u2(x,0) then ui(t,x), i = 1,2 is a subsolution to (1)-(3). According to the comparison principle, for enormous beginning data, the solutions of (1)-(3) blow up in a finite amount of time. The proof is finished.
A 1
/ m(k +1) Theorem 3. If q1 q2 < -
and pi > 1 H— m + —, then every solution ' k / k
m +1
of the problem (1)-(3) is blow-up in finite time.
Proof. It is efficiently enough to show that the problem (1)-(3) without sources fulfilled out conditions. Let, construct
uib(t,x)= t*- gi(^i), £,i = xt-Yi (16)
where gi are two compactly supported functions,
Mi =
m[m(k + 1) + (m + 1)qt]
(m(k +1 ))2 - (m +1 )2qiq3-?Ti (m(k +1 ))2 - (m +1 )2qiqs-i
Now substituting (16) into (1)-(3) and obtain the next:
m[mk(k +1 ) + (mk — 1 )qt] — (m +1 )q1 q2
d
dgk
d^i
m—1
dgk \ dgi d^ ) + — Migi > 0
(17)
dgk m— 1 dgk
d^i d^i
< g3q-i(0)
Now it is time to find self-similar solutions of the problem (17), (18). Let
then
gi(£i) = Bi(bi — £i) mk-1 dgi Bim
d^i mk — 1
(bi — £,i)mk-1
(18)
(19)
dgi
Yi^i^— Migi =
d^i
Bi^ , m i m
-£,i (bi — ^i) mm-T—1 — MiBi (bi — mmT =
mk — 1
B i ^^ mi mi
(bi — £i)mm-T—1 — MiBi (bi — ^i)mm-T—1 (bi — ^i) >
mk — 1
—
biBim mk — 1
MibiBi (bi — £,i)
mk-1
+
m
m
d
dti
dgk
dti
m—1
dti
= Bi
mk
m
mk — 1
m+1
km (bi — Li)
--
m iU C lmk-1 ' —
- biBd Mi +
m
mk - 1
(bi—Li)mk-1
Bmk—1 > Wmkz1
i - km I m
m+1
Mi +
1 m ^ dgk m-1 k d9k
' mk — 1 1 ' dti dti
/ ti
< 93-i(0)
Applying comparison principles to the expressions above it is obtained:
Bf (bi - torn
dgk m- 1 dgk
dti dti
Bk (bi—torn-1
i—
m-1
i—
ti=0
*mk (bi—Li)+mm-1 1)m
ti=0
ti=0
m qim
= Bmkbimk-1 < Bqiib3T—
, M 1
And this show that when pt > ( 1 + — m + —, (17) and (18) are valid. It results from
v k
the comparison concept that (1)-(3) have solutions blowing up in a finite time.
Theorem 4. If q 1 q2 < (m(k +1)) and pt > 1, then every solution of the problem (1)-(3) is blow-up in finite time.
Proof. Theorems can be proved in the same manner as it was done in [8,14].
Conclusion
It is accomplished to acquire the diffusive system's solution of the type Zeldovich-Barenblatt. It is demonstrated that the nonlinear diffusion issue characterized that global solutions may not exist for degenerate parabolic systems coupled via nonlinear boundary conditions for specific values of numerical parameters. Using the comparison approach, it is possible to study the finite speed properties of the problem diffusion with a source.
References
1. Samarskii A. A., Galaktionov V.A., Kurdyumov S.P., Mikhailov A. P. Blow-up in Quasilinear Parabolic Equations. Berlin: Walter de Gruyter, 1995.
2. Wu Z.Q., Zhao J.N., Yin J.X., Li H.L. Nonlinear Diffusion Equations. River Edge NJ: World Scientific Publishing Co Inc., 2001.
3. Kalashnikov A. S. Some problems of the qualitative theory of nonlinear degenerate parabolic equations of second order, Uspekhi Mat Nauk, 1987. vol.42, pp. 135-176.
4. Deng K., Levine H. A. The role of critical exponents in blow-up theorems., The sequel J Math Anal Appl, 2000. vol. 243, pp. 85-126.
5. Galaktionov V. A., Levine H.A.On critical Fujita exponents for heat equations with nonlinear flux conditions on the boundary, Israel J Math, 1996. vol.94, pp. 125-146.
6. Rakhmonov Z. R., Tillaev A. I. On the behavior of the solution of a nonlinear polytropic filtration problem with a source and multiple nonlinearities., Nanosystems: physics chemistry mathematics, 2018. vol. 9, no. 3, pp. 1-7.
m
7. Rakhmonov Z. R., Khaydarov A. T., Urunbaev J. E. Global Existence and Nonexistence of Solutions to a Cross Diffusion System with Nonlocal Boundary Conditions., Mathematics and Statistics, 2020. vol.8, no. 4, pp. 404 - 409.
8. 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. vol. 9, pp. 236-245.
9. Zhaoyin X., Chunlai M., Yulan W. Critical curve of the non-Newtonian polytropic filtration equations coupled via nonlinear boundary flux, Rocky mountain journal of mathematics, 2009. vol. 2, no. 39.
10. Levine H.A.The role of critical exponents in blow up theorems, SIAM Rev., 1990. vol.32, pp. 262288.
11. Li Z.P., Mu C.L., Cui Z.J. Critical curves for a fast diffusive polytropic filtration system coupled via nonlinear boundary flux, Z. Angew Math. Phys. ,2009. vol.60, pp. 284-296.
12. Chen B., Mi. Y., Mu Ch. Global existence and nonexistence for a doubly degenerate parabolic system coupled via nonlinear boundary flux, Acta Mathematica Scientia, 2011. vol.31B(2), pp. 681-693.
13. Yongsheng M., Chunlai M., Botao Ch. Critical exponents for a doubly degenerate parabolic system coupled via nonlinear boundary flux, J. Korean Math, 2011. vol. 48, no. 3, pp. 513-527.
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. vol. 8(3), pp. 317-322.
15. Rakhmonov Z., Urunbaev J., Alimov A. Properties of solutions of a system of nonlinear parabolic equations with nonlinear boundary conditions, AIP Conference Proceedings, 2022. vol. 2637, no. 040008.
16. Rakhmonov Z., Parovik R., Alimov A. Global existence and nonexistence for a multidimensional system of parabolic equations with nonlinear boundary conditions., AIP Conference Proceedings, 2021. vol.2365, no. 060022.
17. Aripov M., Matyakubov A., Bobokandov M. Cauchy problem for the heat dissipation equation in non-homogeneous medium, AIP Conference Proceedings, 2023. vol.2781.
Information about authors
W
Alimov Akram Abdurashidovich A - Senior Lecturer of the Dep. of Information Systems and Mathematical Disciplines, Tashkent branch of the G.V. Plekhanov Russian University of Economics, ORCID 0009-0009-8518-4366.
Rakhmonov Zafar Ravshanovich A - D. Sci. (Phys. & Math.), Dean of the Faculty of Applied Mathematics and Intellectual Technologies of the National University of Uzbekistan, Tashkent, Uzbekistan., ORCID 0000-0002-4190-7069.
Информация об авторах
Алимов Акрам Абдурашидович А - старший преподаватель кафедры Информационные системы и математические дисциплины, Ташкентский филиал Института им. Плеханова Российский экономический университет, ORCID 0009-0009-8518-4366.
Рахмонов Зафар Равшанович математических наук, декан математики, и интеллектуальных университета Узбекистана, г.
ORCID 0000-0002-4190-7069.
А - доктор физико-факультета прикладной технологий Национального Ташкент, Узбекистан,