- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
多孔介质中可压缩可混溶驱动问题的共轭梯度迭代法的L^2模误差估计
维普资讯
2006年3月 高 等 学 校 计 算 数 学 学 报 第28卷第 1期
多孔介质中可压缩可混溶驱动问题的
共轭梯度迭代法的 2模误差估计
马克颖
(山东大学数学与系统科学学院,济南250100)
L2.ERROR ESTIMATEOFCONJUGATE
GRADIENT ITER觚 ’ION M ETH0D F0R M ISCIBLE
C0M PRESSIBLE DISPLACEM ENT PR0BLEM
MaKeying
(SchoolofMath.andSystemSciences,ShandongUniv.,Jinan250100)
Abstract Misciblecompressibledisplacementinporousmediaismodelledbya
nonlinearcoupledsystem oftwoparabolicequations:thepressureequationand
theconcentrationequation.ThepressureisapproximatedbyastandradGalerkin
methodwhiletheconcentrationisapproximatedby acombinationofaGalerkin
methodnadthemethodofcharacteristics.Aconjugategradientiterationmethod
isusedtosolvethealgebraicproblemsarisingfrom thosetwoapproximateschemes.
Bydetailedtheoreticalanalyses,L2~errorestimateisobtainedbetweentheexactSO-
lutionoforiginalproblemandthesolutionofconjugategradientiterationmethod.
Keywords misciblecompressibledisplacement,conjugategradientiteration
method,errorestimate.
AMS(2000)subjectclassifications65M10,65N30
中图法分类号 O241.82
掌 国家自然科学基金项目,国家自然科学基金项 目19972023),
教育部博士点基金资助项目(20030422047).
收稿 日期:2003.10.22.
维普资讯
2006年 3月 高等学校计算数学学报 ·51·
1数学模型
多孔介质中可压缩可混溶驱动问题的模型是两个非线性抛物型方程:压力方程和饱
和度方程.Douglass和Roberts曾提出其数学模型并研究了半离散化方法 [2-3】.袁益让对
此模型研究了特征一有限元方法 8【】和差分法 9【1.
本人对可压缩可混溶驱动问题的模型曾研究了共轭梯度迭代解与原问题真解的最优
阶日 模误差估计 [5J_其中饱和度方程的弥散项为一V·(D(x)Vc),而本:讨论的是D(u)
情况下的 模误差估计.就 模而言,对此模型目前尚未有人讨论过.从本文可看到,
由于饱和度方程中含有6(c)鬻这一项,使得 模误差估计不得不损失一个h因子,为
了得到 0模误差估计还必须利用到 日 模估计这一技巧.
下面给出此非线性抛物型方程组的初边值问题
a()d蟹)窑一 (。c())=d(co瓦p+Vu=‘口,z∈, t∈ o(,,
b()塞+6c() Vc—V.D(u()V e(-c)q,z (1.1)
文档评论(0)