We solve the boundary-value problem of the heat transfer modeling in wet surface grinding, considering a constant heat transfer coefficient over the workpiece surface and a general heat flux profile within the friction zone between wheel and workpiece. We particularize this general solution to the most common heat flux profiles reported in the literature, that is, constant, linear, parabolic, and triangular. For these cases, we propose a fast method for the numerical computation of maximum temperature, in order to avoid the thermal damage of the workpiece. Also, we provide a very efficient method for the numerical evaluation of the transient regime duration (relaxation time).