Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients

JOURNAL OF COMPUTATIONAL PHYSICS

卷:229

Article

Zeman, Jan¹  Vondrejc, Jaroslav¹  Novak, Jan²  Marek, Ivo³ 

[1] Czech Tech Univ, Fac Civil Engn, Dept Mech, Prague 16629 6, Czech Republic

[2] Czech Tech Univ, Fac Civil Engn, Ctr Integrated Design Adv Struct, Prague 16629 6, Czech Republic

[3] Czech Tech Univ, Fac Civil Engn, Dept Math, Prague 16629 6, Czech Republic

关键词: Numerical homogenization;    FFT-based solvers;    Trigonometric collocation method;    Conjugate gradient solvers;   

DOI  :  10.1016/j.jcp.2010.07.010

来源: Elsevier

【 摘 要 】

In this short note, we present a new technique to accelerate the convergence of a FFT-based solver for numerical homogenization of complex periodic media proposed by Moulinec and Suquet [1]. The approach proceeds from discretization of the governing integral equation by the trigonometric collocation method due to Vainikko [2], to give a linear system which can be efficiently solved by conjugate gradient methods. Computational experiments confirm robustness of the algorithm with respect to its internal parameters and demonstrate significant increase of the convergence rate for problems with high-contrast coefficients at a low overhead per iteration. (C) 2010 Elsevier Inc. All rights reserved.

【 授权许可】

Free