FFT-based homogenisation accelerated by low-rank tensor approximations
arXiv:1902.07455 · doi:10.1016/j.cma.2020.112890
Abstract
Fast Fourier transform (FFT) based methods have turned out to be an effective computational approach for numerical homogenisation. In particular, Fourier-Galerkin methods are computational methods for partial differential equations that are discretised with trigonometric polynomials. Their computational effectiveness benefits from efficient FFT based algorithms as well as a favourable condition number. Here these kind of methods are accelerated by low-rank tensor approximation techniques for a solution field using canonical polyadic, Tucker, and tensor train formats. This reduced order model also allows to efficiently compute suboptimal global basis functions without solving the full problem. It significantly reduces computational and memory requirements for problems with a material coefficient field that admits a moderate rank approximation. The advantages of this approach against those using full material tensors are demonstrated using numerical examples for the model homogenisation problem that consists of a scalar linear elliptic variational problem defined in two and three dimensional settings with continuous and discontinuous heterogeneous material coefficients. This approach opens up the potential of an efficient reduced order modelling of large scale engineering problems with heterogeneous material.
arXiv comment27 pages
References in corpus (15)
- A numerical method for computing the overall response of nonlinear composites with complex microstructure
- 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
- Computational Homogenization of Polycrystals
- Fourier-based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields
- Finite strain FFT-based non-linear solvers made simple
- Spectral tensor-train decomposition
- A finite element perspective on non-linear FFT-based micromechanical simulations
- 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
- Competing damage mechanisms in a two-phase microstructure: how microstructure and loading conditions determine the onset of fracture
- Improved guaranteed computable bounds on homogenized properties of periodic media by Fourier-Galerkin method with exact integration
- Energy-based comparison between the Fourier--Galerkin method and the finite element method
- Double-grid quadrature with interpolation-projection (DoGIP) as a novel discretisation approach: An application to FEM on simplexes