Fast Method of Particular Solutions for Solving Partial Differential Equations

Method of particular solutions (MPS) has been implemented in many science and engineering problems but obtaining the closed-form particular solutions, the selection of the good shape parameter for various radial basis functions (RBFs) and simulation of the large-scale problems are some of the challenges which need to overcome. In this dissertation, we have used several techniques to overcome such challenges. The closed-form particular solutions for the Matérn and Gaussian RBFs were not known yet. With the help of the symbolic computational tools, we have derived the closed-form particular solutions of the Matérn and Gaussian RBFs for the Laplace and biharmonic operators in 2D and 3D. These derived particular solutions play an important role in solving inhomogeneous problems using MPS and boundary methods such as boundary element methods or boundary meshless methods. In this dissertation, to select the good shape parameter, various existing variable shape parameter strategies and some well-known global optimization algorithms have also been applied. These good shape parameters provide high accurate solutions in many RBFs collocation methods. Fast method of particular solutions (FMPS) has been developed for the simulation of the large-scale problems. FMPS is based on the global version of the MPS. In this method, partial differential equations are discretized by the usual MPS and the determination of the unknown coefficients is accelerated using a fast technique. Numerical results confirm the efficiency of the proposed technique for the PDEs with a large number of computational points in both two and three dimensions. We have also solved the time fractional diffusion equations by using MPS and FMPS

Hybrid Chebyshev Polynomial Scheme for the Numerical Solution of Partial Differential Equations

In the numerical solution of partial differential equations (PDEs), it is common to find situations where the best choice is to use more than one method to arrive at an accurate solution. In this dissertation, hybrid Chebyshev polynomial scheme (HCPS) is proposed which is applied in two-step approach and one-step approach. In the two-step approach, first, Chebyshev polynomials are used to approximate a particular solution of a PDE. Chebyshev nodes which are the roots of Chebyshev polynomials are used in the polynomial interpolation due to its spectral convergence. Then, the resulting homogeneous equation is solved by boundary type methods including the method of fundamental solution (MFS) and the equilibrated collocation Trefftz method. However, this scheme can be applied to solve PDEs with constant coefficients only. So, for solving a wide variety of PDEs, one-step hybrid Chebyshev polynomial scheme is proposed. This approach combines two matrix systems of two-step approach into a single matrix system. The solution is approximated by the sum of particular solution and homogeneous solution. The Laplacian or biharmonic operator is kept on the left hand side and all the other terms are moved to the right hand side and treated as the forcing term. Various boundary value problems governed by the Poisson equation in two and three dimensions are considered for the numerical experiments. HCPS is also applied to solve an inhomogeneous Cauchy-Navier equations of elasticity in two dimensions. Numerical results show that HCPS is direct, easy to implement, and highly accurate

Numerical Solution of Partial Differential Equations Using Polynomial Particular Solutions

Polynomial particular solutions have been obtained for certain types of partial differential operators without convection terms. In this dissertation, a closed-form particular solution for more general partial differential operators with constant coefficients has been derived for polynomial basis functions. The newly derived particular solutions are further coupled with the method of particular solutions (MPS) for numerically solving a large class of elliptic partial differential equations. In contrast to the use of Chebyshev polynomial basis functions, the proposed approach is more flexible in selecting the collocation points inside the domain. Polynomial basis functions are well-known for yielding ill-conditioned systems when their order becomes large. The multiple scale technique is applied to circumvent the difficulty of ill-conditioning. The derived polynomial particular solutions are also applied in the localized method of particular solutions to solve large-scale problems. Many numerical experiments have been performed to show the effectiveness of the particular solutions on this algorithm. As another part of the dissertation, a modified method of particular solutions (MPS) has been used for solving nonlinear Poisson-type problems defined on different geometries. Polyharmonic splines are used as the basis functions so that no shape parameter is needed in the solution process. The MPS is also applied to compute the sizes of critical domains of different shapes for a quenching problem. These sizes are compared with the sizes of critical domains obtained from some other numerical methods. Numerical examples are presented to show the efficiency and accuracy of the method