In the case of large scale simulations involving complex domain, the task of mesh generation is quite a challenging and tedious task both in terms of time and human effort. In recent years, significant research efforts have been invested in the development of mesh-free methods to solve engineering and scientific problems without mesh generation. These methods also overcome the shortcomings of the mesh-based methods such as element distortion in case of problems involving large deformations and crack propagation. Of these methods, the meshless local Petrov Galerkin (MLPG) method is amongst the noticeable ones.  MLPG method is a promising meshfree method for continuum problems in complex domains, especially for large deformation, moving boundary and phase change problems.

However, most of the applications have been limited to small scale problems which are usually solved using direct solvers. Direct solvers such as the Gauss elimination method or LU decomposition method are of the complexity of O(N^3^), where N is the order of the stiffness matrix. There are large scale problems, which have to be modeled and simulated with high accuracy and results in millions of discretization points. The use of direct solvers is not practical for such large scale problems. Very least number of attempts have been made to solve large scale complex geometry problems using MLPG approach which needs many millions (or even billions) of computing nodes and would require fast iterative solvers. The MLPG method leads to the asymmetric system of algebraic equations. Bi-Conjugate Gradient Stabilised (BiCGSTAB) and Generalized Minimum Residual (GMRES) solvers are the most popular Krylov subspace solvers for an asymmetric linear system of equations. The convergence rate of these methods depends on the preconditioner used. There are basic iterative solvers such as Jacobi and Gauss-Seidel, but these are very slow in convergence. However, these basic methods can still be used as a preconditioner for Krylov subspace solvers.    

 

In the present work, restarted version of the GMRES method, with and without preconditioner, is applied in connection with the Interpolating MLPG method. Two dimensional and three-dimensional linear steady-state heat conduction problems are taken as model test problems. The error has been calculated based on the exact solution of these problems. Its performance is compared with the preconditioned BiCGSTAB method based on computation time and convergence behavior. Jacobi and successive over-relaxation (k) methods are used as a preconditioner in both the solvers. This work is an attempt to develop a robust and fast MLPG solver.


GMRES Solver for MLPG Method Applied to Heat Conduction

Division / Sponsors: 
Aerospace, Bioengineering, Applied Mechanics, Design Engineering Division, Advanced Energy Systems, Manufacturing Engineering Division, Noise Control and Acoustics, Nondestructive Evaluation, Diagnosis, and Prognosis (NDPD), Nuclear Engineering Division, Process Industries Division, Fluids Engineering, Heat Transfer, Materials, Microelectromechanical Systems, Safety Engineering and Risk Analysis, Computers & Information in Engineering


Paper Proceedings: https://asmedigitalcollection.asme.org/IMECE/IMECE2020/volume/84478

IMECE2020 - Virtual Conference

Supporting the need among researchers worldwide for additional online and multimedia content, ASME now offers the ASME Conference Video Collection. This platform will collect and provide long-term access to conference presentation videos, encompassing all technology areas within ASME's technical conference and events program. About 30,000-40,000 engineers across all mechanical engineering and related disciplines participate in ASME conferences each year. This includes authors of approximately 6,000 to 8,000 technical papers published by ASME, the majority of whom also submit video presentations for our virtual conferences. The video collection will be available as an add-on to existing product packages or as a separate subscription.