#### Abstract

The interpolating moving least-squares (IMLS) method is discussed in detail, and a simpler formula of the shape function of the IMLS method is obtained. Then, based on the IMLS method and the Galerkin weak form, an interpolating element-free Galerkin (IEFG) method for two-point boundary value problems is presented. The IEFG method has high computing speed and precision. Then error analysis of the IEFG method for two-point boundary value problems is presented. The convergence rates of the numerical solution and its derivatives of the IEFG method are presented. The theories show that, if the original solution is sufficiently smooth and the order of the basis functions is big enough, the solution of the IEFG method and its derivatives are convergent to the exact solutions in terms of the maximum radius of the domains of influence of nodes. For the purpose of demonstration, two selected numerical examples are given to confirm the theories.

#### 1. Introduction

Conventional computational methods, such as the finite element method (FEM) and the boundary element method (BEM), cannot be applied well to some engineering problems. For the extremely large deformation and crack growth problems, the remeshing technique must be used. Meshless methods have been developed in recent years. The most important common feature of meshless methods is that the trial function is constructed from a set of nodes with no meshing at all. Then some complex problems, such as the large deformation and crack growth problems, can be simulated with the method without the remeshing techniques [1, 2].

Many kinds of meshless methods have been proposed, such as element-free Galerkin (EFG) method , complex variable meshless method , mesh-free reproducing kernel particle Ritz method , radial basis function (RBF) method , finite point method (FPM) , meshless local Petrov-Galerkin (MLPG) method , reproducing kernel particle method (RKPM) , meshless manifold method , boundary element-free method (BEFM) , and local boundary integral equation (LBIE) method [41, 42].

The element-free Galerkin (EFG) method is one of the most powerful meshless methods . The EFG method can obtain a solution with high precision. Various problems have been successfully analyzed by the EFG method. By using the orthogonal function system with a weight function as the basis function, Zhang et al. presented the improved element-free Galerkin method , which has high computational efficiency. By combining the complex variable moving least-squares (CVMLS) approximation and the EFG method, Peng et al. presented the complex variable element-free Galerkin (CVEFG) method . Compared with the conventional EFG methods, the CVEFG method has greater computational precision and efficiency.

The EFG method is constructed based on the moving least-squares (MLS) approximation. The shape function that is formed with MLS approximation can obtain a solution with high precision. However, a disadvantage of the MLS approximation is that its shape function does not satisfy the property of Kronecker function. Then the EFG method based on the MLS approximation cannot apply the essential boundary conditions directly and easily. The essential boundary conditions need to be introduced by additional approaches, such as Lagrange multipliers  and penalty methods . However, for Lagrange multipliers, the corresponding discrete system will introduce additional unknowns which are not directly associated with the solution themselves. Furthermore, the banded structure of the matrix equation system is seriously worsened, as well as the conditioning properties. And, for penalty methods, the optimal value of penalty factor is hard to be set, which always affects the accuracy of the final solution.

To overcome this problem, Most and Bucher designed a regularized weight function with a regularization parameter , by which the MLS approximation can almost fulfill the interpolation with high accuracy . Most and Bucher enhanced the regularized weighting function to obtain a true interpolation MLS approximation . Another possible solution for this problem is the interpolating moving least-squares (IMLS) method presented by Lancaster and Salkauskas . The IMLS method is established based on the MLS approximation by using singular weight functions. The shape function of the IMLS method satisfies the property of Kronecker function. Thus the meshless methods based on the IMLS method can apply the essential boundary condition directly without any additional numerical effort. Based on the IMLS method, Kaljević and Saigal  presented an improved EFG method, in which the boundary condition is applied directly. Ren improved the expression of the shape function of the IMLS method and then presented the interpolating element-free Galerkin (IEFG) method and the interpolating boundary element-free (IBEFG) method for two-dimensional potential and elasticity problems . To overcome the singularity of the weight function in the IMLS method, Wang et al. presented the improved interpolating moving least-squares (IIMLS) method, in which nonsingular weight function is used . In the IEFG method, the essential boundary conditions are applied directly and easily, and the number of unknown coefficients in the trial function of the IMLS method is less than that in the trial function of the MLS approximation. Therefore, the IEFG method based on the IMLS method has high computational efficiency and precision.

Error estimation for meshless method is certainly important to increase the reliability and reduce the cost of numerical computations in many engineering problems. Some error analyses have been done for the MLS approximation and the meshless method based on it . Krysl and Belytschko studied the convergence of the continuous and discontinuous shape functions of the second-order elliptic partial differential equations . Chung and Belytschko proposed the local and global error estimates using the difference between the values of the projected stress and these given directly by the EFG solution . Dolbow and Belytschko studied the integration error . Gavete et al. presented a procedure to estimate the error in elliptic equations and then proposed a posteriori error approximation [78, 79]. R. J. Cheng and Y. M. Cheng studied the error estimate of the finite point method based on the MLS approximation  and the error estimates of element-free Galerkin method for potential and elasticity problems [81, 82]. For the IEFG method, since the essential boundary condition is applied directly, then the error estimate of the IEFG method is no doubt different from that of the EFG method. However, until now the error analysis of the IEFG method has not been seen in the recent literature.

Two-point boundary value problems occur in applied mathematics, theoretical physics, engineering, and optimization theory. Since it is usually impossible to obtain analytical solutions to two-point boundary value problems met in practice, these problems must be attacked by numerical methods. Many numerical methods have been proposed for the solutions of these problems, such as the Galerkin and collocation methods, boundary value method, variational iteration method, and meshless method based radial basis functions.

In this paper, the IMLS method is discussed in detail. The computation of the shape function of this paper is simpler than the corresponding expression presented by Lancaster and Salkauskas. Then based on the IMLS method of this paper and the Galerkin weak form, an IEFG method for two-point boundary value problems is presented. Since the shape function of the IMLS method satisfies the property of Kronecker function, then the IEFG method can apply the essential boundary condition directly. And as the number of the coefficients in the trial function of the IMLS method is less than that in the MLS approximation, then fewer nodes are selected in the entire domain in the IEFG method than in the conventional EFG method. Hence, the IEFG method has high computing speed and precision.

Then the error analysis of the IEFG method for two-point boundary value problems is presented. The convergence rates of the numerical solution and its derivatives of the IEFG method are presented. The theoretical results show that if the exact solution is smooth enough and the order of the polynomial basis functions is big enough, then the solution of the IEFG method and its derivatives are convergent to the exact solutions in terms of the maximum radius of the domains of influence of nodes. For the purpose of demonstration, some selected numerical examples are given to confirm the theory.

#### 2. Interpolating Moving Least-Squares Method

Let be a set of all nodes in the bounded domain , where is the number of nodes. The parameter denotes the radius of the domain of influence of node , and denotes the Euclidean norm. The domain of influence of is defined by

Let be the function of the field variable defined in . The approximation function of is denoted by . In this paper, the following weight function is used: where satisfying for and for . In general, can be chosen to be any weight function of the MLS approximation.

For a given point , the inner product is defined as where is the number of nodes whose compact support domains cover .

Let and let be given basis functions. Then a new set of basis functions will be generated from these given basis. Let where is a linear operator defined as

Define a local approximation as where is the point in the domain of influence of , and () are the unknown coefficients of basis functions.

Then define a functional as where shown in (2) is a weight function with compact support, are the nodes with domains of influence that cover the point , and .

By minimizing the functional , we have

From (4) and (10), we have Then (11) reduces to

In , the unknown parameters are solved from (13). In fact, (13) can be simplified.

If the weight function of (2) is used, , it can be proved that there exists

Then (13) can be simplified as

Equation (15) is simpler than the expression (13) presented in  and can be rewritten as where and is a matrix, and

From (16), we have

Then the local approximation function is obtained as

Thus the global interpolating approximation function of can be obtained as where is the shape function vector as

Equation (25) is the shape function of the IMLS method, and then the IMLS method is presented.

Equations (13) and (15) show that the computation of the shape function of this paper is simpler than the corresponding expression in . The reduction of computational amount is related to the order of the operations. And, in any case, the calculation amount of the shape functions of this paper is at least multiplication operations less than that of the IMLS method in .

#### 3. Interpolating Element-Free Galerkin Method for Two-Point Boundary Value Problems

Consider the following two-point boundary value problem: where , and are known sufficiently smooth functions and . , suppose that problem (29) has a unique solution.

The Galerkin weak form of (29) is where .

From the IMLS method, the unknown solution at arbitrary field point in the interval can be expressed as where is the number of nodes whose compact support domains cover the point .

Substituting (31) into (30) yields Because the nodal test function is arbitrary, the final discretized equation of (29) is obtained as where

Since the shape function of the IMLS method satisfies the property of Kronecker function, then the essential boundary conditions can be applied directly. Substituting the boundary conditions into (33) directly, we can obtain the unknowns at nodes by solving the linear equations (33).

To evaluate the integrals in (34), it is necessary to generate integration cells over the whole domain of the problem. These cells can be defined arbitrarily, but a sufficient number of quadrature points must be used to obtain a well-conditioned and nonsingular system of (33). In one dimension, one example of a cell structure is to set the quadrature cells equal to the intervals between the nodes. Once the cells and corresponding quadrature points are defined, the discrete equations are assembled by looping over each quadrature point.

The numerical procedure of IEFG method for two-point boundary value problems is listed as follows.(1)Looping over background cells to determine all Gauss points to find out its location and weight.(2)Looping over Gauss points for integration of (34) to(a)determine the support domain for specified Gauss point and select neighboring nodes based on a defined criterion;(b)compute shape function and its derivatives for each Gauss point;(c)assemble the contribution of each Gauss point to form system equation.(3)Enforcing essential (displacement) boundary conditions.(4)Solving the system equation to obtain nodal displacements.

Thus the IEFG method is presented for two-point boundary value problems.

#### 4. Error Estimates

In this section, the error analysis of the IEFG method for two-point boundary value problems is presented. The convergence rates of the numerical solution and its derivatives of the IEFG method are presented.

Let be a set of all nodes in the interval , where is the number of nodes. Let and , where is the radius of the domain of influence of node . For a given , denotes the maximum radius of the influence domains of nodes whose compact support domains cover . And suppose that there exist constants and such that and , respectively.

Define the Lebesgue space as where

The Sobolev space is defined as where

Define

Then the variational problem in accordance with (29) is to find such that where

Suppose that the bilinear form on the Sobolev space is bounded and coercive; that is, there exist constants and such that

Since the shape function of the IMLS method satisfies the property of Kronecker function, then the finite-dimensional solution space of the IEFG method for two-point boundary value problems is defined as

Then the IEFG method according to (29) is to find such that

Obviously, there exists . Hence, the IEFG method for two-point boundary value problems has a unique solution. And the following theorem can be obtained.

Theorem 1. Suppose that is the solution of the variational problem (40) and is the solution of the IEFG method (44). Then there exist(a), ;(b);(c).

In fact, the approximation function of the IMLS method provides a linear operator defined as

If , then we have proved that there exist bounded function and constant such that

Let . Then the following error estimates of the energy norm and the norm can be obtained.

Theorem 2. Suppose that . Let and be, respectively, the solutions of the problems (40) and (44). Then there exist and , which are independent with the parameter , such that

Proof. From Theorem 1 and (47) we have Then this theorem holds.
By using the Aubin-Nitsche method, the following error estimates in the norm can be obtained.

Theorem 3. Suppose that . Let and be, respectively, the solutions of the problems (40) and (44). Then there exists a constant , which is independent with the parameter , such that

Proof. For , let be the solution to the adjoint problems If the coefficients of the bilinear form are sufficiently smooth, then there exists the following estimate:
Let . Then we have
For arbitrary , from Theorem 1, we have
It follows from (53) and (54) that
In (55), if we let   and  , then there exists
From (52) and (56), we have Then this theorem is proved.

To study the error estimates of the high derivatives of the numerical solution of the IEFG method, we need to firstly prove the following inverse estimates of the function in the shape function space.

Theorem 4. Suppose that is defined by (25). Then ; there exists a constant , which is independent with the parameter , such that

Proof. From (46), there exists a bounded function independent with such that where .
It is obvious that there exist bounded functions and independent with such that
Then we have
Then it follows from (61) that where is independent with .
From (62), we have And there exists
From (63) and (64), we have
Then from (62) and (65), this theorem holds.
By using the inverse estimates, the following error estimates of the high derivatives of the numerical solution can be obtained.

Theorem 5. Suppose that . Let and be, respectively, the solutions of the problems (40) and (44). Then there exists a constant independent with , such that

Proof. From Theorem 4 we have Then it follows from Theorem 2 and (47) that There certainly exits From (68), (69), and (47), we have Then this theorem holds.

#### 5. Numerical Examples

In this section, two numerical examples are presented to show the applicability and the theoretical error estimates of the IEFG method of this paper. In our numerical computation, the nodes are arranged regularly, and the radius of the domain of influence of node is determined by , where is a positive scalar. The value of must be chosen so that the solution of (22) exists. In our following examples, and . Define the error norms as where and are, respectively, the analytical and numerical solutions. The integration in (71) is obtained numerically by fourth-order Gaussian quadrature. The in (2) is chosen to be the cubic spline weight function.

The first example considered is a linear elastostatics problem. A one-dimensional bar of unit length is subjected to a body force of magnitude . The displacement of the bar is fixed at the left end, and the right end is traction free. The bar has a constant cross sectional area of the unit value, and the elastic modulus is . Then the equilibrium equation and boundary conditions of this problem are

The analytical solution to the above problem is

Let . Under the quadratic basis functions and 21 regular distributed nodes, the analytical and numerical solutions of the displacement and strains along the bar are shown, respectively, in Figures 1 and 2, where the numerical values of the IEFG method are in good agreement with the exact ones.

The absolute errors of the displacements obtained by the IEFG and EFG methods are shown in Figure 3. Here, the essential boundary conditions are enforced by the penalty method in the EFG method, and the penalty factor is chosen to be . The CPU times to obtain these results by using the IEFG and EFG methods are, respectively, 0.0267 s and 0.0253 s. It can be seen that the IEFG method in this paper has higher precision than the EFG method under the similar CPU time.

The error norms of , , and under the quadratic and linear basis functions are, respectively, shown in Figures 4 and 5. The convergence rates of , , and with quadratic basis are, respectively, about 3, 2, and 1, and the convergence rates with linear basis are 2, 1, and . It is also shown that the second derivatives of the numerical solution of the IEFG method are not convergent to the exact values on the radius when the linear basis is used. It can be seen that these numerical results are in excellent agreement with the ones of the theories of the paper.

The second example considers the following equilibrium equation: with the boundary conditions The analytical solution of this example is

When the quadratic basis functions and 21 regular distributed nodes are used, the analytical and numerical solutions of and are shown, respectively, in Figures 6 and 7, where the numerical values of the IEFG method are also in accordance well with the exact ones.

The absolute errors of obtained by the IEFG and EFG methods are shown in Figure 8. In the EFG method, the cubic spline weight function is used, and the essential boundary conditions are enforced by the penalty method. The penalty factor is chosen to be . The CPU times to obtain these results by the IEFG and EFG methods are, respectively, 0.0361 s and 0.0325 s. Again, the IEFG method has higher precision than the EFG method under the similar CPU time.

The convergence rates of , , and are shown in Figures 9 and 10, respectively, with the quadratic and linear basis functions. The convergence rates of , , and with the quadratic basis are, respectively, 3, 2, and 1. And the corresponding rates with linear basis are, respectively, 2, 1, and . Figure 10 also shows that the second derivatives of the numerical solution of the IEFG method are not convergent to the exact values on the radius under the linear basis. It is also evident that these numerical results agree well with the ones of the theories of the paper.

#### 6. Conclusions

In this paper, the IMLS method is discussed in detail. The computation of the shape function of this paper is simpler than the corresponding expression presented by Lancaster and Salkauskas. Then, based on the IMLS method of this paper and the Galerkin weak form, an IEFG method for two-point boundary value problems is presented. Since the shape function of the IMLS method satisfies the property of Kronecker function, then the IEFG method can apply the essential boundary condition directly. And as the number of the coefficients in the trial function of the IMLS method is less than that in the MLS approximation, then fewer nodes are selected in the entire domain in the IEFG method than in the conventional EFG method. Hence, the IEFG method has high computing speed and precision.

The error analysis of the IEFG method for two-point boundary value problems is presented. The convergence rates of the numerical solution and its derivatives of the IEFG method are presented. The theories of this paper show that if the analytical solution is sufficiently smooth and the order of the polynomial basis functions is big enough, then the solution of the IEFG method and its derivatives are convergent to the analytical solutions in terms of the maximum radius of the domains of influence of nodes. For the purpose of demonstration, some selected numerical examples are given to confirm the theories.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 11171208), Shanghai Leading Academic Discipline Project (no. S30106), and the Natural Science Foundation of Ningbo (no. 2013A610103).