Solving the biharmonic equation in irregular domains by the least squares collocation method

This paper addresses a new version of the least squares collocation (LSC) method of high order accuracy proposed and implemented for the numerical solution of the nonhomogeneous biharmonic equation. The differential problem is projected onto a polynomial space of the fourth and eighth degrees by the LSC method. The algorithm implemented is applied in irregular domains. The boundaries of these domains are given by analytical curves, in particular, by splines. The irregular domain is embedded in a rectangle covered by a regular grid with rectangular cells. In this paper we use the irregular cells (i-cells) which are cut off by the domain boundary from the rectangular cells of the initial regular grid. The idea of attaching elongated i-cells to the neighboring ones is used. A separate piece of the analytical solution is constructed in the combined cells. The collocation and matching points located outside the domain are used to approximate the differential equations in the boundary cells. These two approaches allow us to reduce essentially the conditionality of the corresponding system of linear algebraic equations. It is shown that the approximate solutions obtained by the LSC method converge with an increased order and coincide with the analytical solutions of the test problems with high accuracy in the case of the known solution. The numerical results are compared with those found by other authors who used a high order finite difference method (FDM). The nonhomogeneous biharmonic equation is used to model the stress-strain state (SSS) of isotropic thin irregular plates as an application.