两类发展方的部分迎风有限元方法.pdfVIP

两类发展方程的部分迎凰查腿显杰法 马智勇 (山东大学数学与系统科学学院，济南250ioo) 中文摘要 r 』 ＼有限元法是求解值煎坌左盈的一类重要方法，是R．Courant于40年代首先提 出的：随着计算机的迅速发展，有限元方法也得到了巨大的发展，广泛地应用到各 个领域．有限元方法具有网格剖分灵活，适用区域广泛，易于处理第二、三类边值 问题，计算矩阵稀疏，对有限维空间函数光滑性要求低，精度高等诸多优点． 近年来，有限元方法被广泛用于流体力学计算．作为流体力学的一个重要模 型，对流扩散方程已被许多数值分析工作者广泛研究．在众多的方法中，出现了一 法，流线扩散法，特征一Petrov-Galerkin有限元法等著名方法． 另一类重要的方法是迎风有限元法．这类方法易于实现，适用于多维问题，尤 其是它们保留了原问题的两个重要性质：极值原理和质量守恒原理．显然，无论从 物理或是力学方面，还是从数值分析方面，这两个性质都是十分重要的．P 全迎风格式为了保持数值解的稳定性，常常带有过多的人工粘性，从而抹平了 解的变化，降低了数值解的精度．为了使数值解既稳定又有较好的精度，Hughes 等人113]利甩改变数值积分中结点位置的方法xCZ．维等参双线性元提出了部分迎 风有限元格式，这是一种非协调有限元过程，八十年代初，Ikeda综合了他本人早 分迎风格式(16】，获得了比全迎风格式好的结果． 本文采用都分迎风有限元方法，对两类发展方程作出了数值分析，全文共分 两章． 第一章，关于非线性半导体器件睡盔闺壁的郡分迎风有限元格式．半导体器 件的数疽模拟在半导体工业中有重要的实用价值，因而其数值方法的研究受到了 普遍的重视．半导体器件的数值模拟需要研究扩夔攘型的非线性偏微分方程组的 初边值问题． f考虑如下的数学模型： l假设nC群是有界的多边形区域，I，=(0，司．亍是一定常数．未知函数是 电场位势中，电子、空穴浓度e，P和温度函数T —Aq,=o加一e+Ⅳ(z)) 塞=V．{D。扛)Ve-#A王)ev曲】一R(e，pT) 宝=V．{DA坍p饥【z)pV妒卜球，p，T) (T、t)∈n×J 初始条件： o n E e(o，0)=eo(z)，p(卫，0)=；加(z)，T(x，0)=To(z)， 边界条件： r×J 筹j，=嘉I，=赛I，=舅I，=o，(州)e 著名学者Gummel于1964年提出用序列迭代法计算这类问题，开创了半导体 器件数值模拟这一新领域。袁益让教授于1996年用特征差分方法对这类阐题进行 了分析．得到了最优的工2估计。本文第一章对此类非线性问题提出了一种易于实 现的部分迎风有限元格式．在一定的正则性假设下，得到了误差估计． lJ“一“^肌’+0砷一妒^I|Ⅳ≤cp+^s+^铲1)} 第二章，考虑线性Sobolev方程的Dirichht问题。悚u=u(。。￡)，使得 f J (z，￡)∈Q×(o，研 uI锄=0 (z．t)∈an l x(o．T1 l u(茁．o)=u。(卫) 王∈n，t=o Sobolev方程在流体穿过裂缝岩石的渗透理论，土壤中湿气迁移问题，不同介 质间的热传导问题等许多数学物理方面有着广泛的应用．关于Sobolev问题解的存 在唯一性讨论，参见文献[911101．在许多实际应用中，常常涉及到带有对流项的 Sobotev方程，研究对流占优的Sohoiev方程的有限元数值解法有着重要的实际意 义． 本文第二章对此类问题提出了如下的部分迎风有限元格式： { =(，n，讥) v％∈vo“ uo=l^uo(z)