The flash method is a commonly used technique for estimating the thermal diffusivity of solid materials. In this work, we develop a new efficient formula for the thermal diffusivity calculation to process the flash method data in the case of heat losses at the rear and front faces. First, the transient heat equation is analytically solved using the Green functions under time-dependent convection boundary conditions. The property of the descending part of the analytical solution, containing series which converges very quickly for long times, allows to express the thermal diffusivity as a function of the slope of the natural logarithm of the temperature rise and the Biot numbers related to heat losses at the rear and front faces. These Biot numbers are calculated from the ratio between the integrals of the complete temperature trends of the two faces. For this purpose, we developed a novel way to evaluate these integrals using the analytical solutions in Laplace domain. To verify the accuracy of these formulas, we simulated experimental data adding a Gaussian noise to the theoretical temperature results of the flash method data.