Open AccessArticle

Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma

Sherif M. E. Ismaeel

^1,2,*,

Abdul-Majid Wazwaz

Elsayed Tag-Eldin

⁴

and

Samir A. El-Tantawy

^5,6

Department of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia

Department of Physics, Faculty of Science, Ain Shams University, Cairo 11566, Egypt

Department of Mathematics, Saint Xavier University, Chicago, IL 60655, USA

⁴

Faculty of Engineering and Technology, Future University in Egypt, New Cairo 11835, Egypt

⁵

Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt

⁶

Research Center for Physics (RCP), Department of Physics, Faculty of Science and Arts, Al-Mikhwah, Al-Baha University, Al-Baha 1988, Saudi Arabia

Author to whom correspondence should be addressed.

Symmetry 2023, 15(1), 57; https://doi.org/10.3390/sym15010057

Submission received: 30 November 2022 / Revised: 13 December 2022 / Accepted: 21 December 2022 / Published: 26 December 2022

(This article belongs to the Section Mathematics)

Download

Browse Figures

Versions Notes

Abstract

In this work, a damped modified Kawahara equation (mKE) with cubic nonlinearity and two dispersion terms including the third- and fifth-order derivatives is analyzed. We employ an effective semi-analytical method to achieve the goal set for this study. For this purpose, the ansatz method is implemented to find some approximate solutions to the damped mKE. Based on the proposed method, two different formulas for the analytical symmetric approximations are formally obtained. The derived formulas could be utilized for studying all traveling waves described by the damped mKE, such as symmetric solitary waves (SWs), shock waves, cnoidal waves, etc. Moreover, the energy of the damped dressed solitons is derived. Furthermore, the obtained approximations are used for studying the dynamics of the dissipative dressed (modified Kawahara (mK)) dust-ion acoustic (DIA) solitons in an unmagnetized collisional superthermal plasma consisting of inertia-less superthermal electrons and inertial cold ions as well as immobile negative dust grains. Numerically, the impact of the collisional parameter that arises as a result of taking the ion-neutral collisions into account and the electron spectral index on the profile of the dissipative structures are examined. Finally, the analytical and numerical approximations using the finite difference method (FDM) are compared in order to confirm the high accuracy of the obtained approximations. The achieved results contribute to explaining the mystery of several nonlinear phenomena that arise in different plasma physics, nonlinear optics, shallow water waves, oceans, and seas, and so on.

Keywords:

modified Kawahara equation; dissipative modified Kawahara equation solitons; collisional plasmas; superthermal electrons

1. Introduction

The study of differential equations (DEs) (e.g., ordinary, partial, linear, nonlinear, etc.) is interesting and inevitable. The significance of DEs is paramount for understanding how nature works, as it helps in modeling a vast number of problems in engineering, physics, chemistry, biology, medicine, etc. This enables us to understand how scientific phenomena behave and how to manifest the precise details of nature and to control its behavior [1,2]. There are many equations that have greatly succeeded in explaining the mechanism of solitary wave propagation in various physical models [3,4,5,6]. An example is the third- and/or fifth-order Korteweg–de Vries (KdV) equation

(\partial_{τ} ϕ + a ϕ \partial_{ζ} ϕ + b \partial_{ζ}^{3} ϕ = 0)

, which is considered one of the most famous equations of the twenty-first century due to its great role in modeling many nonlinear phenomena that arise in electrical circuits, seas, oceans, nonlinear optics, shallow water waves, in addition to plasma physics [7,8]. The KdV equation succeeded in describing the mechanism of propagation of many nonlinear structures (e.g., solitons) in various fields of science, especially in plasma, seas, and oceans. However, the KdV equation sometimes fails to model the nonlinear structures at some physical values of the parameters related to the model under consideration, which are called ‘critical parameters’. At these critical parameters, the coefficient of the nonlinear term “a” disappears, and the KdV equation becomes unsuitable for modeling nonlinear structures. Accordingly, many researchers have derived another high-order nonlinearity evolution equation to describe the nonlinear structures at the critical values of plasma compositions, which is called the ‘modified KdV (mKdV) equation’:

\partial_{τ} ϕ + c ϕ^{2} \partial_{ζ} ϕ + b \partial_{ζ}^{3} ϕ = 0

, where c represents the coefficient of the cubic nonlinear term. However, sometimes both KdV and mKdV equations fail to describe the nonlinear waves near/close to the critical values of some plasma compositions. Therefore, many authors have come up with another evolution equation, known as the ‘Extended/Gardner KdV (EKdV) equation’:

\partial_{τ} ϕ + (a ϕ + c ϕ^{2}) \partial_{ζ} ϕ + b \partial_{ζ}^{3} ϕ = 0

, which has succeeded in modeling many nonlinear waves near/close to the critical values of some plasma compositions.

It is noted that the aforementioned equations are called the ‘family of the third-order dispersion KdV equation’ or the ‘KdV-type equation’. On the contrary, there is another form of the KdV-type equation of a fifth-order dispersion term, which is called the ‘fifth-order KdV-type (fKdV-type) equation’ or sometimes the ‘Kawahara-type equation’, which includes the following family of the Kawahara equation (KE) [9,10,11,12,13]:

\begin{matrix} damped KE & : \partial_{τ} ϕ + a ϕ \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ = 0, \end{matrix}

(1)

\begin{matrix} damped KE & : \partial_{τ} ϕ + a ϕ \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ + R ϕ = 0, \end{matrix}

(2)

\begin{matrix} nonplanar KE & : \partial_{τ} ϕ + a ϕ \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ + \frac{m}{2 τ} ϕ = 0, \end{matrix}

(3)

\begin{matrix} nonplanar damped KE & : \partial_{τ} ϕ + a ϕ \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ + (R + \frac{m}{2 τ}) ϕ = 0, \end{matrix}

(4)

\begin{matrix} modified KE (mKE) & : \partial_{τ} ϕ + A ϕ^{2} \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ = 0 . \end{matrix}

(5)

Many other forms, such as forced KE, forced mKE, Gardner/Extended KE (EKE), and forced EKE, have been studied and solved using a variety of analytical and numerical approaches. Here, “a” and “B” are the coefficients of the quadratic nonlinear term and the third-order dispersion term, respectively. However, the parameter “d” is the coefficient of the fifth-order dispersion term, and

ϕ \equiv ϕ (ζ, τ) .

The family of the fifth-order dispersion term has many applications in several branches of science, such as capillary-gravity water waves, waves in plasma physics [14,15,16], shallow water waves (WWs), WWs with surface tension, and many other applications [17,18,19,20,21]. Due to the importance of this family in modeling many nonlinear physical problems, it has been analyzed and solved using several analytical and numerical techniques, e.g., the tanh method (TM) [22,23,24], tanh–coth method [22], a direct algebraic method [25], the homotopy analysis method (HAM) [26], sine–cosine [27], mesh-free method [28], septic B-spline collocation method [29], and so on [30]. For instance, Wazwaz [31] used many exact methods, including TM, extended TM, and sine–cosine method, to find a series of analytical solutions to the mKE using the sine–cosine method, the tanh method, and the ansatz method based on hyperbolic functions for analyzing and solving the mKE. Moreover, the homogeneous balance method (HBM) [32] has been employed for solving KE, and the authors [32] compared the obtained results with the (

G^{'} / G

)-expansion method solution and the F-expansion method solutions. Additionally, Gepreel et al. [33] used the improved general mapping deformation technique to find some exact analytic solutions to some nonlinear PDEs, such as the planar mKE. In addition, both the homotopy perturbation method (HPM) and the variational iteration method (VIM) [34] have been professionally applied to analyze the mKE and some related equations of the fifth-order dispersion term. The authors of [34] compared the numerical solutions of the VIM and HPM methods and the exact solutions to the mKE. The authors of [34] found that both VIM and HPM are characterized by their high accuracy and efficiency and ease of processing that many other numerical methods do not furnish. Biazar et al. [35] used both the VIM and Adomian decomposition method (ADM) to find some analytical approximations for both KE and mKE. Additionally, the KE has been solved numerically using the optimal homotopy asymptotic technique (OHAT) [36,37]. Further, the authors of [36] compared the solution obtained via the OHAT and the exact solution to this equation using VIM, HPM, and Variational HPM (VHPM) and found that the solution obtained using OHAT has higher accuracy than some of the indicated methods. Moreover, the OHAT [37] has been applied to determine an approximate solution to the mKE. Recently, El-Tantawy [13] used TM, the ansatz method based on the trigonometric, the Weierrtrass, and the Jacobi elliptic functions for deriving exact solutions to both mKE and EKE.

All the aforementioned studies focused on analyzing the family of differential equations that have a fifth-order dispersion term, including KE, mKE, EKE, etc., that support exact analytical solutions. However, there are many physical effects that if taken into account, such as the viscosity, collision, etc., eventually lead to some evolutionary equations with the fifth-order dispersion term with no exact analytical solution. For example, we report the case of the collision of charged particles with each other or neutral particles within some plasma systems. In this case, the fluid equations of the plasma model under study can be reduced to the family of damped Kawahara-type equations in order to study the properties of the nonlinear dissipative waves that propagate in a plasma model [38,39,40]. Further, if the geometrical effects are considered, then the set of fluid equations are reduced to the family of a completely non-integrable nonplanar (cylindrical and spherical) KE [11]. The families of the damped Kawahara-type equation, nonplanar Kawahara-type equation, and nonplanar damped Kawahara-type equation do not support exact solutions and must be solved using some numerical or semi-analytical methods. For instance, the non-integrable nonplanar damped KE (4) has been solved analytically using the ansatz method to derive semi-analytical solutions with high accuracy in order to study the properties of nonplanar dissipative Kawahara solitary and cnoidal waves in superthermal dusty plasma [12]. Further, the authors of [12] derived a general formula for the nonplanar dressed soliton energy. Sometimes, the collisional force between the charged particles cannot be neglected in some types of plasma compositions [41,42,43,44,45], especially cold plasma. Further, at certain values of plasma parameters, the nonlinear quadratic term disappears, so this forces us to search for higher-order nonlinearity orders to be able to describe the physical phenomenon at these values. Accordingly, the following new evolutionary equation, which is called ’damped mKE’ is obtained:

\partial_{τ} ψ + A ψ^{2} \partial_{ζ} ψ + B \partial_{ζ}^{3} ψ - d \partial_{ζ}^{5} ψ + R ψ = 0,

(6)

where

ψ \equiv ψ (ζ, τ)

, and R indicates the coefficient of the linear term that arises as a result of taking the collisional force between the charged and neutral plasma particles into account.

Equation (6) is used to describe large but finite amplitudes of dissipative nonlinear waves at some critical plasma compositions or any other fluid model described by this equation. To our knowledge and based on a comprehensive and thorough survey, this is the first analysis of the damped mKE using a semi-analytical approach. Thus, the main objective of the present investigation is to derive semi-analytical solutions/analytical approximations to Equation (6) using the ansatz method [10,11,12]. Moreover, the finite difference method (FDM) is applied to find numerical approximations for comparison with the semi-analytical solutions. Further, the global maximum residual error (GMRE) for all obtained approximations is estimated to confirm the effectiveness of all obtained approximations:

L_{R} = {FindMaximum|}_{℧} |\partial_{τ} ψ + A ψ^{2} \partial_{ζ} ψ + B \partial_{ζ}^{3} ψ - d \partial_{ζ}^{5} ψ + R ψ|,

(7)

where

℧ \in [ζ_{i}, ζ_{f}] \times [τ_{i}, τ_{f}]

represents the whole study domain.

The second objective of this investigation is to derive a formula for the energy of the damped modified Kawahara soliton and prove that it is not conserved due to the existence of the damping term. Further, we prove that the obtained formula of the damped modified Kawahara soliton energy can recover the energy of the conserved case. The third objective is to use all obtained semi-analytical solutions in order to study the distinctive characteristics of the damped (dissipative) modified Kawahara solitary waves (SWs) in a super-thermal dusty plasma composed of fluid inertial cold ions and superthermal electrons in addition to the immobile negative dust grains.

The paper is structured as follows: in Section 2, the anstaz method is carried out for analyzing the damped mKE in order to find semi-analytical solutions to the damped mKE (6). All derived semi-analytical solutions are implemented to find the dissipative modified Kawahara SWs, as illustrated in Section 3. For the physical application, the reductive perturbation method (RPM) is used for reducing the fluid equations of a superthermal dusty plasma to a damped mKE (6). After that, the impact of the plasma parameters on the profile (including the width and the amplitude) of the dissipative modified Kawahara solitons is discussed in Section 4. In Section 5, we briefly summarize the most important results we obtained.

2. The Anstaz Method for Analyzing the Damped mKE

Since the damped mKE (6) cannot be solved analytically via exact analytical methods due to the existence of the collisional term, in this section, we try to determine semi-analytical solutions to this equation using one of the semi-analytical methods that give good results with this type of non-integrable equation. Thus, the anstaz method is applied to find symmetric approximations to the damped mKE (6). This method is summarized as follows:

We assume the solution of the damped mKE (6) is in the form

$ψ (ζ, τ) = F_{1} ϕ (ζ F_{2}, F_{3}) \equiv F_{1} ϕ_{F} .$

(8)

Here, $F_{1} \equiv F_{1} (τ)$ , $F_{2} \equiv F_{2} (τ)$ , and $F_{3} \equiv F_{3} (τ)$ are, respectively, time-dependent functions related to the temporal amplitude, inverted width, and wave velocity of the damped structures. Remember that for $F_{1, 2} = 1$ and $F_{3} = t$ , the planar undamped case is recovered. Further, the function $ϕ (ζ, τ)$ indicates the exact solution to the mKE (5).
Inserting Equation (8) into the damped mKE (6), the following values of $\partial_{ζ}^{n} ψ$ and $\partial_{τ} ψ$ $(n = 1, 3, 5)$ are obtained:

$\{\begin{matrix} \partial_{ζ} ψ = F_{1} F_{2} \partial_{ζ} ϕ_{F}, \\ \partial_{ζ}^{3} ψ = F_{1} F_{2}^{3} \partial_{ζ}^{3} ϕ_{F}, \\ \partial_{ζ}^{5} ψ = F_{1} F_{2}^{5} \partial_{ζ}^{5} ϕ_{F}, \end{matrix}$

(9)

and

$\partial_{τ} ψ = {\dot{F}}_{1} ϕ_{F} + F_{1} [ζ {\dot{F}}_{2} \partial_{ζ} ϕ_{F} + {\dot{F}}_{3} (- A ϕ_{F}^{2} \partial_{ζ} ϕ_{F} - B \partial_{ζ}^{3} ϕ_{F} + d \partial_{ζ}^{5} ϕ_{F})] .$

(10)
Substituting the values of $\partial_{ζ}^{n} ψ$ and $\partial_{τ} ψ$ into the damped mKE (6) and reorganizing all terms, the following residual equation is obtained:

$Q_{1} \partial_{ζ}^{5} ϕ_{F} + Q_{2} \partial_{ζ}^{3} ϕ_{F} + Q_{3} ϕ_{F}^{2} \partial_{ζ} ϕ_{F} + Q_{4} \partial_{ζ} ϕ_{F} + Q_{5} ϕ_{F} = 0,$

(11)

with

$(\begin{matrix} Q_{1} \\ Q_{2} \\ Q_{3} \\ Q_{4} \\ Q_{5} \end{matrix}) = (\begin{matrix} (d) F_{1} (- F_{2}^{5} + {\dot{F}}_{3}), \\ B F_{1} (F_{2}^{3} - {\dot{F}}_{3}) \\ A F_{1} (F_{1}^{2} F_{2} - {\dot{F}}_{3}), \\ ζ F_{1} {\dot{F}}_{2}, \\ {\dot{F}}_{1} + F_{1} R, \end{matrix})$

where ${\dot{F}}_{i} \equiv \partial_{τ} F_{i}$ $(i = 1, 2, 3) .$
Now, for $Q_{i} = 0$ $(i = 1, 2, 3, 5)$ , we get the following system of differential equations:

$\begin{matrix} {\dot{F}}_{1} + F_{1} R & = 0, \end{matrix}$

(12)

$\begin{matrix} F_{1}^{2} F_{2} - {\dot{F}}_{3} & = 0, \end{matrix}$

(13)

$\begin{matrix} - F_{2}^{5} + {\dot{F}}_{3} & = 0, \end{matrix}$

(14)

$\begin{matrix} F_{2}^{3} - {\dot{F}}_{3} & = 0 . \end{matrix}$

(15)
Solving Equation (12) by using the initial condition (IC) $F_{1} (τ_{0}) = 1$ , the following value of $F_{1}$ is obtained:

$F_{1} = e^{R (τ_{0} - τ)} .$

(16)
It is clear from Equations (14) and (15) that there are two independent ODEs in $(F_{2}, F_{3})$ in addition to Equation (13) which means that there are two values for $(F_{2}, F_{3})$ that should be determined. By solving Equations (13) and (14) together, we get the first values for $(F_{2}, F_{3}) .$ Further, by solving Equations (13) and (15) together, we get the second values for $(F_{2}, F_{3})$ . Inserting the value of $F_{1},$ given in Equation (16), into Equation (13), yields

$(e^{2 R (τ_{0} - τ)} F_{2} - {\dot{F}}_{3}) = 0 .$

(17)
For the first values of $(F_{2}, F_{3})$ : inserting Equation (14) $({\dot{F}}_{3} = F_{2}^{5})$ into (17) yields

$(e^{2 R (τ_{0} - τ)} F_{2} - F_{2}^{5}) = 0$

(18)

and by solving Equation (18), we get

$F_{2} = e^{\frac{R}{2} (τ_{0} - τ)} .$

(19)

To find the value of $F_{3}$ corresponding to the value of $F_{2}$ given in Equation (19), we must solve Equation (14) using both IC $F_{3} |_{τ = τ_{0}} = 0$ and the value of $F_{2}$ given in Equation (19), which leads to

$F_{3} = \frac{2}{5 R} (1 - e^{\frac{5 R}{2} (τ_{0} - τ)}),$

(20)

It is obvious that the values of $(F_{2}, F_{3})$ given in Equations (19) and (20) can be recovered in many cases, which is discussed later.
For the second values of $(F_{2}, F_{3})$ : Equation (15) $({\dot{F}}_{3} = F_{2}^{3})$ is inserted into (17) to get

$(e^{2 R (τ_{0} - τ)} F_{2} - F_{2}^{3}) = 0$

(21)

and by solving Equation (21), the following new value of $F_{2}$ is obtained

$F_{2} = e^{R (τ_{0} - τ)} .$

(22)

Now, by replacing the new value of $F_{2}$ , given in Equation (22), into Equation (15), and using the IC $F_{3} |_{τ = τ_{0}} = 0$ , we finally get

$F_{3} = \frac{1}{3 R} (1 - e^{3 R (τ_{0} - τ)}),$

(23)
By inserting both the first and second values of $(F_{2}, F_{3})$ as well as the value of $F_{1}$ given in Equation (16) into the ansatz solution (8), the following two general approximate solutions of the damped mKE (6) are obtained:

$ψ_{1} (ζ, τ) = e^{R (τ_{0} - τ)} ϕ [ζ e^{\frac{R}{2} (τ_{0} - τ)}, \frac{2}{5 R} (1 - e^{\frac{5 R}{2} (τ_{0} - τ)})],$

(24)

and

$ψ_{2} (ζ, τ) = e^{R (τ_{0} - τ)} ϕ [ζ e^{R (τ_{0} - τ)}, \frac{1}{3 R} (1 - e^{3 R (τ_{0} - τ)})] .$

(25)

3. Dissipative Modified Kawahara Solitons and Their Energy

To find the dissipative SWs of the damped mKE (6), firstly, we need the exact soliton solution to the mKE (5). As we demonstrated in the introductory section, the mKE (5) has been solved using different analytical and numerical techniques. One of the symmetric exact dressed soliton solutions of the mKE (5) reads [13]

ϕ (ζ, τ) = ϕ_{m a x} {sech}^{2} [W (ζ - \frac{4 B^{2}}{25 d} τ)] .

(26)

Here,

ϕ_{m a x} = \frac{3 B}{\sqrt{10 A d}}

is the maximum amplitude of the SWs, while

W = \frac{1}{2} \sqrt{\frac{B}{5 d}}

gives the inverse width of the SWs. Now, to obtain the dissipative dressed soliton solution of the damped mKE (6), we should insert Equation (26) into the general approximations (24) and (25):

ψ_{1} (ζ, τ) = ϕ_{m a x} e^{R (τ_{0} - τ)} {sech}^{2} [W (ζ e^{\frac{R}{2} (τ_{0} - τ)} - \frac{4 B^{2}}{25 d} \frac{2}{5 R} (1 - e^{\frac{5 R}{2} (τ_{0} - τ)}))],

(27)

and

ψ_{2} (ζ, τ) = ϕ_{m a x} e^{R (τ_{0} - τ)} {sech}^{2} [W (ζ e^{R (τ_{0} - τ)} - \frac{4 B^{2}}{25 d} \frac{1}{3 R} (1 - e^{3 R (τ_{0} - τ)}))] .

(28)

In the following section, the semi-analytical solutions (27) and (28) are applied to study the properties of the dissipative modified Kawahara solitons that can be generated in an electronegative dusty plasmas.

To find the energy of the dissipative modified Kawahara soliton, we first determine the planar modified Kawahara soliton energy and prove that it is conserved (does not change with time). The energy of the modified Kawahara soliton can be obtained using the following conservation law:

\begin{matrix} I & = \int_{- \infty}^{\infty} ϕ \partial_{τ} ϕ \partial_{ζ} = \int_{- \infty}^{\infty} \frac{1}{2} \partial_{τ} ϕ^{2} \partial_{ζ} = \partial_{τ} (\int_{- \infty}^{\infty} \frac{1}{2} ϕ^{2} \partial_{ζ}) = \partial_{τ} E_{m K E} \\ = - \int_{- \infty}^{\infty} ϕ (A ϕ^{2} \partial_{ζ} ϕ + B \partial_{ζ}^{3} ϕ - d \partial_{ζ}^{5} ϕ) \partial_{ζ} \\ = - \int_{- \infty}^{\infty} (A ϕ^{3} \partial_{ζ} ϕ + B ϕ \partial_{ζ}^{3} ϕ - d ϕ \partial_{ζ}^{5} ϕ) \partial_{ζ} \\ = - A \int_{- \infty}^{\infty} (ϕ^{3} \partial_{ζ} ϕ) \partial_{ζ} - B \int_{- \infty}^{\infty} (ϕ \partial_{ζ}^{3} ϕ) \partial_{ζ} + d \int_{- \infty}^{\infty} (ϕ \partial_{ζ}^{5} ϕ) \partial_{ζ} \\ = 0 . \end{matrix}

(29)

Thus, the planar modified Kawahara SWs energy is conserved and is given by

\begin{matrix} E_{m K E} & = \frac{1}{2} \int_{- \infty}^{\infty} ϕ^{2} \partial_{ζ} = \frac{1}{2} ϕ_{m a x}^{2} \int_{- \infty}^{\infty} {sech}^{4} [W (ζ - \frac{4 B^{2}}{25 d} τ)] \partial_{ζ} \\ = \frac{2}{3} \frac{ϕ_{m a x}^{2}}{W} . \end{matrix}

(30)

We shall now proceed to derive the energy of the dissipative modified Kawahara soliton. As we proved above, the modified Kawahara soliton energy is conserved, i.e.,

\partial_{τ} E_{0} = 0

. On the other side, since all physical quantities (amplitude, width, and velocity) associated with these waves change with time, their energy becomes unconserved, i.e.,

\partial_{τ} E_{R} (τ) \neq 0

. The value of

E_{R} (τ)

can be obtained as follows:

\begin{matrix} \partial_{τ} E_{R} (τ) & = \partial_{τ} (\int_{- \infty}^{\infty} \frac{1}{2} ψ^{2} \partial_{ζ}) = \int_{- \infty}^{\infty} ψ \partial_{τ} ψ \partial_{ζ} \\ = - \int_{- \infty}^{\infty} ψ (A ψ^{2} \partial_{ζ} ψ + B \partial_{ζ}^{3} ψ - d \partial_{ζ}^{5} ψ + R ψ) \partial_{ζ} \\ = - \int_{- \infty}^{\infty} (A ψ^{3} \partial_{ζ} ψ + B ψ \partial_{ζ}^{3} ψ - d ψ \partial_{ζ}^{5} ψ + R ψ^{2}) \partial_{ζ} \\ = - 2 R \int_{- \infty}^{\infty} \frac{1}{2} ψ^{2} \partial_{ζ} = - 2 R E_{R} (τ) . \end{matrix}

(31)

By integrating Equation (31) once over

τ

, we obtain

\begin{matrix} \int \frac{\partial_{τ} E_{R} (τ)}{E_{R} (τ)} & = - 2 \int R \partial_{τ} + c c \Rightarrow \\ ln (E_{R} (τ)) & = - 2 R τ + c c . \end{matrix}

(32)

To determine the value of the integration constant

c c

, we use the IC

E_{R} (τ_{0}) = E_{0}

, which leads to

c c = ln (E_{0}) + 2 R τ_{0} .

(33)

Inserting Equation (33) into Equation (32), we get

\begin{matrix} ln (E_{R} (τ)) & = - 2 R τ + ln (E_{0}) + 2 R τ_{0} \\ \Rightarrow ln (\frac{E_{R} (τ)}{E_{0}}) = 2 R (τ_{0} - τ) \\ \Rightarrow E_{R} (τ) = E_{0} e^{2 R (τ_{0} - τ)} . \end{matrix}

(34)

It is noticed that the Kawahara soliton energy is a function of time and decays with the increase of time, which confirms that the energy of the soliton is not conserved. It is noteworthy that for

R = 0

, the conserved soliton energy of the mKE is recovered as

E_{R} (τ) = E_{0} .

(35)

4. Dissipative Modified Kawahara Solitons in a Dusty Plasma

Now, we can apply the analytical approximations (24) and (25) for studying the dissipative dust-ion acoustic modified Kawahara solitons in a collisional super-thermal complex plasma consisting of inertialess superthermal electrons and inertial cold ions in addition to stationary negative dust grains to preserve the neutrality condition. In this model, the ion-neutral collisions are considered. Thus, the set of normalized fluid equations that govern the dynamics of dust-ion acoustic waves (DIAWs) are given by: [46]

The continuity equation of the positive ions:

\partial_{t} n_{i} + \partial_{x} (n_{i} u_{i}) = 0,

(36)

the momentum equation of the positive ions:

\partial_{t} u_{i} + u_{i} \partial_{x} u_{i} + υ u_{i} - \partial_{x} φ = 0,

(37)

and Poission’s equation:

\partial_{x}^{2} φ + n_{i} - n_{e} - α_{d} = 0,

(38)

while the superthermal electrons’ normalized density

n_{e}

is given by:

\begin{matrix} n_{e} & = μ {[1 - \frac{φ}{κ - \frac{3}{2}}]}^{- κ + \frac{1}{2}} \\ = μ (1 + α_{1} φ + α_{2} φ^{2} + α_{3} φ^{3} + \dots), \end{matrix}

(39)

with

\begin{matrix} α_{1} & = \frac{(2 κ - 1)}{(2 κ - 3)}, \\ α_{2} & = \frac{(2 κ - 1) (2 κ + 1)}{2 {(2 κ - 3)}^{2}} \\ α_{3} & = \frac{(2 κ - 1) (2 κ + 1) (2 κ + 3)}{6 {(2 κ - 3)}^{3}} . \end{matrix}

Equations (36)–(39) are written in the dimensionless form, i.e., the dependent quantities

(n_{i}, u_{i}, φ, n_{e})

and independent variables

(υ, x, t)

are written in normalized form. Here,

u_{i}

and

n_{i}

indicate the normalized velocity and density of the positive ions, respectively,

ϕ

expresses the normalized wave potential,

υ

is the normalized ion-neutral collisional frequency, while

(x, t)

represent the scaled space-time independent variables [46]. Here,

κ

indicates the spectral index of the superthermal electron distribution

(1.5 < κ)

μ

and

α_{d}

are, respectively, the electron and negative dust concentrations that satisfy the following neutrality condition:

μ = 1 - α_{d}

The RPM is implemented to reduce the fluid Equations (36)–(39) to an evolution equation governing the propagation of nonlinear structures that can be created in the present plasma model. Accordingly, the following, stretched to the independent variables, are considered:

(ζ, τ, υ) = (ε (x - V_{p h} t), ε^{3} t, ε^{3} υ_{0}),

(40)

while the following expansion of the dependent quantities are as follows:

Y (x, t) = Y^{(0)} + \sum_{j = 1}^{\infty} ε^{j} Y^{(j)} (ζ, τ)

(41)

with

\begin{matrix} Y (x, t) & \equiv [n_{i}, u_{i}, φ], \\ Y^{(j)} (ζ, τ) & \equiv [n_{i}^{(j)}, u_{i}^{(j)}, φ^{(j)}], \\ Y^{(0)} & \equiv [1, 0, 0], \end{matrix}

where

V_{p h}

expresses the normalized phase velocity of DIAWs, and

ε

is a small real parameter (

ε < < 1

Inserting Equations (40) and (41) into Equations (36)–(39) and rearranging all terms, we get reduced equations with different orders of

ε

. For the lowest-orders of

ε

, the first-order quantities

(n_{i}^{(1)}, u_{i}^{(1)})

and the normalized phase velocity

V_{p h}

are obtained:

n_{i}^{(1)} = \frac{1}{V_{p h}} u_{i}^{(1)} = \frac{1}{V_{p h}^{2}} φ^{(1)} and V_{p h} = \frac{1}{\sqrt{α_{1}}} .

(42)

For the next-order to

ε

, we get the second-order quantities

(n_{i}^{(2)}, u_{i}^{(2)})

as follows:

\{\begin{matrix} n_{i}^{(2)} = \frac{1}{V_{p h}^{4}} (\frac{3}{2} φ^{(1) 2} + V_{p h}^{2} φ^{(2)}), \\ u_{i}^{(2)} = \frac{1}{V_{p h}^{3}} (\frac{1}{2} φ^{(1) 2} + V_{p h}^{2} φ^{(2)}) . \end{matrix}

(43)

In addition, Poisson’s equation gives us

\frac{1}{V_{p h}} (2 α_{2} V_{p h}^{4} - 3) φ^{(1) 2} + 2 V_{p h} (α_{1} V_{p h}^{2} - 1) φ^{(2)} = 0 .

(44)

In Equation (44), the coefficient of the second term represents the dispersion relation, i.e.,

(α_{1} V_{p h}^{2} - 1) = 0

, while the coefficient of the first term represents the coefficient of the quadratic nonlinear term in the damped KdV equation, which must vanish (i.e.,

2 α_{2} V_{p h}^{4} - 3 = 0

) to get the damped mKdV equation. Now, by solving

(2 α_{2} V_{p h}^{4} - 3) = 0

, the critical value of dust concentration

α_{d c}

is obtained:

α_{d c} = \frac{4}{3} \frac{(κ - 1)}{(2 κ - 1)} .

(45)

Note here that the value of

α_{d c}

must satisfy the neutrality condition:

μ = 1 - α_{d c}

For the next-higher order to

ε

, we get

\begin{matrix} \partial_{τ} n_{i}^{(1)} + \partial_{ζ} (n_{i}^{(1)} u_{i}^{(2)}) + \partial_{ζ} (n_{i}^{(2)} u_{i}^{(1)}) - V_{p h} \partial_{ζ} n_{i}^{(3)} + \partial_{ζ} u_{i}^{(3)} = 0 \end{matrix}

(46)

\begin{matrix} υ_{0} u_{i}^{(1)} + \partial_{τ} u_{i}^{(1)} + \partial_{ζ} (u_{i}^{(1)} u_{i}^{(2)}) - V_{p h} \partial_{ζ} u_{i}^{(3)} + \partial_{ζ} φ^{(3)} = 0, \end{matrix}

(47)

while Poission’s equation gives us:

n_{i}^{(3)} - α_{3} φ^{(1) 3} - 2 α_{2} φ^{(1)} φ^{(2)} - α_{1} φ^{(3)} + \partial_{ζ}^{2} φ^{(1)} = 0 .

(48)

Solving Equations (46)–(48) and taking the fifth-order derivative/perturbation into consideration, the following damped mKE is obtained:

\partial_{τ} ψ + A ψ^{2} \partial_{ζ} ψ + B \partial_{ζ}^{3} ψ - d \partial_{ζ}^{5} ψ + R ψ = 0,

(49)

with

A = \frac{1}{4 λ^{3}} (15 - 6 α_{3} λ^{6}) and B = \frac{λ^{3}}{2},

(50)

where

R = υ_{0} / 2

denotes the coefficient of the damping term and

ψ \equiv φ^{(1)}

Before proceeding to analyze the effect of plasma parameters on the dissipative modified Kawahara solitons, we should study the impact of the electron spectral index

κ

and the collisional frequency R on the dressed soliton energy

E_{R} (τ)

. Figure 1 demonstrates the relation between the parameters

(κ, R)

and the energy

E_{R} (τ) .

One can see that the energy

E_{R} (τ)

increases (decreases) with the increase in

κ

(R). We can deduce an important result from this analysis: as long as the energy

E_{R} (τ)

of the dressed dissipative soliton decreases with a specific parameter, this, in turn, leads to a decrease in the speed and nonlinearity of the dressed soliton, and therefore, it works to reduce the amplitude of the dressed soliton and vice versa.

Now, we can study the impact of the physical parameters

(κ, R)

on the properties of the symmetric dissipative dressed solitons using the two approximations (27) and (28). Figure 2 exhibits the effect of the ion-neutral collisional frequency R on the profiles of the dissipative dressed solitons using the two approximations (27) and (28). We can observe that by increasing R, the dissipative dressed soliton amplitude decreases. A similar trend is observed as time goes on, as illustrated in Figure 3, i.e., the dissipative dressed soliton amplitude decays over time. It is important to note here that for ultrafine values of

υ_{0} / R

, the nonlinearity of the DIA dissipative dressed solitons increases, which leads to the enhancement of the soliton amplitude, and the profile of the dissipative dressed solitons becomes identical to the planar dressed solitons. Figure 4 is considered for investigating the impact of the electron spectral index

κ

on the profile of the DIA dissipative dressed solitons. It is found that the amplitude of the dressed solitons increases with

κ

, which means that the dressed soliton nonlinearity and energy increase with increasing

κ

(this result is confirmed above during analyzing the soliton energy). Moreover, the GMRE

L_{R}

is estimated for

(κ, R, d) = (2, 0.1, 0.1)

, as shown in Table 1.

It is observed from Table 1 and Figure 5 that there is almost complete harmony between all semi-analytical solutions and numerical solutions using FDM, but there are some minor differences in the accuracy of each them. As shown in Table 1, the accuracy of the first formula for the semi-analytical solution (24) or (27) is better than both the second formula for the semi-analytical solution (25) or (28) and the FDM numerical approximation. Furthermore, the obtained results are compared with numerical solutions using Wolfram Mathematica. It is noticed that the semi-analytical solutions (27) and (28) as well as the FDM numerical approximations are better than the numerical solutions using Wolfram Mathematica. In general, the obtained semi-analytical solutions are characterized by high-accuracy and are more stable over a long period of time. Moreover, we analyzed the damped mKE (49) numerically using HPM, but unfortunately, this method gave a bad result over a long period of time; in addition, the accuracy of its approximations was not good due to the existence of the damping term. Thus, we did not add these results here.

5. Conclusions

The damped modified Kawahara equation (mKE) was analyzed via an ansatz method to find semi-analytical symmetric solutions. After applying the proposed method, two general semi-analytical solutions were formally derived. All derived symmetric approximations can recover all traveling wave solutions to the mKE. As an example to the derived approximations, the dissipative modified Kawahara solitary wave solution was obtained. Moreover, a generic energy formula of the dissipative dressed solitons was derived. Further, the dust-ion acoustic (DIA) dissipative solitons in a superthermal plasma were studied using the derived approximations. The effects of the collisional frequency R and the electron spectral index

κ

on the profile of the dissipative dressed solitons were studied and discussed. It was found that the two parameters have an opposite effect on the soliton profile, i.e., the dissipative dressed soliton amplitude increases (decreases) with increasing

κ

(R). Further, the accuracy of both analytical and numerical approximations was investigated by estimating the global maximum residual error along the whole study domain. We observed that numerical approximation using the finite difference method (FDM) is better than approximation (25) or (28). The current study can contribute to studying the mechanisms of various nonlinear phenomena that arise in several plasma models, nonlinear optics, fluid mechanics, etc.

Future work: Due to the importance of fractional differential equations [47] in explaining many mysterious phenomena, we will employ some effective methods for analyzing the damped fractional Kawahara-type equation and some linked fractional equations with the fifth-order dispersion term.

Author Contributions

Conceptualization, S.M.E.I. and S.A.E.-T.; methodology, A.-M.W. and S.A.E.-T.; software, S.A.E.-T.; validation, A.-M.W., E.T.-E. and S.A.E.-T.; formal analysis, A.-M.W., E.T.-E. and S.A.E.-T. investigation, S.M.E.I. and S.A.E.-T.; resources, S.M.E.I., E.T.-E. and S.A.E.-T.; data curation, A.-M.W. and S.A.E.-T.; writing—original draft preparation, S.M.E.I., E.T.-E. and S.A.E.-T.; writing—review and editing, S.M.E.I., A.-M.W. and S.A.E.-T.; visualization, E.T.-E.; supervision, A.-M.W. and S.A.E.-T.; project administration, S.M.E.I. and S.A.E.-T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IF2/PSAU/2022/01/21730).

Conflicts of Interest

The authors declare no conflict of interest.

References

Ray, S.S. Nonlinear Differential Equations in Physics; Springer: Singapore, 2020. [Google Scholar] [CrossRef]
Muravnik, A.B. Nonlocal problems and functional-differential equations: Theoretical aspects and applications to mathematical modelling. Math. Model. Nat. Phenom. 2019, 14, 601. [Google Scholar] [CrossRef]
Islam, S.M.R.; Kumar, D.; Fendzi-Donfack, E.; Inc, M. Impacts of nonlinearity and wave dispersion parameters on the soliton pulses of the (2+1)-dimensional Kundu–Mukherjee–Naskar equation. Rev. Mex. Fís. 2022, 68, 061301. [Google Scholar] [CrossRef]
Islam, S.M.R.; Khan, S.; Arafat, M.S.Y.; Akbar, M.A. Diverse analytical wave solutions of plasma physics and water wave equations. Results Phys. 2022, 40, 105834. [Google Scholar] [CrossRef]
Arafat, S.M.Y.; Khan, K.; Islam, S.M.R.; Rahman, M.M. Parametric effects on paraxial nonlinear Schrödinger equation in Kerr media. Chin. J. Phys. 2022. [Google Scholar] [CrossRef]
Fatema, K.; Islam, M.E.; Arafat, S.M.Y.; Akbar, M.A. Solitons’ behavior of waves by the effect of linearity and velocity of the results of a model in magnetized plasma field. J. Ocean. Eng. Sci. 2022. [Google Scholar] [CrossRef]
Iniewski, K. Circuits at the Nanoscale: Communications, Imaging, and Sensing, 1st ed.; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
Crighton, D.G. Applications of KdV. Acta Appl. Math. 1995, 39, 39. [Google Scholar] [CrossRef]
Wazwaz, A.-M. A Fifth-Order Korteweg–de Vries Equation for Shallow Water with Surface Tension: Multiple Soliton Solutions. Acta Phys. Pol. A 2016, 130, 679. [Google Scholar] [CrossRef]
Aljahdaly, N.H.; El-Tantawy, S.A. Novel anlytical solution to the damped Kawahara equation and its application for modeling the dissipative nonlinear structures in a fluid medium. J. Ocean Eng. Sci. 2022, 7, 492. [Google Scholar] [CrossRef]
Alharthi, M.R.; Alharbey, R.A.; El-Tantawy, S.A. Novel analytical approximations to the nonplanar Kawahara equation and its plasma applications. Eur. Phys. J. Plus. 2022, 137, 1172. [Google Scholar] [CrossRef]
El-Tantawy, S.A.; El-Sherif, L.S.; Bakry, A.M.; Alhejaili, W.; Wazwaz, A.-M. On the analytical approximations to the nonplanar damped Kawahara equation: Cnoidal and solitary waves and their energy. Phys. Fluids 2022, 34, 113103. [Google Scholar] [CrossRef]
El-Tantawy, S.A.; Salas, A.H.; Alharthi, M.R. Novel analytical cnoidal and solitary wave solutions of the Extended Kawahara equation. Chaos Solitons Fractals 2021, 147, 110965. [Google Scholar] [CrossRef]
El-Shewy, E.K.; Zahran, M.A.; Schoepf, K.; Elwakil, S.A. Contribution of higher order dispersion to nonlinear dust-acoustic solitary waves in dusty plasma with different sized dust grains and nonthermal ions. Phys. Scr. 2008, 78, 025501. [Google Scholar] [CrossRef]
Elwakil, S.A.; El-Shewy, E.K.; Zahran, M.A. Higher-order solution of an electron acoustic solitary waves via vortex electron distribution. Chaos Solitons Fractals 2004, 22, 13. [Google Scholar] [CrossRef]
Alyousef, H.A.; Salas, A.H.; Alharthi, M.R.; El-tantawy, S.A. New Periodic and Localized Traveling Wave Solutions to a Kawahara-Type Equation: Applications to Plasma Physics. Complexity 2022, 2022, 9942267. [Google Scholar] [CrossRef]
Berloff, N.G.; Howard, L.N. Solitary and periodic solutions to nonlinear nonintegrable equations. Stud. Appl. Math. 1997, 99, 1–24. [Google Scholar] [CrossRef]
Bridges, T.; Derks, G. Linear instability of solitary wave solutions of the Kawahara equation and its generalizations. SIAM. J. Math. Anal. 2002, 33, 1356–1378. [Google Scholar] [CrossRef] [Green Version]
Hunter, J.K.; Scheurle, J. Existence of perturbed solitary wave solutions to a model equation for water waves. Phys. D 1998, 32, 253–268. [Google Scholar] [CrossRef]
Kakutani, T.; Ono, H. Weak non-linear hydromagnetic waves in a cold collision-free plasma. J. Phys. Soc. Jpn. 1969, 26, 1305. [Google Scholar] [CrossRef]
Kawahara, T. Oscillatory solitary waves in dispersive media. J. Phys. Soc. Japan 1972, 33, 260–264. [Google Scholar] [CrossRef]
Wazwaz, A.-M. Partial Differential Equations and Solitary Waves Theory; Higher Education Press: Beijing, China, 2009. [Google Scholar]
Wazwaz, A.-M. New solitary wave solutions to the Kuramoto-Sivashinsky and the Kawahara equations. Appl. Math. Comput. 2006, 182, 1642. [Google Scholar] [CrossRef]
Mohamad-Jawad, A.; Slibi, T. Traveling Wave Solutions using Tanh Method for solving Kawahara and Modified Kawahara Equations. J. -Rafidain Univ. Coll. Sci. 2012, 2012, 21–34. [Google Scholar] [CrossRef]
Sirendaoreji. New exact travelling wave solutions for the Kawahara and modified Kawahara equations. Chaos Solitons Fractals 2004, 19, 147. [Google Scholar] [CrossRef]
Kurulay, M. Approximate analytic solutions of the modified Kawahara equation with homotopy analysis method. Adv. Differ. Equ. 2012, 2012, 178. [Google Scholar] [CrossRef] [Green Version]
Yusufoğlu, E.; Bekir, A.; Alp, M. Periodic and solitary wave solutions of Kawahara and modified Kawahara equations by using Sine–Cosine method. Chaos Solitons Fractals 2008, 37, 1193. [Google Scholar] [CrossRef]
Zarebnia, M.; Aghili, M. A new approach for numerical solution of the modified Kawahara equation. J. Nonlinear Anal. Appl. 2016, 2016, 48. [Google Scholar] [CrossRef] [Green Version]
Ak, T.; Karakoc, S.B.G. A numerical technique based on collocation method for solving modified Kawahara equation. J. Ocean Eng. Sci. 2018, 3, 67. [Google Scholar] [CrossRef]
Mancas, C. Traveling wave solutions to Kawahara and related equations. Differ. Equations Dyn. Syst. 2019, 27, 19–37. [Google Scholar] [CrossRef] [Green Version]
Wazwaz, A.-M. New solitary wave solutions to the modified Kawahara equation. Phys. Lett. A 2007, 360, 588. [Google Scholar] [CrossRef]
Taha, W.M.; Ibrahim, I.A.; Noorani, M.S.M. New Application for Generalized Regularized Long Wave (GRLW) Equation, Modified Dispersive Water Wave (MDWW) Equation and Kawahara Equation by Homogeneous Balance Method. J. Eng. Appl. Sci. 2019, 14, 4010. [Google Scholar] [CrossRef]
Gepreel, K.A.; Nofal, T.A. Exact Solutions for Nonlinear Partial Differential Equations and Its Applications. Int. J. Pure Appl. Math. 2013, 86, 241. [Google Scholar] [CrossRef]
Jin, L. Application of variational iteration method and homotopy perturbation method to the modified Kawahara equation. Math. Comput. Model. 2009, 49, 573. [Google Scholar] [CrossRef]
Biazar, J.; Gholamin, P.; Hosseini, K. Variational Iteration and Adomian Decomposition Methods for Solving Kawahara and Modified Kawahara Equations. Appl. Math. Sci. 2008, 2, 2705. [Google Scholar]
Kashkari, B.S. Application of Optimal Homotopy Asymptotic Method for the Approximate Solution of Kawahara Equation. Appl. Math. Sci. 2014, 8, 875. [Google Scholar] [CrossRef] [Green Version]
Ullah, H.; Nawaz, R.; Islam, S.; Idrees, M.; Fiza, M. The optimal homotopy asymptotic method with application to modified Kawahara equation. J. Assoc. Arab. Univ. Basic Appl. Sci. 2015, 18, 82. [Google Scholar] [CrossRef]
Alyousef, H.A.; Salas, A.H.; Matoog, R.T.; El-Tantawy, S.A. On the analytical and numerical approximations to the forced damped Gardner Kawahara equation and modeling the nonlinear structures in a collisional plasma. Phys. Fluids 2022, 34, 103105. [Google Scholar] [CrossRef]
El-Tantawy, S.A.; Salas, A.H.; Alyousef, H.A.; Alharthi, M.R. Novel exact and approximate solutions to the family of the forced damped Kawahara equation and modeling strong nonlinear waves in a plasma. Chin. J. Phys. 2022, 77, 2454. [Google Scholar] [CrossRef]
Aljahdaly, N.H.; Ashi, H.A.; Wazwaz, A.-M.; El-Tantawy, S.A. Exponential time differencing scheme for Modeling the dissipative Kawahara solitons in a two-electrons collisional plasma. Rom. Rep. Phys. 2022, 74, 109. [Google Scholar]
Kashkari, B.S.; El-Tantawy, S.A. Homotopy perturbation method for modeling electrostatic structures in collisional plasmas. Eur. Phys. J. Plus. 2021, 136, 121. [Google Scholar] [CrossRef]
Masood, W.; Rizvi, H.; Hasnain, H.; Batool, N. Dust drift shock waves with non-Maxwellian ion population in nonuniform collisional dusty plasmas in planetary environments. Astrophys Space Sci. 2013, 345, 49. [Google Scholar] [CrossRef]
Shan, S.A. Dissipative electron-acoustic solitons in a cold electron beam plasma with superthermal trapped electrons. Astrophys Space Sci. 2019, 364, 36. [Google Scholar] [CrossRef]
Khan, S.; Rahman, A.U.; Hadi, F.; Zeb, A.; Khan, M.Z. Weakly dissipative dust acoustic solitons in the presence of superthermal particles. Contrib. Plasma Phys. 2017, 57, 223. [Google Scholar] [CrossRef]
Masood, W.; Karim, S.; Shah, H.A.; Siddiq, M. Coupled nonlinear drift and ion acoustic waves in dense dissipative electron-positron-ion magnetoplasmas. Phys. Plasmas 2009, 16, 112302. [Google Scholar] [CrossRef]
Pakzad, H.R.; Nobahar, D. Dust-ion acoustic solitons in superthermal dusty plasmas. New Astron. 2022, 93, 101752. [Google Scholar] [CrossRef]
Ünal, S.Ç.; Daşcıoğlu, A.; Bayram, D.V. New exact solutions of space and time fractional modified kawahara equation. Phys. A 2020, 551, 124550. [Google Scholar]

Figure 1. The damped modified Kawahara soliton energy

E_{R} (τ)

is plotted against the physical parameters

(κ, R)

Figure 1. The damped modified Kawahara soliton energy

E_{R} (τ)

is plotted against the physical parameters

(κ, R)

Figure 2. The profile of the dissipative dressed solitons according to (a) the approximation (27) and (b) the approximation (28) is plotted against the collisional frequency parameter R.

Figure 3. The profile of the dissipative dressed solitons according to (a) the approximation (27) and (b) the approximation (28) is plotted against the time of propagation

τ

Figure 3. The profile of the dissipative dressed solitons according to (a) the approximation (27) and (b) the approximation (28) is plotted against the time of propagation

τ

Figure 4. The profile of the dissipative dressed solitons according to (a) the approximation (27) and (b) the approximation (28) is plotted against the electron spectral index

κ

Figure 4. The profile of the dissipative dressed solitons according to (a) the approximation (27) and (b) the approximation (28) is plotted against the electron spectral index

κ

Figure 5. A comparison between (a) the approximations (27) and (28) as well as FDM approximation and (b) the approximations (27) and (28) as well as the numerical solution using Wolfram Mathematica.

Table 1. In this table, the residual error

L_{R}

is estimated for all obtained approximations.

Table 1. In this table, the residual error

L_{R}

is estimated for all obtained approximations.

Approximation	$L_{R}$
Approx. (27)	0.0205979
Approx. (28)	0.0399192
FDM Approx.	0.0206002
Mathematica	0.0607945

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Ismaeel, S.M.E.; Wazwaz, A.-M.; Tag-Eldin, E.; El-Tantawy, S.A. Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma. Symmetry 2023, 15, 57. https://doi.org/10.3390/sym15010057

AMA Style

Ismaeel SME, Wazwaz A-M, Tag-Eldin E, El-Tantawy SA. Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma. Symmetry. 2023; 15(1):57. https://doi.org/10.3390/sym15010057

Chicago/Turabian Style

Ismaeel, Sherif M. E., Abdul-Majid Wazwaz, Elsayed Tag-Eldin, and Samir A. El-Tantawy. 2023. "Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma" Symmetry 15, no. 1: 57. https://doi.org/10.3390/sym15010057

APA Style

Ismaeel, S. M. E., Wazwaz, A. -M., Tag-Eldin, E., & El-Tantawy, S. A. (2023). Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma. Symmetry, 15(1), 57. https://doi.org/10.3390/sym15010057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Simulation Studies on the Dissipative Modified Kawahara Solitons in a Complex Plasma

Abstract

1. Introduction

2. The Anstaz Method for Analyzing the Damped mKE

3. Dissipative Modified Kawahara Solitons and Their Energy

4. Dissipative Modified Kawahara Solitons in a Dusty Plasma

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI