摘要
本文主要研究了两个问题,在文中分别加以阐述.下面就两部分内容作简单介绍.
文章的第一部分,探讨了稳态温控问题自适应有限元方法的收敛性.我们知道在实际应用中,自适应有限元方法(参见[2][5][6])常被用来数值求解最优控制问题.对本文而言,我们将致力于研究稳态温控问题.在文章的开端,首先讨论了当模型问题采用标准有限元逼近时,分别在L2-模和H1-模下得到的后验误差估计.然后,基于上述结果,我们得到了一个误差缩减率.更进一步,结合数据振荡的缩减率,可以构造模型问题的一个自适应有限元算法.在后续部分,我们证明了该算法是收敛的.
文章的第二部分,围绕气体在块状纳米材料中的压力扩散行为展开.由于任何固体材料都或多或少地能够渗透过一些气体,气体在固体中扩散的规律与气体自扩散公式有相同的形式.描述扩散行为的基本参量,即扩散系数是由气体一固体组合的性质决定.对于特定的气体,它在不同固体材料中的扩散系数是不同的.因此,对于未知固体材料,研究气体在其内部的扩散系数也能从一方面帮助人们了解它的结构、属性,这对工程领域的研究有着实际的意义.本文通过模拟气体在未知纳米材料块中的压力扩散实验,建立对应的偏微分方程模型,试图反演出气体在固体中的压力扩散系数.借助物理学相关知识,本文将反问题最终转化为最优化问题来处理.由于实验测量数据均为离散值,所以主体部分采用了偏微分方程数值解、数值积分(参见[29][30][31][32])的方法来处理.最优化问题拟用黄金分割搜索法(参见[26][28][33])求解.
In present paper, we focus on two problems. They will be discussed in two parts respectively.
In the first part, we concentrate on the Convergence of Adaptive Fi-nite Element Methods for Static Temperature Control Problem. Adaptive Finite Element Methods (see [2][5][6]) for numerically solving optimal control problems are often used in practice. In this part, we concentrate on the Static Temperature Control Problem. At the beginning, we discuss a posteriori error estimates for the standard finite element approxima-tion of the model problem in L2-norm and H1-norm respectively. Then based on the above results, we obtain an error reduction rate. Further-more, together with a reduction rate of data oscillation, we construct an adaptive FEM algorithm for the model problem. Next, we prove that this algorithm is convergent.
In the second part, we work on the Pressure diffusion behavior of the gas in nanomaterials. Gas can penetrate any solid matter more or less, and the diffusion rule of gas in solid is similar to its self-diffusion. Dif-fusion coefficient is the basic parameter, which is decided by the nature of the gas-solid combination. In different solid, the diffusion coefficient of gas is different as well. As an application, for an unknown kind of solid, to study the diffusion coefficient of gas in it can help us understand its structure. By simulating the pressure diffusion behavior of the gas in nanomaterials, we establish the corresponding partial differential equa-tions, and we try to calculate the diffusion coefficient. With knowledge of physics, the problem above can be transformed into optimization prob-lem. Because the data in this experiment is discrete, we have to solve the problem using numerical methods (see [29][30][31][32]). Optimization algorithm (see [26][28][33]) is needed also.
引文
[1]Wenbin Liu, Ningning Yan. Adaptive Finite Element Methods for Optimal Con-trol Governed by PDEs. Beijing:Science Press,2008.
[2]Pedro Morin, Ricardo H. Nochetto, Kunibert G. Siebert. Convergence of Adap-tive Finite Element Methods. SIAM,2002,Vol 44, No.4:631-658.
[3]Scott L R and Zhang S. Finite element interpolation of non-smooth functions satisfying boundary conditions. Math Comp,1992,54:483-493.
[4]Verfurth R. A Review of a Posteriori Error Estimation anf Adaptive Mesh-Refinement Techniques, New York:Wiley-Teubner,1996.
[5]Peter Binev, Wolfgang Dahmen, Ron DeVore. Adaptive Finite Element Methods with convergence rates. Numer. Math.,2004,97:219-268.
[6]A. SCHMIDTA, K.G. SIEBERTA. LBERT-Software for scientific computations and applications, Acta Math. Univ. Comenian.(N.S.),2001,70:105-122.
[7]F. BREZZI, M. FORTIN. Mixed and Hybrid Finite Element Methods. Springer-Verlag, New York,1991.
[8]A. COHEN, W. DAHMEN, R. DEVORE. Adaptive wavelet methods for el-liptic operator equations-Convergence rates. Math. Comp.2001,70:27-75.
[9]S. DAHLKE, W. DAHMEN, K. URBAN. Adaptive Wavelet Methods for Sad-dle Point Problems-Optimal Convergence Rates, preprint.
[10]S. DAHLKE, R. HOCHMUTHA, K. URBAN. Adaptive wavelet methods for saddle point problems. M2AN Math. Model. Numer. Anal.2000,34:1003-1022.
[11]R.B. KELLOGG. On the Poisson equation with intersecting interfaces. Appl. Anal.1975,4:101-129.
[12]W. MITCHELL. A comparison of adaptive refinement techniques for elliptic problems. ACM Trans. Math. Software.1989,15:210-227.
[13]P. MORIN, R.H. NOCHETTO, K.G. SIEBERT. Data oscillation and conver-gence of adaptive FEM, SIAM J. Numer. Anal.2000,38:466-488.
[14]P. MORIN, R.H. NOCHETTO, K.G. SIEBERT. Local problems on stars:A posteriori error estimators, convergence, and performance. Math. Comp. to appear.
[15]Morin,P., Nochetto, R., Siebert, K. Data Oscillation and Convergence of Adap-tive FEM. SIAM J. Numer. Anal.2000,38:466-488.
[16]Binev, P., Dahmen,W., DeVore, R., Petrushev, P.. Approximation classes for adaptive methods. Serdica Math. J.2002,28:391-416.
[17]Binev, P., DeVore, R.. Fast computation in tree approximation. Numer. Math., this issue.
[18]Agouzal, A.. A posteriori error estimator for nonconforming finite element methods. Appl.Math. Lett.1994,7:1017-1033.
[19]Binev, P., Dahmen,W., DeVore, R..Adaptive finite element methods with con-vergence rates. Numer. Math.2004,97:219-268.
[20]Brezzi, F., Fortin, M..Mixed and hybrid finite element methods. Springer, Berlin-Heidelberg-NewYork,1991.
[21]Oswald, P.:On a hierarchical basis multilevel method with nonconforming P1 elements. Numer. Math.1992,62:189-212.
[22]Kanschat, G., Suttmeier, F.-T..A posteriori error estimates for nonconforming finite element schemes. Calcolo.1999,36:129-141.
[23]Morin, P., Nochetto, R.H., Siebert, K.G.. Data oscillation and convergence of adaptive FEM. SIAM J. Numer. Anal.2000,38(2):466-488.
[24]王烈衡,许学军,有限元方法的数学基础,北京:科学出版社,2007.
[25]龙长江,齐欢,张翼鹏.气体储罐泄漏扩散模型研究.数学的实践与认识,2006,Vol36,No.6:111-114.
[26]贺国强,孟泽红.求解热传导反问题的一种正则化Newton型迭代法.应用数学和力学,2007,Vol 28, No.4:479-486.
[27]宋文彪.空气污染控制工程.第二版.北京:冶金工业出版社,1988.14—16.
[28]马明书,付立志,谷淑敏.解三维抛物型方程的一个ADI格式. Pure and Applied Mathematics,2009,Vol.25, No.1:8-13.
[29]宋巨龙,钱富才,基于黄金分割的全局最优化方法.计算机工程与应用,2005.94-130.
[30]梅树立,张森文,徐加初,郭幸福,非线性动力学方程的自适应精细积分.暨南大学学报,2005,Vo1.26:319-323.
[31]周长银,贺国平,李述山,基于随机搜索的一个约束优化问题的全局收敛算法.高校应用数学学报,2005,20(2):197-206.
[32]邱爱保,重积分数值解的MATLAB实现.宜春学院学报,2008, Vol.130:15-16.
[33]陆金甫,关治.偏微分方程数值解.北京:清华大学出版社,2004.