Title
Linear Computational Cost Implicit Solver For Parabolic Problems
Abstract
In this paper, we use the alternating direction method for isogeometric finite elements to simulate transient problems. Namely, we focus on a parabolic problem and use B-spline basis functions in space and an implicit time-marching method to fully discretize the problem. We introduce intermediate time-steps and separate our differential operator into a summation of the blocks that act along a particular coordinate axis in the intermediate time-steps. We show that the resulting stiffness matrix can be represented as a multiplication of two (in 2D) or three (in 3D) multi-diagonal matrices, each one with B-spline basis functions along the particular axis of the spatial system of coordinates. As a result of these algebraic transformations, we get a system of linear equations that can be factorized in a linear O(N) computational cost at every time-step of the implicit method. We use our method to simulate the heat transfer problem. We demonstrate theoretically and verify numerically that our implicit method is unconditionally stable for heat transfer problems (i.e., parabolic). We conclude our presentation with a discussion on the limitations of the method.
Year
DOI
Venue
2020
10.7494/csci.2020.21.3.3824
COMPUTER SCIENCE-AGH
Keywords
DocType
Volume
isogeometric analysis, implicit dynamics, linear computational cost, direct solvers
Journal
21
Issue
ISSN
Citations 
3
1508-2806
0
PageRank 
References 
Authors
0.34
0
4
Name
Order
Citations
PageRank
Grzegorz Gurgul100.34
Marcin Los200.34
Maciej Paszynski319336.89
Victor M. Calo419138.14