Next Article in Journal
Reliability Evaluation for a Stochastic Flow Network Based on Upper and Lower Boundary Vectors
Previous Article in Journal
The Topological Transversality Theorem for Multivalued Maps with Continuous Selections
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Solving Modified Helmholtz Equation in Layered Materials Using the Multiple Source Meshfree Approach

1
Department of Harbor and River Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
2
Center of Excellence for Ocean Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(11), 1114; https://doi.org/10.3390/math7111114
Submission received: 6 November 2019 / Revised: 12 November 2019 / Accepted: 13 November 2019 / Published: 15 November 2019

Abstract

:
This paper presents a study for solving the modified Helmholtz equation in layered materials using the multiple source meshfree approach (MSMA). The key idea of the MSMA starts with the method of fundamental solutions (MFS) as well as the collocation Trefftz method (CTM). The multiple source collocation scheme in the MSMA stems from the MFS and the basis functions are formulated using the CTM. The solution of the modified Helmholtz equation is therefore approximated by the superposition theorem using particular nonsingular functions by means of multiple sources located within the domain. To deal with the two-dimensional modified Helmholtz equation in layered materials, the domain decomposition method was adopted. Numerical examples were carried out to validate the method. The results illustrate that the MSMA is relatively simple because it avoids a complicated procedure for finding the appropriate position of the sources. Additionally, the MSMA for solving the modified Helmholtz equation is advantageous because the source points can be collocated on or within the domain boundary and the results are not sensitive to the location of source points. Finally, compared with other methods, highly accurate solutions can be obtained using the proposed method.

1. Introduction

The solution of the modified Helmholtz equation plays a crucial role in the science and engineering fields, such as for boundary detection problems [1], water wave problems [2], Cauchy problems [3,4], diffusion equations [5], topological sensitivity analyses [6], and boundary value problems [7]. Over the past 10 years, conventional approaches have been widely adopted for solving the modified Helmholtz equation [8,9,10]. In contrast to the conventional approaches, several boundary collocation meshfree methods that have the basis functions satisfying the partial differential equation have also been proposed, for instance, the boundary knot method (BKM) [11,12], the boundary node method [13,14], the singular boundary method (SBM) [15,16,17], and the regularized meshless method [18,19]. In 2003, Chen and Hon [12] adopted the BKM to solve Helmholtz, modified Helmholtz, and the combination of diffusion and convection problems. The convergence of the BKM for the above equations was also investigated. Later, Chen et al. [15] applied the SBM to solve the modified Helmholtz equation. To examine the accuracy of the SBM, several wave numbers were considered in the numerical model [15]. For the modeling of Helmholtz-type equations in composite layered materials, Bin-Mohsin and Lesnic [20] adopted the method of fundamental solutions (MFS).
The MFS [21,22] is categorized into the indirect boundary element method or the boundary collocation meshfree approach, since only the boundary collocation points are used in the solution process. To avoid the singularity of the fundamental solutions, the MFS must collocate the source points outside the physical domain. Since it is a challenge to identify the location of source points, several attempts regarding this issue have been studied [23,24]. In 2007, Young et al. [25] proposed the modified MFS for finding the solution to the Laplace equation with an arbitrary domain. Later, Liu [26] adopted the multiple-length MFS (MLMFS) to investigate the appropriate placement of the sources for the Laplace equation. By using the MLMFS, the condition number can be reduced and the accuracy of the MLMFS may be increased. In recent years, Shigeta et al. [27] presented the multilayer method of fundamental solutions in conjunction with the weighted greedy QR decomposition (WGQRD) to solve the Laplace equation. As for the WGQRD, it can be used to find the necessary source points and unnecessary ones are removed in the numerical model.
Recently, Ku et al. [28] proposed a novel boundary-type meshfree method which is a hybrid of the MFS as well as the collocation Trefftz method (CTM). It adopts the collocation scheme from the MFS and the matrix form can be constructed in terms of superpositioning the nonsingular basis function from the CTM. This approach was then used for solving the Laplace-type subsurface flow problem. Xiao et al. [29] further applied the method to the problems of heat transfer in heterogeneous multilayer materials in two dimensions. Since the hybrid boundary-type meshless methods are only found to solve Laplace-type problems, it is especially interesting to explore this newly developed method for solving other non-Laplace problems. To the authors’ knowledge, pioneering work using the new boundary-type meshfree method to solve the modified Helmholtz equation in layered materials has not yet been proposed. We therefore propose a new boundary-type meshfree method using the basis functions based on the first-kind modified Bessel function capable of solving the two-dimensional modified Helmholtz equation over arbitrary geometries with doubly and multiply connected domains.
Starting from previous studies [28,29], we conducted pioneering work to solve the modified Helmholtz equation in layered materials. Since it is necessary to derive particular nonsingular basis functions for the two-dimensional modified Helmholtz equation, the specific method in this study is named the multiple source meshfree approach (MSMA). The key idea of the MSMA starts with the MFS and the CTM. The multiple source collocation scheme in the MSMA stems from the MFS, and the basis functions are formulated from the CTM. The solution of the modified Helmholtz equation is therefore approximated by the superposition theorem using particular nonsingular functions by means of multiple sources located within the domain. To deal with the two-dimensional modified Helmholtz equation in layered materials, the domain decomposition method (DDM) [30,31] was adopted. Numerical examples were carried out to validate the method.
The article is organized as follows: A brief introduction of the two-dimensional modified Helmholtz equation and the formulation of the MSMA considering the domain decomposition method are given in Section 2. In Section 3, several numerical examples are included to evaluate the accuracy of the proposed approach by comparing it with exact solutions. In Section 4, we discuss the research. Findings are summarized in Section 5.

2. The Methodology

We consider a two-dimensional domain for the problem of interest in 2 and assume that Ω is enclosed by a boundary Ω = Γ . The two-dimensional modified Helmholtz equation and the corresponding boundary conditions are stated as
( 2 λ 2 ) u ( x , y ) = 0 , ( x , y ) Ω .
The Dirichlet and Neumann conditions are
u ( x , y ) = g ( x , y )   on   Γ D
u ( x , y ) n = f ( x , y )   on   Γ f
in which 2 denotes the two-dimensional domain, 2 denotes the Laplacian, u may represent temperature (or acoustic pressure), Ω represents the domain of the problem, λ is the wave number, g and f are the Dirichlet and Neumann boundary data, Γ D and Γ f denote the boundaries where the Dirichlet and Neumann conditions are given, and n represents the outward normal direction.
The idea of this study was originally motivated by the MFS and the CTM, in which the collocation scheme is from the MFS and the basis functions are from the Trefftz formulation. To elaborate the proposed MSMA, the collocation scheme is introduced first. Figure 1 demonstrates the collocation points for the MFS, CTM, and MSMA bounded by simply and multiply connected domains. As illustrated in Figure 1a, a source point is located within the physical boundary enclosed by a simply connected domain in the CTM. Instead of using just a source point, the concept of the MFS is to approximate the solutions via a fundamental solution by means of locating the sources away from the physical boundary, as illustrated in Figure 1b. As for the MSMA, the solution is approximated through the linear combination of the particular nonsingular functions in terms of multiple sources collocated within the physical boundary, as illustrated in Figure 1c. For a multiply connected domain, which includes more than one cavity, multiple sources are placed within the physical boundary. The collocation scheme of the sources is demonstrated in Figure 1d–f for the CTM, MFS, and MSMA. From Figure 1, it is obvious that the proposed MSMA is advantageous because we may place the source and boundary points at the same position. Also, unlike the MFS, the source points can be placed at any position and are not sensitive to the results.
In the MSMA, the solution of the modified Helmholtz equation is approximated by the superposition theorem using nonsingular general solutions by means of multiple sources. Numerical solutions of the modified Helmholtz equation in two dimensions may be expressed as the combination of nonsingular general solutions. We may rewrite the two-dimensional modified Helmholtz equation in the polar coordinate as follows:
2 u r 2 + 1 r u r + 1 r 2 2 u θ 2 λ 2 u = 0
where r represents the radial coordinate and θ is the angular coordinate. For Equation (4), particular nonsingular functions can be derived using the separation of variables. Since the solutions are expressed as the linear superposition of the nonsingular functions, we may obtain the numerical solution of the two-dimensional modified Helmholtz equation represented by the linear superposition of the nonsingular basis functions:
u ( x ) j = 1 S k = 1 P a j , k H j , k ( x , y j )
in which x Ω , y j is the source point, a j , k denotes a vector of unknown coefficients, S denotes the source point number, and P represents the terms of the particular solutions. The particular solutions [4], H j , k ( x , y j ) , can be derived using the separation of variables for the two-dimensional modified Helmholtz equation in a simply connected domain given by
H j , k ( x , y j ) = [ I k ( λ ρ j ) sin ( k θ j ) I k ( λ ρ j ) cos ( k θ j ) ] T
where I k is the modified Bessel function of order k for the first kind; ρ j denotes the distance between the boundary collocation point and the source point, ρ j = | x y j | ; and θ j is the angle of the boundary point and the source point. For the Dirichlet boundary condition, we have
u ( x l ) j = 1 S k = 1 P a j , k H j , k ( x l , y j ) = g ( x l )
in which l = 1 , , Q , Q represents the boundary point number, and g ( x l ) represents the Dirichlet boundary data given on the boundary points. From Equation (3), the boundary of the Neumann condition is rewritten as
u ( x l ) n = ( u ( x l ) ρ ρ x + u ( x l ) θ θ x ) n x + ( u ( x l ) ρ ρ y + u ( x l ) θ θ y ) n y
where n x and n y at x l denote the outward normal direction vector. Using Equation (5), the derivatives of u / ρ and u / θ are given by
u ( x l ) ρ = j = 1 S k = 1 P a j , k [ λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) sin ( k θ j ) λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) cos ( k θ j ) ] T
u ( x l ) θ = j = 1 S k = 1 P a j , k [ k I k ( λ ρ j ) cos ( k θ j ) k I k ( λ ρ j ) sin ( k θ j ) ] T .
Substituting Equations (9) and (10) into Equation (8), the following equation is listed below:
u ( x l ) n = j = 1 S k = 1 P a j , k [ ( I k ( λ ρ j ) sin ( k θ j ) ) / n ( I k ( λ ρ j ) cos ( k θ j ) ) / n ] T
where
( I k ( λ ρ j ) sin ( k θ j ) ) n = j = 1 S k = 1 P [ λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) sin ( k θ j ) cos θ j k I k ( λ ρ j ) cos ( k θ j ) ( sin θ j ρ j ) ] n x + [ λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) sin ( k θ j ) sin θ j + k I k ( λ ρ j ) cos ( k θ j ) ( cos θ j ρ j ) ] n y ,
( I k ( λ ρ j ) cos ( k θ j ) ) n = j = 1 S k = 1 P [ λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) cos ( k θ j ) cos θ j + k I k ( λ ρ j ) sin ( k θ j ) ( sin θ j ρ j ) ] n x + [ λ I ( k + 1 ) ( λ ρ j ) + ( k ρ j ) I k ( λ ρ j ) cos ( k θ j ) sin θ j k I k ( λ ρ j ) sin ( k θ j ) ( cos θ j ρ j ) ] n y ,
Neumann boundary data are also imposed on the domain boundary using Equation (8):
u ( x l ) n j = 1 S k = 1 P a j , k [ ( I k ( λ ρ j ) sin ( k θ j ) ) / n ( I k ( λ ρ j ) cos ( k θ j ) ) / n ] T = f ( x l )
where l = 1 , , Q , and f ( x l ) represents the Neumann boundary data given on the boundary points. Using Equations (7) and (14), the following equation can be obtained as
A α = b
where A = [ I 1 ( λ ρ 1 ) sin ( θ 1 ) I 1 ( λ ρ 1 ) cos ( θ 1 ) I k ( λ ρ 1 ) sin ( k θ 1 ) I k ( λ ρ 1 ) cos ( k θ 1 ) I 1 ( λ ρ 2 ) sin ( θ 2 ) I 1 ( λ ρ 2 ) cos ( θ 2 ) I k ( λ ρ 2 ) sin ( k θ 2 ) I k ( λ ρ 2 ) cos ( k θ 2 ) I 1 ( λ ρ c ) sin ( θ c ) I 1 ( λ ρ c ) cos ( θ c ) I k ( λ ρ c ) sin ( k θ c ) I k ( λ ρ c ) cos ( k θ c ) N j , k = 1 I N j , k = 1 II N j , k = P I N j , k = P I I N j , k = 1 I N j , k = 1 I I N j , k = P I N j , k = P I I N j , k = 1 I N j , k = 1 I I N j , k = P I N j , k = P I I ] ,
α = [ a 1 a 2 a 3 a m 2 a m 1 a m ] T ,   b = [ g 1 g 2 g c f 1 f 2 f d ] T ,
and
N j , k I = j = 1 S k = 1 P [ λ I ( k + 1 ) ( λ ρ j ) + ( k / ρ j ) I k ( λ ρ j ) sin ( k θ j ) cos θ j k I k ( λ ρ j ) cos ( k θ j ) ( sin θ j / ρ j ) ] n x + [ λ I ( k + 1 ) ( λ ρ j ) + ( k / ρ j ) I k ( λ ρ j ) sin ( k θ j ) sin θ j + k I k ( λ ρ j ) cos ( k θ j ) ( cos θ j / ρ j ) ] n y
N j , k I I = j = 1 S k = 1 P [ λ I ( k + 1 ) ( λ ρ j ) + ( k / ρ j ) I k ( λ ρ j ) cos ( k θ j ) cos θ j + k I k ( λ ρ j ) sin ( k θ j ) ( sin θ j / ρ j ) ] n x + [ λ I ( k + 1 ) ( λ ρ j ) + ( k / ρ j ) I k ( λ ρ j ) cos ( k θ j ) sin θ j k I k ( λ ρ j ) sin ( k θ j ) ( cos θ j / ρ j ) ] n y .
In the preceding equations, A is a Q × S matrix, S represents the source point number, α represents a vector with the size of S × 1 for unknown coefficients, b denotes a vector with the size of Q × 1 for boundary values, c and d represent the boundary point number for Dirichlet and Neumann boundary data, respectively, a 1 , a 2 , , a m are unknown coefficients to be determined, and g 1 , g 2 , , g c and f 1 , f 2 , , f d denote the Dirichlet and Neumann boundary data, respectively.
To deal with the two-dimensional modified Helmholtz equation in layered materials, the DDM [30,31] was adopted. It should be noted that the DDM was also used to deal with the problems bounded by doubly and multiply connected domains in this study because the proposed MSMA approximates the solution of the two-dimensional modified Helmholtz equation as the combination of only particular nonsingular basis functions.
Unlike the CTM or the MFS, the MSMA does not consider the negative basis functions or the so-called singular basis functions in the proposed method. As a result, we are not able to tackle the singular solution in the cavity by using only the nonsingular basis functions. It is therefore necessary to split the domain into several simply connected subregions, such that we may guarantee subregions existing only in a nonsingular solution.
In the DDM, the physical boundary is split into several simply connected subregions which are connected at the interface between two consecutive subregions, as depicted in Figure 2. At the interface, the overspecified boundary data need to be satisfied. A doubly connected domain, Ω , is split into two subregions, Ω 1 and Ω 2 . For the subregions, we have Γ 1 , Γ 2 , …, and Γ 8 sub-boundaries, in which Ω 1 is composed of Γ 1 , Γ 2 , Γ 3 , and Γ 4 and Ω 2 is composed of Γ 5 , Γ 6 , Γ 7 , and Γ 8 . To connect Ω 1 and Ω 2 , we may use both the Dirichlet and Neumann boundary conditions to satisfy the continuity condition on the interface. The form of the continuity equations is described as follows:
u | Γ 1 = u | Γ 5 ,   u n | Γ 1 = u n | Γ 5 ,   u | Γ 3 = u | Γ 7 ,   u n | Γ 3 = u n | Γ 7 .
Using the DDM, we can obtain the matrix form, which is expressed as follows:
A D α D = b D
A D = [ A Ω 1 0 Ω 2 A I | Γ 1 , Γ 3 A I | Γ 5 , Γ 7 0 Ω 1 A Ω 2 ] ,   α D = [ α Ω 1 α Ω 2 ] ,   b D = [ b Ω 1 b I b Ω 2 ]
where A Ω 1 and A Ω 2 are the matrix with the size of m 1 × n 1 and m 2 × n 2 for Ω 1 and Ω 2 ; m 1 and m 2 are the boundary collocation points for Ω 1 and Ω 2 , respectively; and n 1 and n 2 are the terms of the particular solutions for Ω 1 and Ω 2 , respectively. At the interface, A I | Γ 1 , Γ 3 has the size of m I × n 1 for the boundary Γ 1 and Γ 3 , and A I | Γ 5 , Γ 7 has the size of m I × n 2 for the boundary Γ 5 and Γ 7 . m I represents the boundary point number at the interface. 0 Ω 1 and 0 Ω 2 are vectors for which all values are zero with the size of m 2 × n 1 and m 1 × n 2 , respectively. α Ω 1 denotes a n 1 × 1 vector of unknown coefficients of Ω 1 , α Ω 2 denotes a n 2 × 1 vector of unknown coefficients of Ω 2 . b I is the vector with the size of m I × 1 at the interface, b I = [ 0 g 0 f ] T , and 0g and 0f are vectors for which all values are zero. Solving the matrix form expressed in Equation (20), unknown coefficients b Ω 1 and b Ω 2 for subdomains Ω 1 and Ω 2 can then be determined.

3. Numerical Examples

To validate the MSMA for solving the modified Helmholtz equation, six examples were conducted to verify the proposed meshfree approach. The absolute error and relative root-mean-square error were used to evaluate the accuracy of the numerical solution as follows:
Absolute   error = | u a n a ( x i ) u n u ( x i ) | ,
Relative   root - mean - square   error = 1 n t i = 1 n t ( u n u ( x i ) u a n a ( x i ) ) 2 / 1 n t i = 1 n t u a n a 2 ( x i )
where n t is the number of inner points in the domain, and u n u ( x i ) and u n u ( x i ) are the analytical and numerical solutions at x i , respectively.

3.1. Modeling of Modified Helmholtz Equation Bounded by a Simply Connected Region

This example is the modeling of the modified Helmholtz equation bounded by a simply connected region with a gear-shaped domain, as depicted in Figure 3. The shape parameter of the physical boundary is given by
Ω = { ( x , y ) | x = ρ ( θ ) cos ( σ ( θ ) ) , y = ρ ( θ ) sin ( σ ( θ ) ) } ,
where ρ ( θ ) = 0.2 ( 2 + 0.5 sin ( 7 θ ) ) , σ ( θ ) = θ + 0.5 sin ( 7 θ ) , and 0 θ 2 π .
We consider the following exact solution:
u = e λ / 2 ( x + y ) .
In this case, the Dirichlet and Neumann boundary data are imposed using the exact solution. The wave number λ is set as 2 . Figure 3 illustrates the maximum absolute error (MAE) versus the terms of the particular solutions, k . The MSMA can obtain highly accurate results in the order of 10 11 while k = 15 . Figure 4 demonstrates the relationship of the MAE with the source point number. When the source points number is more than 100, the proposed MSMA can obtain accurate results with the order of 10 10 . As shown in Figure 5, we found that the accuracy of the results decreased with the increase of the wave number λ . To obtain the MAE, 1211 inner points were uniformly located within the domain. Figure 6 shows the comparison of the numerical and analytical solutions. The results demonstrate that the solutions of the MSMA agree very well with the exact solution.

3.2. Accuracy Comparison of the Proposed Method

The second case is the modeling of the modified Helmholtz equation bounded by an amoeba-shaped boundary, as illustrated in Figure 7. In this example, the results computed by MSMA were compared with those from the SBM [15]. The shape parameter of the irregular physical boundary is described by
Ω = { ( x , y ) | x = ρ ( θ ) cos ( θ ) , y = ρ ( θ ) sin ( θ ) } ,
in which ρ = e ( sin θ sin 2 θ ) 2 + e ( cos θ cos 2 θ ) 2 , and 0 θ 2 π .
Three different exact solutions can be found as follows:
u = e λ x ,
u = e ( 2 λ / 2 ) ( x y ) ,
u = e ( λ / 2 ) ( x + 3 y ) .
In the numerical implementation, the terms of the particular solutions, k , was 15. The Dirichlet boundary data were imposed on the domain boundary using three different analytical solutions, as depicted in Equations (27)–(29). The wave number λ was set as π . The computed results using the MSMA were compared with the SBM [15]. We obtained highly accurate numerical solutions while the boundary and source point number was greater than 100, as shown in Figure 8a. Compared with the results using the SBM [15] as shown in Figure 8b, we found that the accuracy of the MSMA was much more accurate than the SBM. Figure 9a shows the accuracy for the relative root-mean-square error versus the wave number. It was found that the accuracy of the increased wave number may decrease the accuracy of the computed results. Figure 9b shows the relative root-mean-square error versus the wave number computed by the SBM [15]. Figure 9 shows that the accuracy of the MSMA is better than that from the SBM.

3.3. Investigation of the Wave Number

The third case is the modeling of the modified Helmholtz equation bounded by a peanut-shaped boundary, as depicted in Figure 10. To examine the effect of the wave number for the numerical results, different values of the wave number were considered in the example. The shape parameter of the irregular physical boundary is described by
Ω = { ( x , y ) | x = ρ ( θ ) cos ( θ ) , y = ρ ( θ ) sin ( θ ) } ,
where ρ = 0.1 cos ( 2 θ ) + 1.1 si n 2 ( 2 θ ) , and 0 θ 2 π .
To solve the modified Helmholtz equation, the exact solution can be found as
u = e ( 2 λ / 2 ) ( x y ) .
The Dirichlet boundary values were imposed on the domain boundary using the exact solution. We investigated the accuracy of the proposed method by using different values of the wave number ranging from 1 to 50. Figure 10 shows the accuracy for the MAE versus the wave number using the proposed meshfree approach. As illustrated in Figure 10, the accuracy can be found in the order of 10 9 while λ = 50 . The numerical solutions revealed that the MSMA may yield highly accurate results even when the high wave number is considered.

3.4. Solution of the Modified Helmholtz Equation Bounded by a Doubly Connected Region

The fourth case is the modeling of the modified Helmholtz equation bounded by a doubly connected region, as depicted in Figure 11. The shape parameter of the concentric circles is described by
Ω 1 = { ( x , y ) | x = ρ 1 cos θ 1 , y = ρ 1 sin θ 1 } ,   ρ 1 = 1 ,
Ω 2 = { ( x , y ) | x = ρ 2 cos θ 2 , y = ρ 2 sin θ 2 } ,   ρ 2 = 0.5 .
The exact solution of this example can be found as follows:
u = I 3 ( k ρ ) cos ( 3 θ ) + K 3 ( k ρ ) cos ( 3 θ )
where I 3 is the modified Bessel function of the third order for the first kind, and K 3 is the modified Bessel function of the third order for the second kind. In this example, the terms of the particular solutions were set to be 15. The wave number was set as 2 . For this example, the physical boundary, enclosed by a doubly connected region, was split into two simply connected subregions using the DDM [30,31], as depicted in Figure 11. Four hundred boundary and source points were interactively arranged on the physical boundary. Figure 12 shows that the sub-boundaries at the interface contained four sub-boundaries: Γ 1 , Γ 3 , Γ 5 , and Γ 7 . The overspecified boundary conditions enforced on the boundary and source points at the interface are described by
u | Γ 1 = u | Γ 5 ,   u n | Γ 1 = u n | Γ 5 ,   u | Γ 3 = u | Γ 7 ,   u n | Γ 3 = u n | Γ 7 .
To examine the accuracy, 2415 inner points were placed within the domain to obtain the numerical results. We also compared the numerical solutions with the above exact solution, as illustrated in Figure 13. It seems that the numerical solutions perfectly agree with the exact solution. Furthermore, the MAE of the MSMA can be obtained in the order of 10 6 .

3.5. Solution of Modified Helmholtz Equation Bounded by a Multiply Connected Region

The fifth case is the modified Helmholtz equation bounded by a multiply connected region, as depicted in Figure 14. For the multiply connected domain including more than one cavity, we must place multiple source points within the physical boundary. Because the proposed MSMA does not consider the negative basis functions or the so-called singular basis functions, we cannot tackle the singular solution in the cavity by using only the nonsingular basis functions. Accordingly, to deal with problems bounded by multiply connected domains, the DDM was adopted, which split the domain into several simply connected subregions, such that we may guarantee subregions existing only in a nonsingular solution. The shape parameter of the physical boundary is given by
Ω 1 = { ( x , y ) | x = ρ 1 cos θ 1 , y = ρ 1 sin θ 1 } ,   ρ 1 = ( cos ( 3 θ 1 ) + 2 ( sin ( 3 θ 1 ) ) 2 ) 1 / 3 , 0 θ 1 2 π ,
Ω 2 = { ( x , y ) | x = ρ 2 cos θ 2 , y = ρ 2 sin θ 2 } ,   ρ 2 = e ( sin θ 2 sin 2 θ 2 ) 2 + e ( cos θ 2 cos 2 θ 2 ) 2 ,   0 θ 2 2 π ,
Ω 3 = { ( x , y ) | x = ρ 3 cos θ 3 , y = ρ 3 sin θ 3 } ,   ρ 3 = cos ( 2 θ 3 ) + 1.1 ( sin ( 2 θ 3 ) ) 2 , 0 θ 3 2 π .
The exact solution can be found as follows:
u = I 1 ( k ρ ) cos θ + K 1 ( k ρ ) cos θ ,
where I 1 is the modified Bessel function of the first order for the first kind, and K 1 is the modified Bessel function of the first order for the second kind. In this example, the term of the particular solutions was 15. The wave number was set as 2 . For this example, the physical boundary, enclosed by a doubly connected region, was split into two simply connected subregions using the DDM [30,31], as depicted in Figure 14. Four hundred boundary and source points were interactively arranged on the physical boundary. Figure 15 shows that the sub-boundaries at the interface contained six sub-boundaries: Γ 2 , Γ 4 , Γ 6 , Γ 8 , Γ 10 , and Γ 12 . The overspecified boundary conditions enforced on the boundary and source points at the interface are described by
u | Γ 2 = u | Γ 8 ,   u n | Γ 2 = u n | Γ 8 ,   u | Γ 4 = u | Γ 10 ,   u n | Γ 4 = u n | Γ 10 ,   u | Γ 6 = u | Γ 12 ,   u n | Γ 6 = u n | Γ 12 .
To examine the accuracy, 2737 inner points were placed within the domain to obtain the numerical results. We also compared the numerical solutions with the exact solution, as illustrated in Figure 16. The results demonstrate that the numerical solutions perfectly agree with the exact solution. Furthermore, the MAE of the MSMA can be obtained in the order of 10 6 .

3.6. Solution of Modified Helmholtz Equation in Two Layered Materials

The last example is the modeling of the modified Helmholtz equation in two layered materials, as depicted in Figure 17. The wave number in two subdomains were set to be λ 1 and λ 2 , in which λ 1 = 2 and λ 2 = 1 , respectively. The shape parameter of the boundary is described by
Ω 1 = { ( x 1 , y 1 ) | 0 x 1 1 , 0 y 1 0.5 } ,
Ω 2 = { ( x 2 , y 2 ) | 0 x 2 1 , 0 y 2 0.5 } .
The exact solutions are as follows:
u 1 = e x + y ,
u 2 = e x ( 1 + y ) .
In this example, the terms of the particular solutions were 10. For modeling the modified Helmholtz equation in two layered materials, the physical boundary was split into two simply connected subregions using the DDM [30,31], as illustrated in Figure 17. For each subregion, 600 boundary and source points were interactively arranged on the physical boundary. Figure 17 also shows that the sub-boundaries at the interface contained two sub-boundaries: Γ 1 and Γ 5 . The overspecified boundary conditions enforced on the boundary and source points at the interface are described by
u | Γ 1 = u | Γ 5 ,   u n | Γ 1 = u n | Γ 5 .
To evaluate the accuracy of the problem, 900 inner points were arranged within the domain to obtain the numerical results. We also compared the numerical solutions with the exact solution, as depicted in Figure 18. The results demonstrate that our results agree well with the above exact solution. Furthermore, we found that the MAE of the MSMA can be obtained in the order of 10 7 , as depicted in Figure 19.

4. Discussion

In this paper, we proposed a multiple source meshfree approach for solving a two-dimensional modified Helmholtz equation in layered materials. The novel conception of the proposed MSMA started with the MFS and the CTM. The advantages of the proposed method are as follows.
The proposed method for solving the modified Helmholtz equation avoids a complicated procedure for finding the location of the sources. By adopting the nonsingular functions, the source points can be collocated on or within the domain boundary. Compared with other methods, highly accurate solutions can be obtained regardless of the location of source points. We also demonstrated that the proposed method incorporated with the DDM can be used to solve the modified Helmholtz equation in layered materials as well as the domain bounded by simply and multiply connected regions.
We investigated the accuracy of the proposed method by using different values of the wave number ranging from 1 to 50 for the modified Helmholtz equation. The accuracy of the MSMA was found in the order of 10 9 , while the wave number was 50. The computed results revealed that the MSMA may yield highly accurate results even when a high wave number is considered. The proposed method is, however, limited for solving linear governing equations due to the approximation of the solution using the addition theorem.

5. Conclusions

In this study, we presented a meshfree approach for solving the modified Helmholtz equation bounded by simply, doubly, and multiply connected regions using the MSMA in two dimensions. The modified Helmholtz equation in two layered materials was also considered. The method was verified and numerical examples were also performed. We may point out the two major achievements of the paper:
  • The key idea of the MSMA stems from the indirect boundary element method, which adopts multiple source points. Because of the adoption of nonsingular functions, the sources can be collocated on or within the domain boundary without using a complicated searching algorithm to find the appropriate location of the source points.
  • To the best of our knowledge, the MSMA using nonsingular basis functions is newly developed. A pioneering work for solving the modified Helmholtz equation bounded by a multiply connected region was conducted using the MSMA in this study. Furthermore, the MAE of the proposed approach for the modified Helmholtz equation in two layered materials can reach up to the order of 10 7 . From the computed results, we conclude that the MSMA is relatively simple because it avoids a complicated procedure for finding the appropriate location of the source points. Moreover, the MSMA has advantages of highly accurate and boundary collocation only for solving problems with complex geometry.

Author Contributions

Conceptualization, C.-Y.K. and W.Y.; Data curation, C.-Y.L.; Methodology, C.-Y.K. and J.-E.X.; Validation, J.-E.X.; Writing—original manuscript, C.-Y.K.; Writing—review and editing, C.-Y.K. and J.-E.X.

Funding

This research was funded by the Ministry of Science and Technology of R.O.C. (grant number MOST 108-2119-M-019-001).

Acknowledgments

The authors appreciate the editor and the referees for invaluable comments that have enhanced the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, C.M.; Liu, Y.C.; Chan, H.F.; Hsiao, S.S. Solutions of boundary detection problem for modified Helmholtz equation by Trefftz method. Inverse Probl. Sci. Eng. 2014, 22, 267–281. [Google Scholar] [CrossRef]
  2. Chen, K.H.; Chen, J.T.; Chou, C.R.; Yueh, C.Y. Dual boundary element analysis of oblique incident wave passing a thin submerged breakwater. Eng. Anal. Bound. Elem. 2002, 26, 917–928. [Google Scholar] [CrossRef]
  3. Nguyen, H.T.; Tran, Q.V.; Nguyen, V.T. Some remarks on a modified Helmholtz equation with inhomogeneous source. Appl. Math. Model. 2013, 37, 793–814. [Google Scholar] [CrossRef]
  4. Liu, C.S.; Qu, W.; Chen, W.; Lin, J. A novel Trefftz method of the inverse Cauchy problem for 3D modified Helmholtz equation. Inverse Probl. Sci. Eng. 2016, 25, 1278–1298. [Google Scholar] [CrossRef]
  5. Ben-Avraham, D.; Fokas, A.S. The Solution of the modified Helmholtz equation in a wedge and an application to diffusion-limited coalescence. Phys. Lett. A 1999, 263, 355–359. [Google Scholar] [CrossRef]
  6. Jleli, M.; Samet, B.; Vial, G. Topological sensitivity analysis for the modified Helmholtz equation under an impedance condition on the boundary of a hole. J. Math. Pures Appl. 2014, 103, 557–574. [Google Scholar] [CrossRef]
  7. Li, X. On solving boundary value problems of modified Helmholtz equations by plane wave functions. J. Comput. Appl. Math. 2006, 195, 66–82. [Google Scholar] [CrossRef]
  8. Yoneta, A.; Tsuchimoto, M.; Honma, T. An analysis of axisymmetric modified Helmholtz equation by using boundary element method. IEEE Trans. Magn. 1990, 26, 1015–1018. [Google Scholar] [CrossRef]
  9. Marina, L.; Lesnic, D.; Mantič, V. Treatment of singularities in Helmholtz-type equations using the boundary element method. J. Sound Vib. 2004, 278, 39–62. [Google Scholar] [CrossRef]
  10. Wang, R.; Chen, X. A fast solver for boundary integral equations of the modified Helmholtz equation. J. Sci. Comput. 2014, 65, 553–575. [Google Scholar] [CrossRef]
  11. Shen, D.J.; Lin, J.; Chen, W. Boundary knot method solution of Helmholtz problems with boundary singularities. J. Mar. Sci. Technol. 2014, 22, 440–449. [Google Scholar]
  12. Chen, W.; Hon, Y.C. Numerical investigation on convergence of boundary knot method in the analysis of homogeneous Helmholtz, modified Helmholtz, and convection-diffusion problems. Comput. Methods Appl. Mech. Eng. 2014, 192, 1859–1875. [Google Scholar] [CrossRef]
  13. Song, R.; Chen, W. A hybrid boundary node method. Int. J. Numer. Meth. Eng. 2002, 53, 751–763. [Google Scholar]
  14. Mukherjee, Y.X.; Mukherjee, S. A boundary node method for three-dimensional linear elasticity. Int. J. Numer. Meth. Eng. 1999, 46, 1163–1184. [Google Scholar]
  15. Chen, W.; Zhang, J.Y.; Fu, Z.J. Singular boundary method for modified Helmholtz equations. Eng. Anal. Bound. Elem. 2014, 44, 112–119. [Google Scholar] [CrossRef]
  16. Junpu, L.; Chen, W. A modified singular boundary method for three-dimensional high frequency acoustic wave problems. Appl. Math. Model. 2018, 54, 189–201. [Google Scholar]
  17. Ragusa, M.A. Cauchy-Dirichlet problem associated to divergence form parabolic equations. Commun. Contemp. Math. 2004, 6, 377–393. [Google Scholar] [CrossRef]
  18. Sun, L.; Chen, W.; Zhang, C. A new formulation of regularized meshless method applied to interior and exterior anisotropic potential problems. Appl. Math. Model. 2013, 37, 7452–7464. [Google Scholar] [CrossRef]
  19. Liu, L. Single layer regularized meshless method for three dimensional exterior acoustic problem. Eng. Anal. Bound. Elem. 2017, 77, 138–144. [Google Scholar] [CrossRef]
  20. Bin-Mohsin, B.; Lesnic, D. The method of fundamental solutions for Helmholtz-type equations in composite materials. Comput. Math. Appl. 2011, 62, 4377–4390. [Google Scholar] [CrossRef]
  21. Liu, C.S. Improving the ill-conditioning of the method of fundamental solutions for 2D Laplace equation. CMES-Comp. Model. Eng. 2008, 28, 77–93. [Google Scholar]
  22. Li, Z.C.; Huang, H.T.; Lee, M.G.; Chiang, J.Y. Error analysis of the method of fundamental solutions for linear elastostatics. J. Comput. Appl. Math. 2013, 251, 133–153. [Google Scholar] [CrossRef]
  23. Nishimura, R.; Nishimori, K.; Ishihara, N. Determining the arrangement of fictitious charges in charge simulation method using genetic algorithms. J. Electrost. 2000, 49, 95–105. [Google Scholar] [CrossRef]
  24. Karageorghis, A. A practical algorithm for determining the optimal pseudo-boundary in the method of fundamental solutions. Adv. Appl. Math. Mech. 2009, 1, 510–528. [Google Scholar] [CrossRef]
  25. Young, D.L.; Chen, K.H.; Chen, J.T.; Kao, J.H. A modified method of fundamental solutions with source on the boundary for solving Laplace equations with circular and arbitrary domains. CMES-Comp. Model. Eng. 2007, 19, 197–221. [Google Scholar]
  26. Liu, C.S. An equilibrated method of fundamental solutions to choose the best source points for the Laplace equation. Eng. Anal. Bound. Elem. 2012, 36, 1235–1245. [Google Scholar] [CrossRef]
  27. Shigeta, T.; Young, D.L.; Liu, C.S. Adaptive multilayer method of fundamental solutions using a weighted greedy QR decomposition for the Laplace equation. J. Comput. Phys. 2012, 231, 7118–7132. [Google Scholar] [CrossRef]
  28. Ku, C.Y.; Xiao, J.E.; Liu, C.Y.; Fan, C.M. On modeling subsurface flow using a novel hybrid Trefftz-MFS method. Eng. Anal. Bound. Elem. 2019, 100, 225–236. [Google Scholar] [CrossRef]
  29. Xiao, J.E.; Ku, C.Y.; Huang, W.P.; Su, Y.; Tsai, Y.H. A novel hybrid boundary-type meshless method for solving heat conduction problems in layered materials. Appl. Sci. 2018, 8, 1887. [Google Scholar] [CrossRef]
  30. Boubendir, Y.; Midura, D. Non-overlapping domain decomposition algorithm based on modified transmission conditions for the Helmholtz equation. Comput. Math. Appl. 2018, 75, 1900–1911. [Google Scholar] [CrossRef]
  31. Meric, R.A. Domain decomposition methods for Laplace’s equation by the BEM. Commun. Numer. Meth. Eng. 2000, 16, 545–557. [Google Scholar] [CrossRef]
Figure 1. Illustration of the collocation points for the collocation Trefftz method (CTM), method of fundamental solutions (MFS), and multiple source meshfree approach (MSMA).
Figure 1. Illustration of the collocation points for the collocation Trefftz method (CTM), method of fundamental solutions (MFS), and multiple source meshfree approach (MSMA).
Mathematics 07 01114 g001
Figure 2. Diagram of the domain decomposition method (DDM) in a doubly connected region.
Figure 2. Diagram of the domain decomposition method (DDM) in a doubly connected region.
Mathematics 07 01114 g002
Figure 3. Maximum absolute error (MAE) versus the order of k .
Figure 3. Maximum absolute error (MAE) versus the order of k .
Mathematics 07 01114 g003
Figure 4. MAE versus the source point number.
Figure 4. MAE versus the source point number.
Mathematics 07 01114 g004
Figure 5. MAE versus the wave number λ .
Figure 5. MAE versus the wave number λ .
Mathematics 07 01114 g005
Figure 6. Results comparison with the exact solution.
Figure 6. Results comparison with the exact solution.
Mathematics 07 01114 g006
Figure 7. Boundary and source points for an amoeba-shaped boundary.
Figure 7. Boundary and source points for an amoeba-shaped boundary.
Mathematics 07 01114 g007
Figure 8. Comparison of the accuracy with the MSMA and the singular boundary method (SBM): (a) the MSMA and (b) the SBM.
Figure 8. Comparison of the accuracy with the MSMA and the singular boundary method (SBM): (a) the MSMA and (b) the SBM.
Mathematics 07 01114 g008
Figure 9. Accuracy comparison for the relative root-mean-square error versus the wave number: (a) the MSMA and (b) the SBM.
Figure 9. Accuracy comparison for the relative root-mean-square error versus the wave number: (a) the MSMA and (b) the SBM.
Mathematics 07 01114 g009
Figure 10. MAE versus the wave number λ .
Figure 10. MAE versus the wave number λ .
Mathematics 07 01114 g010
Figure 11. Collocation points for a doubly connected region.
Figure 11. Collocation points for a doubly connected region.
Mathematics 07 01114 g011
Figure 12. Configuration of sub-boundaries.
Figure 12. Configuration of sub-boundaries.
Mathematics 07 01114 g012
Figure 13. Results comparison with the exact solution.
Figure 13. Results comparison with the exact solution.
Mathematics 07 01114 g013
Figure 14. Collocation points for the example bounded by a multiply connected domain.
Figure 14. Collocation points for the example bounded by a multiply connected domain.
Mathematics 07 01114 g014
Figure 15. Configuration of sub-boundaries.
Figure 15. Configuration of sub-boundaries.
Mathematics 07 01114 g015
Figure 16. Results comparison with the exact solution.
Figure 16. Results comparison with the exact solution.
Mathematics 07 01114 g016
Figure 17. Boundary and source points for the example with two layered materials.
Figure 17. Boundary and source points for the example with two layered materials.
Mathematics 07 01114 g017
Figure 18. Results comparison with the exact solution.
Figure 18. Results comparison with the exact solution.
Mathematics 07 01114 g018
Figure 19. The absolute error.
Figure 19. The absolute error.
Mathematics 07 01114 g019

Share and Cite

MDPI and ACS Style

Ku, C.-Y.; Xiao, J.-E.; Yeih, W.; Liu, C.-Y. On Solving Modified Helmholtz Equation in Layered Materials Using the Multiple Source Meshfree Approach. Mathematics 2019, 7, 1114. https://doi.org/10.3390/math7111114

AMA Style

Ku C-Y, Xiao J-E, Yeih W, Liu C-Y. On Solving Modified Helmholtz Equation in Layered Materials Using the Multiple Source Meshfree Approach. Mathematics. 2019; 7(11):1114. https://doi.org/10.3390/math7111114

Chicago/Turabian Style

Ku, Cheng-Yu, Jing-En Xiao, Weichung Yeih, and Chih-Yu Liu. 2019. "On Solving Modified Helmholtz Equation in Layered Materials Using the Multiple Source Meshfree Approach" Mathematics 7, no. 11: 1114. https://doi.org/10.3390/math7111114

APA Style

Ku, C. -Y., Xiao, J. -E., Yeih, W., & Liu, C. -Y. (2019). On Solving Modified Helmholtz Equation in Layered Materials Using the Multiple Source Meshfree Approach. Mathematics, 7(11), 1114. https://doi.org/10.3390/math7111114

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