摘要
众所周知,椭圆型偏微分方程Cauchy问题在Hadamard意义下严重不适定,表现在Cauchy数据的微小扰动可导致Cauchy问题解的巨大误差.来源于科学和工程中的许多理论和应用问题可归结为椭圆型方程Cauchy问题,如工程无损探测,地球物理勘查,心脏病学等.Cauchy问题的不适定性给上述问题的研究带来了很大困难,表现在难以构造稳定,高效的算法.一般来说,椭圆型方程Cauchy问题不具有稳定性,但若对该问题的解作先验有界的假设,则可获得稳定性估计,如Holder稳定性,对数稳定性等.
本文研究了三类椭圆型方程柯西问题:Laplace方程柯西问题,Helmholtz方程柯西问题和变系数椭圆方程柯西问题.分析了这些问题的不适定性本质和不适定性程度,给出了几种正则化方法.
本文分为四部分.第一章简要介绍了反问题的概念,反问题的数学特征以及正则化方法.第二章用线方法和谱方法求解了Laplace方程柯西问题,得到了稳定的误差估计.第三章用拟逆方法和拟边界值方法求解了Helmholtz方程和修正的Helmholtz方程,并且得到了Helmholtz方程在一般源条件下的最优误差界.第四章用Fourier方法,修正的Tikhonov正则化方法和小波对偶最小二乘法求解了变系数椭圆方程柯西问题,得到了稳定的误差估计.
此外,我们讨论了所有这些方法的数值实现,给出大量的数值例子来测试所提出的正则化方法各方面的性质.这些测试显示我们的方法是数值可行和有效的.
It is well know that the Cauchy problem of elliptic partial equations are severely ill-posed in the sense of Hadamard:a small noise in the Cauchy data may cause dramatically large error in the solution. The Cauchy problem of elliptic equation arises from many science and engineering problems such as nondestructive testing techniques, geophysics and cardiology. The ill-posedness of Cauchy problem lead to enormous difficulty to research these problems, that is, it is difficult to construct stable, efficient algorithm. In general, the Cauchy problem of elliptic equation has not stability. Under an additional a-priori bound condition, a stable estimate can be obtained, such as Holder stability or logarithmic stability.
In this thesis, we study three kinds of Cauchy problem of elliptic equa-tion:Cauchy problem of Laplace equation; Cauchy problem of Helmholtz equation and Cauchy problem of elliptic equation with variable coefficients.
We analyze the ill-posedness of these inverse problems, and discuss their degree of ill-posedness. For stably computing these problems, we proposed several regularization methods.
This thesis is divided into four parts. The first chapter is preface, the definition of inverse problem, the mathematical characteristic of inverse problem and the regularization method are reported. In the second chap-ter, we use line method and spectral method solve the Cauchy problem of Laplace equation, the stable error estimate is obtained. In the third chap-ter, we apply the quasi-reversibility method and the quasi-boundary value method to Helmholtz equation. Moreover, we get the optimal error bound for the Cauchy problem of Helmholtz equation at general source condition. The fourth chapter is devoted to Fourier method, revised Tikhonov method and Wavelet dual least squares method for the Cauchy problem of elliptic equation with variable coefficients, we obtain stable error estimate.
In addition, we discuss the numerical implementation of all these meth-ods, and give large number of numerical examples to test various properties of the proposed regularization methods. These tests show that our methods are numerically feasible and effective.
引文
[1]Shaw C B Jr, Improvement of the Resolution of an Instrument by Numer-ical Solution of an Integral Equation. J.Math.Anal.Appl,1972,37:83-112.
[2]R. Gorenflo and S. Vessella, Abel Integral Equations:Analysis and Appli-cations, Lectures Notes in Math.1991, Berlin, Springer.
[3]G. Alessandrini, Stable determination of a crack from boundary measure-ments, Proc. R. Soc. Edinburgh A,1993,123,497-516.
[4]U. Tautenhahn, Optimal Stable Solution of Cauchy Problems for Elliptic Equation, Journal for Analysis and its Application,1996,4,961-984.
[5]T. Reginska and K. Reginski, Approximate solution of a Cauchy problem for the Helmholtz equation, Inverse Problems,2006.22,975-989.
[6]L. Elden, F. Berntsson and T. Reginska, Wavelets and Fourier methods for solving the sideways heat equation, SIAM J. Sci. Comput.,2000,21, 2187-2205.
[7]Karen A. Ames, Gordon W. Clark. James F. Epperson, and Seth F. Op-penheimer, A comparison of regularizations for an ill-posed problem, Math-ematics of Computation,1998,224(67),1451-1471.
[8]H. W. Engl, M. Hanke and A. Neubauer, Regularization of Inverse Prob-lems,1996, Kluwer Academic Publisher, Dordrecht Boston London.
[9]Z. Qian, C.L.Fu and X.T. Xiong, Fourth-order modified method for the Cauchy problem for the Laplace equation, Journal of Computational and Applied Mathematics,2006,19,205-218.
[10]X.T. Xiong, Central difference regularization method for the Cauchy prob-lem of Laplace's equation, Applied Mathematics and Computation.2006, 181,675-684.
[11]Ailin Qian, Yujiang Wu, A computational method for a Cauchy problem of Laplace equation, Applied Mathmatics and Computation.2009,207, 478-485.
[12]Ailin Qian, Yujiang Wu, on a quasi-reversibility regularization method for a Cauchy problem of the Helmholtz equation, Journal of Computational and Applied Mathematics.2010,233,1969-1979.
[13]Ailin Qian, A new wavelet method for solving an ill-posed problem, Applied Mathmatics and Computation.2008,203,635-640.
[14]U. Tautenhahn, Optimality for linear ill-posed problems under general source conditions, Numer. Funct. Anal. Optim.,1998,19,377-398.
[15]D. N. Hao and D. Lesnic, The Cauchy for Laplace's equation via the con-jugate gradient method, IMA Journal of Appl. Math.2000,65,199-217.
[16]Ang DD, Nghia NH and Tam NC, Regularized solutions of a Cauchy prob-lem'for the Laplace equation in an irregular layer:a three-dimensional model. Acta Math. Vietnam,1998,23,65-74.
[17]Ang DD, Thanh DN and Thanh VV, Regularized solutions of a Cauchy problem for the Laplace equation in an irregular strip, J. Integral Equations Appl,1993,5,429-441.
[18]Barenblatt GI, Zheltov IP and Kochina IN, Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks (strata), Priklad. Mat. Mekh.,1963,24,852-864.
[19]C. L. Fu et al, Fourier regularization method for a cauchy problem of Laplace equation, Inverse probl. Sci. Eng.2007,16,159-169.
[20]C. Y. Qiu and C. L. Fu, Wavelets and regularization of the Cauchy problem for the Laplace equation, J. Math. Anal. Appl.2008,338,1440-1447.
[21]Baumeister J and Reinhardt HJ, On the approximate solution of a two-dimensional inverse heat conduction problem, Inverse and Ill-posed Prob-lems, ed H G Engl and CWGroetsch (Orlando, FL:Academic),1987,325-344.
[22]Beck JV, Nonlinear estimation applied to the nonlinear inverse heat con-duction problem, Int. J. Heat Mass Transfer,1970.13,703-716.
[23]F. Berntsson, Numerical methods for inverse heat conduction problems, Doctoral thesis, Linkopings University,2001, Linkoping, Sweden.
[24]F. Berntsson and L. Elden, Numerical solution of a Cauchy problem for the Laplace equation, Inverse Problems,2001,17,839-853.
[25]L. Bourgeois, A mixed formulation of quasi-reversibility to solve the Cauchy problem for Laplace's equation, Inverse Problems,2001,21,1087-1094.
[26]Cannon JR, The one-dimensional heat equation, Reading,1984, MA, Addison-Wesley.
[27]Chang JR, Yeih WC and Shieh MH, On the modified tikhonov's regular-ization method for the cauchy problem of the laplace equation, Journal of Marine Science and Technology.2001,9,113-121.
[28]C. L. Fu, H. F. li, Z. Qian and X. T. Xiong,2007, Fourier regularization method for solving a Cauchy problem for the Laplace equation, Inverse Problems in Science and Engineering, in press.
[29]J. Cheng, Y. C. Hon and Y. B. Wang, A numerical method for the dis-continuous solutions of Abel integral equations, In inverse problems and spectral theory,2004,348,233-243.
[30]J. Cheng, Y. C. Hon, T. Wei and M. Yamamoto, Numerical computation of a Cauchy problem for Laplace's equation, ZAMMZ. Angew. Math. Mech, 2001,81,665-674.
[31]J. Cheng, X. Z. Jia and Y. B. Wang, Numerical differnetiation and its applications, Inv. Prob. Sci. Eng.,2007,15,339-357.
[32]F. Berntsson, A spectral method for solving the sideways heat equation, Inverse Probl.1999,15,891-906.
[33]B. D. Coleman and W. Noll, An approximation theorem for functionals, with application in continuum mechanics, Arch. Rational Mech. Anal., 1960,6,355-370.
[34]P. Colli-Franzone and E. Magenes, On the inverse potential problemof elec-trocardiology, Calcolo,1979,16,459-538.
[35]J. Cullum, Numerical differentiation and regularization, SIAM J. Numer. Anal.,1971,8,254-265.
[36]L. Elden, Approximations for a Cauchy problem for the heat equation, Inverse Problems,1987,3,263-273.
[37]L. Elden, Hyperbolic approximations for a Cauchy problem for the heat equation, Inverse Problems,1988,4,59-71.
[38]L. Elden, Numerical solution of the sideways heat equation by difference approximation in time, Inverse Problems,1995,11(4),913-923.
[39]L. Elden, Solving the sideways heat equation by a "method of lines", J. Heat Trans. ASME.,1997,119,406-412.
[40]Falk RS and Monk PB,, Logarithmic convexity for discrete harmonic func-tions and the approximation of the Cauchy problem for Poisson's equation, Math. Comput.,1986,47,135-149.
[41]C. L. Fu, Simplied Tikhonov and Fourier regularization methods on a gen-eral sideways parabolic equation, J. Comput. Appl. Math.,2004,167,449-463.
[42]C. L. Fu and C. Y. Qiu, Wavelet and error estimation of surface heat flux, J. Comput. and Appl. Math.,2003,150,143-155.
[43]C. L. Fu, Y. B. Zhu and C. Y. Qiu, Wavelet regularization for an inverse heat conduction problem, J. Math. Anal. Appl.,2003,288,212-222.
[44]C. L. Fu, X. T. Xiong, H. F. Li and Y. B. Zhu, Wavelet and spectral regular-ization methods for a sideways parabolic equation, Appl. Math. Comput., 2005,160,881-908.
[45]C. L. Fu, X. T. Xiong and Z. Qian, Fourier regularization for a backward heat equation, J. Math. Anal. Appl.,2007,331,472-480.
[46]P. Fu, C. L. Fu, X. T. Xiong and H. F. Li, Two regularization methods and the order optimal error estimates for a sideways parabolic equation, Comput. Math. Appl.2005,49,777-788.
[47]T. Reginska, Application of wavelet shrinkage to solving the sideways heat equation, BIT,2001,41(5),1101-1110.
[48]X.T. Xiong, C.L. Fu, Determining surface temperature and heat flux by a wavelet dual least squares method, J. Coumput. Appl. Math.,2007,201, 198-207.
[49]X.T. Xiong, C.L. Fu, Error estimate on a backward heat equation by a wavelet dual squares method, International Journal of Wavelets, Multires-olution and Information Processing 2007,5(3),389-397.
[50]I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Regional Conference Series in Applied Mathematics. SIAM, Philadelphia, PA,1992.
[51]Golubev G and Khasminskii R,1999, A statistical approach to the Cauchy problem for the Laplace equation, State of the art in probability and statis-tics (Leiden,1999),419-33, IMS Lecture Notes Monogr. Ser.,36, Inst. Math. Statist., Beachwood, OH.
[52]R. Gorenflo, Funktionentheoretische Bestimmung des Aussenfeldes zu einer zweidimensionalen magnetohydrostatischen Konfiguration, Z. Angew. Math. Phys.,1965,16,279-290.
[53]H.H.Qin, T.Wei, Modified regularization method for the Cauchy problem of the Helmholtz equation, Applied Mathematical Modelling,2009,33, 2334-2348.
[54]H.H.Qin, D.W.Wen, Tikhonov type regularization method for the Cauchy problem of the modified Helmholtz equation, Applied Mathematics and Computation,2008,203,617-628.
[55]R. Gorenflo and S. Vessella, Abel integral equations,1991, Springer-Verlag, Berlin.
[56]Groetsch CW, The theory of Tikhono regularization for Fredholm equations of the first kind, Pitman Publishing INC,1984, Toronto.
[57]Groetsch CW, Differentiation of approximately specified functions, Amer. Math. Monthly,1991,98,847-850.
[58]Guo L and Murio DA, A mollified space-marching finite-difference algo-rithm for the two-dimensional inverse heat conduction problem with slab symmetry, Inverse Problems,1991,7,247-259.
[59]J. Hadamard, Lectures on the Cauchy problem in linear partial differential equations, Yale University Press,1923, New Haven.
[60]Han H, Ingham DB and Yuan Y, The boundary element method for the solution of the backward heat conduction equation, J. Comput. Phys.,1995, 116,292-299.
[61]M. Hanke and 0. Scherzer, Error analysis of an equation error method for the identification of the diffusion coefficient in a quasi-linear parabolic differential equation, SIAM J. Appl. Math.,1999,59,1012-1027.
[62]D. N. Hao, A mollification method for ill-posed problems, Numer. Math., 1994,68,469-506.
[63]D. N. Hao and P. M. Hien, Stability results for the Cauchy problem for the Laplace equation in a strip, Inverse Problems,2003,19,833-844.
[64]D. N. Hao, P. M. Hien and H. Sahli, Stability results for a Cauchy problem for an elliptic equation, Inverse Problems,2007,23,421-461.
[65]D. N. Hao and D. Lesnic, The Cauchy problem for Laplace's equation via the conjugate gradient method, IMA J. Appl. Math.,2000,65,199-217.
[66]D. N. Hao, Reinhardt HJ and Schneider A, Numerical solution to a sideways parabolic equation, Int. J. Numer. Meth. Engng.,2001,50,1253-1267.
[67]D.N. Hao, T. D. Van and R. Gorenflo, Towards the Cauchy problem for the Laplace equation, Banach Center Publ.,1992,27,111-128.
[68]M. Imber, Temperature extrapolation mechanism for two-dimensional heat flow, AIAA J.,1974,12,1089-1093
[69]V. Isakov, Inverse problems for partial differential equations,2nd edn,2005, Berlin, Springer.
[70]Johnson CR, Computational and numerical methods for bioelectric field problems, Crit. Rev. Biomed. Eng.,1997,25,1-81.
[71]Keller JB, Inverse problems, Am. Math. Mon.,1976,83,107-118.
[72]Kirsch A, An introduction to the mathematical theory of inverse problems, 1999, Berlin, Springer.
[73]Klibanov MV and Santosa F, A computational quasi-reversibility method for Cauchy problems for Laplace's equation, SIAM J. Appl. Math.,1991, 51,1653-1675.
[74]Knosowski Y, Lieres E and Schneider A, Regularization of a non-characteristic Cauchy problem for a parabolic equation in multiple dimen-sions, Inverse Problems,1999,15,731-743.
[75]Lattes R and Lions JL, The method of quasi-reversibility:Applications to partial differential equations, Elsevier,1969, New York.
[76]Lavrent'ev MM, Romanov VG and Shishatski SP,1986, Ill-Posed Prob-lems of Mathematical Physics and Analysis (Translations of Mathematical Monographs vol 64), Providence, RI:American Mathematical Society.
[77]Lesnic D and Wake GC, A mollified method for the solution of the Cauchy problem for the convection-diffusion equation, Inv. Prob. Sci. Eng.,2007, 4,293-302.
[78]H. F. Li, C. L. Fu and X. T. Xiong, Optimal error bound in Hq-scale of regularized approximation solution for a sideways parabolic equation, Appl. Math. Mech.,2005,26,1238-1244.
[79]刘继军,不适定问题的正则化方法及应用,科学出版社,2005,北京.
[80]Liu L and Murio DA, Numerical experiments in 2-D IHCP on bounded domains part Ⅰ:the "interior" cube problem, Comput. Math. Appl.,1996, 31,15-32.
[81]Manselli P and Miller K, Calculation of the surface temperature and heat flux on one side of a wall from measurements on the opposite side, Annali di Matematica Pura ed Applicata,1980,123,161-183.
[82]Mera NS, Elliott L and Ingham DB, An inversion method with decreasing regularization for the backward heat conduction problem, Numerical Heat Transfer, Part B,2002,42,215-230.
[83]Mera NS, Elliott L, Ingham DB and Lesnic D, An iterative boundary el-ement method for solving the one dimensional backward heat conduction problem, Int. J. Heat Mass Transfer,2001,44,1937-1946.
[84]D. E. Beskos, Boundary element methods in dynamic analysis:Part ⅱ(1986-1996), Applied Mechanics Reviews.1997,50(3),149-197.
[85]Chen, J. T., and Wong, F. C, Dual formulation of multiple reciprocity method for the acoustic mode of a cavity with a thin partition. Journal of sound and vibration.1998,217(1),75-95.
[86]M. R. Bai, Application of bem(boundary element method)-based acoustic holography to radiation analysis of sound sources with arbitrarily shaped geometries. J. Acoust. Soc. Am.1992,92(2),533-549.
[87]H. Cheng, J. Huang, and T. Leiterman, An adaptive fast solver for the mod-ified Helmholtz equation in two dimensions.J.Comput.Phys.2006,211(2), 616-637.
[88]M. Itagaki, Higer order three-dimensional fundamental solutions to the helmholtz and the modified helmholtz equations.Engineering Analysis with Boundary Elements.1995,15(3),289-293.
[89]A. Tikhonov and V. Arsenin, Solutions of ill-posed problems. V. H. Win-stion Sons, Washington, D. c. John Wiley Sons,1977, New York.
[90]L. Marin, L. Elliott, Heggs,P.J.,Ingham,D.B.,Lesnic,D.,and Wen,X.An al-ternating iterative algorithm for the Cauchy problem associated to the Helmholtz equation. Comput.Methods Appl.Mech.Engrg.2003,192(5-6), 709-722.
[91]B. Jin and Y. Zheng, A meshless method for some inverse problems asso-ciated with the Helmholtz equation.Comput.Methods Appl. Mech. Engrg. 2006,195(19-22),2270-2288.
[92]L. Marin and D. Lesnic, The method of fundamental solutions for the Cauchy problem associated with two-dimendional helmholtz-type equa-tions. Computers and Structures.2005,83(4-5),267-278.
[93]B. Jin and Y. Zheng, Boundary knot method for the Cauchy problem as-sociated with the inhomogeneous helmholtz equation. Engineering Analisis with Boundary Elements.2005,29(10),925-935.
[94]Miranker WL, A well posed problem for the backward heat equation, Pro-ceedings of the American Mathematical society.1961,12,243-7.
[95]Murio DA, The mollification method and the numerical solution of ill-posed problems,1993, New York, Wiley.
[96]Murio DA, Mejfa CE and Zhan S, Discrete mollification and automatic numerical differentiation, Computers Math. Appl.,1998,35,1-16.
[97]Z. Qian C. L. Fu and R. Shi, A modified method for a backward heat conduction problem, Appl. Math. Comput.,2007,185,564-573.
[98]Z. Qian, C. L. Fu and X. T. Xiong, Semidiscrete central difference method in time for determining surface temperatures, Int. J. Math. Math. Sci., 2005,3,393-400.
[99]Z. Qian, C. L. Fu and X. T. Xiong, A modified method for a non-standard inverse heat conduction problem, Appl. Math. Comput.,2006,180,453-468.
[100]Z. Qian, C. L. Fu, X. T. Xiong and T. Wei, Fourier truncation method for high order numerical derivatives, Appl. Math. Comput.2006,181,940-948.
[101]C. Y. Qiu and C. L. Fu, Wavelet regularization for an ill-posed problem of parabolic equation, Acta Math. Sci.,2002,22,361-372.
[102]C. Y. Qiu, C. L. Fu and Y. B. Zhu, Wavelets and regularization of the sideways heat equation, Comput. Math. Appl.,2003,46,821-829.
[103]W. L. Miranker, A well posed problem for the backward heat equation, Proceedings of the American Mathematical Society.1961,12(2),243-274.
[104]Qu R, A new approach to numerical differentiation and integration, Math. Comput. Modelling,1996,24,55-68.
[105]Ramm AG and Smirnova AB, On stable numerical differentiation, Math. Comput.,2001,70,1131-1153.
[106]J. J. Reinhardt, H. Han and D. N. Hao, Stability and regularization of discreted approximation to the Cauchy problem for the Laplace's equation, SIAM, Numer. Anal.,1999,36,890-905.
[107]T. Reginska, Sideways heat equation and wavelets, J. Comput. Appl. Math.,1995,63,209-214.
[108]T. Reginska and L. Elden, Soving the sideways heat equation by a Wavelet-Galerkin method, Inverse Problems,1997,13,1093-1096.
[109]T. I. Seidman and L. Elden, An'optimal filtering' method for the sideways heat equation, Inverse Problems,1990,6,681-696.
[110]Tikhonov A N and Arsenin V Y, Solutions of ill-posed problems, Winston and Sons,1977, Washington.
[111]Trong DD, Long NT and Dinh PN, Nonhomogeneous heat equation:Iden-tification and regularization for the inhomogeneous term, J. Math. Anal. Appl.,2005,312,93-104.
[112]Trong DD, Quana PH and Dinh PN, Determination of a two-dimensional heat source, Uniqueness, regularization and error estimate, J. Comput. Appl. Math.2006,191,50-67.
[113]Tu JS and Beck JV, Solution of inverse heat conduction problems using a combined function specification and regularization method and a direct finite element solver, Proc.2nd Ann. Inverse Problems in Engineering Semi-nar, ed J V Beck (East Lansing, MI:College of Engineering, Michigan State University),1989,1-25.
[114]Vani C and Avudainayagam A, Regularized solution of the Cauchy problem for the Laplace equation using Meyer Wavelets, Math. Comput. Model., 2002,36,1151-1159.
[115]Wang YB, Jia XZ and Cheng J, A numerical differentiation method and its application to reconstruction of discontinuity, Inverse Problems,2002, 18,1461-1476.
[116]Weber CF, Analysis and solution of the ill-posed inverse conduction prob-lem, Int. J. Heat Mass Transfer,1981,24,1783-1792.
[117]T. Wei, Numerical methods for inverse problems, Doctoral thesis, City University of Hong Kong.2005.
[118]T. Wei and M. Li, High order numerical derivatives for one-dimensional scattered noisy data, Appl. Math. Comput.,2006,175,1744-1759.
[119]X. T. Xiong, C. L. Fu and H. F. Li, Central difference schemes in time and error estimate on a non-standard inverse heat conduction problem, Appl. Math. Comput.,2004,157,77-91.
[120]X. T. Xiong, C. L. Fu and H. F. Li, Central difference method of a non-standard inverse heat conduction problem for determining surface heat flux from interior observations, Appl. Math. Comput.,2006,173,1265-1287.
[121]X. T. Xiong, C. L. Fu and H. F. Li, Fourier regularization method of a sideways heat equation for determining surface heat flux, J. Math. Anal. Appl.,2006,317,331-348.
[122]X. T. Xiong, C. L. Fu, Z. Qian, and X. Gao, Error estimates of a difference approximation method for a backward heat conduction problem, Int. J. Math. Math. Sci.,2006,1-9.
[123]X. T. Xiong, C. L. Fu and Z. Qian, Two numerical methods for solving a backward heat conduction problem, Appl. Math. Comput.,2006,179, 370-377.
[124]Yildiz B, Ozdemir M, Stability of the solution of backward heat equation on a weak compactum, Appl. Math. Comput.,2000,111,1-6.
[125]T. Hohage, Regularization of exponentially ill-posed problems, Numer. Funct. Anal. Optim.,2000,21,439-464.
[126]Yildiz B, Yetis H, Sever A, A stability estimate on the regularized solution of the backward heat equation, Appl. Math. Comput.,2003,135,561-567.
[127]Yoshimura T and Ikuta K, Inverse heat conduction problem by finite ele-ment formulation, Int. J. Syst. Sci.,1985,16,1365-1376.
[128]Zabaras N and Liu JC, An analysis of two-dimensional linear inverse heat transfer problems using an integral method, Numer. Heat Transfer,1998, 13,527-533.
[129]Songming Hou, Xu-Dong Liu, A Numerical method for solving variable co-efficient elliptic equation with interfaces. Journal of computational physics. 2005,202,411-445.
[130]M.Oeverman, R. Klein, A cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, Journal of computational physics.2006,219,749-769.
[131]W.T.Ang,J.Kusuma,A boundary element method for a second order elliptic partial differential equation with variable coefficients, Engineering Analysis with Boundary Elements.1996,18,311-316.