Abstract : Abstract: In this work, we propose a numerical method to handle discontinuous fluxes arising in transport-like equations. More precisely, we study hyperbolic PDEs with flux transmission conditions at interfaces between subdomains where coefficients are discontinuous. A dedicated finite volume scheme with a limited high order enhancement is adapted to treat the discontinuities arising at interfaces. The validation of the method is done on 1D and 2D toy problems for which exact solutions are available, allowing us to do a thorough convergence study. We then apply the method to a biological model focusing on complex cell dynamics, that initially motivated this study, and illustrates the full potentialities of the scheme.