An h-adaptive local discontinuous Galerkin method for the Navier–Stokes–Korteweg equations

Lulu Tian, Y. Xu, Johannes G.M. Kuerten, Jacobus J.W. van der Vegt

    Research output: Contribution to journalArticleAcademicpeer-review

    9 Citations (Scopus)


    In this article, we develop a mesh adaptation algorithm for a local discontinuous Galerkin (LDG) discretization of the (non)-isothermal Navier–Stokes–Korteweg (NSK) equations modeling liquid–vapor flows with phase change. This work is a continuation of our previous research, where we proposed LDG discretizations for the (non)-isothermal NSK equations with a time-implicit Runge–Kutta method. To save computing time and to capture the thin interfaces more accurately, we extend the LDG discretization with a mesh adaptation method. Given the current adapted mesh, a criterion for selecting candidate elements for refinement and coarsening is adopted based on the locally largest value of the density gradient. Astrategy to refine and coarsen the candidate elements is then provided. We emphasize that the adaptive LDG discretization is relatively simple and does not require additional stabilization. The use of a locally refined mesh in combination with an implicit Runge–Kutta time method is, however, non-trivial, but results in an efficient time integration method for the NSK equations. Computations, including cases with solid wall boundaries, are provided to demonstrate the accuracy, efficiency and capabilities of the adaptive LDG discretizations.
    Original languageUndefined
    Pages (from-to)242-265
    Number of pages24
    JournalJournal of computational physics
    Publication statusPublished - 15 Aug 2016


    • Solid wall boundaries
    • EWI-27381
    • Mesh adaptation
    • Accuracy and stability
    • METIS-319470
    • IR-102077
    • (Non)-isothermal Navier–Stokes–Korteweg equations
    • Local discontinuous Galerkin method

    Cite this