The proposed FEM model describes natural convective air cooling of cheese curd and cheeses with different sizes, chemical composition and initial temperature, including temperature-dependent functions accounting for the variation of specific heat capacity and thermal conductivity of cheeses. Both the calculated convective heat transfer coefficients (from 3.58 to 15.15 W/m(2) K) and the ratio between Grashof and square Reynolds numbers confirmed that heat exchange was natural convective. The model permitted to accurately predict the transient temperature change into the cheese, as shown by the mean RMSEs values (from 0.34 to 2.29 degrees C). Higher RMSEs values (up to 3.29 degrees C) were obtained for cheese curds, because some deviations from the assumptions of the model occurred. These higher RMSEs values for cheese curds quantify the importance of compliance with the model's assumptions to ensure a best fit between the simulated and experimental data. (C) 2011 Elsevier Ltd. All rights reserved.