非线形粘弹性流体稳定化有限元方法的研究.pdfVIP

( ) 2002 年 10 月 四川大学学报 自然科学版 Oct . 2002 第 39 卷第 5 期 Journal of Sichuan University (Natural Science Edition) Vol. 39 　No. 5 文章编号 (2002) 非线形粘弹性流体的稳定化有限元方法研究 1 2 潘　璐 ,游　雄 ( 1. 四川大学数学学院 ,成都 610064 ;2. 南京农业大学理学院 ,南京 210095) 摘要 :对遵循 Oldroyd B 型粘弹性动问题 ,提出了一种基于流线迎风 PetrowGalerkin 方法和最 小二乘法相结合的稳定化有限元方法. 这种方法有效地解决了以往粘弹性流体研究过程中出 现的由于应力方程对流控制占优引起的拟振动现象和有限元空间组合不匹配产生的不稳定现 ( ) 象. 近似应力 ,速度和压力分别是 Pk 连续的 , Pk +1 和 Pk 的 k ≥0 . 假定连续问题满足充分光 滑且小的解 ,则由不动点定理可得近似问题解的存在性及误差分析. 关键词 : 非线性形粘弹性流体 ; 稳定化有限元方法 ; 误差估计 中国分类号:O242. 21 　　　文献标识码 :A 1 　引言 我们的研究遵循 Oldroyd B 型非线形粘弹性流体模型的有限元方法. 其理论以 NS 方程理论作为基础 之一. 发生在粘弹性流中的现象主要与所考虑物质的记忆和伸缩效应有关. 在绝热情形 ,此问题的模型包括 质量守恒方程、动量守恒方程及液变基本方程 ,它们通常是非线性的[1～4 ] . ( ) 采用基于 SU PG 流线迎风 PetrovGalerkin 方法 和最小二乘法相结合的稳定化有限元方法 ,这种方法 有效地解决了以往粘弹性流体研究过程中出现的由于应力方程对流控制占优引起的拟振动现象和有限元空 间组合不匹配而出现的不稳定现象. 我们在现有文献基础之上首次将有限元空间的应用推广到了任意阶 ,即 解决了 P 连续应力, P 连续速度, P 连续压力, k ≥0 的情况. 假定连续问题满足充分光滑且小的解, 则 k k +1 k 由不动点定理可得近似问题解的存在性及误差分析. 2 　问题的提出 我们考虑 R 2 中具有Lipschizian 边界 г的有界连通开集Ω的流体. 使用以下记号 :对标量函数 p , 其梯度 p , ( p) i 记为p , i , 对向量函数 u , 其梯度张量 u , ( u) ij = ui , j , 对向量函数 u , 其散度是标量 . u = d σ σ ( σ) σ ui , i , 且令 u . = ui , 对张量函数 , 其散度是向量 . , . i = ij , j . dx i σ λ( ) σ λβ( σ ) α ( ) Ω 问题 P : + u . + , u - 2 d u = 0 , 　in