对流-扩散-反应方程的变分多尺度解法.pdfVIP

第26卷 ~PJ6JtJ] 工 程 数 学 学 报 Vo1．26No．6 2009q-12)] CHINESEJOURNALOFENGINEERINGMATHEMATICS Dec．2009 文章编~-：1005—3085(2009)06—0997—08 对流一扩散一反应方程的变分多尺度解法术 朱海涛， 欧阳洁 f西北工业大学应用数学系，西安 710072) 摘 要：根据变分多尺度的思想求解了对流项和反应项占优的对流一扩散一反应方程。在变分多尺度思想 的理论框架内，推导了附加于Galerkin变分弱形式的稳定化结构和具体的稳定化系数；阐述了 这种稳定化结构和经典的SUPG稳定化结构之问的关系：数值算例表明，该稳定化系数可以适 应均匀和非均匀的计算网格。通过网格的恰当加密，变分多尺度方法消除了算例中的数值伪振 荡 。 关键词：对流一扩散一反应方程；变分多尺度方法；SUPG；稳定化方法 分类号：AMS(2000180M25 中图分类号：O241 文献标识码：A 1 引言 用标准Galerkin方法求解对流一扩散 一反应方程 ，当对流项或反应项 占优时，数值结 果会伴随着剧烈的数值伪振荡 。通过在Galerkin变分弱形式的左端加上一项基于原偏微 分方程残差的稳定项而产生的稳定化类方法 ，既使该类数值计算方法获得 了相容性 ，又 获得 了稳 定性 1【】，从而保证 了稳定化类方法的收敛性 。这类方法包括：SUPG(Streamline Upwind／Patrov—Galerkin)、GLS(GalerkinLeastSquares)、USFEM(UnusualStabilizedFinite ElementMethods1等。但是，这类方法的缺点是最优的稳定化系数很难确定。上世纪90年 代，Hughes考虑了各种稳定化方法，提出了变分多尺度方法的基本理论框架 2【，3j。本文根据变 分多尺度的思想求解了对流项和反应项 占优的对流一扩散一反应方程，主要的工作是：1) 在 变分多尺度的理论框架内，用泡函数近似 “细”尺度上的解，推导了用变分多尺度方法求解该 方程的附加于Galerkin变分弱形式的稳定化结构及对应的稳定化系数：2) 阐述了在线性和双 线性的情形下，变分多尺度方法的稳定化结构和经典的SUPG稳定化结构之间的关系；3) 基 于变分多尺度思想求解了两个具有代表性的数值算例，并通过网格加密技术与变分多尺度方法 的耦合，消除了算例中由对流项和反应项 占优引起的数值伪振荡。 2 对流．扩散．反应方程 考虑如下的对流一扩散一反应方程 a·Vu—kAu+8U=f 在 内 (1) 乱=0 在 Q 上 收稿 日期：2008-05—09． 作者简介：朱海涛 (1982年4月生)，男，硕士．研究方向：数值仿真 基金项目：国家自然科学基金 ；国家重点基础研究发展计划 (2o05CB321704)． 998 1： 程 数 学 学 报 第26卷 其叶1a为已知的对流速度； 是扩敞系数，k 0；s是已知的反应项系数；f是 知的源项。 (，为求解 域 2【的I规剂分，K ∈^是 If1的任 有限单元，hK是 元K 的直 Inx{ }， ={h∈硎(Q)I^IK∈p1()，V ∈ )， ∈ P (K)是单元K上连续的线性或双线性多项式函数空间。 方程(1)对应的标准Galerkin离散格式为：求 “h∈ C础 (Q)，满足 (a。VUh，)Q+k(Vuh，Vvh)n+(suh，Vh)n=(f，^)2【，VVh∈ C 】喏， (2) 这里 (，u)= uvdx，除非特别声明，D=Q。当 《