A three dimensional Kolmogorov problem with extended forcing term for Navier-Stokes equations is considered in the extended periodic domain. The Galerkin{Fourier method is applied and two problem congurations are considered: the symmetry preserving subset of solutions and the full set of solutions. The bifurcation patterns are revealed through the numerical analysis: the eigenvalues of the linearised perturbed system are analyzed (for stationary and periodic solutions) as well as the phase space trajectories that the system generates. Linear stability of the main solution is analyzed. For cubic domain we detected a fast transition to chaos through double Hopf bifurcation. Then the full bifurcation scenario is analyzed for rectangular parallelepiped with double side stretching. The initial stage of laminar-turbulent transition undergoes supercritical pitchfork bifurcation followed by subcritical Hopf bifurction. The system can either go through the series of cycles in Feigenbaum and Sharkovsky ordering or through the bifurcations on invariant tori up to the three and four dimensional tori. The transition to chaos goes either through the emergence of singular attractors of dimension greater than three or through the resonance of tori.