A numerical model based on computational fluid dynamics is developed to analyze fluid flow and thermal aspects in grinding. The model uses multiphase fluid flow with heat transfer based on the volume-of-fluid method, convection, conduction in solids and a multiple reference frame model of the porous grinding zone. Fluid velocity vectors, useful flow rate, grinding temperatures and energy partition are predicted using the model. In lieu of direct measurements of these quantities, the verification relies on the indirect assessment of surface integrity. The simulation results provide adequate agreement with the measured residual stress, depth of heat-affected zone and full width at half maximum profile with respect to the grinding temperatures.