MEASURING RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY OF SOFT FOODS USING A MODIFIED SQUEEZE-FLOW APPARATUS ADRIANO SUN and SUNDARAM GUNASEKARAN 1 University of Wisconsin-Madison Food & Bioprocess Engineering Laboratory 460 Henry Mall, Madison, WI 53706 Accepted for Publication March 15, 2009 ABSTRACTjts_181 275..287 The rheological characteristics of soft foods such as peanut butter, mayonnaise and margarine were evaluated by the UW Meltmeter, a modified squeezing-flow device, applying a constant load under constant volume and changing area. Sample thickness profiles were measured as a function of time and the yield stress was calculated via a Herschel Bulkley model as a measure of spreadability. Peanut butter exhibited the highest yield stress value (1.31 1.45 kpa) followed by margarine (0.80 kpa) and mayonnaise (0.31 kpa). Calculated values of biaxial extensional viscosity and consistency coefficient values calculated from the UW Meltmeter data followed the expected trends of relative spreadability of peanut butter, margarine, and mayonnaise. These results show that the UW Meltmeter is a relatively simple and reliable device for measuring yield stress and rheological properties of soft foods. PRACTICAL APPLICATIONS A new method has been proposed to measure yield stress of spreadable foods. KEYWORDS Rheology, soft foods, spreadability, squeeze flow, UW meltmeter, yield stress 1 Corresponding author. TEL: 608-262-1019; FAX: 608-262-1228; EMAIL: guna@wisc.edu Journal of Texture Studies 40 (2009) 275 287. 2009, Wiley Periodicals, Inc. 275
276 A. SUN and S. GUNASEKARAN INTRODUCTION Soft, spreadable foods such as cream cheese and peanut butter (PB) are elastoplastic or viscoplastic materials. Consumer acceptance of these foods depends on their textural characteristics such as spreadability a measure of how easily and uniformly they can be deformed and spread at end-use temperatures. The rheological properties of spreadable foods have been studied by a variety of methods. Daubert et al. (1998) and Breidinger and Steffe (2001) used yield stress and yield strain data from vane measurements to construct texture maps of spreadable foods. Citerne et al. (2001) examined the rheological and physical characteristics of PB samples by creep, dynamic oscillatory test, and with equilibrium flow curves by fitting the Casson and Bingham models. However, semiliquid and soft foods like spreadable foods are often difficult to work with conventional rotational rheometers because of the possible wall slip and excessive sample disruption during loading into narrow gaps. The squeezing-flow rheometry has become a very attractive method to measure the rheological properties of semi-liquid and soft foods because wall slip is in fact a prerequisite for a proper test and due to its simplicity in interpreting results from sample thickness as a function of time. The technique was primarily used to measure the flow properties of highly viscous materials such as polymer melts, but has also been a useful method for semiliquid and semisolid food materials such as processed and melted cheeses, butter, dough, PB, etc. (Casiraghi et al. 1985; Wang et al. 1998; Gunasekaran and Ak 2003). The objectives of this research were to: (1) evaluate the rheological characteristics of soft, foods by a modified squeezing-flow device, named the UW Meltmeter; (2) compare yield stress measurements of the UW Meltmeter with traditional techniques such as the vane and shear stress ramp tests; and (3) relate rheological properties measured using UW Meltmeter to relative spreadability of the foods studied. MATERIALS AND METHODS Spreadable Food Materials Three store-bought commercial foods were used: Peanut Butter Delight (Aldi Inc., Batavia, IL), Shedd s Spread Country Crock margarine (Unilever, Englewood Cliffs, NJ) and Hellman s mayonnaise (Unilever, Englewood Cliffs, NJ). These products were chosen because they represent a range of subjective spreadability characteristics such as easy (mayonnaise [MY]), mild
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 277 1 (margarine [MG]) and semi-hard (PB), and are known to be self-lubricating because of their high fat content, which is a requirement of a proper squeezeflow test. PB tested was from six containers: four were labeled as C2, C3, C4 and C5, which were manufactured at the same time (lot date of January 4, 2007); the other two were labeled C1 and C6, and were from lot dates of December 19, 2006 and February 16, 2007, respectively. MY and MG tested were from only one container each. All materials were stored in a refrigerator (~4.0C) prior to testing. MG and MY were tested at refrigerated temperatures (5.0 and 7.5C, respectively). PB was left at room conditions for 2 3 h for temperature equilibration before testing at room temperature (22C). Temperatures of all samples were monitored periodically to ensure that maximum variation was less than 1.0C. When samples reached the upper temperature limit, they were brought back to the refrigerator. Prior to measurements, samples were allowed to rest for 5 min upon loading. In addition, all experiments were performed in triplicates (n = 3) following a balanced randomized run order with respect to replicates and measurement techniques to account for structural and consistency variations throughout testing. Modified Squeezing-Flow Apparatus A modified squeezing-flow device, originally designed and fabricated to measure the melt and flow behavior of cheeses, was used. This device, called the UW Meltmeter, was developed in the Food and Bioprocess Engineering Laboratory of the Biological Systems Engineering Department at the University of Wisconsin-Madison (Wang et al. 1998). A schematic diagram of the UW Meltmeter and photograph of the system are shown in Figs. 1 and 2, respectively. It consists of one or two sample wells (7 mm in depth and 30 mm in diameter) with a stationary inner ring and a moveable outer ring. In each unit, the inner ring is equipped with an electric heater operated by a temperature controller (CN 4400, Omega Engineering, Inc., Stanford, CT). A stationary center piston allows the outer ring to move up and down through a lever. When the outer ring is in the UP position by rising the lever, a sample well is formed in which sample can be easily loaded without excessive structural disruption (Fig. 3a). When the lever is lowered, the outer ring is in the DOWN position exposing a cylindrical sample specimen of the shape and dimensions of the sample well (Fig. 3b). A linear variable differential transformer (LVDT; DC-E 1000, Schaevitz Engineering, Pennsauken, NJ) is attached to a 66-mm-diameter circular top plate that rests on the outer ring. By lowering the lever, the outer ring is lowered and the top of the sample specimen is in contact with the top plate
278 A. SUN and S. GUNASEKARAN To computer data acquisition LVDT Circular top plate 7 mm 30 mm LVDT rod Outer ring (flow platform) Sample Stationary center piston Moveable outer ring-annulus Heater Lever Arm Power supply FIG. 1. SCHEMATIC DIAGRAM OF THE SAMPLE WELL UNIT OF THE UW MELTMETER (GUNASEKARAN AND AK 2003) LVDT, linear variable differential transformer. Circular top plate LVDT rod Lever FIG. 2. PICTURE OF THE MODIFIED SQUEEZING-FLOW DEVICE, THE UW MELTMETER, WITH TWO SAMPLE WELLS LVDT, linear variable differential transformer.
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 279 Outer ring 30-mm 7-mm Inner ring a b FIG. 3. (a) SAMPLE WELL IN THE UW MELTMETER WHEN THE OUTER RING IS IN UP POSITION AND (b) SAMPLE SPECIMEN FORMED WHEN OUTER RING IS IN DOWN POSITION FIG. 4. (a) STARTING POSITION WHILE SAMPLE IS RESTING TO RELEASE SHEAR EFFECTS; (b) LEVER IS LOWERED AND THE UPPER PLATE IS IN CONTACT WITH SAMPLE SPECIMEN; (c) LOAD GOVERNED BY UPPER PLATE SQUEEZES SAMPLE SPECIMEN (GUNASEKARAN AND AK 2003) (Fig. 4). The force applied by the circular top plate, the LVDT rod and any additional weights, causes the specimen to be compressed and squeezed. The UW Meltmeter test is performed applying a constant load under constant volume and changing area. During UW Meltmeter tests PB, MG and MY samples were scooped from their original containers at respective experimental temperatures and were gently placed into the sample well. Special care was taken to load the samples consistently and to prevent the presence of voids or air pockets. A metal spatula
280 A. SUN and S. GUNASEKARAN was used to smooth the top surface and to remove any excess material. The LVDT monitors the diminishing height of the specimen during deformation. A computer data acquisition system (DAS 16G High Speed Analog I/O Board, Metrabyte Corp., Taunton, MA) and software (Easyest LX Software, Asyst Technologies, Inc., Rochester, NY) were used to collect the sample height versus time data. Measurements were performed in three replicates (n = 3) at a range of applied forces (e.g., 0.74 1.62 N) for each spreadable food. Yield Stress Measurements The determination of the yield stress from the UW Meltmeter is based on the residual stress after deformation and takes into consideration that a viscous-plastic fluid (i.e., Bingham plastic and Herschel-Bulkley fluid) flows until it reaches an equilibrium height or limiting thickness (H L ). The yield stress parameter of spreadable foods, denoted here as the apparent extensional yield stress, was determined by applying the analytical expression for Herschel-Bulkley fluids proposed by Ak and Gunasekaran (2000) shown in Eq. (1), using the H L value obtained from the sample height versus time data. σ y L = WH 2 π RH o o (1) From the equation above, s y is the so-called apparent extensional yield stress, W is the load applied on the circular top plate, R o is the initial radius of the cylindrical specimen and H o is the initial height of the specimen sample. The mean value of the last 50 data points (20 s) was taken as the H L value, as long as the recorded height did not change by more than 5% during that time. A dynamic rheometer, Bohlin C-VOR (Malvern Inc., Southborough, MA), was used to compare yield stress values with the vane and stress ramp tests. For the vane test, the rheometer was equipped with a four-blade vane (35-mm high and 10 mm in diameter) and was operated in the controlled-rate mode while collecting 150 data points. Samples were loaded into a standard stainless steel holding cup of 62-mm high and 27 mm in diameter. The tests were run at constant rotational speed of 0.5 rpm and angular velocity of 0.052 rad/s. For the vane test, triplicate measurements were made with the top surface of vane blades aligned with the sample surface taking new samples each time. Triplicate of tests were performed using the single-point determination, based on the total torque required to overcome the yield stress of the fluid, using the following expression (Dzuy and Boger 1983, 1985)
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 281 σ o,vane 2Mo = h 1 + 3 πd d m+ 3 1 (2) where, s o,vane is the vane yield stress, M o is the maximum or peak torque, h and d are the height and diameter circumscribed by the tip of the vane blades, respectively, and m is a constant. The assumption of m = 0 has shown to be reliable in previous studies (Steffe 1996; Daubert et al. 1998; Corradini and Peleg 2005); thus, the yield stress was determined with Eq. (3): σ o,vane 6Mo = 2 πd ( 3h+ d) (3) For the stress ramp test, a 25-mm-diameter serrated-plate geometry (SP 25-mm, Malvern Instruments, Westborough, MA) attached to the rheometer was used with a 1-mm gap. The test was performed by applying gradually increasing stresses onto the sample under controlled-stress mode and collecting 100 data points. The instantaneous (or apparent) viscosity was monitored for the presence of an inflexion point, indicating onset of flow. The yield stress (s o,ramp ) was given by the stress at the inflexion point in the instantaneous viscosity versus shear stress curve. RESULTS AND DISCUSSION Yield Stress Determination A typical height versus time response is shown in Fig. 5 for PB (e.g., PB-C5) under the constant force of 1.62 N. The specimen height diminishes exponentially with time, and the initial rate of decrease of sample thickness is inversely proportional to the applied load. The apparent extensional yield stress (s y ) values measured with the UW Meltmeter are listed in Table 1 along with those measured by the vane and stress ramp tests. The magnitudes of yield stress determined from all three techniques followed similar trend MY exhibited the lowest yield stress followed by MG and PB. The yield stress values of all three spreadable materials determined from the UW Meltmeter and the vane tests were in close agreement. However, this was not true when comparing the yield stress obtained from the stress ramp test, as values were lower by more than an order of magnitude than the corresponding values from the other two tests for all three foods. This suggests that yield stress values obtained from shear stress
282 A. SUN and S. GUNASEKARAN 7.0 Run 1 Run 2 Run 3 Sample height (mm) 6.5 6.0 5.5 5.0 0 5 10 Time (s) 15 FIG. 5. TYPICAL SAMPLE HEIGHT VERSUS TIME RESPONSE FROM THE UW MELTMETER OF PEANUT BUTTER (PB-C5) UNDER SQUEEZING FORCE OF 1.62 N TABLE 1. YIELD STRESS VALUES OF PEANUT BUTTER (PB-C1 THROUGH PB-C6), MARGARINE (MG), AND MAYONNAISE (MY) DETERMINED BY THE UW MELTMETER, VANE AND STRESS RAMP TESTS Product Manufacturing lot date UW Meltmeter Load (N) s y (kpa) Vane method s o,vane (kpa) Stress ramp s o,ramp (kpa) PB-C1 December 19, 2006 1.44 1.27 0.07 1.37 0.04 0.51 0.08 PB-C2 January 4, 2007 1.44 1.31 0.04 1.35 0.12 0.56 0.03 1.62 1.43 0.06 PB-C3 January 4, 2007 1.62 1.43 0.06 1.40 0.06 0.59 0.05 PB-C4 January 4, 2007 1.62 1.40 0.06 1.42 0.10 0.53 0.02 1.80 1.48 0.10 PB-C5 January 4, 2007 1.44 1.31 0.04 1.39 0.07 0.56 0.10 1.62 1.43 0.06 PB-C6 February 16, 2007 2.07 1.79 0.05 1.71 0.03 0.73 0.05 MG 0.74 0.80 0.05 0.86 0.07 0.25 0.02 MY 0.30 0.31 0.02 0.34 0.005 0.12 0.005 Average of three replicates standard deviation.
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 283 7.0 Limiting thickness, HL (mm) 6.0 5.0 4.0 R 2 = 0.93 R 2 = 0.99 3.0 PB * R 2 MG = 0.95 MY 2.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Initial Applied Stress (kpa) FIG. 6. LIMITING THICKNESS VERSUS INITIAL APPLIED STRESS FOR PEANUT BUTTER (PB), MARGARINE (MG) AND MAYONNAISE (MY) * Includes averaged values from all PB containers. ramp curves are to be considered as the static yield stress; the dynamic yield stress tests were performed under excessive sample disruption during loading and measurement, for example, while lowering the upper serrated plate to a 1-mm gap size. The limiting thickness H L decreased with initial stress (s i ) applied, and correspondingly the calculated s y values increased. These trends are presented in Figs. 6 and 7, respectively. The initial applied stress was determined dividing the applied load (W) used by the initial area of the specimen (707 mm 2 ). Under high applied load the material is deformed and squeezed at a great extent such that the material instantly collapses. On the other hand, under low load deformation is rather slow requiring long time to reach H L, which suggests that the material is able to resists the applied stress, and perhaps able to rearrange and rebuild its structure somewhat. Therefore, at higher loads s y may be regarded as the static measurement considering both its magnitude (i.e., higher yield stress) and state of physical structure of the material. We recognize that the magnitude of s y of the materials tested with the UW Meltmeter may be lower than the true yield stress, as s y accounts for the remaining material structure after squeezing.
284 A. SUN and S. GUNASEKARAN Apparent extensional yield stress (kpa) 2.0 1.5 1.0 0.5 R 2 = 0.97 R 2 = 0.99 R 2 = 0.99 PB* MG MY 0.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Initial applied stress (kpa) FIG. 7. APPARENT EXTENSIONAL YIELD STRESS VERSUS INITIAL APPLIED STRESS FOR PEANUT BUTTER (PB), MARGARINE (MG) AND MAYONNAISE (MY) * Includes averaged values from all PB containers. Yield Stress of PB as a Function of Manufacturing Date Results from all the three measurement methods indicate that variations in PB manufactured on different lot dates are some what predictable based on yield stress (seen in Table 1). The yield stress of PB of lot dates January 4, 2007 (PB-C2, C3, C4 and C5) and December 19, 2006 (PB-C1) were similar; PB of lot date February 16, 2007 (PB-C6) had the highest yield stress. The reasons for these deviations on yield stress are unknown and beyond the scope of our study, yet possible physical and chemical changes during storage, transportation, and distribution may have contributed to structural changes of PB, not to mention variations in the manufacturing processes. Biaxial Extensional Strain Rate and Biaxial Extensional Viscosity The data of specimen height as a function of time (Fig. 5) obtained from the UW Meltmeter were converted into biaxial extensional strain rate ( ε B ) and biaxial extensional viscosity (h B ) with Eqs. (4) and (5), respectively, for all three products (Chatraei et al. 1981; Dealy 1995):
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 285 1e+7 Biaxial extensional viscosity (Pa-s) 1e+6 1e+5 1e+4 MY MG PB-C5 1e+3 0.0001 0.001 0.01 0.1 1 Biaxial extensional strain rate (1/s) FIG. 8. PLOT OF h B VERSUS ε B FOR MAYONNAISE (MY), MARGARINE (MG), AND CONTAINER 5 OF PEANUT BUTTER (PB-C5) UNDER SQUEEZING LOADS OF 0.30 N, 0.74 N AND 1.62 N, RESPECTIVELY 1 1 1 dh [] t εb = εt = [] 2 2 Ht dt (4) η B = H 2 2 FH () t dh () t Ao dt o ( ) where e T is the momentary strain rate, F is the squeezing flow force, H(t)isthe height at time t, H o is the initial height of the specimen, and A o is the initial cross-sectional area of the specimen. Figure 8 displays the h B as a function of ε B for MY, MG and PB under squeezing loads of 0.30, 0.74 and 1.62 N, respectively. It was observed that increasing squeezing loads resulted in increasing h B, causing the curve to shift upwards; conversely, low loads resulted in lower h B and curves tended to move downwards. Overall, biaxial extensional viscosity of MY was the lowest followed by those of MG and PB, which follows the expected relative spreadability values for these foods. (5)
286 A. SUN and S. GUNASEKARAN TABLE 2. CALCULATED POWER-LAW CONSTANTS n AND K FOR PEANUT BUTTER (PB), MARGARINE (MG) AND MAYONNAISE (MY) BASED ON THE UW MELTMETER DATA Product Manufacturing lot date Load (N) Flow behavior index, n Consistency coefficient, K (Pa s n ) PB-C1 December 19, 2006 1.44 0.019 0.005 1,796 13.7 PB-C2 January 4, 2007 1.44 0.022 0.009 1,880 82.1 1.62 0.015 0.005 2,030 104.3 PB-C3 January 4, 2007 1.62 0.011 0.002 1,978 112.9 PB-C4 January 4, 207 1.44 0.014 0.004 1,844 46.2 1.62 0.014 0.001 1,978 21.4 PB-C5 January 4, 2007 1.62 0.018 0.004 1,922 13.2 1.80 0.024 0.001 2,195 151.2 PB-C6 February 16, 2007 2.07 0.021 0.006 1,987 97.3 MG 0.73 0.010 0.007 942 34.7 MY 0.30 0.018 0.005 399 8.8 Power-Law Constants The consistency coefficient (K) and flow behavior index (n) were determined assuming the squeezing-flow behavior of the selected spreadable foods follows the power-law equation (Table 2). The n values were in the order of 0.010 to 0.024 indicating highly shear thinning behavior for all three spreadable products. These values were much lower than previously reported in the literature based on shear flow measurements (Campanella and Peleg 1987; Steffe 1996). On the other hand, coefficient K represented well the differences in consistency between products, as values for PB were about two and six times greater than those of MG and MY, respectively. CONCLUSIONS The rheological properties of three soft, spreadable foods were measured with a modified squeezing-flow device, the UW Meltmeter. Information of sample height versus time and the resulting limiting thickness (H L ) obtained were used to calculate the apparent extensional yield stress parameter based on a Herschel Bulkley model. In general, decreasing deformation rate of sample thickness as a function of time was proportional to the increase in applied squeezing force; and decreasing H L resulted in increasing yield values. The magnitude of yield stress increased in the following the order: MY (0.31 kpa); MG (0.80 kpa); and PB (1.31 1.45 kpa). Additionally, yield stress values of selected spreadable foods measured with the UW Meltmeter were in good agreement with values obtained with the vane method. The biaxial extensional
RHEOLOGICAL CHARACTERISTICS AND SPREADABILITY 287 viscosity calculated from UW Meltmeter data followed the expected relative spreadability values for the three foods tested. REFERENCES AK, M.M. and GUNASEKARAN, S. 2000. Simulation of lubricated squeezing flow of a Herschel-Bulkley fluid. Appl. Rheol. 10(6), 274 279. BREIDINGER, S.L. and STEFFE, J.F. 2001. Texture map of cream cheese. J. Food Sci. 66, 453 456. CAMPANELLA, O.H. and PELEG, M. 1987. Squeezing flow viscosimetry of peanut butter. J. Food Sci. 52, 180 184. CASIRAGHI, E.M., BAGLEY, E.B. and CHRISTIANSON, D.D. 1985. Behavior of mozzarella, cheddar, and processed cheese spread in lubricated and bonded uniaxial compression. J. Texture Studies 16(3), 281 301. CHATRAEI, S., MACOSKO, C.W. and WINTER, H.H. 1981. Lubricated squeezing flow: A new biaxial extensional rheometer. J. Rheol. 25, 433 443 CITERNE, G.P., CARREAU, P.J. and MOAN, M. 2001. Rheological properties of peanut butter. Rheologica Acta V40, 86 96. CORRADINI, M.G. and PELEG, M. 2005. Consistency of dispersed food systems and its evaluation by squeezing flow viscometry. J. Texture Studies 36(5 6), 605 629. DAUBERT, C.R., TKACHUK, J.A. and TRUONG, V.D. 1998. Quantitative measurement of food spreadability using the vane method. J. Texture Studies 29, 427 435. DEALY, J.M. 1995. Official nomenclature for material functions describing the response of a viscoelastic fluid to various shearing and extensional deformations. J. Rheol. 39, 253 265. DZUY, N.Q. and BOGER, D.V. 1983. Yield stress measurement for concentrated suspensions. J. Rheol. 27, 321 349. DZUY, N.Q. and BOGER, D.V. 1985. Direct yield stress measurement with the vane method. J. Rheol. 29, 335 347. GUNASEKARAN, S. and AK, M.M. 2003. Cheese Rheology and Texture. CRC Press, New York, NY. STEFFE, J.F. 1996. Rheological Methods in Food Process Engineering. Freeman Press, East Lansing, MI. WANG, Y.C., MUTHUKUMARAPPAN, K., AK, M.M. and GUNASEKA- RAN, S. 1998. A device for evaluating melt/flow characteristics of cheeses. J. Texture Studies 29, 43 55.