论文部分内容阅读
偏微分方程被广泛用于描述现代科学和工程计算中的许多实际问题。对于解析解不存在的微分方程,有效地进行数值求解就尤为重要。随着数学理论的深入研究和计算方法的不断发展,有限差分法,有限体积法以及有限元方法已经能够较好得求解大部分微分方程。然而,对于数值解比较奇异的方程,如果采用均匀网格需要大量的计算资源,尤其是高维问题可能会超出计算机的计算能力。移动网格方法根据数值解的特点对网格进行重新分布,可以在不浪费计算资源的前提下有效减少计算误差。同时,在实际数值计算中,选取均匀的时间步长可能需要较长的计算时间。时间自适应方法可以在计算的过程中不断调整时间步长,从而提高数值计算的效率。 1925年,Einstein预测了在极低温度下气体中的粒子会处于相同的量子态。1995年,在稀薄的碱金属气体中发现了Bose-Einstein凝聚态(BEC)。该问题引起了物理学家和数学家的广泛关注,通常用非线性薛定谔(NLS)方程来描述Bose-Einstein凝聚态的单粒子性。大量的科研工作者在理论和数值方面对非线性薛定谔方程进行了研究,并提出了一系列的数值求解方法。在无穷势阱下,当粒子间存在强相互作用时,Bose-Einstein凝聚的基态解中会出现边界层。因此,利用均匀网格计算该基态解需要大量的计算资源。同时,求解Bose-Einstein凝聚的基态解就是在限制条件下求能量泛函的极小值点,该能量在数值计算初期变化剧烈,而在接近收敛时变化非常缓慢,因此采用均匀时间步长需要较长的计算时间。根据该问题数值解在空间和时间上的特点,在空间上利用移动网格方法,在时间上利用时间自适应方法能够有效提高数值计算的效率。 本文主要介绍一种时空自适应有限元方法来求解Bose-Einstein凝聚态的基态解。首先,本文介绍了自适应方法和非线性薛定谔方程的相关理论知识。其次,介绍了一维问题基于等分布原理的移动网格方法,二维问题基于调和映射的移动网格方法以及时间自适应方法。然后,本文分析了不同势阱下Bose-Einstein凝聚态基态解的数值特点,提出了如何在空间上实现移动网格技术以及在时间上实现自适应。基于时空自适应有限元方法,本文给出了一维和二维情况下Bose-Einstein凝聚基态解的数值算例,分析比较了均匀网格和移动网格的数值结果,并指出时空自适应方法的有效性。