The newly developed unifying discontinuous formulation named the correction procedure viareconstruction (CPR) for conservation laws is extended to solve the Navier-Stokesequations for 3D mixed grids. In the current development, tetrahedrons and triangularprisms are considered. The CPR method can unify several popular high order methodsincluding the discontinuous Galerkin and the spectral volume methods into a more efficientdifferential form. By selecting the solution points to coincide with the flux points,solution reconstruction can be completely avoided. Accuracy studies confirmed that theoptimal order of accuracy can be achieved with the method. Several benchmark test casesare computed by solving the Euler and compressible Navier-Stokes equations to demonstrateits performance.