We consider the Cauchy problem for the incompressible homogeneous Navier-Stokes (NS) equations on a d-dimensional torus, in the C^infinity formulation described, e.g., in [25]. In [22][25] it was shown how to obtain quantitative estimates on the exact solution of the NS Cauchy problem via the a posteriori analysis of an approximate solution; such estimates concern the interval of existence of the exact solution and its distance from the approximate solution. In the present paper we consider an approximate solutions of the NS Cauchy problem having the form u^N(t) = sum_{j=0}^N R^j u_j(t), where R is the mathematical Reynolds number (the reciprocal of the kinematic viscosity) and the coefficients u_j(t) are determined stipulating that the NS equations be satisfied up to an error O(R^{N+1}). This subject was already treated in [24], where, as an application, the Reynolds expansion of order N=5 in dimension d=3 was considered for the initial datum of Behr-Necas-Wu (BNW). In the present paper, these results are enriched regarding both the theoretical analysis and the applications. Concerning the theoretical aspect, we refine the approach of [24] following [25] and use the symmetries of the initial datum in building up the expansion. Concerning the applicative aspect we consider two more (d=3) initial data, namely, the vortices of Taylor-Green (TG) and Kida-Murakami (KM); the Reynolds expansions for the BNW, TG and KM data are performed via a Python program, attaining orders between N=12 and N=20. Our a posteriori analysis proves, amongst else, that the solution of the NS equations with anyone of the above three data is global if R is below an explicitly computed critical value. Our critical Reynolds numbers are below the ones characterizing the turbulent regime; however these bounds have a sound theoretical support, are fully quantitative and improve previous results of global existence.