Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R^2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.

偏微分方程反问题:模型、算法和应用

偏微分方程反问题是一个重要的数学研究领域,覆盖了偏微分方程、泛函分析、非线性分析、优化算法和数值分析等不同的数学分支,在介质成像、遥感遥测和图像处理等当代重要的工程领域有广泛的应用.基于问题的不适定性,求解这类问题需要引进正则化思想.但是由于模型的复杂性和广泛性,很难建立统一的正则化框架.本文旨在对几类重要的偏微分方程反问题的研究给出一个系统的总结.在阐明偏微分方程反问题起源和特点的基础上,对以电阻抗成像、波场逆散射和介质热成像为应用背景的三类重要的偏微分方程反问题,系统阐述了核心研究问题、已有结果和方法、未来重要的研究方向.最后从反演方法有效实现的角度,对影响偏微分方程反问题数值求解精度和误差估计的主要因素给出了分析.

ON SOURCE ANALYSIS BY WAVE SPLITTING WITH APPLICATIONS IN INVERSE SCATTERING OF MULTIPLE OBSTACLES

We study wave splitting procedures for acoustic or electromagnetic scattering problems. The idea of these procedures is to split some scattered field into a sum of fields coming from different spatial regions such that this information can be used either for inversion algo- rithms or for active noise control. Splitting algorithms can be based on general boundary layer potential representation or Green＇s representation formula. We will prove the unique decomposition of scattered wave outside the specified reference domain G and the unique decomposition of far-field pattern with respect to different reference domain G. Further, we employ the splitting technique for field reconstruction for a scatterer with two or more separate components, by combining it with the point source method for wave recovery. Us-ing the decomposition of scattered wave as well as its far-field pattern, the wave splitting procedure proposed in this paper gives an efficient way to the computation of scattered wave near the obstacle, from which the multiple obstacles which cause the far-field pattern can be reconstructed separately. This considerably extends the range of the decomposition methods in the area of inverse scattering. Finally, we will provide numerical examples to demonstrate the feasibility of the splitting method.

A NEW APPROACH AND ERROR ANALYSIS FOR RECONSTRUCTING THE SCATTERED WAVE BY THE POINT SOURCE METHOD

Consider an inverse scattering problem by an obstacle D belong to R^2 with impedance boundary. We investigate the reconstruction of the scattered field u^s from its far field pattern u^∞ using the point source method. First, by applying the boundary integral equation method, we provide a new approach to the point-source method of Potthast by classical potential theory. This extends the range of the point source method from plane waves to scattering of arbitrary waves. Second, by analyzing the behavior of the Hankel function, we obtain an improved strategy for the choice of the regularizing parameter from which an improved convergence rate （compared to the result of [15]） is achieved for the reconstruc- tion of the scattered wave. Third, numerical implementations are given to test the validity and stability of the inversion method for the impedance obstacle.

NUMERICAL SOLUTION OF THE SCATTERING PROBLEM FOR ACOUSTIC WAVES BY A TWO-SIDED CRACK IN 2-DIMENSIONAL SPACE

The wave scattering problem by a crack F in R2 with impedance type boundary is considered. This problem models the diffraction of waves by thin two-sided cylindrical screens. A numerical method for solving the problem is developed. The solution of the problem is represented in the form of the combined angular potential and single-layer potential. The linear integral equations satisfied by the density functions are derived for general parameterized arcs. The weakly singular integrals and the Cauchy singular integral arising in these equations are computed using a highly accurate scheme with a truncation error analysis. The advantage of the scheme proposed in this paper is, in one hand, the fact that we do not need the analyticity property of the crack and we allow different complex valued surface impedances in both sides of the crack. In the other hand, we avoid the hyper-singular integrals. Numerical implementations showing the validity of the scheme are presented.

On fluorescence imaging:The diffusion equation model and recovery of the absorption coefficient of fluorophores

To quantify fluorescence imaging of biological tissues,we need to solve an inverse problem for the coupled radiative transfer equations which describe the excitation and emission fields in biological tissues.We begin by giving a concise mathematical argument to derive coupled diffusion equations with the Robin boundary condition as an approximation of the radiative transfer system.Then by using this coupled system of equations as a model for the fluorescence imaging,we have a nonlinear inverse problem to identify the absorption coefficient in this system.The associated linearized inverse problem is to ignore the absorbing effect on the excitation field.We firstly establish the estimates of errors on the excitation field and the solution to the inverse problem,which ensures the reasonability of the model approximation quantitatively.Some numerical verification is presented to show the validity of such a linearizing process quantitatively.Then,based on the analytic expressions of excitation and emission fields,the identifiability of the absorption coefficient from the linearized inverse problem is rigorously analyzed for the absorption coefficient in the special form,revealing the physical difficulty of the3-dimensional imaging model by the back scattering diffusive system.

SOLVING THE BACKWARD HEAT CONDUCTION PROBLEM BY DATA FITTING WITH MULTIPLE REGULARIZING PARAMETERS

We propose a new reconstruction scheme for the backward heat conduction problem. By using the eigenfunction expansions, this ill-posed problem is solved by an optimization problem, which is essentially a regularizing scheme for the noisy input data with both the number of truncation terms and the approximation accuracy for the final data as multiple regularizing parameters. The convergence rate analysis depending on the strategy of choosing regularizing parameters as well as the computational accuracy of eigenfunctions is given. Numerical implementations are presented to show the validity of this new scheme.

On the Error Estimate of the Harmonic B_z Algorithm in MREIT from Noisy Magnetic Flux Field

Abstract Magnetic resonance electrical impedance tomography （MREIT, for short） is a new medical imaging technique developed recently to visualize the cross-section conduc- tivity of biologic tissues. A new MREIT image reconstruction method called harmonic Bz algorithm was proposed in 2002 with the measurement of Bz that is a single component of an induced magnetic flux density subject to an injection current. The key idea is to solve a nonlinear integral equation by some iteration process. This paper deals with the convergence analysis as well as the error estimate for noisy input data Bz, which is the practical situation for MREIT. By analyzing the iteration process containing the Laplacian operation on the input magnetic field rigorously, the authors give the error estimate for the iterative solution in terms of the noisy level 6 and the regularizing scheme for determining ABz approximately from the noisy input data. The regularizing scheme for computing the Laplacian from noisy input data is proposed with error analysis. Our results provide both the theoretical basis and the implementable scheme for evaluating the reconstruction accuracy using harmonic Bz algorithm with practical measurement data containing noise.