Two flow configurations, involving bilayers, are explored systematically: through an inclined channel comprised of two planar surfaces in parallel alignment and down an inclined plane. Both liquid layers are assumed to be perfectly immiscible and the lower rigid bounding surface contains locally defined steep-sided topographical features - either a step-up/-down or a trench/peak. A common governing equation-set based on the long-wave approximation and depth-averaging is derived, embodying the more general case of a system of n-layers, and solved numerically. Results for the particular case of flow in a vertically aligned channel are found to be indistinguishable from corresponding solutions in the literature; those for the case of a step-up and non-zero Reynolds number having not been reported hitherto. New to this, the case of flow in a inclined channel is investigated as is the situation when, in horizontal alignment, the channel's upper planar bounding surface moves with a constant speed inducing a shear flow. Gravity-driven bilayer film flow reveals a number of interesting features dependent on the fluid properties, the Reynolds number and the ratio of the upper to lower liquid layer thickness, with parallels drawn to the practically relevant cascade/slide-coating mode of multi-layer thin-film deposition. In the limit when both layers have identical properties the corresponding equivalent single layer solution is recovered exactly. (C) 2018 Elsevier Ltd. All rights reserved.