耿發(fā)展

（常熟理工學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，江蘇 常熟 215500）

奇異攝動問題在流體動力學(xué)、海洋大氣層環(huán)流、化學(xué)反應(yīng)、最優(yōu)化控制等領(lǐng)域有著廣泛的應(yīng)用，因此，該問題引起了廣大研究人員的興趣.奇異攝動問題的解具有邊界層的特性，經(jīng)典的數(shù)值方法難以獲得該類問題的精度較高的近似解.具有轉(zhuǎn)向點(diǎn)的奇異攝動問題數(shù)值求解難度更大，尋求這類問題的有效數(shù)值方法更具挑戰(zhàn)性.文獻(xiàn)[1-4]討論了這類問題的數(shù)值方法，本文基于再生核空間理論，提出求解具有兩個(gè)邊界層的奇異攝動轉(zhuǎn)向點(diǎn)問題的有效數(shù)值方法.

1 數(shù)值方法

考慮下面的奇異攝動轉(zhuǎn)向點(diǎn)問題：

其中函數(shù)a(x)，b(x)和 f(x)是在-1＜x＜1內(nèi)充分光滑.

依據(jù)a(x)的取值，方程（1.1）的解表現(xiàn)出邊界層或內(nèi)部層的特性.系數(shù)a(x)=0時(shí)點(diǎn)稱為轉(zhuǎn)向點(diǎn).與非轉(zhuǎn)向點(diǎn)問題相比，具有邊界層或內(nèi)部層的轉(zhuǎn)向點(diǎn)問題處理起來往往更加棘手.這里，我們主要考慮如何解決具有邊界層的奇異攝動轉(zhuǎn)向點(diǎn)問題的情況.

假設(shè)：

在此假設(shè)下，問題（1.1）具有兩個(gè)邊界層.

1.1 問題（1.1）的轉(zhuǎn)化

引入變量變換 x=h(s)=sin(π sin(π sin(πs/2)/2)/2)，

問題（1.1）轉(zhuǎn)化成了下面的新的邊值問題

問題（1.3）的解不再具有邊界層特性，處理起來比（1.1）容易得多，一旦得到（1.3）的解，（1.1）的解即可得到，下面我們將用再生核方法求解（1.3）.

1.2 求解問題（1.3）的再生核方法

為了解決上述問題，引入一個(gè)未知函數(shù)

其中 ?(s)=γ0+γ1s，滿足 ?(-1)=α，?(1)=γ.

定義 Lv(s)= εv″(s)+p(s)v′(s)+q(s)v(s)，問題（1.3）轉(zhuǎn)化為

其中 F(s)=g(s)-L?(s).

為了求解問題（1.5），我們構(gòu)造下面的再生核空間W4[-1,1].

定義1.1 W4[-1,1]={u(x)|u?(x)是絕對連續(xù)函數(shù)，且u(4)(x)∈L2[-1,1],u(-1)=0, }u(1)=0，W4[-1,1]的內(nèi)積和范數(shù)分別表示為

和

定理1.1 W4[-1,1]為再生核空間，它的再生核表示為

定理1.1的證明及再生核的求解方法參見文獻(xiàn)[3，5].

不難證明，L:W4[-1,1]→W1[-1,1]為有界線性算子，從而其共軛算子 L*存在.記 φi(s)=kˉ(s,si)和ψi(s)=L*φi(s)，對實(shí)施施密特正交化，即得

定理1.2 對于問題（1.5），如果在區(qū)間[-1,1]內(nèi)是稠密的，那么為空間W4[-1,1]的完全系.

定理的證明過程見文獻(xiàn)[6-7].

定理1.3 如果在區(qū)間[-1,1]內(nèi)稠密，那么問題（1.5）的解可表示為

證明根據(jù)定理1.2，很顯然有：是再生核空間W4[-1,1]的完全正交基.

進(jìn)一步，我們可以得到

命題得證.

對級數(shù)表示（1.6）進(jìn)行N項(xiàng)截?cái)啵憧傻脝栴}（1.5）的近似解.

定理1.4 問題（1.5）的近似解yN(x)及其導(dǎo)數(shù)y′N(s)一致收斂.

證明由于W4[lε,rε]是一個(gè)希爾伯特空間，顯然有

由于

從而可得

命題得證.

由（1.4）及變量變換 x=h(s)，可得（1.1）的近似解 uN(s)=vN(h-1(x)),其中 vN(s)=yN(s)+?(s).

2 數(shù)值算例

考慮下面的奇異攝動問題

f(x)已知且保證問題的精確解是u(x)=e-2x(2x-1)/ε.使用本文提出的方法，取ε=2-8,2-12.本例所有的計(jì)算通過Mathemati?ca7.0完 成 ，所得的數(shù)值結(jié)果見圖1，圖2.

圖1 ε=2-8時(shí)的絕對誤差

圖2 ε=2-12時(shí)的絕對誤差

[1]Rai P,Sharma K K.Numerical method for singularly perturbed differential-difference equation with turning point[J].Int J Pure Ap?pl Math,2011,73(20):451-470.

[2]Kadalbajoo M K,Arora P,Gupta V.Collocation method using artificial viscosity for solving singularly perturbed turning point prob?lem having twin boundary layers[J].Comput Math Appl,2011,61:1595-1607.

[3]Geng F Z,Qian S P.Reproducing kernel method for singularly perturbed turning point problems having twin boundary layers[J].Ap?plied Mathematics Letters,2013,26:998-1004.

[4]Geng F Z,Qian S P,Li S.A numerical method for singularly perturbed turning point problems with an interior layer[J].Journal of Computational and Applied Mathematics,2013,255:97-105.

[5]Cui M G,Lin Y Z.Nonlinear Numerical Analysis in Reproducing Kernel Space[M].New York:Nova Science Pub Inc,2009:15-50.

[6]Geng F Z,Cui M G.Solving a nonlinear system of second order boundary value problems[J].J Math Anal Appl,2007,327:1167-1181.

[7]Geng F Z,Cui M G.A novel method for solving a class of singularly perturbed boundary value problems based on reproducing ker?nel method[J].Appl Math Comput,2011,218:4211-4215.