--- title: Good Boussinesq Equation url: https://www.emergentmind.com/topics/good-boussinesq-equation type: topic --- # Good Boussinesq Equation The good Boussinesq equation denotes the well-posed sign choice within the Boussinesq family of bidirectional long-wave models. In the four-parameter family \[ u_{TT}-c\,u_{XX}-a_4\,u_{XXXX}-a_2\,(u^2)_{XX}=0, \] the choice \(c=-1\) is called the “good” Boussinesq equation, whereas \(c=+1\) gives the ill-posed classical Boussinesq equation [1504.03206]. In the literature, this designation appears in several shifted and rescaled normalizations, including \[ u_{tt}-u_{xx}+(u^2)_{xx}+u_{xxxx}=0, \qquad u_{tt}+\tfrac43\,(u^2)_{xx}+\tfrac13\,u_{xxxx}=0, \] as well as shifted forms obtained by removing linear terms through \(u=w+\tfrac12\) [2003.02777]. The equation belongs to the Korteweg–de Vries kind of equations, models long-wave propagation, and occupies a central position at the intersection of dispersive PDE, Hamiltonian dynamics, integrable systems, inverse scattering, and structure-preserving numerics. ## 1. Formulations, normalizations, and variable conventions A standard periodic form on the one-dimensional torus \(\mathbb T=\mathbb R/(2\pi\mathbb Z)\) is \[ u_{tt}+u_{xxxx}-u_{xx}+\partial_{xx}(u^2)=0, \qquad u(0,x)=u_0(x),\quad u_t(0,x)=u_1(x), \] with data in \(H^s(\mathbb T)\times H^{s-2}(\mathbb T)\) [1201.1942]. On the line, inverse-scattering treatments often use \[ u_{tt}-u_{xx}+(u^2)_{xx}+u_{xxxx}=0, \qquad x\in\mathbb R,\ t>0, \] or an equivalent first-order system obtained by introducing \[ v(x,t)=\int_{-\infty}^x u_t(x',t)\,dx', \qquad u_t=v_x \] [2003.04789]. A shifted periodic Hamiltonian formulation begins with \[ w_{tt}=-w_{xxxx}+w_{xx}+(w^2)_{xx}, \] then introduces \[ u(x,t)=w(x,t)+\tfrac12, \] so that \[ u_{tt}=-u_{xxxx}+(u^2)_{xx}. \] This form is convenient for Hamiltonian and geometric numerical analysis [1807.05182]. Some expositions also consider a reduced “good” Boussinesq equation obtained after dropping the fourth-order term, \[ u_{tt}-u_{xx}-(u^2)_{xx}=0, \] as a simplified well-posed model associated with the same sign choice \(c=-1\) [1504.03206]. A persistent source of terminological ambiguity is therefore not the sign structure itself but the coexistence of several equivalent or reduced normalizations across the analytical and numerical literature. ## 2. Derivation and physical interpretations One derivation starts from the two-dimensional irrotational Euler equations for a fluid layer of mean depth \(h_0\) under gravity, in the long-wave, small-amplitude regime. Introducing the nondimensional variables \[ X=\delta x,\qquad T=\delta t,\qquad \delta\ll1, \] and expanding the velocity potential and free-surface elevation in powers of \(\delta\), one obtains, after eliminating the potential and retaining terms up to \(O(\delta^2)\), the Boussinesq family \[ u_{TT}-c\,u_{XX}-a_4\,u_{XXXX}-a_2\,(u^2)_{XX}=0 \] [1504.03206]. In this setting, \(u\) represents the leading-order surface elevation. The same source connects the equation to the motion of long waves in two dimensions under gravitation and to nonlinear wave propagation in rods or waveguides interacting with an external environment. In such settings one must account for weak nonlinearity \((u^2)_{xx}\), dispersion due to elasticity or fluid filling, and possible energy leakage at the boundary. For thin rods immersed in a viscous medium, a Hamilton principle followed by asymptotic expansion leads to a third-order corrected Boussinesq equation with two dispersive terms, while the reduced good Boussinesq model captures the leading long-wave balance and predicts solitary pulses and periodic wave trains [1504.03206]. In the shifted periodic formulation, the equation admits a Hamiltonian first-order system \[ u_t=v_x,\qquad v_t=-u_{xxx}+(u^2)_x, \] with Poisson operator \[ J=\partial_x \begin{pmatrix} 0&1\\ -1&0 \end{pmatrix} \] and Hamiltonian functional \[ H[u,v]=\frac12\int_a^b\Bigl(v^2+u_x^2+\tfrac23\,u^3\Bigr)\,dx. \] The variational derivatives are \[ \delta_u H=u^2-u_{xx},\qquad \delta_v H=v, \] and periodic integration by parts yields conservation of \(H\). A second quadratic invariant, \[ M[u,v]=\int_a^b u\,v\,dx, \] is also conserved [1807.05182]. These conservation laws are fundamental in both analysis and geometric discretization. ## 3. Well-posedness, low-regularity theory, and smoothing For the periodic equation on \(\mathbb T\), \[ u_{tt}+u_{xxxx}-u_{xx}+\partial_{xx}(u^2)=0, \] local well-posedness holds in \(H^s(\mathbb T)\times H^{s-2}(\mathbb T)\) for \[ s>-\tfrac38, \] under the mean-zero assumption \(\int_{\mathbb T}u_0=\int_{\mathbb T}u_1=0\). The proof uses a normal form approach that explicitly extracts the rougher part of the solution. The nonlinear remainder then lies in \[ C([0,T];H^{s+a}(\mathbb T)), \qquad 0\le a<\min\{2s+1,\tfrac12\}, \] so the nonlinear part gains up to \(\min(2s+1,1/2)\) derivatives [1201.1942]. For the reduced equation \[ u_{tt}-u_{xx}-(u^2)_{xx}=0, \] an energy computation yields the formally conserved quantity \[ E[u](t)=\int_{\mathbb R}\Bigl(\tfrac12u_t^2+\tfrac12u_x^2+\tfrac13u^3\Bigr)\,dx. \] Because the quadratic form \(\int (u_t^2+u_x^2)\,dx\) is positive definite, one obtains an a priori bound of \(\|u\|_{H^1_x}+\|u_t\|_{L^2_x}\). In the linearization about zero, \(u_{tt}-u_{xx}=0\) is the standard wave operator, and a fixed-point or energy-method argument in \[ \mathcal C([0,T];H^s(\mathbb R))\cap \mathcal C^1([0,T];H^{s-1}(\mathbb R)) \] gives local well-posedness for \(s>3/2\), while conservation of \(E\) together with \(H^1\hookrightarrow L^\infty\) gives global existence. By contrast, when \(c=+1\), the energy is indefinite and the problem is ill-posed [1504.03206]. On the half-line, the initial-boundary-value problem \[ u_{tt}-u_{xx}+u_{xxxx}+(u^2)_{xx}=0, \qquad x>0,\ t>0, \] with data \[ u(x,0)=f(x),\quad u_t(x,0)=g(x),\quad u(0,t)=h_1(t),\quad u_x(0,t)=h_2(t), \] is locally well-posed in \(H^s(\mathbb R^+)\) for \[ s\in(-4,5),\qquad s\neq-\tfrac12,\ \tfrac34. \] The boundary data belong to \[ (h_1,h_2)\in H^{(2s+1)/4}(\mathbb R^+)\times H^{(2s-1)/4}(\mathbb R^+), \] with compatibility conditions when \(s>-\tfrac12\) and \(s>\tfrac34\). The proof uses Bourgain-type restricted-norm spaces \(X^{s,b}\), linear estimates for the free and boundary-forcing operators, and a bilinear estimate for \[ M=(-\partial_{xx}+\partial_{xxxx})^{-1/2}. \] Moreover, the nonlinear part of the solution is smoother than the data: \[ u_{nl}\in C([0,T];H^{s+a}(\mathbb R_x^+)), \qquad 0-4\) and the half-derivative smoothing are sharp [1611.09255]. ## 4. Traveling waves, bilinearization, and integrable structures A traveling-wave reduction for the reduced good Boussinesq equation takes \[ u(x,t)=h(z),\qquad z=x-\mu t. \] Then \[ (\mu^2-1)\,h''-(h^2)''=0. \] After integration, the analysis in the cited paper yields three principal families of exact solutions in Jacobi elliptic functions: a periodic \(\mathrm{sn}\)-wave, a periodic \(\mathrm{cn}^3\)-wave, and a \(\mathrm{dn}\)-family whose \(m\to1\) limit is the solitary \(\sech^2\) wave. The latter can be written as \[ u(x,t)=\sech^2(x-t) \] in a particular normalization, or more generally \[ u(x,t)=\tfrac12(\mu^2-1)\,\sech^2\!\Bigl(\tfrac12\sqrt{\mu^2-1}\,(x-\mu t)\Bigr), \qquad \mu>1 \] [1504.03206]. For the normalized fourth-order form \[ u_{tt}-u_{xx}+(u^2)_{xx}+\tfrac13\,u_{xxxx}=0, \] the binary Bell polynomial method yields a bilinear representation. Introducing \[ u=q_{xx},\qquad q=2\ln F, \] one obtains \[ (D_t^2-D_x^2+\tfrac13 D_x^4)\,F\!\cdot\!F=0. \] The same framework produces the \(n\)-soliton \(\tau\)-function, a bilinear Bäcklund transformation, a Lax pair, an infinite hierarchy of conservation laws, and Wronskian determinant solutions. The authors state that these constructions fully demonstrate the complete integrability of the good Boussinesq equation in this normalization [2305.06853]. A distinct integrable-system development establishes a Miura-type transformation from the second-order equation \[ i\,q_t-\tfrac1{\sqrt3}\,q_{xx}-2\sqrt3\,\bar q\,\bar q_x=0 \] to the fourth-order good Boussinesq equation \[ u_{tt}-u_{xx}+(u^2)_{xx}+u_{xxxx}=0. \] The correspondence is derived by comparing a regular and a singular \(3\times3\) Riemann–Hilbert problem, so the Miura map exists simultaneously at the PDE, Lax-pair, and RH levels [2301.07620]. ## 5. Lax pairs, Riemann–Hilbert formulations, and inverse spectral theory For the rescaled good Boussinesq equation \[ u_{tt}+\tfrac43\,(u^2)_{xx}+\tfrac13\,u_{xxxx}=0, \] a line-based inverse scattering transform is formulated through the Lax pair \[ \mathsf L=\partial_x^3+2u\,\partial_x+u_x+v, \qquad \mathsf A=\partial_x^2+\tfrac43\,u, \] or equivalently through a \(3\times3\) matrix system with spectral parameter \(k\), cubic-root symmetry \(\omega=e^{2\pi i/3}\), and contour \[ \Gamma=\mathbb R\cup\omega\mathbb R\cup\omega^2\mathbb R. \] Under no-soliton and nondegeneracy assumptions, the RH problem is determined by two reflection coefficients, \[ r_1(k)=\frac{s_{12}(k)}{s_{11}(k)},\qquad k>0, \qquad r_2(k)=\frac{s^A_{12}(k)}{s^A_{11}(k)},\qquad k<0, \] and the solution is reconstructed by \[ u(x,t)=-\tfrac32\,\partial_x\lim_{k\to\infty}k\bigl[M_{33}(x,t,k)-1\bigr], \qquad v(x,t)=-\tfrac32\,\partial_t\lim_{k\to\infty}k\bigl[M_{33}(x,t,k)-1\bigr]. \] The formulation is specifically designed to support Deift–Zhou steepest descent analysis [2003.02777]. On the half-line, assuming existence of a sufficiently smooth rapidly decaying solution, the good Boussinesq equation can likewise be recovered from a \(3\times3\) Riemann–Hilbert problem depending only on initial and boundary values. In that setting the jump contour consists of twelve half-lines, the unknown matrix is sectionally analytic in twelve sectors, and the reconstruction formula again uses the large-\(k\) behavior of the \((3,3)\)-entry. The analysis is carried out under “no-soliton” and “generic at \(k=0\)” hypotheses [2603.11951]. In periodic inverse-spectral theory, the good Boussinesq equation is written in Lax form with the third-order operator \[ L\,y=(y''+p\,y)'+p\,y'+q\,y \] on \([0,2]\) under the three-point Dirichlet conditions \[ y(0)=y(1)=y(2)=0. \] Its discrete Dirichlet spectrum constitutes the auxiliary spectrum of the equation, in analogy with the Hill operator for periodic KdV. For \((p,q)\in H^1(\mathbb T)\oplus H^1(\mathbb T)\) near zero, each neighborhood \(D_n\) contains exactly one simple real eigenvalue \(\lambda_n\), and the corresponding norming constants \(h_n\) admit uniform high-energy asymptotics. The spectral data determine \((p,q)\) uniquely near the zero potential through a nonlinear Riemann–Hilbert or nonlinear Gelʹfand–Levitan–Marchenko scheme adapted to third-order operators. This work is described as the first in a series devoted to solving the inverse problem for the Boussinesq equation [2409.10988]. ## 6. Long-time asymptotics and dispersive regimes For Schwartz-class initial data on the line, under no-soliton and generic-pole-at-zero assumptions, the inverse-scattering solution can be analyzed by Deift–Zhou steepest descent. In the region \[ \zeta=\frac{x}{t}>\zeta_0,\qquad k_0(\zeta)=\frac{\zeta}{2}, \] the solution exhibits modulated oscillatory asymptotics of size \(t^{-1/2}\). Writing \[ \nu(\zeta)=-\frac1{2\pi}\ln\bigl(1-|r_1(k_0)|^2\bigr)\ge0, \] one obtains an explicit leading term with amplitude proportional to \[ \frac{k_0\sqrt{\nu}}{\sqrt{t}} \] and phase containing \(-\sqrt3\,k_0^2 t\), \(\nu\ln t\), \(\arg r_1(k_0)\), and an integral correction, with error \[ O\!\bigl(\tfrac{\ln t}{t}\bigr) \] uniformly for \(\zeta\) in compact subsets of \((\zeta_0,\infty)\) [2003.04789]. A distinct asymptotic regime arises in the Painlevé region \[ |x|=O(\sqrt t). \] There the associated \(3\times3\) RH problem deforms to a Painlevé IV model problem, and the modified Boussinesq fields \(p(x,t)\) and \(q(x,t)\) are described by the Clarkson–McLeod solution \(P_{\rm IV}(y)\) of Painlevé IV with parameters \[ \alpha=-\tfrac16,\qquad \beta=-\tfrac23, \qquad y=-\frac{\sqrt3}{2\sqrt t}\,x. \] Via the Miura transformation, the good Boussinesq field \(u(x,t)\) has an explicit leading term of order \(t^{-1}\), with uniform remainder \[ O(t^{-3/2}) \] for \(|x|/\sqrt t\le c_1\). The same work states that the theoretical asymptotic solutions were validated against direct numerical simulations [2511.13382]. Taken together, these results separate at least two analytically distinct dispersive regimes: an oscillatory region governed by stationary-phase analysis around nonzero critical points, and a near-origin Painlevé region governed by the small-\(k\) structure of the \(3\times3\) RH problem. ## 7. Numerical analysis and structure-preserving discretization A fully discrete Fourier pseudospectral method with second-order temporal accuracy has been analyzed for the periodic good Boussinesq equation \[ u_{tt}=-u_{xxxx}+u_{xx}+(u^p)_{xx}, \qquad p\ge2. \] Using a specially designed second-order time-stepping with an auxiliary variable \(\psi\approx u_t\), the analysis proves unconditional nonlinear stability, with no restriction of the form \(\Delta t\le C h^2\), and establishes the convergence estimate \[ \|u_{,h}^n-u^*(\cdot,t^n)\|_{H^2} + \|\psi_{,h}^n-u_t^*(\cdot,t^n)\|_2 \le C(\Delta t^2+h^m) \] under the stated regularity assumptions [1401.6327]. A geometric alternative combines Fourier spectral semi-discretization in space with Hamiltonian Boundary Value Methods in time. In the shifted Hamiltonian formulation, HBVM\((k,s)\) has order \(2s\), is symmetric for all \(k\ge s\), and exactly conserves polynomial Hamiltonians of degree \(\deg H\le\lfloor 2k/s\rfloor\). For the cubic discrete Hamiltonian \(H_N\), exact energy conservation requires \[ k\ge\tfrac32 s. \] The SHBVM variant, with large \(s\), acts as a spectral integrator in time and, in the reported tests, preserved both the Hamiltonian and momentum to round-off [1807.05182]. Exponential-type integrators have also been developed for the periodic equation. One first-order scheme proves convergence in \(H^r\), \[ \|u(t_n)-u^n\|_{H^r}\le C\tau^\gamma, \qquad r>\tfrac12,\quad 0<\gamma\le1, \] and in particular first-order convergence when the exact solution remains in \(H^{r+1}\). A second scheme proves \[ \|u(t_n)-u^n\|_{H^r}\le C\tau^{1+\gamma} \] when the exact solution belongs to \(H^{r+1+2\gamma}\), yielding second order in \(H^r\) for \(u\in H^{r+3}\), and \[ \|z(t_n)-z^n\|_{L^2}+\|z_t(t_n)-z_t^n\|_{H^{-2}}\le C\tau^2 \] for \(H^3\)-solutions. The same paper emphasizes that the regularity requirements are lower than those of classical exponential integrators and states an \(O(N\log N)\) cost per step [1902.07478]. A Deuflhard-type exponential integrator Fourier pseudospectral method rewrites the equation as \[ z_{tt}+Lz+N(z)=0,\qquad L=-\partial_x^2+\partial_x^4,\qquad N(z)=-\partial_{xx}(f(z)), \] and combines a trapezoidal-type exponential integrator in time with FFT-based mode updates in space. For sufficiently smooth solutions, the method satisfies \[ \|z(\cdot,t_n)-z_I^n\|_{H^m(\Omega)} + \|\partial_t z(\cdot,t_n)-v_I^n\|_{H^{m-2}(\Omega)} \le C(\tau^2+h^\sigma), \] so it is quadratically convergent in time and spectrally accurate in space, without any CFL-type condition [1910.03267]. More recently, a filtered Lie–Trotter splitting scheme has been analyzed for low-regularity data on the periodic problem \[ z_{tt}+z_{xxxx}-z_{xx}-(z^2)_{xx}=0. \] After rewriting the PDE as a first-order Schrödinger-type system and introducing the spectral cut-off \[ \Pi_\tau f(x)=\sum_{|k|\le\tau^{-1/2}}\widehat f_k e^{ikx}, \] the analysis in discrete Bourgain spaces \(X_\tau^{s,b}\) yields, for \(s\in(0,2]\), \[ \|z(t_n)-z^n\|_{L^2} + \|z_t(t_n)-z_t^n\|_{H^{-2}} \lesssim \tau^{s/2}. \] The paper states that this extends low-regularity convergence down to \(s>0\), overcoming the classical \(s>1/2\) bilinear barrier in smooth Sobolev spaces [2402.11266]. The numerical literature therefore reflects the same structural themes present in the analysis: dispersive stiffness, low-regularity sensitivity, and the special role of Hamiltonian and integrable structure in designing accurate long-time algorithms. Source: https://www.emergentmind.com/topics/good-boussinesq-equation