All Traveling Wave Exact Solutions of the Kawahara Equation Using the Complex Method

Abstract

In this article, we prove that the $⟨p,q⟩$ condition holds, first by using the Fuchs index of the complex Kawahara equation, and then proving that all meromorphic solutions of complex Kawahara equations belong to the class W. Moreover, the complex method is employed to get all meromorphic solutions of complex Kawahara equation and all traveling wave exact solutions of Kawahara equation. Our results reveal that all rational solutions $u_r(x+\nu t)$ and simply periodic solutions $u_{s,1}(x+\nu t)$ of Kawahara equation are solitary wave solutions, while simply periodic solutions $u_{s,2}(x+\nu t)$ are not real-valued. Finally, computer simulations are given to demonstrate the main results of this paper. At the same time, we believe that this method is a very effective and powerful method of looking for exact solutions to the mathematical physics equations, and the search process is simpler than other methods.

1. Introduction and Main Results

In this article, we study all traveling wave exact solutions of the Kawahara equation [1]:

$$u_t+au{u}_x+c{u}_{xxx}-d{u}_{xxxxx}=0,$$

(1)

where $a,c,d$ are constants, x and t represent spatial and time variables, $u{u}_x$ is the nonlinear disturbance term, and $u_{xxx}$ and $u_{xxxxx}$ are the dispersion terms of order three and five, respectively. The solutions of Equation (1) could be used to analyze and interpret a lot of nonlinear dispersive phenomena that can arise in optical fiber, ocean, plasma physics, etc. [2,3]. Given the importance of this equation, it remains the subject of study for many researchers. Many different techniques were devoted to studying various solutions of Equation (1) and its family [4,5,6,7,8,9,10,11,12,13]. For example, Kudryashov [14], obtained exact meromorphic solutions of the Kawahara equation using the Laurent series. Wazwaz [15], found some different solutions (compacton and solitons) to Equation (1) in terms of trigonometric functions. Khan [16], analyzed the Kawahara equation using the variational approach and derived new conditions for obtaining solitary wave solutions. Using the traveling wave ansatz, Baiswas [17], studied the generalized Kawahara equation and derived a solitary wave solution for the family of the Kawahara equation. A lot of effective methods are applied to study the exact solutions, which makes the research more abundant [18,19,20,21,22,23]. Wang [24] used ansatz method to derive the exact solitary wave solution for the generalized Korteweg–de Vries–Kawahara (GKdV-K) equation. Aiman Zara [25] studied numerical approximation of the modified Kawahara equation using the Kernel smoothing method. El-Tantawy [26] derived a set of novel exact and approximate analytic solutions to the family of the forced damped Kawahara equation (KE) using the ansatz method.

In recent years, Yuan et al. [27,28] introduced the $\langle p,q\rangle$ condition to study a class of constant-coefficient complex algebraic differential equations with dominant terms, and obtained all possible meromorphic solutions which belong to the class W. That is to say, the general representation of rational function solutions, finitely growing simply periodic solutions and elliptic solutions, and both rational function solutions and simply periodic solutions can be obtained by degenerating elliptic general solutions. Futhermore, a complex method for obtaining exact solutions in the mathematical physics is presented. Using the complex method [27,28], we prove that all meromorphic solutions of complex Kawahara equation belong to the class W, and then get all meromorphic solutions of complex Kawahara equation and all traveling wave exact solutions of Equation (1). This method is simpler than other methods in the process of finding solutions, and can also be applied to the solution of other nonlinear differential equations.

Substitute the traveling wave transform

$$u=u(x,t)=w(z),z=x+\nu t$$

(2)

into Equation (1); integrate it to get

$$C_0+\nu w+\frac{1}{2}a{w}^2+c{w}^{″}-d{w}^{(4)}=0,$$

(3)

where $C_0$ is an integral constant.

Multiply Equation (3) by $w^{\prime }$, and integrate it; thus, we get the complex ordinary differential Kawahara equation

$$w^{\prime }{w}^{‴}-\frac{1}{2}{(w^{″})}^2-A{w}^3+B{(w^{\prime })}^2+C{w}^2+Dw+E=0,$$

(4)

where $A=\frac{a}{6d},B=-\frac{c}{2d},C=-\frac{\nu }{2d},D=-\frac{C_0}{d},E=-\frac{C_1}{d}$, $C_0$ and $C_1$ are integral constants.

Now, we give the main results in our paper.

Theorem 1.

Let $A\ne 0$, then all meromorphic solutions $w(z)$ of Equation (4) belong to the class W. Here, class W consists of elliptic functions and their degenerations with the form $R(z)$ or $R(e^{\alpha z})$, $\alpha \in \mathbb{C}$, where R is a rational function.

In addition, the solutions of Equation (4) have the forms:

(1) All elliptic function solutions

$$\begin{gathered} w_d(z)=\frac{35}{2A}{[\frac{{\wp }^{\prime }(z)+M}{\wp (z)-N}]}^4-[\frac{140(\wp (z)+N)}{A}-\frac{70B}{39A}]{[\frac{{\wp }^{\prime }(z)+M}{\wp (z)-N}]}^2 \\ \hspace{1em}\hspace{1em}+\frac{280}{A}{[\wp (z)+N]}^2-\frac{280B}{39A}[\wp (z)+N]+\frac{C}{3A}-\frac{62B^2}{1521A}-\frac{28}{A}{g}_2, \end{gathered}$$

(5)

where $g_3=\frac{31B^3}{593190}-\frac{7B}{390}{g}_2,M^2=4N^3-g_{2}N-g_3,$

$$g_2^2-\frac{2B^2}{507}{g}_2+\frac{1457}{41384889}{B}^4-\frac{AD}{1932}-\frac{C^2}{5796}=0,$$

$$E=\frac{1}{A^2}[(\frac{1148}{69}AD+\frac{1148}{207}{C}^2-\frac{661248}{656903}{B}^4)g_2+\frac{23808}{4826809}{B}^6$$

$$+(\frac{20}{1521}AD+\frac{20}{4563}{C}^2)B^2-\frac{1}{3}ACD-\frac{2}{27}{C}^3],$$

D and N are constants.

(2) All simply periodic solutions

$$\begin{array}{c}{\displaystyle w_{s,1}(z)=\frac{2B^2}{169A}[35{coth}^4\frac{\sqrt{26B}i}{26}(z-z_0)-70{coth}^2\frac{\sqrt{26B}i}{26}(z-z_0)+23]+\frac{C}{3A},}\hfill \end{array}$$

(6)

and

$$\begin{gathered} {\displaystyle w_{s,2}(z)=\frac{B^2}{507A{\gamma }^2}[210{coth}^4\sqrt{\frac{B}{26\gamma }}(z-z_0)-(140\gamma +560){coth}^2\sqrt{\frac{B}{26\gamma }}(z-z_0)} \\ +166\gamma +352]+\frac{C}{3A}, \end{gathered}$$

(7)

where ${\displaystyle z_0\in \mathbf{C};D=\frac{1}{85683}\frac{5184B^4-28561C^2}{A},}$

$$E=\frac{1}{130323843}\frac{746496B^6-2628288B^{4}C+4826809C^3}{A^2}$$

in $w_{s,1}(z)$, and

$$D=\frac{1}{23905557}\frac{(217620\gamma -59724)B^4}{{\gamma }^{4}A}-\frac{C^2}{3A},$$

$$E=\frac{1}{A^2}[\frac{16B^6}{448893237}\frac{3001-7007\gamma }{{\gamma }^6}-\frac{B^{4}C}{71716671}\frac{217620\gamma -59724}{{\gamma }^4}+\frac{C^3}{27}],$$

γ is a root of $31z^2-31z+10=0$ in $w_{s,2}(z)$.

(3) All rational function solutions

$${\displaystyle w_r(z)=\frac{280}{A}\frac{1}{{(z-z_0)}^4}+\frac{C}{3A},}$$

(8)

where $B=0,D=-\frac{C^2}{3A},E=\frac{C^3}{27A^2}$ and $z_0\in \mathbf{C}$.

All traveling wave exact solutions of Equation (1) are obtained by substituting (2) into meromorphic solutions w of Equation (4). So, we get the following Theorem.

Theorem 2.

Let $ad\ne 0$, all traveling wave exact solutions $u(x,t)$ of Equation (1) have the following forms:

(1) All elliptic function solutions

$$\begin{array}{c}\hfill u_d(x+\nu t)=\frac{105d}{a}{[\frac{{\wp }^{\prime }(x+\nu t)+M}{\wp (x+\nu t)-N}]}^4-(\frac{840d(\wp (x+\nu t)+N)}{a}+\frac{70c}{13a}){[\frac{{\wp }^{\prime }(x+\nu t)+M}{\wp (x+\nu t)-N}]}^2+\\ \hfill \frac{1680d}{a}{[\wp (x+\nu t)+N]}^2+\frac{280c}{13a}[\wp (x+\nu t]+N)-\frac{168d}{a}{g}_2-\frac{31c^2}{507ad}-\frac{\nu }{a},\end{array}$$

(9)

where $M^3=4N^3-g_{2}N-g_3,g_3=-\frac{31c^3}{4745520d^3}+\frac{7c}{780d}{g}_2,$

$$g_2^2-\frac{c^2}{1014d^2}{g}_2+\frac{1457c^4}{662158224d^4}+\frac{a{C}_0}{11592d^2}-\frac{{\nu }^2}{23184d^2}=0,$$

and $C_0,N$ are constants.

(2) All simply periodic solutions

$$\begin{array}{c}\hfill u_{s,1}(x+\nu t)=\frac{3c^2}{169ad}[35{coth}^4(\frac{1}{2}\sqrt{\frac{c}{13d}}(x-x_0+\nu (t-t_0)))\\ \hfill -70{coth}^2(\frac{1}{2}\sqrt{\frac{c}{13d}}(x-x_0+\nu (t-t_0)))+23]-\frac{\nu }{a},\end{array}$$

(10)

and

$$\begin{array}{c}\hfill u_{s,2}(x+\nu t)=\frac{c^2}{338ad{\gamma }^2}[210{coth}^4(\frac{1}{2}\sqrt{-\frac{c}{13d\gamma }}(x-x_0+\nu (t-t_0)))\\ \hfill -(140\gamma +560){coth}^2(\frac{1}{2}\sqrt{-\frac{c}{13d\gamma }}(x-x_0+\nu (t-t_0)))\\ \hfill +166\gamma +352]-\frac{\nu }{a},\end{array}$$

(11)

where γ is a root of $31z^2-31z+10=0$.

(3) All rational function solutions

$$u_r(x+\nu t)=\frac{1680d}{a}\frac{1}{{(x-x_0+\nu (t-t_0))}^4}-\frac{\nu }{a},$$

(12)

where $c=0,x_0$ and $t_0$ are arbitrary real constants.

The rest of the article is organized in the following order: In Section 2, we will present the relevant lemmas and methodology. Section 3 gives the detailed proof process of Theorem 1 and concise proof method of Theorem 2. Section 4 illustrates our main results using computer simulations. In the last section, some conclusions are given.

2. Preliminary Lemmas and the Complex Method

We need some definitions and lemmas in order to prove Theorem 1 and present the complex method.

Set $k,n\in \mathbf{N}:=\{1,2,3,...\},r_j\in {\mathbf{N}}_0=\mathbf{N}\cup \{0\}$, $r=(r_0,r_1,...,r_k)$, $j=0,1,...,k$.

$$M_r[w](z):={[w(z)]}^{r_0}{[w^{\prime }(z)]}^{r_1}{[w^{″}(z)]}^{r_2}\cdots {[w^{(k)}(z)]}^{r_k}.$$

(13)

The degree of $M_r[w]$ defined by $d(r):=r_0+r_1+\cdots +r_k$.

Definition 1 ([28]).

A differential polynomial is defined by

$$P[w]:=\sum _{r\in \Lambda }{b}_r{M}_r[w],$$

(14)

where Λ is a finite index set, and $b_r$ are constants. $degP[w]:={max}_{r\in \Lambda }\{d(r)\}$ is called the degree of $P[w]$.

Consider the differential equation

$$E(z,w):=P[w]-a{w}^n=0,$$

(15)

where $a\ne 0$ is a constant.

The dominant part of $E(z,w)$ composes all dominant terms which can determine the multiplicity q of w in $E(z,w)$, and is denoted by $\widehat{E}(z,w)$. $D(q)$ and $D_r(q)$ represent the multiplicity of pole of each term in $\widehat{E}(z,w)$ and the multiplicity of pole of each monomial in $E(z,w)-\widehat{E}(z,w)$, respectively.

Obviously

$$D_r(q)=qd(r)+r_1+2r_2+\cdots +k{r}_k<D(q).$$

(16)

Definition 2 ([28]).

The derivative of $\widehat{E}(z,w)$ with respect to w can be calculated by the following formula, for any $\chi ,$

$${\widehat{E}}^{\prime }(z,w)\chi =\underset{\lambda \to 0}{lim}\frac{\widehat{E}(z,w+\lambda \chi )-\widehat{E}(z,w)}{\lambda }.$$

(17)

Definition 3 ([28]).

Substituting Laurent series

$$w(z)=\sum _{l=-q}^{\infty }{c}_l{z}^l,$$

(18)

into Equation (15), where $c_{-q}\ne 0,q>0$. Then, we can get p different principle

$$\sum _{l=-q}^{-1}{c}_l{z}^l$$

with pole of multiplicity q at $z=0$, Equation (15) is said to satisfy weak $\langle p,q\rangle$ condition. If Equation (15) has p different meromorphic solutions with pole of multiplicity q at $z=0$, Equation (15) satisfies $\langle p,q\rangle$ condition.

Definition 4 ([28]).

Let $T_1,T_2$ be two given complex numbers, such that $Im\frac{T_1}{T_2}>0,L=L[2T_1,2T_2]$ is discrete subset $L[2T_1,2T_2]=\{T|T=2m{T}_1+2n{T}_2,m,n\in \mathbb{Z}\}$, which is isomorphic to $\mathbb{Z}\times \mathbb{Z}$. The discriminant $\Delta =\Delta (c_1,c_2):=c_1^3-27c_2^2$ and

$$s_n=s_n(L):=\sum _{T\in L\backslash \{0\}}\frac{1}{T^n}.$$

(19)

Weierstrass elliptic function $\wp (z):=\wp (z,g_2,g_3)$ is a meromorphic function with double periods $2T_1,2T_2$, which satisfies the following formula

$${({\wp }^{\prime }(z))}^2=4\wp {(z)}^3-g_2\wp (z)-g_3,$$

(20)

where $g_2=60s_4,g_3=140s_6$ and $\Delta (g_2,g_3)\ne 0.$

Definition 5 ([29]).

The Fuchs index of Equation (15) can be defined as the root of the equation

$$P(i)=\underset{z\to 0}{lim}{z}^{-i+D(q)}{\widehat{E}}^{\prime }(z,c_{-q}{z}^{-q})z^{i-q}=0$$

(21)

Lemma 1 ([28,30,31,32]).

Set $p,q,m,n\in \mathbf{N},degP[w]<n$. If Equation (15) satisfies $\langle p,q\rangle$ condition, all meromorphic solutions $w(z)$ of Equation (15) belong to class W.

Any elliptic function solution with pole at $z=0$ is given in the form

$$\begin{gathered} {\displaystyle w(z)=\sum _{i=1}^{m-1}\sum _{j=2}^q\frac{{(-1)}^j{c}_{-ij}}{(j-1)!}\frac{d^{j-2}}{d{z}^{j-2}}(\frac{1}{4}{[\frac{{\wp }^{\prime }(z)+M_i}{\wp (z)-N_i}]}^2-\wp (z))} \\ {\displaystyle +\sum _{i=1}^{m-1}\frac{c_{-i1}}{2}\frac{{\wp }^{\prime }(z)+M_i}{\wp (z)-N_i}+\sum _{j=2}^q\frac{{(-1)}^j{c}_{-mj}}{(j-1)!}\frac{d^{j-2}}{d{z}^{j-2}}\wp (z)+k_0,} \end{gathered}$$

(22)

where $c_{-ij}$ can be determined by (18), $c_{-ij}$ and $k_0$ are constants. $M_i^2=4N_i^3-g_2{N}_i-g_3$ and ${\displaystyle \sum _{i=1}^m{c}_{-i1}=0.}$

Any rational function solution $w:=R(z)$ is expressed as

$${\displaystyle R(z)=\sum _{i=1}^m\sum _{j=1}^q\frac{c_{ij}}{{(z-z_i)}^j}+k_0,}$$

(23)

which has $m(\le p)$ different poles of multiplicity q.

Any simply periodic solution $w:=R(\xi )$ is a rational function of $\xi =e^{\alpha z}(\alpha \in \mathbb{C})$ and can be given in the form

$${\displaystyle R(\xi )=\sum _{i=1}^m\sum _{j=1}^q\frac{c_{ij}}{{(\xi -{\xi }_i)}^j}+k_0,}$$

(24)

which has $m(\le p)$ different poles of multiplicity q.

Lemma 2 ([33,34]).

Weierstrass elliptic functions $\wp (z)$ have the addition formula and two successive degeneracies, as shown below:

(I) Addition formula

$$\wp (z-z_0)=-\wp (z)-\wp (z_0)+\frac{1}{4}{[\frac{{\wp }^{\prime }(z)+{\wp }^{\prime }(z_0)}{\wp (z)-\wp (z_0)}]}^2.$$

(25)

(II) If $\Delta (g_2,g_3)=0$, Weierstrass elliptic functions degenerate to simply periodic functions, which can be expressed as

$$\wp (z,3{\delta }^2,-{\delta }^3)=2\delta -\frac{3\delta }{2}{coth}^2\sqrt{\frac{3\delta }{2}}z.$$

(26)

(III) If $g_2=g_3=0$, Weierstrass elliptic functions degenerate to rational functions of z, which can be expressed as

$$\wp (z,0,0)=\frac{1}{z^2}.$$

(27)

Next, we give the complex method.

• Substituting the transform $T:u(x,t)\to w(z),(x,t)\to z$ into a given PDE yields a nonlinear ODE: Equation (4) here.
• Insert (18) into Equation (4) here to determine that weak $\langle p,q\rangle$ condition holds.
• By (22)–(24), we obtain all meromorphic solutions $w(z)$ of Equation (4) here with pole at $z=0$.
• Get all meromorphic solutions $w(z-z_0)$ by Lemmas 1 and 2.
• Inserting the inverse transform $T^{-1}$ into $w(z-z_0)$, we obtain all exact solutions $u(x,t)$ of the given partial differential equation.

3. Proof of Theorem 1

From balance the order of the poles in Equation (4), yields

$$(q+1)+(q+3)=2(q+2)=3q,$$

and we get $q=4$.

Substitute (18) into Equation (4),and set the cofficients to zero, we have

$$c_{-4}=\frac{280}{A},c_{-3}=0,c_{-2}=\frac{280B}{39A},c_{-1}=0,$$

$$c_0=-\frac{62B^2-507C}{1521A},c_1=0,c_2=\frac{62B^3}{59319A},c_3=0,$$

$$c_4=\frac{1}{425673144}\frac{-11780B^4+257049AD+85683C^2}{A},\cdots .$$

(28)

and then determine ${\displaystyle p=1.}$ Therefore, Equation (4) is said to satisfy weak $\langle 1,4\rangle$ condition. We also get that all meromorphic solutions of Equation (4) belong to W if Equation (4) satisfies $\langle 1,4\rangle$ condition by Lemma 1.

In fact, since (17) and $\widehat{E}(z,w)=w^{\prime }{w}^{‴}-\frac{1}{2}{(w^{″})}^2-A{w}^3$ of Equation (4), we have $D(4)=12$,

$$D_r(4)<12,A{c}_{-4}=280,$$

(29)

and

$$\begin{gathered} {\widehat{E}}^{\prime }(x,c_{-4}{x}^{-4})x^{i-4}=({(c_{-4}{x}^{-4})}^{\prime }\frac{d^3}{d{x}^3}+{(c_{-4}{x}^{-4})}^{‴}\frac{d}{dx}-{(c_{-4}{x}^{-4})}^{″}\frac{d^2}{d{x}^2} \\ -3A{(c_{-4}{x}^{-4})}^2)x^{i-4} \\ =-\left(4(i-1)(i-4)(i-5)+120i+360\right)c_{-4}{x}^{i-12}. \end{gathered}$$

(30)

Thus, from (29), (30) and (21), we infer that the Fuchs index of Equation (4) are zeros of the function

$$P(i)=-\left(4(i-1)(i-4)(i-5)+120i+360\right)c_{-4}.$$

(31)

It is easy to prove that $f(x):=4(x-1)(x-4)(x-5)+120x+360>0$ if $x>-1.$ Hence, if $i\in \mathbf{N}$, then $P(i)\ne 0$. That is to say, the Fuchs index of Equation (4) cannot be a positive integer. This implies that (ref. [29], p. 90): The principle part of the Laurent series of $w(z)$ determines the whole Laurent series of $w(z)$. Therefore, weak $\langle 1,4\rangle$ condition implies that $\langle 1,4\rangle$ condition holds.

From, (23), we get the indeterminant form of rational function solutions of Equation (4) at pole $z=0$

$${\displaystyle R_1(z)=\frac{c_{-4}}{z^4}+\frac{c_{-2}}{z^2}+k_0=\frac{280}{A}\frac{1}{z^4}+\frac{280B}{39A}\frac{1}{z^2}+k_0.}$$

(32)

Inserting (32) into Equation (4), and then setting the coefficients to zero, we get the system of Equations (1) which are given in Appendix A. Computing the system of Equations (1), we get

$$k_0=\frac{C}{3A},B=0,D=-\frac{C^2}{3A},E=\frac{C^3}{27A^2}.$$

So we find that the rational function solutions of Equation (4) are

$${\displaystyle w_{r0}(z)=\frac{280}{A}\frac{1}{z^4}+\frac{C}{3A},}$$

(33)

at pole $z=0$, here $B=0,D=-\frac{C^2}{3A},E=\frac{C^3}{27A^2}$.

Hence, all rational function solutions of Equation (4) are

$${\displaystyle w_r(z)=\frac{280}{A}\frac{1}{{(z-z_0)}^4}+\frac{C}{3A},}$$

(34)

where $B=0,D=-\frac{C^2}{3A},E=\frac{C^3}{27A^2}$ and $z_0\in \mathbf{C}$.

By (22) and (26), we get the indeterminant form of simply periodic solutions of Equation (4) at pole $z=0$

$$\begin{gathered} w_{s0}(z)=\frac{1}{6}{c}_{-4}{\wp }^{″}(z,3{\delta }^2,-{\delta }^3)+c_{-2}\wp (z,3{\delta }^2,-{\delta }^3)+k_0 \\ \hspace{1em}\hspace{1em}=\frac{280}{6A}{\wp }^{″}(z,3{\delta }^2,-{\delta }^3)+\frac{280B}{39A}\wp (z,3{\delta }^2,-{\delta }^3)+k_0, \end{gathered}$$

(35)

where $\wp (z,3{\delta }^2,-{\delta }^3)=2\delta -\frac{3}{2}\delta {coth}^2\sqrt{\frac{3\delta }{2}}z.$

Expanding $\wp (z,3{\delta }^2,-{\delta }^3)$ at $z_0=0$, we have

$$\wp (z,3{\delta }^2,-{\delta }^3)=\frac{1}{z^2}+\frac{3{\delta }^2}{20}{z}^2-\frac{{\delta }^3}{28}{z}^4+o(z^4)$$

(36)

Substituting (36) into (35), we get

$$\begin{array}{c}\hfill w(z)=\frac{280}{A}\frac{1}{z^4}+\frac{280B}{39A}\frac{1}{z^2}+\frac{14{\delta }^2}{A}+k_0+(\frac{14B{\delta }^2}{13A}-\frac{20{\delta }^3}{A})z^2\\ \hfill +(-\frac{10B{\delta }^3}{39A}+\frac{21{\delta }^4}{2A})z^4+o(z^4)\end{array}$$

(37)

Comparing coefficients of (37) and (28), we obtain the system of Equations (2) which are showed in Appendix A. Solving the system of Equations (2), we have

$$\delta =-\frac{B}{39},k_0=-\frac{76B^2}{1521A}+\frac{C}{3A},D=\frac{1}{85683}\frac{5184B^4-28561C^2}{A},$$

(38)

and

$$\delta =\frac{B}{39\gamma },k_0=\frac{B^2}{1521A}\frac{6-62\gamma }{{\gamma }^2}+\frac{C}{3A},$$

$$D=\frac{1}{23905557}\frac{(217620\gamma -59724)B^4}{{\gamma }^{4}A}-\frac{C^2}{3A}..$$

(39)

Substituting (35), (38) and (39) into Equation (4), from the correlation of coefficients we can obtain, respectively,

$$E=\frac{1}{130323843}\frac{746496B^6-2628288B^{4}C+4826809C^3}{A^2}$$

and

$$E=\frac{1}{A^2}[\frac{16B^6}{448893237}\frac{3001-7007\gamma }{{\gamma }^6}-\frac{B^{4}C}{71716671}\frac{217620\gamma -59724}{{\gamma }^4}+\frac{C^3}{27}],$$

where $\gamma$ is a root of $31z^2-31z+10=0.$

Hence, we obtain the simply periodic solutions of Equation (4) at pole $z=0$ are

$$\begin{array}{c}{\displaystyle w_{s,1}(z)=\frac{2B^2}{169A}[35{coth}^4\frac{\sqrt{26B}i}{26}z-70{coth}^2\frac{\sqrt{26B}i}{26}z+23]+\frac{C}{3A},}\hfill \end{array}$$

(40)

and

$$\begin{array}{c}{\displaystyle w_{s,2}(z)=\frac{B^2}{507A{\gamma }^2}[210{coth}^4\sqrt{\frac{B}{26\gamma }}z-(140\gamma +560){coth}^2\sqrt{\frac{B}{26\gamma }}z}\hfill \\ +166\gamma +352]+\frac{C}{3A},\hfill \end{array}$$

(41)

where $\delta =-\frac{B}{39},k_0=-\frac{76B^2}{1521A}+\frac{C}{3A},$ $D=\frac{1}{85683}\frac{5184B^4-28561C^2}{A},$

$$E=\frac{1}{130323843}\frac{746496B^6-2628288B^{4}C+4826809C^3}{A^2}$$

in $w_{s,1}(z)$, and $D=\frac{1}{23905557}\frac{(217620\gamma -59724)B^4}{{\gamma }^{4}A}-\frac{C^2}{3A},$

$$E=\frac{1}{A^2}[\frac{16B^6}{448893237}\frac{3001-7007\gamma }{{\gamma }^6}-\frac{B^{4}C}{71716671}\frac{217620\gamma -59724}{{\gamma }^4}+\frac{C^3}{27}],$$

$\gamma$ is a root of $31z^2-31z+10=0$ in $w_{s,2}(z)$.

Furthermore, all simply periodic solutions of Equation (4) are given by

$$\begin{array}{c}{\displaystyle w_{s,1}(z)=\frac{2B^2}{169A}[35{coth}^4\frac{\sqrt{26B}i}{26}(z-z_0)-70{coth}^2\frac{\sqrt{26B}i}{26}(z-z_0)+23]+\frac{C}{3A},}\hfill \end{array}$$

(42)

and

$$\begin{array}{c}{\displaystyle w_{s,2}(z)=\frac{B^2}{507A{\gamma }^2}[210{coth}^4\sqrt{\frac{B}{26\gamma }}(z-z_0)-(140\gamma +560){coth}^2\sqrt{\frac{B}{26\gamma }}(z-z_0)}\hfill \\ +166\gamma +352]+\frac{C}{3A},\hfill \end{array}$$

(43)

where ${\displaystyle z_0\in \mathbf{C},}$ $D=\frac{1}{85683}\frac{5184B^4-28561C^2}{A},$

$$E=\frac{1}{130323843}\frac{746496B^6-2628288B^{4}C+4826809C^3}{A^2}$$

in $w_{s,1}(z)$, and $D=\frac{1}{23905557}\frac{(217620\gamma -59724)B^4}{{\gamma }^{4}A}-\frac{C^2}{3A},$

$$E=\frac{1}{A^2}[\frac{16B^6}{448893237}\frac{3001-7007\gamma }{{\gamma }^6}-\frac{B^{4}C}{71716671}\frac{217620\gamma -59724}{{\gamma }^4}+\frac{C^3}{27}],$$

$\gamma$ is a root of $31z^2-31z+10=0$ in $w_{s,2}(z)$.

From (22), we infer the indeterminant forms of elliptic function solutions of Equation (4) are

$$\begin{gathered} w_{d0}(z)=\frac{c_{-4}}{6}{\wp }^{″}(z)+c_{-2}\wp (z)+k_0 \\ \hspace{1em}\hspace{1em}=\frac{280}{6A}{\wp }^{″}+\frac{280B}{39A}\wp +k_0. \end{gathered}$$

(44)

with pole at $z=0$.

Expanding $\wp (z,g_2,g_3)$ at $z=0$, we have

$$\wp (z,g_2,g_3)=\frac{1}{z^2}+\frac{g_2}{20}{z}^2+\frac{g_3}{28}{z}^4+o(z^4).$$

(45)

Put (45) into (44) and apply (20), we get

$$\begin{array}{c}\hfill w(z)=\frac{280}{A}\frac{1}{z^4}+\frac{280B}{39A}\frac{1}{z^2}+\frac{14g_2}{3A}+k_0+(\frac{14B{g}_2}{39A}+\frac{20g_3}{A})z^2\\ \hfill +(\frac{10B{g}_3}{39A}+\frac{7g_2^2}{6A})z^4+o(z^4).\end{array}$$

(46)

Comparing coefficients of $w(z)$ and (28), we obtain the system of Equations (3). which are shown in the Appendix. Computing the system of Equations (3), we derive

$$k_0=\frac{C}{3A}-\frac{14}{3A}{g}_2-\frac{62B^2}{1521A},g_3=\frac{31B^3}{593190}-\frac{7B}{390}{g}_2.$$

(47)

and $g_2$ satisfies

$$g_2^2-\frac{2B^2}{507}{g}_2+\frac{1457B^4}{41384889}-\frac{AD}{1932}-\frac{C^2}{5796}=0.$$

Substituting (44) and (47) into Equation (4), from the correlation of coefficients, we can get

$$\begin{array}{c}\hfill E=\frac{1}{A^2}\{(\frac{1148}{69}AD+\frac{1148}{207}{C}^2-\frac{661248}{656903}{B}^4)g_2+\frac{23808}{4826809}{B}^6\\ \hfill +(\frac{20}{1521}AD+\frac{20}{4563}{C}^2)B^2-\frac{1}{3}ACD-\frac{2}{27}{C}^3\},\end{array}$$

then we have

$$w_{d0}(z)=\frac{280}{A}{\wp }^2(z)+\frac{280B}{39A}\wp (z)+\frac{C}{3A}-\frac{62B^2}{1521A}-\frac{28}{A}{g}_2.$$

(48)

Thus, all elliptic function solutions of Equation (4) are

$$w_d(z)=\frac{280}{A}{\wp }^2(z-z_0)+\frac{280B}{39A}\wp (z-z_0)+\frac{C}{3A}-\frac{62B^2}{1521A}-\frac{28}{A}{g}_2,$$

(49)

where $z_0\in \mathbf{C},g_3=\frac{31B^3}{593190}-\frac{7B}{390}{g}_2,$

$$\begin{array}{c}\hfill E=\frac{1}{A^2}\{(\frac{1148}{69}AD+\frac{1148}{207}{C}^2-\frac{661248}{656903}{B}^4)g_2+\frac{23808}{4826809}{B}^6\\ \hfill +(\frac{20}{1521}AD+\frac{20}{4563}{C}^2)B^2-\frac{1}{3}ACD-\frac{2}{27}{C}^3\},\end{array}$$

and $g_2$ satisfies

$$g_2^2-\frac{2B^2}{507}{g}_2+\frac{1457B^4}{41384889}-\frac{AD}{1932}-\frac{C^2}{5796}=0.$$

By using the addition formula, we can get another representation of it as

$$\begin{array}{c}\hfill w_d(z)=\frac{35}{2A}{[\frac{{\wp }^{\prime }(z)+M}{\wp (z)-N}]}^4-[\frac{140(\wp (z)+N)}{A}-\frac{70B}{39A}]{[\frac{{\wp }^{\prime }(z)+M}{\wp (z)-N}]}^2\\ \hfill +\frac{280}{A}{[\wp (z)+N]}^2-\frac{280B}{39A}[\wp (z)+N]+\frac{C}{3A}-\frac{62B^2}{1521A}-\frac{28}{A}{g}_2,\end{array}$$

(50)

where

$$M^2=4N^3-g_{2}N-g_3,g_3=\frac{31B^3}{593190}-\frac{7B}{390}{g}_2,$$

$$g_2^2-\frac{2B^2}{507}{g}_2+\frac{1457}{41384889}{B}^4-\frac{AD}{1932}-\frac{C^2}{5796}=0,$$

$$E=\frac{1}{A^2}[(\frac{1148}{69}AD+\frac{1148}{207}{C}^2-\frac{661248}{656903}{B}^4)g_2+\frac{23808}{4826809}{B}^6$$

$$+(\frac{20}{1521}AD+\frac{20}{4563}{C}^2)B^2-\frac{1}{3}ACD-\frac{2}{27}{C}^3],$$

D and N are constants.

So far, we have completed the proof of Theorem 1.

Substituting (2) into all meromorphic solutions $w(z)$ of Equation (4), we obtained Theorem 2. According to the theorem in ref. [35], we can get the following corollary.

Corollary. All rational solutions $u_r(x+\nu t)$ and simply periodic solutions $u_{s,1}(x+\nu t)$ of Equation (1) are real valued, while simply periodic solutions $u_{s,2}(x+\nu t)$ are not real valued.

4. Computer Simulations

This subsection will show our results through computer simulations of $u_r(x+\nu t)$ and $u_{s,1}(x+\nu t)$, as demonstrated in the following figures.

• By applying the complex method, we are able to achieve the rational solution $u_r(x+\nu t)$ of Equation (4). Figure 1 describes the 3D graphs of solution $u_r(x+\nu t)$ for $a=1$, $\nu =1$, and $d=\frac{1}{1680}$ within the interval $-5\le x,t\le 5$. Figure 2 shows the 2D graphs of solution $u_r(x+\nu t)$ for $a=1$, $\nu =1$, and $d=\frac{1}{1680}$ within the interval $-10\le x\le 10$ when $t=0$. It could be observed that they have one generation pole, which is shown by Figure 1 and Figure 2.
• By employing the complex method, we are able to obtain the simply periodic solutions $u_{s,1}(x+\nu t)$ and $u_{s,2}(x+\nu t)$ of Equation (4). Figure 3 shows the 3D graphs of solution $u_{s,1}(x+\nu t)$ for $a=3$, $c=13$, $d=1$, and $\nu =3$ within the interval $-2\pi \le x,t\le 2\pi$. Figure 4 describes the 2D graphs of solution $u_{s,1}(x+\nu t)$ for $a=3$, $c=13$, $d=1$, and $\nu =3$ within the interval $-2\pi \le x\le 2\pi$ when $t=0$.

5. Conclusions

In this paper, we are the first to utilize the complex method to prove that all meromorphic solutions of complex Kawahara equation belong to the class W, and then we get all meromorphic solutions of complex Kawahara equation and all traveling wave exact solutions of Equation (1). In addition, we find that all rational solutions $u_r(x+\nu t)$ and simply periodic solutions $u_{s,1}(x+\nu t)$ of Equation (1) are solitary wave solutions, which could be used to analyze and interpret a lot of nonlinear dispersive phenomena. This research enriches the methods of solving differential equations. Our results also reveal that the complex method of looking for traveling wave exact solutions is general and feasible, and can be applied to other nonlinear partial differential equations.