\[\Delta T_{p} = -\frac{6}{\rho_{ap} V_p c_p d_p} \left( \overline{h}(T_p - T_a) + \sigma\epsilon(T^{4}_{p} - T^{4}_{a}) \right)\Delta t\]