Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1226-4768(Print)
ISSN : 2288-1247(Online)
Food Engineering Progress Vol.21 No.4 pp.303-311
DOI : https://doi.org/10.13050/foodengprog.2017.21.4.303

Kinetic Data for Texture Changes of Foods During Thermal Processing

Seung Hwan Lee*
School of Food Science and Industry, Kyungil University
Corresponding author: Seung Hwan Lee School of Food Science and Industry, Kyungil University 50 Gamasilgil, Hayangeup, Gyeongsan, Gyeongbuk, 38428, Korea +82-53-600-5745; +82-53-600-5759e.seunghwan@gmail.com
25/10/2017 05/11/2017 08/11/2017

Abstract

To automate cooking processes, quantitative descriptions are needed on how quality parameters, such as texture change during heating. Understanding mechanical property changes in foods during thermal treatment due to changes in chemical composition or physical structure is important in the context of engineering models and in precise control of quality in general. Texture degradation of food materials has been studied widely and softening kinetic parameters have been reported in many studies. For a better understanding of kinetic parameters, applied kinetic models were investigated, then rate constants at 100°C and activation energy from previous kinetic studies were compared. The food materials are hardly classified into similar softening kinetics. The range of parameters is wide regardless of food types due to the complexity of food material, different testing methods, sample size, and geometry. Kinetic parameters are essential for optimal process design. For broad and reliable applications, kinetic parameters should be generated by a more consistent manner so that those of foods could be compared or grouped.


초록


    Introduction

    Texture is one of the most important characteristics of foods and quantifying changes in texture is quite useful to product and process and equipment design (Lund, 1982). For producers, postharvest handlers, processors, marketers and consumer, knowledge of texture properties is crucial because texture measurement is the most common techniques and procedures (Chen & Opara, 2013).

    Understanding of pattern and mechanism of the texture changes is necessary to control, improve and predict the quality of heat processed foods. Many researchers have examined thermal softening of food materials and applied chemical reaction kinetic models to quantify a quality attribute. Foods are unstable in the thermodymic sense and tend to be stable kinetically; therefore, the models describing quality changes are kinetic models (Van Boekel, 2008). Real foods are complex biological system. Various reactions such as chemical oxidation, microbial growth, biochemical denaturation and physical changes occur in series or in parallel during thermal processing. Final food quality is decided by interacting and complex reactions rather than a single reaction (Ling et al., 2015). Due to the complexity, modeled foods are normally regarded as a simple system.

    To predict quality change during a simple cooking process, heat transfer and texture kinetic models should be considered simultaneously, and parameters for the both are required to simulate time-temperature history inside foods (Lee et al., 2007). It might be necessary to obtain objective and compatible texture kinetic parameters for precise quality prediction.

    For the better understanding of texture kinetic parameters, published texture softening data, mainly vegetables, were investigated. Temperature of reaction rate constants was set to 100oC for reasonable comparison. Recently, Ling et al. (2015) reviewed kinetics of food quality changes during thermal processing and presented published kinetic parameters of different food product for color, nutrient content, and texture. However, they showed rate constants at different temperature condition.

    This review describes 1) quantifying textural quality, 2) kinetic models on thermal softening of foods, and 3) previous studies on kinetic parameters for future texture kinetic studies. The better understanding of thermal softening kinetic parameters would lead to better process control.

    Quantifying Textural Quality

    The textural quality of foods changes as a result of physicochemical changes during thermal processing (Rao & Lund, 1986) and the quality is perceived by various sensory inputs, interactions between foods and consumers before and during consumption; so, sensory evaluation is a method for expressing texture. However, sensory analysis is difficult to be conducted during each processing step (Jack et al., 1995). Therefore, mechanical properties of foods by instrumental measurements have been used as indices of texture for predicting effect of processing. Wide range of texture measurement instrument from simple hand-held device to instron machine and texture analyzer are adopted for objective tests (Chen & Opara, 2013).

    Human perception to foods is based on several interacting factors while instrumental measurements measure one factor (Sawyer, 1971); therefore, instrumental measurements could not reflect exact sensory perception by consumers. However, both measurement methods are associated with a certain property. For softening, instrumental and sensory results are reported (Rao & Lund, 1986) to be linearly related in studies on green beans, corns, peas, and asparaguses (Hayakawa et al., 1977), black beans (Silva et al., 1981b), and shrimps (Ma et al., 1983). In addition, Collison et al. (1980) showed significant correlation between rupture load and taste panel score of potatoes during cooking. Böhler et al. (1986, 1987) examined the quality of cooked potatoes using both sensory and instrumental methods and coupled texture properties with appropriate instrumental methods.

    Instrumental measurements are practical methods to quantify, collect, and generalize texture changes of foods to predict effects of processing conditions even though it is impossible to thoroughly replace sensory measurements with instrumental measurements.

    Kinetic Models on Thermal Softening of Foods

    First order kinetics

    Most published texture studies express thermal softening as simple first order kinetic models. The first order kinetic model is simple and useful approach to express the degree of softening of foods (Rao & Lund, 1986). First order kinetic model for a texture property, P, is given by (equation 1) and the temperature dependence of softening rate, k, is described by the Arrhenius relationship (equation 2). Texture degradation of vegetables follows a first order kinetic model and the kinetic parameters are valuable for optimizing product quality (Rizvi & Tong, 1997). Even though the first order kinetic model can not express all of the data for extended thermal process times, it is still useful if the desired final texture value is within the range which the simple model describes.(1)(2)

    d P d t = k P
    (1)

    k = k o exp ( E a R T )
    (2)

    Dual mechanism first order kinetics

    Huang & Bourne (1983) studied kinetics of thermal softening of vegetables for extended process times and proposed a dual mechanism first order kinetic model. They examined the rate of softening in several canned vegetables; beans, peas, beets and carrots during retort process at 100-121oC during which the softening was consistent with two simultaneous first order kinetic mechanisms. They hypothesized that there were two substrates ‘a’ and ‘b’ contributing to the firmness A and B and the sum of two was the total texture property at any time (equation 3). Both degradation of A and B could be expressed by first order kinetics. The rate constants for mechanism 1 (equation 4) acting on substrate ‘a’ is greater than those for mechanism 2 (equation 5) acting on substrate ‘b’, and temperature dependent rate constants in both mechanisms are expressed by Arrhenius equation (equation 2). A few assumptions have been made to explain the dual mechanisms thermal softening although there was no complete explanation for the mechanisms of two substrates. Huang & Bourne (1983) reported that the first rapid decrease is attributed to pectic changes in the interlamella layer but the mechanism of the following slow decrease is unknown. Alvarez et al. (2001) observed that the thermal softening kinetics of potato tissue is composed of two simultaneous pseudo first order reactions at 70 and 80oC while it is one pseudo first order kinetic mechanism at 50, 90 and 100oC, and presumes that the fast decreasing in the beginning of cooking time at 70 and 80oC is due to both starch gelatinization and slight changes of pectic substances, and the following slow decreasing is mainly affected by the pectic changes.

    The dual mechanism first order kinetic model describes well the texture softening of vegetables but it needs further explanation for mechanisms of two substrates and abnormal kinetic parameters such as negative activation energy for second mechanism of green beans (Huang & Bourne, 1983).(3)(4)(5)

    Total property P = A + B
    (3)
    3

    Mechanism 1 d A d t = k A A
    (4)

    Mechanism 2 d B d t = k B B
    (5)

    The two simultaneous first order kinetic models by Huang & Bourne (1983) need rigorous mathematical treatment to obtain behavior of mechanism 1 and 2, which are effective whole cooking time. Feng et al. (1989) and Nourian & Ramaswamy (2003) showed two separate first order kinetic models for fast and slow degradation, respectively. It is easier to obtain kinetic parameters of two separate models than two simultaneous models. However, there is no clue to define the time when the first mechanism ends and the second, slower mechanism begins during cooking at different temperatures.

    Fractional conversion technique

    Foods lose their firmness completely after prolonged cooking in the simple first order kinetic and dual mechanism model but foods have nonzero maximum retainable texture property even after prolonged heating. Rizvi & Tong (1997) reexamined published texture softening data described by dual mechanism model and confirmed that the kinetics of softening of vegetables could be expressed well by a simple first order kinetic model when nonzero equilibrium texture property was considered. The fractional conversion factor, f, is defined as the ratio what has reacted in specified time to what must be reacted for completion (equation 6), and the change of (1-f) follows a simple first order reaction kinetics (equation 7). The rate constants in fractional conversion technique follow Arrhenius relationship (equation 2).

    They showed that fractional conversion technique provided more reasonable and reliable texture kinetic parameters than the dual mechanism model. The simple kinetic model using a fractional conversion technique explained well even the thermal softening data of green beans having the negative activation energy (Huang & Bourne, 1983) that the materials softened primarily during the early stage of the reaction and slowly reached to a nonzero equilibrium texture property at a higher temperature.(6)(7)

    f = P o P t P o P
    (6)

    d ( 1 f ) d t = k ( 1 f )
    (7)

    Previous Studies on Kinetic Parameters

    Thermal softening of foods could be quantified and generalized using kinetic models, and the model parameters such as rate constant and activation energy make possible to predict and control product quality. Ling et al. (2015), Rao & Lund (1986) and Villota & Hawkes (1986) reviewed softening kinetics of foods and summarized published kinetic parameters. Available kinetic parameters on softening of foods from these and other studies are summarized in Table 1 where the reaction rate constants at 100oC are calculated from the activation energies which are given or derived from the experiment results. The calculated rate constants do not always fit to experiment results well but the confidence interval issue is beyond the scope of this study. In most of the previous kinetic studies, vegetables were selected as a subject, deformation force or stress was measured as an index texture parameter, and the changes of texture properties were expressed by simple first order kinetic model. Unusually, Mansfield (1974) and Silva et al. (1981a) used the time to reach a certain texture property as an index texture parameter.

    Lund (1982) reviewed that activation energies of muscle tissues were over 209 kJ/mol and those of plant tissues were 63-167 kJ/mol. Salmon (Kong et al., 2007), beef (Bertola et al., 1994) and canned shrimp (Ma et al., 1983) showed toughening and softening during cooking; so it is hard to compare with kinetics parameters from plant foods. In order to compare muscle foods with plant foods, more kinetic studies are needed.

    Previous studies show wide ranges of kinetic parameters even among the same material. Texture indices are affected by not only testing methods such as instrument types and deformation methods but also sample geometry and size; therefore, it is hard to explain clearly the reason for the variation. Remarkably, the rate constants are similar in each study regardless of testing parameters though magnitude of the values is quite different among the studies. After measuring different texture properties of potatoes, both Alvarez et al. (2001) and Nourian & Ramaswamy (2003) obtained narrow ranges of rate constants, 0.059-0.098 and 0.127-0.169 /min, respectively. Harada et al. (1985) tested maximum shear force of three varieties of potatoes and found 0.279-0.396 /min of rate constants. For carrots, Verlinden et al. (1996) showed 0.161- 0.660 /min of rate constants after measuring different texture properties whereas Paulus & Saguy (1980) obtained 0.041- 0.069 /min from three different sizes and varieties of carrots.

    Stoneham et al. (2000) pointed out that the variance of range of rate constants was due to not only testing parameters but also kinetic data reduction procedure on softening kinetics; so, insisted that the inconsistency of kinetic parameters could be overcome by using the fractional conversion technique because this data reduction procedure normalizes experimental parameters on the basis of total texture changes with nonzero final property. They showed that when the texture degradation of potatoes was measured at 80 °C in different instrument types, cross-head speeds, sample sizes, and texture parameters, the rate constants for in the fractional conversion technique were independent of experimental conditions even though the values of texture properties were affected significantly by the testing parameters. There have been only few kinetic studies providing both rate constants and activation energy for texture degradation using the fractional conversion technique. Solomon & Jindal (2003b) used the fractional conversion technique to describe texture changes of potatoes in various testing methods. In order to validate whether the change of conversion factor is certainly independent of testing parameters, further studies are required. Furthermore, the effects of testing protocols and temperatures on initial and equilibrium texture properties should be investigated.

    Calculated rate constants at 100 °C from Loh & Breene (1981), Feng et al. (1989), and Verlinden et al. (1995) were higher than others. They measured very beginning of cooking time and obtained the rate constants, which are only effective for 6 minutes, 15 seconds, and 5 minutes, respectively. Fast loss of turgor pressure at the initial cooking time might be one of the main reasons of the higher rate constants. Later, Verlinden et al. (1996) and Verlinden & De Baerdemaeker (1997) proposed compartment models to describe texture changes of carrots in which the mechanical properties were expressed by sum of components and the mechanical property contributed by turgor pressure, one of the components, was treated as a step function because it disappeared at the beginning of cooking. Kozempel (1988) suggested two separate first order kinetics for potatoes after measuring softening of potatoes at 74-85.5 °C and explained that the initial faster softening was caused by starch gelatinization. The predicted rate constant for the fast softening at 100 °C is very high as 0.781 /min but it is unknown how long the fast degradation lasts and how much the gelatinization process affects on the softening of potatoes at 100 °C. Alvarez et al. (2001) observed similar biphasic degradation of potatoes only at 70 and 80 °C. They concluded that the contribution of gelatinization to thermal softening was limited even though the initial fast softening mechanism was influenced by the gelatinization process, and the softening was mainly affected by solubilization of pectic material in cell wall and middle lamella after starch was gelatinized.

    Conclusions

    Kinetic parameters strongly depend on measurement methods and cooking conditions, and influence of physicochemical changes on texture is quite complex. Rao & Lund (1986) emphasized that it is impossible to quantify all the changes of foods on account of thermal processing, and thermal lag effect must be considered in studies on softening because it depends on both size of the container and investigation methods.

    The better prediction model could be obtained from the better understanding of kinetic parameters, and more reliable process control would be possible as kinetic parameters are generated in more consistent manner. Furthermore, the sensory perception during consumption has not been related to an objective engineering parameter. The relationship between sensory perception by large deformation during chewing and engineering parameters should be defined in future work.

    List of Symbols

    • A: Texture property due to substrate a

    • B: Texture property due to substrate b

    • C: Compartment model

    • D1: Dual mechanism model; two simultaneous first order kinetics

    • D2: Dual mechanism model; two separate first order kinetics Ea: Activation energy (kJ/mol)

    • Ea: Activation energy (kJ/mol)

    • EaA: Activation energy for texture property A (kJ/mol)

    • EaB: Activation energy for texture property B (kJ/mol)

    • Ea1: Activation energy for fast process (kJ/mol)

    • Ea2: Activation energy for slow process (kJ/mol)

    • f : The ratio what has reacted in specified time to what must be reacted for completion

    • F: First order kinetic model using the fractional conversion technique

    • k: Reaction rate constant at temperature T (/min)

    • kA: Reaction rate constant for texture property A (/min)

    • kB: Reaction rate constant for texture property B (/min)

    • ko: Reaction rate constant at T= (/min)

    • k1: Reaction rate constant for fast process (/min)

    • k2: Reaction rate constant for slow process (/min)

    • Pt: Texture property at time t

    • Po: Texture property at time = 0

    • P: Maximum retainable texture property at time = ∞

    • R: Gas constant (8.314 J/mol·K)

    • S: Simple first order kinetic model

    • T: Absolute temperature (K)

    Figure

    Table

    Published texture kinetic parameters of foods
    S: Simple first order kinetic model
    D1: Dual mechanism model; two simultaneous first order kinetics
    D2: Dual mechanism model; two separate first order kinetics
    C: Compartment model
    F: First order kinetic model using the fractional conversion technique

    Reference

    1. Alvarez M.D. , Canet W. , Tortosa M.E. (2001) Kinetics of thermal softening of potato tissue (cv. Monalisa) by water heating. , Eur. Food Res. Technol., Vol.212 ; pp.588-596
    2. Bertola N.C. , Bevilacqua A.E. , Zaritzky N.E. (1994) Heattreatment effect on texture changes and thermal-denaturation of proteins in beef muscle. , J. Food Process. Preserv., Vol.18 ; pp.31-46
    3. Böhler G. , Escher F. , Solms J. (1986) Evaluation of cooking quality of potatoes using sensory and instrumental methods. 1. Sensory evaluation. , Lebensm. Wiss. Technol., Vol.19 ; pp.338-343
    4. Böhler G. , Escher F. , Solms J. (1987) Evaluation of cooking quality of potatoes using sensory and instrumental methods. 2. Instrumental evaluation. , Lebensm. Wiss. Technol., Vol.20 ; pp.207-216
    5. Chen L. , Opara U.L. (2013) Texture measurement approaches in fresh and processed foods - a review. , Food Res. Int., Vol.51 ; pp.823-835
    6. Collison R. , Johnson K. , Okikiolu O.O. , West A. (1980) Subjective and objective assessments of the degree of cooking of potatoes heated by different methods. , J. Food Technol., Vol.15 ; pp.1-8
    7. De Roeck A. , Mols J. , Duvetter T. , Van Loey A. , Hendrickx M. (2010) Carrot texture degradation kinetics and pectin changes during thermal versus high-pressure/high-temperature processing: a comparative study. , Food Chem., Vol.120 ; pp.1104-1112
    8. Feng B. , Cuvelier G. , Brouard F. (1989) Softening of ground beets during heat-treatment. , Sci. Aliments, Vol.9 ; pp.77-88
    9. Harada T. , Tirtohusodo H. , Paulus K. (1985) Influence of temperature and time on cooking kinetics of potatoes. , J. Food Sci., Vol.50 ; pp.459-462
    10. Hayakawa K. , Timbers G.E. , Stier E.F. (1977) Influence of heattreatment on quality of vegetables - organoleptic quality. , J. Food Sci., Vol.42 ; pp.1286-1289
    11. Huang Y.T. , Bourne M.C. (1983) Kinetics of thermal softening of vegetables. , J. Texture Stud., Vol.14 ; pp.1-9
    12. Jack F.R. , Paterson A. , Piggott J.R. (1995) Perceived texture - direct and indirect methods for use in product development. , Int. J. Food Sci. Technol., Vol.30 ; pp.1-12
    13. Jaiswal A.K. , Gupta S. , Abu-Ghannam N. (2012) Kinetic evaluation of colour, texture, polyphenols and antioxidant capacity of Irish York cabbage after blanching treatment. , Food Chem., Vol.131 ; pp.63-72
    14. Ko W.C. , Liu W.C. , Tsang Y.T. , Hsieh C.W. (2007) Kinetics of winter mushrooms (Flammulina velutipes) microstructure and quality changes during thermal processing. , J. Food Eng., Vol.81 ; pp.587-598
    15. Kong F.B. , Tang J.M. , Rasco B. , Crapo C. (2007) Kinetics of salmon quality changes during thermal processing. , J. Food Eng., Vol.83 ; pp.510-520
    16. Kozempel M.F. (1988) Modeling the kinetics of cooking and precooking potatoes. , J. Food Sci., Vol.53 ; pp.753-755
    17. Lau M.H. , Tang J. , Swanson B.G. (2000) Kinetics of textural and color changes in green asparagus during thermal treatments. , J. Food Eng., Vol.45 ; pp.231-236
    18. Lee S.H. , Rao M.A. , Datta A.K. (2007) How does cooking time scale with size? A numerical modeling approach. , J. Food Sci., Vol.72 ; pp.E1-E10
    19. Ling B. , Tang J. , Kong F. , Mitcham E.J. , Wang S. (2015) Kinetics of food quality changes during thermal processing: a review. , Food Bioprocess Technol., Vol.8 ; pp.343-358
    20. Loh J. , Breene W.M. (1981) The thermal fracturability loss of edible plant-tissue - pattern and within-species variation. , J. Texture Stud., Vol.12 ; pp.456-471
    21. Lund D.B. (1982) Quantifying reactions influencing quality of foods - texture, flavor and appearance. , J. Food Process. Preserv., Vol.6 ; pp.133-153
    22. Ma L.Y. , Deng J.C. , Ahmed E.M. , Adams J.P. (1983) Canned shrimp texture as a function of its heat history. , J. Food Sci., Vol.48 ; pp.360-363
    23. Mallidis C.G. , Katsaboxakis C. (2002) Effect of thermal processing on the texture of canned apricots. , Int. J. Food Sci. Technol., Vol.37 ; pp.569-572
    24. Mansfield T. Karel M. , Fennema OR. , Lund DB. (1974) Principles of food science. Part II. Physical principles of food preservation, Marcel Dekker,
    25. Nisha P. , Singhal R.S. , Pandit A.B. (2006) Kinetic modelling of texture development in potato cubes (Solanum tuberosum L.), green gram whole (Vigna radiate L.) and red gram splits (Cajanus cajan L.). , J. Food Eng., Vol.76 ; pp.524-530
    26. Nourian F. , Ramaswamy H.S. (2003) Kinetics of quality change during cooking and frying of potatoes: Part 1. Texture. , J. Food Process Eng., Vol.26 ; pp.377-394
    27. Paulus K. , Saguy I. (1980) Effect of heat-treatment on the quality of cooked carrots. , J. Food Sci., Vol.45 ; pp.239-242
    28. Rao M.A. , Lund D.B. (1986) Kinetics of thermal softening of foods - a review. , J. Food Process. Preserv., Vol.10 ; pp.311-329
    29. Rao M.A. , Lee C.Y. , Katz J. , Cooley H.J. (1981) A kinetic-study of the loss of vitamin-C, color, and firmness during thermalprocessing of canned peas. , J. Food Sci., Vol.46 ; pp.636-637
    30. Rizvi A.F. , Tong C.H. (1997) A critical review ?" fractional conversion for determining texture degradation kinetics of vegetables. , J. Food Sci., Vol.62 ; pp.1-7
    31. Rodrigo C. , Mateu A. , Alvarruiz A. , Chinesta F. , Rodrigo M. (1998) Kinetic parameters for thermal degradation of green asparagus texture by unsteady-state method. , J. Food Sci., Vol.63 ; pp.126-129
    32. Rodrigo C. , Rodrigo M. , Fiszman S.M. (1997) The impact of high temperature, short time thermal treatment on texture and weight loss of green asparagus. , Z. Lebensm. Unters. F. A., Vol.205 ; pp.53-58
    33. Sawyer F.M. (1971) Interaction of sensory panel and instrumental measurement. , Food Technol., Vol.25 ; pp.247-248
    34. Silva C.A. , Bates R.P. , Deng J.C. (1981) Influence of pre-soaking on black bean cooking kinetics. , J. Food Sci., Vol.46 ; pp.1721-1725a
    35. Silva C.A. , Bates R.P. , Deng J.C. (1981) Influence of soaking and cooking upon the softening and eating quality of black beans (phaseolus-vulgaris). , J. Food Sci., Vol.46 ; pp.1716-1720b
    36. Solomon W.K. , Jindal V.K. (2003) Comparison of mechanical tests for evaluating textural changes in potatoes during thermal softening. , J. Texture Stud., Vol.33 ; pp.529-542a
    37. Solomon W.K. , Jindal V.K. (2003) Modeling thermal softening kinetics of potatoes using fractional conversion of rheological parameters. , J. Texture Stud., Vol.34 ; pp.231-247b
    38. Stoneham T.R. , Lund D.B. , Tong C.H. (2000) The use of fractional conversion technique to investigate the effects of testing parameters on texture degradation kinetics. , J. Food Sci., Vol.65 ; pp.968-973
    39. Suzuki K. , Kubota K. , Omichi M. , Hosaka H. (1976) Kinetic studies on cooking of rice. , J. Food Sci., Vol.41 ; pp.1180-1183
    40. Van Boekel M.A. (2008) Kinetic Modeling of Food Quality: A Critical Review. , Compr. Rev. Food Sci. Food Saf., Vol.7 ; pp.144-158
    41. Verlinden B.E. , De Baerdemaeker J. (1997) Modeling low temperature blanched carrot firmness based on heat induced processes and enzyme activity. , J. Food Sci., Vol.62 ; pp.213-218
    42. Verlinden B.E. , De Barsy T. , De Baerdemaeker J. , Deltour R. (1996) Modelling the mechanical and histological properties of carrot tissue during cooking in relation to texture and cell wall changes. , J. Texture Stud., Vol.27 ; pp.15-28
    43. Verlinden BE , Nicolai BM (1995) The starch gelatinization in potatoes during cooking in relation to the modeling of texture kinetics. , J. Food Eng., Vol.24 ; pp.165-179
    44. Villota R. , Hawkes J.G. Okos M.R. (1986) Physical and chemical properties of food., ASAE, ; pp.266-366
    45. Vu TS , Smout C , Sila DN (2004) Effect of preheating on thermal degradation kinetics of carrot texture. , Innov. Food Sci. Emerg. Technol., Vol.5 ; pp.37-44
    46. Yu K. , Wu Y. , Hu Z. , Cui S. , Yu X. (2011) Modeling thermal degradation of litchi texture: comparison of Well model and conventional methods. , Food Res. Int., Vol.44 ; pp.1970-1976