Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 72, 2017 - Issue 5
151
Views
6
CrossRef citations to date
0
Altmetric
Original Articles

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

&
Pages 377-391 | Received 31 Jul 2017, Accepted 26 Oct 2017, Published online: 27 Nov 2017
 

ABSTRACT

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.

Nomenclature

aj=

nonconstant coefficients in MLS

k=

number of iterations in preconditioner call

m=

number of terms in basis

M=

size of iteration cycle for SRJ method

n=

grid size of iteration cycle

N=

grid size

ns=

number of nodes in support domain

p=

level of SRJ scheme

=

specific heat flux

Qg=

heat generation W/m3

Q=

heat flux at Neumann boundary condition W/m2

rq=

radius of quadrature domain, m

rs=

radius of support domain, m

rw=

radius of weight domain, m

Th=

moving least square approximant

=

specified temperature at Dirichlet boundary

v=

test function for MLPG method

w=

weight function used in MLS approximation

z=

index in SRJ iteration cycle

Greek symbols=
αq=

dimensionless parameter of quadrature domain

αs=

dimensionless parameter of support domain

Γ1=

global domain boundary for Dirichlet boundary condition

Γ2=

global domain boundary for Neumann boundary condition

Γ=

global boundary

κ=

thermal conductivity, W/m°C

ϕI=

MLS shape function

Ω=

two-dimensional domain

ΩQ=

local domain

Ω=

boundary of local domain

ω=

relaxation factor for SOR

Superscript=
h=

approximated variable

Subscripts=
g=

heat generation

i,I j,k,l=

indices

q=

quadrature domain

Q=

local domain

S=

support domain

w=

weight function

Nomenclature

aj=

nonconstant coefficients in MLS

k=

number of iterations in preconditioner call

m=

number of terms in basis

M=

size of iteration cycle for SRJ method

n=

grid size of iteration cycle

N=

grid size

ns=

number of nodes in support domain

p=

level of SRJ scheme

=

specific heat flux

Qg=

heat generation W/m3

Q=

heat flux at Neumann boundary condition W/m2

rq=

radius of quadrature domain, m

rs=

radius of support domain, m

rw=

radius of weight domain, m

Th=

moving least square approximant

=

specified temperature at Dirichlet boundary

v=

test function for MLPG method

w=

weight function used in MLS approximation

z=

index in SRJ iteration cycle

Greek symbols=
αq=

dimensionless parameter of quadrature domain

αs=

dimensionless parameter of support domain

Γ1=

global domain boundary for Dirichlet boundary condition

Γ2=

global domain boundary for Neumann boundary condition

Γ=

global boundary

κ=

thermal conductivity, W/m°C

ϕI=

MLS shape function

Ω=

two-dimensional domain

ΩQ=

local domain

Ω=

boundary of local domain

ω=

relaxation factor for SOR

Superscript=
h=

approximated variable

Subscripts=
g=

heat generation

i,I j,k,l=

indices

q=

quadrature domain

Q=

local domain

S=

support domain

w=

weight function

Acknowledgments

The authors gratefully acknowledge the funding to Rituraj Singh from IIT-Roorkee through MHRD research fellowship and computational resources provided by CFD Laboratory, MIED, IIT-Roorkee.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.