paper

An explicit divergence-free DG method for incompressible flow

arXiv:1808.04669 · doi:10.1016/j.cma.2018.11.012

Abstract

We present an explicit divergence-free DG method for incompressible flow based on velocity formulation only. A globally divergence-free finite element space is used for the velocity field, and the pressure field is eliminated from the equations by design. The resulting ODE system can be discretized using any explicit time stepping methods. We use the third order strong-stability preserving Runge-Kutta method in our numerical experiments. Our spatial discretization produces the {\it identical} velocity field as the divergence-conforming DG method of [Cockburn et al., JSC 2007(31), pp.61-73] based on a velocity-pressure formulation, when the same DG operators are used for the convective and viscous parts. Due to the global nature of the divergence-free constraint, there exist no local bases for our finite element space. We present a key result on the efficient implementation of the scheme by identifying the equivalence of the (dense) mass matrix inversion of the globally divergence-free finite element space to a standard (hybrid-)mixed Poisson solver. Hence, in each time step, a (hybrid-)mixed Poisson solver is used, which reflects the global nature of the incompressibility condition. Since we treat viscosity explicitly for the Navier-Stokes equation, our method shall be best suited for unsteady high-Reynolds number flows so that the CFL constraint is not too restrictive.

accepted in CMAME

Cited by in corpus (10)