A two-phase multiple mapping conditioning/large eddy simulation (MMC-LES) model for turbulent flames with solid fuel pyrolysis is derived, numerically implemented and tested. A hybrid approach is used that consists of an Eulerian LES, a stochastic filtered density function method with MMC closure and an inertial Lagrangian fuel particle scheme. A simultaneous interphase heat and mass transfer model with conserved scalar Spalding transfer numbers and localness in mixture fraction space models the pyrolysis of solid fuel particles and volatile gas mixing in the unresolved diffusive layer around each fuel particle. The MMC-LES model is applied to a laboratory-scale, turbulent, pyrolysing coal flame and the predictions are found to be in good overall agreement with the experimental data and comparable to previously reported LES and carrier-phase direct numerical simulations. This paper focuses on the sensitivity of the predictions to numerical parameters and variations in the models for pyrolysis rate and volatile gas composition. Both single first order reaction rate and competing rate pyrolysis models are used and the volatile composition is modelled as both pure methane and a detailed gaseous mixture. The sensitivity of these variations on the predictions of temperature and reactive species mass fractions is marked. Additionally the sensitivity to the presence of envelope flames around individual coal particles is explored and while there are some minor local variations, the steady-state statistics do not vary much. (C) 2019 The Combustion Institute. Published by Elsevier Inc. All rights reserved.