Engquist-Osher method

The Engquist-Osher method is the finite volume method using the Engquist-Osher flux. Under $(k/d)\sup |f\prime |\le 1$ on the invariant state interval, its update is monotone. It preserves the state range, is contractive for summable perturbations in $L^1$, and decreases discrete total variation.

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$.