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

G.L. Kooij, M.A. Bochev, B.J. Geurts

    Research output: Contribution to journalArticleAcademicpeer-review

    7 Citations (Scopus)

    Abstract

    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.
    Original languageEnglish
    Pages (from-to)229-246
    Number of pages18
    JournalJournal of computational and applied mathematics
    Volume316
    DOIs
    Publication statusPublished - 15 May 2017

    Keywords

    • Parallel computing
    • Exponential integrators
    • Partial differential equations
    • Parallel in time
    • Block Krylov subspace

    Fingerprint Dive into the research topics of 'A block Krylov subspace implementation of the time-parallel Paraexp method and its extension for nonlinear partial differential equations'. Together they form a unique fingerprint.

  • Cite this