二元样条函数办法求偏微分方程数值解.pdfVIP

摘 要 样条函数是一种很好的数值逼近工具，它是逼近理论和计算机结合的优秀 产物．在计算机高速发展的推动下，样条函数以其便于计算机存储、计算稳定、 局部支集以及便于交互控制等优点，成为各类工程计算．计算机辅助制造／设计 和几何建模等大型软件的重要数学工具之一。值得一提的是：除了在几何建模 方面样条函数有着无可替代的地位之外，样条有限元在有限元大家族中也占有 重要的地位，它不仅能用较少的自由度得到比其他类型有限元更高的精度之外， 还能与几何建模软件之间方便地共享数据，因为它们都是采用相同的格式来存 储数据的．另外样条函数结合配置法求解流体闯题也受到比较多的关注f641． 本论文主要讨论的是多元样条函数方法及其在偏微分方程数值解中的应 用，事实上也可以认为是一种特殊形式的样条有限元．Brezzi和Fortin宅E[18]中 给出了有限元方法的一个广义角度的定义：”有限元方法是一种通用的技巧， 它的含义是在Hilbert空间中构造一个有限维子空间，从而在这个子空间中对交 分问题实施Ritz．Galerkin方法”．从计算的角度来讲，这个定义并不是非常的 直观和确切(有限元方法的特点是使用了基于网格的具有局部支集基函数的有 限维子空间)。但是它从一个非常抽象的角度概括了类型广泛的有限元逼近方 法．在这些方法中，我们现在主要关心的是标准有限元方法和本文所要介绍的 样条函数方法．我们通常所说的有限元方法指的是标准的有限元计算方法，即 通过在离散单元上构造一个”有限元”参考对象的途径来构造有限元空间的办 法fl，45，72，801，这是狭义的有限元概念，相应地我们称前者为广义的有限元概 念．本文所要介绍的多元样条函数方法也属于广义有限元概念的范畴，是由学 者来明俊等首先提出来的[42，57，581．它与有限元方法的主要区别在于采用了 多元样条分析方法来构造所需要的有限维子空间，接着通过求解带线性约束的 方程组来完成样条函数解的构造，在此过程中有限维子空间并没有被显式地构 造出来，这样的做法在这个子空间不是那么容易构造的时候是方便的，我们将 在本文中看到这些现象． 和样条有限元一样，样条函数方法可以沿用有限元数值分析的许多结论， 前提条件是所使用的样条函数空间必须是存在并且可表示的．这一前提在多元 样条分析的工作中现在已经可以找到很好的理论保证(参考『631以及其中引用 摘 要 的文献)．因此，这篇论文中对于样条函数方法的讨论均是从实际计算的角度 出发的．值得一提的是，如果样条函数方法所采用的样条函数空间和有限元方 法所采用的有限元空间是同一个函数空间，那么得到的近似解也是相同的．因 此样条函数方法似乎没有很明显的优势．然而通过非协调元在很多问题中的成 功应用我们可以看到一些线索：非协调元(包括间断有限元)总是通过改造经 典变分原理的手段来满足近似问题与原变分问题的相容性，它的优势是放松了 单元边界之间的协调约束(连续或光滑)条件，使得有限元空间尽可能好地逼 近真实解空闻，通过这种方式锝到的近似解显然铯更充分地表现真实解的性质． 这就是很多时候非协调有限元比协调有限元应用更广泛的主要原因．样条函数 方法可以看成是一种协调问题的非协调处理办法，因为它把单元边界上的协调 条件单独提出来了．这完全是一个代数的问题，因此也和有限元方法一样具有 通用性，但从形式上看又象是一种非协调有限元． 论文共分五章．第一章首先介绍了一些本文所涉及到的一些有限元基本理 论以及有限元方法处理问题的基本步骤，由于样条函数方法也是包括在广义有 限元范畴中的，因此变分原理以及有限元方法成为了我们讨论样条函数方法的 起点，其中有非常丰富的理论结果和计算技巧是可以从有限元方法借鉴的．接 着，我们简要地介绍了二十多年来关于二元样条分析的一些主要理论结果．主 要从二元样条函数空间存在性、Sobolev空间中的样条函数逼近程度和多元样条 函数空间表示的实际技巧等方面做了必要的介绍，这些都为后面的章节做准备 工作．需要指出的是关于二元样条分析的文献是非常丰富的，我们这里所提到 的一些结论完全是为了方便讨论样条函数方法而从中提取出来的． 第二章介绍了三角形上的二元B形式及其一些常用的局部计算方法，比 如de Casteljau求值算法，升次、求导、求积以及求内积算法等．这些算法是样 条函数方法中计算单元方程常用的．其中有不少可以加速椭圆边值问题以及高 次方程的离散过程．进一步我们利用B．形式的光滑性条件，将形式与样条函 数之间建立了转换关系，得到了多元