A recently developed parallel three-dimensional flow solver of the turbulent Navier-Stokes equations is presented. The space discretization of the continuous equations is performed by means of the Fluctuation Splitting algorithm, and the implicit integration process of the discrete equations is based on an effectively preconditioned Newton-Krylov subspace algorithm. The high computational performance is demonstrated using a suite of two- and three-dimensional test cases with grids featuring up to 6.5 million nodes.