Next Article in Journal
Finance for the Environment: A Scientometrics Analysis of Green Finance
Next Article in Special Issue
Local Dynamics of Logistic Equation with Delay and Diffusion
Previous Article in Journal
Solution of Moore–Gibson–Thompson Equation of an Unbounded Medium with a Cylindrical Hole
Previous Article in Special Issue
Oscillation of Second-Order Differential Equations with Multiple and Mixed Delays under a Canonical Operator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Infinitely Smooth Polyharmonic RBF Collocation Method for Numerical Solution of Elliptic PDEs

1
Center of Excellence for Ocean Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
2
Department of Harbor and River Engineering, School of Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(13), 1535; https://doi.org/10.3390/math9131535
Submission received: 5 May 2021 / Revised: 23 June 2021 / Accepted: 25 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Recent Advances in Differential Equations and Applications)

Abstract

:
In this article, a novel infinitely smooth polyharmonic radial basis function (PRBF) collocation method for solving elliptic partial differential equations (PDEs) is presented. The PRBF with natural logarithm is a piecewise smooth function in the conventional radial basis function collocation method for solving governing equations. We converted the piecewise smooth PRBF into an infinitely smooth PRBF using source points collocated outside the domain to ensure that the radial distance was always greater than zero to avoid the singularity of the conventional PRBF. Accordingly, the PRBF and its derivatives in the governing PDEs were always continuous. The seismic wave propagation problem, groundwater flow problem, unsaturated flow problem, and groundwater contamination problem were investigated to reveal the robustness of the proposed PRBF. Comparisons of the conventional PRBF with the proposed method were carried out as well. The results illustrate that the proposed approach could provide more accurate solutions for solving PDEs than the conventional PRBF, even with the optimal order. Furthermore, we also demonstrated that techniques designed to deal with the singularity in the original piecewise smooth PRBF are no longer required.

1. Introduction

A large number of multidisciplinary problems may require solutions involving complicated mathematical models and numerical techniques [1,2]. Among these approaches, the meshfree method has gained the attention of researchers from different scientific fields due to its capability of dealing with partial differential equations (PDEs) with complicated and irregular geometry [3,4,5,6,7]. The radial basis function collocation method (RBFCM) is one meshfree approach for analyzing governing equations, where the unknowns are represented by a function approximation [8,9,10]. The RBFCM is also capable of solving linear and nonlinear PDEs requiring solutions involveing sophisticated mathematical models and numerical algorithms [8,9]. Function approximation using the interpolation method of the radial basis function (RBF) was proposed in 1971 by Hardy [11], who used the multiquadric (MQ) RBF for scattered data interpolation to solve problems. In addition to the MQ RBF, several RBFs may be found, such as the inverse multiquadric (IMQ), Gaussian, and polyharmonic spline (PS) functions [12,13,14]. Among them, the PS and the MQ RBFs may provide much more accurate solutions than other RBFs [15,16]. The RBFs are usually categorized into piecewise and infinitely smooth functions for solving PDEs in the RBFCM. For example, the PS is piecewise smooth. On the other hand, the MQ is infinitely smooth, with the shape parameter [17,18]. For the PS, a minimum order of two must be adopted to keep it smooth for solving Laplace-type equations [19]. Despite the great success of the above RBFs as effective numerical techniques for dealing with several kinds of PDEs, there is still growing interest in the application and development of new and advanced RBFs [20]. A significant number of modifications to RBFs have been proposed, such as the pseudo-spectral RBF [21,22], Gaussian RBF [23], RBF QR alternative basis method [24], finite difference RBF [25,26], partition of unity RBF [27,28], stabilized expansion of the Gaussian RBF [29], rational RBF [30,31], and RBF based on partition of unity of Taylor series expansion [32].
An adaptive greedy technique for the Kansa method is adopted for general PDE problems to be solved on complicated domains [5]. The adaptive greedy algorithm prefers exterior centers for large scaling parameters. Similar to the idea of exterior centers, Ku et al. presented a fictitious source collocation scheme where the centers are replaced by fictitious sources collocated outside the domain as well [33,34]. However, different from the greedy algorithm, the role of fictitious sources is emphasized on the prevention of singularity without the shape parameter for the MQ. Since the distance between the source and the interior point is always greater than zero, the MQ and its derivatives are always smooth and globally infinitely differentiable without using the shape parameter.
Motivated by this concept, a novel infinitely smooth polyharmonic radial basis function (PRBF) collocation method for solving elliptic partial differential equations (PDEs) is presented. The seismic wave propagation problem, groundwater flow problem, unsaturated flow problem, and groundwater contamination problem were investigated to reveal the robustness of the proposed PRBF. Comparisons of the conventional PRBF with the proposed method were carried out as well. The formulation of the PRBF is introduced in Section 2. The convergence analysis is presented in Section 3. Section 4 details the results of the investigation into the geomechanical problems listed above. The findings of this research are finally summarized in Section 5.

2. The Polyharmonic Radial Basis Function Collocation Method

2.1. The Governing Equation

The governing equation for elliptic partial differential equations in two dimensions can be expressed as
2 u ( x ) + A u ( x ) + B u ( x ) = C ( x )   in   Ω
u ( x ) = g ( x )   on   Ω
where denotes the divergence operator; denotes the dot product; u ( x ) denotes the unknown, x = ( x , y ) ; Ω denotes a domain with the boundary Ω ; g ( x ) denotes the boundary value; A is the velocity in two dimensions, A = ( A x ,   A y ) ; B is the given function; and C ( x ) is the given function value.
In the RBFCM, several RBFs can be implemented, such as the MQ, Gaussian, IMQ, and PS functions. However, the accuracy of the Gaussian, MQ, and IMQ RBFs may depend on the shape parameter. Thus, most applications of the Gaussian, MQ, and IMQ RBFs utilize optimization approaches or experimental techniques to determine the optimal shape parameter [17,18]. Since the PS may provide good agreement without providing an additional parameter to prevent singularity, we considered the PRBF in this study. The PRBF with natural logarithm is defined as
φ ( r ) = r 2 k ln ( r ) , k = 1 , 2 , 3
where k is the order; r is the radial distance, r = x x s ; x denotes the interior point; and x s denotes the center point. The equation is a kth order PRBF.
Taking the derivative of (3) with respect to x gives
φ ( r ) x = ( x x s ) r 2 k 2 [ 2 k ln ( r ) + 1 ]
where x s is the x-coordinate of the center point. We may then take the derivative of y and obtain
φ ( r ) y = ( y y s ) r 2 k 2 [ 2 k ln ( r ) + 1 ]
where y s is the y-coordinate of the center point. Similarly, we have
2 φ ( r ) x 2 = ( x x s ) 2 r 2 k 4 [ ( 4 k 2 4 k ) ln r + 4 k 2 ] + r 2 k 2 [ 2 k ln ( r ) + 1 ] .
Again, we take the derivative of y and obtain
2 φ ( r ) y 2 = ( y y s ) 2 r 2 k 4 [ ( 4 k 2 4 k ) ln r + 4 k 2 ] + r 2 k 2 [ 2 k ln ( r ) + 1 ] .
Utilizing the PRBF to build up the function approximation, we obtain
u ( x ) = j = 1 n c λ j φ ( r j )
where n c denotes the total number of center points; λ j denotes the unknown coefficient; φ ( r j ) denotes the basis function; r j is the radial distance at the jth center point, r j = x x j s ; and x j s denotes the jth center point, described as x j s = ( x j s , y j s ) .
Inserting the aforementioned equations into (1), the function approximation of the governing equation is as follows:
j = 1 n c λ j r j 2 k 2 ( 4 k 2 ln r j + 4 k ) + j = 1 n c λ j [ A x ( x x j s ) + A y ( y y j s ) ] r j 2 k 2 [ 2 k ln ( r j ) + 1 ] + j = 1 n c λ j B r j 2 k ln ( r j ) = C ( x ) .
From (9), it is clear that ln ( r j ) becomes singular when r j = 0 , where the position of the center point exactly coincides with the position of the interior point. In this study, we propose a novel idea to convert the piecewise smooth PRBF into an infinitely smooth PRBF using source points collocated outside the domain to ensure that the radial distance is always greater than zero, thus avoiding the singularity of the conventional PRBF. Accordingly, the terms in (9) are always smooth, and no extra techniques are required to deal with the singularity. The concept of fictitious sources is elaborated in the following section.

2.2. The Fictitious Sources

In the conventional PRBF, the interior, center, and boundary points have to be placed where the positions of the interior and center points are usually the same, as displayed in Figure 1a. To prevent singularity, the PRBF with natural logarithm is usually implemented as follows [15]:
φ ( r ) = r 2 k 1 ln ( r r ) , k = 1 , 2 , 3 .
Using (8), we may obtain the function approximation for the governing equation again
j = 1 n s λ j r j 2 k 3 [ 4 k 2 ln ( r j r j ) + 4 k r j ] + j = 1 n s λ j [ A x ( x x j s ) + A y ( y y j s ) ] r j 2 k 3 [ 2 k ln ( r j r j ) + r j ] + j = 1 n s λ j B r j 2 k 1 ln ( r j r j ) = C ( x ) .
where n s denotes the number of fictitious sources. To prevent singularity, a minimum order of two must be utilized to ensure that r j 2 k 3 in the above equation is a smooth and completely monotone RBF. Accordingly, the PRBF with natural logarithm is a piecewise smooth function in the conventional RBFCM.
In this study, fictitious sources are adopted. The center points in the PRBF are regarded as the fictitious sources, where the fictitious sources are placed outside the domain, as depicted in Figure 1b.
The following equation is utilized for placing the source points.
x j s = ( η ρ j s cos θ j s , η ρ j s sin θ j s ) , j = 1 , ,   n s
where η denotes the dilation parameter; ρ j s denotes the radius of the fictitious sources at the jth source point; and θ j s denotes the angle of the fictitious sources at the jth source point. Figure 1b depicts the proposed fictitious sources. As shown in Figure 1b, the positions of the source points can, therefore, be collocated with the consideration of different values of the dilation parameter by utilizing (12). Since the radial distance is always nonzero, the PRBF and its derivatives in the governing equation, as depicted in (3) and (9), are always continuous. Based on the proposed fictitious sources, the original piecewise smooth PRBF is then converted into an infinitely smooth PRBF.
To investigate the unknown coefficients, approximations are applied to satisfy (1) with the boundary conditions, as shown in (2), at each point. The following linear system is acquired:
[ A L ] n i × n s [ A B ] n b × n s λ = [ f ] n i × 1 [ g ] n b × 1
where A L denotes a n i × n s matrix for the interior points; A B denotes a n b × n s matrix for the boundary points; λ denotes a n s × 1 vector which includes unknown coefficients λ j to be solved; f denotes a n i × 1 vector of function values for the interior points, written as f = [ C 1 ,   C 2 ,   ,   C n i ] ; g denotes an n b × 1 vector of boundary data, written as g = [ g 1 ,   g 2 ,   ,   g n b ] ; n i denotes the number of interior points; and n b denotes the number of boundary points. The undetermined coefficients can be obtained by solving a system of simultaneous linear equations in matrix form as shown in (13).
Accordingly, the solution of u ( x ) can be obtained using the function approximation of (8). To evaluate the error, the root-mean-square error (RMSE) and the maximum absolute error (MAE) are both adopted. The accuracy of the computed results is compared to the exact solutions using the following equations.
RMSE = i = 1 I ( u E ( x i ) u ( x i ) ) 2 I ,
MAE = maximum   u E ( x i ) u ( x i ) ,   1 i I
where I is the total number of interior points, RMSE is the root-mean-square error, MAE is the maximum absolute error, and u ( x i ) and u E ( x i ) represent the numerical and exact solutions at the ith interior point, respectively.

3. Convergence Analysis

To compare the proposed PRBF with the conventional PRBF, we firstly consider a steady-state groundwater flow problem [33]:
2 u ( x ) = 0 .
An irregular boundary was defined as
Ω = ( x , y ) x = ρ ( θ ) cos θ ,   y = ρ ( θ ) sin θ , ρ ( θ ) = [ cos ( 3 θ ) + 2 sin 2 ( 3 θ ) ] 1 / 3 ,   0 θ 2 π .
The boundary data, as depicted in (18), were assigned on the irregular boundaries.
g ( x ) = e y sin ( x ) + e x cos ( y ) .
We used 321 interior points, 396 fictitious sources, and 75 boundary points. In this example, we consider six cases with different shapes of fictitious sources using (19) to (24) as demonstrated in Figure 2a–f for Cases I, II, III, IV, V, and VI, respectively.
ρ j 1 s ( θ j s ) = [ cos ( 3 θ j s ) + 2 sin 2 ( 3 θ j s ) ] 1 / 3 ,   0 θ j s 2 π ,
ρ j 2 s ( θ j s ) = e ( sin θ j s sin 2 θ j s ) 2 + e ( cos θ j s cos 2 θ j s ) 2 ,   0 θ j s 2 π ,
ρ j 3 s ( θ j s ) = 1 + cos 2 ( 6 θ j s ) ,   0 θ j s 2 π ,
ρ j 4 s ( θ j s ) = 1 + sin 4 ( 2 θ j s ) ,   0 θ j s 2 π ,
ρ j 5 s ( θ j s ) = 6 cos 4 ( 4 θ j s ) ,   0 θ j s 2 π ,
ρ j 6 s ( θ j s ) = 4 + sin 5 ( 5 θ j s ) ,   0 θ j s 2 π .
The fictitious sources were placed on the external domain by utilizing (25). The fictitious source boundaries were defined as [33,34]
Ω s = ( x j s , y j s ) x j s = η ρ j s cos θ j s ,   y j s = η ρ j s sin θ j s
where Ω s denotes the boundaries of fictitious sources. Figure 3 illustrates the MAE values of the conventional and proposed PRBFs with different orders (values of k). In Case I, the MAE tended to increase with an increase in the order k, as displayed in Figure 3a. In the other five cases, it was found that the MAE of the proposed PRBF remained lower than that of the conventional one with different orders, as displayed in Figure 3b–f. The accuracy is less sensitive to the order k. From Figure 3, accurate results with order k in the range of one to ten were acquired. Additionally, the results demonstrate that the proposed PRBF with an order of one provides more accurate results than the conventional one, as depicted in Table 1. The MAE of the conventional PRBF was in the order of 10 1 to 10 4 as the order k ranged from two to ten. It appears that the MAE for the conventional PRBF could only reach the order of 10 4 when the order k was set to seven, as displayed in Figure 3.
The MAE of the proposed PRBF with different dilation parameters was investigated, as displayed in Figure 4. In Cases I, II, III, and IV, the MAE declined with an increase in the dilation parameters—thus, all cases could obtain accurate results. As depicted in Figure 4e,f these results show that accurate results in the order of 10 6 to 10 11 were yielded when an arbitrary dilation parameter in the range of zero to ten was utilized.
It is found that the error may grow with the increase of the dilation parameter. However, results from all cases reveal that we may obtain accurate results with the use of an arbitrary value in the range of 1.5 to 5 for the dilation parameter. Nevertheless, the error may increase while the larger dilation parameter is adopted. However, it reveals that we could still obtain better accuracy than original PS even in the worst case.

4. Numerical Examples

To demonstrate the applicability of the proposed method, several types of geomechanical problems—the seismic wave propagation problem, groundwater flow problem, unsaturated flow problem, and groundwater contamination problem—were investigated; the results are presented in this section.

4.1. The Seismic Wave Propagation Problem

The problem of seismic wave propagation [35] is considered as
( 2 λ 2 ) u ( x , y ) = 0
where λ 2 denotes the wave number depending on the property of the geomaterial. The boundary is
Ω = ( x , y ) x = ρ ( θ ) cos θ ,   y = ρ ( θ ) sin θ , ρ ( θ ) = 2 + 1.3 cos ( θ ) ,   0 θ 2 π .
Data are assigned on the boundaries by the exact solution:
g ( x ) = e ( λ / 2 ) ( x + 3 y ) .
In this wave problem, we assigned λ 2 = 1 and λ 2 = 0.0001 . The conventional PRBF and the proposed PRBF with fictitious sources were utilized to solve this problem. In the conventional PS RBF, the order k 2 was adopted to prevent the singularity. In the proposed PRBF with fictitious sources, the irregular boundaries of the exterior sources, Ω s , defined by (24), were utilized in this example. Figure 5 shows the original collocation points for the conventional case. We used 480 center points and 75 boundary points. For the fictitious sources, demonstrated in Figure 6, the interior points, boundary points, and fictitious sources were placed within, on, and outside the region, respectively. The numbers of source, boundary, and interior points were 480, 75, and 405, respectively.
Figure 7a shows the order k versus the MAE when the wave number was 0.0001. From Figure 7a, the MAE of the conventional PRBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. It appears that the proposed approach with order of one acquired the most accurate results, with MAE in the order of 10 13 —much better than the conventional PS RBF. Figure 7b shows the order k versus the MAE when the wave number was 1. From Figure 7b, the MAE of the conventional PS RBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. However, the MAE of the proposed PRBF utilizing fictitious sources drastically decreased from 10 2 to 10 10 as the order k ranged from one to ten. It is clear that the proposed PRBF with fictitious sources can provide significantly more accurate results than the conventional PRBF. Figure 8 shows the MAE versus the dilation parameter. The results show that the locations of the fictitious sources are of low sensitivity with respect to the solutions.

4.2. The Groundwater Flow Problem

The Poisson equation [36] is referred to for groundwater flow with recharge, defined as
2 u ( x , y ) = p 2 y cos ( p y ) + y sin ( p x )
where p denotes the parameter depending on the recharge; values of p of 1 and 0.0001 were considered in this problem. The boundary is
Ω = ( x , y ) x = ρ ( θ ) cos θ ,   y = ρ ( θ ) sin θ , ρ ( θ ) = sin 3 θ + 5 sin 2 3 θ ,   0 θ 2 π .
The Dirichlet data are
g ( x ) = y sin ( p x ) + x cos ( p y ) .
The conventional PRBF and the proposed PRBF with fictitious sources were adopted to solve this problem. In the conventional PRBF, the order k 2 was adopted to prevent the singularity. In the proposed PRBF with fictitious sources, irregular boundaries of the exterior sources, Ω s , defined by (24), were utilized in this example. Figure 9 depicts the collocation points for the conventional case. We considered 382 center points and 75 boundary points. For the proposed fictitious sources, as demonstrated in Figure 10, we assigned the interior, boundary, and fictitious sources within, on, and outside the region, respectively. The numbers of fictitious source, boundary, and interior points were 457, 75, and 382, respectively.
Figure 11a shows the order k versus the MAE when value of p was set to 0.0001. From Figure 11a, the MAE of the conventional PRBF fluctuated between 10 1 and 10 5 as the order k ranged from two to ten. However, it was found that the proposed method with order of one acquired significantly more accurate results than the conventional method. Figure 11b shows the order k versus the MAE when the value of p was set to 1. From Figure 11b, the MAE of the conventional PRBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. However, the MAE of the proposed PRBF with fictitious sources drastically decreased from 10 1 to 10 10 as the order k ranged from one to ten. It is clear that the proposed PRBF with fictitious sources can provide more accurate results than the conventional one. Figure 12 shows the MAE versus the dilation parameter. The results show that the locations of the fictitious sources may not be sensitive to the solutions.

4.3. The Unsaturated Flow Problem

The third example is the stationary convection–diffusion problem, usually referred to as unsaturated flow in two dimensions [1].
x K r ( h ) h ( x , y ) x + y K r ( h ) ( h ( x , y ) y + 1 ) = 0 .
Here, K r ( h ) denotes the relative hydraulic conductivity, and h denotes the pressure head. The following equation is obtained by utilizing the Gardner model [1]:
2 h ^ ( x , y ) x 2 + 2 h ^ ( x , y ) y 2 + α g h ^ ( x , y ) y = 0
where α g denotes the pore size distribution parameter; h ^ denotes the linearized pressure head, expressed as h ^ = e α g h e α g h d ; and h d denotes the dry pressure head. The irregular boundary is
Ω = ( x , y ) x = ρ ( θ ) cos θ ,   y = ρ ( θ ) sin θ , ρ ( θ ) = cos 2 θ + 2 + sin 2 2 θ ,   0 θ 2 π .
The boundary data are
g ( x ) = x e α g y .
In this example, soil pore size distribution values of α g =   0.0001 and α g = 1 were considered. The conventional PRBF and the proposed PRBF with fictitious sources were implemented to solve this problem. In the conventional PRBF, the order k 2 was adopted to prevent the singularity. In the proposed method with fictitious sources, irregular boundaries of the exterior sources, Ω s , defined by (24), were utilized in this example. Figure 13 shows the collocation points for the conventional case. We used 502 center points and 75 boundary points. As demonstrated in Figure 14, we collocated the interior, boundary, and fictitious sources within, on, and outside the domain, respectively. The numbers of source, boundary, and interior points were 502, 75, and 427, respectively.
Figure 15a shows the order k versus the MAE when α g = 0.0001 . From Figure 15a, the MAE of the conventional PRBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. However, it was found that the proposed method with order of one acquired significantly more accurate results than the conventional method. Figure 15b shows the order k versus the MAE when the pore size distribution parameter was set to 1. From Figure 15b, the MAE of the conventional PS RBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. The MAE of the proposed method with fictitious sources drastically decreased from 10 1 to 10 9 as the order k ranged from one to ten. It is clear that the proposed PRBF with fictitious sources can provide more accurate results than the conventional one. Figure 16 shows the MAE versus the dilation parameter. The results show that the positions of the fictitious sources may not be sensitive to the solutions.

4.4. The Groundwater Contamination Problem

The last example is the convection–diffusion–reaction problem, usually referred to as the groundwater contamination problem [37]:
2 u ( x ) + A u ( x ) + B u ( x ) = C ( x )
where A denotes the velocity field, written as A = ( A x ,   A y ) , where A x = y 2 sin x and A y = x e y ; B = sin ( x ) + cos ( y ) ; and C ( x ) is the given function value describing the reaction process, defined as follows:
C ( x ) = π 2 y sin ( π x ) + x cos ( π y ) + y 2 sin x π y cos ( π x ) + cos ( π y ) + x e y sin ( π x ) π x sin ( π y ) + sin x + cos y y sin ( π x ) + x cos ( π y ) .
In this example, the boundary is
Ω = ( x , y ) x = ρ ( θ ) cos θ ,   y = ρ ( θ ) sin θ , ρ ( θ ) = cos ( 5 θ ) + 5 sin 2 ( 5 θ ) 5 ,   0 θ 2 π .
The boundary data are
g ( x ) = y sin ( π x ) + x cos ( π y ) .
The conventional PRBF and the proposed PRBF with fictitious sources were implemented to solve this problem. In the conventional PRBF, the order k 2 was utilized to prevent the singularity. In the proposed method, irregular boundaries of the exterior sources, Ω s , defined by (24), were utilized in this example. Figure 17 shows the collocation points for the conventional case. We used 369 center points and 75 boundary points. For the fictitious sources, demonstrated in Figure 18, we assigned the interior, boundary, and source points within, on, and outside the domain, respectively. The numbers of source, boundary, and interior points were 444, 75, and 369, respectively.
Figure 19 shows the order k versus the MAE. From Figure 19, the MAE of the conventional PRBF fluctuated between 10 1 and 10 4 as the order k ranged from two to ten. However, the MAE of the proposed method utilizing fictitious sources drastically decreased from 10 1 to 10 8 as the order k ranged from one to ten. It is clear that the proposed PRBF with fictitious sources can provide more accurate results than the conventional PRBF. Figure 20 shows the MAE versus the dilation parameter. The results show that the positions of the fictitious sources may not be sensitive to the solutions.

5. Conclusions

An infinitely smooth polyharmonic radial basis function (PRBF) collocation method for solving elliptic partial differential equations (PDEs) was developed in this study. The significance of this study can be summarized as follows:
(1)
We convert the piecewise smooth PRBF into an infinitely smooth PRBF using source points collocated outside the domain. The distance between the interior point and fictitious source is, therefore, greater than zero. The PRBF and its derivatives in the governing equation always have the characteristics of global infinite differentiability and smoothness. Accordingly, techniques to deal with the singularity in the original PRBF are no longer required.
(2)
To clarify the possible influences for the positions of the sources on the accuracy, an analysis was conducted to evaluate the accuracy versus the shapes and positions of sources. The results show that the error may grow with the increase of the dilation parameter. However, results from all cases reveal that we may obtain accurate results with the use of an arbitrary value in the range of 1.5 to 5 for the dilation parameter. Nevertheless, the error may increase while the larger dilation parameter is adopted. However, it reveals that we could still obtain better accuracy than original PS even in the worst case. Furthermore, it appears that the accuracy may be less sensitive to the order k from the numerical results. It is found that accurate results with order k in the range of 3 to 10 were acquired.
(3)
The seismic wave propagation problem, groundwater flow problem, unsaturated flow problem, and groundwater contamination problem were investigated to show the robustness and accuracy of the proposed approach. From the results, the proposed infinitely smooth PRBF could provide much more accurate solutions than the conventional PRBF even with the optimal order.

Author Contributions

Conceptualization, C.-Y.L. and S.-M.H.; acquisition of data, L.-D.H. and S.-M.H.; methodology, C.-Y.K.; validation, L.-D.H. and C.-Y.L.; drafting the manuscript, C.-Y.L.; revising the manuscript, C.-Y.K. and C.-Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Datasets generated during this work are available on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, C.Y.; Ku, C.Y.; Huang, C.C.; Lin, D.G.; Yeih, W.C. Numerical solutions for groundwater flow in unsaturated layered soil with extreme physical property contrasts. Int. J. Nonlinear Sci. Numer. Simul. 2015, 16, 325–335. [Google Scholar] [CrossRef]
  2. Vu, T.D.; Ni, C.F.; Li, W.C.; Truong, M.H.; Hsu, S.M. Predictions of groundwater vulnerability and sustainability by an integrated index-overlay method and physical-based numerical model. J. Hydrol. 2021, 596, 126082. [Google Scholar] [CrossRef]
  3. Gu, Y. Meshfree methods and their comparisons. Int. J. Comput. Methods 2005, 2, 477–515. [Google Scholar] [CrossRef]
  4. Hon, Y.C.; Chen, W. Boundary knot method for 2D and 3D Helmholtz and convection–diffusion problems under complicated geometry. Int. J. Numer. Methods Eng. 2003, 56, 1931–1948. [Google Scholar] [CrossRef] [Green Version]
  5. Ling, L.; Opfer, R.; Schaback, R. Results on meshless collocation techniques. Eng. Anal. Bound. Elem. 2006, 30, 247–253. [Google Scholar] [CrossRef] [Green Version]
  6. Xiong, J.; Wen, J.; Liu, Y.-C. Localized boundary knot method for solving two-dimensional Laplace and Bi-Harmonic equations. Mathematics 2020, 8, 1218. [Google Scholar] [CrossRef]
  7. Li, P.W. Space–time generalized finite difference nonlinear model for solving unsteady Burgers’ equations. Appl. Math. Lett. 2020, 114, 106896. [Google Scholar] [CrossRef]
  8. Chen, W.; Fu, Z.J.; Chen, C.S. Recent Advances in Radial Basis Function Collocation Methods; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  9. Ma, Z.; Li, X.; Chen, C.S. Ghost point method using RBFs and polynomial basis functions. Appl. Math. Lett. 2021, 111, 106618. [Google Scholar] [CrossRef]
  10. Fornberg, B.; Natasha, F. Solving PDEs with radial basis functions. Acta Numer. 2015, 24, 215. [Google Scholar] [CrossRef] [Green Version]
  11. Chen, C.S.; Karageorghis, A.; Dou, F. A novel RBF collocation method using fictitious centres. Appl. Math. Lett. 2020, 101, 106069. [Google Scholar] [CrossRef]
  12. Hardy, R. Multiquadric equations of topography and other irregular surfaces. J. Geophys. Res. 1971, 76, 1905–1915. [Google Scholar] [CrossRef]
  13. Gu, J.; Jung, J.H. Adaptive Gaussian radial basis function methods for initial value problems: Construction and comparison with adaptive MQ RBF methods. J. Comput. Appl. Math. 2020, 381, 113036. [Google Scholar] [CrossRef]
  14. Beatson, R.K.; Light, W.A. Fast evaluation of radial basis functions: Methods for two-dimensional polyharmonic splines. IMA J. Numer. Anal. 1997, 17, 343–372. [Google Scholar] [CrossRef]
  15. Madych, W.R.; Nelson, S.A. Polyharmonic cardinal splines. J. Approx. Theory 1990, 60, 141–156. [Google Scholar] [CrossRef] [Green Version]
  16. Xiang, S.; Bi, Z.Y.; Jiang, S.X.; Jin, Y.X.; Yang, M.S. Thin plate spline radial basis function for the free vibration analysis of laminated composite shells. Compos. Struct. 2011, 93, 611–615. [Google Scholar] [CrossRef]
  17. Tsai, C.C.; Cheng, A.D.; Chen, C.S. Particular solutions of splines and monomials for polyharmonic and products of Helmholtz operators. Eng. Anal. Bound. Elem. 2009, 33, 514–521. [Google Scholar] [CrossRef]
  18. Sarra, S.A.; Sturgill, D. A random variable shape parameter strategy for radial basis function approximation methods. Eng. Anal. Bound. Elem. 2009, 33, 1239–1245. [Google Scholar] [CrossRef]
  19. Biazar, J.; Hosami, M. An interval for the shape parameter in radial basis function approximation. Appl. Math. Comput. 2017, 315, 131–149. [Google Scholar] [CrossRef]
  20. Segeth, K. Polyharmonic splines generated by multivariate smooth interpolation. Comput. Math. Appl. 2019, 78, 3067–3076. [Google Scholar] [CrossRef]
  21. Fasshauer, G.E. RBF collocation methods as pseudospectral methods. WIT Trans. Model. Simul. 2005, 39, 10. [Google Scholar]
  22. Krowiak, A. Radial basis function-based pseudospectral method for static analysis of thin plates. Eng. Anal. Bound. Elem. 2016, 71, 50–58. [Google Scholar] [CrossRef]
  23. Fasshauer Gregory, E.; Mccourt Michael, J. Stable evaluation of Gaussian radial basis function interpolants. SIAM J. Sci. Comput. 2012, 34, 737–762. [Google Scholar] [CrossRef]
  24. Gholampour, F.; Hesameddini, E.; Taleei, A. A global RBF-QR collocation technique for solving two-dimensional elliptic problems involving arbitrary interface. Eng. Comput. 2020, 1–19. [Google Scholar] [CrossRef]
  25. Gonzalez-Rodriguez, P.; Bayona, V.; Moscoso, M.; Kindelan, M. Laurent series based RBF-FD method to avoid ill-conditioning. Eng. Anal. Bound. Elem. 2015, 52, 24–31. [Google Scholar] [CrossRef] [Green Version]
  26. Bayona, V.; Flyer, N.; Fornberg, B.; Barnett, G.A. On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs. J. Comput. Phys. 2017, 332, 257–273. [Google Scholar] [CrossRef] [Green Version]
  27. Shcherbakov, V.; Larsson, E. Radial basis function partition of unity methods for pricing vanilla basket options. Comput. Math. Appl. 2016, 71, 185–200. [Google Scholar] [CrossRef]
  28. Cavoretto, R.; De Marchi, S.; De Rossi, A.; Perracchione, E.; Santin, G. Partition of unity interpolation using stable kernel-based techniques. Appl. Numer. Math. 2017, 116, 95–107. [Google Scholar] [CrossRef] [Green Version]
  29. Rashidinia, J.; Fasshauer, G.E.; Khasi, M. A stable method for the evaluation of Gaussian radial basis function solutions of interpolation and collocation problems. Comput. Math. Appl. 2016, 72, 178–193. [Google Scholar] [CrossRef]
  30. Wright, G.B.; Fornberg, B. Stable computations with flat radial basis functions using vector-valued rational approximations. J. Comput. Phys. 2017, 331, 137–156. [Google Scholar] [CrossRef] [Green Version]
  31. De Marchi, S.; Martínez, A.; Perracchione, E. Fast and stable rational RBF-based partition of unity interpolation. J. Comput. Appl. Math. 2019, 349, 331–343. [Google Scholar] [CrossRef]
  32. Bawazeer, S.A.; Baakeem, S.S.; Mohamad, A.A. New approach for radial basis function based on partition of unity of Taylor series expansion with respect to shape parameter. Algorithms 2021, 14, 1. [Google Scholar] [CrossRef]
  33. Ku, C.Y.; Liu, C.Y.; Xiao, J.E.; Hsu, S.M. Multiquadrics without the shape parameter for solving partial differential equations. Symmetry 2020, 12, 1813. [Google Scholar] [CrossRef]
  34. Ku, C.Y.; Liu, C.Y.; Xiao, J.E.; Hsu, S.M.; Yeih, W. A collocation method with space–time radial polynomials for inverse heat conduction problems. Eng. Anal. Bound. Elem. 2021, 122, 117–131. [Google Scholar] [CrossRef]
  35. Cao, L.; Qin, Q.H.; Zhao, N. An RBF–MFS model for analysing thermal behaviour of skin tissues. Int. J. Heat Mass Transf. 2010, 53, 1298–1307. [Google Scholar] [CrossRef]
  36. Mai-Duy, N.; Tran-Cong, T. Indirect RBFN method with thin plate splines for numerical solution of differential equations. CMES Comp. Model. Eng. Sci. 2003, 4, 85–102. [Google Scholar]
  37. Li, J.; Zhai, S.; Weng, Z.; Feng, X. H-adaptive RBF-FD method for the high-dimensional convection-diffusion equation. Int. Commun. Heat Mass Transf. 2017, 89, 139–146. [Google Scholar] [CrossRef]
Figure 1. Location of the collocation points: (a) the conventional PRBF; (b) this study.
Figure 1. Location of the collocation points: (a) the conventional PRBF; (b) this study.
Mathematics 09 01535 g001
Figure 2. The fictitious sources at different positions: (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Figure 2. The fictitious sources at different positions: (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Mathematics 09 01535 g002
Figure 3. The computed MAE values for the original PRBF and the proposed PRBF with fictitious sources: (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Figure 3. The computed MAE values for the original PRBF and the proposed PRBF with fictitious sources: (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Mathematics 09 01535 g003
Figure 4. The computed MAE versus the dilation parameter ( k = 1 ) (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Figure 4. The computed MAE versus the dilation parameter ( k = 1 ) (a) Case I; (b) Case II; (c) Case III; (d) Case IV; (e) Case V; (f) Case VI.
Mathematics 09 01535 g004
Figure 5. The conventional radial basis function collocation method (RBFCM) collocation scheme.
Figure 5. The conventional radial basis function collocation method (RBFCM) collocation scheme.
Mathematics 09 01535 g005
Figure 6. The fictitious sources.
Figure 6. The fictitious sources.
Mathematics 09 01535 g006
Figure 7. The computed MAE values of the original PRBF and the proposed method: (a) λ 2 = 0.0001 ; (b) λ 2 = 1 .
Figure 7. The computed MAE values of the original PRBF and the proposed method: (a) λ 2 = 0.0001 ; (b) λ 2 = 1 .
Mathematics 09 01535 g007
Figure 8. Maximum absolute error versus the dilation parameter ( k = 6 ) .
Figure 8. Maximum absolute error versus the dilation parameter ( k = 6 ) .
Mathematics 09 01535 g008
Figure 9. The conventional RBFCM collocation scheme.
Figure 9. The conventional RBFCM collocation scheme.
Mathematics 09 01535 g009
Figure 10. The fictitious sources.
Figure 10. The fictitious sources.
Mathematics 09 01535 g010
Figure 11. The computed MAE values of the original PRBF and the proposed method: (a) p = 0.0001; (b) p = 1.
Figure 11. The computed MAE values of the original PRBF and the proposed method: (a) p = 0.0001; (b) p = 1.
Mathematics 09 01535 g011
Figure 12. Maximum absolute error versus the dilation parameter ( k = 7 ) .
Figure 12. Maximum absolute error versus the dilation parameter ( k = 7 ) .
Mathematics 09 01535 g012
Figure 13. The conventional RBFCM collocation scheme.
Figure 13. The conventional RBFCM collocation scheme.
Mathematics 09 01535 g013
Figure 14. The fictitious sources.
Figure 14. The fictitious sources.
Mathematics 09 01535 g014
Figure 15. The computed MAE values of the original PRBF and the proposed method: (a) α g = 0.0001 ; (b) α g = 1 .
Figure 15. The computed MAE values of the original PRBF and the proposed method: (a) α g = 0.0001 ; (b) α g = 1 .
Mathematics 09 01535 g015
Figure 16. Maximum absolute error versus the dilation parameter ( k = 4 ) .
Figure 16. Maximum absolute error versus the dilation parameter ( k = 4 ) .
Mathematics 09 01535 g016
Figure 17. The conventional RBFCM collocation scheme.
Figure 17. The conventional RBFCM collocation scheme.
Mathematics 09 01535 g017
Figure 18. The fictitious sources.
Figure 18. The fictitious sources.
Mathematics 09 01535 g018
Figure 19. The computed MAE values of the original PS and the proposed method.
Figure 19. The computed MAE values of the original PS and the proposed method.
Mathematics 09 01535 g019
Figure 20. Maximum absolute error versus the dilation parameter ( k = 10 ) .
Figure 20. Maximum absolute error versus the dilation parameter ( k = 10 ) .
Mathematics 09 01535 g020
Table 1. A comparison of the results.
Table 1. A comparison of the results.
Optimum
Order
Dilation
Parameter
MAERMSE
The conventional PRBF7- 5.87 × 10 4 8.21 × 10 6
The proposed methodCase I15 8.91 × 10 7 2.14 × 10 8
Case II15 2.51 × 10 10 5.04 × 10 12
Case III15 3.03 × 10 12 4.95 × 10 14
Case IV15 2.66 × 10 10 4.35 × 10 12
Case V15 5.13 × 10 9 8.98 × 10 11
Case VI15 6.65 × 10 8 1.48 × 10 9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, C.-Y.; Ku, C.-Y.; Hong, L.-D.; Hsu, S.-M. Infinitely Smooth Polyharmonic RBF Collocation Method for Numerical Solution of Elliptic PDEs. Mathematics 2021, 9, 1535. https://doi.org/10.3390/math9131535

AMA Style

Liu C-Y, Ku C-Y, Hong L-D, Hsu S-M. Infinitely Smooth Polyharmonic RBF Collocation Method for Numerical Solution of Elliptic PDEs. Mathematics. 2021; 9(13):1535. https://doi.org/10.3390/math9131535

Chicago/Turabian Style

Liu, Chih-Yu, Cheng-Yu Ku, Li-Dan Hong, and Shih-Meng Hsu. 2021. "Infinitely Smooth Polyharmonic RBF Collocation Method for Numerical Solution of Elliptic PDEs" Mathematics 9, no. 13: 1535. https://doi.org/10.3390/math9131535

APA Style

Liu, C. -Y., Ku, C. -Y., Hong, L. -D., & Hsu, S. -M. (2021). Infinitely Smooth Polyharmonic RBF Collocation Method for Numerical Solution of Elliptic PDEs. Mathematics, 9(13), 1535. https://doi.org/10.3390/math9131535

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop