Time-Discrete Higher-Order ALE Formulations: Stability

Andrea Bonito, Irene Kyza, Ricardo H. Nochetto

    Research output: Contribution to journalArticle

    16 Citations (Scopus)

    Abstract

    Arbitrary Lagrangian Eulerian (ALE) formulations deal with PDEs on deformable domains upon extending the domain velocity from the boundary into the bulk with the purpose of keeping mesh regularity. This arbitrary extension has no effect on the stability of the PDE but may influence that of a discrete scheme. We examine this critical issue for higher-order time stepping without space discretization. We propose time-discrete discontinuous Galerkin (dG) numerical schemes of any order for a time-dependent advection-diffusion-model problem in moving domains, and study their stability properties. The analysis hinges on the validity of the Reynolds' identity for dG. Exploiting the variational structure and assuming exact integration, we prove that our conservative and nonconservative dG schemes are equivalent and unconditionally stable. The same results remain true for piecewise polynomial ALE maps of any degree and suitable quadrature that guarantees the validity of the Reynolds' identity. This approach generalizes the so-called geometric conservation law to higher-order methods. We also prove that simpler Runge-Kutta-Radau methods of any order are conditionally stable, that is, subject to a mild ALE constraint on the time steps. Numerical experiments corroborate and complement our theoretical results.

    Original languageEnglish
    Pages (from-to)577-604
    Number of pages28
    JournalSIAM Journal on Numerical Analysis
    Volume51
    Issue number1
    DOIs
    Publication statusPublished - 2013

    Keywords

    • SCHEMES
    • CONVECTION-DIFFUSION EQUATION
    • discrete Reynolds' identities
    • geometric conservation law
    • dG methods in time
    • material derivative
    • ALE formulations
    • domain velocity
    • moving domains
    • stability
    • FLUID
    • GEOMETRIC CONSERVATION LAW

    Cite this