Abstract We analyze a second-order projection method for the incompressible Navier-Stokes equations on bounded Lipschitz domains. The scheme employs a Backward Differentiation Formula of order two (BDF2) for the time discretization, combined with conforming finite elements in space. Projection methods are widely used to enforce incompressibility, yet rigorous convergence results for possibly non-smooth solutions have so far been restricted to first-order schemes. We establish, for the first time, convergence (up to subsequence) of a second-order projection method to Leray–Hopf weak solutions under minimal assumptions on the data, namely $$u_0 \in L^2_{\textrm{div}}(\Omega )$$ u 0 ∈ L div 2 ( Ω ) and $$f \in L^2(0,T;L^2_{\textrm{div}}(\Omega ))$$ f ∈ L 2 ( 0 , T ; L div 2 ( Ω ) ) . Our analysis relies on two ingredients: A discrete energy inequality providing uniform $$L^\infty (0,T;L^2(\Omega ))$$ L ∞ ( 0 , T ; L 2 ( Ω ) ) and $$L^2(0,T;H^1_0(\Omega ))$$ L 2 ( 0 , T ; H 0 1 ( Ω ) ) bounds for suitable interpolants of the discrete velocities, and a compactness argument combining Simon’s theorem with refined time-continuity estimates. These tools overcome the difficulty that only the projected velocity satisfies an approximate divergence-free condition, while the intermediate velocity is controlled in space. We conclude that a subsequence of the approximations converges to a Leray–Hopf weak solution. This result provides the first rigorous convergence proof for a higher-order projection method under no additional assumptions on the solution beyond those following from the standard a priori energy estimate.