A dual-level parallel incompressible flow solver accelerates turbulent flow computations on GPU clusters. This approach solves the pressure Poisson equation with a full-depth, amalgamated parallel geometric multigrid method, and implements a Lagrangian dynamic model for subgrid-scale turbulence modeling.