Abstract

We construct new exact traveling wave solutions involving free parameters of the nonlinear reaction diffusion equation by using the improved (πΊξ…ž/𝐺)-expansion method. The second-order linear ordinary differential equation with constant coefficients is used in this method. The obtained solutions are presented by the hyperbolic and the trigonometric functions. The solutions become in special functional form when the parameters take particular values. It is important to reveal that our solutions are in good agreement with the existing results.

1. Introduction

Nonlinear evolution equations (NLEEs) describe many problems of solid state physics, nonlinear optics, plasma physics, fluid mechanics, population dynamics and many others which arise in mathematical biology, engineering sciences and other technical arena. In recent years, several methods have been developed to obtain traveling wave solutions for many NLEEs, such as the theta function method [1], the Jacobi elliptic function expansion method [2], the Hirota’s bilinear transformation method [3], the 𝐹-expansion method [4], the Backlund transformation method [5, 6], the generalized Riccati equation method [7, 8], the sub-ODE method [9], the homogeneous balance method [10, 11], the tanh-coth method [12–14], the sine-cosine method [15], the first integral method [16], the Cole-Hopf transformation method [17], the Exp-function method [18–25], and others [26–43].

Recently, Wang et al. [44] presented the (πΊξ…ž/𝐺)-expansion method and implemented to four well-established equations for constructing traveling wave solutions. In this method, the second-order linear ordinary differential equation (ODE) πΊξ…žξ…ž+πœ†πΊξ…ž+πœ‡πΊ=0 is used, where πœ† and πœ‡ are arbitrary constants. Afterwards, many researchers used this method to many nonlinear partial differential equations and obtained many new exact traveling wave solutions. For instance, Malik et al. [45] applied the (πΊξ…ž/𝐺)-expansion method for getting traveling wave solutions of some nonlinear partial differential equations. Bekir [46] concerned about this method to study nonlinear evolution equations for constructing wave solutions. Zayed [47] investigated the higher-dimensional nonlinear evolution equations by using the same method to get solutions. In [48], Naher et al. implemented the method for constructing abundant traveling wave solutions of the Caudrey-Dodd-Gibbon equation. Lately, Hayek [49] extended the method called extended (πΊξ…ž/𝐺)-expansion method to obtain exact analytical solutions to the KdV Burgers equations with power-law nonlinearity whilst Guo and Zhou [50] expand the method and applied to the Whitham-Broer-Kaup-Like equations and Coupled Hirota-Satsuma KdV equations to construct traveling wave solutions. Zayed and Al-Joudi [51] concerned about the method to find solutions of the NLPDEs in mathematical physics and so on.

More recently, Zhang et al. [52] extended the method which is called the improved (πΊξ…ž/𝐺)-expansion method for constructing abundant traveling wave solutions of the nonlinear evolution equations. Then, many researchers implemented the method to construct exact solutions. For example, Hamad et al. [53] solved the higher-dimensional potential YTSF equation by using this powerful and useful method for getting many new exact solutions. In [54], Nofel et al. investigated the higher-order KdV equation via the same method while Zhao et al. [55] applied this method to obtain traveling wave solutions for the variant Boussinesq equations. Tao and Xia [56] executed the method for searching exact solutions of the (3+1)-dimensional KdV equation and so on.

Many researchers studied the nonlinear reaction diffusion equation to obtain traveling wave solutions by using different methods. For instance, Zayed and Gepreel [57] used the (πΊξ…ž/𝐺)-expansion method to solve this equation. To the best of our knowledge, the nonlinear reaction diffusion equation is not investigated by using the improved (πΊξ…ž/𝐺)-expansion method.

In this paper, we apply the improved (πΊξ…ž/𝐺)-expansion method to construct new exact traveling wave solutions of the nonlinear reaction diffusion equation which is very important equation in mathematical biology.

2. Explanation of the Improved (πΊξ…ž/𝐺)-Expansion Method

Suppose the general nonlinear partial differential equation:𝑃𝑒,𝑒𝑑,𝑒π‘₯,𝑒π‘₯𝑑,𝑒𝑑𝑑,𝑒π‘₯π‘₯ξ€Έ,…=0,(2.1) where 𝑒=𝑒(π‘₯,𝑑) is an unknown function. 𝑃 is a polynomial in 𝑒=𝑒(π‘₯,𝑑) and the subscripts indicate the partial derivatives.

The main steps of the improved (πΊξ…ž/𝐺)-expansion method [52] are as follows.

Step 1. Consider the traveling wave variable:𝑒(π‘₯,𝑑)=𝑣(πœ‚),πœ‚=π‘₯βˆ’π‘‰π‘‘,(2.2) where 𝑉 is the speed of the traveling wave. Using (2.2), (2.1) is converted into an ordinary differential equation for 𝑒=𝑣(πœ‚): 𝑄𝑣,π‘£ξ…ž,π‘£ξ…žξ…žξ€Έ,…=0,(2.3) where the superscripts stand for the ordinary derivatives with respect to πœ‚.

Step 2. Suppose that the traveling wave solution of (2.3) can be presented in the following form [52]: 𝑣(πœ‚)=π‘šξ“π‘—=βˆ’π‘šπ‘Žπ‘—ξ‚΅πΊξ…žπΊξ‚Άπ‘—,(2.4) where 𝐺=𝐺(πœ‚) satisfies the second-order linear ODE: πΊξ…žξ…ž+πœ†πΊξ…ž+πœ‡πΊ=0,(2.5) where π‘Žπ‘—(𝑗=0,Β±1,Β±2,…,Β±π‘š),πœ†, and πœ‡ are constants.

Step 3. To determine the integer π‘š, substitute (2.4) along with (2.5) into (2.3) and then take the homogeneous balance between the highest-order derivatives and the highest-order nonlinear terms appearing in (2.3).

Step 4. Substitute (2.4) together with (2.5) into (2.3) with the value of π‘š obtained in Step 3. Equating the coefficients of (πΊξ…ž/𝐺)π‘Ÿ,(π‘Ÿ=0,Β±1,Β±2,…), then setting each coefficient to zero, yields a set of algebraic equations for π‘Žπ‘—(𝑗=0,Β±1,Β±2,…,Β±π‘š),𝑉,πœ†, and πœ‡.

Step 5. Solve the system of algebraic equations with the aid of commercial software Maple and we obtain values for π‘Žπ‘—(𝑗=0,Β±1,Β±2,…,Β±π‘š),𝑉,πœ†, and πœ‡. Then, substituting obtained values in (2.4) along with (2.5) with the value of π‘š, we obtain exact traveling wave solutions of (2.1).

3. Applications of the Method

In this section, we investigate the nonlinear reaction diffusion equation by applying the improved (πΊξ…ž/𝐺)-expansion method for constructing exact traveling wave solutions.

3.1. The Nonlinear Reaction Diffusion Equation

In this work, we consider the nonlinear reaction diffusion equation involving parameters followed by Zayed and Gepreel [57]: 𝑒𝑑𝑑+𝛼𝑒π‘₯π‘₯+𝛽𝑒+𝛾𝑒3=0,(3.1) where 𝛼,𝛽, and 𝛾 are nonzero constants.

Using the traveling wave transformation Equation (2.2), (3.1) is transformed into the ODE: 𝛼+𝑉2ξ€Έπ‘£ξ…žξ…ž+𝛽𝑣+𝛾𝑣3=0,(3.2) where the superscripts indicate the derivatives with respect to πœ‚.

Taking the homogeneous balance between π‘£ξ…žξ…ž and 𝑣3 in (3.2), we obtain π‘š=1. Therefore, the solution of (3.2) is in the form as following:𝑣(πœ‚)=π‘Žβˆ’1ξ€·πΊξ…žξ€Έ/πΊβˆ’1+π‘Ž0+π‘Ž1ξ€·πΊξ…žξ€Έ,/𝐺(3.3) where π‘Žβˆ’1,π‘Ž0, and π‘Ž1 are all constants to be determined.

Substituting (3.3) together with (2.5) into the (3.2), the left-hand side of (3.2) is converted into a polynomial of (πΊξ…ž/𝐺)π‘Ÿ(π‘Ÿ=0,Β±1,Β±2,…). According to Step 4, collecting all terms with the same power of (πΊξ…ž/𝐺) and setting each coefficient of this polynomial to zero yield a set of algebraic equations (which are omitted to display, for simplicity) for π‘Žβˆ’1,π‘Ž0,π‘Ž1,𝑉,πœ†, and πœ‡.

Solving the system of obtained algebraic equations with the aid of algebraic software Maple, we obtain the following.

Case 1. One has π‘Žβˆ’1=0,π‘Ž0ξƒŽ=βˆ“πœ†βˆ’π›½π›Ύξ€·πœ†2ξ€Έβˆ’4πœ‡,π‘Ž1ξƒŽ=Β±2βˆ’π›½π›Ύξ€·πœ†2ξ€ΈξƒŽβˆ’4πœ‡,𝑉=Β±2π›½πœ†2βˆ’4πœ‡βˆ’π›Ό,(3.4) where 𝛼,𝛽,and𝛾 are nonzero constants and πœ†2βˆ’4πœ‡β‰ 0.

Case 2. One has π‘Žβˆ’1ξƒŽ=Β±2πœ‡βˆ’π›½π›Ύξ€·πœ†2ξ€Έβˆ’4πœ‡,π‘Ž0ξƒŽ=βˆ“πœ†βˆ’π›½π›Ύξ€·πœ†2ξ€Έβˆ’4πœ‡,π‘Ž1ξƒŽ=0,𝑉=Β±2π›½πœ†2βˆ’4πœ‡βˆ’π›Ό,(3.5) where 𝛼,𝛽,and𝛾 are nonzero constants and πœ†2βˆ’4πœ‡β‰ 0.
Substituting the general solution Equation (2.5) into (3.3), we obtain two different families of traveling wave solutions of (3.2).

Family 1 (Hyperbolic Function Solutions). When πœ†2βˆ’4πœ‡>0, we obtain𝑣(πœ‚)=π‘Žβˆ’1ξƒ©βˆ’πœ†2+12ξ”πœ†2βˆšβˆ’4πœ‡π΄sinh(1/2)πœ†2βˆšβˆ’4πœ‡πœ‚+𝐡cosh(1/2)πœ†2βˆ’4πœ‡πœ‚βˆšπ΄cosh(1/2)πœ†2βˆšβˆ’4πœ‡πœ‚+𝐡sinh(1/2)πœ†2ξƒͺβˆ’4πœ‡πœ‚βˆ’1+π‘Ž0+π‘Ž1ξƒ©βˆ’πœ†2+12ξ”πœ†2βˆšβˆ’4πœ‡π΄sinh(1/2)πœ†2βˆšβˆ’4πœ‡πœ‚+𝐡cosh(1/2)πœ†2βˆ’4πœ‡πœ‚βˆšπ΄cosh(1/2)πœ†2βˆšβˆ’4πœ‡πœ‚+𝐡sinh(1/2)πœ†2ξƒͺ.βˆ’4πœ‡πœ‚(3.6) Various known solutions can be rediscovered, if 𝐴 and 𝐡 take particular values.
For example:
(i) if 𝐴=0 but 𝐡≠0, we obtain𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2+12ξ”πœ†21βˆ’4πœ‡coth2ξ”πœ†2ξ‚Άβˆ’4πœ‡πœ‚βˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2+12ξ”πœ†21βˆ’4πœ‡coth2ξ”πœ†2ξ‚Ά;βˆ’4πœ‡πœ‚(3.7)(ii) if 𝐡=0 but 𝐴≠0, we obtain𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2+12ξ”πœ†21βˆ’4πœ‡tanh2ξ”πœ†2ξ‚Άβˆ’4πœ‡πœ‚βˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2+12ξ”πœ†21βˆ’4πœ‡tanh2ξ”πœ†2ξ‚Ά;βˆ’4πœ‡πœ‚(3.8)(iii)if 𝐴≠0,𝐴>𝐡, we obtain𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2+12ξ”πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2βˆ’4πœ‡πœ‚+πœ‚0ξ‚Άξ‚Άβˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2+12ξ”πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2βˆ’4πœ‡πœ‚+πœ‚0.ξ‚Άξ‚Ά(3.9)

Family 2 (Trigonometric Function Solutions). When πœ†2βˆ’4πœ‡<0, we obtain𝑣(πœ‚)=π‘Žβˆ’1ξƒ©βˆ’πœ†2+124πœ‡βˆ’πœ†2βˆšβˆ’π΄sin(1/2)4πœ‡βˆ’πœ†2βˆšπœ‚+𝐡cos(1/2)4πœ‡βˆ’πœ†2πœ‚βˆšπ΄cos(1/2)4πœ‡βˆ’πœ†2βˆšπœ‚+𝐡sin(1/2)4πœ‡βˆ’πœ†2πœ‚ξƒͺβˆ’1+π‘Ž0+π‘Ž1ξƒ©βˆ’πœ†2+124πœ‡βˆ’πœ†2βˆšβˆ’π΄sin(1/2)4πœ‡βˆ’πœ†2βˆšπœ‚+𝐡cos(1/2)4πœ‡βˆ’πœ†2πœ‚βˆšπ΄cos(1/2)4πœ‡βˆ’πœ†2βˆšπœ‚+𝐡sin(1/2)4πœ‡βˆ’πœ†2πœ‚ξƒͺ.(3.10) Various known solutions can be rediscovered, if 𝐴 and 𝐡 are taken particular values.
For example,
(iv) if 𝐴=0 but 𝐡≠0, we obtain 𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2+124πœ‡βˆ’πœ†2cot124πœ‡βˆ’πœ†2πœ‚ξ‚Άβˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2+124πœ‡βˆ’πœ†2cot124πœ‡βˆ’πœ†2πœ‚ξ‚Ά;(3.11)(v)if 𝐡=0 but 𝐴≠0, we obtain 𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2βˆ’124πœ‡βˆ’πœ†21tan24πœ‡βˆ’πœ†2πœ‚ξ‚Άβˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2βˆ’124πœ‡βˆ’πœ†21tan24πœ‡βˆ’πœ†2πœ‚ξ‚Ά;(3.12)(vi)if 𝐴≠0,𝐴>𝐡, we obtain 𝑣(πœ‚)=π‘Žβˆ’1ξ‚΅βˆ’πœ†2+124πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚βˆ’πœ‚0ξ‚Άξ‚Άβˆ’1+π‘Ž0+π‘Ž1ξ‚΅βˆ’πœ†2+124πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚βˆ’πœ‚0.ξ‚Άξ‚Ά(3.13)

Family 1 (Hyperbolic Function Solutions). Substituting (3.4) and (3.5) together with the general solution (2.5) into the (3.3), we obtain the hyperbolic function solution Equation (3.6), and then using (3.7), we obtain solutions respectively (if 𝐴=0 but 𝐡≠0),𝑣1ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξ‚΅ξ”βˆ’4πœ‡βˆ’2πœ†+πœ†2ξ‚΅1βˆ’4πœ‡coth2ξ”πœ†2,βˆ’4πœ‡πœ‚ξ‚Άξ‚Ά(3.14) where βˆšπœ‚=π‘₯Β±((2𝛽/(πœ†2βˆ’4πœ‡))βˆ’π›Ό)𝑑,πœ†2βˆ’4πœ‡β‰ 0, and 𝛾≠0. 𝑣2ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξƒ©ξ‚΅βˆ’4πœ‡2πœ‡βˆ’πœ†2+12ξ”πœ†21βˆ’4πœ‡coth2ξ”πœ†2ξ‚Άβˆ’4πœ‡πœ‚βˆ’1ξƒͺ,βˆ’πœ†(3.15) where βˆšπœ‚=π‘₯Β±((2𝛽/(πœ†2βˆ’4πœ‡))βˆ’π›Ό)𝑑,πœ†2βˆ’4πœ‡β‰ 0, and 𝛾≠0.
Again, substituting (3.4) and (3.5) together with the general solution Equation (2.5) into Equation (3.3), we obtain the hyperbolic function solution Equation (3.6), and then using (3.8), our solutions become, respectively (if 𝐡=0 but 𝐴≠0),𝑣3ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξ‚΅ξ”βˆ’4πœ‡βˆ’2πœ†+πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2.βˆ’4πœ‡πœ‚ξ‚Άξ‚Ά(3.16)𝑣4ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξƒ©ξ‚΅βˆ’4πœ‡2πœ‡βˆ’πœ†2+12ξ”πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2βˆ’4πœ‡πœ‚ξ‚Άξ‚Άβˆ’1ξƒͺ.βˆ’πœ†(3.17) Also, substituting (3.4) and (3.5) together with the general solution (2.5) into the (3.3), we obtain the hyperbolic function solution (3.6), and then using (3.9), we obtain the following solutions, respectively (if 𝐴≠0,𝐴>𝐡): 𝑣5ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξ‚΅ξ”βˆ’4πœ‡βˆ’2πœ†+πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2βˆ’4πœ‡πœ‚+πœ‚0,ξ‚Άξ‚Ά(3.18) where βˆšπœ‚=π‘₯Β±((2𝛽/(πœ†2βˆ’4πœ‡))βˆ’π›Ό)𝑑,πœ†2βˆ’4πœ‡β‰ 0,𝛾≠0, and πœ‚0=tanhβˆ’1(𝐡/𝐴). 𝑣6ξƒŽ(πœ‚)=Β±βˆ’π›½π›Ύξ€·πœ†2ξ€Έξƒ©ξ‚΅βˆ’4πœ‡2πœ‡βˆ’πœ†2+12ξ”πœ†2ξ‚΅1βˆ’4πœ‡tanh2ξ”πœ†2βˆ’4πœ‡πœ‚+πœ‚0ξ‚Άξ‚Άβˆ’1ξƒͺ,βˆ’πœ†(3.19) where βˆšπœ‚=π‘₯Β±((2𝛽/(πœ†2βˆ’4πœ‡))βˆ’π›Ό)𝑑,πœ†2βˆ’4πœ‡β‰ 0,𝛾≠0, and πœ‚0=tanhβˆ’1(𝐡/𝐴).

Family 2 (Trigonometric Function Solutions). Substituting (3.4) and (3.5) together with the general solution Equation (2.5) into the (3.3), we obtain the trigonometric function solution Equation (3.10), and then using (3.11), our solutions become respectively (if 𝐴=0 but 𝐡≠0),𝑣7ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†2ξ€Έξ‚΅ξ”βˆ’2πœ†+4πœ‡βˆ’πœ†2cotξ‚΅124πœ‡βˆ’πœ†2πœ‚,ξ‚Άξ‚Ά(3.20) where 4πœ‡βˆ’πœ†2β‰ 0 and𝛾≠0. 𝑣8ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†22πœ‡βˆ’πœ†2+124πœ‡βˆ’πœ†2cotξ‚΅124πœ‡βˆ’πœ†2πœ‚ξ‚Άξ‚Άβˆ’1ξƒͺ,βˆ’πœ†(3.21) where 4πœ‡βˆ’πœ†2β‰ 0 and𝛾≠0.
Also, substituting (3.4) and (3.5) together with the general solution Equation (2.5) into the (3.3), we obtain the trigonometric function solution Equation (3.10), and then using (3.12), our traveling wave solutions become respectively (if 𝐡=0 but 𝐴≠0),𝑣9ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†2ξ€Έξ‚΅ξ”βˆ’2πœ†βˆ’4πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚.ξ‚Άξ‚Ά(3.22)𝑣10ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†22πœ‡βˆ’πœ†2βˆ’124πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚ξ‚Άξ‚Άβˆ’1ξƒͺ.βˆ’πœ†(3.23) Moreover, substituting (3.4) and (3.5) together with the general solution Equation (2.5) into Equation (3.3), we obtain the trigonometric function solution Equation (3.10), and then using (3.13), our obtained solutions (if 𝐴≠0,𝐴>𝐡) are as follows: 𝑣11ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†2ξ€Έξ‚΅ξ”βˆ’2πœ†+4πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚βˆ’πœ‚0,ξ‚Άξ‚Ά(3.24) where πœ‚0=tanβˆ’1(𝐡/𝐴). 𝑣12ξƒŽ(πœ‚)=±𝛽𝛾4πœ‡βˆ’πœ†22πœ‡βˆ’πœ†2+124πœ‡βˆ’πœ†2ξ‚΅1tan24πœ‡βˆ’πœ†2πœ‚βˆ’πœ‚0ξ‚Άξ‚Άβˆ’1ξƒͺ,βˆ’πœ†(3.25) where πœ‚0=tanβˆ’1(𝐡/𝐴).

4. Results and Discussion

It is noteworthy to mention that some of our obtained solutions are in good agreement with the existing results which are shown in Table 1. Furthermore, the graphical presentations of some of obtained solutions are depicted in the Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11.

4.1. Comparison between Zayed and Gepreel [57] Solutions and Our Solutions

Beyond Table 1, we obtain many new exact traveling wave solutions which have not been found in the previous literature.

4.2. Graphical Representations of the Solutions

The graphical illustrations of the solutions are described in the Figures with the aid of Maple.

5. Conclusions

In this paper, we obtain abundant new exact traveling wave solutions for the nonlinear reaction diffusion equation involving parameters by applying the improved (πΊξ…ž/𝐺)-expansion method. The obtained solutions are expressed in terms of the hyperbolic and the trigonometric function forms. The solutions of the nonlinear reaction diffusion equation have many potential applications in biological sciences. The validity of the obtained traveling wave solutions is proved by comparing with the published results. We expect that the used method will be effectively used to construct many new exact traveling wave solutions for other kinds of nonlinear evolution equations which are arising in technical arena.

Acknowledgments

This paper is supported by the USM short-term grant (reference no. 304/PMATHS/6310072) and the authors would like to express their thanks to the School of Mathematical Sciences, USM, for providing related research facilities. The authors are also grateful to the referee(s) for their valuable comments and suggestions.