Systems of nonlinear PDEs arising in multilayer channel flows