A NUMERICAL INVESTIGATION OF MODIFIED BURGERS’ EQUATION IN DUSTY PLASMAS WITH NON-THERMAL IONS AND TRAPPED ELECTRONS Текст научной статьи по специальности «Физика»

NANOSYSTEMS:

PHYSICS, CHEMISTRY, MATHEMATICS

Deka H., Sarma J. Nanosystems: Phys. Chem. Math., 2023,14 (1), 5-12.

http://nanojournal.ifmo.ru

Original article

DOI 10.17586/2220-8054-2023-14-1-5-12

A numerical investigation of modified Burgers' equation in dusty plasmas with non-thermal ions and trapped electrons

Harekrishna Deka1'", Jnanjyoti Sarma2'6

XK. K. HandiquiState Open University, Department of Mathematics, Khanapara, Guwahati, 781022, India 2R.G. Baruah College, Fatasil Ambari, Guwahati, 781025, India

aharekrishnadeka11@gmail.com, bsjnanjyoti@gmail.com Corresponding author: Harekrishna Deka, harekrishnadeka11@gmail.com

PACS 02.70Bf, 52.27Lw, 52.35Fp, 52.35Tc

Abstract In this paper, one-dimensional lower order modified Burgers' equation (MBE) in dusty plasmas having non-thermal ions and trapped electrons is investigated numerically by finite difference explicit method. The numerical results obtained by the finite difference explicit method for various values of the nonlinear and dissipative coefficients have been compared with the analytical solutions. The obtained numerical results are found to have good agreement with the analytical solutions. It is found that the nonlinear and dissipative coefficients have very important effect on the dust acoustic waves in the system. The absolute error between the analytical and the numerical solutions of the MBE is demonstrated. The stability condition is derived in terms of the equation parameters and the discretization using the von Neumann stability analysis. It has been observed that the waves become flatten and steeper when the dissipative coefficient decreases. It can be concluded that the finite difference explicit method is suitable and efficient method for solving the modified Burgers' equation.

Keywords plasma, dusty plasmas, non-thermal ions, reductive perturbation method, modified Burgers' equation, finite difference explicit method, von Neumann stability analysis

For citation Deka H., Sarma J. A numerical investigation of modified Burgers' equation in dusty plasmas with non-thermal ions and trapped electrons. Nanosystems: Phys. Chem. Math., 2023,14 (1), 5-12.

1. Introduction

In the last few decades, many researchers studied the electrostatic and electromagnetic waves propagation in dusty plasmas in various environments such as in the upper part of the Earth atmosphere, planetary rings, comet tails, interstellar space, the solar atmosphere and low temperature plasmas in laboratory [1-3]. Many authors [4,5] investigated the effect of higher order nonlinearity for dusty plasma considering the negative ions and hot isothermal as well as non-isothermal electrons. Asgari et al. [6] derived the nonlinear Burgers' equation with a non-thermal ion in dusty plasma environment. Many authors [7-9] also investigated the nonlinear behaviors of electrostatic waves in a dusty plasma with trapped particles as well as in unmagnetized and magnetized plasmas. Dev et al. [10] have investigated the wave propagation in a non-magnetized and warm dusty plasma containing trapped electrons as well as non-thermal positive and negative ions under the influence of lower order nonlinearity. In this paper, we deal with the nonlinear modified Burgers' equation in dusty plasmas having negative and positive non-thermal ions with trapped electrons. The finite difference method was first developed by Thomas in 1920 to solve nonlinear hydrodynamic equations [11]. Finite difference methods are the first techniques for numerical solving of nonlinear partial differential equations [12]. The most commonly used finite difference methods for the solution of partial differential equations are as follows: Explicit method, Implicit method and Crank Nicolson method. Many authors applied finite difference explicit method to obtain numerical solutions of nonlinear partial differential equations. The modified Burgers' equation is a nonlinear expansion of the Burgers' equation. There are several methods proposed for solving the modified Burgers' equation that can be briefly presented as follows.

Bratsos et al. [13] applied a finite difference scheme for calculating the numerical solution of the modified Burgers' equation. A collocation method based on quantic splines was proposed by Ramadan and El-Danaf [14]. Irk [15] also proposed the sextic B-spline collocation method for numerical solution of the modified Burgers' equation. Aswin and Awasthi [16] have solved the modified Burgers' equation using iterative differential quadrature algorithms. Roshan and Bharma [17] solved the modified Burgers' equation by the Petrov-Galerkin method. Grienwank and El-Danaf [18] have employed a non-polynomial spline based method to obtain numerical solutions of the modified Burgers' equation. Duan et al. [19] proposed Lattice Boltzmann method to obtain numerical solution of modified Burgers' equation. Bashan et al. [20] used quintic B-spline Differential Quadrature method to solve modified Burgers' equation. Ucar et al. [21] proposed the finite difference method for numerical solution of the modified Burgers' equation. Karakoc et al. [22] have proposed the

© Deka H., Sarma J., 2023

quartic B-spline subdomain finite element method (SFEM) for finding the numerical solution of the Burgers' equation and modified Burgers' equation. In the paper, we have used the finite difference explicit method to obtain the numerical solution of the modified Burgers' equation in dusty plasmas having non-thermal ions and trapped electrons. The finite difference approximation to partial derivatives can be obtained from Taylor series expansion using either the backward, forward or central difference approximations. The paper is organized in the following manner. In section 2, we discuss the basic equations governing the proposed dusty plasma system, modified Burgers' equation with trapped electrons and reductive perturbation technique. In section 3, the stability analysis of the finite difference explicit method is presented. We include the solutions and discussion of numerical results in section 4. Conclusion is given in section 5.

2. Basic equations and Modified Burgers' Equation

The basic equations, governing the dynamical system, of the dust particles in a one-dimensional dust acoustic wave for the dusty plasma is as follows [10]:

TF + I 'n") = ", <■>

1

dvd + v dvd +___

dt dx md'nd dx

dpd _ _qd_ d^ + d2Vd md dx ^ dx2 '

dpd dpd dvd -^r + vd^~ + YPd-K-dt dx dx

0,

d2 ^ A .

_ 4ne (ne + nn - Up + ZdUd),

(2)

(3)

(4)

where nd is the number density of the negatively charged dust particles in the plasma, vd is the dust fluid velocity, ^ is the electrostatic potential, zd is the number of electrons residing on the dust surface at equilibrium, pd is the pressure of the dust fluid, e is the electron charge and md is the mass of the dust particle and y = 3 is the adiabatic index. The non-thermal number densities of positive ions np and negative ions nn can be expressed by the using the following relations [23]:

np = np0 (1 + + a<f>2^ exp (—p$), (5)

(6)

Uno 11 + aap0 + a (ap)2 j exp (zri<7p4>)

with

4yi T -, where 71 represents the population of the non-thermal ions, and ap = —p. Also, zp(zn) is the positive

1 + 371 F Tn

(negative) ion's charge state, Te is the electron temperature and Tp'Tn) is the positive ion (negative ion) temperatures. The electron density in the presence of trapped electron can be expressed by using the following relations [7,24] :

Ueo j 1 + (ßft - b (ß0)3/2 + 2 (ß0)2 - ... J

(7)

T 4(1 — 72) Tef

where ft = -f, b = —— and the parameter y2 is as follows y2 = , with Tef and Tet being the temperatures of

Te 3yj n let

free electrons and trapped electrons in plasma, respectively.

The parameter y2 determines the nature of the distribution function, giving a plateau if y2 = 0 and a dip if y2 < 0 and a hump shape if y2 > 0. However, y2 = 1 corresponds to the Maxwellian distribution of the electrons. In the present plasma system, the range of y2 will be considered as 0 < y2 < 1 for non-isothermal (trapped) electrons.

Now, Nd dust number density is normalized to its equilibrium value nd0, Vd dust-fluid velocity is normalized to

ks T

-p, where kB is the Boltzmann constant, the time

Csd _

iZdkß Tp\ 1/2

md

variable T is normalized to w.

, $ is the DA wave's potential normalized to

/ ™ , N1/2

1

md

, the space variable X is normalized to XD1d

\4nndoz'2de2

and pressure pd is normalized to pd = ndokBTd.

The normalized forms of the basic equations (1 - 4) are as follows

w- + rn 'NdVd) = 0,

„ dVd „ t, dVd dPd „r d$ d2Vd Nd + NdVd dX + dx = Nd dX + VdXd' Pd Pd Vd

-di + Vd dx^ + 3pd dx' = °,

-pi$ — p2$3/2 + p3$2 — p4$5/2 + (Nd — 1)

^Anndoz2de2^ 1/2

kß Tp

d20 dX 2

with the overall charge neutrality condition

neo _ npo - Zdndo - nno

(8)

(9) (10) (11)

(12)

n

e

e

and

Mp =

Про

Zd'do '

Td

ad =

ZdTp '

Mn =

ппО

ПеО

Zd'ndo' Zd ndo

= Mp - Mn - i, n =

P3

(Mp - Mn - ^77 h Mn

P2 = b (Mp - Mn - i) в3/2,

/ (Zn)2 h h ( ^

—--h aaZn h a (a) - Mp

(zP)2

M

Pl = {(Mp - Mn - i) в h Mn (Zn h aa) - Mn (a - Zp)} ,

- aZp h a

and

P4 = (Mp - Mn - 1) b^5/2.

The modified Burgers' equation for the propagation of small and finite amplitude dust-acoustic shock waves (DASWs) is derived. Here the independent variables £ and t are the new stretched variables given by £ = e1/2 (x - V0t) and t = e t, where V0 is the phase speed (normalized to Csd) and e is a small parameter (0 < e < 1) which measures the weakness of the dispersion. The dependent variables Nd, Vd, Pd and ^ can be expanded in the power series as follows [10]:

Nd = 1 + eN((1) + e3/2N2) + ..., (13)

Vd = eVd(1) + e3/2Vd(2) + ..., (14)

Pd = 1 + ePd(1) + e3/2P^2) + .. ., (15)

ф = eф(1) h e3/2ф(2) h e2ф(3) h ....

(16)

Substituting the stretched coordinates and the expression for Nd, Vd, Pd and ^ into the normalized basic equations (8 -11) and equating the coefficients of lower order in e, the required lower order modified Burgers' equation is obtained

дф(1)

h A (ф(1))

1/2 дф(1) д2ф(1)

B

дт ' " V ) д£ д£2 '

where the nonlinear coefficient A and dissipative coefficient B are given by the following expressions

3P2 (V2 - 3ad) B = n 1 2 '

A=

(17)

(1S)

4p1Vo

Equation (17) represents the well-known lower order modified Burgers' equation describing the nonlinear propagation of dust-acoustic shock waves in electronegative dusty plasma with non-thermal ions and trapped electrons.

The stationary shock wave solution of the modified Burgers' equation (17) is obtained by transforming independent variables £ and t to £ = Z - U0t' and t = t' where U0 is the speed of the shock waves. Now, the analytical solution of the modified Burgers' equation is given by the formula

ф(1) = { фт1 ^ i - tan^ j^J J

(19)

3M

4B

where = —— and = —— are the amplitude and the width of the shock waves, respectively, and M is the Mach

4A MW

number.

3. Stability analysis of the explicit finite difference method

In this section, the stability of the finite difference explicit method is investigated by using von Neumann stability analysis. The von Neumann stability theory in which the growth factor of a Fourier mode is defined as = £jeifchi = £jwhere I = %/—I, £j is the amplitude at time level k is the wave number and h = Ax.

To investigate the stability of the numerical scheme, the non-linear term u1/2u of the modified Burgers' equation has been linearized by putting u1/2 = L, where L is constant.

From (29):

Ui,j+1 = U®,j h

kAL

Xi,j+1

i

2h 2kB \

ИГ) U

l,j - ui+l,j

kB + "h2"

ui+l,j - 2ui,j h ui_l,j

h

kAL kB

2h

Substitute u1/2 = L in the above equation, we get

+ h2

Ui- lj h

kB kAL\

h2

2h

Ui+l,

£j e^e = £j e

I0í

2kB kAL kB

i - H h -al h -B

h2

e_" h

kB ^h2

-А-Л ^в'

2h

^ . 2kB\ / -AL kB\ i0 ZkB -Am i0

£ =(i - ""И H -2h- h ^J e_Iв H Ж - ,

2

u

* =I1 - HT) +

2kB \

*

2kB\ , kAL ( _

- I 1

2h V W) + ^2h~K

e ie - eiej +

kB ~hh

eI6 + e

-ie

1

kAL ( T \ kB

- 2I sin 9 J 2 cos 9,

2kB kAL 2kB

* =U - -hr) - —sin 9 + hBcos 9

The stability condition for the numerical scheme is as follows

1*1 < 1,

|*|

2kB

1 -

kAL 2kB

sin 9 +—r^ cos 9

h

h2

1.

2kB h2

So the stability condition is -—r- < 1 or k <7-^. And thus restricts the allowable temporal step size. Since not all

h2 2B

choices of spatial and temporal steps lead to convergent, the explicit scheme (29) is called conditionally stable.

4. Numerical results and discussion

It has been shown in Section 2 that the modified Burgers' equation has an analytical solution in the form (19). In this paper, the modified Burgers' equation (17) in dusty plasmas with non-thermal ions and trapped electrons is solved by the explicit finite difference method and the numerical results are compared with the analytical solutions. For simplicity, we consider (Z, t) = u(x, t) = u(iAx,j At) = ui,j. Equation (17) can be expressed as

du 1/2 du -gd^u

dt dx dx2

11 3 4B

For convenience, we put M = -, U0 = -, Sm1 = —- and â1 = —— = 8B in (19).

' F 2' 0 2 Ym1 8A 1 M

The analytical solution of the modified Burgers' equation is given by the following expression

u(x,t) = ' 8A '1 - tanh^ ' x -

(20)

(21)

The boundary and the initial conditions are taken from the exact solution. In this paper, the numerical solutions of (20) will be sought for the following initial and boundary conditions. In the interval 0 < x < 1, t > 0, and with the initial condition

3x

and the boundary conditions

dx, 0) = (1 - tanh—ll ( , ) \8A\ 8B J J

^HsAl1 + tanh^

u(1,t)H8Ai1 --1

(22)

(23)

(24)

. ,.„ , . . , . du , ________du .52u , .

We discretize the modified Burgers equation by replacing —- by the forward difference and — and —by the

dt dx dx2

central difference approximation, i.e. as follows

du

J'i,j + 1

J'i,3

dt k

d2u ^ ui_i,j - 2ui,j + u,i+i,j dx2 h2

du ^ ui+i,j - ui_i,j

dx ~ 2h .

Thus, (20) becomes as follows

i,j+1, ui,j + Auij1/2

ui+1,j - ui_1,

which can be simplified

kA

kA 1/2 ui,j+i = ui,j + ~2hui,j

2h

ui_1,j - ui+1,j

B

+

kB h2

ui+1,j - ui,j + ui_1,. h2

u,i+i,j - 2ui,j + ui_i,j

(25)

(26)

(27)

(28) (29)

2

2

u

By the von Neumann stability condition, we consider the space steps h = 0.01, the time steps k = 0.005, 0.001 and 0.0005 and B = 0.01, 0.05 and 0.1. The nonlinear coefficient A is a function of pi, p2, and V0. During the solution process, we consider A =1, 1.5 and 2. Numerical solutions are obtained for different values of B and A. The obtained figures of the numerical results are compared with the figures obtained from the analytical solutions which are displayed in Figs. 1-3.

The accuracy of the present method is measured using the absolute error which is defined as

Analytical _ Numerical

Fig. 1. Graph of analytical and numerical solutions with absolute error at A = 1, B = 0.01, h = 0.01, k = 0.005

Analytical Solution Numerical Solution

0 0 0 0.2 0.4 0.6 0.8 1

X

X

Fig. 2. Graph of analytical and numerical solutions with absolute error at A = 1.5, B = 0.05, h = 0.01, k = 0.001

Analytical Solution Numerical Solution

Absolute Error лК)"6 Absolute Error

Fig. 3. Graph of analytical and numerical solutions with absolute error at A = 2, B = 0.1, h = 0.01, k = 0.0005

The analytical solutions and the computed numerical results together with their errors are plotted in Figs. 1 to 3 for various values of the nonlinear coefficient and the dissipative coefficient. But the graphs of the errors have been drawn at time t = 5. It can be seen that the maximal error occurs at the left hand boundary when the greater value of the dissipative coefficient B = 0.1 is considered and the maximal error is found around the location where the shock wave has the highest amplitude with the smaller value of the dissipative coefficient B = 0.01. It can be concluded that as the value of x increases, the errors gradually decreases and it is also seen that as the value of the dissipation coefficient increase, the error will increase.

FIG. 4. The numerical solutions with A =1, B = 0.01, h = 0.01, k = 0.005 at t = 0.5, 1, 1.5, 2

0 0.5 1 0 0. S 1

X X

fig. 5. The numerical solutions with A =1.5, B = 0.05, h = 0.01, k = 0.001 at t = 1, 2, 3, 4

n

FIG. 6. The numerical solutions with A = 2, B =0.1, h = 0.01, k = 0.0005 at t = 1, 2, 3, 4

The computed numerical results for various values of the nonlinear coefficient and the dissipative coefficient at different times are plotted in Figs. 4 to 6. From the figures, it has been observed that as the time increases, the curve of the numerical solution decays.

5. Conclusion

In this paper, the finite difference explicit method has been successfully used for obtaining the numerical solution of the modified Burgers' equation in dusty plasmas having non-thermal ions with trapped electrons. It is obtained using finite difference explicit method. Graphs have been plotted to show a comparison between the analytical and the numerical solutions for various values of the dissipative coefficient. The obtained numerical results show good accuracy when comparing it with the analytical results for various values of the dissipative coefficient. The absolute error has been computed and presented in graphical form. The effects of the nonlinear coefficient and the dissipative coefficient on the shock strength and steepness are investigated. It is found that the shock wave steepness depends more on the dissipative coefficient than on the nonlinear coefficient. In this study, it has been seen that when the dissipative coefficient decreases, the shock waves become flatten and the propagation front becomes steeper. In conclusion, it has been observed that the dissipative coefficient in dusty plasmas having non-thermal ions and trapped electrons plays an important role to dissipate the acoustic shock wave while propagating through the system.

References

[1] Goertz C. Dusty plasmas in the solar system. Reviews of Geophysics, 1989, 27 (2), P. 271-292.

[2] Adhikary N.C., Deka M.K., Bailung H. Observation of rarefactive ion acoustic solitary waves in dusty plasma containing negative ions. Physics of Plasmas, 2009,16 (6), 063701.

[3] Deka M., Adhikary A., Misra A., Bailung H., Nakamura Y. Characteristics of ion-acoustic solitary wave in a laboratory dusty plasma under the influence of ion-beam. Physics of Plasmas, 2012,19 (10), 103704.

[4] Tagare S., Reddy R.V. Effect of higher-order nonlinearity on propagation of nonlinear ion-acoustic waves in a collisionless plasma consisting of negative ions. J. of Plasma Physics, 1986, 35 (2), P. 219-237.

[5] El-Labany S. Contribution of higher-order nonlinearity to nonlinear ion-acoustic waves in a weakly relativistic warm plasma. Part 1. Isothermal case. J. of Plasma Physics, 1993, 50 (3), P. 495-504.

[6] Asgari H., Muniandy S., Wong C. Dust-acoustic shock formation in adiabatic hot dusty plasmas with variable charge. Physics of Plasmas, 2011, 18 (1), 013702.

[7] Schamel H. A modified Korteweg-de Vries equation for ion acoustic waves due to resonant electrons. J. of Plasma Physics, 1973, 9 (3), P. 377-387.

[8] Mamun A., Cairns R., Shukla P. Effects of vortex-like and nonthermal ion distributions on non-linear dust-acoustic waves. Physics of Plasmas, 1996, 3 (7), P. 2610-2614.

[9] Javidan K., Pakzad H. Obliquely propagating electron acoustic solitons in a magnetized plasma with superthermal electrons. Indian J. of Physics,

2013, 87 (1), P. 83-87.

[10] Dev A.N., Deka M.K., Sarma J., Adhikary N.C. Shock wave solution in a hot adiabatic dusty plasma having negative and positive non-thermal ions with trapped electrons. J. of the Korean Physical Society, 2015, 67 (2), P. 339-345.

[11] Smith G.D.,Smith G.D., Smith G.D.S. Numerical solution of partial differential equations: finite difference methods. Oxford University Press, 1985.

[12] Dehghan M. Parameter determination in a partial differential equation from the overspecified data. Mathematical and Computer Modelling, 2005, 41 (2-3), P. 196-213.

[13] Bratsos A., Petrakis L. An explicit numerical scheme for the modified burgers' equation. Int. J. for Numerical Methods in Biomedical Engineering, 2011, 27 (2), P. 232-237.

[14] Ramadan M.A., El-Danaf T.S. Numerical treatment for the modified burgers equation. Mathematics and Computers in Simulation, 2005, 70 (2), P. 90-98.

[15] Irk D. Sextic B-spline collocation method for the modified Burgers' equation. Kybernetes, 2009, 38 (9), P. 1599-1620.

[16] Aswin V., Awasthi A. Iterative differential quadrature algorithms for modified burgers equation. Engineering Computations, 2018, 35 (1), P. 235250.

[17] Roshan T., BhamraK. Numerical solutions of the modified Burgers' equation by Petrov-Galerkin method. Applied Mathematics and Computation, 2011, 218 (7), P. 3673-3679.

[18] Griewank A., El-Danaf T.S. Efficient accurate numerical treatment of the modified Burgers' equation. Applicable Analysis, 2009, 88 (1), P. 75-87.

[19] Duan Y., Liu R., Jiang Y. Lattice Boltzmann model for the modified Burgers' equation. Applied Mathematics and Computation, 2008, 202 (2), P. 489-497.

[20] Bashan A., Karakoc S.B.G., Geyikli T. B-spline differential quadrature method for the modified Burgers' equation. Cankaya University J. of Science and Engineering, 2015, 12 (1), P. 1-13.

[21] Ucar A., Yagmurlu N.M., Tasbozan O. Numerical solutions of the modified Burgers' equation by finite difference methods. J. of Applied Mathematics, Statistics and Informatics, 2017, 13 (1), P. 19-30.

[22] Karakoc S.B.G., Bashan A., Geyikli T. Two different methods for numerical solution of the modified Burgers' equation. The Scientific World J.,

2014, 2014, 780269.

[23] Adhikary A., Deka M., Dev A., Sarmah J. Modified Korteweg-de Vries equation in a negative ion rich hot adiabatic dusty plasma with non-thermal ion and trapped electron. Physics of Plasmas, 2014, 21 (8), 083703.

[24] Schamel H. Stationary solitary, snoidal and sinusoidal ion acoustic waves. Plasma Physics, 1972, 14 (10), 905.

Submitted 13 September 2022; revised 10 December 2022; accepted 13 December 2022

Information about the authors:

Harekrishna Deka - K. K. Handiqui State Open University, Department of Mathematics, Khanapara, Guwahati, 781022, India; ORCID 0000-0003-4280-3728; harekrishnadeka11@gmail.com

Jnanjyoti Sarma - R. G. Baruah College, Fatasil Ambari, Guwahati, 781025, India; ORCID 0000-0002-0793-5680; sjnanjyoti@gmail.com

Conflict of interest: There is no conflict of interest in this research paper.