An FFT-based Galerkin method for homogenization of periodic media

In 1994, Moulinec and Suquet introduced an efficient technique for the numerical resolution of the cell problem arising in homogenization of periodic media. The scheme is based on a fixed-point iterative solution to an integral equation of the Lippmann-Schwinger type, with action of its kernel efficiently evaluated by the Fast Fourier Transform techniques. The aim of this work is to demonstrate that the Moulinec-Suquet setting is actually equivalent to a Galerkin discretization of the cell problem, based on approximation spaces spanned by trigonometric polynomials and a suitable numerical integration scheme. For the latter framework and scalar elliptic problems, we prove convergence of the approximate solution to the weak solution, including a-priori estimates for the rate of convergence for sufficiently regular data and the effects of numerical integration. Moreover, we also show that the variational structure implies that the resulting non-symmetric system of linear equations can be solved by the conjugate gradient method. Apart from providing a theoretical support to Fast Fourier Transform-based methods for numerical homogenization, these findings significantly improve on the performance of the original solver and pave the way to similar developments for its many generalizations proposed in the literature.

[1]  Jan Zeman,et al.  Guaranteed bounds on homogenized periodic media by FFT‐based Galerkin method , 2014 .

[2]  Radhi Abdelmoula,et al.  A damage model for crack prediction in brittle and quasi-brittle materials solved by the FFT method , 2012, International Journal of Fracture.

[3]  J. Cea.  Approximation variationnelle des problèmes aux limites , 1964 .

[4]  W. Rudin.  Real and complex analysis , 1968 .

[5]  Lori Graham-Brady,et al.  A stochastic computational method for evaluation of global and local behavior of random elastic media , 2005 .

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

[7]  Renald Brenner,et al.  Computational approach for composite materials with coupled constitutive laws , 2010 .

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

[9]  G. Strang.  VARIATIONAL CRIMES IN THE FINITE ELEMENT METHOD , 1972 .

[10]  A. Rollett,et al.  Fast fourier transform-based modeling for the determination of micromechanical fields in polycrystals , 2011 .

[11]  Hervé Moulinec,et al.  A numerical method for computing the overall response of nonlinear composites with complex microstructure , 1998, ArXiv.

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

[13]  Jiří Němeček,et al.  Methodology for nanoindentation-assisted prediction of macroscale elastic properties of high performance cementitious composites , 2014 .

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

[15]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

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

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

[18]  H. Moulinec,et al.  Multiscale modeling of ice deformation behavior , 2014 .

[19]  Wolfgang H. Müller,et al.  A study of the coarsening in tin/lead solders , 2000 .

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

[21]  V. Zhikov,et al.  Homogenization of Differential Operators and Integral Functionals , 1994 .

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

[23]  L. Gélébart,et al.  Non-linear extension of FFT-based methods accelerated by conjugate gradients to evaluate the mechanical behavior of composite materials , 2013 .

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

[25]  F. Morel,et al.  Experimental and numerical study of the evolution of stored and dissipated energies in a medium carbon steel under cyclic loading , 2013 .

[26]  S. Shtrikman,et al.  On some variational principles in anisotropic and nonhomogeneous elasticity , 1962 .

[27]  Vondřejc Jaroslav,et al.  FFT-based method for homogenization of periodic media: Theory and applications , 2013 .

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

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

[30]  W. Rudin.  Real and complex analysis, 3rd ed. , 1987 .

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

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

[33]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[34]  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.

[35]  Kaushik Bhattacharya,et al.  A model problem concerning recoverable strains of shape-memory polycrystals , 2005, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

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

[37]  Vít S˘milauer,et al.  Microstructure-based micromechanical prediction of elastic properties in hydrating cement paste , 2006 .

[38]  G. Vainikko.  Fast Solvers of the Lippmann-Schwinger Equation , 2000 .

[39]  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 .

[40]  Zdenek P. Bazant,et al.  Identification of Viscoelastic C-S-H Behavior in Mature Cement Paste by FFT-based Homogenization Method , 2010 .

[41]  R. Fedorak,et al.  Trends and challenges , 1996, Proceedings of Nonvolatile Memory Technology Conference.

[42]  A. Finel,et al.  Coupling phase field and viscoplasticity to study rafting in Ni-based superalloys , 2010 .

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