定常線性薛定諤方程的一種高精度數值解法

張紅梅，尹江華

（1.湖南工業大學 理學院，湖南 株洲 412007；2.株洲歐科億數控精密刀具股份有限公司，湖南 株洲 412008）

1 研究背景

薛定諤方程是量子力學的基本方程之一，其地位與牛頓方程在經典力學中的地位相當。薛定諤方程已被廣泛應用于理論物理、核物理、等離子物理、電磁波理論、化學、光學工程、地震學等領域之中。但對于這些實際應用中的大量復雜問題，用經典的解析方法和實驗研究方法大多難以求解。隨著計算機技術的快速發展，數值解法已成為其另一重要且有效的求解方法。薛定諤方程是耦合方程組，它還涉及復函數，采用經典的離散方法在單層網格下離散原方程所得的離散系統，通常規模較大，耗費大量的計算時間。因此，構建快速的數值解法十分必要。

有限體積元算法除了處理邊界靈活、計算簡便外，還能夠在一定程度上保持物理量的局部和全局守恒性。這些性質使得有限體積元算法成為一種重要方法。近年來，學者們分別采用不同的有限體積元算法求解橢圓方程[1]、Navier-Stokes 方程[2]、擴散方程[3]等。

在求非對稱問題或非正定問題的有限元解時，提出了兩網格離散算法[4]，并將該算法應用到求解非線性邊值問題[5]、Navier-Stokes 問題[6]、流體問題[7]等。隨后，兩網格有限元算法又被成功地應用于求解耦合方程組[8]，并進行了一系列的推廣與應用[9-12]。在上述研究的基礎上，C.S.Chien 等[13]采用兩網格差分法求解非線性薛定諤方程組；Wu L.[5]用兩網格混合有限元算法求解非線性薛定諤方程組。

近年來，國內外許多學者將兩網格離散思想與有限體積元算法相結合，用于求解橢圓方程定解問題[14]、非線性雙曲方程定解問題[15]、拋物方程定解問題[16]等。本文將有限體積元算法計算的靈活性與兩網格算法的快速性相結合，構建兩網格有限體積元算法，用于求解耦合方程組——定常線性薛定諤方程邊值問題：

λΔ 為動能算子，且λ為正實數；

u(x)、V(x)、f(x)分別為未知函數、位勢函數、右端項，均為復函數。

在實際計算時要將解的實部與虛部分開成如下等價的耦合方程組：

式中：u1(x)、u2(x)分別為u(x)的實部和虛部函數；

V1(x)、V2(x)分別為V(x)的實部和虛部函數；

f1(x)、f2(x)分別為f(x)的實部和虛部函數。

為了方便，假設

2 線性有限體積元算法及誤差估計

方程（2）的變分形式可描述為：

設Th={ei}、Zh分別是形狀正則的三角形網格剖分所有單元和剖分點的集合[9]，h為剖分步長。定義線性有限元空間

由式（4）和（8）可得正交關系

并在條件（3）成立的情況下，問題（8）的有限元解uh滿足

由文獻[1]可知算子Ih滿足引理1。

為簡化方程組（2）的計算，將方程組（2）的前兩式在控制體bp上積分，并利用 Green 公式可得

式中n為控制體邊界的單位外法向量。

對上兩式等號左邊第二項進行近似計算，即控制體bp上的積分，被積分函數ui(x)取近似值ui(bp),i=1,2，進行計算，得，將上面兩個方程兩邊分別乘以Ihw1和Ihw2，并對bp求和，則有

為了對式（13）和（14）化簡，先給出如下引理2[2]。

利用引理2 可得

由直接驗證可知

將式（15） 和（16） 代入式（13）和（14），可得如下線性有限體積元變分問題：

定義范數容易證得

由引理3 和Lax-Milgram 定理知，變分問題（17）的解存在且唯一。

根據文獻[9]有引理4。

引理5 證畢。

根據文獻[2]中的式（3.10）有

利用上式和式（23）可得式（24）。引理6 證畢。

引理7設u、uh分別為方程（2）和（17）的解，則有

證由式（21）、（17）和（19）得

引理7 證畢。

因而有

證由式（8）和（17）有

當式（22）和（24）中p=q=2 時，再聯合式（25）則有

在式（26）中，令v=uh-uh，再由式（9）可得式（27）。引理8 證畢。

證由式（7）、（10）和（27）得

定理1 證畢。

3 兩網格有限體積元算法及誤差估計

為了簡化式（18），將其拆成兩部分并引入記號：

兩網格有限體積元算法：

因而有

證將式（17）減去式（29），再由式（18）得

由假設條件（3）及式（12）、（19）、（25）和（28）可得

即得式（31）。

再由式（28）和（31）得

定理2 證畢。

4 兩網格有限體積元算法數值實驗

表1 λ=1 時兩種算法的結果比較Table 1 Comparison of results between the two algorithms with λ=1

表2 λ=100 時兩種算法的結果比較Table 2 Comparison of results between the two algorithms with λ=100

5 結語

薛定諤方程的實際應用比較廣泛，方程的離散系統規模龐大。本文在已求出定常線性薛定諤方程兩網格有限元解的基礎上，重新構造了薛定諤方程的一種新的快速數值解法——兩網格有限體積元算法。從理論上證明了該算法達到最優收斂階，并用數值實驗驗證了該算法能節省大量的計算時間，極大地提高了求解效率。