Topology optimization aims to find the best material layout subject to given constraints. The so-called material distribution methods cast the governing equation as an extended or fictitious domain problem, in which a coefficient field represents the design. When solving the governing equation using the finite element method, a large number of elements are used to discretize the design domain, and an element-wise constant function approximates the coefficient field in the considered design domain. This article presents a spectral analysis of the (large) coefficient matrices associated with the linear systems stemming from the finite element discretization of a linearly elastic problem for an arbitrary coefficient field. Based on the spectral information, we design a multigrid method which turns out to be optimal, in the sense that the (arithmetic) cost for solving the related linear systems, up to a fixed desired accuracy, is proportional to the matrix-vector cost, which is linear in the corresponding matrix size. The method is tested, and the numerical results are very satisfactory in terms of linear cost and number of iterations, which is bounded by a constant independent of the matrix size.

Nguyen, Q., Serra-Capizzano, S., Tablino Possio, C., Wadbro, E. (2022). Spectral analysis of the finite element matrices approximating 2D linearly elastic structures and multigrid proposals. NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 29(4) [10.1002/nla.2433].

Spectral analysis of the finite element matrices approximating 2D linearly elastic structures and multigrid proposals

Tablino Possio, C;
2022

Abstract

Topology optimization aims to find the best material layout subject to given constraints. The so-called material distribution methods cast the governing equation as an extended or fictitious domain problem, in which a coefficient field represents the design. When solving the governing equation using the finite element method, a large number of elements are used to discretize the design domain, and an element-wise constant function approximates the coefficient field in the considered design domain. This article presents a spectral analysis of the (large) coefficient matrices associated with the linear systems stemming from the finite element discretization of a linearly elastic problem for an arbitrary coefficient field. Based on the spectral information, we design a multigrid method which turns out to be optimal, in the sense that the (arithmetic) cost for solving the related linear systems, up to a fixed desired accuracy, is proportional to the matrix-vector cost, which is linear in the corresponding matrix size. The method is tested, and the numerical results are very satisfactory in terms of linear cost and number of iterations, which is bounded by a constant independent of the matrix size.
Articolo in rivista - Articolo scientifico
finite element approximations; matrix sequences; spectral analysis;
English
17-feb-2022
2022
29
4
e2433
open
Nguyen, Q., Serra-Capizzano, S., Tablino Possio, C., Wadbro, E. (2022). Spectral analysis of the finite element matrices approximating 2D linearly elastic structures and multigrid proposals. NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 29(4) [10.1002/nla.2433].
File in questo prodotto:
File Dimensione Formato  
Nguyen-2022-Numer Linear Algebra Appl-VoR.pdf

accesso aperto

Descrizione: Research Article
Tipologia di allegato: Publisher’s Version (Version of Record, VoR)
Dimensione 2.95 MB
Formato Adobe PDF
2.95 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10281/394986
Citazioni
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 1
Social impact