Solutions to conservation laws satisfy the monotonicity property: the number of local extrema is a non-increasing function of time, and local maximum/minimum values decrease/increase monotonically in time. This paper investigates this property from a numerical standpoint. We introduce a class of fully discrete in space and time, high order accurate, difference schemes, called generalized monotone schemes. Convergence toward the entropy solution is proven via a new technique of proof, assuming that the initial data has a finite number of extremum values only, and the flux-function is strictly convex. We define discrete paths of extrema by tracking local extremum values in the approximate solution. In the course of the analysis we establish the pointwise convergence of the trace of the solution along a path of extremum. As a corollary, we obtain a proof of convergence for a MUSCL-type scheme being second order accurate away from sonic points and extrema.

We prove the convergence of a discontinuous Galerkin method approximating the 2-D incompressible Euler equations with discontinuous initial vorticity:$omega_0 in L^2(Omega)$. Furthermore, when $omega_0 in L^infty(Omega)$, the whole sequence is shown to be strongly convergent. This is the first convergence result in numerical approximations of this general class of discontinuous flows. Some important flows such as vortex patches belong to this class.

In this paper we introduce a high-order discontinuous Galerkin method for two-dimensional incompressible flow in the vorticity stream-function formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method. The stream function is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability. The method is efficient for inviscid or high Reynolds number flows. Optimal error estimates are proved and verified by numerical experiments.

A new formulation, a gauge formulation of the incompressible Navier-Stokes equations in terms of an auxiliary field $\mathbf{a}$ and a gauge variable $\phi$, $\mathbf{u} =\mathbf{a}+\nabla\phi$, was proposed recently by E and Liu. This paper provides a theoretical analysis of their formulation and verifies the computational advantages. We discuss the implicit gauge method, which uses backward Euler or Crank-Nicolson in time discretization. However, the boundary conditions for the auxiliary field $\mathbf{a}$ are implemented explicitly (vertical extrapolation). The resulting momentum equation is decoupled from the kinematic equation, and the computational cost is reduced to solving a standard heat and Poisson equation. Moreover, such explicit boundary conditions for the auxiliary field $\mathbf{a}$ will be shown to be unconditionally stable for Stokes equations. For the full nonlinear Navier-Stokes equations the time stepping constraint is reduced to the standard CFL constraint ${\triangle t} / {\triangle x} \le C$. We also prove first order convergence of the gauge method when we use MAC grids as our spatial discretization. The optimal error estimate for the velocity field is also obtained.

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J-G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher-order) finite elements. This method can achieve high-order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright