A multiscale spectral generalized finite element method (MS-GFEM) is presented for the solution of large two and three dimensional stress analysis problems inside heterogeneous media. It can be employed to solve problems too large to be solved directly with FE techniques and is designed for implementation on massively parallel machines. The method is multiscale in nature and uses an optimal family of spectrally defined local basis functions over a coarse grid. It is proved that the method has an exponential rate of convergence. To fix ideas we describe its implementation for a two dimensional plane strain problem inside a fiber reinforced composite. Here fibers are separated by a minimum distance however no special assumption on the fiber configuration such as periodicity or ergodicity is made. The implementation of MS-GFEM delivers the discrete solution operator using the same order of operations as the number of fibers inside the computational domain. This implementation is optimal in that the number of operations for solution is of the same order as the input data for the problem. The size of the MS-GFEM matrix used to represent the discrete inverse operator is controlled by the scale of the coarse grid and the convergence rate of the spectral basis and can be of order far less than the number of fibers. This strategy is general and can be applied to the solution of very large FE systems associated with the discrete solution of elliptic PDE.
Mots clés : generalized finite elements, multiscale method, spectral method, heterogeneous media, fiber reinforced composites
@article{M2AN_2014__48_2_493_0, author = {Babu\v{s}ka, Ivo and Huang, Xu and Lipton, Robert}, title = {Machine {Computation} {Using} the {Exponentially} {Convergent} {Multiscale} {Spectral} {Generalized} {Finite} {Element} {Method}}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {493--515}, publisher = {EDP-Sciences}, volume = {48}, number = {2}, year = {2014}, doi = {10.1051/m2an/2013117}, mrnumber = {3177855}, language = {en}, url = {http://archive.numdam.org/articles/10.1051/m2an/2013117/} }
TY - JOUR AU - Babuška, Ivo AU - Huang, Xu AU - Lipton, Robert TI - Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2014 SP - 493 EP - 515 VL - 48 IS - 2 PB - EDP-Sciences UR - http://archive.numdam.org/articles/10.1051/m2an/2013117/ DO - 10.1051/m2an/2013117 LA - en ID - M2AN_2014__48_2_493_0 ER -
%0 Journal Article %A Babuška, Ivo %A Huang, Xu %A Lipton, Robert %T Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2014 %P 493-515 %V 48 %N 2 %I EDP-Sciences %U http://archive.numdam.org/articles/10.1051/m2an/2013117/ %R 10.1051/m2an/2013117 %G en %F M2AN_2014__48_2_493_0
Babuška, Ivo; Huang, Xu; Lipton, Robert. Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 48 (2014) no. 2, pp. 493-515. doi : 10.1051/m2an/2013117. http://archive.numdam.org/articles/10.1051/m2an/2013117/
[1] Subgrid upscaling and mixed multiscale finite elements. SIAM J. Numer. Anal. 44 (2006) 1150-1171. | MR | Zbl
and ,[2] Homogenization and Its Application. Mathematical and Computational Problems. SYNSPADE 1975, Numer. Solution Part. Differ. Eqs. lll, edited by B. Hubbard. Academic Press (1976) 89-116. | MR | Zbl
,[3] Damage analysis of fiber composites, Part I Statistical analysis on fiber scale. Comput. Methods Appl. Mech Engrg. 172 (1999) 27-77. | MR | Zbl
, , and ,[4] Generalized Finite Element Methods-Main Ideas, Results and Perspective. Int. J. Comput. Methods 1 (2004) 67-103. | Zbl
, and ,[5] Stable Generalized Finiter Element Methods (SGFEM). Comput. Meth. Appl. Mech. Eng. 201-204 (2012) 91-111. | MR | Zbl
and ,[6] Special finite element methods for a class of second order elliptic problems with rough coefficients. SIAM J. Numer. Anal. 31 (1994) 945-981. | MR | Zbl
, and ,[7] Optimal local approximation spaces for Generalized Finite Element Methods with application to multiscale problems. Multiscale Model. Simul., SIAM 9 (2011) 373-406. | MR | Zbl
and ,[8] L2 global to local projection: an approach to multiscale analysis. M3AS 21 (2011) 2211-2226. | MR | Zbl
and ,[9] The Partition of Unity Finite Element Method. Internat. J. Numer. Methods Engrg. 40 (1997) 727-758. | MR | Zbl
and ,[10] J, E. Osborn, Generalized finite element methods:Their performance and their relation to the mixed methods. SIAM, J. Numer. Anal. 20 (1983) 510-536. | MR | Zbl
,[11] Eigenvalue Problems. Handbook of Numerical Analysis, Finite Element Methods (Part 1), Vol. II, edited by P.G. Ciarlet and J.L. Lions. Elsevier Science Publishers, Amsterdam (1991). | MR | Zbl
and ,[12] Homogenization Processes in Periodic Media. Nauka, Moscow (1984). | MR | Zbl
and ,[13] Diffusion and permeation in heterogenous media. Diffusion in Polymers, edited by J. Crank, G.S. Park. Academic Press (1968).
,[14] Existence of ℋ-matrix approximants to the inverse FE-matrix of elliptic operators with L∞ coefficients. Numer. Math. 95 (2003) 1-28. | MR | Zbl
and ,[15] Flux norm approach to finite dimensional homogenization approximations with nonseparated length scales and high contrast. Arch. Rat. Mech. Anal. 198 (2010) 177-221. | MR | Zbl
and ,[16] Asymptotic Analysis for Periodic Structures. North Holland Pub., Amsterdam (1978). | Zbl
, and ,[17] Asymptotic distribution of eigenfunctions and eigenvalues of the basic boundary-contact oscillation problems of the classical theory of elasticity. Georgian Math. J. 6 (1999) 107-126. | MR | Zbl
and ,[18] Critique on theories predicting thermoelastic properties of fibrous composites. J. Comput. Mat. 2 (1968) 332-358.
and ,[19] W. E, B. Engquist, The heterogeneous multiscale methods. Commun. Math. Sci. 1 (2003) 87-132. | MR | Zbl
[20] Analysis of the heterogeneous multiscale method for elliptic homogenization problems. J. Amer. Math. Soc. 18 (2005) 121-156. | MR | Zbl
, and ,[21] Multiscale Finite Element Methods. Springer (2009). | MR | Zbl
and ,[22] Convergence of a nonconforming mutiscale finite element method. SIAM J. Numer. Anal. 37 (2000) 888-910. | MR | Zbl
, and ,[23] Multiscale finite element methods for porous media flows and their applications. Appl. Numer. Math. 57 (2007) 577-596. | MR | Zbl
and ,[24] Asymptotic and numerical homogenization. Acta Numer. 17 (2008) 147-190. | MR | Zbl
and ,[25] Multiscale enrichment based on partition unity. Int. J. Num. Mech. Eng. 62 (2005) 1341-1359. | Zbl
and ,[26] Analysis of Structural Composite Materials. Marcel Dekker, New York (1973).
, and ,[27] The AL basis for the solution of elliptic problems in heterogeneous media. Multiscale Model. Simul. 10 (2012) 245-258. | MR | Zbl
, and ,[28] Composite Materials, in vol II of Mechanics of composite materials, edited by G.P. Sendeckyj. Academic Press (1974).
and ,[29] Theory of Fiber reinforced materials, NASA Report CR-1974 (1972) 1-704.
,[30] A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys. 134 (1997) 169-189. | MR | Zbl
and ,[31] Removing the cell resonance error in the multiscale finite element method via a Petrov-Galerkin formulation. Commun. Math. Sci. 2 (2004) 185-205. | MR | Zbl
, and ,[32] Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulaion, subgrid scale models, bubbles and the origins of stabilized methods. Comput. Methods Appl. Mech. Engrg. 127 (1995) 387-401. | MR | Zbl
,[33] The variational multiscale method. A Paradigm for computational mechanics. Comput. Meth. Appl. Mech. Eng. 166 (1998) 3-24. | MR | Zbl
, , and ,[34] Die Electrizitatskonstante naturlicher und kustlicher Mischkorper. Phys. Zeitschr. XXVII (1926) 115-158. | JFM
,[35] Multiscale methods for elliptic problems. Multiscale Model. Simul. 9 (2011) 1064-1086. | MR | Zbl
,[36] Discussione analitica sul influenze che L'azione di mezo dialettrico hu sulla distribuziione dell' electtricita alla superficie di pin corpi ellecttici diseminati in esso. Mem.Di Math. et di Fisica in Modena 24 (1850) 49.
,[37] Trestise on Electricity and Magnetisum, vol. 1. Oxford Univ. Press (1873) 62. | Zbl
,[38] The Partion of Unity Method Basic Theory and Applications, Comput. Meth. Appl. Mech. Eng. 139 (1996) 289-314. | Zbl
and ,[39] On n-widths for elliptic problems. J. Math. Anal. Appl. 247 (2000) 272-289. | MR | Zbl
,[40] The Theory of Composites. Cambridge University Press, Cambridge (2002). | MR | Zbl
.[41] H-convergence, Séminaire d'Analyse Fonctionelle et Numérique de l'Université d'Alger, mimeographed notes, 1978. L. Tartar Cours Peccot, College de France (1977). Translated into English as F. Murat L. Tartar, H- convergence, in Topics in the Mathematical Modeling of Composite Materials, Progress in Nonlinear Differential Equations and their Applications, in vol. 31, edited by A.V. Cherkaev, R.V. Kohn. Birkhäuser, Boston (1997) 21-43. | Zbl
,[42] A framework for adaptive multiscale methods for elliptic problems. Multiscale Model. Simul. 7 (2008) 171-196. | MR | Zbl
, and ,[43] Metric-based upscaling. Commun. Pure Appl. Math. 60 (2007) 675-723. | MR | Zbl
and ,[44] Localized bases for finite-dimensional homogenization approximations with nonseparated scales and high contrast. Multiscale Model. Simul. 9 (2011) 1373-1398. | MR | Zbl
and ,[45] n-Widths in Approximation Theory. Springer-Verlag, Berlin, Heidelberg, New York 7 (1985). | MR | Zbl
,[46] S.D. Poisson, Second mem. sur la theorie de magnetism, Mem. de L Acad. de France (1822) 5.
[47] On the influence of obstacles in rectangular order upon the properties of the medium. Philos. Mag. 50 (1892) 481. | JFM
,[48] Non-Homogeneous Media and Vibration Theory, in vol. 127 of Lecture Notes in Physics. Springer-Verlag (1980). | MR | Zbl
.[49] Sul limite delle soluzioni di problemi di Cauchy relativi all'equazione del calore. Ann Scu. Norm. Pisa 21 (1967) 657-699. | Numdam | MR | Zbl
,[50] Sulla convergenza di soluzioni di equazioni paraboliche ed ellittiche. Ann. Sc. Norm. Sup. Pisa 22 (1968) 517-597. | Numdam | MR | Zbl
,[51] Convergence in Energy for Elliptic Operators, edited by B. Hubbard. Numer. Solutions Partial Differ Eqs. III, (Synspade 1975, College Park, Maryland 1975). Academic Press, New York (1975). | MR | Zbl
,[52] Variational Methods Applied to Problems of Diffusion and Reaction, Springer Tracts in Natural Philosophy. Springer-Verlag (1973). | Zbl
and ,[53] T. Strouboulis, L. Zhang and I Babuška, Generalized finite element method using mesh-based handbooks application to problem in domains with many voids. Comput. Methods Appl. Mechanics Engrg. 192 (2003) 3109-3161. | MR | Zbl
[54] The design and analysis of the generalized finite element method. Comput. Methods Appl. Mech. Engrg. 181 (2001) 43-69. | MR | Zbl
, and ,[55] p-version of generalized FEM using mesh based handbooks with applications to multiscale problems. Int. J. Num. Meth. Engrg. 60 (2004) 1639-1672. | MR | Zbl
, and ,[56] Random Heterogeneous Materials, Microstructure and Macroscopic Properties. Springer, New York (2002). | MR | Zbl
,Cité par Sources :