Научная статья на тему 'Numerical investigation of solutions to a reaction-diffusion system with variable density'

Numerical investigation of solutions to a reaction-diffusion system with variable density Текст научной статьи по специальности «Математика»

CC BY
73
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДВОЙНАЯ НЕЛИНЕЙНОСТЬ / DOUBLE NONLINEARREACTION-DIffUSION SYSTEM / АВТОМОДЕЛЬНОЕ УРАВНЕНИЕ / SELF-SIMILAR SOLUTIONS / АСИМПТОТИКА РЕШЕНИЙ / ASYMPTOTICS / NUMERICALCALCULATIONS / СИСТЕМА РЕАКЦИИ-ДИФФУЗИИ / ЧИСЛЕННОЕ РЕШЕНИЕ

Аннотация научной статьи по математике, автор научной работы — Sadullaeva Shahlo A.

In thispaper we demonstrate thepossibilities of the self-similar and approximately self-similar approaches for studying solutionsofa nonlinear mutualreaction-diffusion system. The asymptoticbehaviorofcompactly supported solutions and free boundary is studied. Based on established qualitative properties of solutions numerical computation is carried out. The solutions are presented in visualization form, which allows observing evolution of the studied process in time.

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

Текст научной работы на тему «Numerical investigation of solutions to a reaction-diffusion system with variable density»

УДК 519.21

Numerical Investigation of Solutions to a Reaction-diffusion System with Variable Density

Shahlo A. Sadullaeva*

Tashkent University of Information Technology Amir Temur, 108, Tashkent, 700084, Uzbekistan

Received 15.10.2015, received in revised form 06.11.2014, accepted 30.12.2014 In this paper we demonstrate the possibilities of the self-similar and approximately self-similar approaches for studying solutions of a nonlinear mutual reaction-diffusion system. The asymptotic behavior of compactly supported solutions and free boundary is studied. Based on established qualitative properties of solutions numerical computation is carried out. The solutions are presented in visualization form, which allows observing evolution of the studied process in time.

Keywords: double nonlinear reaction-diffusion system, self-similar solutions, asymptotics, numerical calculations.

DOI: 10.17516/1997-1397-2016-9-1-90-101.

Introduction

In the domain Q = {(t,x) : t > 0, x € Rw} properties are investigated of the process of a nonlinear diffusion-reaction with variable density described by the following system:

д (p(x)u)

dt

д (p(x)v) dt

= div (\x\n vm1-1 \Vu\p-2 Vu^j + p(x)1(t)uei, = div (\x\n um2-1 \Vv\p-2 Vv) + p(x)j(t)ve2,

u (0, x) = uo (x) ^ 0, v (0, x) = v0 (x) > 0, x e RN, where m1, m2, n e R, в1, ^ 1, P ^ 2 are given positive numbers.

V(-) - grad(-), uo (x) ,vo (x) > 0, x e RN,

(1)

(2)

X

p(x) = \x\-1, l> 0, 0 <7(t) e C(0, to).

System (1) describes different physical processes in two componential inhomogeneous nonlinear media, for example, the processes of mutual reaction-diffusions, heat conductivity, the theory of combustion, the theory of a polytrophic filtration of a liquid and a gas in the presence of a source whose power is equal to p(x)j(t)u^1, p(x)j(t)v^2. In works [1-6] were considered particular cases of the system (1), when 7(t) = 1, n = l = 0, p = 2 and we investigated arising different types of solutions depending on the parameters of system (1).

System (1) is degenerate in the domain, where u = v = 0 and may have no classical solutions. Therefore, the weak solutions of system (1) are studied, having physical sense. Namely functions 0 < u,v € C (Q) such that \x\n vm1-1 \Vu\p-2 Vu, \x\n um2-1 \Vv\p-2 Vv € C(Q) satisfying the

*[email protected] © Siberian Federal University. All rights reserved

integral identity in the sense of distributions [1]. For solutions of system (1) may take place the phenomena of finite velocity of propagation and space localization of disturbance, i.e. there may exist functions h(t), l2(t), such that u (t,x) = 0 and v (t,x) = 0 if \x\ > l1 (t) and \x\ > l2 (t). In the case li(t), l2(t) < for t > 0 a solution to problem (1)-(2) is called a space localization of disturbance. The surfaces \x\ = l1 (t) and \x\ = l2(t) are called a free boundary or a front.

In the present work we suggest a method of construction of self-similar equations system (1) based of splitting of system (1), and study asymptotics of compactly support solutions and a free boundary and asymptotics of self-similar solutions for the quick diffusion case. It is shown that the coefficient of the main member of the asymptotics of the solution satisfies a certain system of nonlinear algebraic system equation. Based on established qualitative properties of the solution, using approximately self-similar solutions, numerical experiments, visualization of processes described by reaction-diffusion system (1) with variable density were carried out.

1. Construction of a self-similar system of equations

Studying different properties of solutions to system (1) is a complicated problem, even for particular cases of system (1) [2,6-9]. In these works for particular cases of system (1)-(2) the effectiveness was shown of the self-similar approach for studying different properties of solutions to problem (1)-(2). Below a new method (nonlinear splitting) of construction of a self-similar and an approximately self-similar system is suggested. This method gives us a more simple way of investigation of qualitative properties of solutions to problem (1)-(2).

For construction of a self-similar and an approximately self-similar system for system (1) solutions u(t,x), v(t,x) to system of equations (1) are searched in the form

u(t, x) = u(t)w (t(t), \x\)), v(t,x) = v(t)z (t (t),f(\x\)),

(3)

where

u(t)

T + Y(t)dt J 0

i

' ßi-i

v(t)

T + i(t)dt

0

i

ß2-1

, T > 0,

and functions t, ^ will be choosen below.

Substituting (3) in system (1) reduces it to the following system of equations

dW = S-' d zm-

ot o<£ (

dw p- 2 dw

d<p ds

dz p 2 dz

d<p ds

(4)

where the functions t, ^ are chosen as

(t)= i vmi-1(n)up-2(n)dn = i um2-1(n)vp-2(n)dV

00

(5)

V(r) = — \r\pi , \r\ = A pi \

N 2 p - (n + l) N - 1

y x2, pi = -, s = p---—, l + n<p.

¿.—/ 1 p -

i=i

p - (n + l) '

It is easy to establish that system (4) has approximately self-similar solutions of the form

w(t,v) = f (î),

z(t,v) =

t

w + w

where £ = , and the functions f, ^ satisfy to the following approximately self-similar system

t 1/p

of equations

d-

(ts-iri— dfp-2 df)+-df+y № (t)u^-i)v-mi-i) (f+fßi )=o, y d- d- J p d-

i v * \ (7)

-s-1fm2-1 tP ' + + Y(t)T {t)ü-m*-1)vß-v-1) ' + ^ ) =0, y d— d— J pd—

-1-s dd—i-'-'r21

It is easy to provej that

Y(t)u^1-(v-1)v-(mi-1)r(t) ^ const, as t ^ <x>, Y(t)u-(mi-1)v^2-(v-1)T(t) ^ const, as t ^^

if 0 < y(t) G H, where H is Hardy's body [2]. In this case system (1) becomes self-similar. Therefore it is possible to say that (7) is an asymptotically self-similar system for system of equations (1).

Let y(t) = const. Then approximately self-similar system (7) has a self-similar form if

(P2 - 1)(p - (mi + 1)) = (Pi - 1)(p - (m2 + 1)) (9)

In this case for the functions f, ^ we have the following self-similar system

-1-s d- (-s-1fm2-1

df p-2 dA + if

d- d- )

d' p-2 d'\ + ^d' P d-

d- d-.

0,

(10)

where

ai = ($i - 1)$ - 1)/[($i - 1)$ - 1) - ((m1 - 1)($i - 1) + (p - 2)032 - 1))],

a2 = ($1 - 1)($2 - 1)/[($i - 1)($2 - 1) - ((m2 - 1)($2 - 1) + (p - 2)($i - 1))].

We notice that in the singular case $1 = (m1 + p - 2), $2 = (m2 + p - 2) positive solutions system (10) in the case of one equation were studied in [9].

2. Asymptotics of self-similar solutions

Now we will study asymptotics of weak compactly supported solutions (c.s.) to system (10) when j(t) = const.

Consider system of equations (10) with the following boundary condition

f (0) = ci > 0, f (b) = 0,

v>(0) = c2 > 0, ^(b) =0, ()

where 0 <b < .

The existence of a self-similar weak solution to problem (10)-(11) for one equation, in the case j(t) =0, n = l = 0, p = 2 was studied in [1] and conditions for the existence of c.s. solutions were obtained. Asymptotic behavior of a self similar c.s. solution ere established for the case

0

2

one equation in [10] and for other system, but with other nonlinear coefficients in [3-6,11] when p=2.

Now we will study the asymptotic of a c.s. solution system (10) in the case p + m1 — 3 > 0, p+m2 —3 > 0. For this goal we will transform system (10) to the convenient form for investigation by using the following transformations

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

f (o = mvi(ri),

= ■>v(î)У2(n),

n = - ln (a - , 0 <£<a1/x,

(12)

where

rn = (a - ey;, m = (a - er:, a =

p -1 *

(p — 1)(P — (mi + 1)) (p — 1)(p — (m2 + 1))

qi =-, 92 =-,

q q

p > mi + 1, p > m2 + 1, a > 0, q = (p — 2)2 — (mi — 1)(m2 — 1).

Case (p — 2)2 — (mi — 1)(m2 — 1) = 0 is a singular case. In this case the asymptotic of a solutions of the system (10) have another behavior. This case is required an additional investigation.

Theorem 1. Let qi > 0, q2 > 0. Then c.s. solutions to system (10) as n ^ ^ a1-p j

have the asymptotic

f (0 = yl№(1 + o(1)),

m = y0°$ (m+o(1)), )

where 0 < y2 < (i = 1, 2), if one of following conditions is satisfied:

1) fii > q-, i = 1, 2, then (y°, y02) are the real roots of the following system of nonlinear

qi

algebraic equations

(y0)mi-1(y0)p-2 = ci, (y0)m2-1(y0o)p-2 = C2,

1

p(Aqi)p-i

(14)

i = 1, 2,

y1 = ( Ap-2qP2-1p \

\qi(Xpqi) p-2 J

1

y0

I qi (Apqi) p-2 \ ^ Ap-2qp-1p )

Apqi

p-2

p-2 1 Aqi '

where qi, i = 1, 2 are the numbers defined above;

2) ¡3i = q-, i = 1, 2, then iy/2, y2^ are the real roots of the the following system of

qi

nonlinear algebraic equations

qp-1(y0)mi-1(y0)p-2 +

qp-1(y0)m2-1(y0)p-2 +

ai(y0 )ßi'

aApqi

pAp

i

a2

(y20)ß

0 ß2 1

aApq2

pAp

-l

c

2

q

1

1

1

Proof. In order to prove Theorem 1 we use transformation (12). Transformations (12) reduces self-similar system (10) to the following form

dnLi(yi,y2) + (- qi) Li(yi,y2) +y1 ^J?1 m + ap tn (yi + tiM1) = o,

-^ЬЛУ^Уз) + ( ^ф(П) - L1(y1,yi) + V2 -pq2iy2 ф(П) + jp Ф(л) (y2 + Ф2(п)У22) = 00,

(16)

where

ф(п) = , ФМ = e-™(e'-1), i =1, 2.

Li(yi,V2)= У?1-1(\У[ - 41Ух\)р-2(y1 - 41 Ух),

L2(yi,y2) = ут2-1(\у2 - 42у2\)р-2(у2 - 42y2),

and X, 4i, aiy i = 1,2 are the numbers defined above.

Such transformation (12) allows us to reduce studying of the asymptotics of solutions to system (10) as n ^ œ to studying those solutions of system (16), which in some neighborhood of +œ satisfy the inequalities

^ - 4iyi = 0, yi(n) > 0, i = 1, 2. an

First we show that solutions yx(n), y2(n) to system (16) have finite limits as n ^ œ. We introduce the notation

vi(n) = Li(yi,y2), i = 1, 2. Then system (16) may be rewritten in the form

vi = - (jt(v) - qi) vi - yiXpqiyi t(n) - jpt(n) (yi + ti(n)yf1), v2 = - (jt(v) - q^ v2 - y2Xpm1y'2 t(n) - jppt(n) (y2 + MvWl2).

To analyze the solutions to this system, we introduce the auxiliary functions

Si(pi,v) = - (jt(v) - qi) vi - yiXpqi^i t(n) - jppt(n) (yi + ti(v)yfl), 02(^2,n) = - (jt(n) - q^ P2 - V2Xpq2y2t(n) - jpt(n) (y2 + MvWl2),

where pi, i = 1,2, are real numbers. For each value of pi the function 0i(pi,n) preserves the sign on some interval [n^, C [n0, (0 < no < V^i) and for all n G [n^, takes place one of the inequalities

9i(pi,n) > 0, 9i(pi,n) < 0. Given the theorem of Bohl [13] for the functions vi(n) there are a limit in n G ,

lim vi(n) < lim vi (n) = 0.

Therefore

lim yi(n) = y0 < lim yi (n) = 0.

Then we have

lim vi = lim ^ - ( sA(n) - qO 1i - ^' - qiy A(n)--1 A(n) (y + a. (n)yßi

r^+o r^+o [ \A ) Ap 1

~A(n) - apA(n) [yi + A'(n)yi

v2 = ri_:moo j - [j^ri - q^j 12 - V2Xpq1y2 A(n) - ApA(n) (y2 +

0.

Note that as n ^ +œ

l im A(n) ^ 0, l im A(n)Ai(n) =

(17)

0, i f ßi > q-1, i = 1, 2,

q qi 1

............. I 1/a, if pi = q-1, i = 1, 2.

qi

With into account the latest limits we will obtain from (17) the following system of algebraic equations:

(y°2)mi-i(y°1)p-2 = ci, (y°1)m2-i(y°2)p-2 = C2,

(18)

if ßi > q—1 and

qi

qp-i(yï)mi-i(yi)p-2 + q2p-i(y0)m2-i(y0)p-2 +

ai(y0)ßi

aApqi

a2(yQ)ß2-aApq2

1

pjp-i ' 1

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

pjp-i,

(18')

if ßi

qi -1 qi '

Therefore we have the asymptotical representation (13). □

In the case p = 2 or m = 1 in (10), the properties of the different solutions as computing aspects of the system equation (10) were studied by many authors [12].

3. Fast diffusion case (qi < 0, i = 1, 2.)

Consider now the case of fast diffusion. In this case we will study asymptotics of the regular solutions of self-similar system (10) as Ç ^ +œ with the boundary condition

We will replace in (10)

f (0) = ci > 0, f (œ) = 0, ■ (0) = c2 > 0, ■(œ) = 0.

f (0 = f(Ç)yi(n), ■(0 = ■>V(0У2(n),

(19)

(20)

where n = ln (a + ), f = (a + )qi, tp = (a + )q2, a > 0, j, q\, q2 are the numbers defined above.

Theorem 2. Let q\ < 0, q2 < 0, Pi > 1, i = 1, 2. Then solutions of system (10) have the asymptotics

f (0 = yf (0(1 + o (1)),

■ (0 = y0$ (0(1 + ° (1))

0

i

i

as n ^ œ when one of the following conditions is satisfied:

I.

and

or

II.

and

(N - l) (mi - 1) (m-i - 1) - (p - 2)2 - (p - n - l)(p - mi - 1) > 0

ai >

p - (mi + 1)

(mi - 1) (m3-i - 1) - (p - 2)2

1, 2,

(N - l) (mi - 1) (m- - 1) - (p - 2) - (p - n - l)(p - mi - 1) < 0

ai <

p - (mi + 1)

(mi - 1) (m3-i - 1) - (p - 2)2

1, 2,

numbers y2, y2 are the roots of the following nonlinear algebraic equations

(s + Aqi) (\qiy0\)p 2 (y0 )mi-i + (s + Aq2) (\q2y°0\)p-2 (y0)m2-i +

1

ai

pjp-2 1

pjp

2

qiAp-i a2

q2Ap

i

0.

Proof. After replacement (20) system (10) take the following form

dnLi (yi,y2) + ((n) + qi) Li (yyi,y2)

+Ypti (n) (vi + t2 (n) y?') =0,

where

+Ai (n){ di+qiyi) +

ßi '

Li(yi,y2) = y?

er

i

dyi

~r + qiyi

dn

p-2

dyi

ym+qiyi)

Ai(n)

, ^2(n) = e

nqi(ßi-i)

Ai(n), i = 1, 2.

(22)

(23)

en - a

Such studying of solutions of a system (10) as n ^ +œ is reduced to studying of those solutions of system (10), each of which in some neighborhood of œ satisfies to inequalities

^ + qiyi = 0, yi(n) > 0, i = 1, 2. an

Passing in (23) to the limit at n ^ +œ we obtain necessary conditions and system of algebraic equations (22). Rest of a proving is similar to above mentioned manner. □

4. Case (1 = (2 = 1.

Now consider the Cauchy problem for system

d (p(x)u) dt

d (p(x)v) dt

= div (\x\n vmi-i \Vu\p-2 Vu^ + p(x)Yi(t)u, = div (\x\n um2-i \Vv\p-2 Vv) + p(x)j2(t)v.

We will prove a condition of localization of the problem (24),(2).

0

Introduce the notation

ï(t)=exp Yi(s)ds, v(t) = exp Y2(s)ds,

J0 J0

i\ w 1 I IP1 p- (n +l)

nm) = — \r\p1, pi =-,

pi p

\x\

\

Vxi, t(t)= f vmi-\n)up-2{n)dn = f um2-i(n)vp-2 u+ (t,x) = u(t)№,

v+(t,x) = v(t>KO,

e = \x\/[t (t)]i/p,

(25)

where

f(e) = A (a - )+ , m = B (a - )+2

Y =

p

p — 1' q (p — 2)2 — (mi — l)(m,2 — 1)' (p — 2)2 — (mi — l)(m2 — 1)' a > 0, A and B are positive constants.

We notice that functions u+(t,x),v+ (t,x) are Zeldovich-Barenblatt type solutions to system (24). Here it constructed for the system (24) at the first. It is hold the following

(p - 1)(p - (mi + 1)),

qi =

(p - 1)(p - (mi + 1))

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

Theorem 3. Let p > n +1, N > 1, exp ^J Yi(p)dpj < i = 1, 2, t(t) < for yt > 0 and

uo(x) ^ u+(0,x), vo(x) ^ v+(0,x), x £ RN. Then the solution of the system (1) is specially localized.

Proof. In order to prove Theorem 3 we use the following change of variables in system (1)

u(t, x) = u(t)w (t(t), y>( \x\)), v(t,x) = v(t)z (t(t),V (\x\)). Then the system reduces to the form

j4 (w, z)

dw

'JT

dz ~~3t

^s-izmi-i dw P-2 dw

d<p d<p d<p

-A. vs-iwm2-i dz P-2 dz

d<p d<p d<p

-- 0, 0,

(26)

N — l

where s = p---p > n + l.

p — (n + l)

Set in (26)

w(t,x)= f (0> e = f\x] z(t,x) = ^ = t 1/p '

Then we have from (26) the self similar system equation

Ls(f^) = ei-sd (es-irii-i

(27)

df p- ~2 df\ + if pde

de de )

d^ p- 2 d'^> ' + pde

de ~d£ )

w, z

0

Notice the functions f (£), tp(£) have properties

0 < f(Ç), ç'-i'mi-i

0 < m, ç'-1 :

'-i vm2-i

df p- -2 df

dÇ dÇ

d'tp p- 2 dtp

f = -Ap-iBmi-i(AYi)p-iÇ'f (Ç) G C(0, œ),

pi

2i

)p-'ç

To complete proof of the Theorem 3 it is sufficient to check conditions

L3( f) < 0, LA(f,') < 0, in \ Ç\ <

-,(p-i)/p

according to a comparison principle [11]. In fact, since

ÇN- i ■ mi - i

çn- i f ?2 -i

-i df \\ p -2 df

dÇ dÇ

d'tp p- 2 d'tp

- i \\ dÇ dÇ

^B^AYiy^e f (Ç),

l(AY2)p-iÇN '(Ç),

(Yi - 1)(p - 1) + (mi - 1)y2 = Yi, ((m2 - 1)yi + (p - 1)Y2 = Y2, we have

d

Çi-N — Ç

C1-

Therefore

Nd {c

N-i'mi-i -i :cm2 — 1

if m

df \\ p -2 df

dÇ dÇ

d'tp p- 2 d'

dÇ d

p- i mi - i

(AYi)p-iN + Çf),

Am2-i(AY2)p-i (N + çdVy

N

N

Ls(f , t ) =--f < 0, L4(f, ■) =--■ < 0, if \Ç\ < a(p-i)/p

pp if numbers A, B are the roots of the system of algebraic equations

AP-iBmi-i(\Yi)p-i = 1/p, Bp-i Am2-i (\j2)p-i = 1/p.

Hence

L3(f,$) < 0, L4(f,$) < 0,

in D = (t,x) : t> 0, \ x \ < l(t), l(t) = a(p-i)/p[T(t)]i/(p-l-n).

Applying comparison principle [11] we have for solution of the problem (24), (2)

u(t,x) ^ u+(t,x), v(t,x) ^ v+(t,x), in Q.

We notice that the functions w(t, x), z(t, x) have the properties

j(t,x) = z(t,x) = 0, when \x\ ^

-,(p-i)/p

[t (t)]

i/(p-l-n)

In this case for the free boundary we have an estimate

\x(t)\ < a(p-i)/p[T(t)]i/(p-l-n) for t > 0.

Therefore if t(t) < to for all t > 0, then there place a space localization of a solution of the system (1) according condition of the theorem (3). Theorem 3 is proved. □

We notice that in the case

w(t,x) = z(t, x) = 0, when \ x \ > a(p-i)/p[T (t)]i/(p-l-n).

and t(t) ^ to when t ^ to there is a finite speed of perturbation.

i

i

5. Results of numerical experiments and visualization of solutions.

For numerical solution of the problem the equation in two dimensions was used method Alternating Direction Implicit (ADI) in a combination to the method of balance. Iterative processes were constructed on the method, with Picard, Newton and a special method. In special linearization the members u3l, v32, in system (1) are presented as u331 = u^-^ uk, v3 = v32-^ v k, where u0, v0, the solutions of the system of ordinary differential equation

du

0 ßi _ _ n I ! 1

dv

0 - „ß2

dt

dt

Results of computational experiments show, that all listed iterative methods are effective for solving nonlinear problems and lead to nonlinear effects if we use as initial approximation the solutions of self-similar equations constructed by the method of nonlinear splitting and by the method of standard equation [1,5,6]. As it was expected, for obtaining some fixed accuracy the method of Newton requires smaller number of iterations, than the methods of Picard and the special method due to a good choice of an initial approximation. We observe that in each of the considered cases Newton's method has the best convergence due to the good initial approximation.

Below numerical results for one of solutions to problem (1)-(2) in the two-dimensional cases (Fig. 1). The results of the numerical experiment gives the effect of a finite speed of perturbation of a solution, and localization of a solution depending on the values of numerical parameters. The computational experiments were carried out for a slow and a quick diffusion cases.

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

parameters ? = 1 f = 20

mx = 0.2 , m2 =0.2, p = 2.8

ft = 5, /, = 2 'Ol

ß, = 7, /.= 3 1 A11

AA>1 mi + ¿>-3 = 0 /Ml - JM

eps =10~3 -

0

0

mx =1.4 , m1 =1.4. p = 2.5 -at

ß1 = 5, /[ = 2

A =5, /,= 3 .'J ft ...

AÄ >i III à -•<

mi + p-3 > 0 M fl '■1

eps =10~3 Übtet'

Fig. 1.

References

[1] M.Aripov, Sh.A.Sadullaeva. To solutions of one nondivergent type parabolic equation with double nonlinearity, Advances and Progress in Analysis, 2010, 1-18.

[2] M.Aripov, Sh.A Sadullaeva. To properties of the equation of reaction diffusion with double nonlinearity and distributed parameters, Journal of Sibirian Federal University. Mathematics & Physics, 6(2013), 157-167.

[3] M.Aripov, Approximate self-similar approach for solving of quasilinear parabolic equation, Experimentation, Modeling and Computation in Flow, Turbulence and Combustion. Wiley, 2(1997), 19-26.

[4] M.Aripov. Asymptotic of the solution of the non-Newton polytropical filtration equation, Zeitschrift fur Angewandte Mathematik und Mechanik, 80(2000), no. 3, 767-768.

[5] M.Aripov, J.Muhammadiev. Asymptotic behavior of automodel solutions for one system of quasilinear equations of parabolic type, Buletin Stiintific-Universitatea din Pitesti, Seria Matematica i Informatica, 3(1999), 19-40.

[6] M.Aripov. Method of the standard equation for the solution of the nonlinear value problem, Tashkent, Fan, 1988.

[7] R.Kersner, J.Reyes, A.Tesei. One a class of nonlinear parabolic equations with variable density and absorption advances, Diff. Equations, 7(2002), 15-176.

[8] I.Kombe, Double nonlinear parabolic equations with singular lower order term, Nonlinear Analysis, 56(2004), 185-199.

[9] S.N.Dimova, M.S.Kastchiev, M.G.Koleva, D.P.Vasileva. Numerical analysis of the blow-up regimes of combustion of two-component nonlinear heat-conducting medium, Journal of applied mathematics and mathematical physics, 35(1995), no. 3, 303-319.

[10] B.H.Gilding, L.A.Pelletier. On a class of similarity of the porous media equation. 2, J. Math. Anal. and Appl, 57(1977), 52-538.

[11] A.A.Samarskii, V.A Galaktionov, S.P.Kurdyumov, A.P.Mikhajlov, Blow-up in quasilinear parabolic equations, Berlin, 4, Walter de Grueter, 1995.

[12] S.P.Kurdyumov, E.S.Kurkina, O.V.Telkovskaya, Regimes with sharpening in two-component media, Matem. Mod., 1(1989), no. 1, 34-50.

[13] B.P.Demidovich. Lectures on mathematical theory of stability, Moscow, Nauka, 1967 (in Russian).

Численное исследование решений системы реакции-диффузии с переменной плотностью

Шахло A. Садуллаева

В статье мы покажем возможности автомодельного и приближенно-автомодельного подхода к исследованию взаимной системы реакции-диффузии нелинейных уравнений. Изучаются асимптотические поведения компактно распространенных решений и свободные границы. На основе полученных качественных свойств проведены численные расчеты. При этом решения представляются в визуальной форме по времени.

Ключевые слова: двойная нелинейность, система реакции-диффузии, автомодельное уравнение, асимптотика решений, численное решение.

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