function out = k_foam_dT(T)
% Derivative with respect of T of the function k_foam.
out = 1.87e-2-10.6e-5*T+14.835e-8*T.^2;