Recently, Ramos and Garcia Lopez (2024) and Gao (2023a, 2023b) have conducted several comprehensive shallow-water studies in International Journal of Numerical Methods for Heat and Fluid Flow, that is, blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in the shallow water (Ramos and Garcia Lopez, 2024), novel hetero-Bäcklund transformations and similarity reductions for a (2 + 1)-dimensional generalized dispersive long-wave system that models certain dispersive and nonlinear long gravity waves on the oceanic shallow water (Gao, 2023a), and novel similarity reductions for a (2 + 1)-dimensional generalized modified dispersive water-wave system for the nonlinear and dispersive long gravity waves in the shallow water of uniform depth (Gao, 2023b).
Those studies (Gao, 2023a, 2023b; Ramos and Garcia Lopez, 2024) have sparked renewed interest in shallow-water research. Consequently, this Letter is dedicated to the Whitham-Broer-Kaup system for the dispersive long waves in the shallow oceanic water (Ahmad et al., 2015; Fan and Bao, 2024; Sabawi and Hamad, 2024; Fan and Bao, 2022; Gao et al., 2020; Cao et al., 2020; Arshad et al., 2017; Ren and Lin, 2018; Kuo, 2017; Fei et al., 2017; Zhou and Lu, 2017; Xu et al., 2017; Lin et al., 2011a, 2011b; Bhrawy et al., 2014; Shan et al., 2012; Wang et al., 2010, 2009; Zhang et al., 2008) (and references therein), i.e.:
where x and t are the scaled space and time coordinates, the subscripts represent the partial derivatives, p(x, t) and q(x, t) are the real functions denoting the horizontal velocity field and height of the deviation from the equilibrium position, respectively, α is the real dispersion coefficient and β is the real diffusion coefficient.
For System (1), Bäcklund transformation, soliton and similar wave solutions have been studied (Fan and Bao, 2024); quartic B-spline collocation method has been used to obtain the numerical solutions (Sabawi and Hamad, 2024); Painlevé integrability, bilinear forms and Bäcklund transformations have been offered (Gao et al., 2020; Zhang et al., 2008; Fan and Bao, 2022); two different groups of variational principles have been constructed (Cao et al., 2020); travelling wave solutions in the form of solitons, dark solitons, bell and anti-bell periodic waves have been obtained (Ahmad et al., 2015; Arshad et al., 2017); nonlocal symmetry, consistent tanh-expansion solvability and power-series solutions have been studied (Ren and Lin, 2018; Zhou and Lu, 2017); modified simplest equation method has been used to derive some soliton solutions (Kuo, 2017); residual symmetries and interaction solutions have been investigated (Fei et al., 2017); double Wronskian solutions and soliton interactions have been studied via the Wronskian technique (Xu et al., 2017; Lin et al., 2011a, 2011b; Wang et al., 2009); solitons, cnoidal waves and snoidal waves have been presented (Bhrawy et al., 2014); as well as Darboux transformation and some soliton solutions have been investigated (Shan et al., 2012; Wang et al., 2010). Other contributions on System (1) have appeared as the references in the aforementioned papers.
However, to the best of our knowledge, for System (1), generalized Darboux transformation (GDT), which can be used to derive the rogue wave, multi-pole and mixed solutions, has not yet been reported. Using symbolic computation (Gao, 2023a, 2024), we aim to construct a GDT for System (1). With respect to the horizontal velocity field and the height of the deviation from the equilibrium position, the GDT supports the derivation of the multi-pole solitons, multi-pole breathers, higher-order rogue waves and mixed waves. We hope this study contributes to a deeper understanding of the physical dynamics of the nonlinear waves in shallow water environments.
Via the transformations:
we transform System (1) into:
Wang et al. (2010) has offered a Lax pair for System (3):
where Φ = (ϕ, ψ)T is the eigenfunction, λ is a real spectral parameter and ϕ1 and ϕ2 are the complex functions of x and t. The compatibility condition Φxt = Φtx leads to System (3).
We consider an n-fold gauge transformation that transforms Lax Pair (4) into:
where [N] represents the N-fold solutions/matrices, Φ[N] = (ϕ[N], ψ[N])T, Δ[N] is the 2 × 2 gauge transformation matrix and N is a positive integer. U[N] and V[N] own the same forms as U and V except that the initial potentials P and Q are replaced with their N-fold versions P[N] and Q[N]. Combining Gauge Transformation (5) with Lax Pair (4), we derive:
Motivated by Wang et al. (2010), the n-fold gauge matrix Δ[n] can be set as:
where j = 0, 1, 2,…, N – 1, aj’s, bj’s, cj’s and dj’s are 4 N functions of x and t, and ω2 = 1 + cN–1. To determine 4N unknown functions aj’s, bj’s, cj’s and dj’s, we suppose that λ1, λ2, …, λ2n are the 2n roots of Det(Δ[N]), i.e.:
where Det denotes the determinant, k = 1, 2, …, n, n ≤ N, mk is a nonnegative integer and . We assume that Φ2k–1 = (ϕ2k–1, ψ2k–1)T and Φ2k = (ϕ2k, ψ2k)T respectively are the eigenfunction of Lax Pair (4) at λ = λ2k–1, and λ = λ2k, and the superscript T denotes the transpose of the matrix.
According to and , we give the following 4n equations:
Nonetheless, the 4n equations presented in equations (9) are inadequate for uniquely determining the 4N unknown functions Aj’s, Bj’s, Cj’s and Dj’s. We consider the Taylor expansions and at ϵ = 0:
where ϵ is a small parameter, ι = 0, 1, …, ζ, ζ = 0, 1, …, mk and .
We truncate and to the mk-th order and obtain:
i.e.:
Combining Transformations (2), equations (6), N-Fold DT Matrix (7) and equations (11), we calculate out the Nth-order solutions for System (1) as:
where Yk is a 2(mk + 1) × 2N matrix, 1 ≤ ς ≤ 2(mk + 1), 1 ≤ τ ≤ 2N, (Yk)ς,τ denotes the element in the ςth row and τth column of Yk, Yb can be obtained via the modified Y with its (N + 1) th column replaced by the vector (υ)2N×1, Yc can be obtained via the modified Y with its 1st column replaced by the vector (µ)2N×1.
In conclusion, N-Fold DT Matrix (7) and The Nth-Order Solutions (12) together establish an N-fold GDT for System (1). In The Nth-Order Solutions (12), the variable n indicates the number of sets of spectral parameters utilized, the variable mk represents that the kth set of spectral parameters is iterated mk times (a total of mk + 1 times) and the value of N represents the total order of The Nth-Order Solutions (12). Notably, when n = N, i.e. mk = 0, each of the n spectral parameters is iterated only once. In this case, the resulting GDT simplifies to the N-fold DT proposed by Shan et al. (2012) and Wang et al. (2010).
In relation to the horizontal velocity field and height of the deviation from the equilibrium position, the previously obtained GDT allows for multiple iterations of the spectral parameters, leading to the generation of the multi-pole solitons, multi-pole breathers, higher-order rogue waves and mixed waves. This study might assist in exploring the complex and dynamic natural mechanisms underlying real-world shallow water waves.
Statements and Declarations
Funding: This work has been supported by the National Natural Science Foundation of China under Grant No.11772017, and the National Scholarship for Doctoral Students under Grant BSY202306241.
Conflict of interest: The authors declare that they have no conflict of interest.
Data availability: Data sharing not applicable to this article, as no datasets were generated or analysed during the current study.
