Hostname: page-component-77c78cf97d-d2fvj Total loading time: 0 Render date: 2026-04-26T06:30:34.991Z Has data issue: false hasContentIssue false

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

Published online by Cambridge University Press:  11 March 2014

Ivo Babuška
Affiliation:
Institute for Computational Engineering and Science and Department of Aerospace Engineering, University of Texas, Austin, TX 78712, USA. babuska@ices.utexas.edu
Xu Huang
Affiliation:
Department of Mathematics, Louisiana State University, Baton Rouge, LA 70803, USA; xhuang4@tigers.lsu.edu
Robert Lipton
Affiliation:
Department of Mathematics and Center for Computation and Technology, Louisiana State University, Baton Rouge, LA 70803, USA; lipton@math.lsu.edu
Get access

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.

Information

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Article purchase

Temporarily unavailable