A diffusion-split method to deal with thermal shocks using standard linear tetrahedral finite elements

Victor D. Fachinotti (1, 2) , Michel Bellet (1)

1 CEMEF - Centre de Mise en Forme des Matériaux
2 CONICET - Consejo Nacional de Investigaciones Científicas y Técnicas [Buenos Aires]

The thermal analysis using linear standard tetrahedral finite elements may be affected by spurious local extrema in the regions affected by thermal shocks, in such a severe way to directly discourage the use of these elements. The present work proposes a slight modification to the discrete heat equation in order to obtain a system matrix in M-matrix form, which assures an oscillation-free solution. The performance of this method is evaluated by means of test case with analytical solution, as well as an industrial application, for which a well-behaved numerical solution is available.

Mots clés: Diffusion; Tetrahedral Finite Elements

Domaines: Sciences de l'ingénieur [physics]; Matériaux

Victor D. Fachinotti, Michel Bellet. A diffusion-split method to deal with thermal shocks using standard linear tetrahedral finite elements. Materials Processing and Design: Modeling, Simulation and Applications - NUMIFORM 2004 - Proceedings of the 8th International Conference on Numerical Methods in Industrial Forming Processes, Jun 2004, Colombus, United States. ISBN : 0-7354-0188-8 - p. 2258-2263, ⟨10.1063/1.1766871⟩. ⟨hal-00531257⟩