Next Article in Journal
Coulomb Solutions from Improper Pseudo-Unitary Free Gauge Field Operator Translations
Previous Article in Journal
Privacy-Enhancing Security Protocol in LTE Initial Attack
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benjamin–Bona–Mahony Equation with Variable Coefficients: Conservation Laws

by
Ben Muatjetjeja
and
Chaudry Masood Khalique
*
International Institute for Symmetry Analysis and Mathematical Modelling, Department of Mathematical Sciences, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho 2735, South Africa
*
Author to whom correspondence should be addressed.
Symmetry 2014, 6(4), 1026-1036; https://doi.org/10.3390/sym6041026
Submission received: 22 August 2014 / Revised: 18 November 2014 / Accepted: 25 November 2014 / Published: 15 December 2014

Abstract

:
This paper aims to construct conservation laws for a Benjamin–Bona–Mahony equation with variable coefficients, which is a third-order partial differential equation. This equation does not have a Lagrangian and so we transform it to a fourth-order partial differential equation, which has a Lagrangian. The Noether approach is then employed to construct the conservation laws. It so happens that the derived conserved quantities fail to satisfy the divergence criterion and so one needs to make adjustments to the derived conserved quantities in order to satisfy the divergence condition. The conservation laws are then expressed in the original variable. Finally, a conservation law is used to obtain exact solution of a special case of the Benjamin–Bona–Mahony equation.

1. Introduction

The Benjamin–Bona–Mahony equation
u t + u x + u u x u x x t = 0
was, for the first time, studied by Benjamin et al. [1]. It is also known as the regularized long-wave equation and is applicable to shallow water waves and to the study of drift waves in plasma or the Rossby waves in rotating fluids. See, for example, [2] and references therein. It should be noted, however, that the BBM equation was also derived by Peregrine [3].
In the last two decades, various versions of the Benjamin–Bona–Mahony equation have been investigated in the literature. A general form of the Benjamin–Bona–Mahony equation is [2]
u t + α u x + β u u x + δ u x x t = 0
where α, β and δ are constants with the nonlinear and dispersion coefficients β ≠ 0 and δ < 0, respectively. For different values of the constants α, β and δ, Equation (2), results in various types of nonlinear equations which are very useful in the study of various physical phenomena. See [2] and references therein. However, in certain cases, the physical situations of the problem dictate us to consider nonlinear equations with variable coefficients [2,48]. Recently, in [2] the variable coefficients version of the Benjamin–Bona–Mahony equation
u t + α ( t ) u x + β ( t ) u u x δ ( t ) u x x t = 0
was investigated and some exact solutions were obtained using the classical Lie group method [9].
In this study, we consider the Equation (3), but with δ being an arbitrary constant, namely
u t + α ( t ) u x + β ( t ) u u x + δ ( t ) u x x t = 0
where α(t) and β(t) are arbitrary functions of t. The objective of the study is to classify the Noether operators and to construct conservation laws for the Equation (4).
Conservation laws are mathematical expressions of the physical laws, such as conservation of energy, mass, momentum and so on. They play a very important role in the solution and reduction of partial differential equations. Conservation laws have been widely used to investigate the existence, uniqueness and stability of solutions of nonlinear partial differential equations. This can be seen in the references [1012]. They have also been employed in the development and use of numerical methods (see for example, [13,14]). Recently, conserved vectors associated with Lie point symmetries have been used to find exact solutions of some partial differential equations [15] and systems of ordinary equations [16]. Noether theorem [17] gives us an elegant way to construct conservation laws provided a Lagrangian is known for an Euler–Lagrange equation. Thus, the knowledge of a Lagrangian is essential in this case.
It is worth mentioning that no Lagrangian exists for Equation (4) and as a result one can not invoke Noether theorem. However, an interesting approach is employed to construct conservation laws for Equation (4). It should be noted that the present approach, which we use here, fails to construct the conservation laws for Equation (3), i.e., when δ is a function of t.
The paper is organized as follows. In Section 2 we briefly give the fundamental notations and relations concerning the Noether symmetry approach. Section 3 obtains the Noether operators and the corresponding conservation laws for the Equation (4). In Section 4, a conservation law is used to obtain exact solution of a special case of the Benjamin–Bona–Mahony equation. Concluding remarks are presented in Section 5.

2. Fundamental Notations and Relations

Here we present some vital features of Noether operators concerning partial differential equations. These results will be used in Section 3. The reader is referred to [17,18] for further details.
Consider the vector field
X = τ ( t , x , z ) t + ξ ( t , x , z ) x + η ( t , x , z ) z
which has the second-order prolongation X[2] given by
X [ 2 ] = τ ( t , x , z ) t + ξ ( t , x , z ) x + η ( t , x , z ) z + ς t 1 z t + ς x 2 z x + ς t t 1 z t t + ς x x 2 z x x + ς t x 1 z t x +
where the expressions for ζ t 1 , ζ x 2 , ζ t x 1 and ζ x x 2 are given in [5]. The Euler–Lagrange operator is defined by
δ δ z = z D t z t D x z x + D t 2 z t t + D x 2 z x x + D x D t z x t
where the total differential operators are given by
D t = t + z t z + z t t z t + z x t z x +
D x = x + z x z + z x x z x + z x t z t +
Consider a partial differential equation of two independent variables, viz.,
E ( t , x , z , z x , z t , z t t , z x x , ) = 0
which has a second-order Lagrangian L, i.e., Equation (10) is equivalent to the Euler–Lagrange equation
δ L δ z = 0
Definition 1. The vector field X, of the form Equation (5), is called a Noether operator corresponding to a second-order Lagrangian L of Equation (10) if
X [ 2 ] ( L ) + { D t ( ξ 1 ) + D x ( ξ 2 ) } L = D t ( B 1 ) + D x ( B 2 )
for some gauge functions B1(t, x, z) and B2(t, x, z).
We now recall the following theorem.
Theorem 1. (Noether [17]) If X, as given in Equation (5), is a Noether point symmetry generator corresponding to a second-order Lagrangian L of Equation (10), then the vector T = (T1, T2) with components
T 1 = τ L + W δ L δ z t + D x ( W ) δ L δ z t x + D t ( W ) δ L δ z t t B 1
T 2 = ξ L + W δ L δ z x + D t ( W ) δ L δ z t x + D x ( W ) δ L δ z x x B 2
is a conserved vector for Equation (10) associated with the operator X, where W = η − ztτ − zxξ is the Lie characteristic function.

3. Conservation Laws of Equation (4)

Consider the Benjamin–Bona–Mahony Equation (4). Here, we note that Equation (4) does not admit a Lagrangian. Nevertheless, we can transform Equation (4) into a variational form by setting u = zx. Thus, Equation (4) becomes a fourth-order equation, namely
z t x + α ( t ) z x x + β ( t ) z x z x x + δ z x x x t = 0
It can easily be verified that a second-order Lagrangian of the Equation (15) is
L = 1 2 { δ z x x z x t α ( t ) z x 2 z t z x 1 3 β ( t ) z x 3 }
The substitution of L from Equation (16) into Equation (12) and splitting with respect to the derivatives of z, yields an overdetermined system of linear PDEs
τ z = 0 , η z = 0 , ξ z = 0 , τ x = 0 , ξ t = 0 , ξ x = 0 ,
η x x = 0 , η x t = 0 , β ( t ) τ + β ( t ) τ t = 0 ,
α ( t ) ξ 1 + β ( t ) η x + α ( t ) τ t = 0 ,
1 2 η t α ( t ) η x = B z 2 , 1 2 η x = B z 1 ,
B t 1 + B x 2 = 0.
After some lengthy calculations, the solution of the above system yields
τ = a ( t )
ξ = c 1
η = c 2 x + d ( t )
B 1 = 1 2 c 2 z + f ( t , x )
B 2 = 1 2 d z α ( t ) c 2 z + e ( t , x )
f t + e x = 0
β ( t ) a ( t ) + β ( t ) a ( t ) = 0
α ( t ) a ( t ) + α ( t ) a ( t ) + β ( t ) c 2 = 0
The analysis of Equations (28) and (29) prompts the following four cases:
Case 1. α(t) and β(t) arbitrary but not of the form contained in Cases 2–4.
In this case, we obtain two Noether point symmetries. These are given below together with their corresponding gauge functions:
X 1 = x , B 1 = f , B 2 = e , f t + e x = 0
X 2 = d ( t ) z , B 1 = f , B 2 = 1 2 d ( t ) z + e , f t + e x = 0
Invoking Theorem 1, and reverting back to the original variables, the two nontrivial conserved vectors associated with these two Noether point symmetries are, respectively,
T 1 1 = u 2 2 + δ 2 u u x x δ 2 u x 2 f
T 1 2 = 1 2 α ( t ) u 2 + 1 3 β ( t ) u 3 + δ u u x t δ 2 u x u t e
and
T 2 1 = 1 2 d ( t ) u δ 2 d ( t ) u x x f
T 2 2 = α ( t ) d ( t ) u 1 2 d ( t ) u t d x 1 2 β ( t ) d ( t ) u 2 δ d ( t ) u x t + 1 2 δ d ( t ) u x + 1 2 d ( t ) u d x e
Here it can be seen that the above conserved vectors do not satisfy the divergence condition, viz., DiTi|(4) = 0, as some excessive terms emerge that require some further analysis. By making a slight adjustment to these terms, it can be shown that these terms can be absorbed into the divergence condition. For,
D t ( T 1 1 ) + D x ( T 1 2 ) = δ 2 u u x x t δ 2 u x u x t f t e x = D t ( δ 2 u u x x f ) D x ( δ 2 u x u t + e )
hence
D t ( T 1 1 δ 2 u u x x + f ) + D x ( T 1 2 + δ 2 u x u t + e ) = 0.
We now redefine the conserved vectors in the parenthesis as:
T ¯ 1 1 = T 1 1 δ 2 u u x x + f = 1 2 u 2 δ 2 u x 2
T ¯ 1 2 = T 1 2 + δ 2 u x u t + e = α ( t ) 2 u 2 + β ( t ) 3 u 3 + δ u u x t
Thus, the modified conserved vectors T ¯ 2 1 and T ¯ 1 2 satisfy the divergence condition. We observe that the conserved vectors Equations (38) and (39) are local conserved vectors. Likewise, we can then construct the second pair of the conserved quantities associated with T 2 1 and T 2 2 as:
T ¯ 2 1 = d ( t ) 2 u δ 2 d ( t ) u x x
T ¯ 2 2 = α ( t ) d ( t ) u 1 2 d ( t ) u t d x β ( t ) 2 d ( t ) u 2 δ 2 d ( t ) u x t + δ 2 d ( t ) u x + 1 2 d ( t ) u d x
Note that the conserved quantities Equations (40) and (41) are nonlocal conserved vectors, and since d(t) is an arbitrary function of t, one obtains infinitely many nonlocal conserved vectors of Equation (4). A special case of Equations (40) and (41), when d(t) = 1, is
T ¯ 2 1 = u 2 δ 2 u x x
T ¯ 2 2 = α ( t ) u 1 2 u t d x β ( t ) 2 u 2 δ 2 u x t
Case 2. α(t) = γ, β(t) = λ, where γ and λ are non-zero constants.
This case provides us with three Noether symmetry generators, namely, X1, X2 given by the operators Equations (30) and (31) and X3 given by
X 3 = t with B 1 = f , B 2 = e , f t + e x = 0
The use of Noether conserved vector Equations (13) and (14) corresponding to X3 yields
T 3 1 = 1 2 γ u 2 1 6 λ u 3 + δ 2 u x x u t d x f
T 3 2 = γ u u t d x + 1 2 [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ u x t u t d x 1 2 δ u x u t t d x 1 2 δ u t 2 e
Again, the above conserved flow fails to satisfy the divergence criterion. Thus, by inheriting the same procedure as in Case 1, the nontrivial conserved flows associated with X1, X2 and X3 are
T ¯ 1 1 = 1 2 u 2 δ 2 u x 2
T ¯ 1 2 = γ 2 u 2 + λ 3 u 3 + δ u u x t
T ¯ 2 1 = d ( t ) 2 u δ 2 d ( t ) u x x
T ¯ 2 2 = γ d ( t ) u 1 2 d ( t ) u t d x λ 2 d ( t ) u 2 δ 2 d ( t ) u x t + δ 2 d ( t ) u x + 1 2 d ( t ) u d x
and
T ¯ 3 1 = 1 2 γ u 2 1 6 λ u 3
T ¯ 3 2 = γ u u t d x + 1 2 [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ u x t u t d x 1 2 δ u t 2
respectively. The conserved vector Equations (47) and (48) is a local conserved vector whereas the conserved vectors Equations (49) and (52) are nonlocal conserved quantities. It should be noted that one can use Equations (49) and (50) to construct infinitely many nonlocal conserved vectors.
Case 3. α(t) = γeσt, β(t) = λeσt, where γ, σ and λ are nonzero constants.
Here we get three Noether point symmetries operators, viz., X1, X2 given by the operators Equations (30) and (31) and X3 given by
X 3 = e σ t t with B 1 = f , B 2 = e , f t + e x = 0
The application of Theorem 1, with X3, gives
T 3 1 = 1 2 γ u 2 1 6 λ u 3 + δ 2 e σ t u x x u t d x f
T 3 2 = γ u u t d x + 1 2 e σ t [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ e σ t u x t u t d x 1 2 δ e σ t u x u t t d x + 1 2 δ e σ t u x u t d x 1 2 δ e σ t u t 2 e
and as before, the modified conserved vectors are given by
T ¯ 3 1 = T 3 1 δ 2 e σ t u x x u t d x = γ 2 u 2 λ 6 u 3
T ¯ 3 2 = T 3 2 1 2 δ σ e σ t u x u t d x + 1 2 δ e σ t u x u t t d x = γ u u t d x + 1 2 e σ t [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ e σ t u x t u t d x 1 2 δ e σ t u t 2
Thus, the nontrivial conserved flows associated with X1, X2 and X3 are, respectively,
T ¯ 1 1 = 1 2 u 2 δ 2 u x 2
T ¯ 1 2 = γ 2 e σ t u 2 + λ 3 e σ t u 3 + δ u u x t
T ¯ 1 2 = 1 2 d ( t ) u 1 2 δ d ( t ) u x x
T ¯ 2 2 = γ e σ t d ( t ) u 1 2 d ( t ) u t d x λ 2 e σ t d ( t ) u 2 δ 2 d ( t ) u x t + δ 2 d ( t ) u x + 1 2 d ( t ) u d x
T ¯ 3 1 = γ 2 u 2 λ 6 u 3
T ¯ 3 2 = γ u u t d x + 1 2 e σ t [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ e σ t u x t u t d x 1 2 δ e σ t u t 2
Case 4. α(t) = γ/a(t), β(t) = λ/a(t), where γ, λ are constants, with γ, λ ≠ 0 and a(t) an arbitrary function of t.
In this case, we obtain three Noether point symmetry generators, viz., X1, X2 given by Equations (30) and (31) and X3 given by
X 3 = a ( t ) t with B 1 = f , B 2 = e , f t + e x = 0
Following the above procedure, the nontrivial local and nonlocal conserved quantities corresponding to X1, X2 and X3, in this case, are
T ¯ 1 1 = δ 2 u 2 δ 2 u x 2
T ¯ 1 2 = γ 2 a ( t ) u 2 + λ 3 a ( t ) u 3 + δ u u x t
T ¯ 2 1 = d ( t ) 2 u δ 2 d ( t ) u x x
T ¯ 2 2 = γ a ( t ) d ( t ) u 1 2 d ( t ) u t d x λ a ( t ) d ( t ) u 2 δ 2 d ( t ) u x t + δ 2 d ( t ) u x + 1 2 d ( t ) u d x
T ¯ 3 1 = 1 2 γ u 2 1 6 λ u 3
T ¯ 3 2 = γ u u t d x + 1 2 a ( t ) [ u t d x ] 2 + 1 2 λ u 2 u t d x + δ a ( t ) u x t u t d x 1 2 δ a ( t ) u t 2
respectively.
Remark 1. Remark: It should be noted that since the Lagrangian Equation (16) is invariant under the spatial translation symmetry, this will give rise to the linear momentum conservation laws. Moreover, if α(t) = γ, β(t) = λ, where γ and λ are non-zero constants, then the corresponding Lagrangian Equation (16) is also invariant under the time translation symmetry and thus the linear momentum and energy are both conserved.

4. Exact Solution of Equation (4) for a Special Case Using Conservation Laws

First we recall a definition and theorem from [15] that we utilize in this Section.
Definition 2. Suppose that X is a symmetry of Equation (10) and T a conserved vector of Equation (10). Then if X and T satisfy
X ( T i ) + T i D j ( ξ j ) T j D j ( ξ j ) = 0 , i = 1 , 2
then X is said to be associated with T.
Define a nonlocal variable v by Tt = vx, Tx = −vt. Then using the similarity variables r, s, w with the generator X = s, we have Tr = vs, Ts = −vr, and the conservation law is then rewritten as
D r T r + D s T s = 0 ,
where
T r = T t D t ( t ) + T x D x ( r ) D t ( r ) D x ( s ) D x ( r ) D t ( s )
T s = T t D t ( s ) + T x D x ( S ) D t ( r ) D x ( s ) D x ( r ) D t ( s )
Theorem 2. An n-th order partial differential equation with two independent variables, which admits a symmetry X that is associated with a conserved vector T, is reduced to an ordinary differential equation of order n − 1; namely, Tr = k, where Tr is defined by Equation (72) for solutions invariant under X.
We now use the above definition and theorem to obtain an exact solution of one special case of Equation (4) by making use of its conservation law. We consider Case 2 of Section 3.
Recall that the Equation (4) with α(t) = γ, β(t) = λ, where γ and λ are non-zero constants, admits (among others)
X 1 = t , X 2 = x
and possesses the conservation law with conserved vector
T = ( 1 2 u 2 δ 2 u x 2 , γ 2 u 2 + λ 3 u 3 + δ u u x t )
that is associated with both X1 and X2. We now define the combination of X1 and X2 by X = ρX1+X2. Thus, the canonical coordinates of X are given by
s = x , r = ρ x t , u = u ( r )
where u = u(r) is the invariant solution under X if u satisfies Equation (4) with α(t) = γ, β(t) = λ. Employing Equation (72), the r-component of T given in Equation (75) is
T r = 1 2 u 2 α γ 2 u 2 δ α 2 2 u r 2 α λ 3 u 3 + δ α 2 u u r r
Using Theorem 2 above, Equation (77) can be written as
1 2 u 2 α γ 2 u 2 δ α 2 2 u r 2 α λ 3 u 3 + δ α 2 u u r r = k
By letting ur = p(u), we have ... Then Equation (78) reduces to the first order ordinary differential equation
p 1 2 u p = 1 p ( λ 3 δ α u 2 1 2 δ α 2 u + γ 2 δ α u + k δ α 2 u )
The integration of Equation (79) leads to the four parameter family of solutions
± ( λ u 3 / 3 δ α 2 k / α 2 δ u 2 / α 2 δ + γ u 2 / δ α + k 1 u 2 ) 1 2 d u = ρ x t + k 2
of Equation (4) invariant under X = ρX1 + X2.

5. Concluding Remarks

In this paper we derived the conservation laws for the Benjamin–Bona–Mahony equation with variable coefficients Equation (4). This equation does not have a Lagrangian. In order to use Noether theorem, we transformed Equation (4) to a fourth-order Equation (15) which admits a Lagrangian. We then employed the Noether theorem to derive the conservation laws for Equation (15). The derived conserved vectors needed further adjustment to satisfy the divergence criterion. By reverting back to the original variable u, we constructed the conserved vectors for our third-order Benjamin–Bona–Mahony equation with variable coefficients Equation (4). The conserved vectors consisted of some local and infinite number of nonlocal conserved vectors. Finally, we obtained an exact solution for the Benjamin–Bona–Mahony equation using the double reduction theory. The importance of finding conservation laws and their physical applications are mentioned in the paper.

Acknowledgments

The authors thank the anonymous referees whose comments helped to improve the paper. Ben Muatjetjeja would like to thank the North-West University, Mafikeng Campus for their financial support.

Author Contributions

Ben Muatjetjeja and Chaudry Masood Khalique worked together in the derivation of the mathematical results. Both authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Benjamin, T.B.; Bona, J.L.; Mahony, J.J. Model equations for long waves in nonlinear dispersive systems. Philos. Trans. R. Soc. Lond. A 1972, 272, 47–78. [Google Scholar]
  2. Singh, K.; Gupta, R.K.; Kumar, S. Benjamin–Bona–Mahony (BBM) equation with variable coefficients: Similarity reductions and Painlevé analysis. Appl. Math. Comput 2011, 217, 7021–7027. [Google Scholar]
  3. Peregrine, D.H. Calculations of the development of an undular bore. J. Fluid Mech 1966, 25, 321–330. [Google Scholar]
  4. Molati, M.; Khalique, C.M. Lie symmetry analysis of the time-variable coefficient B-BBM equation. Adv. Differ. Equ 2012, 212. [Google Scholar] [CrossRef]
  5. Johnpillai, A.G.; Khalique, C.M. Variational approaches to conservation laws for a nonlinear evolution equation with time dependent coefficients. Quaest. Math 2011, 34, 235–245. [Google Scholar]
  6. Wang, M.L.; Wang, Y.M. A new Bäcklund transformation and multi-soliton solutions to the KdV equation with general variable coefficients. Phys. Lett. A 2001, 287, 211–216. [Google Scholar]
  7. Engui, F. Auto-Bäcklund transformation and similarity reductions for general variable coefficient KdV equations. Phys. Lett. A 2002, 294, 26–30. [Google Scholar]
  8. Singh, K.; Gupta, R.K. Lie symmetries and exact solutions of a new generalized Hirota-Satsuma coupled KdV system with variable coefficients. Int. J. Eng. Sci 2006, 44, 241–255. [Google Scholar]
  9. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer-Verlag: New York, NY, USA, 1993. [Google Scholar]
  10. Lax, P.D. Integrals of nonlinear equations of evolution and solitary waves. Commun. Pure Appl. Math 1968, 21, 467–490. [Google Scholar]
  11. Benjamin, T.B. The stability of solitary waves. Proc. R. Soc. Lond. A 1972, 328, 153–183. [Google Scholar]
  12. Knops, R.J.; Stuart, C.A. Quasiconvexity and uniqueness of equilibrium solutions in nonlinear elasticity. Arch. Ration. Mech. Anal 1984, 86, 234–249. [Google Scholar]
  13. LeVeque, R.J. Numerical Methods for Conservation Laws; Birkhäuser: Basel, Switzerland, 1992. [Google Scholar]
  14. Godlewski, E.; Raviart, P.A. Numerical Approximation of Hyperbolic Systems of Conservation laws; Springer: Berlin, Germany, 1996. [Google Scholar]
  15. Sjöberg, A. Double reduction of PDEs from the association of symmetries with conservation laws with applications. Appl. Math. Comput 2007, 184, 608–616. [Google Scholar]
  16. Muatjetjeja, B.; Khalique, C.M. Lie group classification for a generalized coupled Lane-Emden system in dimension one. East Asian J. Appl. Math 2014, 4, 301–311. [Google Scholar]
  17. Noether, E. Invariante Variationsprobleme. König. Gesell. Wissen. Göttingen. Math Phys. Klasse 1918, 2, 235–257. [Google Scholar]
  18. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1994–1996; Volume 1–3. [Google Scholar]

Share and Cite

MDPI and ACS Style

Muatjetjeja, B.; Khalique, C.M. Benjamin–Bona–Mahony Equation with Variable Coefficients: Conservation Laws. Symmetry 2014, 6, 1026-1036. https://doi.org/10.3390/sym6041026

AMA Style

Muatjetjeja B, Khalique CM. Benjamin–Bona–Mahony Equation with Variable Coefficients: Conservation Laws. Symmetry. 2014; 6(4):1026-1036. https://doi.org/10.3390/sym6041026

Chicago/Turabian Style

Muatjetjeja, Ben, and Chaudry Masood Khalique. 2014. "Benjamin–Bona–Mahony Equation with Variable Coefficients: Conservation Laws" Symmetry 6, no. 4: 1026-1036. https://doi.org/10.3390/sym6041026

APA Style

Muatjetjeja, B., & Khalique, C. M. (2014). Benjamin–Bona–Mahony Equation with Variable Coefficients: Conservation Laws. Symmetry, 6(4), 1026-1036. https://doi.org/10.3390/sym6041026

Article Metrics

Back to TopTop