TY - JOUR

T1 - A mass-redistributed finite element method (MR-FEM) for acoustic problems using triangular mesh

AU - He, Z. C.

AU - Li, Eric

AU - Liu, G. R.

AU - Li, G. Y.

AU - Cheng, A. G.

PY - 2016/10/15

Y1 - 2016/10/15

N2 - The accuracy of numerical results using standard finite element method (FEM) in acoustic problems will deteriorate with increasing frequency due to the “dispersion error”. Such dispersion error depends on the balance between the “stiffness” and “mass” of discretization equation systems. This paper reports an improved finite element method (FEM) for solving acoustic problems by re-distributing the mass in the mass matrix to “tune” the balance, aiming to minimize the dispersion errors. This is done by shifting the integration point locations when computing the entries of the mass matrix, while ensuring the mass conservation. The new method is verified through the detailed numerical error analysis, and a strategy is also proposed for the best mass redistribution in terms of minimizing dispersion error. The relative dispersion error of present mass-redistributed finite element method (MR-FEM) is found to be much smaller than the FEM solution, in both theoretical prediction and numerical examination. The present MR-FEM works well by using the linear triangular elements that can be generated automatically, which enables automation in computation and saving computational cost in mesh generation. Numerical examples demonstrate the advantages of MR-FEM, in comparison with the standard FEM using the same triangular meshes and quadrilateral meshes.

AB - The accuracy of numerical results using standard finite element method (FEM) in acoustic problems will deteriorate with increasing frequency due to the “dispersion error”. Such dispersion error depends on the balance between the “stiffness” and “mass” of discretization equation systems. This paper reports an improved finite element method (FEM) for solving acoustic problems by re-distributing the mass in the mass matrix to “tune” the balance, aiming to minimize the dispersion errors. This is done by shifting the integration point locations when computing the entries of the mass matrix, while ensuring the mass conservation. The new method is verified through the detailed numerical error analysis, and a strategy is also proposed for the best mass redistribution in terms of minimizing dispersion error. The relative dispersion error of present mass-redistributed finite element method (MR-FEM) is found to be much smaller than the FEM solution, in both theoretical prediction and numerical examination. The present MR-FEM works well by using the linear triangular elements that can be generated automatically, which enables automation in computation and saving computational cost in mesh generation. Numerical examples demonstrate the advantages of MR-FEM, in comparison with the standard FEM using the same triangular meshes and quadrilateral meshes.

UR - http://www.scopus.com/inward/record.url?scp=84980484872&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2016.07.025

DO - 10.1016/j.jcp.2016.07.025

M3 - Article

AN - SCOPUS:84980484872

VL - 323

SP - 149

EP - 170

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -