Normally hyperbolic invariant manifolds (NHIMs) are well-known organizing centers of the dynamics in the phase space of a nonlinear system. Locating such manifolds in systems far from symmetric or integrable, however, has been an outstanding challenge. Here, we develop an automated detection method for codimension-one NHIMs in autonomous dynamical systems. Our method utilizes Stationary Lagrangian Coherent Structures (SLCSs), which are hypersurfaces satisfying one of the necessary conditions of a hyperbolic LCS, and are also quasi-invariant in a well-defined sense. Computing SLCSs provides a quick way to uncover NHIMs with high accuracy. As an illustration, we use SLCSs to locate two-dimensional stable and unstable manifolds of hyperbolic periodic orbits in the classic ABC flow, a three-dimensional solution of the steady Euler equations. Invariant manifolds play an important role in transport phenomena that range from chemical reactions through fluid mixing to celestial mechanics. Such manifolds, however, are often challenging to locate in multi-dimensional systems that are far from integrable and possess no special symmetries. Among these invariant manifolds, normally hyperbolic invariant manifolds are particularly important due to their persistence under small perturbations. This renders them robust with respect to modeling errors and numerical inaccuracies. Here, we develop a method to detect codimension-one, normally hyperbolic invariant manifolds in general autonomous systems. We approximate such manifolds as zero level sets of a scalar function derived from the recent variational theory of hyperbolic Lagrangian Coherent Structures (LCSs). This approximation converges exponentially fast to a true invariant manifold as longer and longer flow-map samples are used in its construction. We illustrate this method on the classic steady ABC flow, revealing some of its two-dimensional invariant manifolds at a previously unseen level of detail.