AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 963852 10.1155/2014/963852 963852 Research Article Traveling Wave Solutions and Infinite-Dimensional Linear Spaces of Multiwave Solutions to Jimbo-Miwa Equation Zhang Lijun 1, 2 Khalique C. M. 2 Rui Weiguo 1 Department of Mathematics, School of Science Zhejiang Sci-Tech University Hangzhou, Zhejiang 310018 China zstu.edu.cn 2 Department of Mathematical Sciences International Institute for Symmetry Analysis and Mathematical Modelling North-West University Mafikeng Campus, Private Bag X 2046, Mmabatho 2735 South Africa nwu.ac.za 2014 942014 2014 23 01 2014 19 03 2014 9 4 2014 2014 Copyright © 2014 Lijun Zhang and C. M. Khalique. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The traveling wave solutions and multiwave solutions to (3 + 1)-dimensional Jimbo-Miwa equation are investigated in this paper. As a result, besides the exact bounded solitary wave solutions, we obtain the existence of two families of bounded periodic traveling wave solutions and their implicit formulas by analysis of phase portrait of the corresponding traveling wave system. We derive the exact 2-wave solutions and two families of arbitrary finite N-wave solutions by studying the linear space of its Hirota bilinear equation, which confirms that the (3 + 1)-dimensional Jimbo-Miwa equation admits multiwave solutions of any order and is completely integrable.

1. Introduction

Various nonlinear partial differential equations (NLPDEs) have been proposed to model different kinds of phenomena in natural and applied sciences such as fluid dynamics, plasma physics, solid-state physics, optical fibers, acoustics, mechanics, biology, and mathematical finance. Obviously, it is of significant importance to study the solutions of such NLPDEs from both theoretical and practical points of view. However, the solution spaces of nonlinear equations are infinite-dimensional and contain diverse solution structures, so it is usually a difficult job to determine the solutions to nonlinear NLPDEs.

A great idea to generate exact solutions of NLPDEs is to reduce the NLPDEs into some algebraic equations by assuming the solutions to have some special forms or satisfy some solvable simpler equations. This can be seen in, for example, the exp-function method , the tanh function method , the homogeneous balance method [3, 4], the auxiliary function method [5, 6], the sech-function method , the sine-cosine method [8, 9], the tanh-coth method , the Jacobi elliptic function method , the F -expansion method, and the extended F -expansion method . Normally, it is not an easy task to solve these nonlinear algebraic equations reduced from NLPDEs because they involve many parameters. However, the rapid development of symbolic computation makes it relatively easy to solve these algebraic equations . At the same time one should be very careful when applying these methods as different methods might give the same solutions. For example, solutions obtained by using the sech-function method, the tanh-coth method, and the exp-function method are actually the same [14, 15], because f ( sech 2 x ) = f ( 1 - tanh 2 x ) = f ( 4 / ( e x + e - x ) 2 ) , for any function f .

Recently, the planar dynamical system theorem has been employed to study the traveling wave solutions of NLPDEs [16, 17]. The best advantage of this approach is that the boundedness, periodicity, the shapes of the solutions, and even the singular traveling wave solutions of NLPDEs can be recognized easily from the corresponding orbits of their phase portraits under various different parametric conditions. Also the exact solutions can be derived at the same time. It is worth pointing out that Hirota’s bilinear method  is also an amazing method to find exact solutions to some NLPDEs. By some independent variable transformation, various nonlinear equations of mathematical physics are transformed into Hirota’s bilinear equations [18, 19], which possess some specific properties and thus might be applied to study the solution sets of nonlinear differential equations. Recently, even some programs have been designed and some algorithms have been proposed on searching for integrable bilinear equations . Based on Hirota’s bilinear form, soliton solutions were obtained by the Hirota perturbation technique , the multiple exp-function algorithm , and other methods . Even some generalized bilinear form has been proposed recently by Ma .

The ( 3 + 1 ) -dimensional Jimbo-Miwa equation (1) u x x x y + 3 ( u x u y ) x + 2 u y t - 3 u x z = 0 was first proposed and studied by Jimbo and Miwa . Recently, this equation has attracted a great deal of attention which mainly focuses on its solutions, integrability properties and symmetries. Ma and Lee  proposed a direct approach to solve (1) by using rational function transformations. Li and Dai  applied the generalized Riccati equation method to look for its exact solutions. However, in , Kudryashov and Sinelshchikov pointed out some mistakes and proved that some solutions in  could be written in a uniform form and thus they are not new at all.

The ( 3 + 1 ) - dimensional Jimbo-Miwa equation can be transformed into the Hirota bilinear equation (2) ( D x 3 D y + 2 D t D y - 3 D x D z ) f · f = 0 through the dependent variable transformation u = 2 ( ln f ) x .

In this paper, we firstly study the bounded traveling wave solutions of the Jimbo-Miwa equation (1) by investigating the bifurcation and phase portraits of a planar cubic polynomial ordinary differential equation by using the planar dynamical system theory [16, 37]. We then explore the 1-wave and 2-wave solutions and the infinite-dimensional linear spaces of multiwave solutions to Jimbo-Miwa equation by employing Hirota’s bilinear method, thus confirming that Jimbo-Miwa equation possesses multiwave solutions of arbitrary order.

2. Traveling Wave Solutions to the <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10"> <mml:mo mathvariant="bold">(</mml:mo> <mml:mn>3</mml:mn> <mml:mo mathvariant="bold"> </mml:mo> <mml:mo mathvariant="bold">+</mml:mo> <mml:mo mathvariant="bold"> </mml:mo> <mml:mn>1</mml:mn> <mml:mo mathvariant="bold">)</mml:mo> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Dimensional Jimbo-Miwa Equation

To investigate the traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation, we make the traveling wave transformation u ( x , y , z , t ) = u ( ξ ) , ξ = k x + l y + m z + n t , under which (1) is reduced to the nonlinear ordinary differential equation (3) k 3 l d 4 u d ξ 4 + 6 k 2 l d u d ξ d 2 u d ξ 2 + ( 2 l n - 3 k m ) d 2 u d ξ 2 = 0 . Integrating (1) once with respect to ξ gives (4) k 3 l d 3 u d ξ 3 + 3 k 2 l ( d u d ξ ) 2 + ( 2 l n - 3 k m ) d u d ξ = g , where g is an arbitrary constant. Let v = d u / d ξ ; then (1) becomes (5) k 3 l d 2 v d ξ 2 + 3 k 2 l v 2 + ( 2 l n - 3 k m ) v = g , which is a second-order nonlinear ordinary differential equation. We will study the solutions of (5) by planar dynamical system method and thus derive the traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation.

2.1. Bounded Solutions of (<xref ref-type="disp-formula" rid="EEq5">5</xref>)

First, we rewrite (5) in a simpler and more general form, namely, (6) d 2 v d ξ 2 = a v 2 + b v + c , where a = - 3 / k , b = ( 3 k m - 2 l n ) / ( k 3 l ) , and c = g / ( k 3 l ) . We now study the bifurcation and exact solutions of (6). Let v = w ; then (6) is equivalent to the dynamical system (7) v = w , w = a v 2 + b v + c , which has the Hamiltonian (8) H ( v , w ) = w 2 2 - a 3 v 3 - b 2 v 2 - c v . Clearly, the phase orbits defined by the vector fields of system (7) determine all solutions of (6). The bounded solutions of (6) correspond to the bounded phase orbits of system (7), which we now investigate. Along the orbit corresponding to H ( v , w ) = h , (9) ( d v d ξ ) 2 = 2 a 3 v 3 + b v 2 + 2 c v + 2 h . Consequently, the general formula of the solutions of (6) can be expressed as (10) 0 v d v ( 2 a v 3 / 3 ) + b v 2 + 2 c v + 2 h = ± ξ 0 ξ d ξ .

However, it is not easy to know the properties and the shapes of (10) which actually is determined by the parameters a , b , c , and h . Clearly, the abscissas of equilibrium points of system (7) are the zeros of a v 2 + b v + c = 0 . Obviously, the system has no bounded orbits when b 2 - 4 a c < 0 . We suppose that b 2 - 4 a c > 0 in order to study the bounded orbits of system (7). Denoting v ± = ( - b ± b 2 - 4 a c ) / 2 a , then E + ( v + , 0 ) and E - ( v - , 0 ) are two equilibrium points of system (7). By the theory of planar dynamical system, we know that E + is a saddle point and E - is a center. Denote h ± = H ( v ± , 0 ) , and, by careful computation, we get (11) h ± = ( b 2 - 4 a c ) ( - b ± b 2 - 4 a c ) + 2 a b c 12 a 2 . Obviously, h - < h < h + . H ( v , w ) = h + corresponds to homoclinic orbits, and H ( v , w ) = h - corresponds to the center E - and H ( v , w ) = h , where h + < h < h - corresponds to a family of closed orbits surrounding the center E - , which are surrounded by a homoclinic orbit. That is to say, (10) defines bounded solutions if and only if h + h < h - . To be exact, (10) defines a family of periodic solutions when h + < h < h - .

When h = h + , (10) defines a bounded solution which approaches v + as ξ goes to infinity. Actually, (12) 2 a 3 v 3 + b v 2 + 2 c v + 2 h + = 2 a 3 ( v - v + ) 2 ( v - v 0 ) , where v 0 = - ( b + 2 b 2 - 4 a c ) / 2 a , so (10) can be reduced to (13) v 0 v d v ( v - v + ) a ( v - v 0 ) = 2 3 ( ξ - ξ 0 ) , from which we can get the exact solution as (14) v = v + - ( v + - v 0 ) sech 2 ( a ( v + - v 0 ) 6 ( ξ - ξ 0 ) ) .

By further simplification, (14) becomes (15) v = v + - 3 b 2 - 4 a c 2 a sech 2 [ 1 2 ( b 2 - 4 a c ) 1 / 4 ( ξ - ξ 0 ) ] , which is an exact bounded solution of (5).

Thus, we have the following lemma.

Lemma 1.

The general second-order ODE (6) has bounded solutions if and only if b 2 - 4 a c > 0 . The bounded solutions can be expressed as (10) in an implicit form. In fact, provided h - < h < h + , (10) defines a family of bounded periodic solutions and h = h + defines a bounded solution which approaches v + as ξ goes to infinity and can be expressed explicitly as (15), where v + = ( - b + b 2 - 4 a c ) / 2 a and h ± is defined by (11).

2.2. Bounded Traveling Wave Solutions to the <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M69"> <mml:mo mathvariant="bold">(</mml:mo> <mml:mn>3</mml:mn> <mml:mo mathvariant="bold"> </mml:mo> <mml:mo mathvariant="bold">+</mml:mo> <mml:mo mathvariant="bold"> </mml:mo> <mml:mn>1</mml:mn> <mml:mo mathvariant="bold">)</mml:mo> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Dimensional Jimbo-Miwa Equation

According to the analysis and results in Section 2.1, we know that (5) has only two kinds of bounded solutions, among which one is a family of periodic solutions and another is a family of solutions approaching a fixed number as ξ goes to infinity. Note that what we aim to study is the bounded traveling wave solutions to the ( 3 + 1 ) -dimensional Jimbo-Miwa equation which are determined by v = d u / d ξ and v satisfies (5). So we have to investigate how we can get the bounded solutions to (3) from the bounded solution of (5).

Clearly, u ( ξ ) = ξ 0 ξ v ( ξ ) d ξ , and v ( ξ ) can be expressed implicitly as (10). By the geometry meaning of the integral and the properties of the solutions of (5), we get the traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation. To get the bounded solution, we choose the integral constant g to be 0; that is, c = 0 in (15), and so (16) u ( ξ ) = g - 3 | b | a tanh [ | b | 2 ( ξ - ξ 0 ) ] ; that is, (17) u ( ξ ) = g - k | 2 l n - 3 k m k 3 l | × tanh [ 1 2 | 2 l n - 3 k m k 3 l | ( ξ - ξ 0 ) ] , which is a family of exact bounded kink traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation, where g is an arbitrary constant.

However, we may not get bounded solutions from the family of periodic solutions of (5). It is easy to see that if v ( ξ ) is a periodic solution of (5), then u ( ξ ) = ξ 0 ξ v ( ξ ) d ξ is bounded if and only if 0 T v ( ξ ) d ξ = 0 , where T is the period of the function v ( ξ ) . Recall that the period of the function v ( ξ ) which is defined by (10) with h - < h < h + continuously depends on the parameters a , b , c , and h . So 0 T v ( ξ ) d ξ continuously depends on the parameters a , b , c , and h too. Let V ( a , b , c , h ) = 0 T v ( ξ ) d ξ ; then V ( a , b , c , h ) is a continuous function of a , b , c , and h . We now prove the existence of the root of V ( a , b , c , h ) = 0 to get the existence of the bounded periodic traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation.

By the theory of planar dynamical system, when h - < h < h + , (18) 2 a 3 v 3 + b v 2 + 2 c v + 2 h + = 2 a 3 ( v - v 1 ) ( v - v 2 ) ( v - v 3 ) , where v 1 > v + > v 2 > v - > v 3 > v 0 in the case of a > 0 and v 0 > v 1 > v - > v 2 > v + > v 3 in the case of a < 0 . The periodic solution satisfies v + > v 2 > v ( ξ ) > v 3 > v 0 when a > 0 and v 0 > v 1 > v ( ξ ) > v 2 > v + when a < 0 . Clearly, v 0 > 0 when a > 0 and 0 > b > - 4 a c / 3 . So v ( ξ ) > 0 and thus V ( a , b , c , h ) > 0 when a > 0 and 0 > b > - 4 a c / 3 . However, v ( ξ ) < v + < 0 and thus V ( a , b , c , h ) < 0 when a > 0 , c > 0 , and b > 0 . There must exist at least one zero of V ( a , b , c , h ) = 0 in the region a > 0 of the parameter space since V ( a , b , c , h ) is a continuous function of a , b , c , and h . The same happens in the region a < 0 of the parameter space. So, we know that there exist at least two families of bounded periodic traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation. Thus we have the following theorem.

Theorem 2.

The ( 3 + 1 ) - dimensional Jimbo-Miwa equation has two types of bounded traveling wave solutions as given below.

The ( 3 + 1 ) - dimensional Jimbo-Miwa equation has a family of exact bounded kink traveling wave solutions (19) u ( x , y , z , t ) = g - k | 2 l n - 3 k m k 3 l | × tanh [ 1 2 | 2 l n - 3 k m k 3 l | ( k x + l y + m z + n t - ξ 0 ) ] ,

where ξ 0 and g are two arbitrary constants.

The ( 3 + 1 ) -dimensional Jimbo-Miwa equation has at least two families of bounded periodic traveling wave solutions which are determined implicitly by (10) and u ( ξ ) = ξ 0 ξ v ( ξ ) d ξ , where ξ = k x + l y + m z + n t and ξ 0 is an arbitrary constant.

3. <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M146"> <mml:mi>N</mml:mi> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Wave Solutions Linear Subspace of the <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M147"> <mml:mo mathvariant="bold">(</mml:mo> <mml:mn>3</mml:mn> <mml:mo mathvariant="bold"> </mml:mo> <mml:mo mathvariant="bold">+</mml:mo> <mml:mo mathvariant="bold"> </mml:mo> <mml:mn>1</mml:mn> <mml:mo mathvariant="bold">)</mml:mo> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Dimensional Jimbo-Miwa Equation

In this section we study the N - wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation by the linear superposition principle to Hirota bilinear equations which was proposed firstly by Ma and Fan .

Let wave variables η i = k i · x , 1 i N , where k i = { k 1 , i , k 2 , i , k 3 , i , k 4 , i } T , x = { x , y , z , t } T . According to the linear superposition principle to Hirota bilinear equations [30, 31], we can get the following sufficient and necessary conditions for span { e η 1 , , e η N } being a subspace of the Hirota bilinear equation (2).

Theorem 3.

Let k i = { k 1 , i , k 2 , i , k 3 , i , k 4 , i } T , x = { x , y , z , t } T , and f N = i = 1 N ϵ i e k i · x . Then, for any constants { ϵ i , 1 i N } , f N solves the Hirota bilinear equation (2) if and only if (20) ( k 1 , i - k 1 , j ) 3 ( k 2 , i - k 2 , j ) + 2 ( k 4 , i - k 4 , j ) ( k 2 , i - k 2 , j ) - 3 ( k 1 , i - k 1 , j ) ( k 3 , i - k 3 , j ) = 0 for any 1 j < i N .

It follows from Theorem 3 that f N is an N -wave solution to the Hirota bilinear equation (2) if (20) holds for any 1 j < i N . However, the corresponding u ( x , y , z , t ) = 2 ( f N ) x / f N is not necessarily an N -wave solution to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1) even if { η 1 , η 2 , , η N } is linear independent. In fact, (21) u ( x , y , z , t ) = 2 ( f N ) x f N = 2 i = 1 N ϵ i k 1 , i e k i · x i = 1 N ϵ i e k i · x = 2 k 1,1 + i = 2 N k 1 , i ϵ i e ( k i - k 1 ) · x 1 + i = 2 N ϵ i e ( k i - k 1 ) · x = 2 k 1,1 + 2 i = 2 N k 1 , i ϵ i e k i · x 1 + i = 2 N ϵ i e k i · x , where ϵ i = ϵ i / ϵ 1 , k 1 , i = k 1 , i - k 1,1 and k i = k i - k 1 . Obviously, the function above is an N - 1 -wave solution to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1) if { η i ,   i = 2,3 , , N } is independent; that is, { k i ,   i = 2 , , N } is linear independent, solves (20), and satisfies the dispersion relations. Here, k i = { k 1 , i , k 2 , i , k 3 , i , k 4 , i } T satisfies the dispersion relations meaning that k 1 , i 3 k 2 , i + 2 k 4 , i k 2 , i - 3 k 1 , i k 3 , i = 0 holds.

Recall that the goal of this paper is to investigate the multiwave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation. Combining the conclusion of Theorem 3 with the analysis above, we can now state the following theorem on the N -wave solution to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1).

Theorem 4.

Let k i = { k 1 , i , k 2 , i , k 3 , i , k 4 , i } T , x = { x , y , z , t } T , 1 i N . For any constants g and { ϵ i , 1 i N } , suppose that (22) u ( x , y , z , t ) = 2 i = 1 N k 1 , i ϵ i e k i · x 1 + i = 1 N ϵ i e k i · x + g . Then u ( x , y , z , t ) is an N -wave solution to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1), if { k 1 , k 2 , , k N } is a linear independent solution set of equation (20) satisfying the dispersion relations.

To get other N -wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation, we study the linear independent solution sets to (20) which are required to satisfy the dispersion relations. Actually, (20) is a system possessing N ( N + 1 ) / 2 coupled equations (plus the dispersion conditions). It is usually not so easy to get the solutions of (20). Fortunately, the number of the equations is 1 when N = 1 and it is 3 when N = 2 , which might make it easy to get the solutions in these two cases. Clearly, 1-wave solution is the traveling wave solution. Let us check what kinds of traveling wave solutions we can obtain by the Hirota bilinear method first.

For the case N = 1 in (22), we only need to find the independent solutions k 1 of the dispersion relation (23) ( k 1,1 ) 3 k 2,1 + 2 k 4,1 k 2,1 - 3 k 1,1 k 3,1 = 0 .

That is to say, we get the traveling wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation from (20). Suppose that l i , i = 1,2 , 3,4 are any arbitrary constants, l 1 0 and l 3 = l 2 ( l 1 3 + 2 l 4 ) / ( 3 l 1 ) ; then k 1 = { l 1 , l 2 , l 3 , l 4 } T satisfies (23), and thus (24) u ( x , y , z , t ) = 2 l 1 ϵ e k i · x 1 + ϵ e k i · x + g , where ϵ and g are any arbitrary constants, is a 1 -wave solution, that is, traveling wave solution to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation.

By further computation, (24) could be rewritten as (25) u ( x , y , z , t ) = g - 2 l 1 1 + ϵ e k i · x . For ϵ > 0 , (26) 2 l 1 1 + ϵ e k i · x = l 1 [ 1 - tanh 1 2 ( k i · x + ξ 0 ) ] and, for ϵ < 0 , (27) 2 l 1 1 + ϵ e k i · x = l 1 [ 1 - coth 1 2 ( k i · x + ξ 0 ) ] , where ξ 0 = ln | ϵ | is an arbitrary constant. Consequently, (24) can be rewritten as (28) u ( x , y , z , t ) = g - l 1 tanh 1 2 ( k i · x + ξ 0 ) , (29) u ( x , y , z , t ) = g - l 1 coth 1 2 ( k i · x + ξ 0 ) . By parametric transformation, we can find that the traveling wave solution (28) is exactly the same as solution (19) in Section 2. However, (29) is a family of unbounded solutions.

3.1. 2-Wave Solutions to the <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M224"> <mml:mo mathvariant="bold">(</mml:mo> <mml:mn>3</mml:mn> <mml:mo mathvariant="bold"> </mml:mo> <mml:mo mathvariant="bold">+</mml:mo> <mml:mo mathvariant="bold"> </mml:mo> <mml:mn>1</mml:mn> <mml:mo mathvariant="bold">)</mml:mo> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Dimensional Jimbo-Miwa Equation

To obtain the 2-wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation, we now study the solutions to the following coupled algebraic equations: (30) k 1,1 3 k 2,1 + 2 k 4,1 k 2,1 - 3 k 1,1 k 3,1 = 0 , k 1,2 3 k 2,2 + 2 k 4,2 k 2,2 - 3 k 1,2 k 3,2 = 0 , ( k 1,2 - k 1,1 ) 3 ( k 2,2 - k 2,1 ) + 2 ( k 4,2 - k 4,1 ) ( k 2,2 - k 2,1 ) - 3 ( k 1,2 - k 1,1 ) ( k 3,2 - k 3,1 ) = 0 .

From (30), we get (31) k 3,1 = k 1,1 3 + 2 k 4,1 3 k 1,1 k 2,1 , (32) k 3,2 = k 1,2 3 + 2 k 4,2 3 k 1,2 k 2,2 , (33) k 3,2 - k 3,1 = ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) 3 ( k 1,2 - k 1,1 ) × ( k 2,2 - k 2,1 ) . Clearly, from (31)–(33), we have (34) ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) 3 ( k 1,2 - k 1,1 ) ( k 2,2 - k 2,1 ) = k 1,2 3 + 2 k 4,2 3 k 1,2 k 2,2 - k 1,1 3 + 2 k 4,1 3 k 1,1 k 2,1 , from which we get (35) k 2,2 = ( - k 1,2 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) - 1 ( ( k 1,1 3 + 2 k 4,1 ) ( k 1,2 - k 1,1 ) - k 1,1 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) × ( ( k 1,2 3 + 2 k 4,2 ) ( k 1,2 - k 1,1 ) - k 1,2 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) - 1 ) × k 1,1 k 1,2 k 2,1 . Now, substituting (35) into (32), we obtain (36) k 3,2 = ( - k 1,2 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) - 1 ( ( k 1,1 3 + 2 k 4,1 ) ( k 1,2 - k 1,1 ) - k 1,1 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) × ( ( k 1,2 3 + 2 k 4,2 ) ( k 1,2 - k 1,1 ) - k 1,2 [ ( k 1,2 - k 1,1 ) 3 + 2 ( k 4,2 - k 4,1 ) ] ) - 1 ) × ( k 1,2 3 + 2 k 4,2 ) 3 k 1,2 2 k 1,1 k 2,1 . Consequently, we get the 2-wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation and we can now state the following theorem.

Theorem 5.

Let k i = { k 1 , i , k 2 , i , k 3 , i , k 4 , i } T , x = { x , y , z , t } T , i = 1,2 . For any constants g and { ϵ i , i = 1,2 } , (37) u ( x , y , z , t ) = 2 k 1,1 ϵ 1 e k 1 · x + k 1,2 ϵ 2 e k 2 · x 1 + ϵ 1 e k 1 · x + ϵ 2 e k 2 · x + g is a 2 -wave solution to the ( 3 + 1 ) -dimensional Jimbo-Miwa equation (1), where k 1,1 , k 1,2 , k 2,1 , k 4,1 , and k 4,2 are arbitrary constants and the constants k 3,1 , k 2,2 , and k 3,2 are determined by (31), (35), and (36), respectively.

3.2. <bold /> <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M248"> <mml:mrow> <mml:mi>N</mml:mi></mml:mrow> </mml:math></inline-formula>-Wave Solutions to the <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M249"> <mml:mo mathvariant="bold">(</mml:mo> <mml:mn>3</mml:mn> <mml:mo mathvariant="bold"> </mml:mo> <mml:mo mathvariant="bold">+</mml:mo> <mml:mo mathvariant="bold"> </mml:mo> <mml:mn>1</mml:mn> <mml:mo mathvariant="bold">)</mml:mo> <mml:mtext>-</mml:mtext></mml:math> </inline-formula>Dimensional Jimbo-Miwa Equation

To get the N -wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation, we need to investigate the independent solution set { k 1 , k 2 , , k N } of (20), which satisfies the dispersion relations. Generally, it is very difficult to get the solution if N is greater than 3 in which case the number of (20) is N ( N + 1 ) / 2 . However, it might be possible to solve these equations by assuming some special relations between these parameters [30, 31].

Following the idea in [30, 31], suppose that k i = { a 1 k i n 1 , a 2 k i n 2 , a 3 k i n 3 , a 4 k i n 4 } T ,    i = 1,2 , , N , where a j s and n j s , j = 1,2 , 3,4 , are parameters to be determined later. To get the possible solutions to (20), we choose ( n 1 , n 2 , n 3 , n 4 ) to be ( 1,2 , 4,3 ) or ( 1 , - 1,1 , 3 ) .

For the case when ( n 1 , n 2 , n 3 , n 4 ) = ( 1,2 , 4,3 ) , substituting k i = { a 1 k i n 1 , a 2 k i n 2 , a 3 k i n 3 , a 4 k i n 4 } T , i = 1,2 , , N , into (20) gives (38) a 1 3 a 2 ( k i - k j ) 3 ( k i 2 - k j 2 ) + 2 a 4 a 2 ( k i 3 - k j 3 ) ( k i 2 - k j 2 ) - 3 a 1 a 4 ( k i - k j ) ( k i 4 - k j 4 ) = 0 . Obviously, (38) holds for arbitrary values of k i and k j if and only if a 1 , a 2 , a 3 , and a 4 satisfy the two equations (39) 2 a 2 ( a 1 3 + 2 a 4 ) - 3 a 1 a 3 = 0 , a 2 ( a 1 3 - 4 a 4 ) = 0 . Solving (39), we get a 4 = ( 1 / 4 ) a 1 3 , a 3 = a 1 2 a 2 , where a 1 and a 2 are free parameters. It is easy to see that k i = { a 1 k i , a 2 k i 2 , a 3 k i 4 , a 4 k i 3 } T , i = 1,2 , , N , satisfy the dispersion condition. Consequently, we know that the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1) possesses N -wave solution for any arbitrary positive integer N . Thus, we have the following theorem.

Theorem 6.

For any arbitrary positive integer N and for any different k i , i = 1,2 , , N , (40) u ( x , y , z , t ) = 2 a 1 i = 1 N k i ϵ i e a 1 k i x + a 2 k i 2 y + a 1 2 a 2 k i 4 z + ( 1 / 4 ) a 1 3 k i 3 t 1 + i = 1 N ϵ i e a 1 k i x + a 2 k i 2 y + a 1 2 a 2 k i 4 z + ( 1 / 4 ) a 1 3 k i 3 t + g is a family of N -wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1), where ϵ i , i = 1,2 , , N , and a 1 and a 2 are arbitrary constants.

For the case when ( n 1 , n 2 , n 3 , n 4 ) = ( 1 , - 1,1 , 3 ) , substituting k i = { a 1 k i , a 2 k i - 1 , a 3 k i , a 4 k i 3 } T , i = 1,2 , , N , into (20) gives (41) a 1 3 a 2 ( k i - k j ) 3 ( k i - 1 - k j - 1 ) + 2 a 4 a 2 ( k i 3 - k j 3 ) ( k i - 1 - k j - 1 ) - 3 a 1 a 4 ( k i - k j ) ( k i 3 - k j 3 ) = 0 . Clearly, (41) holds for arbitrary values of k i and k j if and only if a 1 , a 2 , a 3 , and a 4 satisfy (42) 2 a 2 ( a 1 3 - a 4 ) - 3 a 1 a 3 = 0 , a 2 ( a 1 3 + 2 a 4 ) = 0 . Solving (42), we obtain a 4 = - a 1 3 / 2 , a 3 = a 1 2 a 2 , where a 1 and a 2 are free parameters. Unfortunately, unlike the case when ( n 1 , n 2 , n 3 , n 4 ) = ( 1,2 , 4,3 ) above, it is easy to see that k i = { a 1 k i , a 2 k i - 1 , a 3 k i , a 4 k i 3 } T , i = 1,2 , , N , do not satisfy the dispersion condition any more. However, if we let k 0 = { a 1 , a 2 , a 3 , a 4 } T , then, for any i = 1,2 , , N , (43) k i = k i - k 0 = { a 1 ( k i - 1 ) , a 2 ( k i - 1 - 1 ) , a 3 ( k i - 1 ) , a 4 ( k i 3 - 1 ) } T satisfies the dispersion relations because { k 0 , k 1 , k 2 , , k N } is a solution set of (20).

Consequently, besides the family of N -wave solution (40), the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1) has another group of N -wave solution for any arbitrary positive integer N and so we have the following theorem.

Theorem 7.

For any arbitrary positive integer N and for any different k i , i = 1,2 , , N , (44) u ( x , y , z , t ) = 2 a 1 i = 1 N ( k i - 1 ) ϵ i e a 1 ( k i - 1 ) x + a 2 ( k i - 1 - 1 ) y + a 1 2 a 2 ( k i - 1 ) z - ( 1 / 2 ) a 1 3 ( k i 3 - 1 ) t 1 + i = 1 N ϵ i e a 1 ( k i - 1 ) x + a 2 ( k i - 1 - 1 ) y + a 1 2 a 2 ( k i - 1 ) z - ( 1 / 2 ) a 1 3 ( k i 3 - 1 ) t + g is a family of N -wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1), where ϵ i , i = 1,2 , , N , and a 1 and a 2 are arbitrary constants.

4. Concluding Remarks

The dynamical system theory was employed to study the traveling wave solutions of the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1). The multiwave solutions were investigated by studying the linear space of the corresponding Hirota bilinear equation. The exact formulas of two families of multiwave solutions of any order were obtained as well. Wazwaz  employed Hirota’s bilinear method and derived 1-wave and 2-wave solutions to this equation and stated that it is completely integrable and it admits multiple-soliton solutions of any order. In this paper we explicitly obtained two different families of N -wave solutions to the ( 3 + 1 ) - dimensional Jimbo-Miwa equation (1) given by (40) and (44) and this confirmed the statement given in  that the Jimbo-Miwa equation (1) admits multiple-wave solutions of any order N .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the Nature Science Foundation of China (no. 11101371). This paper is motivated by nice discussion with Professor Wen-Xiu Ma during the authors’ visit at University of South Florida in March 2013. The authors would like to express their sincere appreciations to Professor Ma and the anonymous reviewers.

He J.-H. Wu X.-H. Exp-function method for nonlinear wave equations Chaos, Solitons & Fractals 2006 30 3 700 708 10.1016/j.chaos.2006.03.020 MR2238695 ZBL1141.35448 Parkes E. J. Duffy B. R. An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations Computer Physics Communications 1996 98 3 288 300 2-s2.0-0030284901 10.1016/0010-4655(96)00104-X Wang M. Zhou Y. Li Z. Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics Physics Letters A 1996 216 1–5 67 75 2-s2.0-30244524644 Feng Z. Comment on ‘on the extended applications of homogeneous balance method’ Applied Mathematics and Computation 2004 158 2 593 596 10.1016/j.amc.2003.10.003 MR2094643 Zhang S. Xia T. A generalized new auxiliary equation method and its applications to nonlinear partial differential equations Physics Letters A 2007 363 5-6 356 360 10.1016/j.physleta.2006.11.035 MR2287814 ZBL1197.35008 Zhang H. New exact travelling wave solutions of nonlinear evolution equation using a sub-equation Chaos, Solitons & Fractals 2009 39 2 873 881 10.1016/j.chaos.2007.01.132 MR2518927 ZBL1197.35257 Ma W. X. Travelling wave solutions to a seventh order generalized KdV equation Physics Letters A 1993 180 3 221 224 10.1016/0375-9601(93)90699-Z MR1238444 Yan C. A simple transformation for nonlinear waves Physics Letters A 1996 224 1-2 77 84 10.1016/S0375-9601(96)00770-0 MR1427895 ZBL1037.35504 Wazwaz A.-M. A sine-cosine method for handling nonlinear wave equations Mathematical and Computer Modelling 2004 40 5-6 499 508 10.1016/j.mcm.2003.12.010 MR2102204 ZBL1112.35352 Lan H. B. Wang K. L. Exact solutions for two nonlinear equations. I Journal of Physics A. Mathematical and General 1990 23 17 3923 3928 MR1077305 10.1088/0305-4470/23/17/021 ZBL0718.35020 Parkes E. J. Duffy B. R. Abbott P. C. The Jacobi elliptic-function method for finding periodic-wave solutions to nonlinear evolution equations Physics Letters A 2002 295 5-6 280 286 10.1016/S0375-9601(02)00180-9 MR1925472 ZBL1052.35143 Cai G. L. Wang Q. C. Huang J. J. A modified F-expansion method for solving breaking soliton equation International Journal of Nonlinear Science 2006 2 122 128 Hereman W. Zhuang W. Symbolic computation of solitons with Macsyma Computational and Applied Mathematics, II (Dublin, 1991) 1992 Amsterdam, The Netherlands North-Holland 287 296 MR1204693 ZBL0765.35048 Kudryashov N. A. Seven common errors in finding exact solutions of nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2009 14 9-10 3507 3529 10.1016/j.cnsns.2009.01.023 MR2509943 ZBL1221.35342 Zhang L. Huo X. On the exp-function method for constructing travelling wave solutions of nonlinear equations Nonlinear and Modern Mathematical Physics: Proceedings of the First International Workshop 2010 1212 Melville, NY, USA American Institute of Physics 280 285 MR2648971 ZBL1216.35109 Li J. B. Singular Traveling Wave Equations: Bifurcations and Exact Solutions 2013 Beijing, China Science Press Zhang L. Chen L.-Q. Huo X. The effects of horizontal singular straight line in a generalized nonlinear Klein-Gordon model equation Nonlinear Dynamics 2013 72 4 789 801 10.1007/s11071-013-0753-7 MR3064500 Hirota R. The Direct Method in Soliton Theory 2004 155 Cambridge, Mass, USA Cambridge University Press xii+200 Cambridge Tracts in Mathematics 10.1017/CBO9780511543043 MR2085332 Hietarinta J. Hirota's bilinear method and soliton solutions Physics AUC 2005 15, part 1 31 37 Ye Y. Wang L. Chang Z. He J. An efficient algorithm of logarithmic transformation to Hirota bilinear form of KdV-type bilinear equation Applied Mathematics and Computation 2011 218 5 2200 2209 10.1016/j.amc.2011.07.036 MR2831494 ZBL05992915 Zhou Z. Fu J. Li Z. An implementation for the algorithm of Hirota bilinear form of PDE in the Maple system Applied Mathematics and Computation 2006 183 2 872 877 10.1016/j.amc.2006.06.034 MR2290840 ZBL1110.65098 Zhou Z. Fu J. Li Z. Maple packages for computing Hirota's bilinear equation and multisoliton solutions of nonlinear evolution equations Applied Mathematics and Computation 2010 217 1 92 104 10.1016/j.amc.2010.05.012 MR2672566 ZBL1205.65281 Yang X.-D. Ruan H.-Y. A Maple package on symbolic computation of Hirota bilinear form for nonlinear equations Communications in Theoretical Physics 2009 52 5 801 807 10.1088/0253-6102/52/5/07 MR2641434 ZBL1186.35194 Hietarinta J. A search for bilinear equations passing Hirota's three-soliton condition. II. mKdV-type bilinear equations Journal of Mathematical Physics 1987 28 9 2094 2101 10.1063/1.527421 MR904425 ZBL0658.35081 Hietarinta J. A search for bilinear equations passing Hirota's three-soliton condition. III. Sine-Gordon-type bilinear equations Journal of Mathematical Physics 1987 28 11 2586 2592 10.1063/1.527750 MR913411 ZBL0658.35082 Hietarinta J. A search for bilinear equations passing Hirota's three-soliton condition. IV. Complex bilinear equations Journal of Mathematical Physics 1988 29 3 628 635 10.1063/1.528002 MR931466 ZBL0684.35082 Ma W.-X. Huang T. Zhang Y. A multiple exp-function method for nonlinear differential equations and its application Physica Scripta 2010 82 6 2-s2.0-78149402626 10.1088/0031-8949/82/06/065003 065003 Ma W.-X. Strampp W. Bilinear forms and Bäcklund transformations of the perturbation systems Physics Letters A 2005 341 5-6 441 449 10.1016/j.physleta.2005.05.013 MR2148263 ZBL1171.37332 Ma W.-X. Zhou R. Gao L. Exact one-periodic and two-periodic wave solutions to Hirota bilinear equations in ( 2 + 1 ) dimensions Modern Physics Letters A 2009 24 21 1677 1688 10.1142/S0217732309030096 MR2549829 ZBL1168.35426 Ma W.-X. Fan E. Linear superposition principle applying to Hirota bilinear equations Computers & Mathematics with Applications 2011 61 4 950 959 10.1016/j.camwa.2010.12.043 MR2770499 ZBL1217.35164 Ma W.-X. Zhang Y. Tang Y. Tu J. Hirota bilinear equations with linear subspaces of solutions Applied Mathematics and Computation 2012 218 13 7174 7183 10.1016/j.amc.2011.12.085 MR2880302 ZBL1245.35109 Ma W. X. Generalized bilinear differential equations Studies in Nonlinear Sciences 2011 2 4 140 144 Jimbo M. Miwa T. Solitons and infinite-dimensional Lie algebras Publications of the Research Institute for Mathematical Sciences 1983 19 3 943 1001 10.2977/prims/1195182017 MR723457 ZBL0571.35104 Ma W.-X. Lee J.-H. A transformed rational function method and exact solutions to the 3 + 1 dimensional Jimbo-Miwa equation Chaos, Solitons & Fractals 2009 42 3 1356 1363 10.1016/j.chaos.2009.03.043 MR2547030 ZBL1198.35231 Li Z. Dai Z. Abundant new exact solutions for the ( 3 + 1 ) -dimensional Jimbo-Miwa equation Journal of Mathematical Analysis and Applications 2010 361 2 587 590 10.1016/j.jmaa.2009.07.040 MR2568721 Kudryashov N. A. Sinelshchikov D. I. A note on “Abundant new exact solutions for the ( 3 + 1 ) -dimensional Jimbo-Miwa equation Journal of Mathematical Analysis and Applications 2010 371 1 393 396 10.1016/j.jmaa.2010.04.072 MR2661082 Chow S. N. Hale J. K. Methods of Bifurcation Theory 1981 New York, NY, USA Springer Wazwaz A.-M. Multiple-soliton solutions for the Calogero-Bogoyavlenskii-Schiff, Jimbo-Miwa and YTSF equations Applied Mathematics and Computation 2008 203 2 592 597 10.1016/j.amc.2008.05.004 MR2458974 ZBL1154.65366