A Discontinuous Galerkin Finite Element Method for Heat Conduction Problems with Local High Gradient and Thermal Contact Resistance

A discontinuous Galerkin (DG) finite element method for the heat conduction problems with local high gradient and thermal contact resistance is presented. The DG formulation Is constructed by employing the stabilization term and the Bassi-Rebay numerical flux term. The stabilization term is defined by a penalization of the temperature jump at the interface. By eliminating the penalization term of the temperature jump in the region of local high gradient and imperfect contact interfaces, the present DG method is applied to solve problems involving local high gradient and thermal contact resistance where the numerical flux is obtained from the definition of the thermal contact resistance. This treatment leads to a novel approach to capture the peak value of the heat flux in local high gradient field and the temperature jump at the imperfect contact interface. Moreover, all iterative procedure and a relaxation technique are also adopted herein to simulate the nonlinear thermoelastic Coupling between thermal contact resistance, temperature and stress field, which can easily avoid numerical instability and gets reasonable results. Several numerical examples are given to demonstrate the accuracy and the reliability of the present DG finite element method. From Our investigations, we find that the DG method is an attractive and competitive approach for solving heat conduction problems with local high gradient and thermal contact resistance, and also call easily avoid numerically instability in dealing with nonlinear Coupling problems caused by thermal contact resistance.