Premium
Quadratic Störmer‐type methods for the solution of the Boussinesq equation by the method of lines
Author(s) -
Tsitouras Ch.,
Famelis I. Th.
Publication year - 2008
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.20317
Subject(s) - ode , mathematics , discretization , ordinary differential equation , method of lines , partial differential equation , quadratic equation , type (biology) , mathematical analysis , differential equation , differential algebraic equation , geometry , ecology , biology
We study the numerical treatment of Boussinesq PDE equation using the method of lines. For the space discretization, we choose either classical finite differences or Fourier pseudospectral methods. Both cases result in a system of second‐order ordinary differential equations (ODEs) that is quadratic. In order to take advantage of this special feature, we choose to solve the ODE system using a new type of hybrid Numerov method specially constructed for such problems. Other efficient ODE solvers taken from the literature are used to solve the system of ODEs as well. By taking all the combinations of space discretization methods and ODE solvers, we discuss the stability and accuracy features revealed from the numerical tests. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008