The quadratic upstream interpolation for convective kinetics (QUICK)
scheme of Leonard (1979) uses a three-point upstream-weighted quadratic
interpolation for cell face values. The face value of φ is obtained from a
quadratic function passing through two bracketing nodes (on each side of the
face) and a node on the upstream side