几类偏微分方程弱有限元方法研究

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:gordonbest
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
弱有限元方法(weak Galerkin finite element methods,简称WG方法)是最近发展起来的求解偏微分方程的有效数值方法.它的主要思想是利用弱微分算子代替传统意义下的微分算子,然后把其应用到通常的变分形式中以数值求解偏微分方程.弱有限元方法的逼近函数为分片间断多项式,逼近函数在单元与单元之间的联系则通过单元边界上的特定多项式实现.自从弱有限元方法在2011年被王军平和叶秀[1]提出以来,它已经被成功应用到求解各种偏微分方程.本文主要提出了几类偏微分方程的新的弱有限元方法.在第三章,通过引入一个反对称矩阵,构造了反应对流扩散方程的一个新的WG格式.该格式对应的线性系统是正定的,从而很容易获得系统的适定性.该格式的WG元取为形式(Pk(T),Pk-1(e)),即在单元的内部,取k次多项式,在单元的边界,取k-1次多项式.这样,使得该新的WG有限元格式产生更少的自由度.通过误差分析,分别给出了在离散H1范数及L2范数下的误差估计,并作了相应的数值实验研究.理论分析和数值实验表明,在以上范数意义下,该新的WG有限元格式均达到最优收敛阶,即对k次元,在离散H1范数下,其收敛比率为O(hk),在L2范数下,其收敛比率为O(hk+1).在第四章,结合间断有限元思想,通过在单元边界上引入均值:{.}和跃值:[·],提出了反应扩散方程的一个新的MWG(modified weak Galerkin)有限元方法.一方面,该方法产生的自由度比传统的WG方法少;另一方面,与现有的一些MWG方法比较,该方法的有限元函数空间和检验函数空间为同一分片间断多项式空间,即均为Vh={v:v|T ∈ Pk(T),T ∈ T,k≥1},这样的好处是使得误差eh=uh-Q0u(Q0为局部L2投影)属于检验函数空间Vh,从而便于对MWG格式进行误差分析.理论分析和数值实验表明,在离散H1范数及L2范数下,该新的MWG格式均达到最优收敛阶,也就是说,对k次元,在离散H1范数意义下,其收敛比率为O(hk),在L2范数意义下,其收敛比率为O(hk+1).在第五章,利用第四章提出的MWG方法,构造了反应对流扩散方程的MWG有限元格式.在该格式中,通过引入一个反对称矩阵,处理了问题中的对流项,从而保证了格式的正定性.通过误差分析,给出了在离散H1范数和L2范数下的误差估计,并作出了相应的数值实验.理论分析和数值实验表明,在以上范数意义下,该MWG格式均达到最优收敛阶,也就是说,对k次元,在离散H1范数意义下,其收敛比率为O(hk),在L2范数意义下,其收敛比率为O(hk+1).在第六章,借助单元边界上的均值:{·},我们构造了弱函数的新形式:v={u,{v},{▽v}},给出了弱Laplace的新定义,并由此提出了重调和方程的MWG有限元格式.该MWG格式的有限元函数空间和检验函数空间为同一分片间断多项式空间,即均为Sh={v:v|T∈Pk(T),T∈τh,k≥2}.通过误差分析,我们给出了在离散H2范数及L2范数意义下的误差估计.理论分析表明,在离散H2范数意义下,该MWG格式达到最优收敛阶,即对k次元,格式收敛比率为O(hk);在L2范数意义下,当MWG元取2次元时,格式达到次最优收敛阶,即收敛比率为O(h2),当MWG元取k(k≥3)次元时,格式达到最优收敛阶,即其收敛比率为O(hk+1).
其他文献
基于Galerkin框架与正交多项式的优势,传统的谱方法被广泛应用于求解具有高正则度解的微分方程。但在诸多数学与科学计算问题中,由于计算区域的角点奇性、边界条件的不匹配与内部系数的不光滑等这类原因,方程的真解往往具有一定的奇性,从而限制了谱方法的实际应用。为了恢复谱方法对奇性问题的高效性,本文我们设计了几类特殊的谱方法来处理一些当前热门的奇性问题。我们将其细分为以下几个部分:一、广义Jacobi函
本论文分为两部分,第二章,第三章和第四章为第一部分,主要围绕几类流体动力学方程组,即带阻尼的可压Euler方程组,可压磁流体动力学(MHD)方程组和稳态的可压液晶流动力学方程组,研究这三个偏微分方程组解的适定性问题。第五章为第二部分,主要围绕种群动态方程组,研究其逼近方程组强解的适定性问题以及平衡态的不稳定性态问题。在第二章中,我们考虑了三维带阻尼的可压Euler方程组的时间周期解。运用正则化逼近
物联网技术和低功耗电子产品的快速发展,促使无线传感器网络广泛应用于环境监测、结构健康监测、交通安全监测、植入式医疗、智慧家庭和智能电网等领域。在这些应用中,最为突出的问题是为无线传感器网络节点供电。在许多情况下,体积大、寿命短、污染环境的传统电源已经无法满足无线传感器网络对小体积、长寿命、绿色供能技术的需求。由此,作为一项面向无线传感器网络应用的独立供能技术—能量收集技术应运而生。这项技术可以收集
三维激光扫描点云可以用于客观描述和数字化再现真实世界的三维场景。分析点云典型场景(如道路环境)中各类目标之间的位置关系,提高遮挡情况下目标检测的准确性,从而明确导致目标遮挡的空间因素,具有重要的学术和应用意义。然而,实现自动化点云场景可见性分析在目前仍面临诸多挑战。首先,现实场景中受遮挡的目标往往同时受到周围相邻目标的干扰,使得这些处于高噪声点云中的目标检测难度加大。其次,观测模型构建以及遮挡计算
核磁共振(Nuclear Magnetic Resonance,NMR)波谱作为一种无损伤的非侵入性检测技术,业已成为物理、化学、生命科学以及医学等领域研究中一种必不可少的检测手段。高分辨二维谱及多维谱方法可以有效缓解常规一维谱中谱峰拥挤等问题,并获得更多的相关信息,已经成为有机小分子检测以及大分子结构研究中重要的分析工具。然而,高分辨谱图的获取需要有高均匀度的磁场,其空间变化率一般要低于1 pp
C60具有优异的性能,如高温超导电性,因此它的发现激发了对碳团簇广泛的理论和实验研究。C60超导的特性来源于电子-声子耦合,而电声子耦合会随着团簇尺度的减小而增大,因此寻找最小的富勒烯引起了极大的关注。在众多碳团簇中,C20是最小、最简单并且是曲率最大的富勒烯结构。由于C20电声子耦合比C60强,因此可以预测C20。是高温超导的更好的潜在材料。因此,本文研究了基于C20的团簇组装结构。另外,自从石
本文我们考虑几类常见的流体方程,研究它们的强解及相关极限问题,也就是,局部解的粘性消失极限和整体解的衰减这两类问题。更确切地说,粘性消失极限问题是指,当粘性系数或扩散系数趋于零时,粘性流体方程的解收敛到无粘性或理想流体方程的解。在有界区域,边界条件将是一个关键,我们主要考虑的是Slip边界下粘性消失极限问题。而整体解的衰减问题是一个大时间行为,是指当时间趋于无穷大时,能量趋于零,本文也包括衰减率和
E.T.A.霍夫曼是德国18世纪末、19世纪初的著名作家、音乐家、指挥、舞台设计。他被誉为歌德和海涅之间最具世界影响的德国作家。他的很多作品都表现了神秘力量对人的控制,使得人产生分裂和异化,因此他很早就被称为“幽灵霍夫曼”。但他的作品中也确实存在体现和谐、统一的因素,而这些因素都可以归到他最有名的童话小说《金罐》里的“阿特兰提斯”这一概念之中。籍由“阿特兰提斯”表达出来的思想内涵贯穿于霍夫曼的所有
作为计算流体力学研究的一个重要内容,双曲守恒律方程的数值解法在流体力学发展过程中占据着非常重要的地位。在层出不穷的数值计算方法之中,高精度、高分辨率的数值计算方法因为其具有良好的特性,在计算流体力学的发展中占据着重要的地位。本文的主要目的是研究几类具有高分辨率、高精度的数值格式。具体内容如下:首先,我们基于有限体积法思想,通过增加光滑因子中非光滑部分的权重,提出了能有效提高CWENO-Z格式分辨率
随着社会经济的飞速发展和科学技术水平的日益提高,人类面临着越来越复杂的实际决策环境,而决策者认知的模糊性和决策因素的不确定性,导致了决策者往往难以获得确定的决策信息;并且,在实际决策过程中,影响决策者的决策信息越来越多,由此产生的决策数据的维数和量级也越来越大;致使决策者难以进行有效且理性的决策。作为现代决策理论与实践的重要组成部分,直觉模糊多属性决策能够有效的模拟复杂决策环境,而原有基于矩阵理论