This flux is needed for boundary conditions, post-processing and Euler-Euler-like sub-models and functions