The purpose of this review is to describe and appraise components of calculation methods, based on the solution of conservation equations in differential form, for the velocity, temperature and concentration fields in turbulent combusting flows. Particular attention is devoted to the combustion models used within these methods and to gaseous-combustion applications. The differential equations are considered first and the implications of conventional (i.e., unweighted) and density-weighted averaging discussed in the contexts of solution methods and physical interpretation. In general, it is concluded that equations should be solved with dependent variables in density-weighted form and the interpretation of measurements requires special care to distinguish between conventionally averaged and density-weighted properties. Finite difference approximations contained within numerical procedures for solving the equations relevant to two- and three-dimensional recirculating flows, such as are to be found in combustion chambers, are considered briefly. It is concluded that computer storage requirements very often preclude the possibility of reducing the numerical error to entirely negligible proportions everywhere. Considerable care must therefore be taken in both specifying a sufficient number and the distribution of mesh points to be used and also in the interpretation of computational results. Turbulence models are also discussed briefly and deficiencies noted. Since many flows are partly controlled by mechanisms other than diffusion and turbulence transport, these deficiencies are of major importance in a limited range of circumstances which are discussed. The various recent methods proposed to represent reaction in turbulent flames are reviewed in relation to diffusion and premixed flames and to flames in which an element of both is present. The application of laminar flame sheet models, chemical equilibrium assumptions, probability density functions of different forms, and truncated series expansion of reaction rate expressions are considered together with the use of probability density function transport equations and their (Monte Carlo) solution. The appraisal is made in relation to presently available results and future requirements and possibilities.