The article presents the results of numerical simulation of gas flows in free fire vortices. The fundamental possibility of physical modeling of the occurrence of concentrated fire vortices is presented in a series of experimental studies carried out at the Joint Institute for High Temperatures of the Russian Academy of Sciences. In the model of a compressible continuous medium for the complete system of Navier-Stokes equations, an initial-boundary-value problem is proposed that describes complex three-dimensional unsteady flows of viscous compressible heat-conducting gas in ascending swirling heat fluxes. Using explicit difference schemes and the proposed initial-boundary conditions, approximate solutions of the Navier-Stokes system of equations are constructed and the thermodynamic characteristics of three-dimensional unsteady gas flows are numerically determined. The calculation results showed that in the process of formation of fiery vortices, it is formally possible to distinguish several stages. The calculations of the thermodynamic characteristics of gas flows during heating of the underlying surface by several local sources showed that the selected mathematical model under the corresponding initial and boundary conditions allows numerical experiments to describe the arising of complex unsteady three-dimensional flows.