An FFT-based Galerkin Method for Homogenization of Periodic Media

Abstract

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 setting, 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.

22 pages, 1 figure

References in corpus (4)

• A numerical method for computing the overall response of nonlinear composites with complex microstructure
• Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients
• Fourier-based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields
• Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method

Cited by in corpus (35)

• Fourier-based schemes for computing the mechanical response of composites with accurate local fields
• Computational Homogenization of Polycrystals
• Finite strain FFT-based non-linear solvers made simple
• A finite element perspective on non-linear FFT-based micromechanical simulations
• Simulation of the Hall-Petch effect in FCC polycrystals by means of strain gradient crystal plasticity and FFT homogenization
• The universal emergence of self-affine roughness from deformation
• Three-dimensional microstructure generation using generative adversarial neural networks in the context of continuum micromechanics
• DBFFT: A displacement based FFT approach for non-linear homogenization of the mechanical behavior
• Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method
• A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media
• Accelerated computational micromechanics
• An algorithm for stress and mixed control in Galerkin based FFT homogenization
• An FFT framework for simulating non-local ductile failure in heterogeneous materials
• Improved guaranteed computable bounds on homogenized properties of periodic media by Fourier-Galerkin method with exact integration
• Shape optimization of phononic band gap structures using the homogenization approach
• A deep learning driven pseudospectral PCE based FFT homogenization algorithm for complex microstructures
• FFT-based homogenisation accelerated by low-rank tensor approximations
• An FFT-based framework for predicting corrosion-driven damage in fractal porous media
• Equivariant geometric learning for digital rock physics: estimating formation factor and effective permeability tensors from Morse graph
• Elimination of ringing artifacts by finite-element projection in FFT-based homogenization
• Numerical analysis of several FFT-based schemes for computational homogenization
• Energy-based comparison between the Fourier--Galerkin method and the finite element method
• A fast algorithm for semi-analytically solving the homogenization boundary value problem for block locally-isotropic heterogeneous media
• A stochastic discrete slip approach to microplasticity: Application to submicron W pillars
• Computation of effective elastic moduli of rocks using hierarchical homogenization
• Effective Response of Heterogeneous Materials using the Recursive Projection Method
• A Framework for FFT-based Homogenization on Anisotropic Lattices
• A stable and accurate X-FFT solver for linear elastic homogenization problems in 3D
• Representative volume elements for matrix-inclusion composites -- a computational study on periodizing the ensemble
• Double-grid quadrature with interpolation-projection (DoGIP) as a novel discretisation approach: An application to FEM on simplexes
• An FFT-based approach for Bloch wave analysis: application to polycrystals
• Adaptation and validation of FFT methods for homogenization of lattice based materials
• FFT-based homogenization on periodic anisotropic translation invariant spaces
• Approximation of Periodic PDE Solutions with Anisotropic Translation Invariant Spaces
• Efficient topology optimization using compatibility projection in micromechanical homogenization