Approximate Analytical and Numeric Solutions to a Forced Damped Gardner Equation

In this paper, some exact traveling wave solutions to the integrable Gardner equation are reported. The ansatz method is devoted for deriving some exact solutions in terms of Jacobi and Weierstrass elliptic functions. The obtained analytic solutions recover the solitary waves, shock waves, and cnoidal waves. Also, the relation between the Jacobi and Weierstrass elliptic functions is obtained. In the second part of this work, we derive some approximate analytic and numeric solutions to the nonintegrable forced damped Gardner equation. For the approximate analytic solutions, the ansatz method is considered. With respect to the numerical solutions, the evolution equation is solved using both the finite different method (FDM) and cubic B-splines method. A comparison between different approximations is reported.

where u ≡ u(x, t) and (a, b, c) represent the coefficients of the nonlinear and dispersion terms which are function of physical parameters related to the system under study. e Gardner equation (GE) has two nonlinear terms in the quadratic and cubic forms and the dissipative term is of third-order derivative. e GE is an integrable system and Miura transformation connects it to the KdV equation. e GE is a useful model to understand the propagation of acoustic waves in different plasma models [17][18][19]. is equation can be used for describing the small but finite amplitude of the nonlinear structures that can propagate with phase velocity. However, this equation and some related equations such as KdV equation and modified KdV equation can be used indirectly to describe waves that propagate with the group velocity such as rogue waves and dark and bright solitons. For example, there are many researchers interested in studying waves in plasma physics who used this equation to describe the rouge waves in different plasma models [20][21][22][23].
Exact solutions to the GE and some related equations are set up by using various methods . Some solutions containing tanh and coth functions are proposed by the extended form of the tanh method [1]. Solitary wave and periodic solutions are constructed by aid of the projective Riccati equations. ese solutions have various terms including trigonometric or hyperbolic functions in rational forms. In most published papers, the authors focused on the integrable GE. However, in many physical models there are many effects cannot be ignored such as the collisions between the charged and neutral particles in plasma physics as well as the collisions between the charged particles themselves. Also, some of the external periodic forces can be impacted on the physical model under consideration. If these effects are considered in this case, we can get a nonintegrable evolution equation (forced damped GE). In this paper, we shall proceed to obtain some approximations to the following forced damped GE: where c(t) gives the coefficient of the damping term which arises due to different collisions that take place in plasma physics and F(t) indicates the external forces. To our knowledge there is no one attempt for studying this equation. us, the main goal of this work is to find some approximations for this equation using different analytical and numerical methods.

Novel Exact Solutions to the Integrable Gardner Equation
In order to obtain some traveling wave solutions (TWSs) to (1), we make the traveling wave transformation Integrating (3) once over ξ and denoting the constant of integration by d, we finally get the following Helmholtz-Duffing equation: Introducing the following ansatz, into (4) and after direct calculations, we obtain where the coefficients F j are given by F 0 � 2aA 3 k + 3A 2 bk + 6Aλ + 6 d , Solving the system S j � 0 will give us the values of (A, B, d, λ) as us, the cnoidal wave solution to GE (1) reads For letting m ⟶ 1, the soliton solution is recovered.
It is easy verified that the following solution is a rational solution: To obtain a traveling wave solution in terms of sn (ξ, m) and tanh(ξ), we make the following substitution: Inserting the ansatz (12) into GE (1) will give us where the coefficients S j (j � 0, 1, 2) are defined by e solution of the algebraic system S j � 0 (j � 0, 1, 2) gives the values of coefficients (A, B, λ):

e Scientific World Journal
Inserting the values (A, B, λ) given in (15) into ansatz (12), the solution in the form of Jacobi elliptic function sn (ξ, m) is obtained: Moreover, for m ⟶ 1, the shock wave solution is obtained: Remark 1. Equation (4) can be written in the following form: e general solution to (18) can be expressed in the terms of Weierstrass elliptic function as following: By substituting ansatz solution (19) into the Helmholtz-Duffing equation (18) and after tedious calculations, we get, e constants A and ξ 0 are determined from the initial In particular, this gives the general solution to both Duffing and Helmholtz equa- We see that the GE (1) admits solutions in terms of the Weierstrass elliptic function. In general, any partial differential equation (pde) can be converted to the ode (18) via the traveling wave transformation, admits TWSs in terms of Weierstrass elliptic function ℘. Also, Weierstrass function ℘ may be expressed in terms of the Jacobian cn function, which leads to cnoidal wave solutions to GE (1). e relation between the Jacobian and Weierstrass elliptic functions reads with On the other hand, with where ζ is the least in magnitude root of the following cubic equation: Note that for g 3 2 − 27g 2 3 � 0, we have only one real root, en, Let g 3 2 − 27g 2 3 ≠ 0, the discriminant of the cubic (26) reads In this case, our analysis depends on the sign of discriminant (29). Below, we will discuss three cases which depends on the sign of Δ.

Second
Case: Δ < 0. For Δ < 0, the cubic (26) has only one real root: 2.3. ird Case: Δ � 0. Here, suppose that g 3 2 − 27g 2 3 � 0, we have two real roots for cubic (26): It follows from (24) that the period of Weierstrass elliptic function reads where a is the greatest real root to the following cubic equation: Using the identity (24) and the approximation (37), the following approximation is obtained: where m is given by (25) and cos m (t): . e Jacobian elliptic functions cn and sn may be approximated by means of the following expressions: In the following Table 1, we can check the accuracy of the above obtained approximations.
We now give approximate expressions for inverse elliptic cosine and inverse elliptic cosine as follows: where − 1 ≤ t ≤ 1.
In Table 2, the numerical values of approximations (40) and (41) are displayed.

Approximate Analytical Solution to the Forced Damped Gardner Equation
Let us suppose that u � v ≡ v(x, t) is a solution to the GE (1): We seek for approximate solution to the forced damped GE (2) in the ansatz form: 4 e Scientific World Journal e functions f, g, h, and φ are to be determined later. Plugging ansatz (43) into (2) and taking the following value of zv/zt into account: we finally obtain e last expression suggests the choices: We will choose the functions f, g, h, and φ so that Solving system (46) and (47), we obtain the following solutions:  (37) and (38).   (40) and (41). (50) us, an approximate analytical solution to the damped and forced GE will be where φ � φ(t) can found from (50) and V(x, t) � v(x, t) is any analytic solution to the integrable GE (42).

FDM for Analyzing the Forced Damped Gardner Equation
To apply FDM for analyzing (2), we first write this equation in the following initial value problem (i.v.p.): where X i ≤ x ≤ X f and 0 ≤ t ≤ T f . e space-time domain are divided into subintervals with uniform size as where m and n are integer numbers. According to the FDM, the following derivatives formulas are introduced: 12Δt , In the case, when i < 0 or i > m or j < 0 or j > n, we define t) is the analytical approximation defined in Section 3. Now, we solve the following system of the nonlinear algebraic equations: with u i,0 � f(x i ) and i � 0, 1, 2, . . . , m and j � 0, 1, 2, . . . , n.
If we already solved system (55), then, we may construct an interpolation function with the data (x i , t j , u i,j ) for i � 0, 1, 2, . . . , m and j � 0, 1, 2, . . . , n. is interpolation function will represent the approximate numerical solution to the evolution equation. Another numerical solution may be obtained using the NDSolve Mathematica command.

Cubic Splines (Odd-order B-splines) for
Analyzing the Forced Damped Gardner Equation e general odd B-splines of (2r − 1)− order are defined as Here, χ s (x) � 1 for a + sh ≤ x < a + (s + 1)h and 0 otherwise and φ i ≡ φ i (x).
Note that for r � 2, we obtain the so-called cubic B-splines as follows: 6 e Scientific World Journal with Assuming that u(x, t) � n+1 k�i− 1 δ k (t)φ k (x) and with the help of Table 3, we get, ese formulas may be employed for solving PDEs like KdV, KdV-Burgers, MKdV, Gardner, and many third-order PDEs arising in different branches of science specially plasma physics. Now, in order to solve the i.v.p. (52), we must solve the following system of nonlinear odes: Table 3: Formulas for cubic splines at nodes. e Scientific World Journal where i � − 1, 0, . . . , n, n + 1, δ j (t) � 0 for j < − 1 or j > n + 1. We must choose the value of h � (b − a)/n in order to get the least residual error as possible.

Analysis and Discussion
We have obtained some analytical and numerical approximations the integrable GE (1) and nonintegrable forced damped GE (2). For analyzing the obtained solutions, we start by an exact solution to integrable GE (1). Some exact solutions to GE (1) such as the cnoidal wave solution given in (8) and the soliton solution given in (9) are introduced during the analysis the approximate solutions to the nonintegrable forced damped GE (2). Also, the following exact soliton solution is introduced to analyze (2).
Now, based on soliton solution (61) and according to the values (a, b, c, c(t), F 0 , ω, k) � (1, 1, 1, 0.2, 0.1, 4, 1), the profile of the approximate analytic soliton solution and the cubic B-splines soliton solution to the forced damped GE (2) is, respectively, presented in Figure 1. Also, the global residual errors in the whole domain for both approximate analytic solution and the numerical approximation using cubic B-splines for n � 20 are, respectively, estimated as L r � 0.630565 and L r � 0.757672. Note that the accuracy of the approximations depends on the values of physical parameters and the chosen exact solution.

Conclusion
In this work, some novel exact solutions to the integrable Gardner equation (GE) and approximations to the nonintegrable forced damped Gardner equations have been obtained. e most important obtained results can be briefly summarized in the following points: (1) In the first part, the integrable GE was reduced to the Helmholtz-Duffing equation using traveling wave transformation. After that, some exact solutions have been derived using the ansatz method. e obtained solutions have been derived in the form of Jacobi and Weierstrass elliptic functions. Moreover, the relation between Jacobi and Weierstrass elliptic functions has been presented. e obtained solutions can be recovered cnoidal waves, solitary waves, and shock waves to the GE. (2) In the second part, general formula for the approximate analytical solution to the forced damped GE has been derived in detail. is solution can be recovered from many nonlinear solutions that can be created and propagated in plasma physics. Based on this formula, the characteristics of many nonlinear structures in plasma physics such as solitary waves, shock waves, and cnoidal waves can be investigated. (3) In the third part, the evolution equation (the forced damped GE) has been analyzed using FDM in order to obtain an approximate numerical solution. (4) In the fourth-part, the cubic splines (Odd-order B-splines) were employed for analyzing the forced damped GE numerically.
Finally, the obtained solutions can help all researchers who are interested by studying the nonlinear structures in fluid mechanics, optical fiber, physics of plasmas, ocean and seas, and water tank waves.