计算物理 ›› 2000, Vol. 17 ›› Issue (4): 372-380.

• 论文 • 上一篇    下一篇

非匹配网格上广义Stokes问题的区域分裂解法及事后误差估算

周春华   

  1. 南京航空航天大学六系, 江苏 南京 210016
  • 收稿日期:1998-12-03 修回日期:1999-08-10 出版日期:2000-07-25 发布日期:2000-07-25
  • 作者简介:周春华(1965~),男,江苏东台,讲师,博士,主要从事科学与工程计算及数值分析的研究.
  • 基金资助:
    留学回国人员科研启动基金资助项目

DOMAIN DECOMPOSITION FOR GENERALIZED STOKES PROBLEM ON NON-COINCIDENT MESHES AND A POSTERIORI ERROR ESTIMATION

ZHOU Chun Hua   

  1. Department 6, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P R China
  • Received:1998-12-03 Revised:1999-08-10 Online:2000-07-25 Published:2000-07-25

摘要: 发展了一种广义Stokes问题的无覆盖区域分裂解法。子域交界面上的约束条件是通过引入一Lagrange乘子而得到弱满足的,在有限元离散子域的交界处网格可以是非匹配的。应用Petrov Galerkin方法解每个子域上的广义Stokes问题,而交界面上的Lagrange乘子则通过共轭梯度法迭代求解,各变量均由线性函数离散。对上述区域分裂解法,还构造了基于求解当地问题的误差事后估算方法。各变量的当地误差估算器定义在二阶非连续鼓包(bump)函数的空间中。最后给出了基于事后误差估算值的自适应网格上的数值结果。

关键词: 区域分裂, 有限元, Stokes问题, 自适应网格, 误差估算

Abstract: A method of domain decomposition for generalized Stokes problem is developed. The condition of compatibility on the interfaces between each subdomain is satisfied weakly by introducing a Lagrange multiplier technique. In the computational domain, the finite element meshes can be non-coincident at the interfaces. The Petrov-Galerkin formulation is applied to solve the generalized Stokes problem in each subdomain, and the Lagrange multiplier problem on the interfaces is solved using an algorithm of conjugated gradient. Velocity, density (or pressure) and Lagrange multiplier are approached in the spaces of continue piecewise linear functions. Associated with the above domain decomposition method, a local a posteriori error estimation has been constructed. The localization of the error estimation is based on solving local problems. The local error estimators respectively for velocity, density and Lagrange multiplier, are defined in the space of discontinue quadratic bump functions. At the end, some numerical results on the adaptive meshes are also given based on local a posteriori error estimation.

Key words: domain decomposition, finite element, Stokes problem, self-adaptive mesh, error estimation

中图分类号: