A numerical study of an ill-posed Boussinesq equation arising in water waves and nonlinear lattices: Filtering and regularization techniques

We consider an ill-posed Boussinesq equation which arises in shallow water waves and nonlinear lattices. This equation has growing and decaying modes in the linear as well as nonlinear regimes and its linearized growth rate a for short-waves of wavenumber k is given by sigma similar to k(2). Previou s numerical studies have addressed numerical difficulties and construction of approximate solutions for ill-posed problems with shortwave instability up to sigma similar to k, e.g. Kelvin-Helmholtz (sigma similar to k) and Ra yleigh-TayIor (sigma similar to root k) instabilities. These same issues ar e addressed and critically examined here for the present problem which has more severe short-wave instability In order to develop numerical techniques for constructing good approximate solutions of this equation, we use a fin ite difference scheme to investigate the effect of this short-wave instabil ity on the numerical accuracy of the exact solitary wave solution of this e quation. Computational evidence is presented which indicates that numerical accuracy of the solutions is lost very quickly due to severe growth of num erical errors, roundoff as well as truncation. We use both filtering and re gularization techniques to control growth of these errors and to provide be tter approximate solutions of this equation. In the filtering technique, nu merical experiments with three types of spectral filters of increasing orde r of regularity are performed. We examine the role of regularity of these f ilters on the accuracy of the numerical solutions. Numerical evidence is pr ovided which indicates that the regularity of a filter plays an important r ole in improving the accuracy of the solutions. In the regularization techn ique, the ill-posed equation is regularized by adding a higher order term t o the equation. Two types of higher order terms are discussed: (i) one that diminishes the growth rate of all modes below a cutoff wavenumber and sets the growth rate of all modes above it to zero; and (ii) the other one dimi nishes the growth rate of all modes and the growth rate asymptotically appr oaches to zero as the wavenumber approaches infinity. We have argued in fav or of the first type of regularization and numerical results using a finite difference scheme are presented. Numerical evidence is provided which sugg ests that regularization in combination with the most regular (C-2 here) sp ectral filter for small values of the regularization parameter can provide good approximate solutions of the ill-posed Boussinesq equation for longer time than possible otherwise. Some of the ideas presented here can possibly be utilized for solving other ill-posed problems with severe short-wave in stabilities and may have an important role to play in numerical studies of their solutions. (C) 1999 Elsevier Science Inc. All rights reserved.