This paper presents a numerical study on the thermal characteristics of a milling process of titanium alloy with nanofluid minimum-quantity lubrication (MQL). The computational fluid dynamics (CFD) approach is introduced for establishing the numerical model for the nanofluid MQL milling process, and estimated temperatures for pure MQL and for nanofluid MQL using both hexagonal boron nitride (hBN) and nanodiamond particles are compared with the temperatures measured by thermocouples in the titanium alloy workpiece. The estimated workpiece temperatures are similar to experimental ones, and the model is validated.