SUMMARYMolecular imaging possesses the ability to characterize and measure biological processes at the cellular and molecular level in vivo. As one of the new molecular imaging modalities, bioluminescence tomography (BLT) is to reconstruct the light distribution inside a small animal from the photon flux measured on its surface. To obtain accurate and robust reconstruction, it needs a good understanding of the propagation of photons in biological tissues, which is referred to as the forward problem in BLT. Because the turbid media is high scattering and low absorption, the propagation process can be described by the steady diffusion equation. In this paper, an hp-adaptivity method for BLT forward problem is presented based on finite elements of high orders and moderate meshes by finite element method (FEM). Both numerical simulation and physical experiment are performed to evaluate the accuracy of solution and the efficiency of computation. The relevant results show that element order is more critical than mesh size to produce an accurate FEM solution efficiently.