反幂势薛定谔方程与分数阶Laplace方程的谱方法

来源 :中国工程物理研究院 | 被引量 : 0次 | 上传用户:FIGOWEN7
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
基于Galerkin框架和正交多项式的优势,谱与谱元方法被广泛应用于求解具有高正则性解的微分方程。然而在许多科学计算问题中,方程本身及其真解往往具有一定的奇性,从而限制了谱方法的实际应用效果。为了恢复谱方法处理奇性问题的高效性,我们需要根据实际问题设计完全克服方程及其真解奇性的谱方法。本文针对反幂势薛定谔方程和分数阶Laplace方程这两类极具代表性的难点奇性问题,开展具有指数收敛阶的高效谱与谱元方法研究,主要内容包括:(1)提出任意维球体与任意多边形区域上反幂势薛定谔方程特征值问题的高效谱与谱元方法。对含反位势的薛定谔特征值问题-△u+c2∣x∣-2u+γ∣x∣-1u=λu,我们首先在任意维球上构造贴合特征函数奇性的Sobolev正交函数系,并设计新型谱方法的数值格式和算法。在此基础上,我们提出了任意多边形区域上该问题的谱元方法:以奇点为中心划出一个适当大小的圆盘,圆内采用能贴合奇性的非多项式Sobolev基函数谱方法,圆外采用(曲边)四边形或者(曲边)三角形网格与C0-协调谱元,圆盘内外单元之间通过mortar元或其他非协调元方法进行信息交换。针对含反立方势的薛定谔特征值问题-△u+c32∣x∣-3u+c22∣x∣-2u+c12∣x∣-1u=λu,我们设计了类似于含反位势的薛定谔特征值问题的新型高效谱方法。最后,本文将反幂势方程谱方法运用到矩形区域上燃料电池数学模型特征值问题-(?)2/(?)x2u-1/x2(?)2/(?)y2=λu的数值求解中。以上问题的谱方法均取得了指数阶的收敛速度。(2)提出任意维全空间上反幂势薛定谔方程的新型谱方法。以逆平方势薛定谔方程为例,从算法设计、数值实验和理论分析上全面探讨了全空间上反幂势薛定谔方程特征值问题与具有rα型奇异解的薛定谔方程源问题的谱方法。为此,我们定义了关于Muntz级数{αn}的Muntz-Hermite函数Lkαn(r2)e-r2/2αn+1-d/2Yn(?)(ξ),并针对特征值问题和源问题分别选择合适的Muntz级数构造基函数,在此基础上设计Galerkin谱逼近格式和算法。得到的离散代数系统的矩阵具有高度稀疏性,对源问题来说右端项积分还可以通过快速Fourier算法来计算。我们给出了特征值问题中数值特征值和特征函数的最优误差估计以及源问题中Galerkin谱方法的数值解逼近真解的收敛性分析,数值实验结果验证了理论分析的正确性。(3)提出并研究了任意维球上分数阶Laplace方程的高效谱方法。采用球调和函数和Jacobi多项式构造球正交多项式/函数作为分数阶Laplace方程的高效基函数,并在对称变分形式的基础上设计相应的Galerkin谱方法逼近格式与求解算法。在数值分析方面,我们给出了分数阶Laplace方程特征值问题中数值特征值的上下确界估计、特征函数的误差估计和代数特征系统矩阵条件数的紧致估计以及源问题中解的正则性分析。发展了球正交多项式逼近理论,并应用该逼近理论证明了分数阶Laplace方程的收敛性分析。
其他文献
扩散方程在工业制造、油藏模拟、天体物理、等离子体物理等领域具有广泛的应用.因此设计高效精确的数值格式求解这一类方程至关重要.在数值格式的设计中,由于网格变形以及扩散系数的各向异性和间断性等因素,使得建立一般网格上满足保物理特性的数值格式一直是当前具有挑战性的一个重要课题.本论文主要针对这一课题展开,包括五部分的内容:(1)针对三维扩散问题单元中心型保正有限体积格式设计中的节点未知量插值消去方法——
为了利用不同量子系统的优点,科学家们提出了混合量子系统的概念。在各种混合系统中,腔自旋波混合系统受到了人们的广泛关注。一般情况下,腔自旋波混合系统由三维微波腔中的微波光子(腔模)和钇铁石榴石(YIG)样品中的自旋波量子(Kittel模)组成。由于自旋之间的海森堡交换作用,YIG材料内部的自旋系综可以同时拥有较高的自旋密度和较低的耗散率,这使得腔自旋波混合系统可以达到强耦合区,甚至超强耦合区。利用腔
计算流体力学对于航空航天、国防工程、汽车制造、生物科学等各个领域都起着重要的作用。其中对于爆轰及燃烧问题的研究,一直以来都是一个热点问题。常见的爆轰模型包括CJ模型与ZND模型两种。其中CJ模型由Chapman以及Jouguet分别提出,在该模型中,爆轰波阵面被看成一个伴随着化学反应的跳跃间断,分离了已燃介质和未燃介质,且化学反应瞬间完成达到热化学平衡。第二种是由Zeldovich、von Neu
关于孤立环境下的原子过程在过去几十年里已经开展了大量的理论和实验工作,然而在许多领域(如天体物理、激光等离子体、惯性约束聚变)中必须考虑等离子体环境效应的影响,目前理论上可以通过屏蔽势来考察这种环境效应对原子结构和动力学过程的影响。在过去的理论工作中,主要是采用Debye-Huckel势来描述弱耦合等离子体中带电粒子之间的相互作用,采用离子球势(不包含温度)描述强耦合等离子体中粒子之间的相互作用。
量子点是一种人工合成的特殊的材料,在其内部电子自旋会受到强烈的束缚,从而形成离散的轨道能级,与自然界中的原子十分类似。其内部的电子自旋具有良好的相干性,可以作为量子比特或者逻辑门,这在自旋电子学,量子信息和功能材料领域有非常广阔的应用前景。然而限制其实际应用的最大阻碍是电子自旋在环境的干扰下会快速地退相干,从而导致量子信息的丢失。因此寻找那些能在环境的干扰下保持相干性的电子自旋是十分重要的研究课题
激光间接驱动惯性约束聚变利用黑腔产生的X射线辐射源,驱动聚变靶丸内爆压缩氘(D)氚(T)燃料,实现中心点火。由于驱动激光器造价高昂,必须通过高温辐射源的烧蚀压和球形聚心收缩内爆,提高点火热斑的压力,以节省实现点火的驱动能量。但目前激光惯性约束聚变领域面临着内爆辐射驱动不对称和流体力学不稳定性两个科学问题的严重挑战。所以这两个问题的研究对于实现间接驱动惯性约束聚变的点火具有重要的物理意义和应用价值。
激光间接驱动惯性约束聚变利用黑腔产生高温高品质的X射线辐射源,以驱动聚变靶丸内爆压缩,实现中心点火,对可控核聚变等高能量密度物理研究具有基础物理意义。当前黑腔构型在驱动对称性、耦合效率及激光等离子体不稳定性等方面尚有提升空间。围绕高性能黑腔设计和黑腔辐射场表征,本论文开展黑腔构型设计及辐射流角分布研究,主要研究内容如下:1.联合研制三维视角因子程序IRAD3D;该程序集三维激光排布建模和物理仿真分
在经受辐照的材料中会有大量的缺陷形成并保留下来,不同类型的缺陷之间发生错综复杂的交互作用,最终会造成材料微观结构和宏观性能的变化。晶界是一种高效的点缺陷吸收阱,晶界-点缺陷交互作用与材料辐照肿胀的动力学密切相关,而且晶界对于提升材料的耐辐照性能作用显著,但是目前对于晶界-点缺陷交互作用问题的认识还需要进一步深入。本论文将围绕金属型核燃料的辐照效应问题,采用实验和理论模拟的方法,研究辐照环境下材料的
本论文的主要内容包括三部分:(1)守恒型并行差分格式设计与理论分析;(2)保正型并行差分格式设计与理论分析;(3)非完美接触界面问题的迭代方法设计与理论分析.在第一部分中,通过分析具有无条件稳定、二阶数值精度的一维并行差分格式,给出了一个推广形式的并行差分格式,首先,对于一维问题提出了一种加权形式的数值流以及权重的选取范围,然后将此格式推广到二维,最后将格式推广到n-维(n ≥3).理论证明了此守
自从核武器诞生以来,人们就充分意识到了核扩散问题的危害,在用核武器威慑对手的同时,又担心受到核武器的威胁。为了限制核扩散的发展趋势,国际社会做了大量工作,其中国际原子能机构的成立是标志性成就之一。中国作为联合国常任理事国和核大国,应该在核不扩散及核安保工作中发挥自己的积极作用。解决核不扩散和核军控问题,除了国际环境因素之外,在技术层面也提出了较高的要求。由于235U自身裂变能力较弱,因此对于铀部件