Variational multiscale element free Galerkin (VMEFG) and local discontinuous Galerkin (LDG) methods for solving two-dimensional Brusselator reaction-diffusion system with and without cross-diffusion

详细信息    查看全文

• 作者：Mehdi Dehghan ; ^{mdehghan@aut.ac.ir" class="auth_mail" title="E-mail the corresponding author} ; ^{mdehghan.aut@gmail.com" class="auth_mail" title="E-mail the corresponding author} ; Mostafa Abbaszadeh ^{m.abbaszadeh@aut.ac.ir" class="auth_mail" title="E-mail the corresponding author}
• 关键词：65M06 ; 65M20 ; 65M60
• 刊名：Computer Methods in Applied Mechanics and Engineering
• 出版年：2016
• 出版时间：1 March 2016
• 年：2016
• 卷：300
• 期：Complete
• 页码：770-797
• 全文大小：5217 K

文摘

The finite element method (FEM) is one of the basic methods for solving deterministic and stochastic partial differential equations. This method is proposed in the 19 decade and after years several modifications for this well-known technique such as mortal FEM, discontinuous Galerkin FEM, extended FEM, least squares FEM, spectral FEM, mixed FEM, immersed FEM, adaptive FEM, etc. have been proposed. Also, for improving the accuracy and for considering complex domain, shape functions of the moving least squares approximation have been replaced with the traditional shape function of the finite element technique. The name of this improved method is element free Galerkin (EFG) method which is classified in the category of meshless methods. The EFG method has been applied for solving many problems in engineering. In the current paper, we select two numerical procedures that are extracted from FEM. The discontinuous Galerkin approach is a useful technique for solving problems that their exact solutions have discontinuity. It should be noted that many enriched ideas have been explained for improving accuracy of EFG method. In the current paper, the local discontinuous Galerkin and variational multiscale EFG methods are applied for obtaining numerical solution of two-dimensional Brusselator system with and without cross-diffusion. The Brusselator system is a theoretical model for a type of autocatalytic reaction. Up to best of authors’ knowledge, the accuracy of the obtained numerical solutions using local discontinuous Galerkin method is very related to the used time-discrete scheme. Therefore, fourth-order exponential time differencing method has been employed for discretizing the time variable. Also, to achieve a full discretization scheme, the local discontinuous Galerkin finite element method is used. Moreover, several test problems are given that show the acceptable accuracy and efficiency of the proposed schemes.