A Monte Carlo code for the calculation of light transport in heterogeneous scattering media is presented together with its validation. Triangle meshes are used to define the interfaces between different materials, in contrast with techniques based on individual volume elements. This approach allows to address realistic problems in a flexible way. A hierarchical spatial organisation enables a fast photon-surface intersection test. The application of the new environment to evaluate the impact of the trabecular structure of bone on its optical properties is demonstrated. A model of the trabecular micro structure recovered from microCT data was used to compute light distribution within tissue. Time-resolved curves across a spherical bone volume were computed. The work presented enables simulation of radiative transport in complex reality-based models of tissue and serves as a powerful, generic tool to study the effect of heterogeneity in the field of biomedical optics.