An approach for numerical calculation of the transonic three dimensional viscous fluid flow, based on finite volume method, is described in this paper. Stability and fast convergence are obtained in this approach by the introduction of second and fourth order time accurate artificial viscosity terms. A system of differential equations is solved by the application of the LU implicit factorization scheme using Jameson-Yoon flux vector splitting.