Local least–squares element differential method for solving heat conduction problems in composite structures

In this article, a completely new numerical method called the Local Least-Squares Element Differential Method (LSEDM), is proposed for solving general engineering problems governed by second order partial differential equations. The method is a type of strong-form finite element method. In this method, a set of differential formulations of the isoparametric elements with respect to global coordinates are employed to collocate the governing differential equations and Neumann boundary conditions of the considered problem to generate the system of equations for internal nodes and boundary nodes of the collocation element. For each outer boundary or element interface, one equation is generated using the Neumann boundary condition and thus a number of equations can be generated for each node associated with a number of element interfaces. The least-squares technique is used to cast these interface equations into one equation by optimizing the local physical variable at the least-squares formulation. Thus, the solution system has as many equations as the total number of nodes of the present heat conduction problem. The proposed LSEDM can ultimately guarantee the conservativeness of the heat flux across element surfaces and can effectively improve the solution stability of the element differential method in solving problems with hugely different material properties, which is a challenging issue in meshfree methods. Numerical examples on two- and three-dimensional heat conduction problems are given to demonstrate the stability and efficiency of the proposed method.