Finite Element Approximation of the Heat Equation under Axisymmetry Assumption

Authors: Raphael Zanella

Abstract:

This works deals with the finite element approximation of axisymmetric problems. The weak formulation of the heat equation under axisymmetry assumption is established for continuous finite elements. The weak formulation is implemented in a C++ solver with implicit time marching. The code is verified by space and time convergence tests using a manufactured solution. An example problem is solved with an axisymmetric formulation and with a 3D formulation. Both formulations lead to the same result but the code based on the axisymmetric formulation is mush faster due to the lower number of degrees of freedom. This confirms the correctness of our approach and the interest of using an axisymmetric formulation when it is possible.

Keywords: Axisymmetric problem, continuous finite elements, heat equation, weak formulation.

References:

[1] Z. P. Bazant, J.-C. Chern and W. Thonguthai, Finite Element Program for Moisture and Heat Transfer in Heated Concrete, Nuclear Engineering and Design 68, pp. 61-70, 1981.
[2] C. Chaboudez, S. Clain, R. Glardon, D. Mari, J. Rappaz and M. Swierkosz, Numerical Modeling in Induction Heating for Axisymmetric Geometries, IEEE Transactions on Magnetics 33(1), pp. 739-745, 1997.
[3] S. V. Litvinov, L. I. Trush and A. A. Avakov, Some Features of Temperature Field Definition in Axisymmetric Problems, 2017 International Conference on Industrial Engineering, Applications and Manufacturing.
[4] V. A. Dobrev, T. E. Ellis, T. V. Kolev and R. N. Rieben, High-order curvilinear finite elements for axisymmetric Lagrangian hydrodynamics, Computers & Fluids 83, pp. 58-69, 2013.
[5] MFEM: Modular Finite Element Methods library, https://mfem.org.
[6] A. Ern and J.-L. Guermond, Theory and Practice of Finite Elements, 1st ed., Springer, New York, 2004.