Engquist-Osher flux ($F_{\mathrm{EO}}(a,b)$)

The Engquist-Osher numerical flux is $F(a,b)=f(0)+{\int }_0^a\max (f\prime (s),0)ds+{\int }_0^b\min (f\prime (s),0)ds$. It sends positive-speed contributions from the left state and negative-speed contributions from the right state. Its consistency $F(u,u)=f(u)$ and opposite monotonicities in the two arguments yield a monotone conservative scheme under the appropriate Courant–Friedrichs–Lewy condition. For the Inviscid Burgers equation, it is $\frac{1}{2}\max (a,0{)}^2+\frac{1}{2}\min (b,0{)}^2$.

For a flux that is a differentiable convex function with unique minimum at $s_{\ast }$, define $f^{+}(u)=f(u)-f(s_{\ast })$ for $u>s_{\ast }$ and zero otherwise, and $f^{-}(u)=f(u)-f(s_{\ast })$ for $u<s_{\ast }$ and zero otherwise. The Engquist-Osher flux is $F(a,b)=f(s_{\ast })+f^{+}(a)+f^{-}(b)$. Thus transsonic rarefaction wave states $a\le s_{\ast }\le b$ give $F=f(s_{\ast })$, while reversed transsonic states give $F=f(a)+f(b)-f(s_{\ast })$. This follows by integrating the positive and negative parts of $f^{\prime }$.