Linear computational cost implicit solver for parabolic problems

Grzegorz Gurgul, Marcin Los, Maciej Paszynski, Victor Calo


In this paper, we use the alternating direction method for isogeometric finite elements to simulate implicit dynamics. Namely, we focus on a parabolic problem and use B-spline basis functions in space and an implicit marching method to fully discretize the problem. We introduce intermediate time steps and separate our differential operator into a summation of the blocks, acting 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 this algebraic transformations, we get a system of linear equations that can be factorized in linear $O(N)$ computational cost in 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.


isogeometric analysis, implicit dynamics, linear computational cost, direct solvers

Full Text:



Y. Bazilevs, L. Beirao da Veiga, J.A. Cottrell, T.J.R. Hughes, and G. Sangalli, Isogeometric analysis: Approximation, stability and error estimates for h-refined meshes, Mathematical Methods and Models in Applied Sciences, 16 (2006) 1031--1090.

Y. Bazilevs, V.M. Calo, J.A. Cottrell, T.J.R. Hughes, A. Reali, G. Scovazzi, Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering 197 (2007) 173-201.

Y. Bazilevs, V.M. Calo, Y. Zhang, T.J.R. Hughes:

Isogeometric fluid-structure interaction analysis with applications to arterial blood flow, Computational Mechanics 38 (2006).

G. Birkhoff, R.S. Varga, D. Young, Alternating direction implicit methods, Advanced Computing 3 (1962) 189–273.

V.M. Calo, N. Brasher, Y. Bazilevs, T.J.R. Hughes, Multiphysics Model for Blood Flow and Drug Transport with Application to Patient-Specific Coronary Artery Flow Computational Mechanics, 43(1) (2008) 161--177.

K. Chang, T.J.R. Hughes, V.M. Calo, Isogeometric variational multiscale large-eddy simulation of fully-developed turbulent flow over a wavy wall Computers and Fluids, 68 (2012) 94-104.

N. Collier, D. Pardo, L. Dalcin, M. Paszy'{n}ski, and V. Calo, The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers, Computer Methods in Applied Mechanics and Engineering (2012), 213, 353-361.

J. A. Cottrell, T. J. R. Hughes, Y. Bazilevs, Isogeometric Analysis: Toward Unification of CAD and FEA, John Wiley and Sons, (2009)

L. Dede, T.J.R. Hughes, S. Lipton, V.M. Calo, Structural topology optimization with isogeometric analysis in a phase field approach USNCTAM2010, 16th US National Congree of Theoretical and Applied Mechanics.

L. Dede, M. J. Borden, T.J.R. Hughes, Isogeometric analysis for topology optimization with a phase field model ICES REPORT 11-29, The Institute for Computational Engineering and Sciences, The University of Texas at Austin (2011).

J. Douglas, H. Rachford, On the numerical solution of heat conduction problems in two and three space variables,

Transactions of American Mathematical Society 82 (1956) 421–439.

R. Duddu, L. Lavier, T.J.R. Hughes, V.M. Calo, A finite strain Eulerian formulation for compressible and nearly incompressible hyper-elasticity using high-order NURBS elements, International Journal of Numerical Methods in Engineering, 89(6) (2012) 762-785.

L. Gao, V.M. Calo, Fast Isogeometric Solvers for Explicit Dynamics, Computer Methods in Applied Mechanics and Engineering, 274 (1) (2014) 19-41.

L. Gao, V.M. Calo, Preconditioners based on the alternating-direction-implicit algorithm for the 2D steady-state diffusion equation with orthotropic heterogeneous coefficients, 273 (1) (2015) 274-295.

L. Gao, Kronecker Products on Preconditioning, PhD. Thesis, King Abdullah University of Science and Technology (2013).

H. G'omez, V.M. Calo, Y. Bazilevs, T.J.R. Hughes, Isogeometric analysis of the Cahn-Hilliard phase-field model, Computer Methods in Applied Mechanics and Engineering 197 (2008) 4333--4352.

H. G'omez, T.J.R. Hughes, X. Nogueira, V.M. Calo, Isogeometric analysis of the isothermal Navier-Stokes-Korteweg} equations. Computer Methods in Applied Mechanics and Engineering 199 (2010) 1828-1840.

J. L. Guermond, P. Minev, A new class of fractional step techniques for the incompressible Navier-Stokes equations using direction splitting, Comptes Rendus Mathematique 348(9-10) (2010) 581–585.

J. L. Guermond, P. Minev, J. Shen, An overview of projection methods for incompressible flows, Computer Methods in Applied Mechanics and Engineering, 195 (2006) 6011–6054.

G. Gurgul, M. Wo'{z}niak, M. L{}o's D. Szeliga, M. Paszy'{n}ski, Open source JAVA implementation of the parallel multi-thread alternating direction isogeometric L2 projections solver for material science simulations, Computer Methods in Material Science (2017)

E. Hairer, G. Wanner, Solving ordinary differential equations II: Stiff and differential-algebraic problems (second ed.), Berlin: Springer-Verlag, section IV.3 (1996)

S. Hossain, S.F.A. Hossainy, Y. Bazilevs, V.M. Calo, T.J.R. Hughes, Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls, Computational Mechanics (2011).

M.-C. Hsu, I. Akkerman, Y. Bazilevs, emph{ High-performance computing of wind turbine aerodynamics using isogeometric analysis, Computers and Fluids, 49(1) (2011) 93-100.

M. L{}o's M. Wo'{z}niak, M. Paszy'{n}ski, L. Dalcin, V.M. Calo, Dynamics with Matrices Possessing Kronecker Product Structure, Procedia Computer Science 51 (2015) 286-295

M. L{}o's M. Paszy'{n}ski, A. Kl{}usek, W. Dzwinel, Application of fast isogeometric L2 projection solver for tumor growth simulations, Computer Methods in Applied Mechanics and Engineering, 316 (2017) 1257-1269.

M. L{}o's M. Wo'{z}niak, M. Paszy'{n}ski, A. Lenharth, K. Pingali, IGA-ADS : Isogeometric Analysis FEM using ADS solver, Computer and Physics Communications, 217 (2017) 99-116.

M. Paszy'{n}ski, Fast solvers for mesh-based computations, Taylor & Francis, CRC Press (2016)

D.W. Peaceman, H.H. Rachford Jr., The numerical solution of parabolic and elliptic differential equations,

Journal of Society of Industrial and Applied Mathematics 3 (1955) 28–41.

L. Piegl, and W. Tiller, emph{The NURBS Book (Second Edition) Springer-Verlag New York, Inc., (1997).

E.L. Wachspress, G. Habetler, An alternating-direction-implicit iteration technique, Journal of Society of Industrial and Applied Mathematics 8 (1960) 403–423.

M. Wo'{z}niak, M. L{}o'{s M. Paszy'{n}ski, L. Dalcin, V. Calo, Parallel fast isogeometric solvers for explicit dynamics, Computing and Informatics (2017)

M. L{}o's, M. Paszy'nski, Applications of Alternating Direction Solver for simulations of time-dependent problems, Computer Science 18(2) (2017) 117-128



  • There are currently no refbacks.