非薄板腐蚀探测问题的数值解法.pdf

第31卷第4期 南京师大学报 (自然科学版) V01．3lNo．4 2008年 12月 JOURNALOFNANJINGNORMALUNIVERSITY(NaturalScienceEdition) Dec，2008 非薄板腐蚀探测问题的数值解法 陈 宇，黄建国 (上海交通大学数学系，上海 200240) [摘要] 腐蚀探测问题是一个数学物理方程反问题 ，它通过外边界上可获知的电场数据反演求解腐蚀系数．通 常所涉及的数据是带有噪声误差的．在无需假设板或管的厚度很薄的条件下 ，提 出了一个基于Dirichlet，Neu— mann条件求解腐蚀系数的变分模型．该模型最终 由最优化领域中的拟牛顿迭代法实现数值求解．给出若干理 论分析 ，并用数值实验结果说明求解方法的可行有效性． [关键词] 有限元 ，拟牛顿法，反问题，正则化 [中图分类号]0242 [文献标识码]A [文章编号]1001-4616(2008)04-0029-04 NumericalM ethodforCorrosionDetectionProblem inNon．SheetCase ChenYu，HuangJianguo (DepartmentofMathematics，ShanghaiJiaotongUniversity，Shanghai200240，China) Abstract：Theproblem ofrecoveringthecorrosioncoefficientinaninaccessibleinteriorpartfrom theelectricinformation inanaccessiblepartofaphysicaldomainisstudied，whichisatypicalinverseproblem inmathematicalphysics．Usual— ly，theprescribeddatahavenoiseerror．A variationalfomr ulationisproposedtoderivethecorrosioncoefficient，based ontheDirichlet—Neumanndataontheaccessiblepart．Thequasi—Newtoniterativemethodinoptimizationisusedtosolve thenumericalsolutionofthisvariationalproblem．Sometheoreticalanalysisisprovided，andthenumerical experiment showsthatthemethodiSeffective． Keywords：finiteelement，Quasi—NewtonMethod，inverseproblems，regularization 探测固体材料 (内部)不可接触区域腐蚀程度的问题具有重要理论意义和直接应用价值．特别地，此 类问题是无损伤探测研究的基本问题．本文所考虑的是具有传导性的金属材料区域．若在可接触区域施加 一 定的电流，这时会在区域内部产生一个 电场，且其电位势满足 Laplace方程．本文给出了一种通过可测 量的电位势和电流数据，判断不可见区域腐蚀情况的数值方法． 1 问题的描述 考虑文献 [1]中的腐蚀探测问题模型： au(，Y)=0， (，Y)∈力， M(，)=gc(，Y)， (，Y)∈厂c， a (，Y)= (，Y)， (，Y)∈Fc， “(，Y)=gD(，Y)， (，Y)∈FD， a (，Y)十 (，Y) (，Y)：0， (，Y)∈厂，， 其中， 是R(d=2，3)中的Lipschitz区域，n是指向