On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems

Yannis Voet

Applied Mathematics and Computation, 2023, vol. 436, issue C

Abstract: The finite element method is a well-established method for the numerical solution of partial differential equations (PDEs), both linear and nonlinear. However, the repeated reassemblage of finite element matrices for nonlinear PDEs is frequently pointed out as one of the bottlenecks in the computations. The second bottleneck being the large and numerous linear systems to be solved arising from the use of Newton’s method to solve nonlinear systems of equations. In this paper, we will address the first issue. We will see how under mild assumptions the assemblage procedure may be rewritten using a completely loop-free algorithm. Our approach leads to a small matrix-matrix multiplication for which we may count on highly optimized algorithms.

Keywords: Finite element method; Matrix assembly; Mass matrix; Stiffness matrix; Vectorization (search for similar items in EconPapers)
Date: 2023
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0096300322005902
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:apmaco:v:436:y:2023:i:c:s0096300322005902

DOI: 10.1016/j.amc.2022.127516

Access Statistics for this article

Applied Mathematics and Computation is currently edited by Theodore Simos

More articles in Applied Mathematics and Computation from Elsevier
Bibliographic data for series maintained by Catherine Liu ().