The construction of discretely conservative finite volume schemes that also globally conserve energy or entropy

This work revisits an idea that dates back to the early days of scientific computing, the energy method for stability analysis. It is shown that if the scalar non-linear conservation law partial derivative n/partial derivative t + partial derivative n/partial derivative x f(u) = 0 is approximated by the semi-discrete conservative scheme duj/dt + 1/Delta x (f(j+1/2) - f(j-1/2)) = 0 then the energy of the discrete solution evolves at exactly the same rate as the energy of the true solution, provided that the numerical flux is evaluated by the formula f(j+1/2) = integral(1)(0) f((u) over cap )d theta, where ((u) over cap)theta = uj + theta(u(j+1) - u(j)). With careful treatment of the boundary conditions, this provides a path to the construction of non-dissipative stable discretizations of the governing equations. If shock waves appear in the solution, the discretization must be augmented by appropriate shock operators to account for the dissipation of energy by the shock waves. These results are extended to systems of conservation laws, including the equations of incompressible flow, and gas dynamics. In the case of viscous flow, it is also shown that shock waves can be fully resolved by non-dissipative discretizations of this type with a fine enough mesh, such that the cell Reynolds number <= 2.