Article ID Journal Published Year Pages File Type
6933878 Journal of Computational Physics 2013 20 Pages PDF
Abstract
In this paper, a compact 5-point stencil for the discretisation of second-order partial differential equations (PDEs) in two space dimensions is proposed. We employ integrated radial basis functions in one dimension (1D-IRBFs) to construct the approximations for the dependent variable and its derivatives over the three nodes in each direction of the stencil. Certain nodal values of the second-order derivatives are incorporated into the approximations with the help of the integration constants. In the case of elliptic PDEs, one algebraic equation is formed at each interior node, and the obtained final system, of which each row has 5 non-zero entries, is solved iteratively using a Picard scheme. In the case of parabolic PDEs discretised with a Crank-Nicolson procedure, a set of three simultaneous algebraic equations is established at each interior node and the three equations are then combined to form two tridiagonal equations through the implicit elimination approach. Linear and non-linear test problems, including lid-driven cavity flow and natural convection between the outer square and the inner cylinder, are considered to verify the proposed stencil.
Related Topics
Physical Sciences and Engineering Computer Science Computer Science Applications
Authors
, ,