复杂区域上粘性%2f非粘性耦合方程的一个谱元方法论文.pdfVIP

前言 谱方法是一种既古老又新兴的求解偏微分方程的数值方法，它起源 的未知数提供非常精确的逼近，主要缺点是不能处理复杂的物理区域． 这个困难通常是通过区域分解加以克服．由Patel’altlJ等人发展起来的 谱元方法即是基于这一技巧，其基本思想是将复杂区域分解成若干个较 小的四边型或六面体区域，然后在每个小区域上使用谱方法． 在一些数学物理及工程实际问题中，常常要求在不同的子区域中选 用不同的数学模型．如气体绕飞行体流动，在边界层附近为粘性流，在 边界层外为无粘流，两者有不同的数学模型，为粘性／非粘性分析．粘 性／非粘性分析是区域分解技术中较为特殊的一类技巧，显然它不同于 一般分解技术之处在于前者在不同分解区域中考虑不同类型的方程，因 此具有更大的灵活性，但也具有更大的难度．粘性／非粘性分析在计算 流体中的研究已有较长的历史，但始终未能得到较大的发展，其根本原 因在于以前的分解计算方法未能显而易见地达到减少计算量的目的．文 【ll弓I进了求解不可压粘性／非粘性耦合方程的一种薪的耦合技术，通过 构造整体变分形式，耦台问题转化成一个整体鞍点问题，类似于由纯粘 性方程产生的鞍点问题．基于这种技巧，改进后的Uzawa方法得以成功 地应用于耦合问题计算，从而有效克服了传统分解技术中子区域迭代造 成的重复计算问题，达到了减少计算量的目的． 本文尝试将上述技巧推广到更一般的计算区域．利用基于谱元法的 四边形剖分，粘性／非粘性方程得以在很大一类复杂计算区域中计算实 现．象在单区域上的谱方法情形。复杂区域上粘性／非粘性耦合方程的 谱元法主要困难在于构造兼容的速度和压力空间以及压力二次项，使得 著名的Babu；ka-Brezzi’8 Inf-Sup条件得以满足，从而保证谱元离散解的 存在性．通过子区域转换和Sobolev空间上的投影算子，并借助在有限 元方法中采用的Boland—Nicolaides技巧【…，我们证明了粘性区域和非 粘性区域上的lnf_Sup条件分别成立．最后利用文[2】2中处理粘性／非 粘性交面条件的变分方法，证明了耦合区域上整体lnbSup条件成立， 从而构造了粘性／非粘性耦合方程的一个适定的谱元离散，证明了数值 解的存在唯一性． 本文轮廓如下。第一部分介绍一些基本的概念和基本定理．第二部 分我们首先描述了基本的粘性／非粘性耦合方程和整体变分问题，提出 了整体变分问题的谱逼近；然后讨论并分析了此离散问题的稳定性，主 要论证了[nLSup条件成立，从而证明了解的存在唯一性 在这里，我对导师许传炬教授三年来对我学业上的精心指导和鼓励， 表示衷心的感谢．导师严谨的治学态度，开拓进取的精神时常激励着我， 导师的每一悉心指导，耳提面命，总使收益匪浅．我真诚地感谢他．同 时我也由衷地感谢师母悉心指导和支持，以及生活上的关心．另外，我 还得感谢其他的老师和同学的帮助，它使我得以顺利完成三年的学业． 2 复杂区域上粘性／非粘性耦合方程 的一个谱元方法 黄凤辉 摘要 本文给出了分析复杂区域上粘性／非粘性耦合问题的一个谱元法． 通过基于谱元法的区域分解技巧，构造并分析了复杂区域上粘性／非粘 性耦合问题的一个高阶算法．通过整体变分方法并借助推广了的L。。一 一种处理粘性非粘性交面条件的变分法，证明了离散解的存在唯一性． 关铡l词：粘性／非粘性耦合问题，谱元法：解的存在唯一性 Inf-Sup条件．。1 3 I．一些基本概念和基本定理 §1．1鞍点问题及其遣近理论 设x，M为两个Hilbert空间，对应范数分别为”恢，”llM．定义 (1)X×X上的连续的双线型形式n(·，·)，(2)X×M上的连续的双线型 形式6{·，·)． 考虑问题t对给定的f∈x’，找u∈x．p∈M使得 』n(叩)+6(唧)=(，，”)，铷∈x，f11)