摘要

In this paper, we propose a new implicit high order two level finite difference method of 0 (k(2) broken vertical bar kh(i) broken vertical bar h(l)(3)) for the solution of special type of 1D fourth order non-linear parabolic partial differential equation of the form u(xxxx) - 2u(xxt) + u(tt) = integral (x, t, u, u(x), u(xx) - u(t), u(xxx) - u(xt)), a < x < b, t >0 and a new three level finite difference method of O(k(2) + h(l)(3)) for the solution of another type of 1D fourth-order non-linear parabolic partial differential equation of the form A(x, t, u, u(xx))u(xxxx) + u(tt) = g(x, t, u, u(t), u(x), u(xx), u(xxx)), a < x < b, t > 0, both on a variable mesh, subject to suitable initial and Dirichlet boundary conditions, where k > 0 and h(l) > 0 are the mesh sizes in time and space coordinates, respectively. The third order variable mesh methods proposed are based on only three spatial grid points x(l),x(l +/- 1), meaning that no fictitious points are required for incorporating the boundary conditions. We discuss how our formulation is able to handle linear singular problem. The proposed two-level finite difference method is shown to be unconditionally stable for a model linear problem. We applied the proposed methods to find the numerical solution of more complex fourth-order non-linear equations like second-order Benjamin-Ono equation and the good Boussinesq equation. It is evident from the numerical experiments that the numerical results agree well with the exact solutions, hence demonstrating the efficiency of the methods derived.

  • 出版日期2016-1-15