On Galilean Invariant and Energy Preserving BBM-Type Equations

Abstract

We investigate a family of higher-order Benjamin–Bona–Mahony-type equations, which appeared in the course of study towards finding a Galilei-invariant, energy-preserving long wave equation. We perform local symmetry and conservation laws classification for this family of Partial Differential Equations (PDEs). The analysis reveals that this family includes a special equation which admits additional, higher-order local symmetries and conservation laws. We compute its solitary waves and simulate their collisions. The numerical simulations show that their collision is elastic, which is an indication of its $S-$integrability. This particular PDE turns out to be a rescaled version of the celebrated Camassa–Holm equation, which confirms its integrability.

1. Introduction

The famous Benjamin–Bona–Mahony (BBM) equation was derived for the first time D. Peregrine [1], reintroduced later in [2]. In scaled variables, it is given by

$$u_t+u_x+u{u}_x-u_{xxt}=0,$$

(1)

where subscripts ${(-)}_x$ and ${(-)}_t$ denote ${\partial }_x$ and ${\partial }_t$—the partial derivative with respect to the spatial x and temporal t variables respectively. One of the main physical drawbacks of the BBM equation is the absence of the Galilean invariance. This question has already been raised in [3] and a remedy was proposed. It consisted in adding higher order terms that help recover this symmetry. However, this approach led to the simultaneous loss of another physically important property—the presence of the energy conservation law. Various investigations around these questions [4] led the authors of this communication to consider the following family of partial differential equations (PDEs):

$$u_t+u_x+u{u}_x-u_{xxt}-\alpha (u{u}_{xxx}+2u_x{u}_{xx})=0,$$

(2)

with $\alpha \in \mathbb{R}$. It is obvious that PDE (2) reduces to the BBM (1) when $\alpha =0$. It is straightforward to see that the parameter $\alpha$ in Equation (2) cannot be removed by a rescaling or another simple equivalence transformation (see, e.g., [5]). Below, we study symmetries, conservation laws, and solitary wave solutions to family (2). The conservation laws of a variable coefficients BBM equation were studied in [6].

In particular, we show that additional structures and exceptional elastic collision-type numerical solution behavior arise for the PDE (2) corresponding to $\alpha =1/3$:

$$u_t+u_x+u{u}_x-u_{xxt}-\frac{1}{3}u{u}_{xxx}-\frac{2}{3}{u}_x{u}_{xx}=0.$$

(3)

We will refer to this integrable PDE as the eBBM${}_{\frac{1}{3}}$ equation.

The present manuscript is organized as follows. In Section 2 the symmetries, Hamiltonian structure, and conservation laws to Equation (2) are studied for generic values of the parameter $\alpha$. The potential symmetries are studied in Section 3. The classification of symmetries and conservation laws with respect to parameter $\alpha$ is done in Section 4. Finally, the solitary wave collisions are briefly discussed in Section 5. The main conclusions of this study are outlined in Section 6.

2. Symmetries, Conservation Laws, and Hamiltonian Structure of the $\mathbf{\alpha }$-Family

The point symmetries of the family of PDEs (2) holding for all $\alpha$ are given by the generators

$${\mathrm{X}}_1={\partial }_t,{\mathrm{X}}_2={\partial }_x,$$

(4)

$${\mathrm{X}}_3=(\alpha -1)t{\partial }_t+\alpha t{\partial }_x+(1+(1-\alpha )u){\partial }_u.$$

(5)

In particular, when $\alpha =0$, the generator ${\mathrm{X}}_3$ yields the scaling symmetry group holding for the BBM Equation (1), and when $\alpha =1$, ${\mathrm{X}}_3$ yields the Galilei group. The latter case is the only Galilei-invariant representative of the PDE family (2). A direct computation shows that there are no higher-order symmetries for any member of the family (2) depending on derivatives of u up to the second-order, i.e., with generators of the form

$$\widehat{\mathrm{X}}=\zeta (x,t,u,u_t,u_x,u_{tt},u_{xt},u_{xx}){\partial }_u,$$

except the point symmetries (4) and (5). Moreover, one can show there are no symmetries holding for an arbitrary $\alpha$ with components depending on $x$, $t$, $u$, and x-derivatives of u up to order five. (It will be shown below that the eBBM${}_{\frac{1}{3}}$ Equation (3) does admit a higher-order symmetry). In the case of an arbitrary $\alpha$, all PDEs in the $\alpha$-family (2) admit local conservation laws with multipliers

$${\mathrm{\Lambda }}_1=1,{\mathrm{\Lambda }}_2=u,{\mathrm{\Lambda }}_3=u^2-2u_{xt}-\alpha (u_x^2+2u{u}_{xx}),$$

(6)

corresponding respectively to the conservation of the momentum and the energy:

$$\mathcal{M}\left(t\right)={\int }_a^b(u-u_{xx})\mathrm{d}x,\mathcal{E}\left(t\right)={\int }_a^b\frac{1}{2}(u^2+u_x^2)\mathrm{d}x,$$

as well as the $\alpha$-dependent conserved quantity

$${\mathcal{N}}_{\alpha }\left(t\right)={\int }_a^b\left(\frac{1}{3}{u}^3+(\alpha u-1)u_x^2\right)\mathrm{d}x.$$

The latter, in a linear combination with energy $\mathcal{E}\left(t\right)$, yields a conserved quantity

$${\mathscr{H}}_{\alpha }\left(t\right)=\frac{1}{2}{\int }_a^b\left(u^2+\frac{1}{3}{u}^3+\alpha u{u}_x^2\right)\mathrm{d}x,$$

which defines a Hamiltonian for the whole $\alpha$-family (2). Indeed, the PDEs (2) can be written in the Hamiltonian form

$$u_t=\mathbb{J}\frac{\delta {\mathscr{H}}_{\alpha }}{\delta u},$$

with the co-symplectic operator $\mathbb{J}$ shared with the original BBM Equation (1):

$$\mathbb{J}≔{(𝟙-{\partial }_x^2)}^{-1}·(-{\partial }_x).$$

3. The Potential $\mathbf{\alpha }$-Family and Its Lagrangian Structure

It is possible to show that similarly to the Korteweg–de Vries (KdV) and the BBM equations, the PDEs of the $\alpha$-family (2), as it stands, does not have a self-adjoint Fréchet derivative, and hence does not directly arise from a variational principle for any Lagrangian functional. This can also be seen from the obvious differences between conservation law multipliers (6) and the evolutionary forms of symmetries (4) and (5), which would contradict the first Noether’s theorem (cf. [7]).

However, the $\alpha$-family (2) has a potential Lagrangian formulation. In fact, upon the introduction of the potential variable $w=w(x,t)$ defined by

$$u\stackrel{\mathrm{def}}{:=}{w}_x,$$

(7)

the $\alpha$-family (2) is written in the potential form

$$w_{tx}+w_{xx}+w_x{w}_{xx}-w_{xxxt}-\alpha (w_x{w}_{xxxx}+2w_{xx}{w}_{xxx})=0,$$

(8)

which, for every $\alpha$, is nonlocally related to the corresponding PDE (2) (see, e.g., [8]). One can show that the linearization of Equation (8) is self-adjoint. The homotopy formula ([7] Theorem 5.92) can be used to construct the Lagrangian, which up to equivalence is readily found to be

$$ℒ\left[w\right]\stackrel{\mathrm{def}}{:=}-\frac{1}{2}{\int }_{\mathbb{R}}\left(\frac{1}{3}{w}_x^3+w_x^2+w_x{w}_t+w_{xx}{w}_{xt}+\alpha w_x{w}_{xx}^2\right)\mathrm{d}x.$$

(9)

It is straightforward to check that the Euler–Lagrange equation $\delta ℒ/\delta w=0$ yields the potential Equation (8) for each value of the parameter $\alpha$.

It is of interest to briefly consider symmetries and conservation laws of the potential $\alpha$-family (8). The first reason for that is that the potential $\alpha$-family Equation (8) are nonlocally related to the PDEs (2), and hence the local symmetries and conservation laws of (2) vs. (8) can differ (cf. [8]). Secondly, since the PDEs (8) arise from a variational principle, the first Noether’s theorem will provide an evident relation between their local symmetries and conservation laws. Seeking local symmetries of the potential $\alpha$-family Equation (8) in evolutionary form

$$\widehat{\mathrm{X}}=\zeta (x,t,w,w_x,w_{xx},w_{xxt},w_{xxx}){\partial }_w,$$

(10)

holding for all $\alpha \in \mathbb{R}$, we find the admitted symmetry generator components:

$$\begin{array}{c}{\zeta }_1=w_t,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \\ {\zeta }_2=w_x,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \end{array}$$

(11)

$$\begin{array}{c}{\zeta }_3=x-(\alpha -1)(w+t{w}_t)-\alpha t{w}_x,\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \end{array}$$

(12)

$$\begin{array}{c}{\zeta }_4=2(x-t{w}_x)+{(\alpha -1)}^{2}t(w_{xx}^2+2w_x{w}_{xxx})\hfill \\ +(\alpha -1)\left(2w+t(2w_{xxt}+w_{xx}^2+2w_x{w}_{xxx}-w_x^2)\right),\hfill \end{array}$$

(13)

$$\begin{array}{c}{\zeta }_5=w_x^2-2w_{xxt}-\alpha (w_{xx}^2+2w_x{w}_{xxx}),\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \end{array}$$

(14)

$$\begin{array}{c}{\zeta }_F=F\left(t\right).\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \end{array}$$

(15)

The generators (10) with ${\zeta }_1$ and ${\zeta }_2$ correspond to the time and space translation point symmetries ${\mathrm{X}}_1={\partial }_t$, ${\mathrm{X}}_2={\partial }_x$ of the $\alpha -$ family (2). The local symmetry with ${\zeta }_3$ (11) corresponds to the point symmetry ${\mathrm{X}}_3=\alpha t{\partial }_x+(\alpha -1)t{\partial }_t+(x-(\alpha -1)w){\partial }_w$ of (2), which includes the Galilei group only when $\alpha =1$. The generator with ${\zeta }_4$ (13) is generally a higher-order symmetry generator, which also degenerates into the Galilei group ${\mathrm{X}}_4=t{\partial }_x+{\partial }_w$ when $\alpha =1$. For $\alpha \ne 1$, the generator ${\zeta }_4$ explicitly involves the potential variable $w$, and hence is a nonlocal higher-order symmetry generator of the corresponding PDE from the $\alpha$-family (2) (cf. [8]). The local generator ${\widehat{\mathrm{X}}}_5$ with symmetry component ${\zeta }_5$ (14) is a higher-order symmetry of the potential equation for each $\alpha$. The local symmetry generator with ${\zeta }_F$ (15) admitted by the potential family (8) is its point symmetry $w\leftarrow w+F\left(t\right)$, with no projection on the space of variables $(x,t,u)$ of Equation (2), related to the non-uniqueness of the definition (7) of the potential up to an arbitrary function of time.

Since the potential family (8) arises from a Lagrangian, by Noether’s theorem, its variational symmetries yield conservation law multipliers. A brief computation yields

$${\mathrm{\Lambda }}_1={\zeta }_1,{\mathrm{\Lambda }}_2={\zeta }_2,{\mathrm{\Lambda }}_5={\zeta }_5,{\mathrm{\Lambda }}_F={\zeta }_F$$

which are local conservation law multipliers for all PDEs in the potential $\alpha$-family (8), whereas the local symmetries with ${\zeta }_3$ and ${\zeta }_4$ are not variational, and thus have no corresponding conservation laws. The multiplier ${\mathrm{\Lambda }}_F$ corresponds to a conservation law

$$\begin{array}{c}{\mathcal{D}}_t\left(F\left(t\right)(w_x-w_{xxx})\right)+{\mathcal{D}}_x(-F^{\prime }\left(t\right)(w-w_{xx})+\hfill \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}F\left(t\right)\left(w_x(1-\alpha w_{xxx})+\frac{1}{2}(w_x^2-\alpha w_{xx}^2)\right))=0,\end{array}$$

(16)

holding for an arbitrary function $F\left(t\right)$, which is a nonlocal conservation law of the $\alpha$-family (2) when $F\left(t\right)\ne 1$. Indeed, the spatial flux in (16) is explicitly dependent on the potential $w=\int u\mathrm{d}x$ itself, and hence is not equivalent to any local expression in terms of u. The conserved density in (16),

$${\rho }_F\stackrel{\mathrm{def}}{:=}F\left(t\right)(w_x-w_{xxx})\equiv F\left(t\right)(u-u_{xx})$$

is, however, a local quantity in terms of the dependent variable u of (8).

An additional important conservation law arises for the potential $\alpha -$ family (8). In fact, the left-hand side of each Equation (8) is a total $x-$ derivative:

$${\mathcal{D}}_x\left(w_t+w_x+\frac{1}{2}{w}_x^2-w_{xxt}-\alpha \left(\frac{1}{2}{w}_{xx}^2+w_x{w}_{xxx}\right)\right)=0,$$

which can be written as

$$w_t+w_x+\frac{1}{2}{w}_x^2-w_{xxt}-\alpha \left(\frac{1}{2}{w}_{xx}^2+w_x{w}_{xxx}\right)=C\left(t\right),$$

(17)

for an arbitrary function $C\left(t\right)$.

In addition to the potential form (8), other potential systems non-locally related to the $\alpha$-family (2) can be constructed, using, for example, the linearly independent conservation law multipliers (6), to obtain independent singlet, couplet, and triplet potential systems [8,9]. Their study may lead to new analytical results of $\alpha$-family (2), including the original BBM Equation (1).

4. Symmetry and Conservation Law Classifications for the $\mathbf{\alpha }$-Family

Since the PDE family (2) and its potential form (8) involve a parameter $\alpha$, its properties can be classified according to $\alpha$ (see e.g., [8] and references therein). In Section 2, symmetries and conservation laws that arise for an arbitrary $\alpha$ were listed; we now consider symmetry and conservation law classifications. First, classifying local conservation laws of the $\alpha$-family (2) with third-order multipliers, the following cases are distinguished:

• In the general case with arbitrary $\alpha$, as reported above, only three common local conservation laws arise, with multipliers (6);
• An additional local conservation law:

  $$\begin{array}{c}{\mathcal{D}}_t\left(\left(t-x+\frac{1}{2}tu\right)(u-u_{xx})\right)+\hfill \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{\mathcal{D}}_x(\frac{1}{3}t{u}^3-\left(\frac{1}{2}x+t(u_{xx}-1)\right)u^2+\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\left((t-x)(u_{xx}-1)+\frac{1}{2}(u_x+t{u}_{tx})\right)u+\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\frac{1}{2}\left(2+t{u}_t+(x-t)u_x\right)u_x)=0\end{array}$$

  (18)

  with multiplier

  $${\mathrm{\Lambda }}_4=x-t(u+1)$$

  arises when $\alpha =1$; this is the case of the so-called extended BBM (eBBM) equation (see [4]);
• The case when $\alpha =\frac{1}{3}$, that is, the eBBM${}_{\frac{1}{3}}$ Equation (3), turns out to be the only other special case when additional conservation laws arise.

In terms of local conservation laws in the specified multiplier ansatz, the eBBM${}_{\frac{1}{3}}$ Equation (3) can be shown to admit and three additional local conservation law multipliers

$$\begin{array}{c}{\mathrm{\Lambda }}_4^{\left(\frac{1}{3}\right)}=5u^3+(9-4u_{xx})u^2\hfill \\ \hspace{1em}\hspace{1em}-(u_x^2+6u_{xt})u+6(u_x{u}_t+3u_{tt}),\hfill \end{array}$$

(19)

$$\begin{array}{c}{\mathrm{\Lambda }}_5^{\left(\frac{1}{3}\right)}={\left(2(u-u_{xx})+3\right)}^{-1/2},\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \end{array}$$

(20)

$$\begin{array}{c}{\mathrm{\Lambda }}_6^{\left(\frac{1}{3}\right)}={\displaystyle \frac{u_{xxxx}-\frac{1}{2}(2u+3)}{{\left(2(u-u_{xx})+3\right)}^{5/2}}}+\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hfill \\ \hspace{1em}\hspace{1em}\frac{5}{2}{\displaystyle \frac{u_{xxx}^2+u_x(u_x-2u_{xxx})}{{\left(2(u-u_{xx})+3\right)}^{7/2}}}.\hfill \end{array}$$

(21)

The local conservation law with the multipliers (19)–(21) have the form:

$${\mathcal{D}}_t{\mathcal{T}}_i\left[u\right]+{\mathcal{D}}_x{\mathcal{F}}_i\left[u\right]=0,i=4,5,6,$$

with the corresponding conserved quantities

$${\mathcal{Q}}_i^{\left(\frac{1}{3}\right)}\left(t\right)\stackrel{\mathrm{def}}{:=}{\int }_a^b{\mathcal{T}}_i\left[u\right]\mathrm{d}x.$$

The conserved densities for the respective conservation laws are given by

$$\begin{array}{cc}\hfill {\mathcal{T}}_4\left[u\right]& \stackrel{\mathrm{def}}{:=}5u^4+12u^3+(26u_x^2+4u_{xx}^2)u^2+\hfill \\ \hfill & (u_x^2{u}_{xx}+24u_x^2+36u_{tt})u-36u_{tt}{u}_{xx},\hfill \\ \hfill {\mathcal{T}}_5\left[u\right]& \stackrel{\mathrm{def}}{:=}\sqrt{2(u-u_{xx})+3},\hfill \\ \hfill {\mathcal{T}}_6\left[u\right]& \stackrel{\mathrm{def}}{:=}{\left(2(u-u_{xx})+3\right)}^{-5/2}·\left(2\right(2(u-u_{xx}){+3)}^2+\hfill \\ \hfill & (2u+3)\left(2(u-u_{xx})+3)+3(u_{xxx}^2-u_x^2)\right).\hfill \end{array}$$

It is also possible to show by a direct verification that the exceptional eBBM${}_{\frac{1}{3}}$ Equation (3) admits a third-order local evolutionary symmetry generator:

$${\widehat{\mathrm{X}}}_4=\frac{u_x-u_{xxx}}{{(2(u-u_{xx})+3)}^{3/2}}{\partial }_u.$$

(22)

The existence of this higher-order symmetry and the above higher-order conservation laws point at the possibility of integrability of the eBBM${}_{\frac{1}{3}}$ Equation (3), in particular, in the light of the Fokas conjecture in [10], and the similarities between the eBBM${}_{\frac{1}{3}}$ and the integrable short pulse equation [4].

5. Solitary Waves Collision

Motivated by the exceptional symmetry and conservation law properties of eBBM${}_{\frac{1}{3}}$ Equation (3), we decided to study the interaction of solitary waves in this equation. For the sake of comparison, we shall also perform the same simulations in the eBBM equation given by (2) with $\alpha =1$. This choice of $\alpha$ is also interesting because it is a Galilei-invariant model. For extended numerical comparisons, see [4].

First of all, the solitary waves to Equation (2) have to be computed. For this purpose we employ the Petviashvili iterative method [11,12,13]. Then, in order to simulate the $\alpha$-family dynamics, we employ a Fourier-type pseudo-spectral discretization with adaptive higher-order embedded Runge–Kutta temporal discretization [14]. We employ $\mathrm{16,384}$ Fourier modes. Since we deal with unidirectional models, we simulate the overtaking collision of two solitary waves. For this purpose we generate two right-going solitary wave profiles, corresponding to propagation velocities $c^{+}=1.4$ and $c^{-}=1.1$. Then, they are placed in a periodic computational domain $[-140,140]$. The larger (and faster) solitary wave is initially placed at $x=-50$ and the smaller one at $x=50$. The initial condition is depicted in Figure 1. Then, this initial configuration is propagated by both models from $t=0$ to $t=500$. The final simulation time is chosen so that the first overtaking collision takes place and solitary waves have enough time to separate again. A zoom on the numerical solution at $t=500$ is shown in Figure 2. We can see that the interaction in the eBBM${}_{\frac{1}{3}}$ Equation (3) is elastic, while the choice of $\alpha =1$ (the energy-preserving Galilei-invariant eBBM equation) yields an inelastic collision.

6. Discussion and Conclusions

In this communication, a family of PDEs (2) with one free parameter was considered. This family was inspired by our quest for Galilei-invariant and energy-preserving higher order analogs of the classical BBM equation [2]. The detailed symmetry and conservation law analysis allowed us to extract one particular member (3) of this family, which admited a higher-order local symmetry (22) and local conservation laws with higher-order multipliers (19)–(21). Moreover, our spectrally accurate numerical simulations showed that the solitary wave interaction in Equation (3) was elastic in contrast to all other values of the free parameter $\alpha$. This fact allowed us to formulate the conjecture of the integrability of eBBM${}_{\frac{1}{3}}$ Equation (3). This conjecture turned out to be true because a close inspection of Equation (3) showed that it was a rescaled version of the celebrated Camassa–Holm (CH) equation [15,16]. Indeed, the eBBM${}_{\frac{1}{3}}$ equation was related to the CH by a local scaling transformation. For example, a time scaling $t=3\tau$ and the change of notation $\tau \to t$ maps the PDE (3) into the CH equation with parameter $\kappa =\frac{3}{2}$:

$$u_t+3u_x+3u{u}_x-2u_x{u}_{xx}-u_{xxt}-u{u}_{xxx}=0.$$

As a result, the present study can be considered as a symmetry-based derivation of a famous PDE inspired by physically sound modeling of the unidirectional nonlinear wave propagation.