Vortex identifications in turbulent flows past arrays of tandem plates are performed by employing the velocity field obtained by high-fidelity large eddy simulations. Lagrangian coherent structures (LCSs) are extracted to examine the evolution and the nonlinear interaction of vortices and to characterize the spatial and temporal characteristics of the flow. LCSs’ identification method is based on the finite-time Lyapunov exponent (FTLE), which is evaluated using the instantaneous velocity data. The simulations are performed in three-dimensional geometries to understand the physics of fluid motion and the vortex dynamics in the vicinity of plates and surfaces at Reynolds number of 50 000. The instantaneous vorticity fields, Eulerian Q-criterion, and LCSs are presented to interpret and understand complex turbulent flow structures. The three-dimensional FTLE fields provide valuable information about the vortex generation, spatial location, evolution, shedding, decaying, and dissipation of vortices. It is demonstrated here that FTLE can be used together with Eulerian vortex identifiers to characterize the turbulent flow field effectively.