Abstract : Radiative transfer plays an important role in turbulent combustion and should be incorporatedin numerical simulations. However, as combustion and radiation are characterized bydifferent time scales and different spatial and chemical treatments, and the complexity of theturbulent combustion flow, radiation effect is often neglected or roughly modelled. Couplinga large eddy simulation combustion solver and a radiation solver through a dedicated languageCORBA is investigated. Four formulations of Monte Carlo method (Forward Method,Emission Reciprocity Method, Absorption Reciprocity Method and Optimized ReciprocityMethod) employed to resolve RTE have been compared in a one-dimensional flame testcase using three-dimensional calculation grids with absorbing and emitting medium in orderto validate the Monte Carlo radiative solver and to choose the most efficient model forcoupling. In order to improve the performance of Monte Carlo solver, two techniques havebeen developed. After that, a new code dedicated to adapt the coupling work has beenproposed. Then results obtained using two different RTE solvers (Reciprocity Monte Carlomethod and Discrete Ordinate Method) applied to a three-dimensional turbulent reactingflow stabilized downstream of a triangular flame holder with a correlated-k distributionmodel describing the real gas medium spectral radiative properties are compared not onlyin terms of physical behavior of the flame but also in computational performance (storagerequirement, CPU time and parallelization efficiency). Finally, the impact of boundary conditionstaking into account the actual wall emissivity and temperature has been discussed.