A method to predict heat transfer coefficients for metal castings during quenching and/or cooling is provided. First, an initial set of HTC data are obtained from the computational fluid dynamics (CFD) simulation based on the metal casting geometry, initial metal casting temperature (distribution), quench bed/tunnel dimensions and set-up, and a given or baseline (standard) quenching condition including, but not limited to, air and/or gas flow velocity, air and/or gas flow direction relative to the work piece, air and/or gas temperature, air and/or gas humidity, etc. The initial HTC values for the entire surface of the work piece calculated from CFD can then be optimized by multiplying scale factors to minimize the error between the predicted temperature-time profiles and the experimental measurements for the given or standard/baseline quench condition. When the HTC values are optimized for a standard/baseline quench condition, a set of semi-empirical equations (or weight functions) can be used to quickly modify the standard/baseline HTC data for different quenching conditions (i.e., variations of quenching conditions from the baseline) without performing complete heat transferring and optimization calculations. A system and article of manufacture are also provided.