Improved guaranteed computable bounds on homogenized properties of periodic media by Fourier-Galerkin method with exact integration

Abstract

Moulinec and Suquet introduced FFT-based homogenization in 1994, and twenty years later, their approach is still effective for evaluating the homogenized properties arising from the periodic cell problem. This paper builds on the author's (2013) variational reformulation approximated by trigonometric polynomials establishing two numerical schemes: Galerkin approximation (Ga) and a version with numerical integration (GaNi). The latter approach, fully equivalent to the original Moulinec-Suquet algorithm, was used to evaluate guaranteed upper-lower bounds on homogenized coefficients incorporating a closed-form double grid quadrature. Here, these concepts, based on the primal and the dual formulations, are employed for the Ga scheme. For the same computational effort, the Ga outperforms the GaNi with more accurate guaranteed bounds and more predictable numerical behaviors. Quadrature technique leading to block-sparse linear systems is extended here to materials defined via high-resolution images in a way which allows for effective treatment using the FFT. Memory demands are reduced by a reformulation of the double to the original grid scheme using FFT shifts. Minimization of the bounds during iterations of conjugate gradients is effective, particularly when incorporating a solution from a coarser grid. The methodology presented here for the scalar linear elliptic problem could be extended to more complex frameworks.

32 pages, 9 figures

References in corpus (6)

• Fourier-based schemes for computing the mechanical response of composites with accurate local fields
• Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients
• An FFT-based Galerkin Method for Homogenization of Periodic Media
• 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
• A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media

Cited by in corpus (7)

• Finite strain FFT-based non-linear solvers made simple
• A finite element perspective on non-linear FFT-based micromechanical simulations
• A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media
• FFT-based homogenisation accelerated by low-rank tensor approximations
• Energy-based comparison between the Fourier--Galerkin method and the finite element method
• A stable and accurate X-FFT solver for linear elastic homogenization problems in 3D
• Double-grid quadrature with interpolation-projection (DoGIP) as a novel discretisation approach: An application to FEM on simplexes