Fast isogeometric solvers for hyperbolic wave propagation problems

Abstract: We use the alternating direction method to simulate implicit dynamics. ur spatial discretization uses isogeometric analysis. Namely, we simulate a (hyperbolic) wave propagation problem in which we use tensor-product B-splines in space and an implicit time marching method to fully discretize the problem. We approximate our discrete operator as a Kronecker product of one-dimensional mass and stiffness matrices. As a result of this algebraic transformation, we can factorize the resulting system of equations in linear (i.e., O(N)) time at each step of the implicit method. We demonstrate the performance of our method in the model P-wave propagation problem. We then extend it to simulate the linear elasticity problem once we decouple the vector problem using alternating triangular methods. We proof theoretically and experimentally the unconditional stability of both methods.