A (Dis)continuous finite element model for generalized 2D vorticity dynamics

E. Bernsen, Onno Bokhove, Jacobus J.W. van der Vegt

    Research output: Book/ReportReportProfessional

    172 Downloads (Pure)


    A mixed continuous and discontinuous Galerkin finite element discretization is constructed for a generalized vorticity streamfunction formulation in two spatial dimensions. This formulation consists of a hyperbolic (potential) vorticity equation and a linear elliptic equation for a (transport) streamfunction. The generalized formulation includes three systems in geophysical fluid dynamics: the incompressible Euler equations, the barotropic quasi-geostrophic equations and the rigid-lid equations. Multiple connected domains are considered with impenetrable and curved boundaries such that the circulation at each connected piece of boundary must be introduced. The generalized system is shown to globally conserve energy and weighted smooth functions of the vorticity. In particular, the weighted square vorticity or enstrophy is conserved. By construction, the spatial finite-element discretization is shown to conserve energy and is $L^2$-stable in the enstrophy norm. The method is verified by numerical experiments which support our error estimates. Particular attention is paid to match the continuous and discontinuous discretization. Hence, the implementation with a third-order Runge-Kutta time discretization conserves energy and is $L^2$-stable in the enstrophy norm for increasing time resolution in multiple connected curved domains.
    Original languageUndefined
    Place of PublicationEnschede
    PublisherUniversity of Twente
    ISBN (Print)169-2690
    Publication statusPublished - 2005

    Publication series

    PublisherDepartment of Applied Mathematics, University of Twente
    ISSN (Print)0169-2690


    • EWI-3607
    • METIS-227276
    • IR-65971
    • MSC-35M10
    • MSC-65M15

    Cite this