廣義非線性Schro¨dinger方程的一個(gè)新的守恒差分格式

1、:1000-4424(2006)01-0079-08O1引言非線性S方程在物理學(xué)(如非線性光學(xué)Q等離子物理學(xué))的應(yīng)用中扮演()chrodingerNPS著重要的角色,因此最近二十年來對(duì)該方程進(jìn)行了大量的研究R在數(shù)值求解方面,也已提出了許多方法,其中包括文S總結(jié)和提出了五個(gè)有限差分格式,提出了一種三層七點(diǎn)格1TS2T式,提出了US3,4T-V格式和加權(quán)格式R本文將討論更一般的廣義非線性Schrodinger方程:()WNPSaY_YYb0,c,0eZgh,dececfXZY(c,Z)bY(c,Z)b0,0eZgh,dfY(c,0)bY(c),c,0dececf這里是實(shí)數(shù),是正的實(shí)數(shù)R,_ai用(

2、式乘以i用(式乘以得到實(shí)部部分,則(式1.1)1.1)1.1)j(1.3)Y得到虛部部分,XZ有著如下的關(guān)系:(1.1)(1.2)(1.3)收稿日期:2004-02-20基金項(xiàng)目:國(guó)家自然科學(xué)基金(教育部高等學(xué)校博士點(diǎn)基金(10471079);20030422049)dGQ=高校應(yīng)用數(shù)學(xué)學(xué)報(bào)e輯|u|dx=constant,p|u|+|u|E=+2RR第21卷第1期(1.4)dx=constant.A2差分格式及其守恒量對(duì)平面區(qū)域B做網(wǎng)格剖分,取空間步長(zhǎng)時(shí)間步長(zhǎng)為,EFBG,=(Ix)J,xxHExxKLCDDC為網(wǎng)比記N為原問題其中x=x+M,=N=G,1,O,KPN=1,2,=,uxLLM

3、DMCNM2L(x)精確解u的近似,我們引入一些符號(hào)U(1.1)S(1.T)(,)xLMNQN+1N+1N+1NI1NI1NI12NI1uu+uuu+u,Vu=.M+1I2MMI1RM+1I2MMI1RxM2Q2Q(x)(x)對(duì)(式提出下面的差分格式U1.1)2N+1Vu=xMN+1NI12N+12NI1u+Vu=pVIxMxMRQLN+1N+1I2NI12N+12NI1I2NI1N+1NI1u|uuu+|uu+u|u.M|+|M|M|+O+|M|M|M|RMMRQQ+2定理X差分格式(關(guān)于離散電荷和離散能是守恒的,即2.1)NNI1GQ=Q=O=Q,NNI1GE=E=O=E.其中,NN+1N

4、Q=(YuY2+YuY2)J2,N+1+2N+2E=p(YuY+YuY)J2+x(|u+|u|).x|MZ+2M=GNN+1xNx這里KNN2MNxKYuY=xu|,YuY=xZ|ZM=GN+1NI1和將(式與做內(nèi)積,然后取2.1)x(u+u)J2分別表示取實(shí)部和虛部,虛部部分得到U證其中,N+1NI1N+1NI1N+1NI1_=(uIu(u+u,YuY2IYuY2MM)MM)=Z2ML2L=GKK2N+12NI1N+1NI1=p(Vu+Vu(u+u=G,MM)xMxM)Z2M=GN+1N+1+2N+12NI1+2NI1|uu+O+|uu+|ua=cM|+|M|M|M|M|RQZ2M+2=GK

5、2+2(1.5)R(2.1)(2.2)(2.T)(2.4)K(2.5)M=GNNM+1M.x_+a=G.(2.b)_左進(jìn)明等:廣義非線性方程的一個(gè)新的守恒差分格式ZYiabeZn+1n-1n+1n-1(u+u(u+u.jj)jj)=01式中-把,式并令:,的表達(dá)式代入(2.6)u表示u的共軛,nn+1nQ=(u2+u2)/2,nn-1則可得Q遞推之,即可得到(式.=Q,2.2)n+1n-1將(式與然后取實(shí)部部分得:2.1)(/2做內(nèi)積,xu-u)+=0.其中,n+1n-1n+1n-1=R(u-u(u-u,jj)jj)=02tj=0(2.7)J2n+12n-1n+1n-1=Rp(u+u(u-u=

6、xjxj)jj)2j=0n+1n-1n+1n-1-Rp(u+uu-u=xx,xx)2j=0n+12n-12-p(u/2,x-ux)JJJn+1n+1-2n-12=R(|uuuj|+|j|j|+2+2j=0n+12n-1-2n-1n+1n-1n+1n-1|uu+|u)(u+u(u-uj|j|j|jj)jj)=n+1+2n-1+2x(|u-|u.j|j|)+2j=0把,式中,并令:,代入(2.7)n+1+2n+2E=p(u+u)/2+x(|u+|u|),j|j+2j=0nn-1則可得到E遞推之,即可得到(式.=E,2.3)nn+1x2nx2JJO3差分格式的收斂性n當(dāng)在(式中u的近似是收斂的,為了

7、證明這一點(diǎn),定P0,P0時(shí),2.1)(,)xtxtj對(duì)ujn義如下算子:22jn+1jn-1jn+1jn-1xxQu(x,t)-p+jn2t23S-222(|u(x,t)|+|u(x,t)|u(x,t)|+|u(x,t)|Tjn+1jn+1jn-1jn+1+2-2|u(x,t)|+|u(x,t)|)jn-1jn-1Qu=iuu+|u|u=0.t-p由U展開式得:VWXYZjn+1jn-1,2S(3.1)(3.2)2(u(x,t)(t)+jntt2u(x,t)=u(x,t)+(u(x,t)t+jn+1jnjnt=2高校應(yīng)用數(shù)學(xué)學(xué)報(bào)>輯34(u(x,t)(t)+O(t),jnttt6u(x,

8、t)=u(x,t)-(u(x,t)t+jn-1jnjnt第21卷第1期2,t)(t)-(u(xjntt234(u(x,t)(t)+O(t),jnttt6u(x,t)=u(x,t)+(u(x,t)t+(u(x,t)x+j+1n+1jnjntjnx2(u(x,t)(t)+jntt222(u(x,t)(x)+(u(x,t)xt+(u(x,t)x(t)+jnxxjnxtjnxtt22233(u(x,t)(x)t+(u(x,t)(x)+(u(x,t)(t)+jnxxtjnxxxjnttt26644(O(x)+O(t),u(x,t)=u(x,t)+(u(x,t)t-(u(x,t)x+j-1n+1jnjnt

9、jnx2(u(x,t)(t)+jntt222(u(x,t)(x)-(u(x,t)xt-(u(x,t)x(t)+jnxxjnxtjnxtt22233,t)(x)t-(,t)(x)+(,t)(t)+(u(xu(xu(xjnxxtjnxxxjnttt26644(O(x)+O(t),t)=u(x,t)-(u(x,t)t+(u(x,t)x+u(xj+1n-1jnjntjnx2(u(x,t)(t)+jntt222(u(x,t)(x)-(u(x,t)xt+(u(x,t)x(t)-jnxxjnxtjnxtt22233(u(x,t)(x)t+(u(x,t)(x)-(u(x,t)(t)+jnxxtjnxxxjnt

10、tt26644(O(x)+O(t),u(x,t)=u(x,t)-(u(x,t)t-(u(x,t)x+j-1n-1jnjntjnx2(u(x,t)(t)+jntt222(u(x,t)(x)+(u(x,t)xt-(u(x,t)x(t)-jnxxjnxtjnxtt22233(u(x,t)(x)t-(u(x,t)(x)-(u(x,t)(t)+jnxxtjnxxxjnttt26644(O(x)+O(t).則差分格式的截?cái)嗾`差為:322244Lu(x,t)-Lu=u(t)-u(t)-u(x)+(O(x)+O(t).jntttxxttxxxx62123故在網(wǎng)格比確定的情況下,當(dāng)該問題是收斂的.0,0時(shí),L(

11、,)-L0,xtuxtujn<4差分格式的穩(wěn)定性構(gòu)造格式時(shí)需要證明計(jì)算過程中的誤差影響,如果計(jì)算過程的誤差保持有界或消失,就U左進(jìn)明等:廣義非線性R方程的一個(gè)新的守恒差分格式cSToViWXeTQ3說明它是穩(wěn)定的,否則就是不穩(wěn)定的.如果格式不穩(wěn)定,即使收斂性是滿足的,這個(gè)方案也是無用的.2讓r則(可轉(zhuǎn)化為:=t/(x),2.1)n+1n-1n+1n+1n+1n-1n-1n-1i(u-u=pr(uu+uu+u)-jj)j+1-2jj-1+uj+1-2jj-1n+1-2n-12n+12n-1-2n-1n+1n-1|uuuu+|u)(u+u,j|j|+|j|j|j|jj)即n+1(t)(|uj

12、|+2(4.1)i+2pr+n+1n+1-2n-12n+12n-1-2(t)(|uuuuu+j|+|j|j|+|j|j|+2n-1n+1)u+-i+2pr+|uj|jn+1n+1-2n-12(t)(|uuuj|+|j|j|+2n+12n-1-2n-1n-1n+1n+1n-1n-1|uu+|u)u-pr(u)=0.j|j|j|jj+1+uj-1+uj+1+uj-1nnixixj代換u并除去e得到=Vej,j(4.2)i+2pr+n+1n+1n-12n+12n-1-2(t)(|V|+|V|-2|V|+|V|V|+2n-1n+1|V|V+-i+2pr+n+1n+1-2n-12(t)(|V|+|V|V

13、|+2n+12n-1n-1n-1n+1n-1|V|-2+|V|V-2prcos(x)(V+V)=0.|V(4.3)n-1n+1對(duì)于項(xiàng)V從(式可以看出它們的系數(shù)是遞減的,由穩(wěn)定性的定義,4.3)V是有界的,可知該差分格式是穩(wěn)定的.HI數(shù)值分析式兩邊同乘以則可轉(zhuǎn)化為:(2.1),tn+1n-1n+1n+1n+1n-1n-1n-1i(u-u-pr(uu+uu+u)+jj)j+1-2jj-1+uj+1-2jj-1n+1(t)(|uj|+2n+1-2n-12n+12n-1-2n-1n+1n-1|uuuu+|u)(u+u=0,j|j|+|j|j|j|jj)那么可得到下面的線性方程組:n+1n+1n+1n+

14、1n+1n+1n+1+LJjuj+1+Kjujjuj-1=Mj.n+1其中,J=-pr,jn+1K=i+2pr+jn+1L=-pr,jn+1n-1M=prui-2pr-jj+1+(I.1)(I.2)n+1n+1-2n-12n+12n-1-2n-1(t)(|uuuuu+|u),j|+|j|j|+|j|j|j|+2n+1n+1-2n-12(t)(|uuuj|+|j|j|+2n+12n-1-2n-1n-1n-1|uu+|u)u+pru.j|j|j|jj-1將(式寫成矩陣的形式I.2)NO=P,(I.3)kZ其中,高校應(yīng)用數(shù)學(xué)學(xué)報(bào)l輯第21卷第1期1BA1C=T=BN-1A0.B2=A=CBNN

15、15;NTn+1n+1n+1n+1,u,dQ=uD=d1,N,1,N.11n+1n+1-2n-12到u然后用本文的差分格式來處理即可.對(duì)于項(xiàng)|,uuujj|+|j|j|+n+12n-1-2n-1n-1nn+1n+1nn-1我們可用u的外插對(duì)u|uu+|uu=2u,j|j|j|,j和ujj進(jìn)行賦值:jj+uj就可以得到理想的結(jié)果.LM數(shù)值算例現(xiàn)在我們將考慮如下算例:取u=-1,=2,=2.(,0)=Q(+10)2(+NOPRSTPRUVWP則該問題的孤立子解為:10),=-P=15,PXYu(P,t)=QRST(P+10-Zt)RUV2W(P+10)-3Wt.(M.1)取取式精確解差的范數(shù),列于

16、=0.1,=0.05,0.01兩種情況來計(jì)算數(shù)值解和(M.1)表1和表2其中包括文中提到的格式,文中提出的一種守恒(1-()2_abWScdQcb)的e格式和本文提出的守恒e對(duì)稱格式)-:TRRdRfRdgRfRVcWhTRRdRfRdgWUVcWh表1=0.1,=0.05nt-_abWScdQc-eTRRdRfRdgRfRVcWh2-eTRRdRfRdgWUVcWh左進(jìn)明等:廣義非線性S方程的一個(gè)新的守恒差分格式chrodinger表2h=0.1,=0.01nt85-CrankNicolson1-ThreelevelSevenpoint2-ThreelevelSixpoint從表中可以看出,

17、本文的格式優(yōu)于上面的C-rankNicolson格式和ThreelevelSeven格式.而且本文的格式是守恒的,對(duì)于非守恒的格式,在計(jì)算中容易出現(xiàn)非線性的"爆point2,5現(xiàn)象炸".參考文獻(xiàn):1T,Ab.AnahaTRlowitzMJalyticalandnumericalaspectsofcertainnonlinearevolution,.,1984,55:203-230.equationsIINumericalNonlinearSchrodingerEquationJJComputPhys2Z,P-.NuhangFeiérezGgarciaVM,Vdzq

18、uezemericalsimulationofnonlinearSchrodinger:.,1995,71:165-177.systemsanewconservativeschemeJApplMathComput3Z,.-hangTiandeCaofingPriceGg,etalApplicationsofthePRschemeforgeneralizednonlinear.,1999,14:SchrodingerequationsinsolvingsolitonsolutionsJApplMathJChinesehnivSeri21-29.4Z,C,.hanTiandeaofingPriceGg,etalApplicationsoftheweightedschemeforGNeSequations.,1998,5k3l:525-541.insolvingsolitonsolutionsJjoreanJComputApplMath5S-anzSernaJM,VerwerJG.Conservativeandnonconservativeschemesforthesolutionofthe.I,1986,6:25-42.nonlinearSchrmdingerequationJMAJ