定常粘弹性流体的最小二乘混合有限元.pdfVIP

第25卷第8期 怀化 学院学报 Vo1．25．No．8 2006年 8月 JOURNALOFHUAIHUA UNIVERSITY Aug．，2006 定常粘弹性流体的最小二乘混合有限元 张宏伟 ， 鲁祖亮 ， 蔡海涛2 (1．长沙理工大学 数学与计算科学学院， 湖南 长沙 410077； 2．中南大学 数学科学与计算技术学院， 湖南 长沙 410075) 摘 要：针对于定常的服从 OldroydB型本构律的粘弹性流体流动建立 了一种最小二乘混合有限元方法 ．应力、 速度分别采用不连续分片k多次项式 ，连续分片 +1次多项式P+l( O)．分析了逼近问题的解的存在性并给出了 逼近解的误差估计 ． 关键词 ：粘弹性流体流动 ； 最小二乘混合有限元 ； 一致三角剖分 中图分类号：0241．5 文献标识码 ：A 文章编号：1671—9743(2006)O8～0017—05 Dirichlet问题的最小二乘方法是由Bramble和 Nitsehe[1J最早开始研究 的，Raviart、Thomas和 Brezzi[zJ， Carey和 Oden-3发展了这一方法 ．近年来 ，最小二乘混合有限元方法引起了国内外学者 的广泛兴趣 ，Lazarov， Pehlivanov_4J Ewing，Wang，Brandts和我 国的羊丹平_5J、罗振东_6J、黄云清和陈艳萍l7—91都对最小二乘混合 ， 有限元的发展做了卓有成效的工作 ．粘弹性流体问题一直是有限元研究的一个热点问题 ．在许多领域，如油 藏稠油渗流、二维漕道的流动、注入过程和挤 出过程 ，粘弹性流体的数值模拟 J ]都是其 中的一个重要研 究课题 ．我们首次将最4-乘混合有限元方法应用于求解定常的服从 OldroydB型本构律的粘弹性流体流动 问题 ，分析 了解 的存在性并给出了逼近解 的误差估计 ． 1 基本概念和定义 考虑在 n ∈兄上边界为 r的粘弹性流体流动模型的一般形式(Oldroyd问题)[10-12]： 口+0+ (¨·7)口+ g(口，7¨)一2ad(¨)=0，在Q ×[0，T]上 (1．1) 12,一7 ·O一2(1一口)7 ·d(21,)+7P=f，在Q×[0，T]上 (1．2) 7 ·12, 0，在Q ×[0，T]上，21,=0，在 ／-,×[0，T]上 (1．3) 当￡=0时，“= “0，= 0．为弹性应力张量，它是对称的即 = ，“是流速矢量，d(“)=÷(v“+ v “)为变形变率张量，p是流体压力， 0为Weissenberg常数，a为延滞参数(当a：0时，问题退化为协旋 转 Maxwell问题 ，当 = 1时，问题退化为对流 Maxwell问题)．这两类方程已有大量的文献专 门研究 ，因此本文 只考虑0口1的情形 ．g ：R ×R 一 R 是双线性映射： g。(口，7¨)：Or．t／(¨)一 (¨)口一0(d(¨)口+口d(¨))+ (口v¨+v¨ )一 (7 + 7¨)， 其中 一101，叫(¨)= 1(v ¨一v ¨)为扭曲张量 ． 本文中，我们使用标准 的Sob。lev空间 Wm(n)，其 中 ll ll， = ∑ llD ll(n)．当n=2时， H (n)=Wm (n)且 ll·ll Il·ll，2，Il·lI= ll·ll。，2． 定义 1 若三角形剖分族均为拟一致剖分，且任何两个相邻单元构成h．近似平行四边形即存在与h无关 收稿 日期 ：2006—06—07 基金项 目：湖南教育厅青年项 目资助 (项 目编号：05B022)． 作者简介 ：张宏伟 (1966一)，男，山东惠民人 ，长沙理工大学副教授 ，硕士，主要研究偏微分方程数值解 · 18 · 怀化 学院学报