Moving-boundary lubrication flux (source code)

= Moving-boundary lubrication flux

The flux of the local Couette-Poiseuille profile is
$$
Q=h(U_0+U_1)/2-h^3p_x/(12\mu).
$$