Isogeometric boundary element method for axisymmetric steady-state heat transfer

Abstract

Isogeometric analysis (IGA) utilizes the Non-Uniform Rational B-Splines (NURBS) basis functions, which are commonly employed in Computer Aided Design (CAD), to both construct the problem’s geometry and approximate the field variables. In axisymmetric scenarios, the 3D problem can be simplified to 2D, requiring only the discretization of the 1D curve boundary for the boundary element method (BEM). This study proposes an isogeometric boundary element method (IGABEM) to simulate steady-state heat transfer in axisymmetric domains. Both the precise geometry of the boundary and the physical variables are approximated with the same NURBS basis functions, yielding higher-order continuity, superior accuracy per degree of freedom (DOF), and continuous nodal gradients with fewer DOFs than the traditional BEM. Additionally, a zoning method is introduced to handle heterogeneities. Five cases including a solid cylinder, a revolution hyperboloid with a coaxial cylindrical tube, an ellipsoid with a spherical cavity, a bimaterial ellipsoid and a complex geometry with a toroidal tube validate the accuracy, convergence, computational efficiency of IGABEM, and the applicability in addressing multiply-connected domains. The method outperforms Finite Element Method (FEM) and conventional BEM in accurately handling steady-state axisymmetric heat transfer issues under Dirichlet, Neumann, and Robin boundary conditions.