J. Dairy Sci. 89:4475-4489
© American Dairy Science Association, 2006.
Heat Treatment of Whole Milk by the Direct Joule EffectExperimental and Numerical Approaches to Fouling Mechanisms
L. Fillaudeau*,1,
P. Winterton*,
J. C. Leuliet
,
J. P. Tissier
,
V. Maury
,
F. Semet
,
P. Debreyne
,
M. Berthou
and
F. Chopard
* Institut National des Sciences Appliquées (INSA), Laboratoire de BiotechnologieBioprocédés, Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 5504, Institut National de la Recherche Agronomique (INRA) UMR 792-135, avenue de Rangueil, F-31077 Toulouse Cedex 4, France
INRALaboratoire de Génie des Procédés et Technologie Alimentaires (LGPTA) 369, rue Jules Guesde, 59650 Villeneuve DAscq Cedex, France
EDF-R&D, Les Renardières, 77818 Moret Sur Loing Cedex, France
Alfa-Laval Vicarb, rue du Rif Tronchard, 38120 Fontanil-Cornillon, France
1 Corresponding author: Luc.Fillaudeau{at}insa-toulouse.fr
 |
ABSTRACT
|
|---|
The development of alternative technologies such as the direct Joule effect to pasteurize and sterilize food products is of great scientific and industrial interest. Our objective was 1) to gain insight into the ability to ensure ultra-high-temperature treatment of milk and 2) to investigate the links among thermal, hydraulic, and electrical phenomena in relation to fouling in a direct Joule effect heater. The ohmic heater [OH; 

(where E is the electrical field and v is the velocity); P (power) = 15 kW] was composed of 5 flat rectangular cells [e (space between the plate and electrode) = 15 mm, w (wall) = 76 mm, and L (length of the plate in plate heat exchanger or electrode) = 246 mm]3 active cells to ensure heating and 2 (at the extremities) for electrical insulation and the recovery of leakage currents. In the first step, the thermal performance of the OH was investigated vs. the flow regimen [50 < Re (Reynolds number) < 5,000], supplied power (0 < P < 15 kW), and electrical conductivity of fluids (0.1 <
20°C < 2 S/m) under clean conditions with model fluids. This protocol enabled a global thermal approach (thermal and electrical balance, modeling of the temperature profile of a fluid) and local analysis of the wall temperature of the electrode. An empirical correlation was established to estimate the temperature gradient, Tw Tb (where Tw is the wall temperature and Tb is the product temperature) under clean conditions (without fouling) and was used to define operating conditions for pure-volume and direct-resistance heating. In the second step, the ability of OH to ensure the ultra-high-temperature treatment of whole milk was investigated and compared with a plate heat exchanger. Special care was taken to investigate the heat transfer phenomena occurring over a range of temperatures from 105 to 138°C. This temperature range corresponds to the part of the process made critical by protein and mineral fouling. The objectives were 1) to demonstrate the ability of an OH to ensure heat treatment of milk, 2) to study the thermal and hydraulic performance with an increasing power and temperature difference between the inlet and outlet of the OH, 3) to define and validate a criterion to follow heat dissipation efficiency, and 4) to compare the fouling propensity with the different configurations. A heat dissipation coefficient, RhCO, was defined and validated to monitor the fouling propensity through global electrical and thermal parameters. Finally, a numerical simulation was developed to analyze heat profiles (wall, deposit, bulk). Because of an increasing Joule effect in the static deposit, the simulation showed how wall overheating would definitively cause fouling to spiral out of control.
Key Words: ohmic heating wall overheating fouling whole milk
 |
INTRODUCTION
|
|---|
Heat treatment remains the oldest and the most frequently used preservation process in the food industry (heating, pasteurization, sterilization, cooking, and cooling). To perform these operations, heat exchangers are still fundamental pieces of equipment throughout the food processing industry. In spite of great improvements in conventional technologies over the last few decades, the heat treatment of dairy products remains a complex operation. First, fouling occurs during heat treatment and reduces the performance as well as the duration of the production cycle. Second, industrial and legal requirements in terms of sterilization level involve higher temperatures to eliminate all heat-resistant germs. In this context, any new technologies used to pasteurize or sterilize food by thermal treatment [geometrical modification, Joule effect heating, ohmic heating (OH)] in a continuous process or by nonthermal treatment (crossflow filtration, high-pressure process, high-intensity light emission, radiation, cell lysis) are of great scientific and industrial interest.
Ohmic heating is based on a simple mechanismthe generation of heat directly inside the product by the Joule effect. It consists of causing an electrical current to flow directly through the medium to be treated, flowing between a pair of electrodes. This technique can be applied to liquids containing free ions, which are therefore electrical conductors. Constant movement of electrical charges creates heat in the product according to Joules law. Heat generation is a direct function of geometry, the electrical conductivity of a given product, and voltage (Aussudre et al., 1998; Ould-El-Moctar, 1992). Ohmic heating is defined as pure-volume and direct-resistance heating, as opposed to heating by convection from the hot surface of a heat exchanger. Consequently, the heat transfer coefficient between the hot wall and the fluid is assumed to be irrelevant because there is theoretically no hot wall. From a fundamental point of view, direct-resistance heating of liquids is a complex physical problem in which strong interactions between heat transfer, hydrodynamics, and electrical phenomena can be observed (De Alwis and Fryer, 1990; Ould-El-Moctar et al., 1993). The nonexistence of a hot wall should constitute a major advantage for food applications, thereby avoiding degradation of heat-sensitive compounds (causing changes in taste, undesirable reactions, or burning) and reducing the fouling of electrode surfaces (Fillaudeau et al., 2001; Ayadi, 2005). Other advantages are related to the use of electricity because energy conversion (electrical into heat energy) is close to 100% and the uniformity of heating is greatly improved. This technology involves a short thermal inertia and fast, simple, and precise regulation is possible. Heat transfer is a function of the electrical and thermal parameters.
Ohmic heating is a technology that appeared in the 19th century when electrical energy was used for heating materials with a capacity to flow (Bhat and Joshi, 1998). A limited number of scientific articles dealing with the OH of milk have been published, mainly at the beginning of the last century (Anderson and Finkelstein, 1919; Getchell, 1935; Prescott, 1927). Unfortunately, this technology did not succeed because of the lack of suitable electrode materials and accurate, automatic electrical controls. During the past 10 yr, new and improved materials and designs for OH have become available (Aussudre et al., 1998; Berthou and Aussudre, 2000). However, the main scientific works and industrial practices are dedicated to the sterilization of low-acid food products containing particles in continuous (Biss et al., 1989; Yang et al., 1997; Eliot-Godereaux et al., 2001) and discontinuous processes (Sastry and Palaniappan, 1992; Zareifard et al., 2003). More recently, some scientific works have investigated the OH of homogeneous liquids in continuous processes (Ould-El-Moctar, 1992; El-Hajal, 1997; Fillaudeau, 2004; Ayadi, 2005) but have rarely investigated the heat treatment (pasteurization, sterilization) of real products (Fillaudeau et al., 2001).
The present work aims to understand the interactions occurring among hydraulic, thermal, and electrical phenomena in relation to fouling in a continuous ohmic apparatus [3 rectangular cells, e = 15 mm, w = 76 mm, and L = 246 mm, where e is the space between the plate or electrode, w is the wall, and L is the length of the plate in a plate heat exchanger (PHE) or electrode; 

, where E is the electrical field and v is the velocity] and to gain insight into the ability to ensure the UHT treatment of whole milk. In the first step, overheating of the electrode surface was investigated with 2 model fluids under clean conditions (without fouling) to analyze and model the temperature gradient between the wall (electrode) and fluid vs. the operating conditions. In the second step, the impact of fouling on thermal and electrical parameters was studied during the sterilization of whole milk. Deposition was compared, as observed with a PHE. A heat dissipation coefficient, RhCO, based on global electrical and thermal parameters, was defined and validated to monitor fouling. Finally, a numerical simulation was developed to analyze thermal profiles (wall, deposit, fluid) along the ohmic heater. Because of an increasing Joule effect in the static deposit, wall overheating can cause fouling to spiral out of control.
 |
MATERIALS AND METHODS
|
|---|
Experimental Setup
The experimental setup consisted of 3 different sections: 1) a preheating zone (heater I), 2) a heating zone with a PHE or ohmic heater (heater II), and 3) a cooling zone. It included the following elements in series: 2 storage tanks (2.5 m3), a constant level tank, a volumetric feed pump (type K320SD; Inoxpa, Banyoles, Spain; Qmax = 800 L/h, where Q is the volume flow rate), a preheater (heater I), a heater (heater II), a tubular heat exchanger as cooler, and a manual valve at the installation outlet.
The ohmic heater was composed of 5 rectangular cells in series (e = 15 mm, w = 76 mm, and L = 246 mm): 3 cells (1, 2, and 3) that ensured heating and 2 (at the extremities) for electrical insulation and the recovery of leakage currents (Fillaudeau et al., 2004; Ayadi, 2005). Electrodes and gaskets were placed between cells to maintain waterproofing and to enable the flow of the product, and 180° junctions ensured a hydraulic connection between cells. The geometry of the ohmic cells imposed an electrical field perpendicular to the flow pattern, 

. Temperature regulation at the ohmic heater outlet was carried out by controlling the electrical power (15-kW generator, 3-phase system, triangular connection). The power was controlled by a phase angle (voltage and current regulation, AC current, 50 Hz) up to the maximum acceptable voltage and intensity (Umax = 100 V, Imax = 170 A, where U is the electrical potential and I is the current). Electrolytic reactions may occur at the electrodesolution interface if continuous or alternating currents flow through an electrolyte. Electrolytic reactions can be prevented if the potential drop at the electrodesolution interface can be kept below the critical electrode potential, which depends on the physical and chemical properties of the electrode, pH values, current frequency, and density (Amatore et al., 1998; Sastry, 2003; Samaranayake and Sastry, 2005). In our working conditions, electrolysis is unlikely because of the current frequency (50 Hz) and density (<5,000 A/m2) used with a specific electrode material named a dimension-stable anode. Dimension-stable anode electrodes have been patented (Berthou et al., 1998) and are made out of titanium with a layer of noble metal oxides (Ru, Ir, Ta, etc.).
With fouling fluids (whole milk), experiments were carried out on 2 experimental setups (Fillaudeau et al., 2001). The first experimental device, used for configuration 1, consisted of a PHE only (V7 type; Vicarb France, Fontanil-Cornillon, France) equipped with straight corrugation plates (corrugations perpendicular to the main flow direction). The second experimental pilot plant (Experiments 2 to 5) consisted of an association between a PHE (V2 type; Vicarb France) equipped with a 60° corrugation angle and an ohmic heater. The design of the pilot plant enabled heat treatment to be carried out with a high degree of autonomy (electrical power control, temperature regulation, and fluid flow). In the PHE, a countercurrent fluid flow was created, and the nitril gaskets allowed operating conditions of up to 1 MPa and 140°C. In configurations 2 to 5, heat power was entirely supplied by the ohmic heater.
Measurement of Process Parameters
Specific instrumentation enabled the measurement of flow rate, temperature, electrical conductivity, voltage, intensity, differential pressure, and relative pressure. All flow rates were measured using an electromagnetic flow meter (type DN25; Khrone, Duisberg, Germany) accurate to within 1% of the full range. Temperatures were measured by means of platinum resistance probes (#2327; type Pt 100
, precision ± 0.2°C; Sensor-Nite, Kleinostheim, Germany) placed at the inlet and outlet of each zone. Relative pressure sensors (manometer, type 4AP30, precision 0.1%; JUMO, Fulda, Germany) were placed at the inlet and outlet of the experimental setup, whereas 2 differential pressure sensors (0 to 25 kPa and 0 to 200 kPa; Schlumberger, La Defense, France) were used on heaters I and II to monitor the exacerbation of a pressure drop when fouling occurs. In addition, the electrical conductivity (conductimeter, type 9111, accurary 0.1%; Keotron, Lyngby, Denmark) was measured at the inlet of the loop and at the outlet of the ohmic heater.
Electrode wall temperatures were measured along the ohmic heater by 12 thermocouples (ThC; sensor, type J electrically insulated ThC, #J05X, accuracy ± 0.5°C, ThC-0 to ThC-12; Corame, C2M, Mont St-Aignan, France) with nonfouling model fluids (Figure 1
). Electrical power was determined from the voltage measured between the second and third electrodes (voltmeter, 0 to 250 V, Sineax U504; Chauvin Arnoux, Paris, France) and the intensity on phase 2 of the electrical transformer (ammeter, 0 to 200 A, type AC22; Camille Bauer, Wohlen, Switzerland). All signals were electrically conditioned (module SCX-1) and collected using a data acquisition card (AT-MOI-16E-10; Agilent Technologies, Palo Alto, CA). A NI-DAQ (National Instruments-Data Acquisition; National Instruments, Austin, TX) software driver enabled configuration and control of the data acquisition system. Measurements were saved on a PC (Pentium, 200 MHz) with LabVIEW software (National Instruments).

View larger version (16K):
[in this window]
[in a new window]
|
Figure 1. Thermocouple (ThC) location at the electrode surface vs. length and height (position = x, y; g = gravity) for experiments with model fluids (dimensions xmax = 738 mm and ymax = 76 mm; positions x = 186, 306, 432, 552, and 678 mm and y = 25, 37.5, and 50 mm). Tbe = product inlet temperature (°C); U = electrical potential (V); Q = volume flow rate (m3/s); g = gravity (m/s2).
|
|
Properties of Fluids and Operating Conditions
With Nonfouling Fluids (Water, Sucrose Solution).
Model fluids were, respectively, 1) water and 2) a sucrose solution (50%, wt/wt). Their physical properties (
, C,
, µ, where
is the volume mass, C is the specific heat capacity,
is the thermal conductivity, and µ is the viscosity) were measured and compared with the literature (Norrish, 1967; Weast, 19831984). Electrical conductivities (
) were adjusted by sodium chloride adjunction (3 conductivities per solution) and described vs. temperature (Table 1
). Electrical conductivity of the solutions was included within the range of values of liquid food products (Fillaudeau, 2004).
Thermal performance of the ohmic heater was investigated vs. the flow regimen (50 < Re < 5,000, where Re is the Reynolds number), the power supplied (0 < P < 15 kW), and the electrical conductivity of fluids (0.1 <
20°C < 2 S/m). This protocol enabled a global thermal approach (thermal and electrical balance, modeling of the temperature profile of the fluid) and a local analysis of the wall temperature of the electrode (Figure 1
). The temperature difference, Tw Tb (where Tw is the wall temperature and Tb is the product temperature), under fouling-free conditions was then described by an empirical correlation.
With a Fouling Fluid (Whole Milk).
The experimental fluid was whole milk (Danone, Bailleul, France). The thermal dependencies of physical properties (density, viscosity, thermal conductivity, specific heat capacity, electrical conductivity, and viscosity) were taken from the literature (Alais, 1975; Bertsch, 1981; El-Hajal, 1997) and are reported in Table 2
.
A perfect knowledge of fouling that occurs within indirect UHT sterilization is essential to improve the design of this type of apparatus, as well as to compare it with other heating technologies. In this work, 5 different experiments were carried out whose experimental configurations and conditions are described in Table 3
. Milk was sterilized according to the following steps:
- Preheating from T0 to T1.
- Homogenization at a pressure close to 20 MPa.
- Heating from T1 to T2 in the heat recovery section (heater I),
- Heating from T2 to T3 in heater II (a PHE with a hot water loop under pressure, or an ohmic heater).
- Cooling to T4 in the cooling section.
A counterpressure was maintained close to 550 kPa in heater II (hottest point of the process) to prevent any boiling. Configuration 1 was carried out on a conventional PHE (V7 type; Vicarb France) and was used as a reference to analyze the deposit distribution (dry weight) and its chemical composition (lipid, protein, and mineral contents) vs. the bulb temperature. The temperature was regulated by increasing the water loop flow rate to modify the temperature profile. In configurations 2 to 5, heater I remained a PHE but with a smaller surface area (V2 type; Vicarb France) and space between the plates (which increased its sensitivity to fouling), and heater II was replaced by an ohmic heater. Electrical (2.7 < Pelec < 9.5 kW), thermal (7.3 < T3 T2 < 44.5°C, and 0.34 < Tw Tb < 1.91°C), and hydraulic (7,350 > Re > 3,540) constraints were progressively increased in the ohmic heater to evaluate its behavior during fouling by a real product (whole milk).
In each case, heaters I and II were studied independently. The operating conditions and experimental measurements are reported in Tables 4
and 5
, respectively. For the PHE, the overall heat transfer coefficient, hg, and pressure drop, DP, were recorded. At the end of each experiment, the PHE was dismantled and the dry weight of the deposit was measured when present in sufficient quantity. For the ohmic heater, the heat dissipation coefficient, RhCO, was monitored to compare the electrical and thermal performance, as was the pressure drop. In each experiment, an initial starting period was carried out using water. With the ohmic heater, the electrical conductivity of water was previously adjusted by sodium chloride adjunction as close as possible to that of milk to avoid over- or underheating. In configurations 4 and 5, the exchange surface of heater I was increased from 1.32 to 2.4 m2 to reduce its sensitivity to fouling. In configuration 4, a long run was carried out, but milk was recycled after 7 h; the objective was to observe whether the deposit properties in the ohmic heater evolved with time under constant thermal and hydraulic conditions. In configuration 5, the objective was to generate heavy fouling of the ohmic heater; to do so the flow rate was strongly reduced to increase the differences of temperature, T3 T2 (44.5°C), wall overheating, Tw Tb (+1.9°C under clean conditions), and residence time (28.3 s) and to work with a high electrical current (133 A).
Numerical Simulation: SPC Version 1.0 Software and Simulated Conditions
The software SPC version 1.0 (INRA, Paris, France) simulates the thermal performance of continuous processes during the stabilization of complex liquid food products. The aim of the software is 1) to simulate the thermal (temperature of the fluid, deposit, and wall) and hydraulic (pressure drop) performance of heat exchangers, 2) to take into account the impact of fouling, 3) to include the qualitative evolution of the product through chemical or biochemical kinetics, and 4) to compare the efficiency and working conditions of several technologies: a) cocurrent and b) countercurrent tubular heat exchangers, c) indirect and d) direct 

, and e)
/
Joule effect heaters. The source code used was based on the theoretical analysis of the different technologies. To analyze the problems, the differential equations governing heat transfer were established. The equations governing the conservation of mass, momentum, and energy for an incompressible fluid were programmed. The differential equations were then solved numerically by one-dimensional resolution vs. the length of the apparatus. A RungeKutta method with constant step was used because it provides a good compromise between accuracy and speed for convergence criteria based on thermal and mass balance.
Simulated conditions were the pasteurization of whole milk (heating from 50 to 92°C, Q = 200 L/h, RT = 15 to 16 s, where RT is mean residence time) with a variable deposit at the electrode surface. The ohmic heater was composed of 3 rectangular cells in series (e = 15 mm, w = 76 mm, and L = 246 mm), and a discretization with 100 steps along the length was performed. The SPC version 1.0 software provides local values of thermal and hydraulic parameters as well as dimensionless numbers and physical properties of the fluid. The hypotheses and boundary conditions used to solve the problem are as follows:
- The geometry of the channel is uniform (L = 246 mm, w = 76 mm, and e = 15 mm) and the 3 active cells are located in series.
- The thickness of the deposit is assumed to be constant along the electrode surface (0 < ed < 4 mm),
- All physical properties are assumed to depend on temperature, except the thermal and electrical conductivities of the deposit (
ther = 0.6 W/m/K,
20°C = 0015 S/m).
- The fluid is considered to form a homogeneous mixture with constant and uniform velocities (no velocity profiles, piston model) along the apparatus.
- The wall of the heater is assumed to be adiabatic.
- At the channel entrance, the temperature and velocity exhibit a uniform distribution.
- The heat transfer coefficient between the deposit (subject to the Joule effect) and the fluid is given by the correlation of Churchill and Ozoë (1973), as cited in Brienza et al. (1983).
- Tw Tb is null under fouling-free conditions (ed = 0 mm).
 |
RESULTS AND DISCUSSION
|
|---|
Study of Wall Temperature Without Fouling
Temperature Profile and Energy Balance.
The thermal performance of the ohmic heater was investigated vs. the flow regimen (50 < Re < 5000), power supplied (0 < P < 15 kW), and electrical conductivity of the fluids (0.1 <
20°C < 2 S/m). This protocol enabled a global thermal approach (thermal and electrical balance, modeling of the temperature profile of the fluid) and a local analysis of the wall temperature of the electrode. For OH, the equation of conservation of energy (Eq. [1]) was solved under several assumptions to determine an analytical expression (Eq. [2]) of the temperature profile (Fillaudeau, 2004):
 | [1] |
 | [2] |
First, the thermal and electrical balance was established with 160 working conditions. The equality between thermal and electrical powers, (Pelec Pther/Pther) < ± 5%, showed that the energy conversion was close to 100%:
 | [3] |
Second, the analytical expression of the temperature profile, Tb, and measurements of the wall temperature, Tw, (ThC-0 to ThC-12) enabled the comparison and discussion of temperature profiles (Tb, Tw) vs. the axial position and the height (Figure 2
). The temperature gradient, Tw Tb
Tw Tblin constitutes an important criterion. The temperature gradient, Tw Tb, increases vs. the heat power and decreases with the Reynolds number. This point is fundamental but requires being checked and modeled if the heat treatment of the fouling fluid is to be considered. Usually, the thermal profiles, Tb(x) and Tw(x), are parallel curves, which indicate a temperature gradient, Tw Tb, that is almost constant along the heater. However, under specific flow and heat conditions, a strong divergence of temperature profiles, Tw(x) and Tb(x), was observed, characterized by different wall temperatures at the bottom and the top of a section (Ould El Moctar et al., 1993; El Hajal et al., 1998). In fact, the thermophysical properties of a liquid in flow are dependent on temperature, and the temperature field in turn depends on the residence time. Natural convection can occur and generate additional movement to the main flow, and implies a mixed convection regimen. In the conditions studied here, we considered that mixed convection occurred if the standard deviation of wall temperatures in the same section (e.g., ThC-5, ThC-6, and ThC-7) was greater than 5 times the ThC precision (Beuf et al., 2001).

View larger version (29K):
[in this window]
[in a new window]
|
Figure 2. Temperature profiles for fluid [Tb(x) Tbe, Tblin(x) Tbe, and Tw(x) Tbe] vs. length (x) and heat power [sucrose solution, 50% (wt/wt); electrical conductivity, G3; Q = 700 L/h; Tbe = 20°C]. Tb(x) = thermal profile of the product; Tbe = product inlet temperature; Tw(x) = thermal profile of the wall; ThC = thermocouple; Tb = product temperature; Tblin = linear thermal profile; Tw = wall temperature; P1, P2, and P3 = heat power at 2.1, 3.9, and 8.3 kw.
|
|
Analysis of Wall Overheating.
The temperature gradient, Tw Tb, should remain as shallow as possible for the successful treatment of fouling liquid food. The main goal is to reach a temperature difference equal to or lower than those commonly generated with a PHE (<+0.5°C, most of the time +0.2 or +0.3°C). Consequently, the quantification and modeling of the temperature gradient appears essential. In Figure 3
, the evolution of the temperature gradient, Tw Tb, is reported for all experimental conditions (flow regimen, heat power, and electrical conductivity), excluding mixed convection phenomena. An increase in the temperature gradient, Tw Tb (maximal value close to +10°C), was observed when the heat power was increased and the flow regimen decreased. For a given flow regimen, the evolution of Tw Tb exhibited a linear shape vs. the heat fouling power. The slope, a = (Tw Tb)/P (Eq. [4]), was calculated and an empirical correlation established for 50 < Re < 4,500, 0 < P < 10+4 W, and 0.1 <
20°C < 2 S/m, which enabled the temperature gradient, Tw Tb, to be estimated and pure-volume and direct-resistance heating (Tw Tb < 0.5°C) to be defined:

View larger version (64K):
[in this window]
[in a new window]
|
Figure 3. Temperature difference (Tw Tb) vs. heat power and flow regimen [water and sucrose solution, 50% (wt/wt); electrical conductivities, G1, G2, and G3; x = 432 mm; Tw = ThC-7]; Tw = wall temperature; Tb = product temperature.
|
|
 | [4] |
UHT Sterilization of Whole Milk
The aim of this work was to gain insight into the ability of OH to ensure the UHT treatment of whole milk compared with a heat exchanger system. Special care was taken to investigate the heat transfer phenomena, that occur over a range of temperatures from 105 to 138°C. This temperature range corresponds to the critical part of the process causing protein and mineral fouling. The objectives were 1) to show the suitability of OH in the heat treatment of milk, 2) to study the thermal and hydraulic performance with increasing power and temperature differences between the inlet and outlet of the ohmic heater, 3) to define and validate a criterion to follow the heat dissipation efficiency, and 4) to compare the fouling propensity with different configurations.
Heat Dissipation Coefficient, RhCO.
The definition of the heat dissipation coefficient in the ohmic heater, RhCO, was based on the following assumption: "The contact between a fluid (food product) and a surface (electrode) implies the occurrence of interactions. These interactions may induce a deposit on the surface whose growth kinetics may be driven by a temperature difference or chemical reactions due to temperature. If the deposition occurs then the thermal and electrical properties of the deposit will have consequences on heat transfer mechanism as well as on process performance" (Fillandeau et al., 2001).
Ohmic heating presents several advantages, but for efficient implementation in the food industry, different factors such as the fouling propensity and electrical conductivity of deposits have to be considered and accurately quantified. When fouling occurs, additional electrical resistance can be created at the surface of the electrodes, resulting in a new electrical system. This layer is motionless and the Joule effect occurs there, just as it does in the flowing product. In any case, the electrical resistance of the layer will determine the performance of the system. In this new formulation (Eq. [5]), the first term corresponds to heat generated in the deposit and the second to that in the flowing product:
 | [5] |
Heat generated in this layer may be transferred to the fluid through heat convection or may induce reactions inside the deposit by a temperature increase. The balance will be a complex equilibrium between the properties of the layer (electrical and thermal conductivity, temperature sensitivity) and the heat transfer coefficient, depending on the hydraulic conditions.
An increase in electrical resistance from fouling induces a potential drop. The fluid is also less well heated, implying an increase in the global voltage required to maintain performance levels. If the deposit constitutes electrical resistance, the electrical consumption should be increased to deliver the same thermal power. To follow the fouling effect in the ohmic heater, a criterion is definedthe heat dissipation coefficient, RhCO (Eq. [6])based on electrical (I, U) and thermal (Q,
, C, T2, T3) power. This ratio is considered under clean conditions (t = 0) vs. time (t). If a deviation of the heat dissipation coefficient, RhCO, is observed, this means that fouling induces additional energy dissipation:
 | [6] |
Configuration 1: Fouling of a PHE.
Plate heat exchangers are commonly used in UHT treatment of milk because of their compactness and high thermal performance. Knowing how fouling occurs in indirect UHT sterilization has been essential for improving the design of this type of apparatus and especially for increasing its working time between cleanings (Lalande et al., 1984; Tissier et al., 1984). The fouling of PHE surfaces during the heat treatment of milk has been studied in terms of kinetics of formation and chemical composition (Lyster, 1965; Burton, 1968; Tissier et al., 1984). Two kinds of deposits can be distinguished: types A and B (Lyster, 1965; Burton, 1968). A type A deposit is created at temperatures close to 100°C. It is composed of proteins (50 to 60%), minerals (30 to 35%), and fat (around 1%). It has a spongy appearance and a creamy white color. The type B deposit appears at temperatures exceeding 105°C. The deposit composition is identified as minerals (70%) and a small proportion of proteins (15 to 20%) and fat (4 to 8%). It strongly adheres to the surface and is hard, dense, and compact. Tissier et al. (1984) determined the following average compositions of protein deposits: For type A, the composition was 62% ß-LG, 9% Ig, 9%
-LA, and 9%
-CN; for type B, it was 50%
-CN and 27%
sl-CN.
Fouling of the heat exchanger was observed during the experiments through the pressure drop and overall heat transfer coefficient measurements. Pressure drops give a useful indication of the geometrical change within flow channels when fouling takes place. However, because the change in pressure drop with time can be completely different among trials performed under the same conditions (Delplace et al., 1994), the overall heat transfer coefficient seemed to be the most representative measured consequence of the fouling layer formed on heat transfer surfaces. In heaters I and II, the results showed that the thermal and hydraulic performance and the analysis of the deposit were similar to previous works (Delplace, 1995). The pressure drop and thermal resistance increase in classical shapes, including an induction, a fouling, and postfouling period (Delplace, 1995; Delplace and Leuliet, 1995). In this configuration, the duration of the experiment is limited by the exacerbation of the pressure drop.
Configurations 2 to 5: Fouling of the Ohmic Heater.
In heater I (PHE), the change in pressure drop with time (Table 4
) exhibits a classical shape for a PHE (configuration 1) and supports data from the literature (Lalande et al., 1984; Delplace, 1995; Delplace and Leuliet, 1995) showing that fouling increases with temperature and decreases with a larger heat exchange surface area and space between plates. The heat transfer resistance, Rh, evolves in the same way as pressure drops because of the effect of fouling intensity. Configurations 2 to 5 confirmed that an efficient energy recovery unit should be used in association with an ohmic heater to carry out long tests (Table 4
). Figure 4
represents the evolution of deposit dry weight along heater I vs. bulk temperature. When the temperature was greater than 100°C, the fouling propensity drastically increased. We did not observe 2 separate peaks, as mentionned by Lalande et al. (1984), but did observe a regular increase in the amount of deposit with temperature, as reported by Burton (1968).
In contrast, in heater II (the ohmic heater) the performance remained constant, as shown in Table 5
. The temperature regulation through electrical parameters (intensity, tension) was nearly perfect, even during the transition from water to milk. Slight overheating (+2°C to the regulated value) was observed for 10 min in configuration 3, which was the consequence of a conductivity change between the water and milk. At this point, the power supplied by the ohmic heater provided acceptable performance. First, the variation in the pressure drop with time in heater II was not significant (Table 5
), with the value remaining smaller than 2 kPa. The constant value of the pressure drop in the ohmic heater may be explained by the large space between the plates (e = 15 mm) in comparison with the PHE (Vicarb, type V2, with e = 2.5 mm). Energy efficiency was remarkable, remaining higher than 90%. The heat dissipation coefficient, RhCO, with a value close to 1 (precision ± 0.1), showed that there was no significant deviation of electrical consumption (except in configuration 5) in spite of the existence of a layer of deposit at the surface of the electrode. As shown in Figures 5A to 5C
, a propensity to fouling was observed in cell 3, in agreement with the operating conditions. Fouling appeared to be homogeneous along the electrode surface. Experiment 4 demonstrated that deposit properties remained constant even during long runs (7 h under fouling conditions, followed by 5 h with recycled milk) when process parameters were well controlled. In contrast, configuration 5 demonstrated a rapid and irreversible fall of the thermal and electrical parameters (Figure 6
) linked to fouling intensity and thermal degradation at the electrode surface (Fillaudeau et al., 2001). Figure 5D
(cell 3, configuration 5) shows thermal degradation of the deposit, with the appearance of small, blackened craters (burning). The following explanations are proposed concerning 1) the observation of fouling at the electrode surface and 2) the deviation of RhCO in the last configuration. The existing literature (Lalande et al., 1984; Delplace and Leuliet, 1995; Delplace et al., 1997) has demonstrated that for PHE, the fouling propensity increases with the bulk temperature and the temperature difference, Tw Tb. Wall overheating under fouling-free conditions was deduced from Equation 4. The values calculated (Table 5
) ranged from 0.34 to 1.91°C, constituting a significantly high value compared with conventional PHE. When fouling occurs, a Joule effect takes place in the static deposit. The heat generated in this layer should be transferred to the fluid by convection, but in fact, it induced an increase in the wall temperature (Fillaudeau et al., 2001; Ayadi et al., 2004). A physical threshold was reached when the heat generated inside the deposit could not be evacuated, and a drift of the electrical and thermal parameters was observed (RhCO collapse), which indicated thermal degradation of the deposit (burning due to the Joule effect).

View larger version (136K):
[in this window]
[in a new window]
|
Figure 5. Observation of fouling generated by whole milk in cell 3 of the ohmic heater after Experiments 2 (A), 3 (B), 4 (C), and 5 (D).
|
|

View larger version (31K):
[in this window]
[in a new window]
|
Figure 6. Heat dissipation coefficient vs. time in the ohmic heater (configurations 2 to 5). RhCO = heat dissipation coefficient; Q = volume flow rate (m3/s).
|
|
Configurations 1 and 4: Observation of the Deposit.
Configurations 1 and 4, which were similar long runs with 2 different technologies (PHE and OH), were compared in the temperature range of 105 to 138°C. Several parameters were held constant, such as the flow rate:surface area ratio (potential fouling area, potential amount of deposit) and the velocity of the temperature increase (+0.7°C/s) in heater I. Flow regimens (Re
6,000) were similar in heater II, but they were not really comparable because of the large geometrical differences. The reduced pressure drop, reduced heat transfer resistance, and heat dissipation coefficient in heater II for configuration 1 increased drastically, whereas in configuration 4, they remained constant. This illustrates that during UHT treatment of milk, the PHE was much more sensitive to fouling than the ohmic heater. In the ohmic heater, the heat transfer resistance and pressure drop remained constant and stable with time in spite of deposits being observed at the surface of the electrodes.
In configuration 4, the deposit in the ohmic heater was compared to classical fouling occurring in the PHE (configuration 1). Between 105 and 120°C (ohmic cells 1 and 2), the deposit formed appeared similar to a type A deposit. It was voluminous and spongy like a protein deposit. Its average thickness was 0.5 mm and its outer surface seemed to form curls, parallel to the fluid flow. The color of the deposit was creamy white. From 120°C up to 135°C (ohmic cell 3, Figure 5C
), the appearance of the deposit changed markedly. It corresponded to a type B deposit (i.e., it was hard, dense, and had cracks). It adhered strongly to the wall, as would a fine mineral deposit called "milkstone," with a color varying from gray to yellow. Consequently, major fouling occurred in the ohmic heater without significant deviation of its hydrodynamic, thermal, or electrical parameters.
Numerical Investigation of Electrode and Deposit Overheating with Fouling
Overheating of poorly conducting liquids adjacent to hot surfaces not only causes burning or degradation of the liquids, but also leads to heat exchanger fouling. The local temperature in a direct Joule effect heater is then fundamental and can be determined either by physical measurement (i.e., magnetic resonance imaging; Ye et al., 2003) or by numerical simulation using a control volume code (Ould El Moctar et al., 1993).
The temperature profiles of fluid, Tb, deposit, Td, and electrode surface, Tw, vs. heat exchanger length with a deposit thickness ranging from 0 to 4 mm was simulated with the SPC version 1.0 software. For an element dx, the joules generated inside the deposit are proportional to its thickness: Heat generated in this layer should be transferred to the fluid by convection. The thermal balance among the fluid, the deposit, and the electrode results from the physical (
, C,
, µ) and electrical (
) properties of the fluid and the deposit. In addition, the flow conditions at the depositfluid interface define the magnitude of the heat convection coefficient. The mean temperature gradients between the electrode surface and bulk milk, Tw Tb, and between the deposit surface and bulk milk, Td Tb, plotted for an ohmic heater vs. the deposit thickness (Figure 7
) exhibit an exponential shape. The Joule effect inside the deposit causes overheating of the wall that makes fouling spiral out of control, because the fouling mechanism undergoes a "snowball effect." This phenomenon appears to be specific to the ohmic heater (direct Joule effect and 

) and is related to a strong interaction between the Joule effect inside the deposit and the convective transfer at its surface.

View larger version (29K):
[in this window]
[in a new window]
|
Figure 7. Evolution of the mean temperature differences, Tw Tb and Td Tb, along the ohmic heater vs. the deposit thickness (fluid, whole milk; Q = 200 L/h; Tbe = 50°C; Tbs = 92°C). Software used was SPC, version 1.0 (INRA). DT = Td Tb; Td = deposit temperature; Tb = product temperature; Tw = wall temperature; DJE = direct Joule effect.
|
|
Delplace et al. (1997) reported that the denaturation kinetics of whey protein were temperature dependent and demonstrated that in a PHE, denatured protein in the thermal boundary layer is involved in deposit formation, with a contribution of local hydrodynamics. In a realistic configuration with an ohmic heater, a velocity difference between the fluid flowing in the core and the fluid flowing close to the electrode (boundary layer) does exist; thus, the temperature gradient in clean conditions (without fouling; Eq. [4]) can generate a driving force to initiate fouling. Numerical simulation consolidates the experimental results obtained during heat treatment of whole milk, showing that instability of the thermal and electrical parameters is systematically created by fouling at the electrode surface because of an increase in the deposit and wall temperatures.
 |
CONCLUSIONS
|
|---|
The development of alternative (thermal or nonthermal) technologies to pasteurize or sterilize food in a continuous process is of great scientific and industrial interest. Moreover, the fouling of apparatus during pasteurization and sterilization (e.g., of dairy products) remains a severe and recurrent problem in the food industry. Heat treatment by the direct Joule effect exhibits numerous advantages because rapid heating kinetics or homogeneous heat treatment is required. Ohmic heating needs to be investigated with model fluids (under fouling-free conditions) and with real fluids (with fouling). In this work, the thermal, electrical, and hydraulic performance of a continuous ohmic apparatus (3 rectangular cells, e = 15 mm, w = 76 mm, and L = 246 mm, 

) is reported under fouling-free conditions and with fouling caused by whole milk. Fouling at the electrode surface was observed and compared with the conventional type A and B deposits identified in a PHE. In addition, a numerical calculation of temperature profiles (bulk, deposit, and wall) along the ohmic heater was performed for different fouling thicknesses.
The technology investigated shows that several phenomena should be avoided during OH of fouling liquid food:
- Appearance of mixed convection. Under a laminar flow regimen, a velocity difference between the fluid flowing in the core and the fluid flowing close to the electrode (boundary layer) exists, and the resultant temperature gradient can generate free convection superimposed on the main stream. The results show acceleration of the fluid close to the electrodes, emphasizing the effect of residence time and the mixed convection mechanism (Ould El Moctar et al., 1993; El-Hajal et al., 1998). Consequently, mixed convection increases the heterogeneity of the heat treatment and generates wall overheating along a section. The temperature gradient, Tw Tb, constitutes a driving force of fouling and should be maintained as low as possible, as opposed to mixed convection.
- A steep temperature gradient between the fluid and electrode, Tw Tb. Control of the temperature difference is essential because overheating of poorly conducting liquids adjacent to hot surfaces not only causes burning or degradation of the liquids, but exacerbates fouling. An empirical correlation was established to estimate the temperature gradient, Tw Tb, under clean conditions. It defines the operating conditions needed to obtain pure-volume and direct-resistance heating.
- Growth of the fouling deposit. The numerical simulation of working conditions with a deposit demonstrates the difficulty of heating a fouling liquid food by the Joule effect. This analysis highlights a strong interaction between the Joule effect inside the deposit and convective transfer at its surface. Any wall overheating would make the fouling snowball.
- Deviation of the heat dissipation coefficient, RhCO. During the sterilization of whole milk, we observed that the deposit formed inside the ohmic heater appeared similar to classical type A and B fouling in a PHE, although chemical analysis could not be performed. Performance of the ohmic heater remained constant even with a deposit, whereas in a PHE, the heat transfer coefficient decreased and the pressure drop became greater. However, a physical threshold is reached when the heat generated inside the deposit can no longer be evacuated. A drift of the electrical and thermal parameters is then observed (RhCO collapse), which indicates thermal degradation of the deposit (burning due to the Joule effect).
In conclusion, the ohmic heater investigated exhibits a noticeable high conversion ratio from electrical to thermal power. However, its application with fouling liquid food should be restricted 1) to a turbulent flow regimen to avoid mixed convection phenomena and to reduce wall overheating, 2) to heat treatment of fluids with a small propensity to fouling, and 3) to situations in which the deposits formed are thin and have a low electrical conductivity.
 |
APPENDIX: NOMENCLATURE
|
|---|
| Latin Letters |
|
| C |
Specific heat capacity (J/kg per K) |
| DP |
Pressure drop (Pa) |
| e |
Space between plate or electrode (m) |
| E |
Electrical field (V/m) |
| EE |
Energy efficiency (%) |
| hg |
Overall heat transfer coefficient (W/K per m2) |
| g |
Gravity (m/s2) |
| I |
Current (A) |
| L |
Length of plate (PHE) or electrode (m) |
| LM |
Lipid content (% wt/wt) |
| m |
Temperature factor (/°C) |
| P |
Power (W) |
| Q |
Volume flow rate (m3/s) |
| r |
Electrical resistance ( ) |
| Re |
Reynolds number (no units) |
| Rh |
Reduced heat transfer coefficient (no units) |
| RhCO |
Heat dissipation coefficient (no units) |
| RT |
Mean residence time (s) |
| S |
Surface area (m2) |
| SP |
Section (m2) |
| t |
Time (s) |
| T |
Temperature (°C) |
| U |
Electrical potential (V) |
| v |
Velocity (m/s) |
| w |
Width of plate or electrode (m) |
| x, y |
Lengthwise and widthwise position (m) |
|
| Greek Letters |
|
 |
Slope (°C/W) |
| µ |
Viscosity (Pa·s) |
 |
Thermal conductivity (W/m per K) |
 |
Volume mass (kg/m3) |
 |
Electrical conductivity (S/m) |
V |
Viscous dissipation (W/m3) |
|
| Subscripts |
|
| 0 |
Initial condition |
| b |
Product (bulk) |
| d |
Deposit |
| e |
Inlet |
| elec |
Electrical |
| lin |
Linear |
| s |
Outlet |
| ther |
Thermal |
| w |
Wall |
|
Received for publication February 23, 2006.
Accepted for publication June 14, 2006.
 |
REFERENCES
|
|---|
Alais, C. 1975. Sciences du lait: Principes des techniques laitières. 3rd ed. Ed. Lavoisier Technique & Documentation, Paris, France.
Amatore, C., M. Berthou, and S. Herbert. 1998. Fundamental Principles of Electrochemical Ohmic Heating of Solutions. J. Electroanalyt. Chem. 457:191203.
Anderson, A. K., and R. Finkelstein. 1919. A study of the electro-pure process of treating milk. J. Dairy Sci. 2:374406.[Abstract/Free Full Text]
Aussudre, C., M. Berthou, and P. Terrien. 1998. Présentation des techniques de chauffage des fluides par effet Joule dans lindustrie. Ed. DER-EDF, Moret Sur Loing, France.
Ayadi, M. A. 2005. Traitement thermique des fluides alimentaires encrassants par la technologie de chauffage ohmique en géométrie rectangulaire. Ph.D. Thesis, Université Henri Poincarré, Nancy, France.
Ayadi, M. A., J. C. Leuliet, F. Chopard, M. Berthou, and M. Lebouché. 2004. Continuous ohmic heating unit under whey protein fouling. Innov. Food Sci. Emerg. Technol. 5:465473.
Berthou, M., and C. Aussudre. 2000. Panorama sur le chauffage ohmique dans lindustrie agro-alimentaire. Ind. Alim. Agric. 7/8:3138.
Berthou, M., M. H. Laurent, and S. Herbert. 1998. Procédé de détermination du comportement electrochimique dun matériau délectrode et de mesure de chute de potentiel interfaciale électrodes/solution dans les conditions de chauffage ohmique, dispositif et applications à cet effet. French Pat. No. 2,776,073.
Bertsch, A. 1981. Propriétés physiques du lait et des crèmes au dessus de 100°C. Ph.D. Thesis, Ecole Nationale Supérieure des Industries Agricoles et Alimentaires, Massy, France.
Beuf, M., A. Legrand, L. Fillaudeau, J. C. Leuliet, M. Berthou, and P. Terrien. 2001. Critère dapparition de la convection mixte en régime laminaire dans un échangeur tubulaire à flux de chaleur constant (TPC). Vol. 15. Pages 6774 in Récent Progrès en Génie des Procédés. A. Strock, J. Boudrant, and D. Tondeur, ed. Ed. Lavoisier Technique & Documentation, Paris, France.
Bhat, A., and V. K. Joshi. 1998. Ohmic processing of foods: The concept, application, present status and future outlook. Alimentaria (EnneroFebrero):8388.
Biss, C. H., S. A. Coombes, and P. J. Skudder. 1989. The development and application of ohmic heating for the continuous heating of particulate food stuffs. Pages 1725 in Proc. Process Eng. Food Ind. Elsevier, Essex, England.
Brienza, B. M., J. B. Gandy, and L. Lackenback, ed. 1983. Fluid mechanics and heat transfer. In Heat Exchange Design Handbook. Vol. 2. Hemisphere Publ. Corp., Washington, DC.
Burton, H. 1968. Reviews of the progress of dairy science. J. Dairy Res. 35:317330.
De Alwis, A. A. P., and P. J. Fryer. 1990. The use of direct resistance heating in the food industry. J. Food Eng. 11:327.
Delplace, F. 1995. Identification des échangeurs de chaleur à plaques. Application à létude de lencrassement par les produits laitiers. Ph.D. Thesis, Université Henri Poincarré, Nancy, France.
Delplace, F., and J. C. Leuliet. 1995. Modelling fouling of a plate heat exchanger with different flow arrangements by whey protein solutions. Trans. IChemE. 73:112120.
Delplace, F., J. C. Leuliet, and D. Levieux. 1997. A reaction engineering approach to the analysis of fouling by whey proteins of a six-channels-per-pass plate heat exchanger. J. Food Eng. 34:91108.
Delplace, F., J. C. Leuliet, and J. P. Tissier. 1994. Fouling experiments of a PHE by whey protein solutions. Trans. IChemE. 72:163169.
El-Hajal, J. 1997. Etude expérimentale et numérique de la convection mixte dans un écoulement de Poiseuille en présence dune dissipation volumique dénergie par conduction électrique directe. Ph.D. Thesis, Institut des Sciences de lIngénieur en Thermique, Energétique et Matérieu, Nantes, France.
El-Hajal, J., A. Ould El Moctar, and H. Peerhossaini. 1998. Convection mixte dans un écoulement de Poiseuille vertical: Étude comparative de chauffage volumique et pariétal à flux constant. Int. Comm. Heat Mass Transfer 25:309319.
Eliot-Godereaux, S., F. Zuber, and A. Goullieux. 2001. Processing and stabilisation of cauliflower by ohmic heating technology. Innov. Food Sci. Emerg. Technol. 2:279287.
Fillaudeau, L. 2004. Chauffage des fluides agroalimentaires par effet Joule direct: Conductivité électrique et expressions analytiques des profils de température. Ind. Alim. Agric. 26:918.
Fillaudeau, L., G. Delaplace, J. C. Leuliet, J. P. Tissier, M. Berthou, and F. Chopard. 2001. Ohmic heating to achieve UHT sterilization of milkComparison with plate heat exchanger and definition of a heat dissipation coefficient. Pages 713718 in ExHFT-5: Experimental Heat Transfer, Fluid Mechanics and Thermodynamics 2001. Vol. 1. G. P. Celata, P. Di Marco, A. Goulas, and A. Mariani, ed. Edizioni ETS, Pisa, Italy.
Fillaudeau, L., V. Maury, and P. Debreyne. 2004. Echangeur de chaleur à effet Joule direct en agroalimentaire: Mesure et modélisation des températures pariétales et du fluide. Pages 135140 in Rencontres AGORAL (Nantes, France), Nov. 30Dec. 1, 2004. Ed. Lavoisier Technique & Documentation, Paris, France.
Getchell, B. E. 1935. Electric pasteurization of milk. Agric. Eng. 16:408410.
Lalande, M., J. P. Tissier, and G. Corrieu. 1984. Fouling of a plate heat exchanger used in UHT sterilisation of milk. J. Dairy Res. 51:557568.
Lyster, R. L. J. 1965. The composition of milk deposit in UHT plant. J. Dairy Res. 32:203209.
Norrish, R.S. 1967. Scientific and Technical Surveys, Selected Tables of Physical Properties of Sugar Solutions. Leatherhead/British Food Manufacturing Industries Research Association, London, UK.
Ould-El-Moctar, A. 1992. Etude des phénomènes physiques couplés lors du chauffage volumique dun liquide ionique en écoulement par conduction électrique directe. Ph.D. Thesis, Ecole centrale de Nantes, Nantes, France.
Ould-El-Moctar, A., H. Peerhossaini, P. Le Peurian, and J. P. Bardon. 1993. Ohmic heating of complex fluids. Int. J. Heat Mass Transfer 36:31433152.
Prescott, S. C. 1927. The treatment of milk by an electrical method. Am. J. Public Health 17:221223.
Samaranayake, C. P., and S. K. Sastry. 2005. Electrode and pH effects on eletrochemical reactions during ohmic heating. J. Electrochem. Chem. 577:125135.
Sastry, S. K. 2003. Ohmic heating. Pages 707711 in Encyclopedia of Agricultural, Food, and Biological Engineering. D. R. Heldman, ed. Marcel Dekker, New York, NY.
Sastry, S. K., and S. Palaniappan. 1992. Mathematical modeling and experimental studies on ohmic heating of liquidparticle mixtures in a static heater. J. Food Process Eng. 15:241261.
Tissier, J. P., M. Lalande, and G. Corrieu. 1984. A study of milk deposit on heat exchanger surface during UHT treatment. Pages 4958 in Engineering and Food: Engineering Sciences in the Food Industry. Vol. 1. B. McKenna, ed. Elsevier Applied Science, London, UK.
Weast, R. C., ed. 19831984. Handbook of Chemistry and Physics. 64th ed. CRC Press, Inc. Boca Raton, FL.
Yang, T. C. S., J. S. Cohen, R. A. Kluter, P. Tempest, C. Manvell, S. J. Blackmore, and S. Adams. 1997. Microbiological and sensory evaluation of six ohmically heated stew type foods. J. Food Qual. 20:303313.
Ye, X., R. Ruan, P. Chen, K. Chang, K. Ning, I. A. Taub, and C. Doona. 2003. Accurate and fast temperature mapping during ohmic heating using proton resonance frequency shift MRI thermometry. J. Food Eng. 59:143150.
Zareifard, M. R., H. S. Ramaswamy, M. Trigui, and M. Marcotte. 2003. Ohmic heating behaviour and electrical conductivity of two-phase food system. Innov. Food Sci. Emerg. Technol. 4:4555.