A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media

In this paper, we assess the performance of four iterative algorithms for solving non-symmetric rank-deficient linear systems arising in the FFT-based homogenization of heterogeneous materials defined by digital images. Our framework is based on the Fourier-Galerkin method with exact and approximate integrations that has recently been shown to generalize the Lippmann-Schwinger setting of the original work by Moulinec and Suquet from 1994. It follows from this variational format that the ensuing system of linear equations can be solved by general-purpose iterative algorithms for symmetric positive-definite systems, such as the Richardson, the Conjugate gradient, and the Chebyshev algorithms, that are compared here to the Eyre-Milton scheme - the most efficient specialized method currently available. Our numerical experiments, carried out for two-dimensional elliptic problems, reveal that the Conjugate gradient algorithm is the most efficient option, while the Eyre-Milton method performs comparably to the Chebyshev semi-iteration. The Richardson algorithm, equivalent to the still widely used original Moulinec-Suquet solver, exhibits the slowest convergence. Besides this, we hope that our study highlights the potential of the well-established techniques of numerical linear algebra to further increase the efficiency of FFT-based homogenization methods.

[1]  R. Varga,et al.  Chebyshev semi-iterative methods, successive overrelaxation iterative methods, and second order Richardson iterative methods , 1961 .

[2]  Matti Schneider,et al.  Use of composite voxels in FFT-based homogenization , 2015 .

[3]  E. Garboczi,et al.  Finite Element and Finite Difference Programs for Computing the Linear Electric and Elastic Properties of Digital Images of Random Materials | NIST , 1998 .

[4]  Luc Dormieux,et al.  Corrigendum to “Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites” [Comput. Methods Appl. Mech. Engrg. 217–220 (2012) 197–212] , 2015 .

[5]  Jan Novák,et al.  Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients , 2010, J. Comput. Phys..

[6]  Zdenek Strakos,et al.  Composite convergence bounds based on Chebyshev polynomials and finite precision conjugate gradient computations , 2014, Numerical Algorithms.

[7]  M. Schneider,et al.  Computational homogenization of elasticity on a staggered grid , 2016 .

[8]  Graeme W. Milton,et al.  A fast numerical scheme for computing the response of composites using grid refinement , 1999 .

[9]  Frederic Legoll,et al.  Periodic homogenization using the Lippmann--Schwinger formalism , 2014, 1411.0330.

[10]  Jaroslav Vondrejc,et al.  An FFT-based Galerkin method for homogenization of periodic media , 2013, Comput. Math. Appl..

[11]  Guy Bonnet,et al.  Effective properties of elastic periodic composite media with fibers , 2007 .

[12]  Noboru Kikuchi,et al.  Digital image-based modeling applied to the homogenization analysis of composite materials , 1997 .

[13]  C. Lanczos.  Chebyshev polynomials in the solution of large-scale linear systems , 1952, ACM '52.

[14]  L. Dormieux,et al.  Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites , 2012 .

[15]  Graeme W. Milton,et al.  An accelerated FFT algorithm for thermoelastic and non‐linear composites , 2008 .

[16]  D Stoyan,et al.  Reconstruction of random heterogeneous media , 2015, Journal of microscopy.

[17]  Bernard Sonon,et al.  An advanced approach for the generation of complex cellular material representative volume elements using distance fields and level sets , 2015 .

[18]  Jan Zeman,et al.  Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method , 2014, 1404.3614.

[19]  Hervé Moulinec,et al.  A computational scheme for linear and non‐linear composites with arbitrary phase contrast , 2001 .

[20]  V. G. Kouznetsova,et al.  Multi-scale computational homogenization: Trends and challenges , 2010, J. Comput. Appl. Math..

[21]  Dieter H. Pahr,et al.  From high-resolution CT data to finite element models: development of an integrated modular framework , 2009 .

[22]  Vincent Monchiet,et al.  Combining FFT methods and standard variational principles to compute bounds and estimates for the properties of elastic composites , 2015 .

[23]  David M. Young,et al.  On Richardson's Method for Solving Linear Systems with Positive Definite Matrices , 1953 .

[24]  M. Geers,et al.  A finite element perspective on nonlinear FFT‐based micromechanical simulations , 2016, 1601.05970.

[25]  Gennadi Vainikko,et al.  Periodic Integral and Pseudodifferential Equations with Numerical Approximation , 2001 .

[26]  P. Donato,et al.  An introduction to homogenization , 2000 .

[27]  L. Richardson.  The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[28]  Benoît Bary,et al.  Numerical analysis of linear viscoelastic 3D concrete specimens: Comparison between FE and FFT methods , 2014 .

[29]  Lionel Gélébart,et al.  Filtering material properties to improve FFT-based methods for numerical homogenization , 2014, J. Comput. Phys..

[30]  Luc Dormieux,et al.  FFT-based methods for the mechanics of composites: A general variational framework , 2010 .

[31]  Jonathan J. Hu,et al.  Parallel multigrid smoothing: polynomial versus Gauss--Seidel , 2003 .

[32]  J. Daniel.  The Conjugate Gradient Method for Linear and Nonlinear Operator Equations , 1967 .

[33]  H. Moulinec,et al.  A fast numerical method for computing the linear and nonlinear mechanical properties of composites , 1994 .

[34]  Gene H. Golub,et al.  A generalized conjugate gradient method for the numerical solution of elliptic partial differential equations , 2007, Milestones in Matrix Computation.

[35]  Ricardo A. Lebensohn,et al.  Simulation of micromechanical behavior of polycrystals: finite elements versus fast Fourier transforms , 2009 .

[36]  Z. Strakos,et al.  Distribution of the discretization and algebraic error in numerical solution of partial differential equations , 2014 .

[37]  Martin Diehl,et al.  Numerically robust spectral methods for crystal plasticity simulations of heterogeneous materials , 2013 .

[38]  G. Bonnet,et al.  A polarization‐based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast , 2012 .

[39]  J. Michel,et al.  Effective properties of composite materials with periodic microstructure : a computational approach , 1999 .

[40]  Julien Yvonnet,et al.  A critical comparison of several numerical methods for computing effective properties of highly heterogeneous materials , 2013, Adv. Eng. Softw..

[41]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[42]  M. Geers,et al.  Finite strain FFT-based non-linear solvers made simple , 2016, 1603.08893.

[43]  P. Withers,et al.  Quantitative X-ray tomography , 2014 .

[44]  Hervé Moulinec,et al.  Comparison of three accelerated FFT‐based schemes for computing the mechanical response of composite materials , 2014 .

[45]  F. Willot,et al.  Fourier-based schemes for computing the mechanical response of composites with accurate local fields , 2014, 1412.8398.

[46]  Elio Sacco,et al.  Variational methods for the homogenization of periodic heterogeneous media , 1998 .

[47]  R. Varga,et al.  Chebyshev semi-iterative methods, successive overrelaxation iterative methods, and second order Richardson iterative methods , 1961 .

[48]  Hervé Moulinec,et al.  A Computational Method Based on Augmented Lagrangians and Fast Fourier Transforms for Composites with High Contrast , 2000 .

[49]  G. Milton.  The Theory of Composites , 2002 .

[50]  M. Schneider.  Convergence of FFT‐based homogenization for strongly heterogeneous media , 2015 .

[51]  Jaroslav Vondřejc,et al.  Improved guaranteed computable bounds on homogenized properties of periodic media by the Fourier–Galerkin method with exact integration , 2014 .

[52]  F. Willot,et al.  Fourier‐based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields , 2013, 1307.1015.

[53]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[54]  A. Wathen,et al.  Chebyshev semi-iteration in preconditioning for problems including the mass matrix. , 2008 .

[55]  Michael A. Groeber,et al.  Automated serial sectioning methods for rapid collection of 3-D microstructure data , 2011 .

[56]  Petr Hlaváček,et al.  Inorganic foams made from alkali-activated fly ash: Mechanical, chemical and physical properties , 2015 .

[57]  Martin H. Gutknecht,et al.  The Chebyshev iteration revisited , 2002, Parallel Comput..

[58]  Dierk Raabe,et al.  Comparison of finite element and fast Fourier transform crystal plasticity solvers for texture prediction , 2010 .

[59]  Hervé Moulinec,et al.  Comparison of FFT-based methods for computing the response of composites with highly contrasted mechanical properties , 2003 .

[60]  Jaroslav Vondvrejc.  Improved guaranteed computable bounds on homogenized properties of periodic media by Fourier-Galerkin method with exact integration , 2014, 1412.2033.

[61]  Jaroslav Vondrejc,et al.  Analysis of a Fast Fourier Transform Based Method for Modeling of Heterogeneous Materials , 2011, LSSC.

[62]  C. Robert,et al.  A comparison between different numerical methods for the modeling of polycrystalline materials with an elastic-viscoplastic behavior , 2015 .