This article proposes a numerical method for solving the problem of an axisymmetric methane jet propagation in an infinite wake air flow. The dimensionless equations of the turbulent boundary layer of reacting gases in von Mises coordinates and the k-e turbulence model were used in modeling. The equations for N components of the gas mixture were reduced to two equations by introducing the Schwab-Zeldovich functions. To solve the problem in von Mises coordinates, a two-layer, six-point implicit finite-difference scheme was used, which provided the second order of accuracy of the approximation in coordinates. An iterative process was realized due to the nonlinearity of the equations for the conservation and transfer of substances. The effect of the radius of fuel nozzle on the indices of turbulent jet and flame was investigated. It was found that in an infinite wake flow of fuel with a decrease in the nozzle radius, the rate of chemical reaction and the highest temperature in the calculation area decrease, and the amount of unburnt fuel increases.