×

zbMATH — the first resource for mathematics

Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients. (English) Zbl 1197.65191
Summary: We present a new technique to accelerate the convergence of a FFT-based solver for numerical homogenization of complex periodic media proposed by H. Moulinec and P. Suquet [C. R. Acad. Sci., Paris, Sér. II 318, No. 11, 1417–1423 (1994; Zbl 0799.73077)]. The approach proceeds from discretization of the governing integral equation by the trigonometric collocation method due to G. Vainikko [in: Direct and inverse problems of mathematical physics. Papers presented at special sessions of the ISAAC’97 congress, University of Delaware, Newark, DE, USA, June 2–7, 1997. Dordrecht: Kluwer Academic Publishers. Int. Soc. Anal. Appl. Comput. 5, 423–440 (2000; Zbl 0962.65097)] to give a linear system which can be efficiently solved by conjugate gradient methods. Computational experiments confirm robustness of the algorithm with respect to its internal parameters and demonstrate significant increase of the convergence rate for problems with high-contrast coefficients at a low overhead per iteration.

MSC:
65N35 Spectral, collocation and related methods for boundary value problems involving PDEs
PDF BibTeX Cite
Full Text: DOI
References:
[1] Moulinec, H.; Suquet, P., A fast numerical method for computing the linear and nonlinear mechanical properties of composites, Comptes rendus de l’académie des sciences. Série II, Mécanique, physique, chimie, astronomie, 318, 11, 1417-1423, (1994) · Zbl 0799.73077
[2] Vainikko, G., Fast solvers of the lippmann – schwinger equation, (), 423-440 · Zbl 0962.65097
[3] Greengard, L.; Lee, J., Electrostatics and heat conduction in high contrast composite materials, Journal of computational physics, 211, 1, 64-76, (2006) · Zbl 1129.78005
[4] Eyre, D.J.; Milton, G.W., A fast numerical scheme for computing the response of composites using grid refinement, The European physical journal applied physics, 6, 1, 41-47, (1999)
[5] Michel, J.C.; Moulinec, H.; Suquet, P., A computational method based on augmented Lagrangians and fast Fourier transforms for composites with high contrast, CMES-computer modeling in engineering and sciences, 1, 2, 79-88, (2000)
[6] Moulinec, H.; Suquet, P., Comparison of FFT-based methods for computing the response of composites with highly contrasted mechanical properties, Physica B: condensed matter, 338, 1-4, 58-60, (2003)
[7] Saranen, J.; Vainikko, G., Trigonometric collocation methods with product integration for boundary integral equations on closed curves, SIAM journal on numerical analysis, 33, 4, 1577-1596, (1996) · Zbl 0855.65119
[8] Vinogradov, V.; Milton, G.W., An accelerated FFT algorithm for thermoelastic and non-linear composites, International journal for numerical methods in engineering, 76, 11, 1678-1695, (2008) · Zbl 1195.74302
[9] Golub, G.; Van Loan, C.F., Matrix computations, (1996), The Johns Hopkins University Press Baltimore and London · Zbl 0865.65009
[10] Hestenes, M.R.; Stiefel, E., Methods of conjugate gradients for solving linear systems, Journal of research of the national bureau of standards, 49, 6, 409-463, (1952) · Zbl 0048.09901
[11] Fletcher, R., Conjugate gradient methods for indefinite systems, (), 73-89 · Zbl 0326.65033
[12] J. Vondřejc, Analysis of heterogeneous materials using efficient meshless algorithms: one-dimensional study, Master’s thesis, Czech Technical University in Prague (2009). <http://mech.fsv.cvut.cz/∼vondrejc/download/ING.pdf>.
[13] J. Novák, Calculation of elastic stresses and strains inside a medium with multiple isolated inclusions, in: M. Papadrakakis, B. Topping (Eds.), Proceedings of the Sixth International Conference on Engineering Computational Technology, Stirlingshire, UK, 2008, pp. 16, paper 127. doi:10.4203/ccp.89.127. · Zbl 1171.93304
[14] Kanaun, S., Fast calculation of elastic fields in a homogeneous medium with isolated heterogeneous inclusions, International journal of multiscale computational engineering, 7, 4, 263-276, (2009)
[15] Moulinec, H.; Suquet, P., A numerical method for computing the overall response of nonlinear composites with complex microstructure, Computer methods in applied mechanics and engineering, 157, 1-2, 69-94, (1998) · Zbl 0954.74079
[16] Saad, Y., Iterative methods for sparse linear systems, (2003), Society for Industrial and Applied Mathematics Philadelphia, PA, USA · Zbl 1002.65042
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.