非线性对流扩散问题的交替方向差分流线扩散法.pdfVIP

高 等 学 校 计 算 数 学 学 报 第29卷第 4期 非线性对流扩散问题的交替方向 差分流线扩散法 张阳 (南开大学数学科学学院，核心数学与组合数学实验室(LPMC)，天津 300071) AI ERN TING．DIRECTIoN DIFFERENCE STREAMLINE DIFFUSIoN METHoD FoR NoNLINEAR CoNVECTIoN．DoMIN TED DIFFUSIoN PRoBLEMS Zhang Yang (School of Mathematics and LPMC，Nankai University，Tianjin 300071，China) Abstract Alternating—Direction method iS combined with finite difierence streamline—diffusion(FDSD)method to create the Alternating—Direction FDSD procedures for two-dimensional nonlinear convection—dominated diffusion equa- tion．Realization of the scheme and the analysis of error are presented．This scheme not only realizes the purpose of lowering the dimension of solving the two-dimension nonlinear convection diffusion equation，but also keeps the favorable stability and high precision of FDSD method．Finially,numerical examples are presented to clar— ify the method． Key words alternating—direction method，finite difference streamline diffusion method，error estimate． AMS(2000)SUbject classifications 65M12，65M60 中图法分类号 0241．82 国家自然科学基金(103016)资助项目 2007年 12月 高 等学校计算 数学学报 367 1引言 Peaceman，Douglas等人于1955年提出了差分格式的交替方向法．随后，Douglas， Dupont于1972年又提出了有限元格式的交替方向法 [1】_其基本思想是：对两个或三个空 间变量的二阶抛物型和双曲型问题，将交替方向法与Galerkin方法相结合，通过算子分裂 技术，把高维问题转化为一系列低维问题，交替地沿各空间变量的方向求解．[2】、[3】和[4】 给出了对更一般扩散问题(带对流项的抛物方程)的数值求解和误差分析． 流线扩散法(以下简称SD方法)是由Hughes和Johnson等人于1982年提出的一 种数值求解对流占优扩散问题的有限元方法 [6】_它以其良好的稳定性和高阶精度而为 数学界和工程界所关注．1998年，孙澈教授等在此基础上提出了有限差分流线扩散法(简 称FDSD)[ 】l【8]．与传统的sD法相比，FDSD法计算更简便且具有良好的稳定性和高阶精 度．但对于高维问题，每一时间层都要解一个大型线性代数方程组．本文针对二维非线性 对流占优扩散问题，将交替方向法与FDSD法相结合，构造了交替方向有限差分流线扩散 法(AD—FDSD)格式，并对其进行了理论分析．此格式不但减少了计算的工作量且仍保留 了FDSD法良好的高阶精度性质． 设Q∈R 为具有逐片光滑边界 Q的有界区域，用W ，p(Q)， ’p(