The present paper shows an in-depth numerical characterisation of the Jet in Hot Co-flow (JHC) configuration using the Reynolds Averaged Navier-Stokes (RANS) modelling with detailed chemistry. The JHC burner emulates the MILD combustion by means of a hot and diluted co-flow and high speed injection. The current investigation focuses on the effect of turbulent combustion models, turbulence model parameters, boundary conditions, multi-component molecular diffusion and kinetic mechanisms on the results. Results show that the approaches used to model the reaction fine structures, namely as Perfectly Stirred Reactors (PSR) or Plug Flow Reactors (PFR), do not have a major impact on the results. Similarly, increasing the complexity of the kinetic mechanism does not lead to major improvements on the numerical predictions. On the other hand, the inclusion of multi-component molecular diffusion helps increasing the prediction accuracy. Three different Eddy Dissipation Concept (EDC) model formulations are compared, showing their interaction with the choice of the C-1 epsilon constant in the k - epsilon turbulence model. Finally, two approaches are benchmarked for turbulence-chemistry interactions, the EDC model and the Partially Stirred Reactor (PaSR) model. (C) 2017 Elsevier Ltd. All rights reserved.