# Nonlinear Waveforms for Ion-Acoustic Waves in Weakly Relativistic Plasma of Warm Ion-Fluid and Isothermal Electrons

Authors: S. A. El-Wakil; Essam M. Abulwafa; E. K. El-Shewy; H. G. Abdelwahed; Hamdi M. Abd-El-Hamid
Journal: Advances in Mathematical Physics
 (2012)

---

## Abstract

The reductive perturbation method has been employed to derive the Korteweg-de Vries (KdV) equation for small- but finite-amplitude electrostatic ion-acoustic waves in weakly relativistic plasma consisting of warm ions and isothermal electrons. An algebraic method with computerized symbolic computation is applied in obtaining a series of exact solutions of the KdV equation. Numerical studies have been made using plasma parameters which reveal different solutions, that is, bell-shaped solitary pulses, rational pulses, and solutions with singularity at finite points, which called “blowup” solutions in addition to the propagation of an explosive pulses. The weakly relativistic effect is found to significantly change the basic properties (namely, the amplitude and the width) of the ion-acoustic waves. The result of the present investigation may be applicable to some plasma environments, such as ionosphere region.


---

## Body

## 1. Introduction

Nonlinear evolution equations are widely used as models to describe complex physical phenomena and have a significant role in several scientific and engineering fields [1, 2]. The propagation of solitary waves is important as it describes characteristic nature of the interaction of the waves and the plasmas. In the case where the velocity of particles is much smaller than that of light, ion-acoustic waves present the nonrelativistic behaviors, but in the case where the velocity of particles approaches that of light, the relativistic effect becomes dominant [3]. Actually high-speed and energetic streaming ions with the energy from 0.1 to 100 MeV are frequently observed in solar atmosphere and interplanetary space. Nevertheless, relativistic ion-acoustic waves have not been well investigated. When we assume that the ion energy depends only on the kinetic energy, such plasma ions have to attain very high velocity of relativistic order. Thus, by considering the weakly relativistic effect where the ion velocity is about 1/10 of the velocity of light, we can describe the relativistic motion of such ions in the study of nonlinear interaction of the waves and the plasmas [4]. It appears that the weakly relativistic and ion temperature effects play an important role in energetic ion-acoustic waves propagating in interplanetary space [5, 6].Washimi and Taniuti [7] were the first to use reductive perturbation method to study the propagation of a slow modulation of quasi-monochromatic waves through plasma. And then the attention has been focused by many authors [8–11].The evolution of small-but finite-amplitude solitary waves, studied by means of the Korteweg-de Vries (KdV) equation, is of considerable interest in plasma dynamics.Many powerful methods have been established and developed to study nonlinear evolution equations (NLEEs). These methods include the inverse scatting method [12], the tanh function method [13], extended tanh method [14–16], the Exp-function method [17], the extended F-expansion method [18], the Jacobi elliptic function expansion method [19], the homogeneous balance method [20], sech-function method [21], F-expansion method [22], and the multiple exp-function method [23]. There is no unified method that can be used to deal with all types of NLEEs. Fan [24] developed a new algebraic method with computerized symbolic computation, which greatly exceeds the applicability of the existing tanh, extended tanh methods, and Jacobi function expansion method in obtaining a series of exact solutions of nonlinear differential equations.Recently, the ion-acoustic solitary wave in collisionless unmagnetized plasma consisting of warm ions fluid and isothermal electrons is studied using the time-fractional KdV equation by El-Wakil et al. [25]. They showed that the time fractional can be used to modulate the electrostatic potential wave.The major topic of this work is to study the ion-acoustic solitary and other type waves in relativistic warm plasma. This paper is organized as follows. In Section2, we present the basic set of fluid equations governing our plasma model. In Section 3, an algorithm describing the computerized symbolic computation method is presented. In Section 4, explicit solutions for KdV equation are obtained. Finally, some discussions and conclusions are given in Section 5.


## 2. Basic Equations and KdV Equation

Consider collisionless ionization-free unmagnetized plasma consisting of a mixture of warm ion-fluid and isothermal electrons. Assume that the ion flow velocity has a weak relativistic effect, and therefore there exist streaming ions in an equilibrium state when sufficiently small- but finite-amplitude waves propagate one-dimensionally. Such a system is governed by the following normalized equations [26]:(2.1a)∂∂tn(x,t)+∂∂x[n(x,t)u(x,t)]=0,(2.1b)[∂∂t+u(x,t)∂∂x][γ(x,t)u(x,t)]+σn(x,t)∂∂xp(x,t)+∂∂xϕ(x,t)=0,(2.1c)[∂∂t+u(x,t)∂∂x]p(x,t)+3p(x,t)∂∂x[γ(x,t)u(x,t)]=0,(2.1d)∂2∂x2ϕ(x,t)+n(x,t)-ne(x,t)=0.The electron temperatureTe is much larger than the ion temperature Ti and in this case for simplicity one can neglect the inertia of the electrons relative to that of the ions, that is, the high-frequency plasma oscillations are neglected. Since it is interested with the regime of density and velocity fluctuations near the ion plasma frequency, so the isothermal electrons density is given by
(2.1e)ne(x,t)=exp⁡[ϕ(x,t)].For weakly relativistic effects, the relativistic factorγ(x,t)=1/1-u2(x,t)/c2 is approximated by
(2.1f)γ(x,t)≈1+u2(x,t)(2c2),where c is the velocity of light.In (2.1a)–(2.1f), n(x,t) and ne(x,t) are the densities of ions and electrons, respectively, u(x,t) is the ion flow velocity, p(x,t) is the ion pressure, ϕ(x,t) is the electric potential, x is the space coordinate, and t is the time variable. σ=Ti/Te≪1 is the ratio of the ion temperature to the electron temperature. All of these quantities are dimensionless and are normalized in terms of the following characteristic quantities: n(x,t) and ne(x,t) by the unperturbed electron density n0, u(x,t) and c by the sound velocity kBTe/mi, p(x,t) and ϕ(x,t) by n0kBTi and kBTe/e, respectively; t and x by the inverse of the plasma frequency ωpi-1=1/4πe2n0/mi and the electron Debye length λD=kBTe/(4πe2n0), respectively, kB is Boltzmann’s constant, and mi is the mass of plasma ion.According to the general method of reductive perturbation theory, the stretched variables are introduced as [7]
(2.2)τ=ε3/2t,ξ=ε1/2(x-λt),
where λ is the phase velocity and ε represents the amplitude of the perturbation. All the physical quantities that appeared in (2.1a)–(2.1f) are expanded as power series in ε about the equilibrium values as(2.3a)n(ξ,τ)=1+εn1(ξ,τ)+ε2n2(ξ,τ)+⋯,(2.3b)u(ξ,τ)=u0+εu1(ξ,τ)+ε2u2(ξ,τ)+⋯,(2.3c)p(ξ,τ)=1+εp1(ξ,τ)+ε2p2(ξ,τ)+⋯,(2.3d)ϕ(ξ,τ)=εϕ1(ξ,τ)+ε2ϕ2(ξ,τ)+⋯.The boundary conditions of this problem are imposed as |ξ|→∞, n=ne=p=1, u=u0, and ϕ=0.Substituting (2.2) and (2.3a)–(2.3d) into the system of (2.1a)–(2.1f) and equating the coefficients of like powers of ε, then from the lowest order(2.4a)n1(ξ,τ)=ϕ1(ξ,τ),(2.4b)u1(ξ,τ)=(λ-u0)ϕ1(ξ,τ),(2.4c)p1(ξ,τ)=3γ1ϕ1(ξ,τ),
with the transcendental equation of λ that is given as
(2.4d)1-(λ2-3σ)γ1=0⟹λ=±(3σγ1+1)γ1,
where
(2.4e)γ1=1+3γ022,γ0=u0c.Using second-order equations in ε and eliminating the second-order perturbed quantities n2, u2, p2, and ϕ2, the following KdV equation for the first-order perturbed potential is obtained: (2.5a)∂∂τϕ1(ξ,τ)+Aϕ1(ξ,τ)∂∂ξϕ1(ξ,τ)+B∂3∂ξ3ϕ1(ξ,τ)=0.
Therefore, the nonlinear coefficient A and the dispersion coefficient B are represented by
(2.5b)A=B(g1-g2γ2γ13),B=1g2γ1,
with
(2.5c)g1=3σγ1(3γ1+1)+2,g2=23σγ1+1,γ2=3γ0(2c).


## 3. Computerized Symbolic Computation Method

An algebraic method with computerized symbolic computation has been developed by Fan [24], which can be used to solve a given partial differential equation in ϕ1(ξ,τ) of the form(3.1a)H(ϕ1,∂ϕ1∂τ,∂ϕ1∂ξ,∂2ϕ1∂ξ2,…)=0.This equation may be transformed into an ordinary differential equation of the form
(3.1b)H(Φ,dΦdη,d2Φdη2,…)=0
using a traveling frame of reference
(3.1c)ϕ1(ξ,τ)=Φ(η),η=ξ-vτ,where v is the traveling wave propagation velocity.The computational technique used to solve the NLEE is described as follows [24].Step 1.
Reduce partial differential equation (3.1a) to the ordinary differential equation (3.1b) by considering the traveling wave transformation (3.1c).Step 2.
Expand the solution of (3.1b) in the form(3.2a)Φ(η)=∑i=0naiφi(η),
where the new variable φ(η) is a solution of the following ordinary differential equation:
(3.2b)dφ(η)dη=∈∑j=0rcjφj(η),∈=±1.Step 3.
Substituting (3.2b) into (3.1b) and balancing the highest derivative term with the highest nonlinear term lead to a relation between n and r, from which the different possible values of n and r can be obtained. These values lead to the different series expansions of the solutions.Step 4.
Substituting the expansions (3.2a)–(3.2b) into (3.1b) and setting the coefficients of all powers of φi and φidφ/dη to zero will give a system of algebraic equations, from which the parameters ai  (i=0,1,…,n) and cj  (j=0,1,…,r) can be found explicitly.Step 5.
Substituting the parameterscj  (j=0,1,…,r) obtained in Step 4 into (3.2b) gives all the possible solutions φ(η).
It is remarked here that the solutions of (3.1a) depend on the explicit solvability of (3.2b). The solutions of (3.2b) will get a series of fundamental solutions such as polynomial, exponential, soliton, rational, triangular periodic, Jacobi, and Weierstrass elliptic doubly periodic solutions.


## 4. Explicit Solutions for the KdV Equation

For KdV equation (2.5a), the traveling wave transformation (3.1c) leads to
(4.1)-vddηΦ(η)+AΦ(η)ddηΦ(η)+Bd3dη3Φ(η)=0.Balancing the highest derivative term with the highest nonlinear term leads to a relation betweenn and r as r=n+2. Taking n=2 gives r=4 and leads to(4.2a)Φ(η)=a0+a1φ(η)+a2φ2(η),(4.2b)dφdη=∈c0+c1φ(η)+c2φ2(η)+c3φ3(η)+c4φ4(η).Substituting (4.2a)–(4.2b) into (4.1) and equating coefficients of all powers of φi and φidφ/dη to zero will get a system of algebraic equations, from which the parameters ai  (i=0,1,2) and cj  (j=0,1,…,4) can be found explicitly. Substituting cj into (4.2b) and using the symbolic software package Maple give explicit solutions of (4.2b).Substituting the coefficientsai and (4.2a) into (4.1) and using the symbolic software package Maple, we obtain, for KdV equation (4.1), the following solutions:
(4.3)Φ(η)=a0-3(Aa0-v)Asech2(12(v-Aa0)Bη),c0=c1=c4=0,(v-Aa0)B>0,(4.4)Φ(η)=a0-3(Aa0-v)Asec2(12(Aa0-v)Bη),c0=c1=c4=0,(v-Aa0)B<0,(4.5)Φ(η)=a0+32(Aa0-v)Atan2(12(v-Aa0)2Bη),c0≠0,c1=c3=0,(v-Aa0)B>0,(4.6)Φ(η)=a0-32(Aa0-v)Atanh2(12(Aa0-v)2Bη),c0≠0,c1=c3=0,(v-Aa0)B<0,(4.7)Φ(η)=a0+3(Aa0-v)Acsch2(12(v-Aa0)Bη),c0=c1=c3=0,(v-Aa0)B>0,(4.8)Φ(η)=a0-3(Aa0-v)Acsc2(12(Aa0-v)Bη),c0=c1=c3=0,(v-Aa0)B<0,(4.9)Φ(η)=vA-12BAη2,c0=c1=c3=0,(v-Aa0)B=0,(4.10)Φ(η)=a0-3(Aa0-v)Am2(2m2-1)cn2(12Aa0-vB(1-2m2)η),c0≠0,c1=c3=0,v-a0AB>0,(4.11)Φ(η)=a0+3(Aa0-v)A1(m2-2)dn2(12Aa0-vB(m2-2)η),c0≠0,c1=c3=0,v-a0AB>0,(4.12)Φ(η)=a0-3(Aa0-v)Am2(m2+1)sn2(12Aa0-vB(m2+1)η),c0≠0,c1=c3=0,v-a0AB<0(4.13)Φ(η)=vA-34c323BA℘(c343η;-4c33c1,-c0),c2=c4=0,c3>0.
where m is the modulus of the Jacobi elliptic functions, a0, c0, c1, and c3 are arbitrary constants, and  ℘ is the Weierstrass elliptic doubly periodic function.


## 5. Results and Discussion

Solutions (4.3) and (4.6) are hyperbolic wave solutions. In Figure 1, a profile of the bell-shaped solitary pulse is obtained for solution (4.3). Figures (1(a) and 1(b)) show that the soliton amplitude and width are sensitive to the relativistic factor γ0 and the temperatures ratio σ. Also, the arbitrary value a0 gives the same effect of adding a higher-order perturbation correction in increasing the amplitude and decreasing the width as shown in Figure 1(c).Profile of localized pulses for expression (4.3).
(a)
(b)
(c)Solutions (4.4) and (4.5) are triangular solutions that develop solitons with singularity at finite points, which are called “blowup” solutions [27] as in Figure 2(a). In Figure 2(b), for some values of a0 the pulse does not vanish to infinity, so it has a localized form. Therefore, the arbitrary value a0 plays a role in obtaining a pulse-shaped localized solution [28]. Solutions (4.7) and (4.8) lead to the propagation of an explosive pulses [29]. In Figure 3, the profiles of explosive (divergent) pulses are depicted for expression (4.7). For a0=0, the solution gives a localized explosive pulse, as in Figure (2.3a) while a periodic solution is obtained for a0=0.4, as shown in Figure (2.3b). The rational solution (4.9) may be helpful to explain certain physical phenomena. Because a rational solution is a disjoint union of manifolds, particle systems describing the motion of a pole of rational solutions for a KdV equation were analyzed [30]. Equations (4.10)–(4.12) are three Jacobi elliptic doubly periodic wave solutions. When m→1, solutions (4.10) and (4.11) reduce to (4.3) while (4.11) reduces to (4.6) [31]. In Figure 4, a profile of triangular periodic wave solution for expression (4.10) for m=2 is shown. On the other hand, (4.13) gives the Weierstrass elliptic doubly periodic type solution as depicted in Figure 5.Profile of pulses for expression (4.5) for v=0.04, σ=0.001, and γ0=0.01.
(a)
(b)The explosive solution for expression (4.7) for v=0.04, σ=0.001, and γ0=0.01.
(a)
(b)Figure 4
The asymptotical property of Jacobi elliptic doubly periodic wave for expression (4.10) with v=0.04, σ=0.001, a0=0, and m=2.Figure 5
Weierstrass elliptic doubly periodic type solution (4.13) for v=0.04, γ0=0.01, c0=-1, c1=1, and c3=0.1.In summary, it has been found that the amplitude and the width of the ion-acoustic waves as well as parametric regime where the solitons can exist are sensitive to the relativistic factorγ0 and the ratio of the ion to the electron temperatures σ. Moreover, solutions for KdV equation have been obtained. It may be important to explain some physical phenomena in some plasma environments, such as ionosphere region.



---

*Source: 529121-2012-07-15.xml*