Optical tomography (OT) is a promising non-intrusive characterization technique of absorbing and scattering media that uses transmitted and/or reflected signals of samples irradiated with visible or near-infrared light. The quality of OT techniques is directly related to the accuracy of their forward models due to the use of inversion algorithms. In this paper, forward models for transient OT approaches are investigated. The system under study involves a one-dimensional absorbing and scattering medium illuminated by a short laser pulse; this problem is solved using a discrete ordinates–finite volume (DO–FV) method in both time and frequency domain. Previous works have shown that time-domain approaches coupled with first order spatial interpolation schemes cannot represent the physics of the problem adequately as transmitted fluxes emerge before the minimal physical time required to leave the medium. In this work, the Van Leer and Superbee flux limiters, combined with the second order Lax–Wendroff scheme, are used in an attempt to prevent this. Results show that despite significant improvement, flux limiters fail to completely eliminate the physically unrealistic behaviour. On the other hand, results for transmittance obtained from the frequency-based method are accurate, without physically unrealistic behaviours at early time periods. The frequency-dependent approach is however computationally expensive, since it requires approximately five times more computational time than its temporal counterpart when used as a forward model for transient OT. On the other hand, the great advantages of the frequency-based approach is that limited windows of temporal signals can be calculated efficiently (in transient OT), and it can also be used as a forward model for steady-state, frequency-based and transient OT techniques.