Stability analysis of generalized mass formulation in dynamic heat transfer

In this article, the generalized mass formulation is developed in an explicit analysis of transient transport problems. It has been well known that the time step is typically smaller in explicit analysis than in implicit analysis when the same size mesh is used. Further, the over-stiffness of conventional finite-element model may result in poor accuracy with linear triangular or tetrahedral elements. In order to improve the computational efficiency and numerical accuracy, this article proposes a generalized mass formulation by matching the mass matrix to the smoothed stiffness matrix using linear triangular elements in 2-D problems. The proposed mass matrix can be obtained by simply shifting the integration points from the conventional locations. Without loss of generality, several 2-D examples, including conduction, convection, and radiation heat transfer problems, are presented to demonstrate that the generalized mass formulation allows a larger time step in explicit analysis compared with the lumped and consistent mass matrices. In addition, it is found that the maximum allowable time step is proportional to the softened effect of the discretized model in an explicit analysis.