Monolithic discretization of linear thermoelasticity problems via adaptive multimesh hp-FEM

  • Authors:

  • P. Solin;J. Cerveny;L. Dubcova;D. Andrs

  • Affiliations:

  • University of Nevada, Reno, NV 89557, USA and Institute of Thermomechanics, Academy of Sciences of Czech Republic, Dolejskova 5, Prague, Czech Republic;Institute of Thermomechanics, Academy of Sciences of Czech Republic, Dolejskova 5, Prague, Czech Republic;Institute of Thermomechanics, Academy of Sciences of Czech Republic, Dolejskova 5, Prague, Czech Republic;University of Nevada, Reno, NV 89557, USA and Institute of Thermomechanics, Academy of Sciences of Czech Republic, Dolejskova 5, Prague, Czech Republic

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2010

Quantified Score

Hi-index 7.31

Visualization

Abstract

In linear thermoelasticity models, the temperature T and the displacement components u"1,u"2 exhibit large qualitative differences: while T typically is very smooth everywhere in the domain, the displacements u"1,u"2 have singular gradients (stresses) at re-entrant corners and edges. The mesh must be extremely fine in these areas so that stress intensity factors are resolved sufficiently. One of the best available methods for this task is the exponentially-convergent hp-FEM. Note, however, that standard adaptive hp-FEM approximates all three fields u"1,u"2 and T on the same mesh, and thus it treats T as if it were singular at re-entrant corners as well. Therefore, a large number of degrees of freedom of temperature are wasted. This motivates us to approximate the fields u"1,u"2 and T on individual hp-meshes equipped with mutually independent hp-adaptivity mechanisms. In this paper we describe mathematical and algorithmic aspects of the novel adaptive multimesh hp-FEM, and demonstrate numerically that it performs better than the standard adaptive h-FEM and hp-FEM.