无界区域一上维线性薛定谔方程紧致差分数值模拟

发布时间：2019-05-14 04:02

【摘要】：本文考虑如下的薛定谔方程初值问题的数值解,其中,h表示普朗克常量,m为粒子的质量,φ(x,t)为波函数,i=(?)为虚数单位。初始函数φ0(x)和源项f(x,t)具有紧致性。区域的无界性给上述问题的数值求解带来很大的困难。目前,人工边界方法是解决此困难的有效方法之一。引入人工边界之后,无界区域就被划分为两部分：有界的计算区域和无界的外部区域。然后,就可以在人工边界上添加合适的人工边界条件,从而,在有限计算区域上的带人工边界条件的边值问题是原问题的很好的近似。因此,本文采用了两种构造人工边界的方法；一是基于帕德逼近和逆拉普拉斯变换,本文第二章第二节中采用了线性薛定谔方程的局部吸收边界条件；二是考虑有正厚度的”人工边界”,一个包围计算区域的人工层,使大部分外行的波都能被吸收,即完美匹配层(PML)方法。第三章第二节中,我们采用复坐标延伸法,使用了线性薛定谔方程的完美匹配层。由此,无界区域上的薛定谔方程就简化为有界的计算区域上的初边值问题。有限差分方法是数值求解初边值问题的有效方法之一。采用传统的一阶和二阶有限差分格式很难得到高精度的数值解,除非使用大量的网格节点。然而,这会大大增加计算量和计算时间。为解决传统有限差分格式的缺点,一种自然的方法就是构造高阶的紧致有限差分格式。这种格式不仅可以在不增加网格节点的前提下提高数值解的精度,而且可以减少计算量。因此,本文第二章第三节、第三章第三节中,分别介绍了相应初边值问题的四阶紧致差分格式。在时间层上,用中心差商代替时间导数,得到时间上的二阶精度。为得到时间上的四阶精度,第二章第四节介绍四阶龙格-库塔方法。最后,我们给出了数值算例,证明了紧致差分格式的有效性。
[Abstract]:......
【学位授予单位】：山东师范大学
【学位级别】：硕士
【学位授予年份】：2016
【分类号】：O241.82

【参考文献】