High Accuracy Compact Operator Methods for Two-Dimensional Fourth Order Nonlinear Parabolic Partial Differential Equations

Abstract

AbstractIn this study, we develop and implement numerical schemes to solve classes of two-dimensional fourth-order partial differential equations. These methods are fourth-order accurate in space and second-order accurate in time and require only nine spatial grid points of a single compact cell. The proposed discretizations allow the use of Dirichlet boundary conditions only without the need to discretize the derivative boundary conditions and thus avoids the use of ghost points. No transformation or linearization technique is used to handle nonlinearity and the obtained block tri-diagonal nonlinear system has been solved by Newton’s block iteration method. It is discussed how our formulation is able to tackle linear singular problems and it is ensured that the methods retain their orders and accuracy everywhere in the solution region. The proposed two-level method is shown to be unconditionally stable for a class of two-dimensional fourth-order linear parabolic equation. We also discuss the alternating direction implicit (ADI) method for solving two-dimensional fourth-order linear parabolic equation. The proposed difference methods has been successfully tested on the two-dimensional vibration problem, Boussinesq equation, extended Fisher–Kolmogorov equation and Kuramoto–Sivashinsky equation. Numerical results demonstrate that the schemes are highly accurate in solving a large class of physical problems.