抛物型方程非齐次边值问题的推广型LOD有限差分及有限元格式.pdfVIP

维普资讯 2006年 6月 高 等 学 校 计 算 数 学 学 报 第 28卷第 2期 抛物型方程非齐次边值问题的推广型 LOD有限差分及有限元格式 王彩华 王同科 (天津师范大学数学科学学院，天津300384) EXTENDED LoCALI oNE—DIM ENSIoNAL FIN ITE DIFFERENCE AND FIN ITE ELEM ENT SCHEM ES FoR NoNHoM oGENEoUS PARABoLIC DIFFERENTIALEQUATIONSW ITH NoNHoM oGENEoUSBoUNDARY CoNDITIoNS W angCaihua WangTongke (SchoolofMathematicalScience，TianjinNormalUniversity,Tianjin300384) Abstract Thispaperpresentssomeextendedlocallyone—dimensionalfinitedif- ferenceandfiniteelementschemesfornonhomogeneousparabolicdifferentialequa- tionswith nonhomogeneousboundary conditions．ItiSproved thatthesefinite differenceschemeshavesecondorderaccuracywithrespecttodiscreteL2norm． Finally,theschemesareillustratedbynumericalexamplesandtheresultsarevery satisfactory． Key words twoand threedimensionalparabolicdifferentialequation，nonh~ mogenerousboundarycondition，finitedifferencescheme，finiteelementmethod， extendedlocallyone—dimensiona lscheme．errorestimate． AMS(2000)subjectclassifications 65M06，65M25 中图法分类号 O241．82 国家自然科学基金项目资助． 收稿 日期：2004—03—01． 维普资讯 2006年 6月 高 等 学 校计 算 数 学 学 报 ．139． 1引 言 本文考虑区域Q=0【，1】d(d=2，3)上的非齐次抛物型方程第一边值问题 Ou — Ca△u+ u： ，￡)， ∈Q，￡∈(0，】， (1．1) 瓦 t(‘，0)=t正o()， ∈Q， (1．2) t(‘，t)= (，￡)， ∈0fl，t∈(0，】， (1．3) 其中Cl， 为常数且 Cl0，C2 0． 对于以上问题，可以使用有限差分方法及有限元方法进行离散，并采用交替方向方法 求解．交替方向方法能够将高维问题转化为一系列的一维问题进行计算，具有计算量少，计 算稳定且易于并行实现等优点，在大规模科学计算中起着非常重要的作用，一直是计算数 学研究的热点问题 (1【—10】)．关于交替方向差分方法，二维情形主要有P—R格式，Douglas格 式及局部一维 (LOD)格式