The method of fundamental solutions for a biharmonic inverse boundary determination problem

In this paper, a nonlinear inverse boundary value problem associated to the biharmonic equation is investigated. This problem consists of determining an unknown boundary portion of a solution domain by using additional data on the remaining known part of the boundary. The method of fundamental solutions (MFS), in combination with the Tikhonov zeroth order regularization technique, are employed. It is shown that the MFS regularization numerical technique produces a stable and accurate numerical solution for an optimal choice of the regularization parameter. A. Zeb on study leave visiting the University of Leeds.     

Use of non-negative constraint in Tikhonov regularization for particle sizing based on forward light scattering

Abstract

The set of linear equations in the inversion of particle size distribution (PSD) based on forward light scattering is an ill-posed problem. In order to solve the inverse problem of this kind, a number of inversion algorithms have been proposed. The regularization algorithm can reconstruct the PSD, but in usual case, the solution may contain negative values and is strongly oscillating. Owing to the natural reason, the solution should be non-negative and smooth. In this paper, a simple non-negative constraint (NNC) is used with a combination of the Tikhonov regularization. Simulations and experiments show that the regularization with NNC can achieve more reasonable results.     

The method of fundamental solutions for a time-dependent two-dimensional Cauchy heat conduction problem

We investigate an application of the method of fundamental solutions (MFS) to the time-dependent two-dimensional Cauchy heat conduction problem, which is an inverse ill-posed problem. Data in the form of the solution and its normal derivative is given on a part of the boundary and no data is prescribed on the remaining part of the boundary of the solution domain. To generate a numerical approximation we generalize the work for the stationary case in Marin (2011) [23] to the time-dependent setting building on the MFS proposed in Johansson and Lesnic (2008) [15], for the one-dimensional heat conduction problem. We incorporate Tikhonov regularization to obtain stable results. The proposed approach is flexible and can be adjusted rather easily to various solution domains and data. An additional advantage is that the initial data does not need to be known a priori, but can be reconstructed as well.     

The dual reciprocity boundary element method for solving Cauchy problems associated to the Poisson equation

This paper presents an application of the dual reciprocity method (DRM) to a class of inverse problems governed by the Poisson equation. Here the term inverse refers to the fact that the boundary conditions are not fully specified, i.e. they are not known for the entire boundary of the solution domain. In order to investigate the ability of the DRM to reconstruct the unknown boundary conditions using overspecified conditions on the accessible part of the boundary we consider some test problems involving circular, annular and square domains. Due to the ill-posed nature of the problem, i.e. the instabilities in the solution of these problems, the DRM is combined with the Tikhonov regularization method.     

A coupled complex boundary method for the Cauchy problem

Considered in this paper is a Cauchy problem governed by an elliptic partial differential equation. In the Cauchy problem, one wants to recover the unknown Neumann and Dirichlet data on a part of the boundary from the measured Neumann and Dirichlet data, usually contaminated with noise, on the remaining part of the boundary. The Cauchy problem is an inverse problem with severe ill-posedness. In this paper, a coupled complex boundary method (CCBM), originally proposed in [Cheng XL, Gong RF, Han W, et al. A novel coupled complex boundary method for solving inverse source problems. Inverse Prob. 2014;30:055002], is applied to solve the Cauchy problem stably. With the CCBM, all the data, including the known and unknown ones on the boundary are used in a complex Robin boundary on the whole boundary. As a result, the Cauchy problem is transferred into a complex Robin boundary problem of finding the unknown data such that the imaginary part of the solution equals zero in the domain. Then the Tikhonov regularization is applied to the resulting new formulation. Some theoretical analysis is performed on the CCBM-based Tikhonov regularization framework. Moreover, through the adjoint technique, a simple solver is proposed to compute the regularized solution. The finite-element method is used for the discretization. Numerical results are given to show the feasibility and effectiveness of the proposed method.     

Numerical boundary identification for Helmholtz-type equations

We study the identification of an unknown portion of the boundary of a two-dimensional domain occupied by a material satisfying Helmholtz-type equations from additional Cauchy data on the remaining known portion of the boundary. This inverse geometric problem is approached using the boundary element method (BEM) in conjunction with the Tikhonov first-order regularization procedure, whilst the choice of the regularization parameter is based on the L-curve criterion. The numerical results obtained show that the proposed method produces a convergent and stable solution     

Numerical solution of a Cauchy problem for Laplace equation in 3-dimensional domains by integral equations

A numerical method based on integral equations is proposed and investigated for the Cauchy problem for the Laplace equation in 3-dimensional smooth bounded doubly connected domains. To numerically reconstruct a harmonic function from knowledge of the function and its normal derivative on the outer of two closed boundary surfaces, the harmonic function is represented as a single-layer potential. Matching this representation against the given data, a system of boundary integral equations is obtained to be solved for two unknown densities. This system is rewritten over the unit sphere under the assumption that each of the two boundary surfaces can be mapped smoothly and one-to-one to the unit sphere. For the discretization of this system, Weinert’s method (PhD, Göttingen, 1990) is employed, which generates a Galerkin type procedure for the numerical solution, and the densities in the system of integral equations are expressed in terms of spherical harmonics. Tikhonov regularization is incorporated, and numerical results are included showing the efficiency of the proposed procedure.     

不适定热传导反问题的两种求解方法及其比较

本文以不适定热传导反问题为对象，采用两种方法进行了求解。一种方法基于对具有测量误差的边界条件进行适当的微扰，使之化为适定问题；另一种方法基于Ｔｉｋｈｏｎｏｖ的正则平滑思想，对反问题中的输入数据进行平滑处理，以便使函数及其一阶导数均实现一致逼近。通过计算与求解表明，两种方法均能得到具有一定精度与稳定性的结果，其中以正则化法更为理想     

基于自适应步长双参数正则化算法的超声波过程层析成像图像重建

提出了一种新的自适应步长双参数正则化算法,对超声波层析成像系统检测浆体浓度分布进行图像重建。该算法利用转换矩阵将超定解作为先验信息,嵌入到正则化泛函中,避免重建图像被过度平滑,不仅成像速度较快且重建图像具有较高分辨率。仿真实验结果表明,相比于Tikhonov正则化算法以及Landweber算法,自适应步长双参数正则化算法重建图像的相关系数有明显提高并且边界信息更加可靠。     

A comparative study of Domain Embedding Methods for regularized solutions of inverse Stefan problems

In this paper, various Domain Embedding Methods (DEMs) for an inverse Stefan problem are presented and compared. These DEMs extend the moving boundary domain to a larger, but simple and fixed domain. The original unknown interface position is then replaced by a new unknown, which can be a boundary temperature or heat flux, or an internal heat source. In this way, the non-linear identification problem is transformed into a linear one in the enlarged domain. Using different physical quantities as the new unknown leads to different DEMs. They are analysed from various points of view (accuracy, efficiency, etc.) through two test problems, by a comparison with a common Front-Tracking Method (FTM). The first test has a smooth temperature field and the second one has some singularities. The advantage of the DEMs in solving the inverse problem and in computing the corresponding direct mapping is shown. In the direct problem, high-order accurate schemes could be obtained more easily with the DEMs than with the FTM. In the inverse problem, an iterative regularization and a Tikhonov regularization have been employed. For the FTM, the iterative regularization is not efficient—the solution oscillates when the data are noisy. As for the Tikhonov regularization, it requests special care to choose an adequate penalty term. In contrast, both the regularizations give good results with all the considered DEMs, except for the second test problem at the beginning (t=0⁺) when the value of the heat flux and the heat source tends to ∞. Slightly different regularization effects have been obtained when using different DEMs. Finally, an automatic choice of the optimal regularization parameter is also discussed, using data with different noise levels. We propose the use of the curve of the residual norm against the regularization parameter. © 1997 John Wiley & Sons, Ltd.     

The Spring-Damping Regularization Method and the Lie-Group Shooting Method for Inverse Cauchy Problems

The inverse Cauchy problems for elliptic equations, such as the Laplace equation, the Poisson equation, the Helmholtz equation and the modified Helmholtz equation, defined in annular domains are investigated. The outer boundary of the annulus is imposed by overspecified boundary data, and we seek unknown data on the inner boundary through the numerical solution by a spring-damping regularization method and its Lie-group shooting method (LGSM). Several numerical examples are examined to show that the LGSM can overcome the ill-posed behavior of inverse Cauchy problem against the disturbance from random noise, and the computational cost is very cheap.     

Simultaneous determination for a space-dependent heat source and the initial data by the MFS

In this paper we propose a numerical algorithm based on the method of fundamental solutions for recovering a space-dependent heat source and the initial data simultaneously in an inverse heat conduction problem. The problem is transformed into a homogeneous backward-type inverse heat conduction problem and a Dirichlet boundary value problem for Poisson's equation. We use an improved method of fundamental solutions to solve the backward-type inverse heat conduction problem and apply the finite element method for solving the well-posed direct problem. The Tikhonov regularization method combined with the generalized cross validation rule for selecting a suitable regularization parameter is applied to obtain a stable regularized solution for the backward-type inverse heat conduction problem. Numerical experiments for four examples in one-dimensional and two-dimensional cases are provided to show the effectiveness of the proposed algorithm.     

Boundary knot method for some inverse problems associated with the Helmholtz equation

The boundary knot method is an inherently meshless, integration‐free, boundary‐type, radial basis function collocation technique for the solution of partial differential equations. In this paper, the method is applied to the solution of some inverse problems for the Helmholtz equation, including the highly ill‐posed Cauchy problem. Since the resulting matrix equation is badly ill‐conditioned, a regularized solution is obtained by employing truncated singular value decomposition, while the regularization parameter for the regularization method is provided by the L‐curve method. Numerical results are presented for both smooth and piecewise smooth geometry. The stability of the method with respect to the noise in the data is investigated by using simulated noisy data. The results show that the method is highly accurate, computationally efficient and stable, and can be a competitive alternative to existing methods for the numerical solution of the problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Modified boundary Tikhonov-type regularization method for the Cauchy problem of a semi-linear elliptic equation

In this paper, we consider a Cauchy problem for the semi-linear elliptic equation and use a modified boundary Tikhonov-type regularization method to overcome its ill-posedness. The existence, uniqueness and stability of the regularization solution are proven. Under an a-priori bound assumption for the exact solution, a convergence estimate of Hölder type for this method is obtained. Finally, an iterative scheme is proposed to calculate the regularization solution, numerical results show that this method works well.     

Inverse wave scattering in 2-D elastic solids

Summary This paper is concerned with an inverse problem for two-dimensional elastic solids. It seeks to recover the subsurface density profile based on the measurements obtained at the boundary. The method considers a temporal interval for which time dependent measurements are provided. It formulates an optimal estimation problem which seeks to minimize the error difference between the given data and the response from the system. It uses a boundary regularization term to stabilize the inversion. The method leads to an iterative algorithm which, at every iteration, requires the solution to a two-point boundary value problem. Several numerical results are presented which indicate that a close estimate of the unknown density function can be obtained based on the boundary measurements only.     

数值求解Symm积分方程的正则化MINRES方法

Symm积分方程在位势理论中具有重要应用,它是Hadamard意义下的不适定问题。离散该方程将产生对称线性不适定系统。基于GCV准则,并应用截断奇异值分解,本文提出数值求解Symm积分方程的正则化MINRES方法。与Tikhonov正则化方法相比,在数据出现噪声的情况下,新方法能有效地求得Symm积分方程的数值解。     

The MFS for numerical boundary identification in two-dimensional harmonic problems

In this study, we briefly review the applications of the method of fundamental solutions to inverse problems over the last decade. Subsequently, we consider the inverse geometric problem of identifying an unknown part of the boundary of a domain in which the Laplace equation is satisfied. Additional Cauchy data are provided on the known part of the boundary. The method of fundamental solutions is employed in conjunction with regularization in order to obtain a stable solution. Numerical results are presented and discussed.     

Boundary knot method for the Cauchy problem associated with the inhomogeneous Helmholtz equation

In this paper, the boundary knot method is extended to the solution of inhomogeneous equations, and it is applied to the Cauchy problem associated with the inhomogeneous Helmholtz equation. Here, we assume that the boundary condition is specified only on a part of the boundary, and the boundary conditions on the remaining part of the boundary are to be determined with the assistance of additional data. Since the resulting matrix equation is highly ill-conditioned, a regularized solution is obtained by employing the truncated singular value decomposition to solve the matrix equation arising from the boundary knot method, with the regularization parameter determined by the L-curve method. Numerical results are presented for several examples with smooth and piecewise smooth boundaries. The numerical verification shows that the proposed numerical scheme is accurate, stable with respect to data noise, and convergent with respect to decreasing the amount of noise in the data.     

Reconstruction of multiplicative space- and time-dependent sources

This paper presents a numerical regularization approach to the simultaneous determination of multiplicative space- and time-dependent source functions in a nonlinear inverse heat conduction problem with homogeneous Neumann boundary conditions together with specified interior and final time temperature measurements. Under these conditions a unique solution is known to exist. However, the inverse problem is still ill-posed since small errors in the input interior temperature data cause large errors in the output heat source solution. For the numerical discretisation, the boundary element method combined with a regularized nonlinear optimization are utilized. Results obtained from several numerical tests are provided in order to illustrate the efficiency of the adopted computational methodology.     

On a certain class of incorrectly stated problems in analytical theory of thermal conductivity and in the theory of indirect measurements

The variational problem of stationary thermal conductivity in an inhomogeneous solid is formulated. It is assumed that the boundary conditions on the boundary of the body are unknown. In order to obtain a unique and stable solution one requires measurement of the temperature at one point and correct selection of the regularization parameter.