We perform a systematic analysis of flow-like solutions in theories ofEinstein gravity coupled to multiple scalar fields, which arise as holographicRG flows as well as in the context of cosmological solutions driven by scalars.We use the first order formalism and the superpotential formulation to classifysolutions close to generic extrema of the scalar potential, and close to"bounces," where the flow is inverted in some or all directions and thesuperpotential becomes multi-valued. Although the superpotential formulationcontains a large redundancy, we show how this can be completely lift bysuitable regularity conditions. We place the first order formalism in thecontext of Hamilton-Jacobi theory, where we discuss the possibility ofnon-gradient flows and their connection to non-separable solutions of theHamilton-Jacobi equation. We argue that non-gradient flows may be useful in thepresence of global symmetries in the scalar sector.