计算数学
       首页 |  期刊介绍 |  编委会 |  投稿指南 |  期刊订阅 |  下载中心 |  留言板 |  联系我们 |  重点论文 |  在线办公 | 
计算数学  2019, Vol. 41 Issue (3): 259-265    DOI:
论文 最新目录 | 下期目录 | 过刊浏览 | 高级检索 Previous Articles  |  Next Articles  
求解定常不可压Stokes方程的两层罚函数方法
李世顺1, 祁粉粉1, 邵新平2
1. 河南理工大学数学与信息科学学院, 焦作 454003;
2. 杭州电子科技大学理学院, 杭州 310027
TWO-LEVEL PENALTY METHOD FOR THE STEADY INCOMPRESSIBLE STOKES EQUATIONS
Li Shishun1, Qi FenFen1, Shao Xinping2
1. School of Mathematics & Information Science, Henan Polytechnic University, Jiaozuo 454003, China;
2. School of Science & Hangzhou Dianzi University, Hangzhou 310027, China
 全文: PDF (322 KB)   HTML (1 KB)   输出: BibTeX | EndNote (RIS)      背景资料
摘要 借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程(离散后的线性方程组具有相同的对称正定系数矩阵).收敛性分析表明粗网格空间相对于细网格空间可以选择很小,并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果,并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.
服务
把本文推荐给朋友
加入我的书架
加入引用管理器
E-mail Alert
RSS
作者相关文章
关键词两层罚函数方法   定常不可压Stokes方程   罚参数   有限元方法     
Abstract: In this paper, we present a two-level penalty method for the steady incompressible Stokes equations by employing two finite element spaces. This method involves solving one small Stokes equation on the coarse space and two penalty equations on the fine space (the linear systems with same symetric and positive coefficient matrices). The convergence shows that the coarse space can be chosen very small. Moreover, the penalty parameter is only dependent on the coarse mesh size and the regularity of the problem. Therefore, the resulting solution still achieves asymptotically optimal accuracy when the penalty parameter is chosen "not very small". The numerical results confirm the convergence analysis, and the numerical comparison also shows that this method is efficient for solving the steady incompressible Stokes equations.
Key wordstwo-level penalty method   steady incompressible Stokes equations   penalty parameter   finite element method   
收稿日期: 2017-09-02; 出版日期: 2019-08-21
基金资助:

国家自然科学基金项目(No.11401177,11701133),浙江省教育厅科研项目(No.Y201533698).

通讯作者: 李世顺,Email:lss6@sina.com     E-mail: lss6@sina.com
引用本文:   
. 求解定常不可压Stokes方程的两层罚函数方法[J]. 计算数学, 2019, 41(3): 259-265.
. TWO-LEVEL PENALTY METHOD FOR THE STEADY INCOMPRESSIBLE STOKES EQUATIONS[J]. Mathematica Numerica Sinica, 2019, 41(3): 259-265.
 
[1] Brenner S C, Scott L R. The Mathematical Theory of Finite Element Methods[M]. New York:Springer, 2002.
[2] Brezzi F, Fortin M. Mixed and Hybrid Finite Element Methods[M]. New York:Springer, 1991.
[3] Carey G F, Krishnan R. Penalty approximation of Stokes flow[J]. Comput. Methods Appl. Mech. Engrg., 1992, 35:169-206.
[4] Cheng X L, Shaikh A W. Analysis of the iterative penalty method for the Stokes equations[J]. Appl. Math. Lett., 2006, 19:1024-1028.
[5] Dai X X, Cheng X L. The iterative penalty method for Stokes equations using Q1-P0 element[J]. Appl. Math. Comput., 2008, 201:805-810.
[6] Elman H C, David J S, Andrew J W. Finite Elements and Fast Iterative Solvers:With Applications in Incompressible Fluid Dynamics[M]. New York:Oxford University Press, 2005.
[7] He Y. Two-level method based on finite element and Crank-Nicolson extrapolation for the timedependent Navier-Stokes equations[J]. SIAM J. Numer. Anal., 2003, 41(4):1263-1285.
[8] Huang P, He Y, Feng X. Convergence and stability of two-level penalty mixed finite element method for stationary Navier-Stokes equations[J]. Front. Math. China, 2013, 8:837-854.
[9] Kheshgi H, Luskin M. Analysis of the finite element variable penalty method for Stokes equations[J]. Math. Comp., 1985, 45(172):347-363.
[10] Layton W, Tobiska L. A two-level method with backtracking for the Navier-Stokes equations[J]. SIAM J. Numer. Anal., 1998, 35:2035-2054.
[11] Layton W, Ye X. Two-level discretizations of the stream function form of the navier-stokes equations[J]. Numer. Funct. Anal. Optim., 1999, 20:909-916.
[12] Lin S Y, Chin Y S, Wu T M. A modified penalty method for Stokes equations and its applications to Navier-Stokes equations[J], SIAM J. Sci. Comput., 1995, 16(1):1-19.
[13] Oden J T, Kikuchi N, Song Y J. Penalty finite element methods for the analysis of Stokesian flows[J]. Comput. Methods Appl. Mech. Engrg., 1992, 31:297-329.
[14] Xu J. Two-grid discretization techinques for linear and nonlinear PDEs[J]. SIAM J. Numer. Anal., 1996, 33(5):1759-1777.
[1] 王俊俊, 李庆富, 石东洋. 非线性抛物方程混合有限元方法的高精度分析[J]. 计算数学, 2019, 41(2): 191-211.
[2] 邓维山, 徐进. 一种泊松-玻尔兹曼方程稳定算法的高效有限元并行实现[J]. 计算数学, 2018, 39(2): 91-110.
[3] 武海军. 高波数Helmholtz方程的有限元方法和连续内罚有限元方法[J]. 计算数学, 2018, 40(2): 191-213.
[4] 葛志昊, 吴慧丽. 体积约束的非局部扩散问题的后验误差分析[J]. 计算数学, 2018, 40(1): 107-116.
[5] 余涛, 张镭. 线性弹性问题的局部正交分解方法[J]. 计算数学, 2018, 39(1): 10-19.
[6] 杨建宏. 定常Navier-Stokes问题低次等阶稳定有限体积元算法研究[J]. 计算数学, 2017, 38(2): 91-104.
[7] 周宇, 李秋齐. 基于降基多尺度有限元的PGD方法及其在含参数椭圆方程中的应用[J]. 计算数学, 2017, 38(2): 105-122.
[8] 李宏, 杜春瑶, 赵智慧. 反应扩散方程的连续时空有限元方法[J]. 计算数学, 2017, 39(2): 167-178.
[9] 李晓翠, 杨小远, 张英晗. 一类随机非自伴波方程的半离散有限元近似[J]. 计算数学, 2017, 39(1): 42-58.
[10] 曹济伟, 葛志昊, 刘鸣放. Stokes方程基于多尺度函数的稳定化有限元方法[J]. 计算数学, 2017, 38(1): 68-80.
[11] 曹济伟. 求解二维时谐Maxwell方程的一种混合有限元新格式[J]. 计算数学, 2016, 38(4): 429-441.
[12] 赵智慧, 李宏, 罗振东. Sobolev方程的连续时空有限元方法[J]. 计算数学, 2016, 38(4): 341-353.
[13] 王军平, 叶秀, 张然. 弱有限元方法简论[J]. 计算数学, 2016, 38(3): 289-308.
[14] 郑权, 高玥, 秦凤. Helmholtz方程外边值问题的基于修正的DtN边界条件的有限元方法[J]. 计算数学, 2016, 38(2): 200-211.
[15] 石东洋, 史艳华, 王芬玲. 四阶抛物方程H1-Galerkin混合有限元方法的超逼近及最优误差估计[J]. 计算数学, 2014, 36(4): 363-380.

Copyright 2008 计算数学 版权所有
中国科学院数学与系统科学研究院 《计算数学》编辑部
北京2719信箱 (100190) Email: gxy@lsec.cc.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发
技术支持: 010-62662699 E-mail:support@magtech.com.cn
京ICP备05002806号-10