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

DOI10.1016/J.JCP.2016.05.041zbMATH Open1349.94072arXiv1508.02045OpenAlexW2298685810MaRDI QIDQ726919FDOQ726919

Jan Zeman, Nachiketa Mishra, Jaroslav Vondřejc

Publication date: 5 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: 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.

Full work available at URL: https://arxiv.org/abs/1508.02045

Recommendations

• Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients
• Analysis of a fast Fourier transform based method for modeling of heterogeneous materials
• An FFT-based Galerkin method for homogenization of periodic media
• Fourier-accelerated nodal solvers (FANS) for homogenization problems
• Locally-synchronous, iterative solver for Fourier-based homogenization

zbMATH Keywords

conjugate gradient algorithmfast Fourier transformnumerical homogenizationRichardson iterationFourier-Galerkin methodChebyshev semi-iterative methodeyre-milton schemeguaranteed bounds

Mathematics Subject Classification ID

Signal theory (characterization, reconstruction, filtering, etc.) (94A12) Numerical methods for discrete and fast Fourier transforms (65T50)

Cites Work

• Methods of conjugate gradients for solving linear systems
• Parallel iterative methods for sparse linear systems
• Title not available (Why is that?)
• Effective properties of composite materials with periodic microstructure: A computational approach
• Chebyshev semi-iterative methods, successive overrelaxation iterative methods, and second order Richardson iterative methods. I, II
• Multi-scale computational homogenization: trends and challenges
• The Theory of Composites
• Effective properties of elastic periodic composite media with fibers
• Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients
• An accelerated FFT algorithm for thermoelastic and non-linear composites
• Title not available (Why is that?)
• Chebyshev semi-iteration in preconditioning for problems including the mass matrix
• Title not available (Why is that?)
• The Conjugate Gradient Method for Linear and Nonlinear Operator Equations
• Digital image-based modeling applied to the homogenization analysis of composite materials
• Comparison of three accelerated FFT-based schemes for computing the mechanical response of composite materials
• Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites
• Parallel multigrid smoothing: Polynomial versus Gauss--Seidel
• A polarization-based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast
• Fourier-based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields
• Filtering material properties to improve FFT-based methods for numerical homogenization
• Distribution of the discretization and algebraic error in numerical solution of partial differential equations
• Composite convergence bounds based on Chebyshev polynomials and finite precision conjugate gradient computations
• The Chebyshev iteration revisited
• Variational methods for the homogenization of periodic heterogeneous media
• Title not available (Why is that?)
• On Richardson's Method for Solving Linear Systems with Positive Definite Matrices
• Analysis of a Fast Fourier Transform Based Method for Modeling of Heterogeneous Materials
• Convergence of FFT‐based homogenization for strongly heterogeneous media
• Finite strain FFT-based non-linear solvers made simple
• Use of composite voxels in FFT-based homogenization
• Combining FFT methods and standard variational principles to compute bounds and estimates for the properties of elastic composites
• An FFT-based Galerkin method for homogenization of periodic media
• Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method
• Improved guaranteed computable bounds on homogenized properties of periodic media by the Fourier-Galerkin method with exact integration
• An advanced approach for the generation of complex cellular material representative volume elements using distance fields and level sets

Cited In (26)

• Computing the effective response of heterogeneous materials with thermomechanically coupled constituents by an implicit fast Fourier transform-based approach
• Deformation patterning in finite-strain crystal plasticity by spectral homogenization with application to magnesium
• Fourier-accelerated nodal solvers (FANS) for homogenization problems
• An FFT-based fast gradient method for elastic and inelastic unit cell homogenization problems
• On polarization-based schemes for the FFT-based computational homogenization of inelastic materials
• A review of nonlinear FFT-based computational homogenization methods
• A memory saving fast A-EFIE solver for modeling low-frequency large-scale problems
• Finite strain FFT-based non-linear solvers made simple
• FFT based numerical homogenization method for porous conductive materials
• FFT-based homogenisation accelerated by low-rank tensor approximations
• FFT-based homogenization of hypoelastic plasticity at finite strains
• A mixed FFT-Galerkin approach for incompressible or slightly compressible hyperelastic solids under finite deformation
• Discrete Fourier transforms and their application to stress—strain problems in composite mechanics: a convergence study
• Locally-synchronous, iterative solver for Fourier-based homogenization
• FFT-based homogenization on periodic anisotropic translation invariant spaces
• Learning Markovian Homogenized Models in Viscoelasticity
• Voxel‐based finite elements with hourglass control in fast Fourier transform‐based computational homogenization
• Eigendecomposition-based convergence analysis of the Neumann series for laminated composites and discretization error estimation
• Energy-based comparison between the Fourier-Galerkin method and the finite element method
• Improved guaranteed computable bounds on homogenized properties of periodic media by the Fourier-Galerkin method with exact integration
• A finite element perspective on nonlinear FFT-based micromechanical simulations
• A reduced order model for geometrically parameterized two-scale simulations of elasto-plastic microstructures under large deformations
• Lippmann-Schwinger solvers for the explicit jump discretization for thermal computational homogenization problems
• Two-stage data-driven homogenization for nonlinear solids using a reduced order model
• Learning constitutive models from microstructural simulations via a non-intrusive reduced basis method: extension to geometrical parameterizations
• Convergence of iterative methods based on Neumann series for composite materials: theory and practice

Uses Software

• FFTHomPy

This page was built for publication: A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q726919)