# Determining the Optimal Source Points in MFS by Minimizing an Energy Gap Functional for 3D Laplace Operator

Authors: M. Sajjadmanesh; S. Shakib Khanghah; H. Aydi
Journal: Journal of Function Spaces
 (2020)

---

## Abstract

In this paper, an extended version of the method of minimizing an energy gap functional for determining the optimal source points in the method of fundamental solutions (MFS) is applied to the 3D Laplace operator subject to the Dirichlet and Neumann boundary conditions. As we know, the MFS is a more popular meshless method for solving boundary or initial-boundary value problems due to its simplicity and high accuracy. However, the accuracy of the MFS depends strongly on the distribution of the source points. Finally, some of the numerical experiments are carried out to express the simplicity and effectiveness of the presented method.


---

## Body

## 1. Introduction

In the last decades, the fixed point results have been improved and generalized in different directions for solving boundary value problems (see [1–5]), even the fixed point theorems have been extended for establishing the existence of solutions for fractional differential equations, as well as for integral equations (see [6, 7]). These results have been usually obtained by analytic techniques and various fixed point theorems. The method of a fundamental solution (MFS) is one of the meshless numerical methods, which was first proposed by Kupradze and Aleksidze [8] in 1964. This method approximates the solution of a boundary value problem (BVP) by a linear combination of fundamental solutions of the governing partial differential operator, provided that the fundamental solution is known. The MFS is a very popular boundary meshless method due to its simplicity and high accuracy. However, this method has a serious drawback that the resulting algebraic equations system may be highly ill-conditioned when the number of source points is increased [9], or when the distances of source points are increased [10]. The method of minimizing an energy gap functional is an effective and efficient way to determine the location of source points, which was first proposed by Wang et al. in [11] for a mixed boundary value problem of 2D Laplace equation. In this paper, we extend the same recent method for the 3D Laplace equation with Dirichlet and Neumann boundary conditions. Finally, several numerical examples are provided to illustrate the efficiency and simplicity of the method.This paper is organized as follows. In the next section, the MFS for the 3D Laplace operator is described. In Section3, the energy gap functional is constructed and finally, in Sections 4 and 5, the numerical algorithm and some of the numerical experiments, respectively, are carried out for proposed method.


## 2. Statement of Problem and the MFS

Consider the mixed boundary value problem of the 3D Laplace equation in spherical coordinates:(1)Δur,θ,φ=0,inΩ,(2)uρ,θ,φ=fθ,φ,onΓD,(3)∂u∂nρ,θ,φ=gθ,φ,onΓN,where Ω⊂ℝ3 is a bounded simply connected domain with the boundary ∂Ω=ΓD∪ΓN:
(4)ΓD=r=ρθ,φ,0≤θ≤απ,0≤φ≤π,ΓN=r=ρθ,φ,απ<θ<2π,0≤φ≤π.The functionρθ,φ is used to describe the boundary shape, the constant α∈0,2 is given, and n is the outward unit normal at ∂Ω.In the MFS, the numerical solution ofu at the field point z=rcosθsinφ,sinθsinφ,cosφ can be expressed as a linear combination of the fundamental solutions Uz,sj:
(5)uz=∑j=1nscjUz,sj,where the source points
(6)sj=Rjcosθjsinφj,sinθjsinφj,cosφj,1≤j≤nsuniformly located on the complementary set Ωc with
(7)Rj=D+ρθj,φj,1≤j≤ns,and D is an unknown offset to be determined such that the error of the numerical solution obtained by the MFS is minimized.The fundamental solutionUz,sj of 3D Laplace equation (1) is given by [12]
(8)Uz,sj=14πrj,rj=z−sj.By collectingnc points
(9)zi=ρθi,φicosθisinφi,sinθisinφi,cosφi,1≤i≤ncuniformly located on the boundary ∂Ω and satisfying the boundary conditions (2) and (3) and considering (5), a linear equation system can be achieved to determine the coefficients cj.


## 3. Energy Gap Functional

To determine the optimal source points, we have the following result which is an extension of Theorem1 in [11] to 3D Laplace operator.Theorem 1.
The energy gap functional with respect to the mixed BVP (1)-(3) is given by
(10)G=∬ΓDunρθ,φ,θ,φfθ,φdσ+∬ΓNuρθ,φ,θ,φgθ,φdσ−∭Ω∥∇ur,θ,φ∥2ds=0,where dσ=∣∂r/∂θ×∂r/∂φ∣dθdφ and ds=ρ2sinφdρdθdφ.Proof.
Multiplying both sides of Eq. (1) by u and integrating it over domain Ω, we can derive
(11)∭ΩuΔuds=0,which can be rewritten as
(12)∭Ω∇·u∇u−∥∇u∥2ds=0.Applying the Gaussian divergence theorem [13] on Eq. (12) yields that
(13)∭Ω∥∇u∥2ds=∬∂Ωu∇u·ndσ,where dσ=∣∂r/∂θ×∂r/∂φ∣dθdφ, and n is the outward unit normal at ∂Ω=ΓD∪ΓN. Finally, by using the boundary conditions (2) and (3) in Eq. (13), the energy gap functional (10) can be derived, and the proof is completed.The numerical solutionu of the problem (1)-(3) in MFS usually does not satisfy the energy gap in Eq. (10), i.e., G≠0. Then, we attempt to minimize the energy gap G by finding the optimal value of D, i.e.,
(14)minD∈D1,D2G,where D1,D1 is an interval of which we attempt to seek the optimal value of D.


## 4. A Numerical Algorithm

Suppose thatM is the number of subintervals of D1,D2 which we attempt to minimize the energy gap functional ∣G∣ in Eq. (10). Moreover, nc and ns are the number of collocation and source points, respectively. For L=0,1,⋯,M, we do the following steps:
(I)
Let(15)DL=D1+LD2−D1Mand compute
(16)sj=DL+ρθj,φjcosθjsinφj,sinθjsinφj,cosφj,Uzi,sj=14πxic−xjs2+yic−yjs2+zic−zjs2−1/21≤i≤nc,1≤j≤ns,where zi=xic,yic,zic and sj=xjs,yjs,zjs are, respectively, the collocation and source points uniformly disturbed on ∂Ω and ∂D with
(17)∂D=DL+ρθ,φcosθsinφ,sinθsinφ,cosφ,0≤θ<2π,0≤φ≤π(II)
Inserting the collocation and source pointszi,sj into Eq. (5) and enforcing it to satisfy the boundary conditions (2) and (3), the following linear equations system can be derived
(18)Ac=b,where c=c1,c2,⋯,cnsT, and the nc×ns matrix A is in general ill-conditionedThe conjugate gradient method (CGM) can be used to solve the overdetermined system (18) as follows:
(III)
Let(19)R0=1ns∑k=1ns∑j=1ncajk21/2,where ajk is the jkth component of the coefficient matrix ACompute(20)s:,k=R0N:,k,where N:,k denotes the norm of the kth column of matrix A and let
(21)B=AP,where P=diags denotes the diagonal matrix with components s1,s2,⋯,sns. Then, the following overdetermined algebraic system with well-conditioned coefficient matrix B can be derived:
(22)Bc1=b,where the all of rows or columns of matrix B have the same norm R0.
(IV)
Let(23)D=BTB,b1=BTbThen, the following linear equation system with the positive definite matrixD can be obtained:
(24)Dc1=b1.(V)
Solvec1 from system (24) by the CGM and denote it by c1L and let cL=s.c1L, where s is a column matrix s=s1,s2,⋯,snsT, and "." denotes the inner product(VI)
Compute∣G∣ by inserting Eq. (5) for u into Eq. (10)(VII)
Consider the criterion of stopping and if it is achieved, thenD=DL is the optimal offset and inserting D and c=cL, the optimal solution u can be obtained


## 5. Numerical Experiments

In this section, we give two examples of mixed boundary value problems to check the effectiveness of the presented method. The implementation of the algorithm is based on the MATLAB software. We denote the analytical and numerical solutions byuana and unum, respectively. Define the maximum error (ME) and the root-mean-square-error (RMSE) as follows:
(25)MEu=maxuanaxl−unumxl,l=1,2,⋯,nc,RMSEu=1nc∑l=1ncuanaxl−unumxl21/2.In Tables1 and 2, we list the minimum energy gap ∣G∣, ME, and RMSE for different values of nc with ns=50. The curves of errors ME and RMSE and also, the optimal value Dopt in D1,D2, from the minimization problem (14) with G given by Eq. (10) are given in Figures 1 and 2 to show the accuracy and efficiency of the proposed method.Example 2.
Consider the mixed boundary value problem (1)-(3) on the unit sphere
(26)Ω=r,θ,φ;0<r<1,0≤θ<2π,0≤φ≤π,with
(27)ΓD=1,θ,φ;0≤θ≤π,0≤φ≤π,ΓN=1,θ,φ;π<θ<2π,0≤φ≤π.
The exact solution is given byux1,x2,x3=1/2x12+x22−x32. Here, D1,D2=10,15 is an interval we attempt to search the optimal value of D, of which Dopt=13 is obtained as shown in Figure 1.Example 3.
Consider the mixed boundary value problem (1)-(3) on the bumpy sphere
(28)Ω=r,θ,φ;r=ρ1+16sinφ,0≤ρ<1,0≤θ<2π,0≤φ≤π,with
(29)ΓD=x∈∂Ω;0≤θx≤π,0≤φx≤π,ΓN=x∈∂Ω;π<θx<2π,0≤φx≤π.Table 1
Comparing the accuracy with differentnc for Example 2.
ncMERMSEMinimum energy gap509.2605×10−34.237×10−31.477×10−22008.1481×10−42.5906×10−41.6706×10−34507.3568×10−42.1993×10−41.6488×10−3Table 2
Comparing the accuracy with differentnc for Example 3.
ncMERMSEMinimum of energy gap501.3657×10−35.0442×10−41.634×10−32004.5631×10−41.5347×10−43.159×10−44503.0346×10−41.0369×10−41.864×10−4Figure 1
Determining the optimal offset in the MFS by minimizing the energy gap for Example2.Figure 2
Determining the optimal offset in the MFS by minimizing the energy gap for Example3.The exact solution is given byux1,x2,x3=x1x2+x1x3+x2x3. The optimal value of D in D1,D2=1,5 which is Dopt=3.4 as shown in Figure 2.


## 6. Conclusion

In this study, the mixed boundary value problem, which consists of determining the optimal source points in MFS, has been investigated by minimizing an energy gap functional for 3D Laplace operator subject to the Dirichlet and Neumann boundary conditions. The method was first proposed by Wang et al. in [11] for a mixed boundary value problem involving the 2D Laplace equation, and we extended it for the 3D case. Inserting the collocation and source points in MFS and enforcing the boundary conditions yield to a system of linear equations with general ill-conditioned coefficient matrix which can be solved by one of the regularization methods such as the conjugate gradient method (CGM). Finally, the obtained numerical solutions via the MFS can be placed in an energy gap functional until its minimum value is obtained for the appropriate source points. Two examples for the unit spherical and the bumpy spherical boundaries have been considered. The numerical results show that minimizing of an energy gap functional is a simple and efficient method for determining source points in the MFS.



---

*Source: 6633908-2020-12-21.xml*