王建云,　田智鯤

(1.湖南工業(yè)大學(xué) 理學(xué)院, 湖南 株洲 412007; 2.湖南工程學(xué)院 理學(xué)院, 湖南 湘潭 411104)

超收斂分析是提高有限元方法數(shù)值精度和效率的一種強有力的工具, 國內(nèi)外許多學(xué)者做了大量的研究[1-4]. 薛定諤(Schr?dinger) 方程是量子力學(xué)最基本的方程, 在原子、分子、非線性光學(xué)、等離子物理、電磁波理論、核物理等領(lǐng)域中被廣泛應(yīng)用. 關(guān)于其數(shù)值求解方法有許多的研究[5-9], 而在其超收斂方面的研究不是很多[10-12]. 在本文中, 考慮如下二維定常非線性薛定諤方程：

(1)

式中Ω?R2為矩形區(qū)域, 未知函數(shù)u(x) 和右端項函數(shù)f(x) 都為復(fù)數(shù)值, 勢能函數(shù)V(x)∈L(Ω) 為實數(shù)值且非負, 即存在某實數(shù)V0>0 使得V(x)≥V0.

(2)

式中a(u,v)=(u,v)+(Vu,v).

a(uh,vh)+(|uh|2uh,vh)=(f,vh),vh∈Sh.

(3)

a(Phw,vh)=a(w,vh),?vh∈Sh.

(4)

1　超收斂誤差估計

引理1[6]若函數(shù)w(x)∈H2(Ω), 則其投影Phw(x) 有如下估計

‖w-Phw‖≤Ch2‖w‖2.

(5)

引理2[5]對任意函數(shù)v(x)∈Sh, 有如下逆不等式成立

‖v‖≤Ch-1‖v‖.

(6)

引理3設(shè)uh為 (3) 的解, 則有如下誤差估計

‖uh‖≤C,‖uh‖≤C.

(7)

證明在 (3) 中取vh=uh, 有(uh,uh)+(Vuh,uh)+(|uh|2uh,uh)=(f,uh).由于 (uh,uh)≥0,(|uh|2uh,uh)≥0, 有

V0‖uh‖2≤(Vuh,uh)≤|(f,uh)|≤C‖f‖‖uh‖,

因此‖uh‖≤C‖f‖. 即式(7)的第一式得證. 另外, 由式(6)和式(5) 可得

‖uh-Phuh‖≤Ch-1‖uh-Phuh‖≤Ch‖uh‖2,

注意到‖uh‖≤‖Phuh‖+‖uh-Phuh‖, 得到‖uh‖≤‖Phuh‖+Ch‖uh‖2≤C. 因此, (7) 的第二式得證.

定理1設(shè)u和uh分別為式(2)和式(3)的解, 且u∈H2(Ω), 則有如下誤差估計

‖uh-Phu‖≤Ch2,‖uh-Phu1‖≤Ch2.

(8)

證明由 (2)和(3)可得a(u-uh,vh)+(|u|2u-|uh|2uh,vh)=0,?vh∈Sh.令u-uh=ρ-θ, 其中ρ=u-Phu,θ=uh-Phu, 則有

a(ρ-θ,vh)+(|uh|2(ρ-θ),vh)+((|u|2-|uh|2)u,vh)=0.

由 (4) 有a(ρ,vh)=0, 這樣

a(θ,vh)+(|uh|2θ,vh)=(|uh|2ρ,vh)+((|u|2-|uh|2)u,vh),

取vh=θ有(θ,θ)+(Vθ,θ)+(|uh|2θ,θ)=(|uh|2ρ,θ)+((|u|2-|uh|2)u,θ).由于 (θ,θ)≥0,(|uh|2θ,θ)≥0, 則有

V0‖θ‖2≤(Vθ,θ)≤|(|uh|2ρ,θ)|+|((|u|2-|uh|2)u,θ)|,

因此V0‖θ‖≤‖|uh|2ρ‖+‖(|u|2-|uh|2)u‖.由 (5) 和 (7) 的第二式得 ‖|uh|2ρ‖≤‖uh‖‖ρ‖≤Ch2‖u‖2,另外

‖(|u|2-|uh|2)u‖≤‖u‖(‖u‖+‖uh‖)‖u-uh‖≤

‖u‖(‖u‖+‖uh‖)(Ch2‖u‖2+‖θ‖),

若記γ0=‖u‖(‖u‖+‖uh‖), 則得 (V0-γ0)‖θ‖≤Ch2‖u‖2. 假定γ0

即式(8) 的第一式得證, 類似可證明式(8)的第二式.

定理2設(shè)u和uh分別為式(2)和式(3)的解,uI∈Sh為u的雙線性插值函數(shù), 且u∈H3(Ω), 則有如下誤差估計

‖uh-uI‖1≤Ch2.

(9)

證明[11] 中已經(jīng)證得

(10)

2　數(shù)值實驗

算例求解非線性薛定諤方程式(1), 其中Ω=[0,1]×[0,1], 勢能函數(shù)V=1, 右端函數(shù)f(x) 選取其滿足精確解為：

u=x(1-x)y(1-y)+ix(1-x)y2(1-y).

表1　數(shù)值結(jié)果

[1]陳傳淼, 黃云清. 有限元高精度理論[M]. 長沙: 湖南科學(xué)技術(shù)出版社, 1995.

[2]林群, 嚴寧寧. 高效有限元構(gòu)造與分析[M]. 保定: 河北大學(xué)出版社, 1996.

[3]WAHLBIN L B. Superconvergence in Galerkin finite element methods[M]. Berlin: Springer, 1995.

[4]HUANG Y Q, YANG W, YI N Y. A posteriori error estimate based on the explicit polynomial recovery[J]. Nat Sci J Xiangtan Univ, 2011, 33(3):1-12.

[5]AKRIVIS G D, DOUGALIS V A, KARAKASHIAN O A. On fully discrete Galerkin methods of second-order temporal accuracy for the nonlinear Schr?dinger equation[J]. Numer Math, 1991, 59:31-53.

[6]JIN J C, WU X N.Convergence of a finite element scheme for the two-dimensional time-dependent Schr?dinger equation in a long strip[J]. J Comput Appl Math, 2010, 234(3): 777-793.

[7]KARAKASHIAN O, MAKRIDAKIS C. A space-time finite element method for the nonlinear Schr?dinger equation: the continuous Galerkin method[J]. SIAM J Numer Anal, 1999, 36(6):1779-1807.

[8]LU W Y, HUANG Y Q, LIU H L. Mass preserving discontinuous Galerkin methods for Schr?dinger equations[J]. J Comput Phys, 2015, 282:210-226.

[9]WANG J Y, HUANG Y Q. Fully discrete Galerkin finite element method for the cubic nonlinear Schr?dinger equation[J]. Numer Math Theor Meth Appl, 2017, 10(3):670-687.

[10]LIN Q, LIU X Q. Global superconvergence estimates of finite element method for Schr?dinger equation[J]. J Comput Math, 1998, 16(6):521-526.

[11]TIAN Z K, CHEN Y P, WANG J Y. Superconvergence analysis of bilinear finite element for the nonlinear Schr?dinger equation on the rectangular mesh[J]. Adv Appl Math Mech, 2018, 10(2):468-484.

[12]WANG J Y, HUANG Y Q, TIAN Z K, et al. Superconvergence analysis of finite element method for the time-dependent Schr?dinger equation[J]. Comput Math Appl, 2016, 71:1960-1972.