dyna DYNA Dyna rev.fac.nac.minas 0012-7353 Universidad Nacional de Colombia 10.15446/dyna.v86n208.72603 Articulos Influencia de la formulación en las propiedades termofísicas y parámetros de calidad de un helado de crema de leche Influence of the formulation on the thermophysical properties and the quality parameters of dairy ice cream Zambrano-Mayorga Luis F. a Ramírez-Navas Juan Sebastián a b Ochoa-Martínez Claudia I. a Department of Food Engineering, GIPAB, Universidad del Valle, Cali, Colombia. zambrano.luis@correounivalle.edu.co, juan.sebastian.ramirez@correounivalle.edu.co, claudia.ochoa@correounivalle.edu.co Universidad Santiago de Cali Department of Food Engineering GIPAB Universidad del Valle Cali Colombia zambrano.luis@correounivalle.edu.co juan.sebastian.ramirez@correounivalle.edu.co claudia.ochoa@correounivalle.edu.co Department of Chemistry, GIEMA, Universidad Santiago de Cali, Cali, Colombia Universidad Santiago de Cali Department of Chemistry GIEMA Universidad Santiago de Cali Cali Colombia

© The author; licensee Universidad Nacional de Colombia.

Jan-Mar 2019 86 208 117 125 01 06 2018 15 06 2019 28 01 2018 This is an open-access article distributed under the terms of the Creative Commons Attribution License Abstract

Thermophysical properties and the quality parameters of vanilla dairy ice cream were evaluated as a function of formulation. The sucrose were substituted for alternative sweeteners, and whey powder was used to adjust the nonfat solid content. A mixture design for four ingredients (sucrose, fructose, stevia, and whey powder) was applied. The ice cream composition was correlated with the thermophysical properties (density, viscosity, heat capacity, thermal diffusivity, and thermal conductivity). The optimal formulation, based on the quality parameters (time until first drop falls, melting time, and overrun), was determined. Additionally, the quality parameters were correlated with the thermal properties. The whey powder had a positive effect on the variables measured, improving the quality parameters. The best formulation was 12.1782 % whey powder and 0.0218 % stevia. No significant sensory instead difference was found between this formulation and a commercial ice cream sample. Sucrose was totally replaced with stevia and whey powder in the final product.

 Resumen

Se evaluaron las propiedades termofísicas y los parámetros de calidad de un helado de vainilla en función de la formulación. La sacarosa se sustituyó por edulcorantes alternativos, y se usó lactosuero en polvo para ajustar el contenido sólido lácteos no grasos. Se aplicó un diseño de mezcla para cuatro ingredientes (sacarosa, fructosa, estevia y suero en polvo). La composición del helado se correlacionó con las propiedades termofísicas (densidad, viscosidad, capacidad de calor, difusividad térmica y conductividad térmica). Se determinó la formulación óptima, basada en los parámetros de calidad (tiempo de caída de primera gota, el tiempo de derretimiento y overrun). Además, los parámetros de calidad se correlacionaron con las propiedades térmicas. El lactosuero en polvo evidenció un efecto positivo sobre las variables medidas, mejorando los parámetros de calidad. La mejor formulación fue 12.1782% de suero de leche en polvo y 0.0218% de estevia. No se encontraron diferencias sensoriales significativas entre esta formulación y una muestra comercial de helado. La sacarosa fue reemplazada por la estevia y el suero en polvo en el producto final.

 <bold> <italic>Keywords</italic>:</bold> thermal properties stevia whey mixture design <bold> <italic>Palabras clave</italic>:</bold> propiedades térmicas estevia suero diseño de mezcla 1. Introduction

Ice cream is a pasteurized, homogenized, frozen food product characterized by a complex matrix that contains proteins, fat, air, minerals, additives, and sweeteners [1]. Among the many available sweeteners, the most commonly used in the production of ice cream are sucrose, glucose, lactose, and fructose. These sweeteners are intended to provide energy and a sweet taste to the product and to reduce the freezing point of the mixture; however, they also act as antifreeze agents and increase the nonfat solids content of the mixture [2]. Because of the increasing number of people with metabolic disorders, there is a corresponding increase in the demand for low-calorie foods, which is the motivation for this study of alternatives in ice cream formulation. Stevia is a sweetener that has gained increasing interest recently because of its high sweetening power (SP) [3].

Moreover, consumer demand has forced the industry to offer food with better nutritional and sensory quality through innovative ingredients [4] Whey (a byproduct of cheese manufacturing), whose primary components are lactose and nitrogenous protein compounds [5,6], can be used as a substitute for nonfat solids in the manufacturing of ice cream. Whey improves the sensory characteristics, such as a colder feeling, and the textural characteristics, such as the creaminess and the resistance to melting. Whey proteins also improve flavor absorption because they have a high molecular weight, increase the consistency of the ice cream by promoting the formation of smaller ice crystals, and improve the agglomeration of fat globules [2,7].

In the preparation of ice cream, in addition to composition, the thermal properties are important. These properties are essential to estimate freezing times and, of particular importance, to determine the quality and stability of the ice cream as affected by processing and storage temperatures [8].

Most studies aim to replace fats and sugars that have high caloric values with substitutes that do not add many calories to the product [9,10]. Therefore, new ingredients or changes in formulation, such as the addition of whey and gums, have been proposed to improve sensory quality and increase creaminess and melt stability [6,11,12]. These studies were limited to investigating the effects of formulation on quality. We have not yet found any studies on the relation between the quality parameters and the thermal properties of the mixture.

The aim of this study was to assess the effects of replacing the sucrose and the nonfat solids with sweeteners and whey powder, respectively, on the thermophysical properties and the quality parameters of dairy ice cream.

 2. Materials and methods <italic>2.1. Ingredients</italic>

Homogenized (UHT) whole milk (ALQUERÍA S.A, Colombia), nonfat dry milk powder (COLANTA, Colombia), homogenized and pasteurized milk cream (ALQUERÍA S.A., Colombia), whey powder (TECNAS S.A., Colombia), sucrose (MANUELITA S.A., Colombia), fructose (Laboratorio San Jorge Ltda., Colombia), stevia (Shanghai Tianjia Biochemical Co. Ltda., China), emulsifier-stabilizer (CIMPA S.A., Colombia), and vanilla essence (Royal Alimentos Ltda, Colombia) were the ingredients used for the ice cream formulation. Table 1 shows the formulation and the fat and nonfat solids content for each test sample.

Ice cream formulation

a Depending on the restriction on nonfat solids.

b Depending on the mixture design, sucrose (SP: 1), fructose (SP: 1.33) or stevia (SP: 280) were used

Source: The Authors

<italic>2.2. Ice cream preparation</italic>

One liter of the mixture was prepared for each test. The liquid ingredients were mixed and heated to 50°C. The combined dry ingredients were added and mixed using a stirrer (DLS model, Velp Scientifica, Europe) at 800-1000 rpm. The mixture was pasteurized at 80°C for 30 s, cooled to 4°C and matured at 4ºC for 24 hours. Finally, the mixture was frozen at -22°C for 50 min using an ice cream maker (Princess, Model 282600, Germany) at 1800 rpm and stored at -17°C [13]. The prepared ice cream met the minimum requirements of 26 % nonfat solids (NFS) and 10 % fat solids (FS), according to the Colombian Technical Standard [14].

Mixture design for four components

WP: Whey powder, SU: Sucrose, FR: Fructose, ST: Stevia

Source: The Authors

<italic>2.3. Experimental design and statistical analysis</italic>

An extreme vertices mixture design method for 4 components (whey powder (WP), sucrose (SU), fructose (FR), and stevia (ST)) was used to determine the optimal proportions of the components (g/100g). The experimental design is shown in Table 2.

Equations (1) to (6) show the linear restrictions for each component with respect to the SP and the quantity of NFS that they contribute to the mixture. The response variables in the ice cream mixture after maturation were the thermal properties (specific heat capacity, thermal diffusivity and thermal conductivity), the physical properties (density and consistency coefficient), and the parameters of ice cream quality (overrun, time until the first drop falls and melting time). An analysis of variance (ANOVA) was performed to find the optimal treatment, maximizing the quality parameters. The statistical software Minitab (version 17.1.0, Minitab, Inc.) was used.

<italic>2.3.1. General restrictions</italic>

ICONTEC [14] requires ice cream to have a minimum of 26 % NFS and 10 % FS. With this restriction, and including all ingredients, a system of equations was developed. The formulation with the least sucrose content was chosen from among the possible combinations, shown in Eq. (1) as follows:

Where WP, SU, FR, and ST are the concentrations of the sweeteners (g/100 g).

2.3.2. Restrictions on the NFS

The restrictions on the NFS based on sucrose substitution are described by Eqs. (2) to (6) as follows:

 2.3.3. Restriction on the SP

The SP restriction (Eq. (6)) determines the quantity of sweetener that is added to the mixture, which should not be greater than 12.2 g/100 g.

 <italic>2.4. Thermal properties</italic> 2.4.1. Specific heat capacity (Cp)

The C p (J kg-1 K-1) of the mixture was determined using a differential scanning calorimeter (DSC) (Mettler Toledo). Approximately 20 mg of the sample was added to an aluminum crucible. The following protocol was then performed: the sample was a) cooled from 25°C to -80°C at 10°C/min, b) maintained at -80°C for 10 min to promote icing, and c) heated from -80°C to 80°C at 10°C/min [9]. The C p was determined by the relation between the derivatives of enthalpy (H, mJ) and temperature (T, K) with respect to time (t, s), according to Eq. (7). The C p was calculated as the mean of the values obtained during the freezing stage (the linear section of the thermogram H vs. t).

 2.4.2. Thermal conductivity (k)

The k (W m-1 K-1) of the sample was determined by applying Eq. (8) to the data from the melting curve obtained while measuring the C p , where S is the slope of the linear section of the ice-melting peak, h is the height of the aluminum crucible (1.6mm) and A is the radial area of the crucible (28.27mm2) [15] as follows:

 2.4.3. Thermal diffusivity (𝛂)

The α (mm2 s-1) was calculated using Eq. (9) as follows [16]:

 <italic>2.5. Physical properties</italic> 2.5.1. Density (ρ)

The gravimetric-volumetric method, which relates the mass of the ice cream mixture at 4°C to the occupied volume of the sample, was used to calculate ρ (Eq. 10), where m is the mass of the mixture (kg), and V is the volume occupied by the sample (m3), as follows:

 2.5.2. Consistency coefficient (K)

A rotational viscometer (Brookfield LVDV-III-U), coupled to a cooling recirculating bath at a constant temperature (4°C), was used to determine the rheological behavior of the dairy ice cream sample after maturation. The consistency coefficient (K, Pa·sn) and the flow behavior index (n) were calculated from the shear stress (σ, Pa·s) and the shear rate (γ, s-1), using Eq. (11) [9] as follows:

 2.6. Quality parameters 2.6.1. Overrun (O)

Overrun (%) is defined as the volume change resulting from the aeration, i.e., caused by the capacity of air incorporated into the ice cream mixture, and described by Eq. (12) as follows [17]:

Where V m is the volume of the mixture and V i is the volume of the ice cream.

2.6.2. First drop time (D) and melting time (M)

Approximately 50 g of ice cream was placed in a mesh (56 holes/cm2) and held at room temperature, controlled at 24.7±0.96°C.

Thermophysical properties and quality parameters

C p : Specific heat capacity, α: Thermal diffusivity, k: Thermal conductivity, ρ: Density, K: Consistency coefficient, n: flow behavior index (dimensionless), O: Overrun, D: First drop time, M: Melting time.

Source: The Authors

The first drop time (D) was determined and the weight of the melted ice cream was recorded as a function of time. The melting time (M) was defined as the time required to melt 30 g of the product [18].

 2.6.3. Sensory analysis

Descriptive statistics were applied to assess the difference between the optimal formulation and a commercial ice cream sample. The assessment was performed with 44 untrained evaluators. All samples were served in cylindrical shapes 3 cm high with 2 cm diameters in 30 ml plastic containers with lids [19].

 3. Results and discussion

The results for each response variable, according to the experimental design, are shown in Table 3.

<italic>3.1. Effects of the mixture ingredients on the response variables</italic>

The contour graphs of the proportions of the ingredients WP, SU, and FR, considering the restrictions related to NFS for each experimental response are shown in Figs. 1 to 8. The vertices correspond to each of the pure components in the mixture, i.e., the vertex of WP represents 12.2 % WP and 0 % sweetener). The stevia (ST) was omitted because the mass ratio of ST with respect to the others is very low, and the effect of ST on the response variables is reflected in the ratio of WP used to meet the required NFS level.

 3.1.1. Specific heat capacity

The results for C p are shown in Fig. 1. The value of C p for the mixture decreases as the proportions of SU and FR increase and the amount of WP decreases; values from 1646 to 1937 J/kg K were measured.

Agrawal, et al. [20] calculated values between 2516 J/kg K and 2542 J/kg K from the C p of each ingredient for ice cream with ginger juice. The composition in their study was 26 % NFS and 12 % FS (the C p of fat is 2.09 kJ/kg K and the C p of NFS is 1.93 kJ/kgK). It is possible that the differences between values are from the methods used for the measurements of C p (indirect, depending on the composition vs. direct, using DSC). In these analyses, it is considered that the air incorporated in the ice cream mixture has a negligible effect on the thermal and physical properties [21]. In general, a low C p indicates a short ice cream hardening time [22].

Contour graph for specific heat capacity (J kg<sup>-1</sup> K<sup>-1</sup>) in function of the formulation Source: The Authors

Contour graph for thermal conductivity (W m<sup>-1</sup> K<sup>-1</sup>) in function of the formulation Source: The Authors

3.1.2. Thermal conductivity

The results of the measurements of k are shown in Fig. 2. The measured k for the ice cream mixture varied between 0.237 W/m K and 0.374 W/m K.

Cogné, et al. [21] mentioned that the ice cream mixture can be represented as a material with two phases, a solid phase (ice crystals) dispersed in a continuous fluid phase composed of dry matter and unfrozen water at -20°C. These authors reported k values from 0.466 to 1.51 W/m K in commercial ice cream with 9.3 % FS and 31.83 % NFS. The difference in k is probably from the state of the sample because the k of the liquid is lower than the k of the solid, e.g, k water is 0.58 W/m K and k ice is 2.18 W/mK. Additionally, the high proportion of NFS in the composition of the food increases the thermal conductivity. Sudhir and Ashis [23] reported k of 0.543 to 0.580 W/m K for ice cream with 27.95 % NFS and 13 % FS and noted that the air (k air is 0.02 W/m°C) inside the food matrix has a negative effect on the thermal conductivity; however, this influence is not addressed in current literature. Additionally, the properties of a food depend primarily on its composition. A lower C p implies that less heat must be eliminated from the material, whereas a lower k indicates that the rate at which heat is eliminated is also lower [24].

3.1.3. Thermal diffusivity

The results in Fig. 3 show the thermal diffusivity varying from 0.1259 to 0.1774 mm2/s. These results are consistent with the results obtained by Ben-Yoseph and Hartel [25] in vanilla ice cream with 12 % FS and 27.5 % NFS, who reported values of 0.01667 mm2/s and 0.15 mm2/s at -5°C and -20°C, respectively.

The value of α indicates how quickly heat can be eliminated during ice cream hardening under storage. The value of α is primarily determined by the composition of the product and the quantity of air inside the mixture. This property is not fully understood in ice cream. It is likely that α is also influenced by the overrun, the cell size, the air distribution, the ice quantity and the sizes of the ice crystals. Increasing the quantity and the dispersion of the air cells and the ice crystals should reduce the thermal diffusivity, resulting in slower melting, i.e., a higher melting time [2].

Contour graph for thermal diffusivity (mm<sup>2</sup> s<sup>-1</sup>) in function of the formulation Source: The Authors

Contour graph for density (kg m<sup>-3</sup>) in function of the formulation Source: The Authors

3.1.4. Density

The results in Fig. 4 show the influence of the ingredients on the density (ρ) of the ice cream mixture.

The density changes depending on the composition, i.e., ρ increases as the amount of NFS increases and as the fat content decreases [2]. The density values ranged from 1118 to 1143 kg/m3. As the WP content increased, the density of the mixture decreased. This trend is related to the emulsifying and foaming properties of WP. It is also related to the density of each of the components, i.e., ρ of FR and SU are 1690 kg/m3 and 1590 kg/m3, respectively, whereas ρ of WP is from 0.7 to 0.85 kg/m3 [26]. Warren and Hartel [27] reported density values for an ice cream mixture from 1070 to 1160 kg/m3. Those values are comparable to the density values found in this study.

 3.1.5. Rheological properties

The results for the consistency coefficient (K) of the ice cream mixture are shown in Fig. 5. This variable is a key attribute

Contour graph for consistency coefficient (Pasn) in function of the formulation Source : The Authors

because it influences the body and texture of the final product [28]. Values of K equal to 2.022 Pa·sn and to 11.267 Pa·sn were obtained for formulations containing 12.2 % SU and 12.1782 % WP with 0.0218 % ST, respectively. Other studies have reported similar K for ice cream mixtures [29-31].

Dertli, et al. [29] found consistency coefficients of 0.03 Pa·sn and 24.07 Pa·sn for ice cream mixtures without and with fermentation with Streptococcus thermophilus strains. In the case of the fermented ice cream, the consistency coefficient was higher because of the exopolysaccharides. A higher proportion of WP is positively related to K in the matured dairy ice cream mixture.

The flow behavior index (n) is between 0.419 and 0.71, which indicates that the ice cream mixtures are non-Newtonian fluids that exhibit a pseudoplastic behavior (Table 3). Non-Newtonian behavior is generally associated with intermolecular forces and with the effect of the addition of WP on the protein molecule sizes [32]. The n values obtained are consistent with published studies. Toker, et al. [31] studied the effect of processing temperatures in the range of 5 to 35°C on the K of the ice cream mixture (7.036 to 12.883Pa·sn) and on n (0.328 to 0.360). Also Dogan, et al. [30] analyzed the variations in K (0.277 to 74.6Pa·sn) and in n (0.07 to 0.527) caused by the effect of the xanthan gum concentration in ice cream mixtures.

 3.1.6. Overrun

The results of the overrun measurements are shown in Fig. 6. The overrun was between 48 % (treatment with 12.2 % SU) and 126 % (treatments with 12.1782 % WP - 0.0218 % ST). Muse and Hartel [33] reported values from 40 to 70 % overrun for ice cream stored in a batch freezer. The overrun results were good because of the high quantity of air incorporated into the mixture, i.e., the result was a soft ice cream with small crystals that is more resistant to recrystallization during storage [34].

The overrun increased as the WP was increased, and decreased as the proportion of SU increased. This finding could be related to the viscosity and to the difficulty of trapping air in less consistent mixtures (low K). Goff, et al. [35] indicated that whey proteins possess the ability to stabilize emulsions, reducing the interfacial tension and forming a membrane that allows partial coalescence, and therefore improving the binding between fat globules. In addition, stable foams can be formed when heat and the casein micelles interact, resulting in shear stress. Therefore, the use of WP improves the ice cream properties, giving a mild sweetness, increasing the nutritional value and providing good emulsifying properties, structure and texture. This is consistent with the higher values of overrun that are a direct function of the proportion of WP added to the mixture.

Contour graph for overrun in function of the formulation Source: The Authors

3.1.7. First drop time and melting time

The results for the D and the M of the ice cream samples are shown in Figs. 7 and 8, respectively. These characteristics are particularly important to the consumer. Values of D and M for the treatment with sucrose at 12.2 % were 8.11 min and 19.5 min, respectively, while for the treatment with 12.1782 % WP - 0.0218 % ST, the D and M values were 13.76 min and 33.83 min, respectively.

Bahram-Parvar, et al. [36] reported melting times between 37.40 min and 50.80 min when adding gums to the mixture. However, the range of melting rates is consistent with the study of Karaca, et al. [37], who reported a melting rate of 1.17 to 1.91 g/min for ice creams that contain different fat substitutes, values slightly lower than those found in the present study (1.54 g/min and 3.69 g/min). The melting resistance in samples with high proportions of WP is from the increased viscosity of the mixture. In addition, when the overrun is high, the melting time is also high from the effect of the air bubbles in the mixture [27].

Contour graph for first drop time (min) in function of the formulation Source: The Authors

Contour graph for melting time [min (30 g)<sup>-1</sup>] in function of the formulation Source: The Authors

The behavior of the mixture can be associated with the rheological properties of the frozen phase of the ice cream. Ice crystals melt and the water diffuses in the frozen phase. The formulations of ice cream with high consistency coefficients and high viscosities have high melting resistance, i.e., the ice cream melts slowly [33].

 <italic>3.2. Statistical analysis</italic>

Polynomial models describe the correlation between the responses and the variables, as shown in Table 4. The coefficients of determination (R2 and adjusted R2) indicate that the selected models have a good fit of the data. A linear model was chosen for C p and k, quadratic models were adjusted for α, K, O and D and a cubic model was adjusted for the variables ρ and M. In the regression models, only the sources of variation with p ≤ 0.05 were included. The positive sign of a coefficient in the adjusted model indicates the magnitude of an increase in the variable, while a negative sign shows a decrease.

 <italic>3.3. Optimization of the formulations</italic>

The optimal proportions of each ingredient to satisfy the quality parameters shown in Table 5 were determined. Consumers prefer ice cream with a low melting rate, i.e., both D and the M should be high.

The industry prefers to maximize the overrun. Hence, these three quality parameters were the response variables to maximize, giving to all of them the same weighing. The optimal proportions were determined to be 12.1782 % WP, 0 % SU, 0 % FR, and 0.0218 % ST. This result indicates that sucrose can be completely substituted by sweeteners with low-caloric value and high sweetening power by taking advantage of WP to produce a dairy ice cream with vanilla flavor that has very good quality features for high consumer acceptability. This promising result can be important for people that cannot consume ice cream because of its high sugar content.

Mathematical models for determination thermophysical properties and quality parameters

C p : Specific heat capacity (J kg-1 K-1), α: Thermal diffusivity (mm2 s-1), k: Thermal conductivity (W m-1 K-1), ρ: Density (kg m-3), K: Consistency coefficient (Pa.sn), O: Overrun (%), D: First drop time (min), M: Melting time (min/30g), WP: Whey powder (g/12.2 g sweetener), SU: Sucrose (g/12.2 g sweetener), FR: Fructose (g/12.2 g sweetener) y ST: Stevia (g/12.2 g sweetener).

Source: The Authors

Conditions for optimization

Source: The Authors

Properties calculated based on the optimal formulation

Source: The Authors

A sample of ice cream with the optimized formulation was prepared to evaluate it and validate the models. The properties of the sample are shown in Table 6. The results presented in Table 6 show moderate errors in most cases.

 <italic>3.4. Correlation of the response variables</italic>

A nonlinear regression was used to mathematically describe the relation that exists among the quality parameters as a function of the thermal properties of the dairy ice cream. With these correlations, the ice cream behavior can be predicted, depending on the thermal properties of the ice cream mixture. The mathematic correlations and the coefficients of determination R2 and adjusted R2 are listed in Table 7. The properties of the sample are shown in Table 8. Small differences were found between the observed and the predicted values, verifying the optimization process.

Correlations for determining quality parameters depending on the thermal properties

O: Overrun (%), D: First drop time (min), M: Melting time (min/30g), Cp: Specific heat capacity (J kg-1 K-1), α: Thermal diffusivity (mm2 s-1), k: Thermal conductivity (W m-1 K-1), ρ: Density (kg m-3).

Source: The Authors

Quality parameters calculated based on the thermal properties

Source: The Authors

 <italic>3.5. Sensory analysis</italic>

The results of the sensory evaluation of pair-wise comparisons to differentiate between the optimal experimental sample and the commercial sample showed no significant difference (p ≤ 0.05), with a proportion of the total population of evaluators able to differentiate between both samples p d of 50 % and a type II error β of 0.01 No differences between the optimal formulation and the commercial sample were found by consumers, verifying that the optimized formulation shows similar sensory properties to the commercial product [19, 38].

 4. Conclusions

The formulation of the mixture has influenced the thermophysical properties and quality parameters of a dairy ice cream with vanilla flavor. Whey powder and stevia were the most influencing ingredients. An extreme vertices mixture design was applied to find the formulation with the greatest overrun, the highest first drop time, and the highest melting time to satisfy an industrial need and the consumers’ requirements. The optimal formulation replaced all the sucrose with stevia and added whey powder. This formulation exhibited sensory characteristics that were not significantly different from a commercial sample. Mathematical correlations were developed to predict the parameters of ice cream quality as a function of the thermal properties of the ice cream mixture.

 References [1] Ozdemir, C.,Dagdemir, E.,Ozdemir, S.andSagdic, O., The effects of using alternative sweeteners to sucrose on ice cream quality. J Food Qual. 31(4), pp. 415-428, 2008. DOI: 10.1111/j.1745-4557.2008.00209.x. Ozdemir C. Dagdemir E. Ozdemir S. Sagdic O The effects of using alternative sweeteners to sucrose on ice cream quality J Food Qual 31 4 415 428 2008 10.1111/j.1745-4557.2008.00209.x [2] Goff, H.D. and Hartel, R.W., Ice Cream vol. 7. New York, USA: Springer Science Business Media, 2013. Goff H.D. Hartel R.W. Ice Cream 7 New York, USA: Springer Science Business Media 2013 [3] Soukoulis, C.,Fisk,I.D. and Bohn, T., Ice Cream as a vehicle for incorporating health-promoting ingredients: conceptualization and overview of quality and storage stability. Compr Rev Food Sci Food Safety. 13(4), pp. 627-655, 2014. DOI: 10.1111/1541-4337.12083. Soukoulis C. Fisk I.D. Bohn T. Ice Cream as a vehicle for incorporating health-promoting ingredients: conceptualization and overview of quality and storage stability. Compr Rev Food Sci Food Safety 13 4 627 655 2014 10.1111/1541-4337.12083 [4] Ozdemir, C.,Arslaner,A.,Ozdemir,S.andAllahyari, M., The production of ice cream using stevia as a sweetener. J Food Sci Technol. 52(11), pp. 7545-7548, 2015. DOI: 10.1007/s13197-015-1784-5. Ozdemir C. Arslaner A. Ozdemir S. Allahyari M. The production of ice cream using stevia as a sweetener. J Food Sci Technol. 52 11 7545 7548 2015 10.1007/s13197-015-1784-5. [5] Posada, K.,Terán,D.M.yRamírezNavas, J.S., Empleo de lactosuero y sus componentes en la elaboración de postres y productos de confitería. Aliment Latinoam. 292, pp. 66-75, 2011. Posada K. Terán D.M. RamírezNavas J.S. Empleo de lactosuero y sus componentes en la elaboración de postres y productos de confitería. Aliment Latinoam. 292 66 75 2011 [6] Rodríguez, T. and M’Boumba, A., Utilización del suero del queso en helados. Cienc Tecnol Alim. 21(2), pp. 9-12, 2011. Rodríguez T. M’Boumba A. Utilización del suero del queso en helados Cienc Tecnol Alim 21 2 9 12 2011 [7] Young, S., Whey products in ice cream and frozen dairy desserts - Applications Monograph -Houston, Texas, USA, U.S. Dairy Export Council, 2007. Young S. Whey products in ice cream and frozen dairy desserts - Applications Monograph -Houston, Texas, USA, U.S. Dairy Export Council 2007 [8] Renaud, T., Briery, P., Andrieu, J. and Laurent, M., Thermal properties of model foods in the frozen state. J Food Eng. 15(2), pp. 83-97, 1992. DOI: 10.1016/0260-8774(92)90027-4. Renaud T. Briery P. Andrieu J. Laurent M. Thermal properties of model foods in the frozen state. J Food Eng 15 2 83 97 1992 10.1016/0260-8774(92)90027-4. [9] Soukoulis, C.,Rontogianni, E.andTzia, C., Contribution of thermal, rheological and physical measurements to the determination of sensorially perceived quality of ice cream containing bulk sweeteners. J Food Eng. 100(4), pp. 634-641, 2010. DOI: 10.1016/j.jfoodeng.2010.05.012. Soukoulis C. Rontogianni E. Tzia C. Contribution of thermal, rheological and physical measurements to the determination of sensorially perceived quality of ice cream containing bulk sweeteners J Food Eng. 100 4 634 641 2010 10.1016/j.jfoodeng.2010.05.012. [10] Whelan, A.P.,Vega, C.,Kerry, J.P.andGoff, H.D., Physicochemical and sensory optimisation of a low glycemic index ice cream formulation. Int J Food Sci Technol. 43(9), pp. 1520-1527, 2008. DOI: 10.1111/j.1365-2621.2007.01502.x. Whelan A.P. Vega C. Kerry J.P. Goff H.D. Physicochemical and sensory optimisation of a low glycemic index ice cream formulation. Int J Food Sci Technol. 43 9 1520 1527 2008 10.1111/j.1365-2621.2007.01502.x. [11] Bazmi, A., Launay, B., Cuvelier, G. and Relkin, P., Impact of crystalline milk fat on rheological properties of Ice cream mix emulsions during aging time at 4°C. J Texture Stud. 39(4), pp. 309-325, 2008. DOI: 10.1111/j.1745-4603.2008.00145.x. Bazmi A. Launay B. Cuvelier G. Relkin P. Impact of crystalline milk fat on rheological properties of Ice cream mix emulsions during aging time at 4°C. J Texture Stud. 39 4 309 325 2008 10.1111/j.1745-4603.2008.00145.x [12] Klesment, T., Stekolštšikova, J. and Laos, K., Influence of guar gum/furcellaran and guar gum/carrageenan stabilizer systems on the rheological and sensorial properties of ice cream during storage. P Estonian Acad Sci., 63(2), pp. 193-198, 2014. Klesment T. Stekolštšikova J. Laos K. Influence of guar gum/furcellaran and guar gum/carrageenan stabilizer systems on the rheological and sensorial properties of ice cream during storage. P Estonian Acad Sci. 63 2 193 198 2014 [13] Villalva, F., Cravero, A. y Ramón, A., Elaboración y composición química de un helado sabor durazno, con características prebióticas. Helad. Panad. Latinoam. 224, pp. 58-61, 2013. Villalva F. Cravero A. Ramón A. Elaboración y composición química de un helado sabor durazno, con características prebióticas. Helad. Panad. Latinoam. 224 58 61 2013 [14] ICONTEC, NTC1239, Helados y mezclas para helados, vol. 1239, ed. Bogota, Colombia: Instituto Colombiano de Normas Tecnicas y Certificacion, 2002, 21 P. ICONTEC NTC1239, Helados y mezclas para helados 1239 Colombia Instituto Colombiano de Normas Tecnicas y Certificacion 2002 21 [15] Riesen, R., Simple determination of the thermal conductivity of polymers by DSC. Thermal Analysis Aplication. UserCom. 22(2), pp. 19-23, 2005. Riesen R. Simple determination of the thermal conductivity of polymers by DSC. Thermal Analysis Aplication 22 2 19 23 2005 [16] Ibarz, A. and Barbosa- Cánovas, G., Operaciones unitarias en ingeniería de alimentos. Madrid: Mundi-Prensa, 2005. Ibarz A. Barbosa- Cánovas G. Operaciones unitarias en ingeniería de alimentos Madrid Mundi-Prensa 2005 [17] Drewett, E.M. and Hartel, R.W., Ice crystallization in a scraped surface freezer. J Food Eng. 78(3), pp. 1060-1066, 2007. DOI: 10.1016/j.jfoodeng.2005.12.018. Drewett E.M. Hartel R.W. Ice crystallization in a scraped surface freezer. J Food Eng. 78 3 1060 1066 2007 10.1016/j.jfoodeng.2005.12.018. [18] Granger, C., Leger, A., Barey, P., Langendorff, V. and Cansell, M., Influence of formulation on the structural networks in ice cream. Int Dairy J. 15(3), pp. 255-262, 2005. DOI: 10.1016/j.idairyj.2004.07.009. Granger C. Leger A. Barey P. Langendorff V. Cansell M. Influence of formulation on the structural networks in ice cream. Int Dairy J. 15 3 255 262 2005 10.1016/j.idairyj.2004.07.009. [19] ICONTEC, NTC2680, Análisis sensorial. Metodología. Prueba de comparación pareada, vol. 2680, ed. Bogota, Colombia: Instituto Colombiano de Normas Tecnicas y Certificacion , 2013, 32 P. ICONTEC NTC2680, Análisis sensorial. Metodología. Prueba de comparación pareada 2680 Bogota, Colombia Instituto Colombiano de Normas Tecnicas y Certificacion 2013 32 32 [20] Agrawal, A., Karkhele, P., Sandey, K., Sahu, C. and Sinha, G., Effect of incorporation of ginger juice in various rates on the freezing and thermal properties of ice cream. Asian J. Dairy Food Res. 34(2), pp. 92-97, 2015. Agrawal A. Karkhele P. Sandey K. Sahu C. Sinha G. Effect of incorporation of ginger juice in various rates on the freezing and thermal properties of ice cream. Asian J. Dairy Food Res. 34 2 92 97 2015 [21] Cogné, C., Andrieu, J., Laurent, P., Besson, A. and Nocquet, J., Experimental data and modelling of thermal properties of ice creams. J Food Eng. 58(4), pp. 331-341, 2003. DOI: 10.1016/s0260-8774(02)00396-5. Cogné C. Andrieu J. Laurent P. Besson A. Nocquet J. Experimental data and modelling of thermal properties of ice creams. J Food Eng. 58 4 331 341 2003 10.1016/s0260-8774(02)00396-5 [22] Clarke, C., The science of ice cream: Royal Society of Chemistry, 2012. Clarke C. The science of ice cream: Royal Society of Chemistry 2012 [23] Sudhir, K.S. and Ashis, K.D., Thermal properties of frozen peas, clams and ice cream. Can I Food Sci Tech J., 17(4), pp. 242-246, 1984. DOI: 10.1016/s0315-5463(84)72565-x. Sudhir K.S. Ashis K.D. Thermal properties of frozen peas, clams and ice cream. Can I Food Sci Tech J. 17 4 242 246 1984 10.1016/s0315-5463(84)72565-x [24] Clarke, C., The science of ice cream, 2 ed. Cambridge The Royal Society of Chemistry:, 2004. Clarke C. The science of ice cream 2 Cambridge The Royal Society of Chemistry 2004 [25] Ben-Yoseph, E. and Hartel, R.W., Computer simulation of ice recrystallization in ice cream during storage. J Food Eng. 38(3), pp. 309-329, 1998. DOI: 10.1016/s0260-8774(98)00116-2. Ben-Yoseph E. Hartel R.W. Computer simulation of ice recrystallization in ice cream during storage J Food Eng 38 3 309 329 1998 10.1016/s0260-8774(98)00116-2. [26] Chegini, G.R. and Taheri, M., Whey powder: process technology and physical properties: A review. Middle East J Sci Res. 13(10), pp. 1377-1387, 2013. DOI: 10.5829/idosi.mejsr.2013.13.10.1239. Chegini G.R. Taheri M. Whey powder: process technology and physical properties: A review. Middle East J Sci Res. 13 10 1377 1387 2013 10.5829/idosi.mejsr.2013.13.10.1239. [27] Warren, M.M. and Hartel, R.W., Structural, compositional, and sensorial properties of United States commercial ice cream products. J Food Sci. 79(10), pp. E2005-E2013, 2014. DOI: 10.1111/1750-3841.12592. Warren M.M. Hartel R.W. Structural, compositional, and sensorial properties of United States commercial ice cream products. J Food Sci. 79 10 E2005 E2013 2014 10.1111/1750-3841.12592. [28] Stanley, D.W., Goff, H.D. and Smith, A.K., Texture-structure relationships in foamed dairy emulsions. Food Res Int. 29(1), pp. 1-13, 1996. DOI: 10.1016/0963-9969(95)00063-1. Stanley D.W. Goff H.D. Smith A.K. emulsions. Texture-structure Food Res Int. 29 1 1 13 1996 10.1016/0963-9969(95)00063-1 [29] Dertli, E., Toker, O.S., Durak, M.Z., Yilmaz, M.T., Tatlısu, N.B., Sagdic, O. and Cankurt, H., Development of a fermented ice-cream as influenced by in situ exopolysaccharide production: rheological, molecular, microstructural and sensory characterization. Carbohyd Polym. 136, pp. 427-440, 2016. DOI: 10.1016/j.carbpol.2015.08.047. Dertli E. Toker O.S. Durak M.Z. Yilmaz M.T. Tatlısu N.B. Sagdic O. Cankurt H. Development of a fermented ice-cream as influenced by in situ exopolysaccharide production: rheological, molecular, microstructural and sensory characterization. Carbohyd Polym. 136 427 440 2016 10.1016/j.carbpol.2015.08.047 [30] Dogan, M., Kayacier, A., Toker, Ö.S., Yilmaz, M.T. and Karaman, S., Steady, dynamic, creep, and recovery analysis of ice cream mixes added with different concentrations of xanthan gum. Food Bioprocess Technol. 6(6), pp. 1420-1433, 2012. DOI: 10.1007/s11947-012-0872-z. Dogan M. Kayacier A. Toker Ö.S. Yilmaz M.T. Karaman S. Steady, dynamic, creep, and recovery analysis of ice cream mixes added with different concentrations of xanthan gum. Food Bioprocess Technol. 6 6 1420 1433 2012 10.1007/s11947-012-0872-z. [31] Toker, O.S., Karaman, S., Yuksel, F., Dogan, M., Kayacier, A. and Yilmaz, M.T., Temperature dependency of steady, dynamic, and creep-recovery rheological properties of ice cream mix. Food Bioprocess Technol. 6(11), pp. 2974-2985, 2012. DOI: 10.1007/s11947-012-1005-4. Toker O.S. Karaman S. Yuksel F. Dogan M. Kayacier A. Yilmaz M.T. Temperature dependency of steady, dynamic, and creep-recovery rheological properties of ice cream mix. Food Bioprocess Technol. 6 11 2974 2985 2012 10.1007/s11947-012-1005-4 [32] Morison, K.R. and Mackay, F.M., Viscosity of lactose and whey protein solutions. Int J Food Properties. 4(3), pp. 441-454, 2001. DOI: 10.1081/jfp-100108647. Morison K.R. Mackay F.M. Viscosity of lactose and whey protein solutions. Int J Food Properties. 4 3 441 454 2001 10.1081/jfp-100108647 [33] Muse, M.R. and Hartel, R.W., Ice cream structural elements that affect melting rate and hardness. J Dairy Sci. 87(1), pp. 1-10, 2004. DOI: 10.3168/jds.S0022-0302(04)73135-5. Muse M.R. Hartel R.W. Ice cream structural elements that affect melting rate and hardness. J Dairy Sci. 87 1 1 10 2004 10.3168/jds.S0022-0302(04)73135-5 [34] Bahram-Parvar, M., A review of modern instrumental techniques for measurements of ice cream characteristics. Food Chem. 188, pp. 625-631, 2015. DOI: 10.1016/j.foodchem.2015.05.017. Bahram-Parvar M. A review of modern instrumental techniques for measurements of ice cream characteristics. Food Chem. 188 625 631 2015 10.1016/j.foodchem.2015.05.017 [35] Goff, H.D., Kinsella, J.E. and Jordan, W.K., Influence of various milk protein isolates on ice cream emulsion stability. J Dairy Sci. 72(2), pp. 385-397, 1989. DOI: 10.3168/jds.S0022-0302(89)79120-7. Goff H.D. Kinsella J.E. Jordan W.K. Influence of various milk protein isolates on ice cream emulsion stability. J Dairy Sci. 72 2 385 397 1989 10.3168/jds.S0022-0302(89)79120-7 [36] Bahram-Parvar, M., Tehrani, M.M., Razavi, S.M.A. and Koocheki, A., Application of simplex-centroid mixture design to optimize stabilizer combinations for ice cream manufacture. J Food Sci Technol. 52(3), pp. 1480-1488, 2015. DOI: 10.1007/s13197-013-1133-5. Bahram-Parvar M. Tehrani M.M. Razavi S.M.A. Koocheki A. Application of simplex-centroid mixture design to optimize stabilizer combinations for ice cream manufacture. J Food Sci Technol. 52 3 1480 1488 2015 10.1007/s13197-013-1133-5 [37] Karaca, O.B., GÜVen, M., Yasar, K., Kaya, S. and Kahyaoglu, T., The functional, rheological and sensory characteristics of ice creams with various fat replacers. Int J Dairy Technol. 62(1), pp. 93-99, 2009. DOI: 10.1111/j.1471-0307.2008.00456.x. Karaca O.B. GÜVen M. Yasar K. Kaya S. Kahyaoglu T. The functional, rheological and sensory characteristics of ice creams with various fat replacers. Int J Dairy Technol. 62 1 93 99 2009 10.1111/j.1471-0307.2008.00456.x [38] Ramírez-Navas, J.S., Análisis sensorial: pruebas orientadas al consumidor. Rev. RECITEIA. 12(1), pp. 83-102, 2012. Ramírez-Navas J.S. Análisis sensorial: pruebas orientadas al consumidor. Rev. RECITEIA. 12 1 83 102 2012

:Zambrano-Mayorga, L.F., Ramírez-Navas, J.S. and Ochoa-Martínez, C.I., Influence of the formulation on the thermophysical properties and the quality parameters of dairy ice cream.. DYNA, 86(208), pp. 117-125, January - March, 2019

received his BSc. Eng in Food Engineering in 2013 from the Universidad Técnica de Ambato, Ambato, Ecuador, and his MSc. degree in Food Engineering in 2016 from the Universidad del Valle, Cali, Colombia. Currently, he is an Instructor in the Servicio Nacional de Aprendizaje SENA (Buga, Colombia). His research interests include simulation, modeling and food technology process. ORCID: 0000-0001-8934-6021

received his BSc. Chemical Engineering in 2003 from the Universidad Central del Ecuador, Quito, Ecuador, and his PhD in Engineering in 2012 from Universidad del Valle, Cali, Colombia. He has worked in programs and projects in the food, chemistry and nutrition areas since 2012 at the Universidad del Valle, Universidad Santiago de Cali and Corporación Biotec. He is currently an assistant professor in the Food Engineering School, Engineering Faculty, Universidad del Valle, and professor in the Basic Science Faculty, Universidad Santiago de Cali. His research interests include: physical properties of food, process engineering, dairy science, natural ingredients and functional food. He has several publications in scientific journals. ORCID: 0000-0002-6731-2784

received her BSc. Chemical Engineering in 1989, her MSc. in Chemical Engineering in 2001, and her PhD in Engineering in 2006. She has worked in programs and projects in the food area since 2006 at the Universidad del Valle. She is currently a full professor in the Food Engineering School, Facultad de Ingeniería, Universidad del Valle, Cali, Colombia. Her research interests include: modeling, simulation and drying; process engineering. She has several publications in scientific journals and is currently coordinator of the post-graduate program at the Food Engineering School. ORCID:0000-0002-2666-1726