Beltrami场的数值计算方法

来源 :中国科学院研究生院 中国科学院大学 | 被引量 : 0次 | 上传用户:cjc013
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Beltrami场是三维空间中的一种向量场,该向量场与其旋度的向量积处处为零。Beltrami场的数值计算在太阳物理学、流体力学、核聚变学等领域有着十分广泛的应用。本文旨在研究计算Beltrami场的数值方法。由于Beltrami场所满足的方程是一个带约束的非线性偏微分方程,并且其解不唯一,所以设计一种求解该方程的有效数值方法是一个极具挑战性的研究课题。虽然目前天体物理学界已提出了少量计算Beltrami场的方法,但是它们都缺乏数值和理论上的可靠性,而在计算数学界几乎还没有这方面的研究工作。本文的工作包含两个部分:对有界区域和无界区域两种情况,分别提出了计算Beltrami场的高效迭代算法,并得到了相应的收敛性结果。   本文第一部分研究有界区域上Beltrami场的计算。通过引入一个辅助变量,我们把原来的非线性问题转化成一个含两个变量的最优化问题。为了求解该最优化问题,本文提出了一种新的迭代算法。这种迭代算法可以使磁场变量和所引入的辅助变量在不同的迭代步独立地求解,并且比已有方法具有明显的优势。进一步,本文严格证明了该迭代算法的收敛性,并得到了有限元解的最优误差估计。数值结果验证了理论结果的正确性,并且说明了所提出的算法是高效的、可靠的。   本文第二部分研究无界区域上Beltrami场的计算。根据太阳物理学与电磁学的基本知识,本文导出了一个全日面日冕层磁场的近似模型。该模型在有界区域上就是原来的Beltrami场方程,而在无界区域上退化成一个Laplace方程。利用有限元和边界元方法,我们得到了该近似模型对应的离散耦合系统。为了求解该离散耦合系统,本文提出了一种解耦迭代算法,并严格证明了该迭代算法的收敛性。数值结果表明所提出的算法是稳定的和高效的。
其他文献
非阿基米德动力系统是目前国际上广受关注的新方向,涉及(复)动力系统、数论和代数等多个方向.全纯非阿基米德动力系统始于1981年菲尔茨奖得主Yoccoz和Herman等人的研究,目前
流动性是证券市场的生命力所在,关于流动性的各种学术研究也是金融市场微观结构理论中的热点问题。本文主要考察资产流动性、盈余质量与股票流动性的关系,首先,公司的资产流
基于统计学习理论的机器学习的方法——支持向量机(SVM),有着很好的数学理论基础和泛化性,已经广泛应用到模式识别、数据挖掘等研究领域。它通过引入核函数的方法,把线性不可
郑煤集团老君堂煤矿煤田地处滑动构造带,其首采21121采煤工作面采用单体液压支柱配2.4m长Π型钢梁支护,采取放顶煤一次采全高采煤方法炮采落煤。工作面俯采角为5°~13°,伪倾
寄生物与宿主之间协同进化的研究早已在生命科学领域引起广泛关注.现有研究寄生物与其宿主协同进化的模型几乎都是基于寄生物将会导致宿主种群减少的前提建立的.然而,寄生物在
新形势下如何使预备役部队现役干部尽快成才?衡水某预备役团党委在实践中对这一问题做出了较好的回答。他们以高度的责任感和无私的爱心关心干部的成长,想方设法创造条件,千
多尺度现象在科学和工程问题中普遍存在,典型的例子包括材料力学性质的研究和多孔介质流传播性质的研究。拟连续体方法就是一类典型的多尺度方法。本论文分为两个部分。第一
本文主要研究有限维纯跳Lévy过程的分部积分公式,并利用分部积分公式得到了纯跳Lévy过程驱动的非线性随机微分方程对应转移半群的强Feller性和Harnack不等式。   第一部
矩阵的Fan积和Hadamard积是矩阵理论及其应用的重要问题之一,对它们的特征值的界的估计是目前研究的热点之一.本文对两个非奇异M-矩阵A和B的Fan积的最小特征值q(A★B)的下界,
对于系数为非Lipschitz连续的随机微分方程,一般的Euler方法和Milstein方法可能会出现发散的情形.Higham,Mao,Szpruch[Discrete and ContinuousDynamicalSystems B,2013,pp2083-21