A compact high order finite volume scheme for advection-diffusion-reaction equations

We present a new integral representation for the flux of the advection-diffusion-reaction equation, which is based on the solution of a local boundary value problem for the entire equation, including the source term. The flux therefore consists of two parts, corresponding to the homogeneous and particular solution of the boundary value problem. Applying Gauss-Legendre quadrature rules to the integral representation gives the high order finite volume complete flux scheme, which is fourth order accurate for both diffusion dominated and advection dominated flow