Extension of the complete flux scheme to systems of conservation laws

We present the extension of the complete ¿ux scheme to advection-diffusion-reaction systems. For stationary problems, the ¿ux approximation is derived from a local system boundary value problem for the entire system, including the source term vector. Therefore, the numerical ¿ux vector consists of a homogeneous and an inhomogeneous component, corresponding to the advection-diffusion operator and the source term, respectively. For time-dependent systems, the numerical ¿ux is determined from a quasi-stationary boundary value problem containing the time-derivative in the source term. Consequently, the complete ¿ux scheme results in an implicit semidiscretisation. The complete ¿ux scheme is validated for several test problems. Keywords: Advection-diffusion-reaction systems, ¿ux (vector), ¿nite volume method, integral representation of the ¿ux, Green’s matrix, numerical ¿ux, matrix functions, Péclet matrix.