A new boundary meshfree method for calculating the multi-domain constant coefficient heat conduction with a heat source problem

A new high-precision boundary meshfree method, namely virtual boundary meshfree Galerkin method (VBMGM), for calculating the multi-domain constant coefficient heat conduction with a heat source problem is given. In the paper, the radial basis function interpolation is used to solve the virtual source function of virtual boundary and the heat source within each subdomain. Simultaneously, the equation of VBMGM for multi-domain constant coefficient heat conduction with a heat source problem is obtained by the Galerkin method. Therefore, the proposed method has common advantages of the boundary element method, meshfree method, and Galerkin method. Coefficient matrix of this specific expression is symmetrical and the specific expression of VBMGM for the multi-domain constant coefficient heat conduction with a heat source problem is given. Two numerical examples are given. The numerical results are also compared with other numerical methods. The accuracy and feasibility of the method for the multi-domain constant coefficient heat conduction with a heat source problem are proved.     

A new boundary-type mesh-free method for solving the multi-domain heat conduction problem

ABSTRACT

In this article, a new high-precision boundary-type mesh-free method, called the virtual boundary mesh-free Galerkin method (VBMGM), is given for analyzing the heat conduction problem consisting of multiple media. A control equation of the proposed method is established by the the Galerkin method of the weighted residual methods. The virtual source function is approached through the radial basis function interpolation in the mesh-free method. Thus, the main feature of the proposed technique has the advantages of the mesh-free method, the boundary-element method, and the Galerkin method, a symmetrical coefficient matrix. The numerical meanings of the weighted values in the VBMGM equation, such as the partial derivatives the temperature, the heat flux, the temperature connection condition, and the heat flux connection condition, are clear. Numerical examples of a thick-walled cylinder consisting of two media and a two-layered rectangular plate are provided. The feasibility and the accuracy of the proposed technique are proved.     

Virtual boundary element method in conjunction with conjugate gradient algorithm for three-dimensional inverse heat conduction problems

In this article, the virtual boundary element method (VBEM) in conjunction with conjugate gradient algorithm (CGA) is employed to treat three-dimensional inverse problems of steady-state heat conduction. On the one hand, the VBEM may face numerical instability if a virtual boundary is improperly selected. The numerical accuracy is very sensitive to the choice of the virtual boundary. The condition number of the system matrix is high for the larger distance between the physical boundary and the fictitious boundary. On the contrary, it is difficult to remove the source singularity. On the other hand, the VBEM will encounter ill-conditioned problem when this method is used to analyze inverse problems. This study combines the VBEM and the CGA to model three-dimensional heat conduction inverse problem. The introduction of the CGA effectively overcomes the above shortcomings, and makes the location of the virtual boundary more free. Furthermore, the CGA, as a regularization method, successfully solves the ill-conditioned equation of three-dimensional heat conduction inverse problem. Numerical examples demonstrate the validity and accuracy of the proposed method.     

On preconditioned BiCGSTAB solver for MLPG method applied to heat conduction in complex geometry

Meshless local Petrov–Galerkin (MLPG) method is a promising meshfree method for continuum problems in complex domains, especially for large deformation, moving boundary and phase change problems. For large-scale problems, iterative methods for solving the discretized equations are more suitable than direct methods. Krylov subspace solvers of conjugate gradient type are the most preferred iterative solvers. The convergence rate of these methods depends on preconditioner used. Recently, proposed schedule relaxation Jacobi (SRJ) method can be used as a stand-alone solver and as a preconditioner. In the present work, the SRJ method is tested as a stand-alone solver and as a preconditioner for BiCGSTAB solver in the MLPG method, and its performance has been compared with successive overrelaxation (k) preconditioner. Two-dimensional linear steady-state heat conduction in complex shape geometry has been used as the model test problem.