The proposed work is a detailed numerical simulation of three-dimensional subsonic turbulent flow in a channel, with the main focus on symmetric perpendicular jets arising from the walls. The solution of the Favreaveraged Navier-Stokes equations closed by the turbulence model is carried out using an algorithm based on the ENO scheme. To accelerate the convergence of the iterative process, a preconditioning method is used and a transition to a vector of primitive variables is performed. The results of the study are essential for a better understanding of subsonic turbulent flows and may find applications in various fields, including engineering and scientific research. The relevance of the work is highlighted by the development of efficient numerical algorithms capable of solving subsonic three-dimensional Navier-Stokes equations using high-order accuracy schemes, as well as the application of robust turbulence models to analyze supersonic multicomponent flow. The scientific novelty of the work is the successful use of the preconditioning method to accelerate the convergence of the iterative process