The version of simulating fans in mine ventilation network analysis by using the combination of FEM with the laminar solution method is considered. Two ways of fan insertion are discussed: in the form of a separate active element and in the form of combined element (an active element inserted into a passive element). The latter approach allows changing a position of the fan in a design stage without renumbering of a grid. The method of chords is used to linearize in general nonlinear fan characteristic. It is shown that the dependence of air discharge on the upstream-downstream pressure difference is an ascending curve for the fan. This provides a convergence for the iterative procedure of solving the FEM nonlinear equations system. A method of calculation of convective heat transfer between the air flow and the rock wall based on Lyon’s solution for a pipe with the constant density of boundary heat flow is described. As the density of heat flow and the air properties are slowly varying functions along a mine working (a branch of a network), then discretizing each branch into small segments permits to consider these parameters as constant along segments. Calculation of nonhomogeneous temperature field and treatment of gravity allow modelling the natural draught effect. The results of solution of the test problems validating used algorithms are given. Calculation results of flow mass rate in a mine network within the year are also represented which make it possible to estimate the influence of natural draught on venting capability.