Discontinuous Galerkin method in numerical simulation of two-dimensional thermoelasticity problem with single stabilization parameter. (August 2018)

• Record Type: Journal Article
• Authors: Jaśkowiec, Jan; Pluciński, Piotr
• Is Part Of: Advances in engineering software. Volume 122(2018)
• Page Start: 62
• Page End: 80
• Publication Date: 2018-08
• Subjects: Discontinuous Galerkin -- Thermoelasticity -- Polygonal finite elements -- Stabilization parameter; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2018.04.015
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: Coupled 2D non-stationary thermomechanical problem is solved by discontinuous Galerkin method with finite difference rules (DGFD). Continuities of the displacements and temperature fields are enforced with the help of numerical fluxes using single stabilization parameter. Boundary conditions are enforced with the help of high-order finite difference rules. The Chebyshev basis functions are used for approximations. Zienkiewicz–Zhu smoothing procedure is adopted for DGFD method. Untypical polygonal finite elements are used, for example fish-like cells. Very high-order finite elements are used in which p > 10 . Very high-order finite elements can be easily combined with low-order finite elements. Abstract: The aim of the paper is the development of discontinuous Galerkin with finite difference rules (DGFD) to a two-dimensional stationary and non-stationary thermoelasticity problem. Displacement and temperature fields are approximated on the same mesh frame but with various approximation orders, which are set independently for each of the fields. Because the DGFD method does not use nodes, special attention needs to be paid to applying boundary conditions. Various types of thermal and mechanical boundary conditions are considered. In the presented approach only one stabilization parameter for the coupled problem needs to be evaluated in the DGFD method. The same parameter used in thermal and in mechanical part. The considered domain is discretized by a polygonal mesh in which the polygonal elements may have arbitrary shapes, such as e.g. a fish shape, as well as typical rectangular shapes. The orthogonality of Chebyshev basis functions may be utilized for rectangular elements. Very high-order approximate solution can be obtained in such case. In the coupled problem, the same element may be high-order for displacement field while low-order to approximate temperature. The argument contained in the paper is illustrated with few examples.