A block Krylov subspace implementation of the time-parallel Paraexp method and its extension for nonlinear partial differential equations

详细信息    查看全文

• 作者：G.L. Kooij^a ; ^{g.l.kooij@utwente.nl} ; M.A. Botchev^b ; ^c ; ^{m.a.botchev@utwente.nl} ; B.J. Geurts^a ; ^d ; ^{b.j.geurts@utwente.nl}
• 关键词：Parallel computing ; Exponential integrators ; Partial differential equations ; Parallel in time ; Block Krylov subspace
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2017
• 出版时间：15 May 2017
• 年：2017
• 卷：316
• 期：Complete
• 页码：229-246
• 全文大小：894 K
• 卷排序：316

文摘

A parallel time integration method for nonlinear partial differential equations is proposed. It is based on a new implementation of the Paraexp method for linear partial differential equations (PDEs) employing a block Krylov subspace method. For nonlinear PDEs the algorithm is based on our Paraexp implementation within a waveform relaxation. The initial value problem is solved iteratively on a complete time interval. Nonlinear terms are treated as source terms, provided by the solution from the previous iteration. At each iteration, the problem is decoupled into independent subproblems by the principle of superposition. The decoupled subproblems are solved fast by exponential integration, based on a block Krylov method. The new time integration is demonstrated for the one-dimensional advection–diffusion equation and the viscous Burgers equation. Numerical experiments confirm excellent parallel scaling for the linear advection–diffusion problem, and good scaling in case the nonlinear Burgers equation is simulated.