We develop a probabilistic interpretation of local mild solutions of the three dimensional Navier-Stokes equation in the L p spaces, when the initial vorticity field is integrable. This is done by associating a generalized nonlinear diffusion of the McKean-Vlasov type with the solution of the corresponding vortex equation. We then construct trajectorial (chaotic) stochastic particle approximations of this nonlinear process. These results provide the first complete proof of convergence of a stochastic vortex method for the Navier-Stokes equation in three dimensions, and rectify the algorithm conjectured by Esposito and Pulvirenti in 1989. Our techniques rely on a fine regularity study of the vortex equation in the supercritical L p spaces, and on an extension of the classic McKean-Vlasov model, which incorporates the derivative of the stochastic flow of the nonlinear process to explain the vortex stretching phenomenon proper to dimension three.