A rational spectral method for the KdV equation on the half line (Q2390007)

The authors consider the Korteweg-de Vries (KdV) equation on the half line subject to initial condition and boundary conditions \[ U(0,t)=\lim_{x \rightarrow \infty} U(x,t)=\lim_{x \rightarrow \infty} U_x(x,t)=0. \] For the space discretization Legendre rational functions \[ \mathbb{R}_N^0:= \mathbb{R}_N\cap H_0^1(0,\infty), \,\, \mathbb{R}_N:=\{R_k(x):=\frac{\sqrt{2}}{x+1}L_k \Big(\frac{x-1}{x+1}\Big), \, k=0,\dots N\} \] are used as trial space, where \(L_k\) denotes the Legendre polynomial of degree \(k\). The test space is \(\mathbb{R}_{N-1}^0\) and the semi-discrete scheme has the form \[ (\partial_t u_N,v)-\frac{1}{2}(\Pi_N^C u_N^2,\partial_x v)-(\partial_x^2 u_N,\partial_x v)=(\Pi_N^C f,v), \] where \(\Pi_N^C u:= \omega^{1,0}P_N^C (\omega^{-1,0}u)\) with \(\omega^{1,0}:=2(1+x)^{-1}\), \(\omega^{-1,0}:=\frac{1}{2}(1+x)\) and \(P_N^C\) a certain interpolation operator in \(\mathbb{R}_{N}\) using Chebyshev-Gauss points. Also a fully-discrete scheme for the time discretization involving the Crank-Nicolson/leap-frog scheme with the linear part treated implicitely and the nonlinear part explicitly is considered. The schemes are shown to be convergent with order \(\mathcal{O}(N^{2-r})\) in space and second order in time. Two numerical examples are provided.