DISCONTINUOUS GALERKIN FINITE ELEMENT METHODS FOR NUMERICAL SIMULATIONS OF THERMOELASTICITY

摘要

In this article, the discontinuous Galerkin (DG) finite element methods are extended to numerical simulations of uncoupling thermoelasticity. A thermoelastic formulation of interior penalty DG (IP-DG) method is presented and aspects of the numerical implementation are discussed in matrix form. The feasibility of the method for general thermoelastic simulations is validated through some typical test cases. Investigations are performed in which the method is applied to problems with thermal stress discontinuities induced by non-homogeneous inclusions, illustrating that reasonably ranged magnitude of IP term can well deal with the problems, and IP parameter increasing can effectively suppress the non-physical oscillations around discontinuities although its excess may affect solution accuracy in smooth region to some extent.