期刊
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING
卷 67, 期 1, 页码 89-107出版社
WILEY
DOI: 10.1002/nme.1628
关键词
mixed finite element; mass lumping; finite volume; nonconforming finite element
Diffusion-type partial differential equation is a common mathematical model in physics. Solved by mixed finite elements, it leads to a system matrix which is not always an M-matrix. Therefore, the numerical solution may exhibit unphysical results due to oscillations. The criterion necessary to obtain an M-matrix is discussed in details for triangular, rectangular and tetrahedral elements. It is shown that the system matrix is never an M-matrix for rectangular elements and can be an M-matrix for trianaular an tetrahedral elements if criteria on the element's shape and on the time step length are fulfilled. A new mass lumping scheme is developed which leads to a less restrictive criterion: the discretization must be weakly acute (all angles less than pi/2) and there is no constraint on the time step length. The lumped formulation of mixed hybrid finite element can be applied not only to triangular meshes but also to more general shape elements in two and three dimensions. Numerical experiments show that, compared to the standard mixed hybrid formulation, the lumping scheme avoids (or strongly reduce) oscillations and does not create additional numerical errors. Copyright (c) 2006 John Wiley & Sons, Ltd.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据