李 恬，王彩華，鄭尚昆

（天津師范大學(xué) 數(shù)學(xué)科學(xué)學(xué)院，天津300387）

二維奇異擾動問題的非等距有限差分格式

李 恬，王彩華，鄭尚昆

（天津師范大學(xué) 數(shù)學(xué)科學(xué)學(xué)院，天津300387）

通過對一維非等距中心差分格式引入擬合因子，構(gòu)造了一類新型非等距中心差分格式，將其推廣到二維情形，得到一類針對二維奇異擾動問題的新型非等距五點差分格式，對該格式進(jìn)行了截斷誤差估計.數(shù)值實驗部分采用4種非等距網(wǎng)格進(jìn)行處理，結(jié)果表明該非等距差分格式對含邊界層的奇異擾動問題有很好的實用性.

奇異擾動問題；非等距網(wǎng)格；有限差分格式；擬合因子；截斷誤差估計

奇異擾動問題的解具有多尺度特性，由于擾動參數(shù)的引入，導(dǎo)致解在局部區(qū)域的改變很劇烈，而在其他區(qū)域平穩(wěn)變化，即出現(xiàn)邊界層性狀.Pearson［1］首先嘗試將有限差分格式結(jié)合網(wǎng)格寬度來解決這類問題，之后，這種方法被Abrahamsson［2］做了深入探討，闡明奇異擾動問題差分方法的精度不會超過O（ε）.許多學(xué)者進(jìn)一步對奇異擾動問題進(jìn)行了研究.文獻(xiàn)［3］和文獻(xiàn)［4］針對奇異擾動問題得到了精度為O（h）和O（h2）的差分近似；文獻(xiàn) ［5］改進(jìn)了對流擴(kuò)散方程的擾動差分方法，該方法對邊界層高Reynolds數(shù)效應(yīng)有極高的分辨能力；文獻(xiàn)［6］對于一維奇異擾動邊值問題描述了10種產(chǎn)生Il′in差分格式的方法，并分析了相關(guān)方法的特點；文獻(xiàn)［7］建立了二維對流擴(kuò)散方程的變步長擾動有限差分格式，該格式具有精度高、穩(wěn)定性與收斂性好的特點；文獻(xiàn)［8］提出了非均勻網(wǎng)格上一維對流擴(kuò)散方程的高階緊致差分方法；文獻(xiàn)［9］對一維定常對流擴(kuò)散方程構(gòu)造了二階精度的非等距網(wǎng)格差分格式，并給出格式的截斷誤差.對于二維奇異擾動問題，文獻(xiàn)［10-13］給出了利用Taylor級數(shù)展式構(gòu)造高精度指數(shù)型五點差分格式的方法；文獻(xiàn)［14］給出了二維對流擴(kuò)散方程的緊致差分格式；文獻(xiàn)［15］研究了用正交梯度網(wǎng)格構(gòu)造二維高精度有限差分格式的方法.上述二維問題的差分格式都是基于等距網(wǎng)格的.在實際的流場計算中，非等距網(wǎng)格可以在邊界層進(jìn)行細(xì)網(wǎng)格處理，因而相比采用相同結(jié)點數(shù)的等距網(wǎng)格，可以得到更精確的數(shù)值結(jié)果，且能反映出邊界層變化.本研究對于一維非等距中心差分格式，通過增加擬合因子，結(jié)合系數(shù)常數(shù)化，構(gòu)造了一種指數(shù)型差分格式，并將此差分格式推廣到二維情形，得到一類針對二維奇異擾動問題的新型非等距五點差分格式，對該格式進(jìn)行了截斷誤差估計.數(shù)值實驗在4種不同的非等距網(wǎng)格上進(jìn)行，實驗結(jié)果表明本研究格式對于邊界層問題具有很好的實用性.

1 一維非等距差分格式的構(gòu)造

考慮區(qū)間I=［0，1］上一維奇異擾動邊值問題

其中：u=u（x）為待求函數(shù)；ε為一正常數(shù)，0＜ε?1；a、b為給定的邊界值；p（x）為對流項系數(shù)，f（x）為源匯項，且p（x）、f（x）充分光滑.一般地，當(dāng)p（x）＜0時，邊界層在x=0附近；當(dāng)p（x）＞0時，邊界層在x=1附近.

將區(qū)間I=［0，1］剖分為N份，用xi表示節(jié)點坐標(biāo)，設(shè)0=x0＜x1＜…＜xN-1＜xN=1，定義h=xi+1-xi，，并簡記.在等距情形下有.

記u0，u1，…，uN為方程（1）在剖分節(jié)點的精確值，U0，U1，…，UN為待構(gòu)造的離散差分格式在剖分節(jié)點的計算值.

一維奇異擾動問題（1）在常系數(shù)情形下的精確解為

對等距中心差分格式引入擬合因子σi，設(shè)差分格式在節(jié)點處滿足

將精確解代入，求得

其中：ρi=pih/ε，σi=σ（ρi）.即得到Il’in差分格式

其中：

針對非等距網(wǎng)格引入擬合因子σi，設(shè)差分格式滿足

則離散方程在內(nèi)點相應(yīng)的差分格式為

將精確解

代入式（7），可得擬合因子

故可得新的非等距中心差分格式

其中：

2 二維非等距五點差分格式

考慮二維定常奇異擾動問題

求解區(qū)域為Ω=｛（x，y）|a＜x＜b，c＜y＜d｝，函數(shù)p、q、f、γ在Ω內(nèi)充分光滑.對于五點差分格式，記（xi，yj）為0點，（xi+1，yj）為1點，（xi-1，yj）為3點，（xi，yj+1）為2點，（xi，yj-1）為4點，如圖1所示.

圖1 矩形網(wǎng)格的五點Fig.1 Five spots of rectangular mesh

二維奇異擾動問題（9）可分解為如下2個方程

對其在x、y方向上應(yīng)用差分格式，引入擬合因子σh、σk，則得x、y方向上的一維非等距格式

結(jié)合以上2式可得非等距五點新差分格式

其中：

下面考慮非等距五點差分格式（13）的截斷誤差.對格式（13），利用Taylor級數(shù)展開式可得

于是

3 數(shù)值算例

對常系數(shù)和變系數(shù)帶邊界層的奇異擾動問題進(jìn)行數(shù)值實驗，利用Gauss-Seidel方法進(jìn)行計算機(jī)求解，為便于比較誤差，所給問題已知精確解.

實驗在4種非等距網(wǎng)格情形下進(jìn)行.

（1）伸縮變換網(wǎng)格1.使用如下網(wǎng)格生成函數(shù)剖分［0，1］區(qū)間

其中：參數(shù)λ為網(wǎng)格伸縮系數(shù)，滿足|λ|≤1.當(dāng)0＜λ≤1時，網(wǎng)格點隨著λ值的增大集中分布在x=1附近；當(dāng)-1≤λ＜0時，網(wǎng)格點隨著λ值的減小集中分布在x= 0附近；當(dāng)λ=0時，網(wǎng)格點均勻分布在［0，1］中.

（2）伸縮變換網(wǎng)格2.使用如下網(wǎng)格生成函數(shù)剖分［0，1］區(qū)間

其中：參數(shù)c為網(wǎng)格伸縮系數(shù)，滿足c＞0.當(dāng)c較大時，網(wǎng)格點集中分布在x=1附近；當(dāng)c→0時，網(wǎng)格點相對均勻地分布在［0，1］中.

（3）等比網(wǎng)格.當(dāng)公比r＜1時，得到逐漸稠密的網(wǎng)格；當(dāng)公比r＞1時，得到逐漸稀疏的網(wǎng)格.

（4）分塊均勻網(wǎng)格.將區(qū)間I=［0，1］均勻剖分為n份，然后再將其中1個或幾個子區(qū)間剖分為m等份.

例1 考慮二維常系數(shù)對流擴(kuò)散方程

該問題的精確解為

計算區(qū)域為（x，y）∈［0，1］×［0，1］，邊界為

圖2給出了當(dāng)Re=10、102、103、104時的精確解的曲面圖.表1給出了4種不同非等距網(wǎng)格下數(shù)值解與精確解的最大絕對誤差.前3種網(wǎng)格區(qū)域進(jìn)行64×64的非等距剖分，方程為齊次情形，當(dāng)Re=1時沒有出現(xiàn)邊界層，隨著Re的增大在x=1和y=1處出現(xiàn)邊界層.分塊等距網(wǎng)格中將區(qū)間先等距剖分51份，然后將最后一個區(qū)間等距剖分14份.

圖2 當(dāng)Re=10、102、103、104時精確解曲面圖（例1）Fig.2 Exact solutions for Re=10，102，103，104（Example 1）

表1 4種非等距網(wǎng)格下的精確解與數(shù)值解的最大絕對誤差（例1）Tab.1 Maximum absolute errors of exact solutions and numerical solutions on four kinds of non-uniform meshes（Example 1）

例2 考慮變系數(shù)二維對流擴(kuò)散方程

該問題的精確解為

計算區(qū)域為（x，y）∈［0，1］×［0，1］，邊界為

圖3給出了當(dāng)P=10、102、103、104時精確解的曲面圖.表2給出了4種不同非等距網(wǎng)格下數(shù)值解與精確解的最大絕對誤差.前3種網(wǎng)格區(qū)域進(jìn)行64×64的非等距剖分，方程為非齊次情形，當(dāng)P=1時沒有出現(xiàn)邊界層，隨著P的增大在x=1和y=1處出現(xiàn)邊界層.分塊等距網(wǎng)格中將區(qū)間先等距剖分51份，然后將最后一個區(qū)間等距剖分14份.

表2 4種非等距網(wǎng)格下的精確解與數(shù)值解的最大絕對誤差（例2）Tab.2 Maximum absolute errors of exact solutions and numerical solutions on four kinds of non-uniform meshes（Example 2）

由表1和表2可知，對于不同的奇異擾動問題，本研究格式對含小邊界層模型的數(shù)值結(jié)果比較理想，隨著參數(shù)的增大，數(shù)值解和精確解的誤差逐漸變小，當(dāng)網(wǎng)格點在邊界層內(nèi)密集時，該格式能反映解在邊界層內(nèi)的變化，因此本研究的非等距五點差分格式對邊界層問題具有很好的實用性.

圖3 當(dāng)P=10、102、103、104時精確解曲面圖（例2）Fig.3 Exact solutions for P=10，102，103，104（Example 2）

［1］PEARSON C E.On a differential equation of boundary layer type［J］.Journal of Mathematics and Physics，1968，47：134-154.

［2］ABRAHAMSSON L R，KELLER H B，KREISS H O.Difference approximations for singular perturbations of systems of ordinary differential equations［J］.Numerische Mathematik，1974，22（5）：367-391.

［3］KELLOGG R B，TSAN A.Analysis of some difference approximations for a singular perturbation problem without turning points［J］.Mathematics of Computation，1978，32（144）：1025-1039.

［4］BERGER A E，SOLOMON J M，CIMENT M.Uniformly accurate difference methods for a singular perturbation problem［J］.Computational and Asymptotic Methods，1980（3）：14-28.

［5］陳國謙，楊志峰.對流擴(kuò)散方程的指數(shù)型攝動差分法［J］.計算物理，1993，10（2）：197-207.CHEN G Q，YANG Z F.A perturbational finite difference scheme of exponential type for the convection-diffusion equation［J］.Chinese Journal of Computational Physics，1993，10（2）：197-207（in Chinese）

［6］ROOS H G.Ten ways to generate the Il＇in and related schemes［J］.Journal of Computational and Applied Mathematics，1994，53（1）：43-59.

［7］李桂波，李明軍，高智.對流擴(kuò)散方程的變步長攝動有限差分格式［J］.水動力學(xué)研究與進(jìn)展：A輯，2005，20（3）：293-299.LI G B，LI M J，GAO Z.The perturbational finite difference（PFD）scheme of variational step of convection-diffusion equation［J］.Journal of Hydrodynamics，2005，20（3）：293-299（in Chinese）.

［8］田芳，田振夫.非均勻網(wǎng)格上求解對流擴(kuò)散問題的高階緊致差分方法［J］.寧夏大學(xué)學(xué)報：自然科學(xué)版，2009，30（3）：209-212.TIAN F，TIAN Z F.A high order compact finite difference method for convection diffusion problems on non uniform grids［J］.Journal of Ningxia University：Natural Science Edition，2009，30（3）：209-212（in Chinese）.

［9］曹廣滿，王彩華，齊海濤.對流擴(kuò)散方程的非一致網(wǎng)格有限差分方法［J］.天津師范大學(xué)學(xué)報：自然科學(xué)版，2010，30（1）：7-10.CAO G M，WANG C H，QI H T.Non-uniform finite difference scheme method for convection-diffusion equations［J］.Journal of Tianjin Normal University：Natural Science Edition，2010，30（1）：7-10（in Chinese）.

［10］TIAN Z F，DAI S Q.High-order compact exponential finite difference methods for convection-diffusion type problems［J］.Journal of Computational Physics，2007，220（2）：952-974.

［11］WANG C H.A new way to generate an exponential finite difference scheme for 2D convection-diffusion equations［J］.Journal of Applied Mathematics，2014（5）：63-76.

［12］SANYASIRAJU Y，MISHRA N.Exponential compact higher order scheme for nonlinear steady convection-diffusion equations［J］.Communications in Computational Physics，2011，9（4）：897-916.

［13］RADHAKRISHNA PILLAI A C.Fourth-order exponential finite difference methods for boundary value problems of convective diffusion type［J］.International Journal for Numerical Methods in Fluids，2001，37（1）：87-106.

［14］王峰峰，王彩華，齊海濤.二維對流擴(kuò)散方程的緊致差分格式［J］.天津師范大學(xué)學(xué)報：自然科學(xué)版，2011，31（1）：11-16.WANG F F，WANG C H，QI H T.Compact difference scheme for twodimentional convection-diffusion equations［J］.Journal of Tianjin Normal University：Natural Science Edition，2011，31（1）：11-16（in Chinese）.

［15］SPOTZ W F.Formulation and experiments with high-order compact schemes for nonuniform grids［J］.International Journal of Numerical Methods for Heat&Fluid Flow，1998，8（3）：288-303.

（責(zé)任編校 馬新光）

Finite difference scheme on non-uniform mesh for two-dimensional singularly perturbed problems

LI Tian，WANG Caihua，ZHENG Shangkun
（College of Mathematical Science，Tianjin Normal University，Tianjin 300387，China）

A new five-point finite difference scheme on non-uniform mesh by adding fitting factors to the central difference scheme is presented for two-dimensional singularly perturbed problems，and the truncation error estimate is given.Numerical experiments are implemented on four kinds of non-uniform mesh，and the results show that the scheme is well applicable to singularly perturbed problems with boundary layers.

singularly perturbed problems；non-uniform mesh；finite difference scheme；fitting factor；truncation error estimate

O241.82

A

1671-1114（2016）03-0008-05

2015-10-15

國家自然科學(xué)基金資助項目（11071123）.

李 恬（1992—），女，碩士研究生.

王彩華（1973—），女，副教授，主要從事微分方程數(shù)值解方面的研究.