In this paper, we consider an inverse source problem for elliptic partial differential equations with Dirichlet and Neumann boundary data. The unknown source term is to be determined from additional boundary conditions. Unlike the existing methods found in the literature, which usually use some of the boundary conditions to form a boundary value problem for the elliptic partial differential equation and the remaining boundary conditions in the objective functional for optimization to determine the source term, the novel method that we propose here has coupled complex boundary conditions. We use a complex elliptic partial differential equation with a Robin boundary condition coupling the Dirichlet and Neumann boundary data, and optimize with respect to the imaginary part of the solution in the domain to determine the source term. Then, on the basis of the complex boundary value problem, Tikhonov regularization is used to obtain a stable approximate source function and the finite element method is used for discretization. Theoretical analysis is given for both the continuous model and the discrete model. Several numerical examples are provided to show the usefulness of the proposed coupled complex boundary method.
This study considers the inverse source problem for elliptic partial differential equations with both Dirichlet and Neumann boundary data. The unknown source term is to be determined by additional boundary conditions. Unlike the existing methods found in the literature, which usually employ the firstorder in time gradient-like system (such as the steepest descent methods) for numerically solving the regularized optimization problem with a fixed regularization parameter, we propose a novel method with a second-order in time dissipative gradient-like system and a dynamical selected regularization parameter. A damped symplectic scheme is proposed for the numerical solution. Theoretical analysis is given for both the continuous model and the numerical algorithm. Several numerical examples are provided to show the robustness of the proposed algorithm.
In this paper we propose a new class of iterative regularization methods for solving ill-posed linear operator equations. The prototype of these iterative regularization methods is in the form of second order evolution equation with a linear vanishing damping term, which can be viewed not only as a extension of the asymptotical regularization, but also as a continuous analog of the Nesterov's acceleration scheme. New iterative regularization methods are derived from this continuous model in combination with damped symplectic numerical schemes. The regularization property as well as convergence rates and acceleration effects under the Hölder-type source conditions of both continuous and discretized methods are proven.The second part of this paper is concerned with the application of the newly developed accelerated iterative regularization methods with a posteriori stopping rule to the diffusionbased bioluminescence tomography, which is modeled as an inverse source problem in elliptic partial differential equations with both Dirichlet and Neumann boundary data. Several numerical examples, as well as a comparison with the state-of-the-art methods, are given to show the accuracy and the acceleration effect of the new methods. ‡ Corresponding Author
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.