Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method

Ivo Babuška; Xu Huang; Robert Lipton

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2014)

  • Volume: 48, Issue: 2, page 493-515
  • ISSN: 0764-583X

Abstract

top
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.

How to cite

top

Babuška, Ivo, Huang, Xu, and Lipton, Robert. "Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method." ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique 48.2 (2014): 493-515. <http://eudml.org/doc/273253>.

@article{Babuška2014,
abstract = {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.},
author = {Babuška, Ivo, Huang, Xu, Lipton, Robert},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {generalized finite elements; multiscale method; spectral method; heterogeneous media; fiber reinforced composites},
language = {eng},
number = {2},
pages = {493-515},
publisher = {EDP-Sciences},
title = {Machine Computation Using the Exponentially Convergent Multiscale Spectral Generalized Finite Element Method},
url = {http://eudml.org/doc/273253},
volume = {48},
year = {2014},
}

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 - Modélisation Mathématique et Analyse Numérique
PY - 2014
PB - EDP-Sciences
VL - 48
IS - 2
SP - 493
EP - 515
AB - 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.
LA - eng
KW - generalized finite elements; multiscale method; spectral method; heterogeneous media; fiber reinforced composites
UR - http://eudml.org/doc/273253
ER -

References

top
  1. [1] T. Arbogast and K.J. Boyd, Subgrid upscaling and mixed multiscale finite elements. SIAM J. Numer. Anal.44 (2006) 1150–1171. Zbl1120.65122MR2231859
  2. [2] I. Babuška, 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. Zbl0346.65064MR502025
  3. [3] I. Babuška, B. Anderson, P. Smith and K. Levin, Damage analysis of fiber composites, Part I Statistical analysis on fiber scale. Comput. Methods Appl. Mech Engrg. 172 (1999) 27–77. Zbl0956.74048MR1685902
  4. [4] I. Babuška, U. Banerjee and J. Osborn, Generalized Finite Element Methods-Main Ideas, Results and Perspective. Int. J. Comput. Methods1 (2004) 67–103. Zbl1081.65107
  5. [5] I. Babuška and U. Banerjee, Stable Generalized Finiter Element Methods (SGFEM). Comput. Meth. Appl. Mech. Eng. 201-204 (2012) 91–111. Zbl1239.74093MR2862999
  6. [6] I. Babuška, G. Caloz and J.E. Osborn, Special finite element methods for a class of second order elliptic problems with rough coefficients. SIAM J. Numer. Anal.31 (1994) 945–981. Zbl0807.65114MR1286212
  7. [7] I. Babuška and R. Lipton, Optimal local approximation spaces for Generalized Finite Element Methods with application to multiscale problems. Multiscale Model. Simul., SIAM 9 (2011) 373–406. Zbl1229.65195MR2801210
  8. [8] I. Babuška and R. Lipton, L2 global to local projection: an approach to multiscale analysis. M3AS 21 (2011) 2211–2226. Zbl1242.65234MR2860673
  9. [9] I. Babuška and J. Melenk, The Partition of Unity Finite Element Method. Internat. J. Numer. Methods Engrg.40 (1997) 727–758. Zbl0949.65117MR1429534
  10. [10] I. Babuška, J, E. Osborn, Generalized finite element methods:Their performance and their relation to the mixed methods. SIAM, J. Numer. Anal. 20 (1983) 510–536. Zbl0528.65046MR701094
  11. [11] I. Babuška and J.E. Osborn, 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). Zbl0875.65087MR1115240
  12. [12] N.S. Bakhvalov and G. Panasenko, Homogenization Processes in Periodic Media. Nauka, Moscow (1984). Zbl0692.73012MR797571
  13. [13] R.M. Barrer, Diffusion and permeation in heterogenous media. Diffusion in Polymers, edited by J. Crank, G.S. Park. Academic Press (1968). 
  14. [14] M. Bebendorf and W. Hackbusch, Existence of ℋ-matrix approximants to the inverse FE-matrix of elliptic operators with L∞ coefficients. Numer. Math. 95 (2003) 1–28. Zbl1033.65100MR1993936
  15. [15] L. Berlyand and H. Owhadi, Flux norm approach to finite dimensional homogenization approximations with nonseparated length scales and high contrast. Arch. Rat. Mech. Anal.198 (2010) 177–221. Zbl1229.35009MR2721592
  16. [16] A. Besounssan, J.L. Lions and G.C. Papanicolau, Asymptotic Analysis for Periodic Structures. North Holland Pub., Amsterdam (1978). Zbl0404.35001
  17. [17] T. Burchuladze and R. Rukhadze, 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. Zbl0926.74032MR1675486
  18. [18] C.C. Chams and G.P. Sendeckij, Critique on theories predicting thermoelastic properties of fibrous composites. J. Comput. Mat.2 (1968) 332–358. 
  19. [19] W. E, B. Engquist, The heterogeneous multiscale methods. Commun. Math. Sci. 1 (2003) 87–132. Zbl1093.35012MR1979846
  20. [20] Weinan E, P. Ming and P. Zhang, Analysis of the heterogeneous multiscale method for elliptic homogenization problems. J. Amer. Math. Soc.18 (2005) 121–156. Zbl1060.65118MR2114818
  21. [21] Y. Efendiev and T.Y. Hou, Multiscale Finite Element Methods. Springer (2009). Zbl1163.65080MR2477579
  22. [22] Y.R. Efendiev, T.Y. Hou and X.H. Wu, Convergence of a nonconforming mutiscale finite element method. SIAM J. Numer. Anal.37 (2000) 888–910. Zbl0951.65105MR1740386
  23. [23] Y. Efendiev and T. Hou, Multiscale finite element methods for porous media flows and their applications. Appl. Numer. Math.57 (2007) 577–596. Zbl1112.76046MR2322432
  24. [24] B. Engquist and P.E. Souganidis, Asymptotic and numerical homogenization. Acta Numer.17 (2008) 147–190. Zbl1179.65142MR2436011
  25. [25] J. Fish and Z. Huan, Multiscale enrichment based on partition unity. Int. J. Num. Mech. Eng.62 (2005) 1341–1359. Zbl1078.74637
  26. [26] S.K. Garg, V. Svalbonas and G.A. Gurtman, Analysis of Structural Composite Materials. Marcel Dekker, New York (1973). 
  27. [27] L. Grasedyck, I. Greff and S. Sauter, The AL basis for the solution of elliptic problems in heterogeneous media. Multiscale Model. Simul.10 (2012) 245–258. Zbl1250.65140MR2902606
  28. [28] L.J. Gurtman and R.H. Krock, Composite Materials, in vol II of Mechanics of composite materials, edited by G.P. Sendeckyj. Academic Press (1974). 
  29. [29] Z. Hashin, Theory of Fiber reinforced materials, NASA Report CR-1974 (1972) 1–704. 
  30. [30] T.Y. Hou and Xiao-Hui Wu, A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys.134 (1997) 169–189. Zbl0880.73065MR1455261
  31. [31] T.Y. Hou, Xiao-Hui Wu and Yu Zhang, Removing the cell resonance error in the multiscale finite element method via a Petrov-Galerkin formulation. Commun. Math. Sci.2 (2004) 185–205. Zbl1085.65109MR2119937
  32. [32] T.J.R. Hughes, 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. Zbl0866.76044MR1365381
  33. [33] T.J.R. Hughes, G.R. Feijoo, L. Mazzei and J.B. Quincy, The variational multiscale method. A Paradigm for computational mechanics. Comput. Meth. Appl. Mech. Eng. 166 (1998) 3–24. Zbl1017.65525MR1660141
  34. [34] K. Lichtenecker, Die Electrizitatskonstante naturlicher und kustlicher Mischkorper. Phys. Zeitschr. XXVII (1926) 115–158. Zbl52.0914.17JFM52.0914.17
  35. [35] A. Malquist, Multiscale methods for elliptic problems. Multiscale Model. Simul.9 (2011) 1064–1086. Zbl1248.65124MR2831590
  36. [36] G.F. Masotti, 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. [37] G. Maxwell, Trestise on Electricity and Magnetisum, vol. 1. Oxford Univ. Press (1873) 62. Zbl1049.01021
  38. [38] J. Melenk and I. Babuška, The Partion of Unity Method Basic Theory and Applications, Comput. Meth. Appl. Mech. Eng.139 (1996) 289–314. Zbl0881.65099
  39. [39] J.M. Melenk, On n-widths for elliptic problems. J. Math. Anal. Appl.247 (2000) 272–289. Zbl0963.35047MR1766938
  40. [40] G.W. Milton. The Theory of Composites. Cambridge University Press, Cambridge (2002). Zbl0993.74002MR1899805
  41. [41] F. Murat, 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. Zbl0920.35019
  42. [42] J. Nolen, G. Papanicolaou and O. Pironneau, A framework for adaptive multiscale methods for elliptic problems. Multiscale Model. Simul.7 (2008) 171–196. Zbl1160.65342MR2399542
  43. [43] H. Owhadi and L. Zhang, Metric-based upscaling. Commun. Pure Appl. Math.60 (2007) 675–723. Zbl1190.35070MR2292954
  44. [44] H. Owhadi and L. Zhang, Localized bases for finite-dimensional homogenization approximations with nonseparated scales and high contrast. Multiscale Model. Simul.9 (2011) 1373–1398. Zbl1244.65140MR2861243
  45. [45] A. Pinkus, n-Widths in Approximation Theory. Springer-Verlag, Berlin, Heidelberg, New York 7 (1985). Zbl0551.41001MR774404
  46. [46] S.D. Poisson, Second mem. sur la theorie de magnetism, Mem. de L Acad. de France (1822) 5. 
  47. [47] J.W. Rayleigh, On the influence of obstacles in rectangular order upon the properties of the medium. Philos. Mag. 50 (1892) 481. JFM24.1015.02
  48. [48] E. Sanchez-Palencia. Non-Homogeneous Media and Vibration Theory, in vol. 127 of Lecture Notes in Physics. Springer-Verlag (1980). Zbl0432.70002MR578345
  49. [49] S. Spagnolo, Sul limite delle soluzioni di problemi di Cauchy relativi all’equazione del calore. Ann Scu. Norm. Pisa21 (1967) 657–699. Zbl0153.42103MR225015
  50. [50] S. Spagnolo, Sulla convergenza di soluzioni di equazioni paraboliche ed ellittiche. Ann. Sc. Norm. Sup. Pisa22 (1968) 517–597. Zbl0174.42101MR240443
  51. [51] S. Spagnolo, 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). Zbl0347.65034MR477444
  52. [52] W. Streider and R. Aris, Variational Methods Applied to Problems of Diffusion and Reaction, Springer Tracts in Natural Philosophy. Springer-Verlag (1973). Zbl0273.49067
  53. [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. Zbl1054.74059MR2007029
  54. [54] T. Strouboulis, I. Babuška and K. Copps, The design and analysis of the generalized finite element method. Comput. Methods Appl. Mech. Engrg.181 (2001) 43–69. Zbl0983.65127MR1734667
  55. [55] T. Strouboulis, L. Zhang and I. Babuška, p-version of generalized FEM using mesh based handbooks with applications to multiscale problems. Int. J. Num. Meth. Engrg.60 (2004) 1639–1672. Zbl1059.65106MR2069141
  56. [56] S. Torquato, Random Heterogeneous Materials, Microstructure and Macroscopic Properties. Springer, New York (2002). Zbl0988.74001MR1862782

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.