10,884
Views
105
CrossRef citations to date
0
Altmetric
Original Articles

Yield Stress in Foods: Measurements and Applications

&
Pages 70-101 | Received 06 Apr 2007, Accepted 29 Jun 2008, Published online: 08 Jul 2010

Abstract

Though the presence of true yield stress has been debated, it has been accepted as an engineering reality. Now, yield stress is routinely measured and used in the food industry not only for basic process calculations and manufacturing practices, but also as a test for sensory and quality indices and to determine the effect of composition and manufacturing procedures on structural and functional properties. In this article, we present a comprehensive review of various measurement methods and the applications of yield stress in many practical situations. The measurement methods have been grouped under traditional and recent methods. An emphasis was placed on the development of measurement methods applied to a variety of foods, and the use of yield stress as a parameter to evaluate certain food material attributes: flow behavior, end-use quality, processability, composition, and so on.

INTRODUCTION

Yield stress is a well recognized physical and rheological property for liquid and solid materials. The yield stress is classically defined as the minimum shear stress that must be applied to the material to initiate flow.[Citation1] Materials that exhibit yield stress are typically multiphase systems such as cements, soils, paints, pastes, printing inks, greases, and a variety of food products such as salad dressings, sauces, and spreads.

The idea of yield stress was introduced in the 1920s by the work of Bingham, Schwedoff, and Green, and since then, its concept and definition have been continuously debated, especially after Barnes and Walters[Citation2] challenged its existence in their article “The Yield Stress Myth?” The controversy hinges on the early statements describing solid materials as those exhibiting “no continuous deformation (or flow) upon decreasing stresses below the yield stress value.” The debate continued thereafter with the works by Evans,[Citation3] Astarita,[Citation4] Schurz,[Citation5] and Spaans and Williams.[Citation6] The existence of yield stress in fluids still remains controversial; however, there is much to validate that yield stress is an “engineering reality”.[Citation7] For instance, flocculated suspensions and dispersions usually exhibit yield stress as individual particles tend to interact and aggregate resulting in flocs and consequently form a three-dimensional network structure.[Citation8,Citation9] Therefore, yield stress may also relate to the strength of the coherent network structure throughout the volume of the material as the force per unit area required to break down the structure.[Citation1,Citation8,Citation10,Citation11]

In addition, yield stress is generally denoted as the transition stress at which a material behaves either as elastic solid-like or viscous liquid-like. This transition usually occurs within a range of stresses, in which the material exhibits viscoelastic behavior.[Citation1] This transition behavior may be distinguished based on the Deborah number (De), which is defined as the ratio of the characteristic time of the material to the characteristic time of the deformation process. Thus, a material can behave as either a solid-like at high De (i.e., elastic deformation with short characteristic time) or liquid-like at low De (i.e., viscous flow with long characteristic time). Consequently, the timeframe of yield measurement is also a key parameter.[Citation12] Cheng[Citation12] also suggested that the magnitude of the measured yield values is closely related to creep, stress growth, thixotropic breakdown and recovery, and characteristic times of the transient responses. These controversial debates notwithstanding, the quantitative knowledge of yield stress is used extensively in diverse areas of study—colloidal science, food engineering, biomedical sciemces, and geology. The idea of yield stress is now a well accepted physical reality. However, as cautioned by Barnes,[Citation9] defining yield stress and using mathematical models to characterize material behavior should only be used over an appropriate range of shear rates and under proper applicability.

In the food industry, the definition of yield stress as the minimum shear stress required to initiate flow,[Citation8,Citation13,Citation14] is generally accepted due to the limited time-scale of most food processing operations (e.g., emptying a vat) and other food-related activities (e.g., eating).[Citation9] The idea of yield stress as the transition stress and as a measure of the structural strength of food dispersions like mayonnaise and tomato puree are also well known.[Citation11,Citation15,Citation16,Citation17] Cheng[Citation12] introduced the concept of static and dynamic yield stresses, which is useful for food processing applications (). The former represents the yield stress of an undisturbed sample; and the latter, a broken down sample. In general, the magnitude of static yield stress is significantly higher than that of the dynamic yield stress. This concept is important considering that most foods exhibit time-dependent rheological properties and undergo physicochemical changes during storage and/or processing.[Citation14,Citation15]

Figure 1 Static and dynamic yield stress (from[Citation15]).

Figure 1 Static and dynamic yield stress (from[Citation15]).

A precise, quantitative understanding of yield stress in foods is important for several practical applications—from process calculations to product development. For example, the knowledge of yield stress is helpful in estimating how well a fluid food, such as tomato ketchup, drains from a bottle.[Citation18] There are several other applications of yield stress in foods: thickness of coating layer (chocolate on ice cream bar, glazing on doughnuts and cakes), settling of particles in fluids (spices in salad dressing, chocolate in chocolate milk), spreadability (cream cheese, mayonnaise), mouthfeel (creaminess of yoghurt), and process design calculations (flow and velocity profiles through pipelines).

MEASURING YIELD STRESS IN FOODS

The most common method to obtain yield stress is to extrapolate the shear stress versus shear rate curve back to zero shear rate. However, extrapolated yield stress values are strongly influenced by the rheological model chosen.[Citation15, Citation19] An alternative to this is to measure yield stress directly, which relies on an independent assessment of the stress at which the material starts to flow.[Citation1] Direct measurements with conventional rheometers involve testing a material with creep/recovery, stress relaxation, and stress growth techniques. Other methods such as cone penetration, vane method, and dynamic oscillatory measurements have also been studied.[Citation8,Citation12,Citation13,Citation20] A modified squeezing flow apparatus,[Citation21] slump tests,[Citation22] and measurements of static yield stress of suspensions using a slotted-plate device,[Citation23] are among the recently developed methods.

Two major problems encountered in the determination of the yield stress lie on the reproducibility of the experimental data and on the vastly different yield values obtained from different techniques.[Citation15] As indicated by Cheng,[Citation12] yield stress is a time-dependent property; hence, a proper time scale of the experiment must be selected to minimize discrepancy among measurements. Furthermore, most suspensions show thixotropic characteristics, which cause development of a more ordered structure or structural recovery if the applied strain rate is small (i.e., long measurement time).[Citation19]

To make matters worse, fluids that exhibit yield stress are often more difficult to measure and the type of viscometer may impose additional limitations due to possible wall slip.[Citation24] These effects are associated with the displacement of the dispersed phase away from solid boundaries due to surface, hydrodynamic, viscoelastic, chemical, and gravitational forces. As a result, the liquid near the walls are depleted of particles and have lower viscosity, leading to higher measured shear rates at fixed shear stresses.[Citation19] Improvements to existing methods and new methods have been developed to obtain reliable and precise yield stress measurements in foods.

Traditional and Common Methods

Extrapolation

Extrapolation of the shear stress versus shear rate data obtained from conventional rheometers is the most common indirect technique to measure yield stress. The experimental data, also referred to as equilibrium flow curves, can be interpreted with or without a rheological model such as the Bingham, the Herchel-Bulkley or the Casson. Extrapolating the experimental data back to zero shear rate and obtaining the yield stress value at the shear stress intercept is an easy and quick procedure; however, it can only be done for linear data.[Citation24] An alternative numerical technique is to plot apparent viscosity (η a ) versus shear stress and determine the yield value at the point where η a tends to infinity.[Citation15] Similarly, Nguyen and Boger[Citation24] pointed out that a slope of −1 for a logarithm plot of η a versus shear rate curve over a range of low shear rate often indicates the presence of yield stress (e.g., Bingham fluids). The use of rheological models is helpful for proper extrapolation. However, it is important to reiterate that the selection of a particular model depends on the shear rate range and on the intended application.[Citation9,Citation24]

Overall, the goodness of the extrapolated yield stress value relies on the accuracy of the experimental data in the low shear rate range. In addition, as mentioned earlier, wall slip effects, fracture, and expulsion of the fluid sample are some of the problems encountered while using conventional rheometers, which may affect not only the extrapolated yield stress value, but also the rheological flow behavior. Nguyen and Boger[Citation24] described in detail the main features and problems with most common types of viscometers, and proposed some methods to correct these problems. However, the yield stress obtained from extrapolated technique cannot be regarded as an absolute material property or “true yield stress”.[Citation24] The extrapolated yield stress is often used as a quick estimate and to compare results from other methods.

Stress relaxation test

In this method, the fluid material is first sheared at either constant shear rate or constant shear stress in a conventional rotational viscometer, followed by bringing the material to rest either gradually or suddenly. The yield stress, also called as the relaxed stress, is then measured as the limiting or residual stress remaining in the fluid upon cessation of flow.[Citation1] During the test, it is important to ensure that the initial applied stress is higher than the yield stress. When the shear stress in the relaxed state decays to a constant level it represents the yield value. This method has been primarily applied to study linear viscoelasticity of food materials, with little work done for yield stress measurements. One of the reasons is the pre-shearing required, which damages the material structure. Therefore, stress relaxation technique is suitable only for fluids whose yield stress is not dependent on shear history.[Citation24]

Creep/recovery response

The creep/recovery response consists of applying a constant shear stress in steps to the fluid material as illustrated in . If stresses applied are below the yield stress, the test material behaves as an elastic solid with a complete recovery upon removal of stress and it will not flow. This is observed by a horizontal curve with slope nearly zero (τ1 and τ2). When the applied stresses are greater than the yield value, the deformation or strain increases indefinitely with time and viscous flow ensues (τ4). The yielding point can be detected from a drastic change of slope in the time vs. shear strain curve; sometimes, it is difficult to identify the critical shear stress as it approaches the yield stress. It can be seen in that the yield stress lies between τ2 and τ3. Overall, creep/recovery test is a more sensitive and less destructive than other methods for measuring yield stress. However, measurement duration can be very long to determine the steady creep compliance; influence the accuracy of the yield values due to shear history and thixotropic; and possibly cause some physicochemical effects.

Figure 2 Typical strain-time curves obtained during yield stress measurement by the creep-recovery test (from[Citation1]).

Figure 2 Typical strain-time curves obtained during yield stress measurement by the creep-recovery test (from[Citation1]).

Shear Stress Ramp Test

This method is an alternative to creep/recovery technique to shorten the measurement duration. The test is performed by applying gradually increasing stresses onto the sample. The instantaneous (or apparent) viscosity is monitored for the presence of an inflexion point, which indicates onset of flow, and the yield stress is given at the corresponding shear stress (). The measurement of the stress ramp curves are often performed with parallel plate or cone-plate geometries. Possible errors involved are slipping of the sample and effects due to gap size.

Figure 3 Stress ramp curve for margarine (from[Citation21]).

Figure 3 Stress ramp curve for margarine (from[Citation21]).

Cone Penetrometer

Cone penetration is a popular technique used for testing strength properties of soils and other materials of high consistency. In this test, a metal cone of specific dimensions (angle α and height (H) and weight (W) is forced into a flat and smooth horizontal surface of test specimen (). The penetration can be performed in three modes of operation[Citation25]: (1) a metal cone is let to sink and the depth of penetration is measured after a fixed time; (2) a cone is released into the sample and the depth of penetration is measured when the cone is at rest due to the yield stress of the material; and (3) a cone is pushed into the sample at a constant speed and the force required for the penetration is recorded. In the constant force mode, the cone penetrates quickly at first, but it gradually slow down until it comes to rest. The “apparent yield stress” (σ app ) is calculated as a function of the penetration depth h when the cone stops[Citation25]:

Figure 4 Illustration of a cone penetrometer assembly (from[Citation25]).

Figure 4 Illustration of a cone penetrometer assembly (from[Citation25]).

(1)

The cone penetration is considered an empirical method for rapid evaluation of consistency of a variety of solid and semi-solid foods. In fact, this method has been employed extensively for measuring the stiffness of dairy products.[Citation24] As for measuring the yield stress, the method is also considered suitable for fluids with high yield stress and low viscosities as the viscous resistance during cone penetration is minimal. Nguyen and Boger[Citation24] points out that viscous effects retard cone penetration resulting in over estimation of the yield value. On the other hand, it is disadvantageous when measuring highly viscous materials, such as butter, since it might take long time for the cone to reach an equilibrium condition on its own weight.

Other methods based on similar principle of forcing an object through the material under test have been investigated. For example, the falling needle viscometer, which is used to measure the steady shear viscosity of Newtonian and non-Newtonian viscoelastic fluids, has been attempted for yield stress measurements.[Citation26] This method can provide an indirect measure of yield stress, once the density difference between the needle and the liquid is corrected for. However, the needle is not stationary in the liquid.

Uhlherr et al.[Citation27] proposed a new method for static measurement of yield stress using a cylindrical penetrometer. The principle is similar to that of the falling needle viscometer—the penetrometer falls under gravity toward an equilibrium position after being gently released in a fluid. Once the penetrometer is partially immersed, it reaches static equilibrium at the position where the yield stress is achieved. The yield values measured from this technique are comparable to those values obtained from other tests such as vane, extrapolation, stress ramp, and creep for Carbopol and titatnium dioxide solutions.

Dynamic oscillatory testing

Dynamic mechanical measurement is an alternative method for the determination of yield stress as it provides useful information about the structural and deformation behavior of viscoplastic fluids at below or near the yield stress.[Citation24] During oscillatory testing, the material is subjected to a sinusoidal strain and the resulting stress is measured as a function of both time and frequency. Depending on the strain amplitude (γ o ), small deformations occur within the linear viscoelastic region, whereas large deformations will be in the non-linear region.

Dynamic measurements in the linear viscoelastic region of multiphase and composite materials exhibit a low frequency plateau, where the linear viscoelastic terms are independent of frequency.[Citation24] The presence of this plateau has been suggested to be due to the existence of a network structure, which can be correlated to the yield stress (), whereas plateaus are absent for non-yield stress fluids.[Citation24]

Figure 5 Dynamic oscillatory measurement (from[Citation15]).

Figure 5 Dynamic oscillatory measurement (from[Citation15]).

Dynamic measurements are not popularly used for yield stress determination, because it is difficult to achieve the extremely small amplitude required so that the stress and strain sinusoidal waveforms are completely in phase.[Citation1] On the other hand, it might be an advantageous procedure as the small deformation causes minimum disturbances to the structure of the fluid material.[Citation24] Furthermore, Liddell and Boger[Citation1] pointed out that measurement in the non-linear viscoelastic region can yield superior results, once the plateau where the non-linear viscoelastic functions are independent of frequency.

Static stability on an inclined plate

Nguyen and Boger[Citation24] described a technique to measure the yield stress of viscoplastic fluids on an inclined plate. This method takes into consideration the shear stress of the fluid perpendicular to an inclined plane of angle α at a position y measured from the free surface end. This relation is given as[Citation24]:

(2)

where, ρ is the density of the fluid material and g is the acceleration of gravity. Under the above conditions, flow only occurs when shear stress is greater than the yield stress τ (y) > τ o . The yield stress is determined by placing a uniform fluid layer on an initially horizontal plane, followed by progressively increasing the angle of inclination of the plane until a critical value is reached, whereby the fluid starts to flow (). An alternative method is to keep the angle of inclination α fixed and measure the final fluid depth when flow stops after drainage. Although these techniques seem to be relatively simple, accuracy of the measurements is strongly influenced by the precision in measuring the critical angle and the fluid thickness of the.

Figure 6 Schematic of the flow in an inclined plate (from[Citation28]).

Figure 6 Schematic of the flow in an inclined plate (from[Citation28]).

A very similar but different method has also been attempted—withdrawing a plate immersed in a fluid bath at constant speed and measuring the thickness of the material adhering to the plate (). The plate remains in contact with the liquid bath throughout the process, and the amount of material coated to the plate is related to the yield stress.

Figure 7 Illustration of post-withdrawal drainage process (from[Citation29]).

Figure 7 Illustration of post-withdrawal drainage process (from[Citation29]).

Lang and Rha[Citation30] further investigated this method and concluded that the yield stress determined from the thickness of the fluid coated on the vertical plate is not comparable to those obtained with other techniques. De Kee et al.[Citation29] examined the post-withdrawal drainage behavior of Newtonian and viscoelastic fluids and viscoelastic suspensions with a computer-driven laser interferometry, which was used to investigate the transient film thickness. Their results also indicated that the plate-withdrawing technique ineffective for calculating yield stress. Some problems encountered were related to interfacial interactions, rate of withdrawal from the fluid, and draining time of material from the plate.[Citation14]

Vane technique

The vane method was originally employed to eliminate slipping of the sample while studying the coagulation of sodium clay gels[Citation31]; it was further developed to measure flow properties and shear strength of soils.[Citation8,Citation31,Citation32] The vane geometry consists of a vane spindle with typically four to eight thin blades arranged at equal angles, centered on a narrow cylindrical shaft of a rheometer or viscometer. shows a schematic of a typical four-bladed vane and how it is used during measurement.

Figure 8 Schematic of a typical four-bladed vane (left) and how it is used during measurement (right). [H = vane height, D = vane diameter, DT = container diameter, Ω = angular velocity, (a): immersed in the sample (Z1 and Z2 height of sample above and below the vane, respectively); (b): top surface in level with sample] (from[Citation25]).

Figure 8 Schematic of a typical four-bladed vane (left) and how it is used during measurement (right). [H = vane height, D = vane diameter, DT = container diameter, Ω = angular velocity, (a): immersed in the sample (Z1 and Z2 height of sample above and below the vane, respectively); (b): top surface in level with sample] (from[Citation25]).

The vane can be operated under controlled shear-rate mode or controlled shear-stress mode; in both cases the yield stress is the minimum stress required for continuous rotation of the vane. In the controlled shear-stress mode, a constant stress is applied in steps to the vane immersed in the sample. The resulting data corresponds to the strain (or angle of rotation) as a function of time for each increased stress step. The yield stress is determined as the stress corresponding to the point at which there is a rapid increase in the creep angle or strain for example between τ4 and τ5 in . Vane method also minimally disturbs the sample upon insertion of the vane spindle, which is critical when measuring weakly and delicately structured materials such as foams and emulsions.[Citation8,Citation25,Citation31,Citation32] Sample disturbance may also be minimized by using the original product containers such as for foods.[Citation33] The vane method has been suggested as a standard method for determining the yield stress of foods.[Citation34]

Figure 9 Typical response curve for yield stress measurement with the vane method under controlled shear-stress mode. Shear stress τ1 through τ5 are applied each for certain duration.

Figure 9 Typical response curve for yield stress measurement with the vane method under controlled shear-stress mode. Shear stress τ1 through τ5 are applied each for certain duration.

In the controlled shear-rate mode, a constant rotational speed is applied to the vane, and the resulting torque is measured as a function of time (). Studies have shown that rotational speeds between 0.1 and 8 rpm produce superior results, and yield stress values measured were essentially constant.[Citation14] The yield stress can therefore be obtained through a quick single-point determination, based on the total torque to overcome the yield stress (σ o ) of the fluid. A simplified equation is given as:

Figure 10 Typical torque versus time response curve for yield stress measurement with the vane method under controlled shear-rate mode (from[Citation25]).

Figure 10 Typical torque versus time response curve for yield stress measurement with the vane method under controlled shear-rate mode (from[Citation25]).

(3)

where, Tmax is the maximum torque; h is the height of the vane; d is the diameter circumscribed by the tip of the vane blades; and m is a constant. Two key assumptions are made: 1) the shear stress is uniformly distributed over the end surfaces of the vane; and 2) the material yields along a cylindrical surface outlined by the edges of the vane blades. If the first assumption is true, then m = 0. Steffe[Citation14] showed that errors can be up to nearly 15% if this assumption is violated.

An alternative solution to eliminate errors associated to the non-uniformity of stress distribution over the vane ends is to measure the response when the top of the vane blades are aligned evenly with the top surface of the sample (i.e., Z1 = 0, see ). Thus, the following equation can be used to calculate the yield stress:

(4)

In general, the assumption of m = 0 is usually accepted[Citation15]; thus one could rather measure the torque response when the blades are entirely immersed in the sample:

(5)

Steffe[Citation14, Citation15] suggested a procedure to evaluate the validity of m = 0 assumption by calculating the yield stress as the slope of Tmax vs. h curves, using vanes of same diameter but different lengths. With respect to the second assumption, the actual sheared surface can be up to 5% larger than the vane dimensions.[Citation1,Citation14] Barnes and Nguyen[Citation31] pointed out that yield stress values may vary up to 10% if diameter of the sheared surface were greater than the circumscribed circle defined by the blade tips. Steffe[Citation14] examined the results obtained using different vane dimensions (h and d) and container diameter (D) and recommended the following: 1.5 ≤ h/d ≤ 4.0; D/d > 2.0; Z1 = 0 or Z1/d > 1.0 if the vane is fully immersed in the sample; and Z2/d > 0.5. Genovesse and Rao[Citation11] used the vane technique to determine the static and dynamic yield stresses of six commercial foods. The procedure used was to measure the peak torque at undisrupted conditions (static yield stress) followed by applying a continuous cyclical rotation to disrupt the sample and measuring the peak torque afterwards (dynamic yield stress).

Recent and Improved Methods

Squeezing flow

Squeezing flow has been primarily used to measure the flow properties of highly viscous materials such as polymer melts. These flow arrangements were originally tested with industrial viscometers known as “parallel plate plastometers”.[Citation35] Since then, the technique has been applied for rheological measurements of a wide range of materials, including purely viscous liquids, yield stress fluids, viscoelastic solids and liquids, and purely elastic solids.[Citation36] Introduced to foods in the mid-1980s, squeezing flow has been a useful method to measure the rheological properties of semi-liquid and semi-solid food materials such as processed and melted cheeses, butter, dough, peanut butter, etc.[Citation25,Citation27]

The method has became very attractive for measuring rheological properties of very viscous materials and fluids with yield stress because it minimizes the two major problems encountered with conventional rotational viscometers: slip and extensive structural disruption upon insertion of sample into narrow viscometer gap.[Citation38] In addition, the squeeze flow geometry allows producing a wide range of strains and internal deformations of the material.[Citation36]

There are basically three test geometries ( and ): 1) specimens with a constant area and changing volume; 2) specimens with constant volume and changing area, and 3) “imperfect squeezing-flow.” In the first geometry, the specimen radius is usually equal to that of the plates, which guarantees a known contact area at all time and minimization of errors. However, as the sample flows and accumulates outside the plates during compression, an additional variable, pressure that is difficult to predict, is created.[Citation36]

Figure 11 Geometry of squeezing flow test under constant area and constant volume (from[Citation38]).

Figure 11 Geometry of squeezing flow test under constant area and constant volume (from[Citation38]).

Figure 12 Schematic of imperfect squeezing flow apparatus (from[Citation40]).

Figure 12 Schematic of imperfect squeezing flow apparatus (from[Citation40]).

In the constant volume geometry, the plates are larger than the specimen and the radial interface is free. Engmann et al.[Citation36] points out that stresses at the sample edges are more clearly defined, reducing additional pressure build-up at the edges. In imperfect squeezing flow, the bottom plate is replaced by a shallow container.[Citation39] However, not only the boundary conditions at the edges of the plates are unknown, the contact areas between plates and fluid of different size can cause an asymmetric flow.[Citation36]

The flow pattern during the squeezing of a fluid can be classified as being lubricated (perfect slip) and unlubricated (frictional or no slip). Lubrication of the parallel plates can be achieved intentionally (e.g., use of paraffin oil) or by the specimen itself (e.g., peanut butter and melted cheese). The flow pattern can be distinguished by noticing whether the expelled fluid is flat (lubricated) or rounded (non-lubricated). The adhesion and friction between the sample and the plates creates a barrier effect resulting in shear flow, and consequently, during frictional squeezing flow, squeezed specimen is both elongated and sheared. The flow behavior for lubricated, non-lubricated, and partially lubricated conditions have been theoretically analyzed and reviewed by Engmann et al.[Citation36]

Another factor that makes squeezing flow an attractive method is its simplicity. The technique can be performed in two forms: 1) constant or controlled deformation rate; and 2) constant force or stress. At constant deformation rate, the cylindrical specimen is compressed between two parallel plates at constant speed, whereby the force-time or force-height is measured, On the other hand, in constant stress, a creep test is performed and the response is a height-time relationship.

The determination of the yield stress from squeezing flow is based on the residual stress after shear and relaxation. The evaluation takes into consideration that a viscous-plastic fluid (i.e., Bingham plastic and Herschel-Bulkley fluid) when subjected to constant compressive stress flows until it reaches an equilibrium height (HL ) as shown in . The yield stress obtained from the squeezing flow rheometry is sometimes referred to as the apparent extensional yield stress [Citation41] or the dynamic yield stress.[Citation36]

Figure 13 Typical profile in squeeze flow testing for two types of fluids: Newtonian (full line) and Bingham (dashed line) (from[Citation36]).

Figure 13 Typical profile in squeeze flow testing for two types of fluids: Newtonian (full line) and Bingham (dashed line) (from[Citation36]).

Campanella and Peleg[Citation42] evaluated the yield stress of tomato ketchup, mustard and mayonnaise with a non-lubricated squeezing flow apparatus under constant force and constant area using the following analytical expression:

(6)

where, W is the load applied on the upper plate; HL the equilibrium height; and Ro is the initial radius of the cylindrical specimen. As for tests under lubricated squeezing flow (LSF), Engmann et al.[Citation36] points out that dynamic yield stresses cannot be obtained under constant area conditions due to the fact that the material will either continue to deform indefinitely or not deform at all.

Yang[Citation43] presented an expression for LSF of Herschel-Bulkley fluids relating the squeezing force and the fluid thickness. Based on this Ak and Gunasekaran[Citation41] derived an alternative analytical equation to determine the apparent extensional yield stress of a Herschel-Bulkley fluid under constant force and constant volume:

(7)

where, Ho is the initial height of the specimen. While yield stress measurements of semi-liquid foods by squeezing flow are often difficult to perform at constant deformation rate, an estimate of the yield stress can be obtained by determining the residual apparent yield stress with imperfect squeezing flow. Suwonsichon and Peleg[Citation40] evaluated this technique, in which the force exerted by the fluid (the residual force or stress) correlates to the yield stress value. For fluids without yield stress, the recorded force will drop to the buoyancy force instead to zero.

Moreover, an alternative to controlled deformation rate and imperfect squeezing flow geometries is the creep test. Sun[Citation21] measured the yield stress of concentrated spreadable foods using a modified squeezing-flow apparatus, known as the UW Meltmeter, which has been developed and used to measure the melt and flow behavior of cheeses by Gunasekaran and his co-workers.[Citation44] This device consists of one or two sample wells (7-mm deep and 30-mm in diameter) with an inner and outer ring, which makes it easier during sample loading. A stationary center piston allows the outer ring to move up and down through a lever, whereby when it is at the lower position, a cylindrical sample specimen is formed taking the shape and dimensions of the sample well ().

Figure 14 Picture of the dual sample well UW Meltmeter (Left) and sample formed with the lever in the down position (from[Citation21]).

Figure 14 Picture of the dual sample well UW Meltmeter (Left) and sample formed with the lever in the down position (from[Citation21]).

The UW Meltmeter performes squeezing flow test under constant force, constant volume and changing area. Store-bought peanut butter, margarine and mayonnaise were used in this study due to their textural and subjective spreadability (semi-hard, mild and easy, respectively). The “apparent extensional yield stress” of these selected spreadable foods were determined by applying the Herschel-Bulkley model proposed by Ak and Gunasekaran,[Citation41] shown in EquationEquation 7. The yield stress values measured with the UW Meltmeter and vane technique were in good agreement for all three spreadable food products. Overall, the UW Meltmeter is a simple, quick, and reliable technique to measure the yield stress parameter of a variety of semi-liquid foods, including spreadable foods.

Plate method and slotted-plate device

The plate method was developed to overcome some of the problems associated with the assumptions of the vane method that are usually overlooked[Citation19,Citation28]: sample yields along the cylinder described by the blades and uniform stress distribution along the blades. The plate method was specifically designed to produce a uniform stress distribution on the plate surface and to be applicable for both high and low concentration suspensions. De Kee et al.[Citation28] used this static method to determine the yield stress of several liquid and semi-liquid foods (condensed milk, corn syrup, molasses, ketchup, mayonnaise and tomato paste). The experimental procedures consist of immersing the plate into the glass container filled with sample material (). The platform is then lowered at controlled and constant speed (from 0.0127 to 2.54 cm/min). Initially, the yield stress fluid remains fixed to the plate and as the platform descends a force is exerted on the balance beam causing it to pivot.

Figure 15 Illustration of the plate method in measuring yield stress (from[Citation28]).

Figure 15 Illustration of the plate method in measuring yield stress (from[Citation28]).

This stage is represented between time zero and point A in the graphical response shown in . Once the fluid starts moving with respect to the plate, the force acting on the plate surface exceeds the yield stress (shown by the change in slope after point A), following an increase in force to a steady-state value for a given rate of descent (point C). Point B represents the stage in which the platform stops and the balance beam cannot return to its original horizontal position due to the yield stress. The yield stress is then determined as[Citation28]:

Figure 16 (A) Typical response of the plate method; and (B) yield stress obtained from extrapolation of measured shear stress versus platform speed (from[Citation28]).

Figure 16 (A) Typical response of the plate method; and (B) yield stress obtained from extrapolation of measured shear stress versus platform speed (from[Citation28]).

(9)

where, Fb is the buoyant force (g); Fr is the recording force at point B (g); g is the acceleration due to gravity (cm/s2); and s is the plate surface area from both sides and without the edges (cm2).

An alternative procedure is to extrapolate the yield stress value from a shear stress vs platform speed curve (). Although, the yield stress for the liquid and semi-liquid foods measured using the plate apparatus agreed closely with values determined from other methods, different suspensions and the type of plate used may cause discrepancy among the results, not to mention that yielding which may occur at the plate-sample interface as opposed to within the sample.[Citation19,Citation45]

As a result, slotted-plate device has been developed to combine the advantages of both the vane and the plate methods, while providing a uniform stress distribution at the plate surface and eliminating the wall effect.[Citation45] This new device is similar to the plate method apparatus with a thin stainless steel wire connecting the slotted plate with a balance (). A platform holds a double-walled beaker and a computer controls the speed of the platform. The slotted-plate device is housed in a Plexiglas box; two infrared lamps maintain constant air temperature.

Figure 17 (A) Schematic diagram of the slotted-plate device; and (B) typical time response curve (from[Citation45]).

Figure 17 (A) Schematic diagram of the slotted-plate device; and (B) typical time response curve (from[Citation45]).

The concept of the slotted-plate device lies in the fact that slots are filled once the slotted-plate is inserted into the test fluid. As the platform is lowered, the fluid exerts forces on the plate, which is continuously recorded by a balance. Flow begins once the forces exceed the yield stress, characterized by a change in slope in the response data (). Consequently, the apparent yield stress can be calculated as:

(10)

where, F is the measured force recorded by the balance; Fi is the initial force (weight of plate and wire without the buoyant force); Fr is the net force produced (difference between force F and Fi); and S is the total surface area of the plate (including those of the slots). Zhu et al.[Citation19] confirmed that the slotted-plate technique is reliable for yield stress measurement when compared to other indirect and direct methods and recommended using a slotted plate with slot ratio, the ratio between the slot area and the total plate area, greater than 0.5. The further applications of the plate and slotted-plate devices are still under development. In addition to the materials previously tested, fluid and semi-fluid foods,[Citation28] and mineral suspensions,[Citation19] the slotted-plate device has been used to measure the yield behavior of cohesionless powders.[Citation23]

Slump tests

The slump test has been developed as a quick and easy method to measure yield stress of thick, flocculated suspensions. Traditionally used to evaluate the “workability” of fresh concrete, Pashias et al.[Citation22] called it “a fifty-cent rheometer,” as it offers an inexpensive way to determine the yield stress. The test consists of filling a cylindrical mold with the test fluid and lifting the mold off to allow the material to collapse and flow under its own weight (). The slump height, which is the difference between the initial and final heights is measured. The analytical model relating yield stress to the slump height proposed by Murata[Citation46] was corrected by Christensen.[Citation47] Pashias et al.[Citation22] adapted it for a cylindrical geometry as follows:

Figure 18 Schematic diagram of the initial and final state of stress distribution (from[Citation22]).

Figure 18 Schematic diagram of the initial and final state of stress distribution (from[Citation22]).

(11)

where, is the dimensionless yield stress; and s′ is the dimensionless slump value, which is calculated as:

(12)

where, is the dimensionless form of the height of the non-yielded region; while is the dimensionless height of the yielded area ().

Figure 19 Schematic diagram of the slump test (from[Citation22]).

Figure 19 Schematic diagram of the slump test (from[Citation22]).

The standard method defined by the American Society for Testing and Materials for slump test is to use a frustum of a cone, in which the initial height is 30 cm.[Citation48,Citation49] With similar assumptions from Pashias et al.,[Citation22] a relation was developed for a conical geometry[Citation50]:

(13)

The theory behind the analytical model shown above is explained as follow: above a critical height along the undeformed cylinder, the material experiences a shear stress that is smaller than the yield stress, and thus no flow is observed; while below this point the stress created by the weight of the material above the critical point is larger than the yield stress. After the material has slumped, flow of the material takes place in the yielded region until the cross-sectional area increases to such an extent that the stress required to support the weight reaches the yield stress. Meanwhile, the unyielded region remains intact as when the material was underformed.[Citation22,Citation49] Clayton et al.[Citation51] found that yield stress measurements of mineral suspensions were more accurate by using a cylindrical model than the conical model.

From the studies thus far, it is possible to conclude that the slump test is a robust method of measuring yield stress. It was observed that slump height is in fact dependent on yield stress, once increasing yield stress values were found for decreasing slump heights. Studies have also shown that the method is not affected by the type of flocculated material, not to mention that it is sensitive enough to distinguish between a disturbed and undisturbed material based on the slump height. Furthermore, the velocity in which the frustum is lifted off and the surface of the material does not affect measurements.[Citation22]

A possible drawback lies in the range of yield stress values. For high yield stress materials (e.g., greater than 800 Pa), the difference between the initial and final heights is very minimal.[Citation22] However, this problem was overcome by placing additional weights on top of the sample.[Citation52] Omura and Steffe[Citation48] used a centrifugal slump test to predict the yield stresses of semi-solid foods with high yield stresses. In this test, a sample is allowed to slump after it has been subjected to centrifugal force by rotating it on a circular plane. A general agreement was found between yield values of cream cheese, peanut butter, butter, margarine and Neufchatel cheese using the centrifugal slump test and those measured by the vane method. As with the conventional slump test, the centrifugal slump test is also performed by measuring slumped height and original height of the material along with the acceleration of the sample. A simplified equation is given as function of dimensionless slump height[Citation48]:

(14)

where: z is slump height (m); H is original height of the sample (m);  = yield stress (Pa).

Tube viscometer and magnetic imaging resonance (MRI)

The determination of yield stress of fluids using a tube viscometer is done by inducing “plug flow” through a tube. In this method, the fluid flows at constant velocity through a tube of length L and radius R under a known pressure drop (ΔP). The yield stress is calculated from force balance on the fluid:

(15)

The flow is assumed to be steady and fully developed, not to mention absent of both end and wall effects.[Citation24] Thus, overestimation of the yield stress value are due to pressure losses as fluid enters from a reservoir and exits to open air resulting in end effects and wall effects due to apparent reduction in viscosity and dependence on tube dimensions in measuring fluid properties. In addition, slippage at the tube walls is also another major source of error which can be corrected by using a tube with rough walls.

The magnetic resonance imaging (MRI) has been integrated as part of the tube viscometer apparatus in order to ensure fully developed flow through MR velocity images and evaluate rheological parameters, including the yield stress. The MRI allows to directly image the velocity distribution across the tube and to measure the boundary conditions, so that assumptions of velocity at the wall can be avoided (e.g., slip or stick boundary).[Citation53] is a schematic diagram of the MRI-integrated tube viscometer.

Figure 20 Schematic diagram of the MRI apparatus integrated into tube viscometer (from[Citation54]).

Figure 20 Schematic diagram of the MRI apparatus integrated into tube viscometer (from[Citation54]).

Yoon and Mccarthy[Citation55] investigated the flow behavior of yogurt during pipe flow with MRI and found that the average yield stress and consistency index (K) were higher at 25°C (14 Pa and 5.6 Pa sn) than at 35°C (10 Pa and 1.1 Pa sn); whereas flow behavior indices were the opposite (0.36 and 0.88). They also studied the volumetric flow rates of three commercial processed cheeses during melting with the objective for further mixing and pumping of melted processed cheese.[Citation56]

APPLICATIONS OF YIELD STRESS IN FOODS

Concentrated Suspensions

Most foods are dispersions of solids in a liquid or one fluid in another. Concentrated suspensions are dispersions of solids suspended in an aqueous medium with concentration varying from 5 to 24% by weight. Yoo et al.[Citation57] measured the static and dynamic yield stress of structured and unstructured commercial food suspensions with the vane method at controlled shear rate and shear stress. The yield stress has also been used to evaluate the structural recovery of a disrupted sample to its original consistency.[Citation38]

Additionally, the role of particle size and particle concentration on the yield stress of dispersions has received considerable attention. In general, the magnitude of yield stress of food suspensions increases with increasing particle volume fraction, decreasing particle size, and increasing interparticle forces.[Citation58] Qiu and Rao[Citation59] demonstrated that the yield stress of apple sauces increased with increase in pulp content and with decrease in average particle size. However, Missaire et al.[Citation10] showed that differences between structured and unstructured suspensions can cause divergence of results. The yield stress is considered as a measure of interacting particles behavior and the microstructure. Several theoretical models have been attempted to relate yield stress with particle-particle interactions, volume fraction, and particle size.[Citation58]

Foams

Foams are widely used in food applications, either as finished products such as meringues and whipped toppings, or as inclusions during production such as in ice creams and batters for baked goods.[Citation60] Foams consist of a colloidal system containing discontinuous air phase in the form of tiny air bubbles dispersed in an aqueous continuous phase.[Citation61, Citation62] In general, foams are classified as either dry or wet, as well as according to their rheological behavior, solid-like or liquid-like. In both cases, the distinction is based on the gas phase volume (φ), as shown in .

Figure 21 (A) Structure of dry; and (B) wet foams (on Left). Schematic diagram of the solid-like and liquid-like behaviors of foams (on right) (from[Citation64]).

Figure 21 (A) Structure of dry; and (B) wet foams (on Left). Schematic diagram of the solid-like and liquid-like behaviors of foams (on right) (from[Citation64]).

At low φ’s (below approximately 70%), bubbles maintain their spherical structure, are called “bubbly foams” and behave as a viscous fluid; and if it is high, bubbles are pressed and packed together forming polyhedral foams and take the form of dry foams, thus displaying a solid-like behavior.[Citation61–65] Due to their lyophobic characteristics, foams are intrinsically unstable, and as a result, their structure tends to collapse within few seconds to days.[Citation62,Citation64] The addition of surfactants, such as proteins, to colloidal systems enhances foam stability. The main food proteins used are whey, egg white, and soy.

The most common methods to evaluate protein-based foams are measurements of overrun and stability. Solution properties such as ionic strength, pH, ion types and concomitants strongly influence foaming functionality and characteristics.[Citation61,Citation63] A common practice to examine foam formation is to determine how well it will hold up on a spatula or maintain its structure.[Citation14, Citation60] Thus, observing how well the foam resists flow under the influence of gravity is a subjective measure of of yield stress. Objectively, however, the yield stress of foams can be measured with conventional rotational rheometers such as parallel-plate, cup-and-bob, and cone-and-plate, and with capillary flow. Prud'homme and Khan[Citation66] obtained good rheological data by adding sandpaper to a parallel-plate rheometer to prevent wall slip. Similarly, the results from squeezing-flow tests were useful to monitor consistency of foams and to detect gum addition.[Citation38]

Pernell et al.[Citation60] demonstrated that the vane technique can be successfully employed in measuring the yield stress of protein-based foams made at various concentrations and whipping times. In addition, the point method (controlled rate mode) showed superior results compared to slope method (controlled stress mode). The yield stress values were similar for egg white protein (EWP) foams in both testing modes, but not for WPI foams. It was suggested that the latter create different properties when measured in different ways (i.e., difference in mixer models), while EWP foams were more stable and isotropic. Vane yield stresses for EWP foams were greater than those of WPI foams (100–150 Pa and 55–80 Pa, respectively), despite similar proteins volume fraction, bubble size, and surface tension. Based on this, many investigators had focused on determining factors that influence not only the stability and overrun of protein-based foams, but also their yield stresses.

Overall, yield stress values depend on the effects between the protein absorption at the surface and the protein distribution in the lamellae between air bubbles, rather than jus volume fraction or bubble size.[Citation62] Additionally, interfacial measurements of surface tension and interfacial dilatational elasticity suggest that variation in the yield stress among protein foams are related to the ability of the protein to form strong and elastic interfaces throughout the surface.[Citation65] Luck et al.[Citation63] observed that the yield stress of WPI foams is also dependent on pH, salt type, and salt concentration, since these factors alter the foam structure and behavior at the surface.

From the theoretical model of Princen and Kiss,[Citation67] the increase in yield stress of foams is proportional to the cubic root of its air phase volume (i.e., τo∝φ), as large φ’s result in tighter packing of the air bubbles. However, this relation is not entirely true when comparing foams made of different proteins. Although it is intuitive to directly correlate foam stability with yield stress, previous studies have shown that yield stress of foams may be inversely related to their stability. For instance, it is well known that adding sugars (e.g., lactose) to foams increases stability due to an increase in viscosity and decrease in the rate of liquid drainage[Citation61]; thus, yield stress and overrun of foams containing lactose decreases.[Citation63] Foams of hydrolysate proteins result in highest yield stress values and overrun abilities, though they present poor foaming properties.[Citation61,Citation63]

Mleko et al.[Citation68] demonstrated that pH-induced refolding of egg albumin produces substantially firmer foam than untreated egg albumin. Foam firmness was evaluated by the static and dynamic yield stress measured with steady shear experiments using a special cross-hatched geometry that prevents slippage. The authors pointed out that their findings can produce more stable and firmer products where egg whites are employed such as in meringue. In addition, pH treatments may improve overall performance and functionalities (e.g., texture) of protein-based foams in foods.

Concentrated Emulsions

The common measurements of emulsion stability are average particle size of emulsion droplets and evaluation of the thermodynamics of the system, especially at the interface of dispersed and continuous phases. Additionally, viscosity and yield stress are also provide information on how the spherical droplets are packed and interact among each other.[Citation69] While the maximum packing of uniform spheres consist of 74% of the space with the remaining 26% of the space is void, at volume concentrations between 60 and 70%, the contact between droplets affects sensible properties that may cause phase inversion (e.g., from O/W emulsion to W/O). Morrison and Ross[Citation69] pointed out that at concentrations greater than 70% by volume, emulsion takes the characteristics of polyhedral foams due to the packed structure features.

Moreover, the flow of concentrated emulsions is restricted by the interparticle interactions among the droplets. As the concentration increases (≥ 70%), droplets are tightly packed together and flow only occurs after a large shearing stress is applied to overcome the structure built up. Thus, the increase in both yield stress and apparent viscosity are associated with increasing particle volume fraction, decreasing droplet size, and increasing magnitude of interparticle forces.[Citation58,Citation69]

Mayonnaise is an oil-in-water emulsion (vegetable oil in aqueous solutions of mainly egg yolk, sugars and salts) with a volume fraction of approximately 65% in the dispersed phase. The yield stress relates to stability of mayonnaise in low-stress situations such as during storage and transportation. Ma and Barbosa-Cánovas[Citation70] reported that the yield stress of mayonnaise increased significantly with increase in added oil and xanthan gum, from 23 to 235 Pa and 55 to 195 Pa, respectively. The added oil promoted a more compact 3-D network formation between the egg protein molecules and absorbed droplets, while added xanthan gum increased the stability of mayonnaise emulsion through the formation of large aggregates. The structural stability of light and regular-fat content mayonnaise products has been compared by measuring their yield stress.[Citation71] The light mayonnaise exhibited higher yield stress, which indicated that the regular mayonnaise can be pumped more easily. The light mayonnaise requires an added artificial dispersant to prevent destabilization. Salad dressings are another common O/W emulsions with a minimum 30% volume fraction of the dispersed particles.[Citation69] Wendin and Hall[Citation72] investigated the effect of fat, thickener, and emulsifier content on salad dressings from creep-recovery tests. Higher yield stress values were obtained for samples containing high fat and thickener contents, and yield stress and viscosity parameters correlated well with sensory panel results.

Chocolate Coating

The yield stress of foods, such as melted chocolate, directly affects the coating thickness. For instance, if the yield stress is very low, only a small amount of coated layer will remain on the product, which may be undesirable to the consumer. On the other hand, if the material exhibits high yield stress, the coated layer may be too thick, which may be undesirable to the manufacturer. The coating thickness can be measured by studying the material flow on an inclined plane described earlier. Flow occurs when the shear stress exceeds the yield stress, while the fluid remains on the plane if yield stress is much higher than the maximum shear stress.[Citation15] As a result, the maximum coating thickness (h max) can be calculated as[Citation15]:

(16)

where, σ o and ρ are the yield stress and density of the coating material respectively; g is the acceleration due to gravity; and θ is the angle of the inclined plane. EquationEquation 15 assumes that the coating material does not flow off of the surface. Similarly, another expression for h max is given based on post-withdrawal drainage conditions, where θ = 90°[Citation15]:

(17)

In the chocolate manufacturing and confectionary industry, coating is performed with the enrobing technology, in which the molten chocolate flows in a sheet above a moving belt and the product is coated as it moves through the sheet (). The performance of enrobing depends both on the physics of coating flow and the rheological properties of chocolate.[Citation54] As for the latter, in addition to the yield stress, the apparent viscosity is also important.

Figure 22 Schematic of chocolate coating application with enrobing technology (from[Citation54]).

Figure 22 Schematic of chocolate coating application with enrobing technology (from[Citation54]).

The International Office of Cocoa, Chocolate, and Confectionary (IOCCC) has adopted the Casson model as a standard to describe the flow behavior of chocolate.[Citation15] The IOCCC recommends the use of rotational viscometer with a concentric cylinder probe for flow measurements; several studies have suggested additional recommendations such as measuring shear stresses at specific shear rates (e.g., 2, 5, 10, 20, and 50 s−1).[Citation54] Recognizing this possibility, Wichchukit et al.[Citation54] evaluated the rheological properties of molten chocolate during steady pipe flow using a magnetic resonance imaging (MRI) tube viscometric method, which provides velocity profile and flow behavior during pipeline flow. The yield stress values from the MRI tube viscometry were compared to measurements from rotational viscometer as recommended by the IOCCC and from the vane method. The researchers measured the yield stress values as a function of emulsifier concentration.

The Casson model appeared to fit well the experimental data for all three measurement methods. The MRI images showed that the velocity profile was maximum at the center of the tube, which is expected. The maximum velocity was reached with increasing amount of emulsifiers, whereas both yield stress (15 to 1.9 Pa) and apparent viscosity (14.6 to 6 Pa-s) decreased. Comparing the MRI results to those from rotational viscometer, the latter were lower due to possible wall slip effects, while vane results were in good agreement. Wichchukit et al. [Citation54] suggested that the tube viscometer measures static yield stress, as the data correspond to the shear stress to initiate flow. Another aspect of their work was to simulate the average film thickness of milk chocolate during enrobing from an unsteady state mass balance. The numerical solutions were in agreement with experimental results, and as emulsifier level increased from 0 to 0.3%, the predicted thickness decreased by 60%.

Spreadable Foods

Spreadable foods such as butter, margarine, peanut butter, and cream cheese are elastoplastic or viscoplastic and exhibit yield stress. Early studies by Kokini and Dickie[Citation73] on psychophysical mechanisms and fluid mechanics suggested that subjective spreadability of foods could be assessed through the back and forth motion of a knife, and consequently, by estimating the shear stress on the knife surface. This relation may be illustrated as the reciprocal of the maximum shear stress required to create a uniform distribution over a surface[Citation74]:

(18)

In fact, several studies indicate that yield stress is a better measure of spreadability of semi-solid foods than apparent viscosity.[Citation33,Citation41,Citation75] According to Daubert et al.[Citation74], yield stress can be associated with quality and consumer acceptance of spreadable foods, and it provides additional parameter for quality control. A variety of products have been tested and their spreadability values were presented in yield stress-yield strain texture maps ().[Citation33,Citation74] Though yield stress is the major contributor to the spreadability of soft foods, the yield strain provides information relative to the extent of deformation and uniformity upon yielding (i.e., flowing). Consequently, while low yield stress indicates good spreadability, the ability of a soft food to resist deformation at low strains can result in poor spread uniformity.[Citation33,Citation74]

Figure 23 Spreadability map for elastoplastic foods: (M) margarine spread; (PCS) processed cheese spread; (PB) peanut butter; (RFPB) reduced fat peanut butter; (GJ) grape jelly; (CC) cream cheese; (FCC) free cream cheese; (B) touch of butter; (FY) fat free plain yogurt; (SC) sour cream; (NSC) non-fat sour cream; (WT) whipped cream; and (FWT) fat free whipped cream (from[Citation74]).

Figure 23 Spreadability map for elastoplastic foods: (M) margarine spread; (PCS) processed cheese spread; (PB) peanut butter; (RFPB) reduced fat peanut butter; (GJ) grape jelly; (CC) cream cheese; (FCC) free cream cheese; (B) touch of butter; (FY) fat free plain yogurt; (SC) sour cream; (NSC) non-fat sour cream; (WT) whipped cream; and (FWT) fat free whipped cream (from[Citation74]).

Peanut butter is a classic viscoplastic food, which consists of a thick concentrated suspension of small non-colloidal peanut particles in peanut oil[Citation38]; other vegetable oils and stabilizers may be added into the oil matrix to disperse solid particles. Citerne et al.[Citation76] examined the rheological and physical characteristics of two types of peanut butter. One of the products (“100% peanuts”) was an unstabilized suspension with solids and particles in peanut oil, and the second product, referred to as “smooth,” was stabilized with a vegetable oil and contained additives such as salts and sugars. The “apparent yield stress” of both types was obtained by fitting the flow curves, which was measured with parallel plates geometry attached to stress-controlled rheometer, to Casson and Bingham models.

The presence of additives was clearly reflected on the measured yield stress data. The “100% peanut butter” exhibited much smaller yield values, 21.5 and 27 Pa, than the “smooth” type, 363 and 374 Pa, based on Casson and Bingham models, respectively. Similar results were obtained from creep experiments. Transition stress values, denoted by a rapidly increasing viscosity, of 250 Pa and 10 Pa were reported for “smooth” and “100% peanut butter” samples, respectively. Below the transition stress, the sample behaved as viscoelastic solid, and at or above the transition stress as liquid-like material. In addition, Citerne et al.[Citation76] observed, from time-dependent investigation with dynamic oscillatory tests in the non-linear region, an elastic plateau at stresses near the apparent yield stress.

Cheese

Truong and Daubert[Citation77] used vane rheometry as an alternative method for torsion and uniaxial compression tests to evaluate the texture of different cheeses (from hard Cheddar to soft processed cheese). A four-bladed vane was employed to obtain data of stress and deformation from the typical torque-time response at gel fracture. The vane stress was calculated using the relation originally proposed by Dzuy and Boger[Citation2,Citation32] shown in EquationEquation 3, and was used to describe the textural properties of cheeses. Meanwhile, the angular deformation was related to the strain characteristic of the material such as in the torsion test.

Comparing the results from the vane and torsion tests, shear stresses obtained from the vane test were lower than those from the torsion test. Similarly, the vane angular deformation at failure were lower than fracture shear strain values in torsion, but both tests showed Mozzarella cheese the most elastic cheese type with highest deformation. Truong and Daubert[Citation77] recommended the vane technique for stress and deformation measurements of a wide range of cheeses, especially because it overcomes usual problems with torsion and uniaxial compression tests, which are difficult to perform on soft or sticky and elastic materials, respectively.

Furthermore, the yield stress measured with vane has been used to distinguish ingredients used in cheese. Mleko and Foededing[Citation78,Citation79] measured yield stress of rennet casein gels and processed cheese analogs containing various amounts of whey protein isolate (WPI). The yield stress of the analogs increased with addition of WPI. Additionally, substituting part of the rennet casein in processed cheese analogs by whey protein through heated WPI solutions, resulted in stronger cheese with higher yield values.

Ice Cream

Ice cream manufacturers constantly modify and develop new products to meet consumers demand such as for reduced-fat ice creams, flavors, and creaminess levels. Briggs et al.[Citation80] investigated the viability of testing ice cream with the vane technique and measured the “scoopability” of ice cream based on the yield stress at typical scooping temperatures. They showed that the vane yield stress can be potentially used as a quality control tool, as the yield stress is strongly correlated to body, texture, and scoopability. Yield stress values were also used to identify differences in composition. For example, a less expensive product of brand “B” with lower amount of total solids and higher overrun exhibited higher yield stress than a more expensive product of brand “A.” Furthermore, the correlation between temperature and yield tress was stronger for brand “A” than for brand “B.”

The yield stress was also used as a parameter to study the rheological and structural properties of ice cream mix containing stabilizers locust bean gum (LBG) and guar gum (GG) during sub-zero temperature fluctuations.[Citation81] These stabilizers are commonly added to slow down ice crystal growth, especially during recrystallization. The influence of other ingredients (milk solids non-fat (MSNF), sucrose, emulsifiers) was also taken into consideration so that a more realistic system could be mimicked. The yield stress was measured by extrapolation of flow curve data.

Patmore et al.[Citation81] indicated that LBG solutions developed viscoelastic weak gels with temperature cycling, especially in the presence of MSNF, as opposed to GG solutions. The formation of weak gel networks around ice crystals reduced recrystallization rates, because the networks formed entraps melted water and prevented it from freezing onto other ice crystals. The addition of MSNF increased the yield stress overall, though did not cause GG to gel. Meanwhile, fat droplets seemed to interfere with the intermolecular interactions in MSNF/LBG emulsions, which were shown by a decrease in the yield stress.

Yogurts

The yield stress in yogurts is associated to the structural characteristics of yogurt gels. In some types of yogurts, such as fruit-in-the-bottom set types, a considerable firmness is desired to prevent the disruption and breakdown of the structure, whey separation, and better overall texture. Harte et al.[Citation82] examined the initial firmness of yogurt based on measurements of yield stress with the vane technique and responses from trained sensory panelists. Yield stress significantly correlated with sensory initial firmness for retail yogurts and laboratory-made yogurts with different amounts of added gum. Yield stress measurements and sensory responses showed that increasing the gum concentration increased initial yogurt firmness.

Yield stress was also used as a parameter to evaluate the characteristics of yogurts as a function of manufacturing process. Harte et al.[Citation83] identified differences on the gel firmness for yogurts made from raw milk and from milk subjected to thermal processing (85°C for 30 min), or high hydrostatic pressure process (193 or 676 MPa for 5 or 30 min). Yogurt gels made from milk treated at higher-pressure for 30 min and thermally processed had similar yield stress values, and both formed firmer gels in contrast to yogurts made from milk at lower-pressure and not thermally processed. In fact, the yield stress values of yoghurts made using milk from the latter two processes were similar. Notably, the effect of variations in milk treatment on the rearrangement of milk proteins and the yogurt gel microstructure was reflected in the yield stress.

In addition to milk processing, the effects of inoculation rate and incubation, not to mention the combination of all three process stages, affect the structure and physical properties of yogurt gels. Lee and Lucey[Citation84] examined the effects on the initial gel strength and structural breakdown of intact and stirred yogurts made at different pre-heating temperatures and incubation temperatures. The yogurt gels made from high heating temperature (e.g., 85°C) and low incubation temperatures (e.g., 32°C) resulted in firmer gels with higher apparent viscosity and yield stress values. Low yield stress values of intact yogurt gels made with low milk heating temperature (75°C for 30 min) and high incubation temperature (44°C) indicated that these gels had weaker networks. Moreover, stirred yogurt gels had much less resistance than intact gel, as it would be expected. Thus, higher heating temperature and lower incubation temperature result in higher overall attributes in stirred yogurts.

Similarly, Lee and Lucey[Citation85] investigated the effects on the structural network of yogurt gels made at various inoculation rates and incubation temperatures. Yogurt gels made with 0.5% (wt/wt) inoculation rate of starter culture were less stiff than yohurt made with higher inoculation rates (e.g., 2 to 4%). Overall, combinations of intermediate to high inoculation rates along with lower incubation temperatures increase the yield stress of yogurt gels as highly cross-linked and uniformly distributed protein networks are formed.

On the other hand, at lower inoculation rates (e.g., 0.5%), solubilization of colloidal calcium phosphate occurs earlier and at higher pH; thus, the arrangement and interaction of casein micelles in the gel network are changed. Yogurt gels made with lower inoculation rates had lower yield stress values due to weaker interaction between casein particles, causing rearrangement and whey separation more likely to occur.

Hot Cereals and Starch Dispersions

The yield stress values measured by the vane method was used to characterize the flow behavior of wheat and oat hot breakfast cereals.[Citation86] Mimicking the influence of thermal and chemical effects during consumer use, long holding periods of cereals in hot water were determined to increase yield stress, which is caused by the hydration and gelatinization of starch in the cereals. When hot water is added to the mix, the starch granules absorb water and swell. As starch dispersions are gelatinized, the amylose molecules solubilize and leach out of the granules.[Citation86] The swollen starch granules become dispersed in a continuous amylose network.[Citation17] The vane yield stress values were used as a measure to evaluate the strength of network bonds of cross-linked waxy maize (CWM), tapioca, and Amioca starch dispersions.[Citation17] After each of the dispersions has been prepared, sample was let to equilibrate for one hour to recover its structure and reach an equilibrium temperature.

The static yield stress was measured immediately, and the dynamic yield stress after a cyclical flow test had been performed to disrupt the samples' structure. As expected, dynamic yield stress values were lower than the static yield stress values. From the typical torque versus time response, sharpness at the peak stress was more induced with increasing rotation speed (0.05 versus 2.07 rpm). The authors suggested that at very low shear rates, the vane yield stress would be closer to the “true yield stress.” In addition, at low shear rates (e.g., 0.4 s−1), the contribution of viscous drag to the total stress was very minimal.

The stress required to break the network structure decreased in the order of cross-linked waxy maize (CWM), followed by tapioca and Amioca. The authors explained using a texture map that the breaking point for CWM dispersion took place at much lower angular deformations, indicating a brittle behavior due to the amylopectin cross-linking present, which was responsible for restricted swelling. Meanwhile, both tapioca and Amioca behaved more rubber-like once they swelled, respectively, very well and moderately well.

CONCLUSION

Yield stress is a key parameter in various food manufacturing processes and an important tool for product development and quality control. It has been shown that the magnitude of yield stress is strongly correlated to sensory attributes like firmness of yogurt gels, spreadability characteristics, mouthfeel, among many others. Additionally, variation of ingredient sources and texture/structure changes can be detected by yield stress measurement, as yield stress can quantify the level of interacting particles and the microstructure governed by the interparticle network links. Several techniques have been proposed to determine yield stress and many others are being developed and improved continually. Techniques are often classified either as indirect or direct methods, but they may also be categorized as imitative, empirical or fundamental tests. There is no single or “best” method for different foods, and it is not unusual to obtain vastly different yield values obtained for a given product from different techniques. The most appropriate method to measure yield stress depends not only on the availability of the required equipment, but also on the applicability of the measured values.

REFERENCES

  • Liddell , P.V. and Boger , D.V. 1996 . Yield Stress Measurements with the Vane. , 63 ( 2–3 ) : 235 – 261 .
  • Barnes , H.A. and Walters , K. 1985 . The Yield Stress Myth? , V24 ( 4 ) : 323 – 326 .
  • Evans , I.D. 1992 . Letter to the Editor: On the Nature of the Yield Stress . 36 ( 7 ) : 1313 – 1318 .
  • Astarita , G. 1990 . Letter to the Editor: The Engineering Reality of the Yield Stress . 34 ( 2 ) : 275 – 277 .
  • Schurz , J. 1992 . Letter to the Editor: A Yield Value in a True Solution? . 36 ( 7 ) : 1319 – 1321 .
  • Spaans , R.D. and Williams , M.C. 1995 . Letter to the Editor: At Last, a True Liquid-Phase Yield Stress. , 39 ( 1 ) : 241 – 246 .
  • Hartnett , J.P. and Hu , R.Y.Z. 1989 . Technical Note: The Yield Stress—an Engineering Reality. , 33 ( 4 ) : 671 – 679 .
  • Dzuy , N.Q. and Boger , D.V. 1983 . Yield Stress Measurement for Concentrated Suspensions. , 27 ( 4 ) : 321 – 349 .
  • Barnes , H.A. 1999 . The Yield Stress–a Review or `[Pi][Alpha][Nu][Tau][Alpha] [Rho][Var Epsilon][Iota]’–Everything Flows? , 81 ( 1–2 ) : 133 – 178 .
  • Missaire , F. , Qiu , C.G. and Rao , M.A. 1990 . Yield Stress of Structured and Unstructured Food Suspensions. , 21 ( 4 ) : 479 – 490 .
  • Genovese , D.B. and Rao , M.A. 2005 . Components of Vane Yield Stress of Structured Food Dispersions . 70 ( 8 ) : e498 – e504 .
  • Cheng , D.C.H. 1986 . Yield Stress: A Time-Dependent Property and How to Measure It . V25 ( 5 ) : 542 – 554 .
  • Keentok , M. 1982 . The Measurement of the Yield Stress of Liquids . V21 ( 3 ) : 325 – 332 .
  • Steffe , J.F. 1992 . “ Yield Stress: Phenomena and Measurement ” . In Advances in Food Engineering , Edited by: Singh , R.P. and Wirakartakusumah , M.A. 363 London : CRC Press .
  • Steffe , J.F. 1996 . Rheological Methods in Food Process Engineering Freeman Press , MI : East Lansing .
  • Yoo , B. and Rao , M.A. 1994 . Effect of Unimodal Particle Size and Pulp Content on Rheological Properties of Tomato Puree. , 25 ( 4 ) : 421 – 436 .
  • Genovese , D.B. and Rao , M.A. 2003 . Vane Yield Stress of Starch Dispersions. , 68 ( 7 ) : 2295 – 2301 .
  • Barbosa-Cánovas , G.V. , Kokini , J.L. , Ma , L. and Ibarz , A. 1996 . The Rheology of Semiliquid Foods . 39 : 1 – 69 .
  • Zhu , L. , Sun , N. , Papadopoulos , K. and De Kee , D. 2001 . A Slotted Plate Device for Measuring Static Yield Stress. , 45 ( 5 ) : 1105 – 1122 .
  • Yoshimura , A.S. , Prud'homme , R.K. , Princen , H.M. and Kiss , A.D. 1987 . A Comparison of Techniques for Measuring Yield Stresses. , 31 ( 8 ) : 699 – 710 .
  • Sun , A. 2007 . Measuring the Yield Stress of Spreadable Foods with the Uw Meltmeter - a Modified Squeezing-Flow Apparatus Thesis , Madison : Type University of Wisconsin .
  • Pashias , N. , Boger , D.V. , Summers , J. and Glenister , D.J. 1996 . A Fifty Cent Rheometer for Yield Stress Measurement. , 40 ( 6 ) : 1179 – 1189 .
  • Zhu , H. and De Kee , D. 2006 . Slotted-Plate Device to Measure the Yield Behavior of Powders: Fem Analysis. , 45 ( 15 ) : 5375 – 5383 .
  • Nguyen , Q.D. and Boger , D.V. 1992 . Measuring the Flow Properties of Yield Stress Fluids. , 24 ( 1 ) : 47 – 88 .
  • Gunasekaran , S. and Ak , M.M. 2003 . Cheese Rheology and Texture , New York : CRC Press .
  • Park , N.A. , Cho , Y.I. and Irvine , T.F. 1990 . Steady Shear Viscosity Measurements of Viscoelastic Fluids with the Falling Needle Viscometer. , 34 ( 3 ) : 351 – 357 .
  • Uhlherr , P.H.T. , Guo , J. , Fang , T.-N. and Tiu , C. 2002 . Static Measurement of Yield Stress Using a Cylindrical Penetrometer. , 14 ( 1 ) : 17 – 23 .
  • De Kee , D. , Turcotte , G. , Fildey , K. and Harrison , B. 1980 . New Method for the Determination of Yield Stress. , 10 ( 3 ) : 281 – 288 .
  • De Kee , D. , Schlesinger , M. and Godo , M. 1988 . Postwithdrawal Drainage of Different Types of Fluids. , 43 ( 7 ) : 1603 – 1614 .
  • Lang , E.R.B. and Rha , C. 1981 . Determination of the Yield Stress of Hydrocolloid Dispersions. , 12 ( 1 ) : 47 – 62 .
  • Barnes , H.A. and Nguyen , Q.D. 2001 . Rotating Vane Rheometry – a Review. , 98 ( 1 ) : 1 – 14 .
  • Dzuy , N.Q. and Boger , D.V. 1985 . Direct Yield Stress Measurement with the Vane Method. , 29 ( 3 ) : 335 – 347 .
  • Breidinger , S.L. and Steffe , J.F. 2001 . Texture Map of Cream Cheese. , 66 ( 3 ) : 453 – 456 .
  • Rao , M.A. and Steffe , J.F. 1997 . Measuring Yield Stress of Fluid Foods , 51 ( 2 ) : 50 – 52 .
  • Leider , P. and Bird , R.I. 1974 . Theoretical Analysis . Squeezing Flow between Parallel Disks , 13 ( 4 ) : 336 – 341 .
  • Engmann , J. , Servais , C. and Burbidge , A.S. 2005 . Squeeze Flow Theory and Applications to Rheometry: A Review. , 132 ( 1–3 ) : 1 – 27 .
  • 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. , 16 ( 3 ) : 281 – 301 .
  • Corradini , M.G. and Peleg , M. 2005 . Consistency of Dispersed Food Systems and Its Evaluation by Squeezing Flow Viscometry. , 36 ( 5–6 ) : 605 – 629 .
  • Lee , S.-J. and Peleg , M. 1992 . Imperfect Squeezing Flow Viscosimetry with a Wide Plate and a Shallow Container. , 23 ( 3 ) : 267 – 278 .
  • Suwonsichon , T. and Peleg , M. 1999 . Imperfect Squeezing Flow Viscometry of Mustards with Suspended Particulates. , 39 ( 2 ) : 217 – 226 .
  • Ak , M.M. and Gunasekaran , S. 2000 . Simulation of Lubricated Squeezing Flow of a Herschel-Bulkley Fluid. , 10 ( 6 ) : 274 – 279 .
  • Campanella , O.H. and Peleg , M. 1987 . Determination of the Yield Stress of Semi-Liquid Foods from Squeezing Flow Data. , 52 ( 1 ) : 214 – 215 .
  • Yang , F. 1998 . Exact Solution for Compressive Flow of Viscoplastic Fluids under Perfect Slip Wall Boundary Conditions . V37 ( 1 ) : 68 – 72 .
  • Wang , Y.C. , Muthukumarappan , K. , Ak , M.M. and Gunasekaran , S. 1998 . A Device for Evaluating Melt/Flow Characteristics of Cheeses. , 29 ( 1 ) : 43 – 55 .
  • Zhu , L. , Papadopoulos , K.D. and De Kee , D. 2004 . Slotted-Plate Apparatus to Measure Yield Stress of Suspensions. , 75 ( 3 ) : 737 – 740 .
  • Murata , J. 1984 . Flow and Deformation of Fresh Concrete. , V17 ( 2 ) : 117 – 129 .
  • Christensen , G. 1991 . Modelling the Flow of Fresh Concrete: The Slump Test. Thesis , Type Princeton University .
  • Omura , A.P. and Steffe , J.F. 2001 . Centrifugal Slump Test to Measure Yield Stress. , 66 ( 1 ) : 137 – 140 .
  • Roussel , N. and Coussot , P. 2005 . “Fifty-Cent Rheometer” for Yield Stress Measurements: From Slump to Spreading Flow , 49 ( 3 ) : 705 – 718 .
  • Schowalter , W.R. and Christensen , G. 1998 . Toward a Rationalization of the Slump Test for Fresh Concrete: Comparisons of Calculations and Experiments. , 42 : 865 – 870 .
  • Clayton , S. , Grice , T.G. and Boger , D.V. 2003 . Analysis of the Slump Test for on-Site Yield Stress Measurement of Mineral Suspensions. , 70 ( 1–4 ) : 3 – 21 .
  • Baudez , J.-C. , Chabot , C. and Coussot , P. 2002 . Rheological Interpretation of the Slump Test. , 12 : 133 – 141 .
  • Hills , B. 1998 . Magnetic Resonance Imaging in Food Science , John Wiley & Sons .
  • Wichchukit , S. , Mccarthy , M.J. and Mccarthy , K.L. 2005 . Flow Behavior of Milk Chocolate Melt and the Application to Coating Flow . 70 ( 3 ) : E165 – E171 .
  • Yoon , W.B. and Mccarthy , K.L. 2002 . Rheology of Yogurt During Pipe Flow as Characterized by Magnetic Resonance Imaging. , 33 ( 5 ) : 431 – 444 .
  • Yoon , W.B. and Mccarthy , K.L. 2003 . Flow Behavior of Processed Cheese Melts. , 26 ( 6 ) : 559 – 576 .
  • Yoo , B. , Rao , M.A. and Steffe , J.F. 1995 . Yield Stress of Food Dispersions with the Vane Method at Controlled Shear Rate and Shear Stress. , 26 ( 1 ) : 1 – 10 .
  • Genovese , D.B. , Lozano , J.E. and Rao , M.A. 2007 . The Rheology of Colloidal and Noncolloidal Food Dispersions . 72 ( 2 ) : R11 – R20 .
  • Qiu , C.G. and Rao , M.A. 1988 . Role of Pulp Content and Particle Size in Yield Stress of Apple Sauce. , 53 ( 4 ) : 1165 – 1170 .
  • Pernell , C.W. , Foegeding , E.A. and Daubert , C.R. 2000 . Measurement of the Yield Stress of Protein Foams by Vane Rheometry. , 65 ( 1 ) : 110 – 114 .
  • Damodaran , S. 2005 . Protein Stabilization of Emulsions and Foams . 70 ( 3 ) : R54 – R66 .
  • Pernell , C.W. , Foegeding , E.A. , Luck , P.J. and Davis , J.P. 2002 . Properties of Whey and Egg White Protein Foams. , 204 ( 1–3 ) : 9 – 21 .
  • Luck , P.J. , Bray , N. and Foegeding , E.A. 2002 . Factors Determining Yield Stress and Overrun of Whey Protein Foams. , 67 ( 5 ) : 1677 – 1681 .
  • Höhler , R. and Cohen-Addad , S. 2005 . Rheology of Liquid Foam . 17 ( 41 ) : R1041
  • Davis , J.P. , Foegeding , E.A. and Hansen , F.K. 2004 . Electrostatic Effects on the Yield Stress of Whey Protein Isolate Foams. , 34 ( 1 ) : 13 – 23 .
  • Prud'homme , R.K. and Khan , S.A. 1996 . “ Experimental Results on Foam Rheology ” . In Foams , Edited by: Prud'homme , R.K. and Khan , S.A. 217 – 241 . New York : Marcell-Dekker Inc .
  • Princen , H.M. and Kiss , A.D. 1989 . Rheology of Foams and Highly Concentrated Emulsions : Iv . An Experimental Study of the Shear Viscosity and Yield Stress of Concentrated Emulsions. , 128 ( 1 ) : 176 – 187 .
  • Mleko , S. , Kristinsson , H.G. , Liang , Y. and Gustaw , W. 2007 . Rheological Properties of Foams Generated from Egg Albumin after Ph Treatment. , 40 ( 5 ) : 908 – 914 .
  • Morrison , I.D. and Ross , S. 2002 . “ Emulsions ” . In Colloidal Dispersions : Suspensions, Emulsions, and Foams , 420 New York : Wiley-Interscience .
  • Ma , L. and Barbosa-Cánovas , G.V. 1995 . Rheological Characterization of Mayonnaise . Part Ii: Flow and Viscoelastic Properties at Different Oil and Xanthan Gum Concentrations , 25 ( 3 ) : 409 – 425 .
  • Tabilo-Munizaga , G. and Barbosa-Cánovas , G.V. 2005 . Rheology for the Food Industry. , 67 ( 1–2 ) : 147 – 156 .
  • Wendin , K. and Hall , G. 2001 . Influences of Fat, Thickener and Emulsifier Contents on Salad Dressing: Static and Dynamic Sensory and Rheological Analyses. , 34 ( 4 ) : 222 – 233 .
  • Kokini , J.L. and Dickie , A.N.N. 1982 . A Mode of Food Spreadability from Fluid Mechanics. , 13 ( 2 ) : 211 – 227 .
  • Daubert , C.R. , Tkachuk , J.A. and Truong , V.D. 1998 . Quantitative Measurement of Food Spreadability Using the Vane Method. , 29 ( 4 ) : 427 – 435 .
  • Mortensen , B. and Danmark , H. 1982 . Consistency Characteristics of Butter . 27 : 530 – 532 .
  • Citerne , G.P. , Carreau , P.J. and Moan , M. 2001 . Rheological Properties of Peanut Butter. , V40 ( 1 ) : 86 – 96 .
  • Truong , V.D. and Daubert , C.R. 2001 . Textural Characterization of Cheeses Using Vane Rheometry and Torsion Analysis. , 66 ( 5 ) : 716 – 721 .
  • Mleko , S. and Foegeding , E.A. 2000 . Physical Properties of Rennet Casein Gels and Processed Cheese Analogs Containing Whey Proteins. , 55 ( 9 ) : 513 – 516 .
  • Mleko , S. and Foegeding , E.A. 2001 . Incorporation of Polymerized Whey Proteins into Processed Cheese Analogs. , 56 ( 11 ) : 612 – 615 .
  • Briggs , J.L. , Steffe , J.F. and Ustunol , Z. 1996 . Vane Method to Evaluate the Yield Stress of Frozen Ice Cream. , 79 ( 4 ) : 527 – 531 .
  • Patmore , J.V. , Goff , H.D. and Fernandes , S. 2003 . Cryo-Gelation of Galactomannans in Ice Cream Model Systems. , 17 ( 2 ) : 161 – 169 .
  • Harte , F. , Clark , S. and Barbosa-Cánovas , G.V. 2007 . Yield Stress for Initial Firmness Determination on Yogurt. , 80 ( 3 ) : 990 – 995 .
  • Harte , F. , Amonte , M. , Luedecke , L. , Swanson , B.G. and Barbosa-Cánovas , G.V. 2002 . Yield Stress and Microstructure of Set Yogurt Made from High Hydrostatic Pressure-Treated Full Fat Milk. , 67 ( 6 ) : 2245 – 2250 .
  • Lee , W.J. and Lucey , J.A. 2006 . Impact of Gelation Conditions and Structural Breakdown on the Physical and Sensory Properties of Stirred Yogurts. , 89 ( 7 ) : 2374 – 2385 .
  • Lee , W.J. and Lucey , J.A. 2004 . Structure and Physical Properties of Yogurt Gels: Effect of Inoculation Rate and Incubation Temperature. , 87 ( 10 ) : 3153 – 3164 .
  • Dejongh , J.J. and Steffe , J.F. 2004 . Yield Stress of Hot Cereals by the Vane Method. , 35 ( 1 ) : 1 – 9 .

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.