Exact Periodic Wave Solutions for the Perturbed Boussinesq Equation with Power Law Nonlinearity

: In this paper, exact periodic wave solutions for the perturbed Boussinesq equation with power law nonlinearity are obtained for different nonlinear strengths n . When n = 1, the periodic traveling wave solutions can be found by the definition of the Jacobian elliptic function. When n ≥ 1, we construct a transformation to solve for the power law nonlinearity, and the periodic traveling wave solutions can be obtained by applying the extended trial equation method. In addition, we consider the limiting case where the periodicity of the periodic traveling wave solutions vanishes, and we obtain the soliton solution for n = 1. Numerical simulations show the periodicity of the solution for the perturbed Boussinesq equation.


Introduction
The Boussinesq equations are most commonly used to describe the propagation of small amplitude, plane, and long waves over water.The Boussinesq equations are perturbations of the linear wave equation that consider the small effects of nonlinearity and dispersion.In 1872, Joseph Boussinesq first derived the Boussinesq Equation [1] where µ(x, t) represents the non-perturbed surface, with x and t indicating the spatial and temporal coordinates, respectively.b 0 is the gravitational constant, and c 0 is the depth averaged over the surface.In non-dimensional units, the Boussinesq equations can be written as follows: From (1), it can be seen that both µ tt and µ xxxx have the same sign.The Boussinesq equations with this characteristic are often called the good Boussinesq equations, and equations with the opposite characteristic are called the bad Boussinesq equations, as follows: It is worth noting that the terms µ tt and µ xxxx have the opposite signs.
In this study, we consider the following perturbed Boussinesq equation with power law nonlinearity where u(x, t) denotes the waveform, n is the strength of nonlinearity, and d, a, b, β, and ρ are real-valued constants (−d 2 ̸ = β, b ̸ = ρ, and a ̸ = 0).Equation (3) was initially introduced by Anjan Biswas et al. [2].The coefficient β affects the dissipative term, while the coefficient ρ provides higher-order stability.The terms β and ρ represent perturbations, and their presence affects the positivity of u  [3,4].However, the situation is the opposite for the bad Boussinesq equation, which is almost always likely to be not well posed [5].
The good Boussinesq equation describes a two-dimensional non-rotational flow of a viscous fluid in a homogeneous rectangular channel.In recent years, it has been studied in many different aspects.Xu et al. studied the Cauchy problem of the solutions for the good Boussinesq equation with three different initial energy levels [6].Charlier et al. established the asymptotic behavior of the solutions for the good Boussinesq Equation [7].In [8], Luigi obtained numerical solutions for the good Boussinesq equation.The bad Boussinesq equation describes the two-dimensional flow of small amplitude shallow water waves.It is closely related to the dynamics of the anharmonic lattice in the Fermi-Pasta-Ulam (FPU) problem [9].In [5], the local generalized solutions of the bad Boussinesq equation have been shown under a mild condition for the initial data.The non-existence of weak solutions for the initial boundary value problem was investigated in [10].Dai et al. studied the explicit homoclinic orbit and the solitons for the bad Boussinesq equation with periodic boundary conditions [11].The blow-up solutions for the bad Boussinesq equation were found in [12].
Solving nonlinear partial differential equations (PDEs) is an important area of research within the field of nonlinear evolution equations.This branch of research is valuable for helping scientists comprehend complex physical phenomena.There are a large number of methods for obtaining exact solutions of nonlinear evolution Equations [13][14][15][16][17][18][19][20].In [13], the dynamical systems method was proposed for the modified Zakharov equations with a quantum correction.In order to obtain new types of solutions, Benjamin investigated the positive-operator method [14].The generalized Zakharov-Kuznetsov equation was solved using the F-expansion method in [15].The Gardner equation was solved using the Exp-function method by Wang et al. in [16].In [17], Hulya presented the trigonometric type of solutions to nonlinear Schr ödinger's equation using the (G ′ /G) expansion method.
The authors obtained the exact solutions of the nonlinear evolution equation utilizing the B äcklund method [18].Sivalingam et al. studied the fractional differential equations by the neural network scheme in [19,20].
In recent years, little research has been conducted on the solutions of the perturbed Boussinesq Equation (3) [2,21,22].Anjan Biswas obtained the solitary wave solutions for (3) with β = ρ = 0 by the ansatz method [2].The mapping method was employed to derive the exact 1-soliton solution and the singular solution for (3) with n = 1 [21].In [22], the authors also studied the soliton solutions for (3) with n = 1 using the modified auxiliary equation method.However, they only obtained solitary wave solutions for (3).No relevant work has been carried out on the periodic wave solutions for the perturbed Boussinesq Equation (3).As is well-known, the periodic wave solutions of nonlinear evolution equations are widely found in daily life, such as tidal phenomena, animal reproduction, food supply, plant growth, seasonal changes, and so on.The investigation of the periodic solutions to (3) not only extends the form of the exact traveling wave for the perturbed Boussinesq equation but also lays the foundation for researchers to study the qualitative analysis for the perturbed Boussinesq equation in the future, such as the stability.The main contributions of this paper are as follows: • The exact periodic wave solutions for the perturbed Boussinesq equation with power law nonlinearity are studied.To the best of the authors' knowledge, it is the first time to explore periodic wave solutions for the perturbed Boussinesq equation.Moreover, the periodic traveling wave solutions of the perturbed Boussinesq equation are obtained for general n, not just for a specific value of n [21,22].

•
Different from the existing periodic wave solutions [15][16][17] for nonlinear evolution equations, which are expressed in trigonometric functions [15,17] and exponential functions [16].The periodic solutions for the perturbed Boussinesq equation in this paper are expressed in terms of Jacobian elliptic functions.The Jacobian elliptic function is more suitable for engineering applications than other functions, such as the Duffing system, which is often used to describe oscillations in circuit systems [23].• Furthermore, we investigate the limiting case where the periodicity of the periodic traveling wave solution vanishes and derive the soliton solution with a singlepeaked waveform.
The remainder of this paper is arranged as follows.Section 2 introduces the definition of the Jacobian elliptic function and its related properties.Section 3 shows the periodic traveling wave solutions of the perturbed Boussinesq Equation (3) in two cases (n = 1 and n ≥ 1).The simulation results are presented to illustrate the periodicity of three types of Boussinesq equations in Section 4. The paper is concluded in Section 5.

Preliminaries
In this section, we introduce Jacobian elliptic functions.We refer the reader to the book Jaime [24] as the reference to this theory.There are three basic Jacobian elliptic functions, cn(u, h), dn(u, h), and sn(u, h), where h is the elliptic modulus.We start with the definition of the normal elliptic integral of the first kind where 0 < h 2 < 1 and ϕ = am(u, h) = am(u) is the Jacobi amplitude.We denote Further steps can be obtained: Based on the literature [24], the basic properties of the above three types of Jacobian elliptic functions are as follows: where h ′ = √ 1 − h 2 is a replenishment modulus and K(h) = F(π/2, h) is the complete elliptic integral of the first type.

Exact Periodic Wave Solutions for the Perturbed BE
We use the following traveling wave transformation: Inserting (4) into Equation ( 3) and integrate the equation with respect to ξ twice, we obtain with ϕ ′′ = ∂ 2 ϕ ∂ξ 2 , and where A ϕ is the integration constant.The nonlinear term ϕ 2n in ( 5) is difficult to derive directly, and we discuss it separately for the two cases n = 1 and n ≥ 1 using both direct and indirect methods.In the next sub-sections, we discuss the periodic solutions for two cases: n = 1 and n ≥ 1.

n = 1 in Equation (3)
In this sub-section, we discuss the case of the equation at n = 1.When n = 1, Equation ( 5) can be rewritten as follows: We multiply ( 6) by the integrating factor ϕ ′ , and we integrate once more for ξ where B ϕ is another constant of integration.This simplifies to Formula (8) can be rewritten as with , and α 3 be the non-zero roots of the polynomial F(ϕ), such that the polynomial is obtained.Without loss of generality, we assume that We denote φ = ϕ α 3 .Then, substituting it into (8) and (10), we obtain where η i = α i α 3 , i = 1, 2.Moreover, we define a new variable ψ, which satisfies φ = 1 + (η 2 − 1)sin 2 ψ.Substituting it into (12), after some more calculations, we obtain We can further rewrite (13) in the form of the following elliptic function definition where The left-hand side of ( 14) is the standard elliptic integral of the first kind.From this, we can immediately obtain From φ = ϕ α 3 and η i = α i α 3 , i = 1, 2, we have where Based on the properties of the Jacobian elliptic function sn 2 u + cn 2 u = 1, we derive where . According to cn(u + 2K(h)) = −cn(u), it follows that cn 2 has a fundamental period 2K(h).Therefore, the fundamental wave period of the cnoidal wave solution ϕ(ξ) is L n=1 with the value of Remark 1. Indeed, we consider that α 1 and α 2 converge to zero through positive values, h 2 → 1 − and α 3 → − 3(c 2 −d 2 −β)

2a
. In this limit, the elliptic function and its period are also simplified, sn(u, 1 − ) ∼ tanh u and K(h) → 1 + ∞ for h → 1 − .In this limit, the periodicity of the cnoidal wave is lost, and we obtain a waveform with a single hump and an "infinity period" of the form which is exactly the soliton solution of the perturbed Boussinesq Equation (3).Next, we consider the limiting case

n ≥ 1 in Equation (3)
In this section, we first introduce the extended trial equation method and then apply it to construct the solutions for the perturbed Boussinesq equation Equation (3).

Details of the Extended Trial Equation Method
Assume that the nonlinear development equation takes the following form: where u = u(x, t) is the unknown function, and G is the polynomial with respect to u and its partial derivatives.The main steps to obtain the exact solution of Equation ( 19) using the extended trial equation method are as follows [25]: Step 1. Substitute the following transformation into ( 19) This transformation turns (19) into the following form of an ordinary differential Equation (ODE): Step 2. Denote the trial equation in the form where the function θ satisfies the following trial differential equation where s i , i ∈ {0, 1, . . .δ} and ζ j , j ∈ {0, 1, . . ., p} are constants to be determined later.Using the relations ( 22) and ( 23), it can be deduced that where Φ(θ) and Υ(θ) are polynomials with respect to θ.
Step 3. Balancing the highest derivative term with the nonlinear term, the relations between σ, p, and δ can be derived.
Step 5. Equation ( 23) can be simplified to the elementary integral form where ξ 0 is the integration constant.As in the plane, the finite number of singularities of a planar dynamical system depends on the number of real roots of the polynomial F(θ), denoted by the discriminant of the real roots of F(θ) = 0. Classification of the roots of Φ(θ) using the complete discriminant system of polynomials can be completed, with the help of software (Maple2018.0 (WindowsX64), Maple Soft, Inc., Waterloo, ON, Canada), to solve (28) and classify the solutions of ( 21).The following Lemma 1 needs to be introduced.
Lemma 1.For a cubic polynomial F(θ) = θ 3 + d 2 θ 2 + d 1 θ + d 0 , the discriminant of the roots of the polynomial F(θ) is 3 ) 3 , and To determine which cases of F(θ) have three different real roots, there are three cases to be discussed: Case 1: ∆ = 0, D < 0, F(θ) = 0 has two distinct real roots with multiplicities one and two, respectively; that is, = 0 has only one real root with multiplicity three; that is, F(θ) = (θ − α) 3 ; Case 3: ∆ < 0, F(θ) = 0 has only one real root; that is, Remark 2. Lemma 1 gives some applications of polynomial complete discriminant systems to the classifications of the integral Equation (28).It is also presented that the form of the traveling wave solution is different in different cases.These solutions include trigonometric solutions, soliton solutions, periodic elliptic function solutions, etc.Using the extended trial equation method and Lemma 1, it is possible to solve not only the soliton solution of (3) but also the periodic solution.
Cases 1-3 show the conditions that need to be satisfied to obtain solitary wave solutions.However, the soliton solution of (3) has been solved using other methods [2,21,22].In this paper, we focus on periodic solutions, only case 4 needs to be satisfied, and no other solutions are discussed.Further explanation can be found in [26].

Application
From (5), the nonlinear term ϕ 2n is difficult to solve directly, which forces us to perform a transformation to solve the problem.First of all, the integration constant A ϕ = 0 for (5) can be assumed.Then, we multiply ϕ 4n−3 (ξ) on both sides of Equation (5).Applying the transformation and after some calculations, Equation ( 5) can be rewritten as Substituting ( 24) and ( 25) into Equation (31), and applying the balance principle different σ can be obtained by choosing different p and δ; therefore, Equation (32) has infinite solutions.In (32), we choose σ = 3, p = 0, and δ = 1.Then, where s 3 ̸ = 0, ζ 0 ̸ = 0. We substitute (33) into Equation (31) and denote each coefficient in the polynomial over θ to be zero.Furthermore, by solving the algebraic equations for these coefficients, we obtain , where ζ 0 , τ 0 , and τ 1 are parameters.Substituting these results into ( 23) and (28), where For T(θ) and q in Equation ( 36), the values are related to the parameters a, n, τ 1 , etc.If q > 0, we take T(θ); if q < 0, we take −T(θ).From Lemma 1, Now, T(θ) has three different roots: Combining this with (36), we have Then, we let θ = m 3 + (m 2 − m 3 )sin 2 φ.We can obtain where Furthermore, m 1 , m 2 , and m 3 are the roots of the polynomial equation T(θ) = 0.Then, we can obtain from (39) that According to sn(u + 2K(h)) = −sn(u), it follows that sn 2 has the fundamental period 2K(h).Therefore, the fundamental period of θ is L θ , with the value of Theorem 1.When the perturbed Boussinesq equation satisfies (37), Equation (3) has the following periodic wave solutions: where h = m 2 −m 3 m 1 −m 3 .
Remark 3. From sn 2 u + cn 2 u = 1, (42) can be rewritten in the form of a function about the cnoidal wave Remark 5.According to sn(u + 2K(h)) = −sn(u), the fundamental periods of sn and sn 2 are 4K(h) and 2K(h).Moreover, for n ≥ 1, the period of the solution of Equation (3) is determined by n.See [27] for details.

Numerical Simulations
In this section, we plot the simulation of the periodic wave solutions for the perturbed Boussinesq Equation (3) with two cases (n = 1 and n = 2).The nature of (3) depends on the values of the selected parameters d, β, b, and ρ.

Conclusions
Exact periodic wave solutions for the perturbed Boussinesq equation with power law nonlinearity are obtained for different nonlinear strengths n.Equation (3) contains three types of equations due to the differences in the positivity and negativity of the parameters, namely, the good Boussinesq equation, the bad Boussinesq equation, and neither a good nor bad Boussinesq equation.For the first time, the analytic periodic wave solutions of (3) were obtained using direct and indirect methods.For n = 1 in (3), the direct method was used to obtain the periodic solution.For n ≥ 1 in (3), the periodic solution can be found using the extended trial equation method.In addition, the simulation results are presented to illustrate the periodicity of three types of Boussinesq equations.
xx , and u xxxx .When d = 1, a = 1, b = 1, β = ρ = 0, and n = 1 in Equation (3), it represents the classical Boussinesq equation.This paper studies the Boussinesq equation with general n and non-zero β and ρ.On the left side of Equation (3), a(u 2n ) xx represents the nonlinear term.Equation (3) contains three types of equations: (1) when d 2 + β > 0 and b − ρ > 0, Equation (3) represents a good Boussinesq equation; (2) when d 2 + β > 0 and b − ρ < 0, Equation (3) represents a bad Boussinesq equation; (3) when d 2 + β < 0, Equation (3) represents neither a good nor bad Boussinesq equation.The perturbed Boussinesq equations have distinct mathematical properties due to the positive and negative nature of d 2 + β and b − ρ.The good Boussinesq equation has been widely studied in the existing literature, covering the local and global well posedness to the solutions of initial value problems

Figure 1
depicts the shape of the periodic solutions(17) of the three types of Boussinesq equations for n = 1.Figure2illustrates the shape of the periodic solutions (42) of the three types of Boussinesq equations for n = 2. (a) shows the shape of the periodic solution for the good Boussinesq equation when values for d, β, b, and ρ are selected to satisfy conditions d 2 + β > 0 and b − ρ > 0. (b) illustrates the shape of the periodic solution for the bad Boussinesq equation when values for d, β, b, and ρ are chosen to satisfy conditions d 2 + β > 0 and b − ρ < 0. Finally, (c) displays the shape of the periodic solution for neither a good nor bad Boussinesq equation when values for d and β are selected to satisfy the condition d 2 + β < 0.

)
Moreover, (44) (at least formally) includes two basic solutions of the perturbed Boussinesq Equation (3), which are obtained as approximate values of periodic solutions.
Remark 4. We note that (44) involves a large number of periodic solutions for Equation (3).Comparing (17) with (44), it can be realized that the solutions are of the same form but not exactly equal.It is worth noting that α i , i ∈ {1, 2, 3} in (17) and m j , j ∈ {1, 2, 3} in (44) are roots of the polynomials F(ϕ) and T(θ), respectively.