A cheap preconditioner based on fast diagonalization method for matrix-free weighted-quadrature isogeometric analysis applied to nonlinear transient heat transfer problems

被引:3
作者
Fuentes, Joaquin Cornejo [1 ]
Elguedj, Thomas [1 ]
Dureisseix, David [1 ]
Duval, Arnaud [1 ]
机构
[1] Univ Lyon, INSA Lyon, CNRS, LaMCoS,UMR5259, F-69621 Villeurbanne, France
关键词
IsoGeometric Analysis; Weighted Quadrature; Matrix-Free; Fast Diagonalization; FINITE-ELEMENTS; NURBS;
D O I
10.1016/j.cma.2023.116157
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
In this paper, IsoGeometric Analysis is applied to a transient nonlinear heat conduction problem using a predictor-corrector scheme. To obtain a better approximation of the solution, we propose the k-refinement method coupled to a Matrix-Free Weighted Quadrature approach which speeds up matrix-vector operations and also reduces memory consumption. The main goal of this work is to introduce a preconditioning strategy to improve the convergence rate of the iterative solver. This method is based on Fast Diagonalization technique which exploits the tensor structure of basis functions. We show that our preconditioner is inexpensive, easy to code and as effective as other more sophisticated methods. It is also robust with respect to polynomial degree, mesh size and time-integration scheme parameters. Numerical experiments for different geometries and materials are presented, which support our theoretical results.& COPY; 2023 Elsevier B.V. All rights reserved.
引用
收藏
页数:22
相关论文
共 31 条
[1]  
Anderson Edward, 1999, LAPACK Users' Guide, Vthird, DOI DOI 10.1137/1.9780898719604
[2]   Efficient matrix computation for tensor-product isogeometric analysis: The use of sum factorization [J].
Antolin, P. ;
Buffa, A. ;
Calabro, F. ;
Martinelli, M. ;
Sangalli, G. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2015, 285 :817-828
[3]  
Barrett R., 1994, TEMPLATES SOLUTION L, DOI DOI 10.1137/1.9781611971538.CH2
[4]  
Bergheau J., 2008, Finite Element Simulation of Heat Transfer, V1st, DOI [10.1002/9780470611418, DOI 10.1002/9780470611418]
[5]   An isogeometric locking-free NURBS-based solid-shell element for geometrically nonlinear analysis [J].
Bouclier, Robin ;
Elguedj, Thomas ;
Combescure, Alain .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2015, 101 (10) :774-808
[6]   A SCALABLE AND ROBUST VERTEX-STAR RELAXATION FOR HIGH-ORDER FEM [J].
Brubeck, Pablo D. ;
Farrell, Patrick E. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2022, 44 (05) :A2991-A3017
[7]   Fast formation of isogeometric Galerkin matrices by weighted quadrature [J].
Calabro, F. ;
Sangalli, G. ;
Tani, M. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2017, 316 :606-622
[8]   The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers [J].
Collier, Nathan ;
Pardo, David ;
Dalcin, Lisandro ;
Paszynski, Maciej ;
Calo, V. M. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2012, 213 :353-361
[9]  
Cottrell J.A., 2009, Isogeometric analysis: toward integration of CAD and FEA, DOI DOI 10.1002/9780470749081
[10]  
Duval A., 2022, YETI YET ANOTHER IGA