|
|
||||||||
Department of Chemistry, University College Dublin, Belfield, Dublin 4
Corresponding author: V. Buckin; e-mail: buckinv{at}eircom.net.
| ABSTRACT |
|---|
|
|
|---|
Key Words: high resolution ultrasonic spectroscopy concentrated milk heat stability ultrasonic velocity
Abbreviation key: HCT = heat coagulation time, HRUS = high resolution ultrasonic spectroscopy, REM = recombined evaporated milk, SMP = skim milk powder, SMUF = simulated milk ultrafiltrate, WPI = whey protein isolate
| INTRODUCTION |
|---|
|
|
|---|
Numerous studies have been carried out in the field of assessment of the heat stability of milks and their ingredients. Sommer and Hart (1919) developed a heat stability test as early as 1919. Various tests, similar in principle, have been subsequently devised (Sommer and Hart, 1922; Webb and Holm, 1932; Pyne and McHenry, 1955; Rose and Tessier, 1959; Belec and Jenness, 1962) until the 1960s when Davies and White (1966) published their subjective method that forms the basis of the 2 tests prevailing today: the Irish Dairy Board method and the Standard Association of Australia method. Both methods are based on the visual determination of heat coagulation of a milk sample at 120°C. In the Irish method (Irish Dairy Board, 1974), a small milk sample (typically 2.5 mL) is placed in a rocking sealed glass tube immersed in an oil bath maintained at constant temperature; the time required for particles of coagulated protein to become visible is recorded. Assay temperature, agitation rate, degree of tube fill, and angle of tilt influence the coagulation time determined by this method (Davies and White, 1966). The Australian method (Standard Association of Australia, 1978) uses bigger samples (30 mL) and operates on a discontinuous principle. The bottles containing the milk samples are taken out of the oil bath every 3 min, shaken to check for flocculation, and then returned to the heating bath.
Other techniques have been suggested in the past for the determination of heat coagulation times (HCT). White and Davies (1966) determined HCT by following the percentage of total nitrogen remaining in dispersion after centrifugation. The use of this technique is favored if a detailed description of the heat coagulation process is required. Foissy and Kneifel (1984) developed an automatic method using an electromagnetic sensor device, and de Wit et al. (1986) measured HCT by increases in viscosity using a computer-linked, falling-ball viscometer called the Klaro-Graph. These automated techniques can be left unsupervised, but are destructive. Therefore, they have not been widely adopted in studies of the heat stability of milk.
In this work, we describe the application of novel high resolution ultrasonic spectroscopy (HRUS) for the study of heat stability of 2 REM samples submitted to different preheat treatments and of milk ingredients such as casein and whey proteins. High resolution ultrasonic spectroscopy is a novel analytical technique with high potential for analyses of a wide range of samples and molecular processes (Buckin and ODriscoll, 2002). It presents various advantages such as being fast, noninvasive, highly sensitive, and able to perform measurements in optically opaque systems. It also allows continuous monitoring of the changes happening in milk while submitted to a high temperature, which is difficult using other techniques. This technique is based on the precise measurements of parameters of ultrasonic waves propagating through the samples. Ultrasonic parameters are extremely sensitive to the molecular organization and molecular interactions in the samples and can be exploited in the analysis of a broad range of molecular processes (Povey et al., 1999; Kudryashov et al., 2000; Sigfusson et al., 2004). Until now, the use of ultrasonic spectroscopy for material analysis was restrained to certain areas because of a limited resolution and the need for large sample volumes (
10 mL). However, it is now possible to measure ultrasonic parameters to a very high resolution (down to 0.2 mm/ s for ultrasonic velocity, which is up to 1000 times better than traditional ultrasonic techniques) in small volumes (down to 0.03 mL, 1 mL typically) and a broad range of sample types with modern commercial HRUS spectrometers. Monitoring both the ultrasonic attenuation,
, and velocity, u, of milk samples submitted to high temperature allows the detection and analysis of molecular and structural transformations in the heat coagulation process.
In various studies, ultrasonic spectroscopy has been used to monitor and analyze different types of milk coagulation. Examples include the study of acidic coagulation as used in the process of making yogurt (Kudryashov et al., 2000; Buckin and Kudryashov, 2001) or rennet coagulation used for cheese making (Gunasekaran and Ay, 1996; Povey, 1997; Nassar et al., 2001). A study of heat-induced coagulation of calcium-fortified milk was previously carried out, leading to the development of calcium-fortified milk that is stable at temperatures <100°C (Smyth et al., 1999). However, no study of ultrasonic parameters of milk has ever been achieved at temperatures >100°C. This is due to the fact that until recently ultrasonic devices did not provide the required temperature control at such high temperatures. Another problem was the requirement of sample pressurization to avoid boiling at temperatures >100°C. Introduction by Ultrasonic Scientific Ltd. (Dublin, Ireland) of the HRUS 102 spectrometer provides a new flexible tool perfectly adapted for these conditions. The measuring cells can easily be put under pressure and withstand a temperature of 120°C. Under these conditions, the HRUS 102 spectrometer provides high precision measurements of both ultrasonic attenuation and velocity, thus allowing monitoring of phase transitions occurring in the sample. Processes such as enzymatic reactions, which are expected to produce only small changes in ultrasonic velocity (in the order of a few cm/ s) require resolutions down to 0.5 mm/s to be properly analyzed. This order of resolution can be achieved by HRUS, making it a very powerful analytical technique. In the present work, we used HRUS as a new method for the measurement and characterization of the heat stability and the coagulation curves of milk powder solutions and milk ingredients.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The REM samples were then split into smaller samples, and their pH was adjusted to values covering the range from 6.0 to 7.0 using dilute HCl or dilute NaOH. These samples were left overnight at room temperature to reach pH equilibrium, and their pH was measured before introducing them in the measuring cell to determine their coagulation time.
All of the milk ingredient solutions and milk model solutions were prepared using simulated milk ultrafiltrate (SMUF) at concentrated milk level (SMUFx2) as a buffer. The SMUFx2 is a solution designed to maintain the salt composition at concentrated milk levels, which should aid in maintaining the integrity of the casein micellar system. The SMUFx2 was prepared following the method described by Jenness and Koops (1962), and its composition is as follows (mg/100 mL): Ca, 71.8; Na, 84.0; K, 308.0, Mg, 15.6; P, 72.0; Cl, 230.0; citrate, 370.2; sulfate, 19.8; and CO2, 19.0.
Phosphocasein
Phosphocasein powder (84% proteins, 4 to 5% moisture, 1 to 2% lactose, 1 to 2% fat, and 7 to 8% ash) was obtained by following the procedure described by Kelly et al. (2000) and provided by The Dairy Products Research Centre (Moorepark, Co. Cork, Ireland). The samples were prepared by dissolving the correct amount of phosphocasein in SMUFx2 to obtain a 5% (by weight) protein solution (similar to the amount of phosphocasein present in 20% solid reconstituted SMP); 0.04% sodium azide was added as a preservative.
The samples were then split into several smaller samples, and their pH was adjusted in the same manner as for the REM samples.
Whey Protein
Whey protein isolate (WPI) containing 89.6% protein (4 to 5% moisture and 5% ash) was obtained from Davisco International Inc. (MN). The WPI solutions were prepared in the same manner as REM and phosphocasein. The sample contained 1.5% protein by weight, which is similar to the amount of whey protein present in 20% solid reconstituted SMP.
Milk Model
Milk model solutions were prepared by dissolution of phosphocasein and WPI in SMUFx2 at the same concentrations as previously used.
Ultrasonic Measurements
Ultrasonic parameters in the samples were measured using HR-US 102 high-resolution spectrometer (Ultrasonic Scientific, Dublin, Ireland) fitted with two 1-mL ultrasonic cells. The spectrometer can work in differential and single cell modes. In differential mode, one of the cells contains the sample, and the second is filled with the reference liquid; the difference in ultrasonic velocity between the cells is measured. The temperature controller, Haake F8 heating bath provided by the manufacturer of HR-US 102 device, allowed maintenance of the required temperature within ± 0.01°C. Tetraethylene glycol was used as a heating liquid. The limiting resolution for measurements of ultrasonic velocity in differential mode was ±0.2 mm/s. As our measurements did not require such high resolution, single cell mode was used in the current project. In this mode, the reference cell is not used, and ultrasonic velocity in the sample cell only is monitored. Single cell mode provides resolution at ±1.5 cm/s for ultrasonic velocity and down to 0.2% for ultrasonic attenuation. The device was calibrated with water, for which ultrasonic parameters are well known over a broad temperature range (Del Grosso and Mader, 1972).
We monitored several frequencies in the frequency range from 2.5 to 14 MHz. The ultrasonic velocity u and attenuation
were automatically recorded during the coagulation process at 120°C.
Ultrasonic velocity is the distance traveled by the wave through the sample per unit of time. Ultrasonic attenuation,
, is defined by the exponential decay of ultrasonic wave A = A0e
x, where A is the amplitude the of the wave traveled a distance x in the sample and A0 is the initial amplitude of the wave.
Pressure Control
Because the experiments took place at a temperature exceeding the boiling point of milk (120°C), the pressurization of the measuring cell was required. A TP1 hand-held SI Pressure Instruments pneumatic pump (Birmingham, UK) was fitted to the measuring cell and used to maintain a constant pressure of 2 bar throughout the duration of the experiment. Once the cell was pressurized, the measurements were started. The times were adjusted afterward so that on every graph, t = 0 corresponds to the time of injection of the sample into the measuring cell.
Particle Sizing
In homogeneous materials, velocity and attenuation are determined by intrinsic properties of the medium. Heterogeneous dispersions produce an additional scattering contribution to ultrasonic velocity and attenuation, which is a function of particle size. In the long wavelength limit [i.e. when the wavelength of ultrasound (
) is much greater than the particle radius (r)], explicit expressions for the ultrasonic scattering in dispersions have been derived (Epstein and Carhart, 1953; Waterman and Truell, 1961; Allegra and Hawley, 1972; Povey, 1997). The basic mechanism of interaction of ultrasonic wave with particles in dispersions in this regime consists of 2 major contributors, thermoelastic and viscoinertial scattering. Thermoelastic scattering results from oscillations of temperature caused by the adiabatic compressions and decompressions in ultrasonic waves. These oscillations in temperature create a heat flow between the particle and the surroundings. This flow heats and cools the boundary layer between the particle and the surroundings with an associated expansion and contraction of the boundary layer, thus turning the particle into a secondary source of ultrasonic waves. In addition, a thermal wave will flow away from the particle and attenuate. The source of viscoinertial scattering lies in the difference in density between the scatterer and its surroundings. The oscillating forces associated with the deformation in the ultrasonic wave result in a motion of particles relative to their surroundings. This motion increases with increasing density difference between the particles and the surrounding fluid. The relative motion is a source of new waves and a shear field propagating away from the border between the particle and surrounding fluid, which attenuates away over a distance.
In the most comprehensive theories of ultrasonic scattering, the ultrasonic properties of an ensemble of scatterers are characterized in terms of a complex wave vector K. The Psize289 software provided by Ultrasonic Scientific Ltd. is based on the most widely used multiple scattering theoretical approaches (Epstein and Carhart, 1953; Waterman and Truell, 1961; Allegra and Hawley, 1972; Povey, 1997). This software utilizes theoretical relationships between the complex wave vector and the scattering coefficients, which are determined by the physical properties of the particles and the continuous medium.
Subjective Method
Analysis of the heat stability of REM was carried out according to the Irish Dairy Board (1974) method in the laboratory of The Dairy Products Research Centre (Moorepark, Co. Cork, Ireland) for comparison purposes.
| RESULTS |
|---|
|
|
|---|
|
The ultrasonic attenuation profiles of the different samples show a variety of shapes depending on the nature of the sample and its pH, as can be seen in Figures 2
and 3
. This is expected as the scattering contribution to ultrasonic attenuation is highly dependent on particles sizes, which in turn are highly pH dependent (Devold et al., 2000). In samples containing casein micelles, a sharp initial increase in attenuation is followed by a further, slower increase. The heat coagulation (detected as a sharp decrease in ultrasonic velocity) is then marked by a sudden decrease in ultrasonic attenuation (Figure 2
), which can sometimes be preceded by a peak in attenuation (Figure 3
). In whey protein samples, attenuation increases very slowly initially. A peak in attenuation profile coincides with a sharp decrease in ultrasonic velocity (not shown on the figure) and can be attributed to a coagulation process.
|
|
The following physical parameters for dispersed and continuous phases were used in the calculations: 1) dispersed phase, density,
= 1076 kg/m3; specific heat capacity, Cp = 3818 J/kg of K; thermal conductivity,
= 0.521 J/ ms per K, ultrasonic velocity, u = 1563 m/ s; shear modulus, µ = 0.298 GPa; thermal expansion coefficient, ß = 3.1 x 104 K1, and the attenuation at 2.5 MHz,
= 19.37 m1; 2) continuous phase, viscosity,
0 = 2.83 x 104 Pa/s, density,
o = 943.4 kg/m3, specific heat capacity, Cp0 = 4246 J/kg per K, thermal conductivity,
0 = 0.675 J/ms per K, ultrasonic velocity, u0 = 1516 m/s, thermal expansion coefficient, ß0 = 9.0 x 104 K1, and the attenuation at 2.5 MHz,
0 = 0.156 m1. This set of parameters is based on data obtained by Griffin and Griffin (1990) for a temperature of 20°C. To correct the parameters of continuous medium for temperature difference, we have modified the original values given by Griffin and Griffin (1990) by the amount typically observed in liquids for an increase in temperature of 100°C (Smithsonian Institution, 1920; Schaaffs, 1967; American Institute of Physics, 1972; Yaws, 1977; Brandrup and Immergut, 1989; CRC Handbook of Chemistry and Physics, 2000). The parameters for dispersed phase were not corrected, as typically solids show small variation of the parameters with temperature (
20% between room temperature and 120°C).
To evaluate the effect of possible inaccuracy in physical parameters on our estimations, we have varied reasonable combinations of the parameters by a ± 50% change of Griffins values (for both dispersed and continuous phases) and observed that the maximum on attenuation vs. particle size curve (Figure 2
) always remained between 400 and 1000 nm and exceeds the size of casein micelles by several times (Carroll et al., 1971; McMahon and Brown, 1984). In addition, we have found that viscosity (viscous scattering contribution) does not have a significant effect on attenuation in our systems.
Overall, the uncertainty in extrapolated to 120°C physical properties of continuous medium and the particles does not allow exact calculation of particle size in our system. However, regardless of this, qualitative analysis of the aggregation processes can be performed based on the shape of ultrasonic attenuation profiles. Therefore, the position of the maximum of the attenuation vs. particle size curve above 400 nm allows us to conclude that initially, the evolution of ultrasonic attenuation follows the left (rising) side of the curve and that the initial increase in attenuation could then be interpreted as an increase in the size of casein micelles. In milk, Freeman and Mangino (1981) showed that initial heat-induced micellar growth was followed by micellar disintegration, inducing complex changes in particle sizes. In concentrated milk, however, Singh and Creamer (1991a,b) reported that heat treatments caused an increase in average size of the casein micelles and reduced the number of smaller micelles, which is coherent with our results. They also observed increased dissociation of casein from the micelles with heating time and interpreted this by the fact that in concentrated milk, micellar aggregation, and dissociation occur simultaneously.
Upon coagulation, fast formation of large aggregates exceeding 1000 nm effective average diameter induces a slide to the right (falling) part of the curve. Therefore, a further increase in size then causes a decrease in attenuation.
The frequency dependence of ultrasonic attenuation is illustrated in Figure 4
. An increase in ultrasonic attenuation is observed with increasing frequency. At all frequencies, coagulation of the sample is marked by a drop in attenuation. Quantitative analysis of the frequency dependence of ultrasonic attenuation is the object of ongoing work and will be presented in a future paper.
|
| DISCUSSION |
|---|
|
|
|---|
|
![]() |
where u0 is the average velocity before transition (initial velocity),
u is the amplitude of the velocity change caused by the transition, tt is the time of transition (inflection point), v is the transition rate, b is controlled by where the inflexion point of the transition occurs (nearer the lower or upper asymptote), the slope of the velocity profile before transition is given by (m + n) and n is the slope after transition.
The parameters obtained by applying this algorithm to the experimental velocity data allow for the determination of the HCT for the various samples studied. Therefore, it is possible to determine their HCT-pH profiles.
The analysis of the ultrasonic attenuation profiles provides another way of determining the HCT of samples. Because the shape of attenuation profiles varies with the nature and pH of the samples, no systematic analysis was carried out. However, in cases where the attenuation profile shows a similar shape as the velocity profiles, the same algorithm was used for its analysis (Figure 5a
). In other cases (Figure 5b
), no mathematical analysis was carried out. As can be seen in Figure 5(a and b)
however, the transition is visible and happens simultaneously on both velocity and attenuation profiles.
In the case of fast coagulating samples (<5 min), the transition happens during the thermal equilibration stage and is not clearly visible on the velocity profile, as it is "masked" by the sharp initial decrease in velocity. The transition can generally be seen on the attenuation profile, however. This illustrates the benefit of simultaneously monitoring both ultrasonic velocity and attenuation.
The HCT-pH profiles determined by ultrasonic measurements are shown on Figure 6
and 7
for REM and milk ingredients, respectively. The stability profiles of REM samples obtained by ultrasonic measurements were compared with the profiles obtained by the Irish Dairy Board method (Figure 6
). The same pH dependence of the HCT is observed in the profiles obtained by those 2 methods. Overall, there is a good correlation between HRUS and the Irish Dairy Board method. However, exact positions of some points are slightly different. This could be attributed to the different conditions in which milk is coagulating in the 2 methods. For example, while the milk was submitted to intense agitation and contact with air in the Irish Dairy Board method, it was kept still during ultrasonic measurements. Interactions between the milk sample and air-surface effects present in the first case may shift slightly the coagulation time compared with coagulation in the bulk of the sample as it occurs in the ultrasonic method. In the subjective method, the lumps are detected visually on the walls of the tube while, in the ultrasonic method, the overall gelation process is detected in the bulk of the sample. However, the good correlation between the HCT profiles obtained by different methods show that the HRUS technique is an adequate new tool for the measurement of heat stability in milk.
|
|
6.5. This discrepancy can be explained by the fact that Rattray and Jelen (1997) prepared their samples using whey protein concentrate while we use WPI. Their samples would therefore contain a significant amount of lactose that our samples lack. When exposed to high temperature, lactose undergoes several reactions causing the formation of organic acid and is therefore a major source of heat-induced acidity, which could decrease the coagulation times observed in their samples compared to ours.
The stability profile obtained by ultrasonic measurements for phosphocasein (Figure 7
) is in good agreement with the results obtained by Kelly et al. (2000) with the subjective Irish Dairy Board method. The pH dependence observed with both methods is similar, but the absolute values of the coagulation times are different. This can easily be explained by the fact that Kelly carried out the coagulations experiment at 140°C, but the results presented in this work were obtained at 120°C.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication July 16, 2004. Accepted for publication May 30, 2005.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |