¹⁰The most advanced models currently used are based on the solution of the BTE using a spherical harmonic expansion method [16]. This avoids the time-consuming Monte-Carlo simulations of non-equilibrium distribution functions, and increases the accuracy.