% Forced convection from a cylinder
function [h, NuD, ReD] = CylNuT(U,D,mat,Tf)
  [ReD,Pr,kf] = Re(U,D,mat,Tf);                                       
  if min(ReD*Pr) < 0.2 warning = 'CylNuT: ReD*Pr < 0.2' 
  end
  NuD = 0.3+0.62*ReD.^0.5*Pr^(1/3).*(1+(ReD/282000).^0.625).^0.8 ...
       ./(1+(0.4/Pr)^(2/3))^0.25;                                     
  h   = NuD*kf./D;                                                    
end
