Научная статья на тему 'MODELING CROSS-DIFFUSION PROCESSES IN MULTIDIMENSIONAL AREAS'

MODELING CROSS-DIFFUSION PROCESSES IN MULTIDIMENSIONAL AREAS Текст научной статьи по специальности «Математика»

46
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
cross-diffusion systems of equation / asymptotic solutions / new effects / self-similar and approximately self-similar solution / kross-diffuziya tenglamalar sistemasi / asimptotik yechimlar / yangi effektlar / avtomodel va taqribiy avtomodel yechimlar

Аннотация научной статьи по математике, автор научной работы — Abrorjon Mamatov

In this paper, we study modeling cross-diffusion processes in multidimensional areas. It is proved that for a given equation, the parameter values exist and have a numerical solution. The system of equations considered in this paper is based on most physical processes, for example, the cross-diffusion process in this system, thermal conductivity, polytrophic filtration of gas and liquid in nonlinear media is described by a source. There are a lot of partial solutions to this equation. One of the main methods of studying the problem under consideration is the construction of an integral self-similar solution. To do this, initially, when constructing a system of self-similar equations, a nonlinear subtraction method was used. The following results were obtained from this work: the front for the equation of nonlinear heat generation with doubled energy was estimated, the localization process was observed, new effects were observed, an algorithm was constructed in accordance with the obtained self-similar solution, a program code was created in the programming language, and the process modeling was visualized

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

KO‘P O‘LCHOVLI SOHALARDA KROSS-DIFFUZIYA JARAYONLARINI MODELLASHTIRISH

Bu ishda ko'p o'lchovli sohalarda kross-diffuziya jarayonlarini modellashtirish tadqiq qilingan. Ko’rsatilgan tenglama uchun parametr qiymatlarining mavjudligi va sonli yechimga ega ekanligi isbotlangan. Ushbu ishda ko'rib chiqilgan tenglamalar sistemasi ko'pgina fizik jarayonlarga asoslangan bo’lib, masalan, bu sistemada krossdiffuziya jarayoni, issiqlik o'tkazuvchanlik, nochiziqli muhitda gaz va suyuqlikning politropik filtrlanishi manba bilan tasvirlanadi. Bu tenglamaning juda ko’plab xususiy yechimlari mavjud. Ko'rib chiqilayotgan muammoni o'rganishning asosiy usullaridan biri taqribiy avtomodel yechimini qurishdir. Buning uchun dastlab avtomodel tenglamalar sistemasini tuzishda chiziqsiz ajratish usulidan foydalanildi. Ushbu ishdan quyidagi natijalar olindi: ikki karra chiziqsiz issiqlik tarqalish tenglamasi uchun front baholandi, lokalizatsiya jarayoni kuzatildi, yangi effektlar kuzatildi, olingan avtomodel yechimga mos ravishda algoritm qurildi, dasturlash tilida dastur kodi yaratildi va jarayon vizual modellashtirildi

Текст научной работы на тему «MODELING CROSS-DIFFUSION PROCESSES IN MULTIDIMENSIONAL AREAS»

MODELING CROSS-DIFFUSION PROCESSES IN MULTIDIMENSIONAL

AREAS

Abrorjon Mamatov

National University of Uzbekistan, 100174, Tashkent mmtovabrorion1995@gmail.com

ABSTRACT

In this paper, we study modeling cross-diffusion processes in multidimensional areas. It is proved that for a given equation, the parameter values exist and have a numerical solution. The system of equations considered in this paper is based on most physical processes, for example, the cross-diffusion process in this system, thermal conductivity, polytrophic filtration of gas and liquid in nonlinear media is described by a source. There are a lot of partial solutions to this equation. One of the main methods of studying the problem under consideration is the construction of an integral self-similar solution. To do this, initially, when constructing a system of self-similar equations, a nonlinear subtraction method was used. The following results were obtained from this work: the front for the equation of nonlinear heat generation with doubled energy was estimated, the localization process was observed, new effects were observed, an algorithm was constructed in accordance with the obtained self-similar solution, a program code was created in the programming language, and the process modeling was visualized.

Keywords: cross-diffusion systems of equation, asymptotic solutions, new effects, self-similar and approximately self-similar solution.

KO'P O'LCHOVLI SOHALARDA KROSS-DIFFUZIYA JARAYONLARINI

MODELLASHTIRISH

ANNOTATSIYA

Bu ishda ko'p o'lchovli sohalarda kross-diffuziya jarayonlarini modellashtirish tadqiq qilingan. Ko'rsatilgan tenglama uchun parametr qiymatlarining mavjudligi va sonli yechimga ega ekanligi isbotlangan. Ushbu ishda ko'rib chiqilgan tenglamalar sistemasi ko'pgina fizik jarayonlarga asoslangan bo'lib, masalan, bu sistemada kross-diffUziya jarayoni, issiqlik o'tkazuvchanlik, nochiziqli muhitda gaz va suyuqlikning politropik filtrlanishi manba bilan tasvirlanadi. Bu tenglamaning juda ko'plab xususiy yechimlari mavjud. Ko'rib chiqilayotgan muammoni o'rganishning asosiy usullaridan biri taqribiy avtomodel yechimini qurishdir. Buning uchun dastlab avtomodel tenglamalar sistemasini tuzishda

March, 2023

503

ISSN: 2181-1385 ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

chiziqsiz ajratish usulidan foydalanildi. Ushbu ishdan quyidagi natijalar olindi: ikki karra chiziqsiz issiqlik tarqalish tenglamasi uchun front baholandi, lokalizatsiya jarayoni kuzatildi, yangi effektlar kuzatildi, olingan avtomodel yechimga mos ravishda algoritm qurildi, dasturlash tilida dastur kodi yaratildi va jarayon vizual modellashtirildi.

Kalit so'zlar: kross-diffuziya tenglamalar sistemasi, asimptotik yechimlar, yangi effektlar, avtomodel va taqribiy avtomodel yechimlar.

INTRODUCTION

Consider in the domain q = {(t, x): t e r+ , x e r} parabolic system of two cross-diffusion equations:

where > 1,m2 > 1, p > 2, qx, q2, r2 > 1,k, i e r is the parameters, u0(x),v0(x) is the initial conditions, |x|", |x|k is the density of the medium, o <r(t)e c(o,<») is the specified function.

Before constructing a self-similar solution of the system of equations (1), let us consider some cases of diffusion, for example: m, + p —3 > o, i = 1,2 -the state of slow diffusion, mt + p — 3 = 0, i = 1,2 -the critical state(the asymptotic is summed up depending on its solutions), m, + p—3 < 0, i = 1,2 is called the state of fast diffusion. An asymptotic solution is usually understood as a solution of a system of nonlinear equations that can satisfy certain conditions. Let us consider the reaction-diffusion system (1) for the case with a double nonlinear variable density and study its numerical solutions (by calculation methods). The equation (1) represents a number of physical processes [1]: the reaction diffusion process in a nonlinear environment, the heat dissipation process in a nonlinear environment, the filtration of liquid and gas in a nonlinear environment, they represent the existence of the law of polypore and other nonlinear displacements [2-3].

The Cauchy problem and boundary value problems for the equation were observed by many authors in one-dimensional and multi-dimensional cases [4-5]. The equation (1) in the processes represented by the phenomenon of finite distribution of temperature occurs [6]. In the presence of an absorption coefficient, the phenomenon of the "rear" front can occur, that is, the left front can stop after a certain time and move along the medium [7].

(1)

with initial (Cauchy) condition

w(0, x) = u0 (x)> 0, v(0, x) = v0 (x)> 0,

(2)

March, 2023

MAIN RESULTS AND SOLUTION METHODS

We can translate the system of equations (1) into a system of radial-symmetric equations so that we can find a solution to a self-similar or an approximately self-similar. To do this, we first introduce the notation as r = |x|, so that we can translate the system of equations (1) into a radial-symmetric system:

k CU rk — = div dt

k Cv rk — = div dt

rn+N ~lum-1

rn+ N-lvm2-1

du

p-2 du_' ^ dr

+ r(t )uqi vr rk

p-2 v ^

dr

(3)

+ r(t )uq vri rk

After performing the substitution (3), to find a self-similar solution of the system of equations (1) and the solution of the approximately self-similar, we use the following method:

fu(t, r ) = u (t )-®(r(t \p(r )) |v(t, r ) = v (t )• z(r(t ),y(r ))

Now we calculate the initial part of the system of equations (1), as required, as follows:

(4)

du dt dv dt

^ = r(t )• u* vri

— = r(t )• Uq2 vr2

u (t ) = A1 t To +\r(l)dri _ 0 _

v (t )= A2 t T0 + jr(rldr _ 0 _

1 - r2 + rj

1 (qi -j)(r2- j)- riq2 x _ 1 - qj + q2

2 (qi- j)(r2-j)- riq2

where it is equal to a1 =

1-r,

a a

(qj -l)(r2 -l)-riq2

a2 =

a1 qi a1 rj .After performing the

calculations, the system of equations (3) takes the following form:

da i - s d -= pl s —

dr dp

f

s-1 m, -1

p a 1

da

dp

dz i-s d — = ( s —

dr dp

f

ps-1zm1-1

dz

dp

P-2 * l^ da

dp

P-2 . l"\

dp

Uq1 (m1+p+l-3)vr1 (aq1 zr1 -a) K7q2l7 r2 (m2 + p + l-3)L q27 r2 _ 7)

(5)

+ r(t)uq2vr2(m2+p+l-3)(aq2z'2 -z

From the system of equations (4)-(5), (1)-(2) there are important considerations on the question: if mt + p+1 * 4; i = 1,2

t t r(t)=\[u{i1)}mi+p+l-4di = \[v{n)]m2+p+l 4dv

0 0

or if mt + p+1 = 4; i * 1,2 is equal to r(t)=t+1. From these considerations, it follows that the values of arbitrary and non-zero variables:

March, 2023

a

a

\

v

y

505

ISSN: 2181-1385

ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

(p{r ) =

p+k-n . P

p + k - n ln (r ) ; n = p + k

; n ^ p + k

s = p

N + k

f®(r,()=f (4)

[z(),() = g(4)

4 = (

) p ^ \

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

- 1/

4 =t 7 p

p + k - n ( = 4T 1 p

4,=- ! 4 p T

(6)

After substituting (6), the system of equations (3) takes the following form

P -2 jW^

__1 4 = ^i-s 1 _d_ p ) d4 t d4

1 dg = gi-s d

p d4 d4

^s-1 j' m1 -1

df

d4

df_

d4

f 4s-1 gm2 -1 V dg p 2 dgL 1

d4 d4

+ y(t)uqi-(m+p+l-3)vri (fqgri - f)

,"792,7 r2-(m2 + p + 1-3)(f q^r, r2 _ „)

(7)

+ r(f ))Uq2 -(m2 + p + l-3)(fq2 g2 - g

(7) by reducing (reducing) the system of equations, we form a system of equations of a new form:

1-s

4 d4

(

^s-1y mi-1

df

1-s

4 d4

(

4s-1 gm2 -1

d4 dg

p - 2 Í A p

d4

d4

p-v,

d4

+ 4r(t Tu9 -(mi+p+l-3 )vr1 (fq1 gr - f) = 0 p

+ 4r(t)uq2 V'2 -(m2 + p+l-3)(fq2 gr2 - g )= 0 p

(8)

To find a solution to this system of equations (8), we introduce another repeating self-similar pattern:

\f (#)=f [g(#) = B(a2)4

where a, > 0; i = 1,2, yt > 0; i = 1,4 is equal.

To facilitate the calculations, we calculate the first equation of the system of equations (8):

^s-1 J1 m1 -1

df

d4

p 2 f = 4s-1 ^m,-1 « _4r. f2-1)

Ar2r4r1 -1(«1 -4

(«1 -4n -1

p-2

{- A'^r^r -1 («1 -4» /2 (l-1)) =

= -Am1 + p +l-3(r)l + p-1(rJ + p-1 '/4s-1 + (n-1)(p-2)+(n-1)(a (m1-1)+ (r2-1)(p-2)+r2l-1

!1+p+i-3fv V-^ V-^s-1+(n-1)(p-1)(

= - Am1+p+l-3 (r2 )p-1(r1)p-114s-1+(r1-1)( p-1)(fl1 -4r1 r

(l+m1 + p-3)+1-p

1-s d

s-1 r m1 -1

d4

= -4

m1 + p+l-3

4s-1 f

df

d4

p - 2

f.

d4

+ 4-J + d1(f91 gr1 - f )= p d4

(9)

Am1 + p+l-3(r2)p-1(r1 )p-1i(s - 1(r1 - 1)(p - 1))4s-2+(n-1)(p-1)(a1 -4r1 r2(l+m1 +p-3)-p ■ Am1+p+l-3 (r1 )p (r 2 )p-1l (r2 (l + m + p - 3)+1 - p)4s-1+(r 1-1) p («1 -4 r1 )2 (l+m1+p-3)- p" i A r1 r24r1 («1 - 4r1 r2 + d f Aq Bq1 («1 - 4r1 r2 91 («2 - 4 r3 )l - A« -4r1)21 = 0

March, 2023

I

y

V

V

V

y

/

506

Now we will find the unknown parameters from equality (9):

(ri " 1)p = 71 = p_r p -1

72 (l + mi +1 - 3)p = 72 - 1 ^72 =

p-1

l + m1 + p - 4

+p+l3 (71 )p (72 )p-1l (1 -p) = 1.47172 (l + m1 + p - 4) p

(10)

A =

(7 )1-p (7 )2- p (l + m1 + p - 4) (71) (72) pl(1 -p) _

1

m1 + p+l-4

As can be seen from the calculation equations, the following parameters are found from the second equation of the system of equations (9):

73 =

p

74 =

p -1

4

p -1 1+m+p - 4

; B =

(7 )1-p (7 )2-P (l + m2 + p - 4) 7) (74) pl(1 - p) ,

1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

m+p+l-4

Combining all the calculated equalities, we get the integral self-similar solution we are looking for:

1-r,

uA (t, x) =

Va (t, x)=[a1-]

(?1 -1Xr2-1)- rq 2

To +\7(il)dv

0 +\7(li)dl 0

B(a2 -#73 )-

a(«1 -#71 )2

(11)

Since all the parameters in equality (11) for which an approximately self-similar of the solution is found, we now see an asymptotic process for some special cases. We have described the asymptotic in detail above. If, in the above equation (8), the following change is made to the process of finding a calculated self-similar solution:

7(t)zuq1 (m1+p+l 3)vr1 ^ const., t —y ^ 7(t )zUq2 Vr2-(m2 + p+l-3)-

(12)

^ const., t

If equality (12) holds, then we can imagine how the solutions of the self-similar

(11) we have found will change. The condition (12) introduced by us is now considered the asymptotic state of equality in case r(t )= const. First, we check the system of equations

(12) through a new limit condition:

f (0)= ci > 0, f (d) = 0, a3.

g(0) = C2 > 0, g(d) = 0. ( )

where is 0 < d <+» . (12)-(13) For the problem in r(t)=0, n = 0,1 = 0, p = 2 cases, that

the solutions have a trivial self-similar solution and existence properties [4], [9], [11] it is

quoted in the works.

To find out in which cases fast and slow diffusion occurs in the

system of equations (12), we perform a substitution in the form:

March, 2023

a

a

0

507

^SSN:2181-138^^^^SI:0,96Tj^ite-Factor:0,89j SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

ff (4) = f {4)yM

(

g(4) = g (4K (v)

^ V = - ln

a-4 p -1

(14)

were

f (4)=a

N 7 2

ai -4

P-1

= Ae - 7 2V

g (4)=b

/

N74

(15)

a2 -4

P-1

= Be' 74

or it can also be replaced with another form:

if (#)= cj (#) [g (#)= c2 g {4)

Combining both methods, we form the following system of equations:

(16)

s-1 ¿-m, -1

4s-1 f

4 s-1 gm2 -1

df p- -2 -; df

d4 d4

dg p- 2 dgl

d4 d4

f + 4s (72 )P-1f e C(0, »)

(17)

or it can also be written in another view

d f 4s- V 1 y m1 -1 df p-2 - 7 P df

d4 d4 d4

d / 4s-1 gm-1 dg P-2 j-i dg

~d4 d4 d4

= -(72 )P-1

sf+4 f

^ )P-1| sf + 4 %

(18)

We will carry out this system of equations (18) separately, putting two equalities. After calculating the first equality, the following result is obtained:

^s-1y mx-1

P-2 l

f = 7 p-14f (y1 )

d4

h (y1 ) = ym1 1

dy1

dv

2 y1

P-2 l

dv

2 y1

From the second equality, the following result is obtained:

s-1 ^ m1 -1

4s-1g

dg

d4

L2 (y2 ) = ym2-1

p-2 , i

^ = 7 P-14sgh2 (y2 ) d4

dy2 dv

" 74 y 2

P-2f J l \

dy2 l

-f2 - 74 y1 dv

Combining the obtained results, we obtain a new system of equations:

March, 2023

P

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

P

4

1—s

1—s

\

/

508

ISSN: 2181-1385

ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

1_s d

(

rx fn

df

1_s d

(

£5-1 gm2 _1

dg

-2 - i

p_2 i— i dg

= WP_7

y

s _T2Ti

= (r)^

y

s _7471

fl _£71 £71

hy1 + 71

fl _£71 dq

a

L2 y2 + 71

a _£71 dq

L1y1

L2 y 2

(19)

If we also consider the newly formed system of equations (19), divided into two parts, and then the following result will be obtained in the first part:

d(Vi)+f— Vy + — r-P^i(lidyL)-7-rv r2 1-r-p^i(1)yi +

di Ui ) p I. di J (qi-1)(r2-1)-ri?2

,(_72 +174 +fl7-1)q

+ 7-

y? y2 = 0

Then the following result is formed in the second part:

d (l2 y 2 )+í-^ 9 (q)_74 )l2 y 2 +—71 >1 (qf-^ _ 74 y2 dq í 71 ) p I. dq

?2 _ +1

7-TV?-TV-71 >1 (q)y2 +

(?1 _ 1)(r2 _ 1)_ r1?2

,(_/4 + r274 + fl272 _1)q

+ 7_

,_q

y?2 y22 = 0

where <p1h)

91

(q)=-

a _ e

this will be equal to. Instead of concluding, we can say that the function <P1h) (12), (13) serves as a property of the asymptotic solution of the problem. We present two theorems that combine all the calculated results. Theorem 1. k1 > 0, k2 > 0 let be given. Next

[f (£) = h°f (#)(1 + o(1))

L + o

[g (£) = h0 g (^)(1 + o(1)) the system of equations (20) will have an asymptotic at i^+x

1_ a p

v

(20) point.

j

Where it will be equal to 0 < h < +x if one of the following conditions is met:

1) If condition st > —m is satisfied, then will h0, ) be the solution of the following

p -1

nonlinear system of algebraic equations with roots (h1, h2)

h0 Y11 (hf)- 2=C1

(h? Y2 1 (h0 )p- 2 = C2

(21)

were c, =

p(zk, )p 1

(21), we obtain the following system of equations:

= 1,2 . By reducing the given values to the system of equations

March, 2023

d

d

V

n

a _ e

p

n

e

n

1

509

ISSN: 2181-1385

ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

h = 1

( 1 ^

xh f xpk

p - 2

xp - 2 pkp-1 1_

p -2

k1(xk1p)

m1 - 2 p

h2 =

k1(Xk1p)

p - 2

2) Given that 2 = of algebraic equations

k -1

i = 1,2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2, (h0, h0) are the roots of the following nonlinear system

kp-1(h10 ^+p-3 +

1(h10 f -1(h0 f-1

«Xpk1

px

p-1

kp-1(h2

(Äo )»2+p-3 + «2 (h10 f -1 (h0 f -1

(22)

px

p-1

Proof. To prove Theorem 1, we use the substitution (14). When using the self-similar finding of the solution (6), the following results are obtained similarly to the substitution (14):

— L(h1, h2 )+|- ((V)-k1 lL1(h1, h2 ) + dv VX )

t

(,)+«_+(1^ )= 0,

where p(v)=-

(23)

(h1 h2h|- v(v)-k2 1l2 (h1 h2) + dv Vz )

t

+p(v)+^ pvifh* +^2 (vh2 )= 0.

zp zp

p (v)=2-1), i = 1,2. after performing the calculations, the system

of equations (23) is compressed as follows:

Afo, h2 )=v1-1

dh1 - k1h1 dv

dT - k1h1

dv

L2 (h1, h2 ) = h2m2 1

dh2 - k h

j k2h2 dv

\p-2

h 2 k2h2 dv

As can be seen from the substitution (14) that we calculated, the properties of the asymptotic solution in (12) v will be appropriate when there is an arbitrary in the system of equations (23), since the following conditions are appropriate

- kfr * 0, h (v)> 0, i = 1,2.

dv

First, we show that h1(v), h2 (v) are the finite roots (solutions) of the system of equations (23). We enter the markup as follows. At

v

March, 2023

m-. - 2

p

k

1

1

n

e

n

a - e

2

510

ISSN: 2181-1385

ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

b, (i) = L, (h1, h2), i = 1,2.

Now we can write the system of equations (23) in its new form

¿1 = fS l k1 )b + ^^ 1 + ¿1 (l)hf1 )= 0

X'

b/ =í— 9q)_ ^2 V + 9(q) + fl^ +92 (qh2 )= 0

Vx J xp 1 Xp

X

To make the system of equations more compact, we will make changes with new functions

V^,?)=-[-1k1 y -^^¿1)--^11 +^1(l)hf1 )= 0

X

X

wi("iihA-¿(i)-k21^1'-h2 ^¿(i)-^^¿(i)(h2 + ¿2(ii)h<22)=0

VX ) X X

where y,, = 1,2 is a real number. The functions ui given in the above equation will be appropriate in any of the following conditions given for each of the

ly ,+x) c 1i0 ,+x)(0 < 10 < ly ) intervals 1 e ly ,+x)

W, ,1)> 0, W2(y2,1)<

In the system of equations (23), the boundary of the function b1 (1) is located on the interval 1 e 1 ,+x) (taking into account the theorem), and the following data are relevant:

lim b, (1) < +x, lim b, (1) = 0 .

1—1—

From the above conclusions, we get the following result:

lim hj (1)= h0 < +x, lim h, (1)= 0 .

1—+x 1—+x

Now we will put the obtained results in a system of substituted equations, from where we will begin the proof of the theorem:

lim b1 (q)= lim

q^+x q^+x

lim b2 (1)= lim

9(q)_ k)n _ p(q) 9(q)(h1 + 91 (q]^ )

J xp xp

= 0

q^+w

q^+w

_| — 9(q)_ k2 n

' h-, _ knh n _

pí2^ 9(q)_ —p 9(q)(h2 +92 (q>22)

X X

(24)

=0

Taking into account the obtained limit intervals, we pass from the system of equations (24) to the system of algebraic equations (21):

f(h2° )m1 -1 (h10 )p-2 = C1

, 0W -

(h10 )m2 _1 (h20 T

kp_1 (h2 )T _1 (< )p_2 + fl1(h10 ^ 1

flxpk1 px

kp_1 (h

1 (h0 )">_1 (h20 )p_2 + fl2 (h20 l22_1

(25)

flXPk2

px

p_1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

March, 2023

2

1

511

where is s =

k -1

i = 1,2 . It can be seen from the results obtained that we have

obtained the asymptotic form (20), that is, the theorem is proved. From Theorem 1, we can conclude that in case k, < o, i=1,2, rapid diffusion occurs. Therefore, we investigate the regular asymptotic of the fast diffusion state in the self-similar solution (20) and

determine the following limiting conditions

j/(o)=ki > 0, /(») = 0 (26)

[g(0) = k2 > 0, g(«) = 0 V 7

In (12), we perform the following replacement

j/(c)=/ (fab) (27)

U (C) = g (fa2 b) ( )

all the coefficients listed here have been determined above. Theorem 2. Solutions of the system of equations (20) under the given conditions k1 <0 ,k2 <0 ,St > 1, i = 1,2 , b ^ +œ (C will have the following asymptotic form

(/(c)=h0/(C)(i+0(1)) U (c)=h0 g (C)(1+0(1))

(28)

here it will be equal to 0 < h0 < +<» if one of the following conditions is met: 1) If condition (N -1)[(m, - 1)(m3_i -1)-(p - 2)2 J-(p - n -1)(p - mt -1)> 0

and

a, >

P -(m, +1)

,1 = J

1,2 are satisfied, then (h0, h0) will be the solution of the

' (mi - 1)(m3-i - 1)-(P - 2)2

following system of nonlinear algebraic equations with roots (h1, h2)

a, <

2) If

p -(m, +1)

ffc )m1 -1 (h? )P-2 = C1 (h? )m2-1 )P-2 = C2

Condition (N -1)[(m, - 1)(m3-i -1) - (p - 2)2 J- (p - n -1)(p - mt -1) < 0

(29)

and

,1 = 1

1,2 are satisfied, then (h0, h0) will be the solution of the

' (mt - 1)(m3-; - 1)-(p - 2)2

following system of nonlinear algebraic equations with roots (h1, h2)

(s+#1 )(kA°| )p-2 (h20 J"1 -1 +

(s + #k2 ^ )p-2 (h? )* -1

p-2

p-1

p-2

k2#

p-1

= 0

= 0

(30)

Proof. We substitute (12) into the system of equations (27) initially, and we get the following form:

py

d L (h1 > h2 ) ++ | - (v) - k (hu h2 ) + dv vy )

1pT^UK -k1h1 l+O-<P1 (v)(h1 +v2lVh)=0,

v ) y-^

(31)

March, 2023

k

1

a

1

a

+

+

512

where is 1v(i)=-

a - e

V(l) = e-ik(S-1W(i), i = 1,2, and

L, (h1, h2 )= h3-i-1

dht di

\p-

+ kh

dh

di

L + kh

As a consequence, the study of the properties of the solution of the asymptotic system of equations (14) in (31) is the same as the study of the system of equations (2) around, since the following condition is always appropriate

dh, di

+kfr * o, h (i)> o, i=1,2.

Thus, we pass from the system of equations (31) in 1—+x to the system of algebraic equations (30) with the necessary conditions. The proof of Theorem +x occurs in this way.

CONCLUSION

The results of computational experiments show that the iterative methods listed above would be effective in solving nonlinear problems. Nonlinear effects result if the nonlinear division method and linear self-similar solutions are used as the initial approximation of the solution, in which the standard equation is constructed in a functional way. As expected, to achieve the same accuracy, the Newton method requires less iteration compared to the Picard methods and the special method due to the successful choice of the initial approach. In each of the cases considered, the Newton method has the best approximation by choosing a good initial approximation. In some cases, the total number of iterations is almost two times less, and the maximum number of iterations is almost 4 times less than in other methods. The results of numerical calculations show the influence of the numerical speed of distortion propagation, and the localization of the resolution depends on the values of the numerical parameters. All the results of digital experiments are visualized.

ACKNOWLEDGMENTS

We are very grateful to experts for their appropriate and constructive suggestions to improve this template.

REFERENCES

1. Angar, J. (2017) Cross-Diffusion systems with entropy structure. arXiv: 1710.01623v1 [math.AP] 4 Oct 2017. Proceedings of equadiff.

2. Aripov, M., Sadullaeva, Sh. (2013) To properties of solutions to reaction diffusion equation with double nonlinearity with distributed parameters. Jour. of Siberian Fed. Univer. Math. & Phys. 6, 157-167.

March, 2023

n

513

ISSN: 2181-1385 ISI: 0,967 | Cite-Factor: 0,89 | SIS: 1,9 | ASI: 1,3 | SJIF: 5,771 | UIF: 6,1

3. Galaktionov, V., Kurdyumov, S. (1979) On comparison of solutions of parabolic equations. DAN an SSSR. 248(3), 586-589.

4. Kurdyumov, S., Zmitrenko, N. (1977) N-and S-modes of compression of the final plasma mass and features of modes with sharpening, PMTF, 5(1), 3-23.

5. Mamatov, A., Nurumova, A. (2022) Modeling of nonequilibrium cross-diffusion systems, J. Phys.: Conf. Ser. (16.12.2022 2388, 012027).

6. Mamatov, A., Nurumova A. (2022) Investigate the Interaction of Nonlinear Thermal Conductivity Process with Density and Source in a Nonhomogeneous Environment. ICPPMS-2021, International conference, Scopus & Web of Science indexed (16.06.2022 2432, 050031)

7. Mamatov, A. (2020) Numerical simulation of the heat conduction process with double nonlinearity in the two-dimensional case in the presence of transfer and source (absorption). Scientific Journal of Samarkand University, ISSN: 2091-544, No. 1, 41-44.

8. Martinson, L. (1983) Evolution of a heat pulse in a nonlinear medium with volumetric heat absorption. TVT, 21(4), 801-803.

9. Martinson, L., Pavlov, K. (1972) On the spatial localization of thermal disturbances in the theory of nonlinear thermal conductivity. ZhVM and MF. 12(4), 1048-1053.

10. Mamatov, A. (2020) Modeling the effect of the two-fold non-linear heat dissipation equation on biological population with ambient density, Scientific Journal of Samarkand University. ISSN: 2091-5446, 5(123). 37-40.

11. Rakhmonov, Z., Urunbayev, J. (2019) On a problem of cross-diffusion with nonlocal boundary conditions. J. Sib.Fed. Univ.- Math. Phys.5(7), 123-128.

12. Samarsky, A., Zmitrenko, I. (1975) Effect of metastable heat localization in a medium with nonlinear thermal conductivity. DAN USSR, 223(6), 134.

13. Samarsky, A., Elenin, G. (1977) Gorenje nonlinear environment in the form of complex structures. DAN USSR, 237(6), 1330-1333.

14. Tedeev, A. (2004) Conditions for the time global existence and nonexistence of a compact support of solutions to the Cauchy problem for quasilinear degenerate parabolic equations. Siberian Math. J. 45(2), 155-164.

15. Wang, M., Wei, Y. (2008) Blow-up properties for a degenerate parabolic system with nonlinear localized sources. J. Math. Anal. Appl. 343, 621-635.

March, 2023

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.