Journal of Mathematics, Physics and Mechanics

INTRODUCTION

Nonlinear dispersive wave equations describe a broad class of phenomena in which the steepening effect of nonlinearity is balanced by the spreading effect of dispersion [1]. The Korteweg–de Vries (KdV) equation is the archetypal model of this class. It was derived for shallow-water waves and later found to govern solitary waves in plasma, optical, and lattice settings, and it is notable for its complete integrability and infinite family of conservation laws [2].

Despite this success, the KdV equation suffers from a well-documented weakness: its linear dispersion relation makes high-wavenumber components poorly behaved, so the model is not uniformly stable in its short-wave content. To remedy this defect, Benjamin, Bona, and Mahony proposed in 1972 an alternative model, now known as the Benjamin–Bona–Mahony equation, or the regularized long-wave equation [3]. The BBM equation replaces the third-order spatial derivative of the KdV equation with a mixed third-order derivative, $u_t-u_{xxt}+u_x+u{u}_x=0$, which yields a bounded phase velocity and removes the short-wave instability [3]. Because of this structural change, the equation is better suited than the KdV model for describing unidirectional propagation of small-amplitude long waves in shallow water, and it has become a standard testbed for the analysis and numerical simulation of nonlinear dispersive systems [4,5].

Over the past five decades, the BBM equation has generated a large and diverse literature. Rigorous results establish local and global well-posedness in classical Sobolev spaces [6,7], identify its three independent conservation laws [8], and characterize the stability of its solitary waves [9,10]. At the same time, the equation has motivated an enormous body of work on exact and numerical solution methods, fractional and stochastic generalizations, and applications ranging from undular bores and tsunamis to ion-acoustic and ocean waves [11-13].

A notable feature of the recent literature is the proliferation of reported exact solutions obtained by heuristic ansatz methods, many of which fail to satisfy the governing equation when checked by direct substitution. This situation has prompted calls for systematic validation and has motivated new work using symbolic computation to derive and certify exact solutions [14]. Such activity, together with the rapid adoption of machine-learning techniques and the continued growth of fractional and stochastic formulations, makes the present an opportune moment for a critical review.

The purpose of this article is threefold. First, we synthesize the mathematical foundations of the BBM equation, from its physical derivation to the modern theory of well-posedness and stability. Second, we provide a critical evaluation of analytical and numerical methods, identifying persistent pitfalls and methodological trends. Third, we highlight open problems and promising directions for future research, with particular attention to literature published in 2025 and 2026. Throughout, we cite the primary sources so that readers can verify claims directly. The organization of the review is as follows. Section 2 covers derivation and mathematical structure. Section 3 addresses solitary and periodic waves. Sections 4 and 5 survey exact and numerical methods, respectively. Sections 6 and 7 treat generalized, fractional, and stochastic variants. Section 8 summarizes applications, and Section 9 concludes with trends and open problems.

Derivation, Model Fidelity, and Mathematical Structure

Physical derivation

The BBM equation originates in the theory of long waves on shallow water. Peregrine's early computations of undular bore development already indicated that corrections to the linearized long-wave model are needed to reproduce the observed waveform [1]. Benjamin, Bona, and Mahony derived their model within the same asymptotic framework as the KdV equation, using small-amplitude and long-wavelength expansions, but with a different ordering of the dispersive terms [3]. The resulting equation,

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

has the linear dispersion relation $\omega =\frac{k}{1+k^2}$, which stays bounded and therefore avoids the instability of the KdV relation ω = k −k³ at large wavenumbers [3]. A useful one-parameter generalization, $u_t+u_x+bu{u}_x-u_{xxt}=0$, allows the nonlinearity to be tuned: the coefficient $b<0$ corresponds to the focusing regime and $b>0$ to the defocusing regime, and the two regimes display a structural duality in their exact solution families [14-16]. The modified long-wave equation, in which the quadratic nonlinearity is replaced by a cubic one, is usually called the modified BBM (MBBM) equation and shares the same structural properties [17].

Model Fidelity and Comparison with Kdv

A central question is whether the BBM equation is as faithful a model of water waves as the KdV equation. Comparison studies show that the two models agree to the same order of approximation for smooth, long-wave data, but differ in their high-frequency behavior [18]. The BBM equation is sometimes described as only approximately integrable, in the sense that it possesses a restricted set of conservation laws and is not solvable by the classical inverse scattering transform in general [18]. Attempts have nevertheless been made to treat the regularized long-wave equation within an inverse scattering framework [19]. The symmetric regularized long-wave equation, which is isotropic in the spatial variable, arises in the description of weakly nonlinear ion-acoustic and space-charge waves and supports hyperbolic-secant-squared solitary waves [12]. Related symmetric models for shallow-water waves were examined in later work [20].

Conservation Laws and Integrability

A defining structural result is that the BBM equation admits exactly three independent local conservation laws, in contrast to the infinite family of the KdV equation. Olver proved this fact using the theory of Euler operators [8]. The cohomological analysis of Duzhin and Tsujishita confirmed that no further nontrivial conservation laws exist [21], and independent derivations reached the same conclusion [22]. These three invariants, commonly associated with mass, momentum, and a quadratic energy functional, play a central role in both the stability theory and the design of conservative numerical schemes. The fact that the BBM equation is not completely integrable, but is close to integrable, makes it a valuable test case for perturbation and asymptotic methods [18].

Well-Posedness

The initial value problem for the BBM equation is well-posed in classical Sobolev spaces. Building on the foundational analysis of the KdV initial value problem [23], early work established existence and uniqueness of periodic solutions [7], global existence in higher dimensions [4], and global well-posedness in arbitrary spatial dimensions under suitable conditions [6]. Decay of solutions in the energy space was proved by Amick, Bona, and Schonbek [24] and refined in subsequent studies [25,26]. The long-period limit, relating periodic and nonperiodic solutions, was studied by Chen [27]. Recent work has sharpened the regularity theory, proving unique continuation properties for the linearized equation [28,29] and, in higher dimensions, Carleman-type estimates that imply strong unique continuation [30]. Hong and Ponce revisited the regularity and continuation properties of solutions to the BBM equation, confirming that the model is well suited to the modern theory of dispersive equations [31]. For the quartic generalized BBM equation, scattering of small solutions has recently been established [32].

Solitary Waves: Existence, Stability, and Instability

Existence and Analyticity

The BBM equation supports solitary waves that persist in shape during propagation. Their existence follows from variational and phase-plane arguments, and their qualitative properties were characterized in the influential work of Bona, Souganidis, and Strauss on solitary waves of KdV type [33]. Solitary-wave solutions of the generalized regularized long-wave equation are smooth and decay exponentially; indeed, Li and Bona proved that they are real analytic [34]. The generalized RLW equation was shown to possess both stable and unstable families of solitary waves, depending on the nonlinearity exponent [35], and stability criteria for model equations for long waves were formulated by Albert [36] and by Bona and Soyeur [37].

Stability Theory

The stability of solitary waves of the BBM equation is governed by the interplay between its conserved quantities. Albert, Bona, and Henry derived sufficient conditions for stability of solitary-wave solutions [9]. A major advance was the proof of asymptotic stability of BBM solitary waves by El Dika and Martel [10]. For the generalized BBM equation, the stability of periodic traveling waves was analyzed by Johnson, who derived necessary conditions for spectral instability [38]. Chen and co-workers established existence and stability results for solitary and periodic waves in perturbed generalized BBM equations [39]. Stability properties for fractional KdV and BBM equations were obtained by Angulo Pava and Natali [40], and periodic waves of the fractional BBM equation were shown to be stable in a range of parameter regimes [41]. Scattering of solutions and stability of solitary waves for the generalized BBM equation were analyzed by Esfahani [42], and orbital stability of solitary waves of the coupled Schrödinger–BBM system was established recently [43].

Instability and Related Phenomena

Not all solitary waves are stable. For large nonlinearity exponents, the generalized BBM equation admits unstable waves, and Jia and Wu recently proved a sharp instability result for the generalized BBM equation [44]. Natali, Pelinovsky, and Wang studied the instability of peaked traveling waves in a local model for shallow-water waves, a phenomenon with no analogue in the classical KdV theory [45]. In addition to solitary waves, the modified BBM equation supports breather-type solutions, which were constructed analytically by Adamashvili using a perturbation reduction to coupled nonlinear Schrödinger equations [46].

Exact and Explicit Solutions: Methods and Their Reliability

Classical Ansatz Methods

A large literature is devoted to constructing exact travelling-wave solutions of the BBM and related equations. Early approaches included the extended tanh method and its variants, expansions in terms of Weierstrass elliptic functions, and Jacobi elliptic function methods [15,47]. Nickel developed a unified approach to elliptic solutions of a generalized BBM equation and expressed the conditions for periodic and solitary waves in terms of the invariants of the Weierstrass function [15]. Mancas and co-workers derived Weierstrass-type travelling-wave solutions for the dissipative BBM equation using dynamical systems theory [16]. Periodic and soliton solutions of generalized BBM equations with variable coefficients were constructed with variational and exp-function techniques [48,49]. The generalized (G′/G)-expansion method was applied to the BBM and MBBM equations [50], and generalized (G′/G)- expansion approaches were used for the KP–BBM equation [51]. Wazwaz obtained exact solutions of compact and noncompact structures for the KP–BBM [52] and ZK–BBM [53] equations. Solitary-wave solutions for variants of the (3+1)-dimensional Wazwaz– BBM equations were derived with the Sardar subequation method [54], and the generalized exponential rational function method produced families of solutions for the generalized BBM equation [55].

More recent contributions continue this tradition. Analytical and numerical studies of the BBM equation were reported by Yokuş and co-workers [56], and soliton solutions of the BBM and modified RLW equations were constructed using classical algebraic methods [57]. New solitary-wave solutions of a generalized BBM equation with distributed delay were obtained in the framework of dynamical systems [58], and perturbed generalized BBM equations were analyzed by Zhu and co-workers [59]. Soliton solutions of the (3+1)-dimensional KdV–BBM equation were presented by Tariq and Seadawy [60], and soliton solutions of the conformable time-fractional generalized BBM equation were derived by Eslami and co-workers [61].

The pace of this activity has not slowed. In 2026 alone, exact solutions were reported for the modified BBM equation through an extended direct algebraic method [62], for the delayed KP–BBM equation [63], for the modified simplified Camassa–Holm and modified BBM equations [64], and for a (2+1)-dimensional BBM–Burgers equation describing wave interactions [65]. Symmetry-based studies of nonlocal structures in the modified BBM equation were reported by Bazghandi [66], and improved versions of the modified simplest equation method were applied to BBM-type systems [67].

The Validation Problem and Symbolic Computation

A serious methodological concern is the reliability of the reported solutions. Many solutions produced by heuristic ansatz methods are not checked by substitution into the governing equation, and independent verification has repeatedly revealed inconsistencies [14]. This has led to a growing emphasis on rigorous validation. In particular, a recent study revisited the travelling-wave reduction of the BBM equation within a symbolic computation framework, derived explicit solutions in terms of Weierstrass elliptic functions for both focusing and defocusing regimes, and verified every solution by symbolic substitution, thereby establishing certified benchmarks for nonlinear dispersive dynamics [14]. The same study demonstrated that certain previously published solitary-wave solutions, including one commonly attributed to the classical literature on the equation, fail to satisfy the equation, underscoring the need for systematic verification practices in this field [14]. These findings echo the broader consensus that reproducible and verifiable exact solutions are essential if analytical benchmarks are to guide numerical experimentation and physical interpretation.

Numerical Methods

Finite Difference Methods

The numerical study of the regularized long-wave equation began with the systematic investigation of Eilbeck and McGuire, who compared several finite difference schemes and assessed their accuracy for solitary waves [5]. Later finite difference treatments include a linearized implicit scheme for the one-dimensional RLW equation [68] and, for the generalized BBM–Burgers equation, three-level difference schemes with proven convergence [69]. Recent work has focused on structure-preserving discretizations: a linearized conservative compact difference scheme for the BBM–Burgers equation was proposed by Su and co-workers [70], a compact finite difference treatment of the BBM–Burgers equation was developed by Kaur and co-workers [71], and a fourth-order difference scheme for a generalized nonlinear time-fractional BBM-type equation was constructed recently [72]. For unbounded domains, local artificial boundary conditions for the BBM–Burgers equation were derived and analyzed by Deng and Li [73].

Galerkin and Finite Element Methods

Finite element methods are well suited to the BBM equation because its mixed spatial–temporal derivative leads to a variational structure that conserves the invariants of the equation. Fully discrete Galerkin methods for KdV-type equations were analyzed by Bona, Dougalis, and Karakashian [74], and the computational study of blow-up and decay for periodic solutions was pursued by the same group [75]. Conservative, high-order numerical schemes for generalized KdV equations were developed in a series of influential papers [76,77]. B-spline Galerkin methods for the RLW equation were introduced by Gardner, Gardner, and Dağ [78], and spline-based and splitting methods were studied by Jain and co-workers [79] and by Bhardwaj and Shankar [80]. For the BBM–Burgers equation with a dissipative term, an adaptive moving-mesh finite element method was developed by Lu and co-workers [81]. Unconditional superconvergence of a nonconforming BDF2 mixed finite element method for the BBM–Burgers equation was established by Xu and Shi [82], and a full discrete mixed finite element method for the BBM equation on a rectangular grid was recently analyzed [83].

Collocation, Meshless, and Spectral Methods

Collocation methods provide efficient alternatives for long-time simulations. Discontinuous Galerkin discretizations of the RLW equation were studied by Hozman [84], and solitary-wave solutions of the RLW equation were computed with a cubic B-spline collocation method by Irk [85]. Cubic B-spline methods were extended to the time-fractional BBM–Burgers equation by Kamran and co-workers [86]. Meshless approaches have gained attention for multi-dimensional problems: Oruç combined an extrapolated second-order BDF method with local radial point interpolation for multi-dimensional RLW equations [87], used Strang splitting with barycentric Floater–Hormann interpolants [88], and a superconvergent finite point method was applied to the generalized BBM– Burgers equation by Hou and Li [89]. A uniform hyperbolic polynomial B-spline collocation method for the RLW equation was reported by Ishwarya and Bhatia [90], and a hybrid shifted Gegenbauer integral-pseudospectral method for time-fractional BBM-type equations was introduced by Elgindy [91].

Analytical Approximation and Machine Learning

Semi-analytical approaches continue to complement fully discrete methods. He's variational method was applied to the BBM and Kawahara equations [92], and homotopy analysis methods were used for the RLW equation [93]. The BLUES function method was recently applied to the KdV–Burgers and BBM–Burgers equations [94]. A significant recent trend is the use of physics-informed neural networks (PINNs), which encode the governing equation into the training loss and thereby produce solutions that respect the model without a mesh. PINN-based solutions for the BBM–Burgers equation, in both forward and inverse formulations, were reported by Sabooni and co-workers [95]. These data-driven methods are promising for inverse problems and uncertainty quantification, though they currently require careful treatment of conservation properties and training data quality.

Structure Preservation and Error Control

Across all of these discretizations, two themes recur. The first is structure preservation: because the BBM equation conserves three independent quantities [8], numerical schemes that preserve discrete analogues of these invariants suppress the secular drift that plagues dissipative time-stepping [76,77]. The second is rigorous error control. Optimal a priori estimates in standard norms have been derived for Galerkin, mixed finite element, and difference methods [69,70,82], and superconvergence results are now available for conservative discretizations of the BBM–Burgers equation [82]. For exterior problems, artificial boundary conditions that are both accurate and efficiently implemented have been developed, allowing the BBM–Burgers equation to be simulated on unbounded domains without excessive computational cost [73]. These developments collectively make the BBM equation one of the best-understood nonlinear PDEs from the numerical standpoint, and a natural testbed for benchmarking new algorithms, including machine-learning-based ones [89,95].

Generalized, Modified, and Fractional Variants

Generalized BBM and BBM–Burgers Equations

The generalized BBM equation, obtained by replacing the quadratic nonlinearity with a power-law nonlinearity, retains the essential dispersive structure of the original model while exhibiting richer wave dynamics, including a threshold between stable and unstable solitary waves [35,44]. When a Burgers-type dissipative term is added, the resulting BBM–Burgers equation interpolates between hyperbolic conservation laws and dispersive regularization. Large-time asymptotic behavior of the BBM–Burgers equation was analyzed by Hayashi, Kaikina, and Naumkin [96], and global existence for the Cauchy problem of the generalized BBM–Burgers equation was established by Geng and Chen [97]. The large-time behavior of solutions converging to the self-similar Burgers diffusion wave was described by Fukuda and Ikeda [98]. Time-periodic solutions of the generalized BBM–Burgers equation with periodic forcing were constructed by Wang [99], and regularization–convergence results for the generalized BBM–Burgers equation were obtained by Kondo and Webler [100]. Nonlinear waves and modulation instability in the generalized Burger–BBM equation were recently investigated by Kumar [101].

Stabilization and control results for the BBM–Burgers equation have also advanced. Exponential input-to-state stability of the delayed BBM–Burgers equation was proved by Polwang and co-workers [102], and local null-controllability of the generalized BBM– Burgers equation with a fourth-order dissipative term was established by Bapu [103].

Modified and Higher-Dimensional Variants

The modified BBM (MBBM) equation with cubic nonlinearity supports breathers and solitons [46,62]. Higher-dimensional extensions, including the KP–BBM, ZK–BBM, and (3+1)-dimensional WBBM equations, generalize the BBM balance of nonlinearity and dispersion to two or three spatial dimensions and have been studied intensively for their exact and periodic wave solutions [52-54,63]. On metric networks, exact travelling waves of the BBM equation were constructed by Mugnolo and Rault, with conditions relating the wave to the geometry of the network [104].

Fractional Formulations

Fractional generalizations of the BBM equation have become a major theme, motivated by the desire to model memory effects and nonlocal dispersion. Fractional BBM–Burgers equations were solved with the homotopy analysis method [105] and with fractional homotopy analysis transform methods [106]. Exact solutions of the time-fractional BBM–Burger equation were obtained with a modified (G′/G)-expansion method [107] and with the residual power series method [108]. Conformable fractional derivatives were used to construct a variety of soliton solutions for the fractional Wazwaz–BBM equations [109]. The generalized fractional BBM equation, which models waves in nonlocally and nonlinearly elastic media, was analyzed both analytically and numerically [110]. A fractal variational approach to the Burgers–BBM equation was developed by Wang [111], and generalized Atangana fractional BBM– Burgers equations were solved with sub-equation and expansion methods [112]. Fractional models of the regularized long-wave equation were analyzed with q-homotopy analysis transform methods [113] and with fractional operators involving Mittag-Leffler kernels [114]. Very recent work includes analytical solutions of the time-fractional symmetric RLW equation [115], Lie symmetry, bifurcation, and stability analysis of the time-regularized long-wave equation [116], and chaos identification in the time-fractional symmetric RLW equation [117]. The long-time existence of a fractional KdV–BBM-type equation was recently established [118], and a fractional-order Kadomtsev–Petviashvili–BBM equation was analyzed for ocean engineering applications [13].

Dissipation, Control, Stochasticity, and Attractors

Damping and Control

The control-theoretic properties of the BBM equation have attracted sustained interest because the equation serves as a simplified model for the management of water waves and other dispersive systems. Micu proved that the linearized BBM equation is approximately controllable but not spectrally controllable [119]. Rosier introduced localized damping and boundary feedback mechanisms, proving global well-posedness and convergence of the damped BBM equation [120]. Very recently, controllability of the semilinear BBM dynamic equation on homogeneous time scales was established [121].

Stochastic Perturbations and Attractors

The long-time behavior of the BBM equation is governed by global and random attractors. Finite-dimensional behavior and strong attractors for the BBM equation were studied by Wang and Yang [122], and attractors for the generalized BBM equation were constructed by Çelebi, Kalantarov, and Polat [123]. Asymptotic smoothing and attractors for the generalized BBM equation in three dimensions were established by Stanislavova, Stefanov, and Wang [124]. For the stochastic BBM equation on unbounded domains, Wang proved the existence of compact random attractors [125]. Invariant manifolds of the stochastic BBM equation were analyzed by Kitaeva [126], and optimal control of the stochastic BBM equation was studied by Dhayal [127]. Weak uniform attractors for the three-dimensional BBM equations were recently characterized [128].

APPLICATIONS

Shallow-Water Waves, Undular Bores, and Tsunamis

The most direct physical applications of the BBM equation concern long waves in shallow water. Undular bores, which form when a tide or flood advances into a river, were modeled in the foundational computations of Peregrine [1] and in the analysis of shallow-water waves on a viscous fluid by Johnson [129]. Tsunamis ascending rivers develop into undular bores, a phenomenon documented by Tsuji and co-workers [11]. Modern treatments continue to employ regularized long-wave models for bore evolution [130]. The symmetric regularized long-wave equation describes weakly nonlinear ion-acoustic and space-charge waves in plasmas [12], and related models for shallow-water waves were given by Shivamoggi [20].

Elastic, Optical, and Ocean Engineering Settings

The generalized fractional BBM equation was derived for longitudinal waves in nonlocally and nonlinearly elastic media [110]. In optics and fluid mechanics, coupled BBM-type systems exhibit periodic, quasi-periodic, and chaotic wave regimes [131], and generalized variational principles for a BBM–Peregrine–Burgers equation were established recently [132]. Asymptotic analysis of the regularized long-wave equation in connection with the van der Pol oscillator was performed by Li [133]. For ocean engineering, fractional-order KP–BBM equations have been used to analyze nonlinear wave propagation [13]. These applications illustrate the versatility of the BBM framework across physical scales, and they explain why the equation remains a preferred benchmark model in applied mathematics.

TRENDS AND FUTURE DIRECTIONS

Several developments define the current trajectory of research on the BBM equation.

First, the reliability of exact solutions has become a pressing concern. The proliferation of unverified ansatz solutions, and the demonstration that some published solutions fail to satisfy the governing equation, has motivated a culture of validation based on symbolic substitution [14]. Future work should adopt standardized verification protocols, ideally by publishing the symbolic checks that certify each solution. Extending the validated Weierstrass-elliptic benchmarks of the BBM equation to generalized, fractional, and stochastic variants is a natural next step [14-16].

Second, machine learning is reshaping the numerical landscape. Physics-informed neural networks and related data-driven solvers have been applied to the BBM–Burgers equation [95], and their extension to fractional and stochastic formulations is an active frontier. Open questions include the conservation of invariants by learned solutions, rigorous error control, and the combination of symbolic computation with neural architectures for certified discovery of exact solutions.

Third, fractional and stochastic models are expanding in scope. Fractional BBM-type equations with Caputo, conformable, Mittag-Leffler, and fractal operators have been studied extensively [109,111,114], but rigorous existence, uniqueness, and stability results for these nonlocal models remain incomplete [118]. Stochastic BBM equations pose similar open problems regarding invariant measures and long-time behavior [125,127].

Fourth, higher-dimensional and networked geometries are under-explored. While KP–BBM, ZK–BBM, and WBBM extensions are well studied analytically [52,54], numerical methods that preserve the structure of these multi-dimensional models are comparatively rare. Travelling waves on networks [104] and wave interactions in (2+1)-dimensions [65] point toward applications in engineered systems, such as coastal protection and metamaterials.

Fifth, the interface between theory and data is growing. The comparison between model equations and measured wave fields, initiated by the classical fidelity studies [18], is being revisited with high-resolution observational data, and the BBM equation is increasingly used as a benchmark for surrogate modeling and uncertainty quantification.

Finally, we emphasize the value of critical synthesis. Reviews such as the present one serve to separate results that have been verified from those that remain conjectural, and to identify the methods, from total-positivity stability criteria [9] to energy-preserving discretizations [76], that have proven their worth over five decades of research.

CONCLUSION

The Benjamin–Bona–Mahony equation occupies a singular position in the theory of nonlinear dispersive waves. Its derivation as a regularized long-wave model addressed a genuine deficiency of the KdV equation, and its structural properties—finite set of conservation laws, well-posedness in classical spaces, and stable solitary waves—make it both physically relevant and mathematically tractable [3,8,10]. Fifty years of research have produced a mature theory of well-posedness, stability, and attractors, a rich collection of exact and numerical solution methods, and a broad range of applications. The most recent literature signals a transition toward rigorous solution validation, machine-learning-based computation, and fractional and stochastic extensions. We believe that the BBM equation will continue to serve as a canonical benchmark for nonlinear dispersive dynamics, and we hope this review provides a reliable map of the field for both newcomers and specialists.

DECLARATIONS

FUNDING

No funding.

COMPETING INTERESTS

The author declares no competing interests.

DATA AVAILABILITY

No new data were generated in support of this review. All references are cited with their digital object identifiers for direct verification.