Next Article in Journal
Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints
Next Article in Special Issue
Multistability of the Vibrating System of a Micro Resonator
Previous Article in Journal
Existence and Uniqueness of Mild Solution Where α ∈ (1,2) for Fuzzy Fractional Evolution Equations with Uncertainty
Previous Article in Special Issue
Recent Advances and Applications of Spiral Dynamics Optimization Algorithm: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain

1
Faculty of Basic Science, University of Bonab, Bonab P.O. Box 55517-61167, Iran
2
Institut Supérieur d’Informatique et des Techniques de Communication, Université de Sousse, Sousse 4000, Tunisia
3
China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
4
Department of Mathematics and Applied Mathematics, Sefako Makgatho Health Sciences University, Ga-Rankuwa 0204, South Africa
5
Department of Mathematics, Taiz University, Taiz 6803, Yemen
6
Research Institute for Natural Sciences, Hanyang University, Seoul 04763, Korea
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2022, 6(2), 66; https://doi.org/10.3390/fractalfract6020066
Submission received: 11 December 2021 / Revised: 8 January 2022 / Accepted: 9 January 2022 / Published: 27 January 2022
(This article belongs to the Special Issue Advances in Optimization and Nonlinear Analysis)

Abstract

:
In this paper, we are interested in an inverse geometric problem for the three-dimensional Laplace equation to recover an inner boundary of an annular domain. This work is based on the method of fundamental solutions (MFS) by imposing the boundary Cauchy data in a least-square sense and minimisation of the objective function. This approach can also be considered with noisy boundary Cauchy data. The simplicity and efficiency of this method is illustrated in several numerical examples.

1. Introduction

The inverse geometry problems, as an important subclass of inverse problems, can be subdivided into two subclasses, depending on the location of the unknown boundary. In the first kind, the portion of the outer boundary of the solution domain is unknown, whilst in the second kind, the inner boundary is unknown.
There are many methods for solving the inverse geometry problems, such as the boundary element regularisation method by Lesnic et al. [1], the method of fundamental solutions and moving pseudo-boundary method by Karageorghis et al. [2,3,4], the boundary function method by Wang et al. [5], the conjugate gradient method (CGM) and the boundary element technique by Huang et al. [6,7].
Bin-Mohsin and Lesnic in 2012 utilised the method of fundamental solutions (MFS) to the modified Helmholtz inverse geometry problem on an annular domain [8].
The purpose of this paper is to extend the aforementioned current approach to the three-dimensional Laplace equation based on the method of fundamental solutions. Finally, two examples are presented to show the simplicity and efficiency of this method.

2. Formulation of the Inverse Geometric Problem

Let D R 3 be a simply connected domain with an unknown boundary D which is compactly contained in a simply connected domain Ω R 3 with the boundary Ω .
Let us consider the following inverse problem:
Δ u = 0 in Ω \ D ¯ ,
subject to the boundary conditions,
u = f on Ω ,
u n = g on Ω ,
u = h on D ,
where f H 1 / 2 ( Ω ) and g H 1 / 2 ( Ω ) are given functions and n is an outward unit normal vector on Ω . Moreover, the function h H 1 / 2 ( D ) is given on the unknown boundary D . Without loss of generality, we can suppose that Ω is the unit disk B ( 0 ̲ ; 1 ) ; otherwise we can conformally map the exterior of the simply connected domain Ω onto the exterior of the unit disk.
The unknown boundary D can be expressed in spherical coordinates as
D = { r ( θ , φ ) ( cos θ sin φ , sin θ sin φ , cos φ ) ; θ [ 0 , 2 π ) , φ [ 0 , π ] }
where r ( θ , φ ) is a 2 π -periodic and π -periodic smooth function with respect to θ and φ , respectively, with values in the interval ( 0 , 1 ) .
The inverse problem we are concerned with is to determine geometrically the domain boundary D by utilising the method of fundamental solutions.

3. The Least-Square Problem Based on the MFS

In the classic MFS, the solution of a homogeneous linear partial differential equation (PDE) is approximated by a linear combination of the fundamental solutions with the set of sources located outside the problem domain and a set of points on the domain boundary. The linear combination coefficients are determined by collocation or, alternatively, with a least-squares fit of the boundary conditions.
Based on the MFS, one can approximate the solution of (1) by a linear combination of its fundamental solutions, which is given by [9]
U ( z , s ) = 1 4 π r ; r = z s .
i.e.,
u ( z ) = j = 1 n s c j U ( z , s j ) ,
where the collocation points z i and z i + M are uniformly located on Ω and D , respectively, i.e.,
z i = ( cos θ ^ i sin φ ^ i , sin θ ^ i sin φ ^ i , cos φ ^ i ) , i = 1 , M ¯
z i + M = r i ( cos θ i sin φ i , sin θ i sin φ i , cos φ i ) , i = 1 , N ¯
Further, the n s : = M + N source points s j and s j + M are uniformly located on the outside of Ω ¯ and the inside of D, respectively, i.e.,
s j = R 1 ( cos θ ^ j sin φ ^ j , sin θ ^ j sin φ ^ j , cos φ ^ j ) , j = 1 , M ¯
s j + M = r j R 2 ( cos θ j sin φ j , sin θ j sin φ j , cos φ j ) , j = 1 , N ¯
where R 1 , R 2 > 1 .
The coefficients vector c = ( c j ) j = 1 , M + N ¯ in linear combination (6) and also, the radial vector r = ( r j ) j = 1 , N ¯ can be determined by imposing the boundary conditions (2)–(4) in a least-square sense, which recasts into minimising the objective function
T ( c , r ) = u f L 2 ( Ω ) 2 + u n g L 2 ( Ω ) 2 + u h L 2 ( D ) 2 .
Upon discretisation, Equation (11) yields
T ( c , r ) = i = 1 M j = 1 M + N c j U ( z i , s j ) f ( z i ) 2 + i = M + 1 2 M j = 1 M + N c j U n ( z i M , s j ) g ( z i M ) 2 + i = 2 M + 1 2 M + N j = 1 M + N c j U ( z i M , s j ) h ( z i M ) 2 .
In general, the boundary data F { f , g , h } are measured noisy data satisfying
F i δ = F i + δ r a n d ( i ) F i ,
where δ is the percentage noise and the number r a n d ( i ) is a random number drawn from the standard uniform distribution on the interval [ 1 , 1 ] generated by the MATLAB code 1 + 2 r a n d ( i ) .
Imposing noise on all measured data implies
T δ ( c , r ) = i = 1 M j = 1 M + N c j U ( z i , s j ) f δ ( z i ) 2 + i = M + 1 2 M j = 1 M + N c j U n ( z i M , s j ) g δ ( z i M ) 2 + i = 2 M + 1 2 M + N j = 1 M + N c j U ( z i M , s j ) h δ ( z i M ) 2 .
The minimisation of (12) or (14) imposes 2 M + N nonlinear equations for the 2 N + M unknowns ( c , r ) , and for a unique solution, it is necessary that M N .

4. Error Analysis and the Regularisation

The accuracy of the presented method is evaluated by the normalised relative root mean square error (RMSE) and L -error:
RMSE = 1 N i = 1 N | r i ( a n ) r i ( n u m ) | 2 1 2 max 1 i N | r i ( a n ) | , L - error = max 1 i N | r i ( a n ) r i ( n u m ) |
where r i ( a n ) and r i ( n u m ) denote the analytical and numerical radial vectors, respectively, at the i t h collocation point on the boundary D .
The obtained numerical radial vectors from the presented method are unstable, especially when noise is added to the boundary data, and so the regularisation is needed. For this, we can add the following regularisation terms via standard zeroth- and first-order Tikhonov’s regularisation with parameters λ 1 , λ 2 0 to the functional (14), i.e.,
R e g ( a ̲ , r ̲ ) = j = 2 M + N + 1 3 M + 2 N λ 1 a j 2 M N 2 + j = 3 M + 2 N + 1 3 M + 3 N 1 λ 2 ( r j 3 M 2 N + 1 r j 3 M 2 N ) 2 ,

5. Numerical Examples

In this section, we give some examples to check the effectiveness of the presented method. We consider a three-dimensional annular domain with an outer boundary as the unit sphere Ω = B ( 0 ̲ ; 1 ) , R 1 = R 2 = 2 and M = N { 25 , 50 } in (7)–(10). Moreover, the percentage noise δ = 5 % is added to every measured boundary data.
The minimisation of functional (12) or (14) is carried out using the MATLAB optimisation toolbox routine lsqnonlin, which solves nonlinear least-squares problems.
Example 1.
Consider a three-dimensional annular domain with an unknown inner boundary D = B ( 0 ̲ ; r ( a n ) ) of radius r ( a n ) = 0.7 . The boundary data are given as follows:
u | Ω = f ( θ , φ ) = 1 2 sin 2 φ ( cos θ + sin θ ) + sin 2 θ sin 2 φ , u n | Ω = g ( θ , φ ) = sin 2 φ ( cos θ + sin θ ) + sin 2 θ sin 2 φ , u | D = h ( θ , φ ) = 49 200 sin 2 φ ( cos θ + sin θ ) + sin 2 θ sin 2 φ .
The exact solution for these input boundary data is u ( x ) = x 1 x 2 + x 1 x 3 + x 2 x 3 .
Table 1 gives the values of the objective functions and the corresponding errors obtained using the optimal initial guess r 0 and M = N { 25 , 50 } without using regularisation parameters. It can be seen that the values of the corresponding errors increase with the number of collocation points and so regularisation is needed.
In Table 2 and Table 3, we present the values of the objective functions and the corresponding errors with initial guess, r 0 , obtained using the regularisation parameters λ 1 , λ 2 { 0 , 10 6 , 10 3 , 10 1 } with M = N { 25 , 50 } and so, in Table 4, we give the minimal objective functions and the corresponding errors with initial guess r 0 .
Example 2.
Consider a three-dimensional annular domain with an unknown inner boundary of radius r ( a n ) = 1 4 ( 1 + cos θ sin 2 φ ) . The boundary data are given as follows:
u | Ω = f ( θ , φ ) = 3 sin 2 φ 2 , u n | Ω = g ( θ , φ ) = 6 sin 2 φ 4 , u | D = h ( θ , φ ) = 1 16 ( 3 sin 2 φ 2 ) ( sin 2 φ cos θ + 1 ) 2 .
The exact solution for these input boundary data is u ( x ) = x 1 2 + x 2 2 2 x 3 2 .
Table 5 gives the values of the objective functions and the corresponding errors obtained using the optimal initial guess r 0 , M = N { 25 , 50 } without using regularisation parameters, whilst Table 6 and Table 7 are obtained using the regularisation parameters λ 1 , λ 2 and so, in Table 8 we give the minimal objective functions and the corresponding errors with initial guess r 0 .

6. Conclusions

In this paper, we extended the aforementioned method presented in [8], based on the method of fundamental solutions to solve numerically the three-dimensional inverse geometry problem on an annular domain. To obtain the stable and accuracy results, Tikhonov’s regularisation parameters were used combined with the problem of the minimising an objective function. From the examples, we can see that our proposed method is effective and stable, even for the boundary data added with noise.

Author Contributions

Validation, H.A.; formal analysis, E.A., C.P.; investigation, M.S., H.A., E.A., C.P.; writing—original draft preparation, M.S.; writing—review and editing, H.A.; supervision, M.S.; funding acquisition, C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work does not receive any external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article, as no data set was generated or analysed during the current study.

Conflicts of Interest

The authors declare no confilct of interests.

References

  1. Lesnic, D.; Berger, J.R.; Martin, P.A. A boundary element regularization method for the boundary determination in potential corrosion damage. Inverse Probl. Eng. 2002, 10, 163–182. [Google Scholar] [CrossRef]
  2. Karageorghis, A.; Lesnic, D.; Marin, L. The method of fundamental solutions for the detection of rigid inclusions and cavities in plane linear elastic bodies. Comput. Struct. 2012, 106, 176–188. [Google Scholar] [CrossRef]
  3. Karageorghis, A.; Lesnic, D.; Marin, L. A moving pseudo-boundary method of fundamental solutions for void detection. Numer. Methods Partial Diff. Equat. 2013, 29, 935–960. [Google Scholar] [CrossRef] [Green Version]
  4. Karageorghis, A.; Lesnic, D.; Marin, L. The method of fundamental solutions for three-dimensional inverse geometric elasticity problems. Comput. Struct. 2016, 166, 51–59. [Google Scholar] [CrossRef] [Green Version]
  5. Wanga, F.; Huaa, Q.; Liub, C.S. Boundary function method for inverse geometry problem in two-dimensional anisotropic heat conduction equation. Appl. Math. Lett. 2018, 84, 130–136. [Google Scholar] [CrossRef]
  6. Huang, C.H.; Yan, J.Y. An inverse problem in simultaneously measuring temperature dependent thermal conductivity and heat capacity. Int. J. Heat Mass Transf. 1995, 38, 3433–3441. [Google Scholar] [CrossRef]
  7. Huang, C.H.; Yeh, C.Y.; Orlande, H.R.B. A non-linear inverse problem in simultaneously estimating the heat and mass production rates for a chemically reacting fluid. Chem. Eng. Sci. 2003, 58, 3741–3752. [Google Scholar] [CrossRef]
  8. Bin-Mohsin, B.; Lesnic, D. Determination of inner boundaries in modified Helmholtz inverse geometric problems using the method of fundamental solutions. Math. Comput. Simul. 2012, 82, 1445–1458. [Google Scholar] [CrossRef]
  9. Vladimirov, V.S. Equations of Mathematical Physics; Mir Publishers: Moscow, Russia, 1981. [Google Scholar]
Table 1. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors with M = N { 25 , 50 } and no regularisation parameters for Example 1.
Table 1. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors with M = N { 25 , 50 } and no regularisation parameters for Example 1.
M = N r 0 Objective Functions L -ErrorRMSE
Without noise25 0.5 ̲ 4.9013 × 10 3 4.38841.4418
50 0.6 ̲ 6.16825.01141.7985
With noise 5%25 0.6 ̲ 2.17992.1210 7.1903 × 10 1
50 0.6 ̲ 1.8326 × 10 1 1.0113 9.1243 × 10 1
Table 2. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 with M = N = 25 for Example 1.
Table 2. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 with M = N = 25 for Example 1.
r 0 Objective Functions L -ErrorRMSE
Without noise λ 1 = 0 λ 2 = 10 6 0.7 ̲ 1.8102 × 10 2 2.61001.0249
λ 2 = 10 3 0.5 ̲ 2.5783 × 10 3 1.2932 × 10 1 1.3835 × 10 1
λ 2 = 10 1 0.3 ̲ 2.2118 × 10 3 1.0442 × 10 1 1.3354 × 10 1
λ 1 = 10 6 λ 2 = 0 0.6 ̲ 5.8822 × 10 3 2.9436 × 10 1 1.5451 × 10 1
λ 1 = 10 3 0.5 ̲ 1.24702.8593 8.8216 × 10 1
λ 1 = 10 1 0.7 ̲ 5.8688 5.1408 × 10 3 1.4688 × 10 3
With noise 5% λ 1 = 0 λ 2 = 10 6 0.9 ̲ 4.1621 × 10 2 3.5650 × 10 1 1.1951 × 10 1
λ 2 = 10 3 0.6 ̲ 1.3166 × 10 2 1.6224 × 10 1 1.3880 × 10 1
λ 2 = 10 1 0.5 ̲ 1.2386 × 10 2 1.0866 × 10 1 1.1611 × 10 1
λ 1 = 10 6 λ 2 = 0 0.4 ̲ 3.1242 × 10 2 6.5845 × 10 1 2.8637 × 10 1
λ 1 = 10 3 0.2 ̲ 1.2775 8.0410 × 10 1 3.8164 × 10 1
λ 1 = 10 1 0.2 ̲ 5.98973.36462.0861
Table 3. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 with M = N = 50 for Example 1.
Table 3. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 with M = N = 50 for Example 1.
r 0 Objective Functions L -ErrorRMSE
Without noise λ 1 = 0 λ 2 = 10 6 0.6 ̲ 3.40924.69981.726
λ 2 = 10 3 0.7 ̲ 2.4650 8.8865 × 10 1 4.5447 × 10 1
λ 2 = 10 1 0.7 ̲ 1.6998 × 10 5 9.8441 × 10 3 1.1249 × 10 2
λ 1 = 10 6 λ 2 = 0 0.6 ̲ 4.3927 × 10 2 4.73631.8483
λ 1 = 10 3 0.5 ̲ 9.7102 × 10 1 4.65461.9311
λ 1 = 10 1 0.2 ̲ 1.4145 × 10 1 1.29831.5821
With noise 5% λ 1 = 0 λ 2 = 10 6 0.6 ̲ 1.26001.5259 4.6191 × 10 1
λ 2 = 10 3 0.8 ̲ 1.2761 × 10 2 1.0640 4.4252 × 10 1
λ 2 = 10 1 0.6 ̲ 6.2499 × 10 3 6.2367 × 10 2 3.1760 × 10 2
λ 1 = 10 6 λ 2 = 0 0.6 ̲ 1.16944.81641.6659
λ 1 = 10 3 0.7 ̲ 9.7398 × 10 1 3.20811.0049
λ 1 = 10 1 0.6 ̲ 1.0836 × 10 1 4.40801.9817
Table 4. The values of the minimal objective functions and the corresponding errors with initial guess r 0 , obtained (with/no) selecting the optimal regularisation parameters with M = N { 25 , 50 } for Example 1.
Table 4. The values of the minimal objective functions and the corresponding errors with initial guess r 0 , obtained (with/no) selecting the optimal regularisation parameters with M = N { 25 , 50 } for Example 1.
M = N With/No NoiseWith/No Regularisation λ 1 λ 2 r 0 Objective Function L -ErrorRMSE
25nono 0.5 ̲ 4.9013 × 10 3 4.38841.4418
with0        10 1 0.3 ̲ 2.2118 × 10 3 1.0442 × 10 1 1.3354 × 10 1
withno 0.6 ̲ 2.17992.1210 7.1903 × 10 1
with0        10 1 0.5 ̲ 1.2386 × 10 2 1.0866 × 10 1 1.1611 × 10 1
50nono 0.6 ̲ 6.16825.01141.7985
with0        10 1 0.7 ̲ 1.6998 × 10 5 9.8441 × 10 3 1.1249 × 10 2
withno 0.6 ̲ 1.8326 × 10 1 1.0113 9.1243 × 10 1
with0        10 1 0.6 ̲ 6.2499 × 10 3 6.2367 × 10 2 3.1760 × 10 2
Table 5. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors for Example 2 with M = N { 25 , 50 } and no regularisation parameters.
Table 5. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors for Example 2 with M = N { 25 , 50 } and no regularisation parameters.
M = N r 0 Objective Functions L -ErrorRMSE
Without noise25 0.5 ̲ 5.0884 × 10 1 4.1355 × 10 1 6.0305 × 10 1
50 0.6 ̲ 1.0837 × 10 2 8.1024 × 10 1 1.2263
With noise 5%25 0.8 ̲ 3.2489 × 10 1 1.92511.9305
50 0.7 ̲ 5.4166 × 10 1 4.51301.8355
Table 6. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 { 0 , 10 6 , 10 3 , 10 1 } with M = N = 25 for Example 2.
Table 6. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 { 0 , 10 6 , 10 3 , 10 1 } with M = N = 25 for Example 2.
r 0 Objective Functions L -ErrorRMSE
Without noise λ 1 = 0 λ 2 = 10 6 0.5 ̲ 5.3972 × 10 1 5.0754 × 10 1 6.3467 × 10 1
λ 2 = 10 3 0.4 ̲ 4.1825 × 10 1 4.0134 × 10 1 4.7861 × 10 1
λ 2 = 10 1 0.8 ̲ 2.0433 × 10 1 4.0078 × 10 1 4.5618 × 10 1
λ 1 = 10 6 λ 2 = 0 0.6 ̲ 4.0430 × 10 1 1.56341.3413
λ 1 = 10 3 0.1 ̲ 5.9128 4.9809 × 10 1 5.6395 × 10 1
λ 1 = 10 1 0.4 ̲ 2.5867 × 10 1 1.96642.9192
With noise 5% λ 1 = 0 λ 2 = 10 6 0.8 ̲ 3.3672 × 10 1 1.92511.9396
λ 2 = 10 3 0.2 ̲ 5.4494 5.3183 × 10 1 4.8568 × 10 1
λ 2 = 10 1 0.1 ̲ 1.9182 × 10 1 4.3568 × 10 1 4.7239 × 10 1
λ 1 = 10 6 λ 2 = 0 0.8 ̲ 3.42574.60552.1547
λ 1 = 10 3 0.2 ̲ 5.9130 6.4875 × 10 1 6.2476 × 10 1
λ 1 = 10 1 0.4 ̲ 2.5843 × 10 1 1.96542.9119
Table 7. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 { 0 , 10 6 , 10 3 , 10 1 } with M = N = 50 for Example 2.
Table 7. The values of the optimal initial guess, r 0 , objective functions and the corresponding errors using the regularisation parameters λ 1 , λ 2 { 0 , 10 6 , 10 3 , 10 1 } with M = N = 50 for Example 2.
r 0 Objective Functions L -ErrorRMSE
Without noise λ 1 = 0 λ 2 = 10 6 0.4 ̲ 8.3494 × 10 1 8.2159 × 10 1 1.2749
λ 2 = 10 3 0.2 ̲ 1.4682 × 10 4 8.0200 × 10 1 1.0192
λ 2 = 10 1 0.8 ̲ 3.8174 9.5145 × 10 1 1.3288
λ 1 = 10 6 λ 2 = 0 0.5 ̲ 2.07326.44082.9474
λ 1 = 10 3 0.7 ̲ 3.63525.21013.879
λ 1 = 10 1 0.4 ̲ 4.7401 × 10 1 1.1726 × 10 1 7.6963
With noise 5% λ 1 = 0 λ 2 = 10 6 0.8 ̲ 3.6080 × 10 1 8.5674 × 10 1 1.1744
λ 2 = 10 3 0.9 ̲ 2.1617 × 10 1 5.2427 × 10 1 6.5482 × 10 1
λ 2 = 10 1 0.7 ̲ 3.2234 × 10 2 4.7546 × 10 1 4.9726 × 10 1
λ 1 = 10 6 λ 2 = 0 0.6 ̲ 1.0459 × 10 1 6.7322 × 10 1 6.447 × 10 1
λ 1 = 10 3 0.6 ̲ 7.0388 9.6421 × 10 1 1.1551
λ 1 = 10 1 0.7 ̲ 7.5548 × 10 1 3.4477 × 10 3 1.2448 × 10 3
Table 8. The values of the minimal objective functions and the corresponding errors with initial guess, r 0 , obtained (with/no) selecting the optimal regularisation parameters λ 1 , λ 2 with M = N { 25 , 50 } for Example 2.
Table 8. The values of the minimal objective functions and the corresponding errors with initial guess, r 0 , obtained (with/no) selecting the optimal regularisation parameters λ 1 , λ 2 with M = N { 25 , 50 } for Example 2.
M = N With/No NoiseWith/No Regularisation λ 1 λ 2 r 0 Objective Function L -ErrorRMSE
25nono 0.5 ̲ 5.0884 × 10 1 4.1355 × 10 1 6.0305 × 10 1
with0        10 1 0.8 ̲ 2.0433 × 10 1 4.0078 × 10 1 4.5618 × 10 1
withno 0.8 ̲ 3.2489 × 10 1 1.92511.9305
with0        10 1 0.1 ̲ 1.9182 × 10 1 4.3568 × 10 1 4.7239 × 10 1
50nono 0.6 ̲ 1.0837 × 10 2 8.1024 × 10 1 1.2263
with0        10 3 0.2 ̲ 1.4682 × 10 4 8.0200 × 10 1 1.0192
withno 0.7 ̲ 5.4166 × 10 1 4.51301.8355
with0        10 1 0.7 ̲ 3.2234 × 10 2 4.7546 × 10 1 4.9726 × 10 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sajjadmanesh, M.; Aydi, H.; Ameer, E.; Park, C. The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain. Fractal Fract. 2022, 6, 66. https://doi.org/10.3390/fractalfract6020066

AMA Style

Sajjadmanesh M, Aydi H, Ameer E, Park C. The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain. Fractal and Fractional. 2022; 6(2):66. https://doi.org/10.3390/fractalfract6020066

Chicago/Turabian Style

Sajjadmanesh, Mojtaba, Hassen Aydi, Eskandar Ameer, and Choonkil Park. 2022. "The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain" Fractal and Fractional 6, no. 2: 66. https://doi.org/10.3390/fractalfract6020066

APA Style

Sajjadmanesh, M., Aydi, H., Ameer, E., & Park, C. (2022). The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain. Fractal and Fractional, 6(2), 66. https://doi.org/10.3390/fractalfract6020066

Article Metrics

Back to TopTop