Takemi Shigeta

Takemi Shigeta, D. L. Young

The purpose of this study is to propose a high-accuracy and fast numerical method for the Cauchy problem of the Laplace equation. Our problem is directly discretized by the method of fundamental solutions (MFS). The Tikhonov regularization method stabilizes a numerical solution of the problem for given Cauchy data with high noises. The accuracy of the numerical solution depends on a regularization parameter of the Tikhonov regularization technique and some parameters of MFS. The L-curve determines a suitable regularization parameter for obtaining an accurate solution. Numerical experiments show that such a suitable regularization parameter coincides with the optimal one. Moreover, a better choice of the parameters of MFS is numerically observed. It is noteworthy that a problem whose solution has singular points can successfully be solved. It is concluded that the numerical method proposed in this paper is effective for a problem with an irregular domain, singular points, and the Cauchy data with high noises.

Takemi Shigeta, D. L. Young

The method of fundamental solution (MFS) is an efficient meshless method for solving a boundary value problem in an exterior unbounded domain. The numerical solution obtained by the MFS is accurate, while the corresponding matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper basis functions is proposed by the introduction of the modified Trefftz method (MTM). The concrete expressions of the corresponding condition numbers and the solvability by these methods are mathematically proven. Thereby, the optimal parameter minimizing the condition number is also mathematically given. Numerical experiments show that the condition numbers of the matrices corresponding to the MTM and the MMFS are reduced and that the numerical solution by the MMFS is more accurate than the one by the conventional method.

Takemi Shigeta

The purpose of this study is to show some mathematical aspects of the adjoint method that is a numerical method for the Cauchy problem, an inverse boundary value problem. The adjoint method is an iterative method based on the variational formulation, and the steepest descent method minimizes an objective functional derived from our original problem. The conventional adjoint method is time-consuming in numerical computations because of the Armijo criterion, which is used to numerically determine the step size of the steepest descent method. It is important to find explicit conditions for the convergence and the optimal step size. Some theoretical results about the convergence for the numerical method are obtained. Through numerical experiments, it is concluded that our theories are effective.

Takemi Shigeta

An efficient numerical method is proposed for computing the Dirichlet-to-Neumann (DtN) map associated with the exterior Dirichlet problem for the two-dimensional Helmholtz equation with an inhomogeneous term. The exterior solution is approximated by the method of fundamental solutions (MFS). When the source and collocation points are equally spaced on concentric circles, the coefficient matrices arising in the discretization become circulant, which enables efficient evaluation of the discrete DtN map by the fast Fourier transform (FFT). By incorporating the boundary condition defined by the DtN map into a finite element formulation, the original exterior problem posed on an unbounded domain is reduced to an equivalent problem on a bounded computational domain, which can then be solved by the finite element method (FEM). Numerical examples show that the proposed MFS-based approach computes the DtN map with high accuracy, and that the FEM incorporating the boundary condition defined by the DtN map yields accurate solutions for exterior Helmholtz problems with an inhomogeneous term posed on unbounded domains.