Local radial basis function collocation method for linear thermoelasticity in two dimensions

Abstract

Purpose

The purpose of this paper is to upgrade our previous developments of Local Radial Basis Function Collocation Method (LRBFCM) for heat transfer, fluid flow and electromagnetic problems to thermoelastic problems and to study its numerical performance with the aim to build a multiphysics meshless computing environment based on LRBFCM.

Design/methodology/approach

Linear thermoelastic problems for homogenous isotropic body in two dimensions are considered. The stationary stress equilibrium equation is written in terms of deformation field. The domain and boundary can be discretized with arbitrary positioned nodes where the solution is sought. Each of the nodes has its influence domain, encompassing at least six neighboring nodes. The unknown displacement field is collocated on local influence domain nodes with shape functions that consist of a linear combination of multiquadric radial basis functions and monomials. The boundary conditions are analytically satisfied on the influence domains which contain boundary points. The action of the stationary stress equilibrium equation on the constructed interpolation results in a sparse system of linear equations for solution of the displacement field.

Findings

The performance of the method is demonstrated on three numerical examples: bending of a square, thermal expansion of a square and thermal expansion of a thick cylinder. Error is observed to be composed of two contributions, one proportional to a power of internodal spacing and the other to a power of the shape parameter. The latter term is the reason for the observed accuracy saturation, while the former term describes the order of convergence. The explanation of the observed error is given for the smallest number of collocation points (six) used in local domain of influence. The observed error behavior is explained by considering the Taylor series expansion of the interpolant. The method can achieve high accuracy and performs well for the examples considered.

Research limitations/implications

The method can at the present cope with linear thermoelasticity. Other, more complicated material behavior (visco-plasticity for example), will be tackled in one of our future publications.

Originality/value

LRBFCM has been developed for thermoelasticity and its error behavior studied. A robust way of controlling the error was devised from consideration of the condition number. The performance of the method has been demonstrated for a large number of the nodes and on uniform and non-uniform node arrangements.

Citation

Mavrič, B. and Šarler, B. (2015), "Local radial basis function collocation method for linear thermoelasticity in two dimensions", International Journal of Numerical Methods for Heat & Fluid Flow, Vol. 25 No. 6, pp. 1488-1510. https://doi.org/10.1108/HFF-11-2014-0359