The numerical solution of Cahn–Hilliard (CH) equation in one, two and three-dimensions via globally radial basis functions (GRBFs) and RBFs-differential quadrature (RBFs-DQ) methods

The present paper is devoted to the numerical solution of the Cahn–Hilliard (CH) equation in one, two and three-dimensions. We will apply two different meshless methods based on radial basis functions (RBFs). The first method is globally radial basis functions (GRBFs) and the second method is based on radial basis functions differential quadrature (RBFs-DQ) idea. In RBFs-DQ, the derivative value of function with respect to a point is directly approximated by a linear combination of all functional values in the global domain. The main aim of this method is the determination of weight coefficients. GRBFs replace the function approximation into the partial differential equation directly. Also, the coefficients matrix which arises from GRBFs is very ill-conditioned. The use of RBFs-DQ leads to the improvement of the ill-conditioning of interpolation matrix RBFs. The boundary conditions of the mentioned problem are Neumann. Thus, we use DQ method directly on the boundary conditions, which easily implements RBFs-DQ on the irregular points and regions. Here, we concentrate on Multiquadrics (MQ) as a radial function for approximating the solution of the mentioned equation. As we know this radial function depends on a constant parameter called shape parameter. The RBFs-DQ can be implemented in a parallel environment to reduce the computational time. Moreover, to obtain the error of two techniques with respect to the spatial domain, a predictor–corrector scheme will be applied. Finally, the numerical results show that the proposed methods are appropriate to solve the one, two and three-dimensional Cahn-Hilliard (CH) equations.