A key issue in dimension reduction of dissipative dynamical systems with
spectral gaps is the identification of slow invariant manifolds. We present
theoretical and numerical results for a variational approach to the problem of
computing such manifolds for kinetic models using trajectory optimization. The
corresponding objective functional reflects a variational principle that
characterizes trajectories on, respectively near, slow invariant manifolds.