Next Article in Journal
Treatment of Cheese Whey Wastewater Using an Expanded Granular Sludge Bed (EGSB) Bioreactor with Biomethane Production
Next Article in Special Issue
Residue Char Derived from Microwave-Assisted Pyrolysis of Sludge as Adsorbent for the Removal of Methylene Blue from Aqueous Solutions
Previous Article in Journal
Hydrogenation of Trans,Trans-Muconic Acid to Bio-Adipic Acid: Mechanism Identification and Kinetic Modelling
Previous Article in Special Issue
Discrete Output Regulator Design for the Linearized Saint–Venant–Exner Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Feedback Control for a Pasteurization Process

1
Department of Chemical Engineering, Biotechnology and Environmental Technology, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
2
Dipartimento di Ingegneria Meccanica, Chimica e dei Materiali, Università degli Studi di Cagliari, Via Marengo 3, 09123 Cagliari, Italy
*
Author to whom correspondence should be addressed.
Processes 2020, 8(8), 930; https://doi.org/10.3390/pr8080930
Submission received: 1 June 2020 / Revised: 16 July 2020 / Accepted: 22 July 2020 / Published: 2 August 2020
(This article belongs to the Special Issue Redesign Processes in the Age of the Fourth Industrial Revolution)

Abstract

:
The milk pasteurization process is nonlinear in nature, and for this reason, the application of linear control algorithms does not guarantee the obtainment of the required performance in every condition. The problem is here addressed by proposing an adaptive algorithm, which was obtained by starting from an observer-based control approach. The main result is the obtainment of a simple PI-like controller structure, where the control parameters depend on the state of the system and are adapted online. The proposed algorithm was designed and applied on a simulated process, where the temperature dependence of the milk’s physical properties was considered. The control strategy was tested by simulating different situations, particularly when time-varying disturbances entered the system. The use of the adaptive rule reduces the variance generally introduced by the PI or PID controller.

1. Introduction

Pasteurization has been used industrially to produce safe and stable dairy products since the late 1800s. The method is based on the observation that fermentation could be prevented through heating at a specific temperature for a defined amount of time. In particular, the method targets the destruction of vegetative pathogenic bacteria, such as Mycobacterium tuberculosis, Coxiella burnetii, salmonellas, Listeria monocytogenes, campylobacters, and Escherichia coli O157:H7 [1].
During the 1900s, technology progressed fast in developing different pasteurization techniques [2]. The dominant technology in the first half of the 1900s was the holding method or vat pasteurization. This technique is based on heating the milk in batch mode at about 63 °C for 30 min. However, high-temperature short-time (HTST) continuous techniques were rapidly accepted at an industrial scale already in the early 1940s [3].
Over the years, different temperature-application times have been proposed, following the technological evolution of the process. Ohmic heating, microwave heating, pulse electric field, and ultrasound are some of the indirect or nonthermal alternative methods to heating processes [4,5,6,7]. Nevertheless, not all of them were industrially accepted compared to traditional heating procedures, and in some cases, their efficiency in reducing the bacterial count was lower than standard pasteurization [8].
Nowadays, HTST is a widely accepted technology for milk pasteurization and it is usually performed in plate heat exchangers.
A modern pasteurization plant consists of different heating and cooling sections that ensure that the milk reaches the required temperatures. Generally, the milk flows from the storage tank into the regeneration section where it is heated at 57–68 °C using pasteurized milk as a heating agent, allowing a heat recovery of up to 93% [9]. The preheated milk is then transferred to the heating section where it is heated up and reaches the final pasteurization temperature. A coiled pipe, called the holding tube, is placed after the heating section and is used to ensure that the milk is held at the pasteurization temperature for the required time [10]. The pasteurized milk passes through a cooling section where an external cooling agent is used to return it to the storage temperature.
Two important research branches have been developed within HTST processes applied to milk pasteurization. The first is related to the optimization of the energy consumption. This branch is related to the availability of process simulators with models that are able to include the milk properties, nutritional indexes, and the thermal inactivation of the microorganisms [11]. This problem was assessed by Arteaga et al. [12], through the study of the optimal time-temperature profiles of the milk heat treatment. The second research branch is related to the control of the process. According to Eurostat (2020), 172.2 million tons of raw milk were produced in 2018. This massive production must result in products that are safe for consumption [13].
The use of a proper control scheme can be a valid tool to improve process performances, reduce costs, and reject disturbances that may significantly move the system from the desired conditions. The control of a HTST pasteurization system includes typical behaviors of many industrial processes, such as complex dynamical models with nonlinearities, making the development of the controller a difficult task. The regulation of both utilities and milk temperatures using a model predictive approach was reported in Niamsuwan et al. [14], where good results were obtained with respect to other control techniques, such as the cascade generic model control (GMC) strategy. Other studies propose different model predictive control (MPC) approaches, with the objective of minimizing the energy consumption [15,16]. Kayalvizh et al. [17] developed individual cascade controllers to maintain the temperature of the milk in the heating and cooling stages of the milk pasteurization process at optimum levels. Pour et al. [18] presented three strategies to design MPC controllers using linear parameter varying models. More recently, the same authors proposed an improved health-aware economic MPC strategy aimed at minimizing the damage of components in a pasteurization plant [19].
In the present paper, the model proposed by Niamsuwan et al. [14] was extended in order to take into account temperature dependence of the density and heat capacity of the milk, using the relationships proposed by Munir et al. [20]. The main variables of the process, which are temperature at the exit of the holding tube and the milk temperature at the end of the pasteurization plant, are controlled by means of two adaptive PI (proportional-integral) controllers, designed using an observer-based structure [21]. The main characteristic of this approach is the systematic tuning design and antiwindup. The performance of the proposed control design was assessed by considering set-point changes and disturbance rejections. The process was designed taking into account the EU legislation that requires a heat treatment of at least 72 °C for 15 s [22].

2. Process Model

On the basis of the milk pasteurization system described in Niamsuwan et al. [14,16], a HTST pasteurization system consisting of a sequence of four heat exchangers was considered. According to the scheme presented in Figure 1, in the heat exchanger sequence it is possible to identify the general sections within a pasteurization unit: regeneration, heating, and cooling. To assure proper heating and cooling of the milk, hot and cold utilities were provided in the form of a boiler and a cooling tower, respectively. Moreover, a ripper plate cooling system was used for the final cooling step. Water was used as heating or cooling media.
The flow arrangement of the heat exchangers is countercurrent for heat exchanger 1, cocurrent for heat exchanger 2, and countercurrent for the two cooling steps in heat exchangers 3 and 4. In heat exchanger 2, a cocurrent flow was chosen to reduce the thermal stress of the milk.
On the basis of the system depicted in Figure 1, first order differential equations describing the temperature change within each unit of the pasteurization system were defined. These are based on the general mass and energy balances of each unit.
The models are based on the following assumptions:
  • Perfectly mixed conditions;
  • Physical and chemical properties of liquids, such as density and heat capacity, are constant, apart from the properties of milk, which are modeled as function of temperature and milk composition;
  • Negligible variation of volume on both sides of the plate heat exchanger, the heating coil of the water tank, the boiler, the ripple plate, and the cooling tower;
  • Well-insulated piping for hot and iced water of the boiler and ripple plate, meaning omitted temperature drop at the piping surface throughout the length of the piping.
Heat exchanger:
dT m , o dt = F m ( T m , i T m , o ) V pp , i ± U pp , i A pp , i Δ T lm ρ m C p , m V pp , i
Boiler:
dT h , i dt = F h ( T h , i T h , o ) V b + m · f H f U b A b ( T h , o T a ) ρ h C p , h V b
dT h , o dt = F h ( T h , i T h , o ) V pp , 2 U pp , 2 A pp , 2 Δ T lm ρ h C p , h V pp , 2
Cooling tower:
dT t , i dt = F t ( T t , i T t , o ) V ct h v E C p , t V ct + ρ w C p , w F w T w ρ t C p , t V ct L d F c T t , o + h A ( T t , o T a ) ρ t C p , t V ct
dT t , o dt = F t ( T t , i T t , o ) V pp , 3 U pp , 3 A pp , 3 Δ T lm ρ t C p , t V pp , 3
Ripple plate:
dT i , i dt = F i ( T i , o T i , i ) V rp + π n s d f h fg A f 60 ρ i C p , i V rp v U rp A rp ( T i , i T a ) ρ i C p , i V rp
dT i , o dt = F i ( T i , i T i , o ) V pp , 4 U pp , 4 A pp , 4 Δ T lm ρ i C p , i V pp , 4
Holding tube:
T m , o = ρ m C p , m F m T m , i + T a U p 2 π d p L p ρ m C p , m F m U p 2 π d p L p
The temperature dependence of the density and heat capacity of the milk were considered. The models used to determine the properties were based on the correlations reviewed in Munir (2015).
ρ m ( T , X f a t ) = ( 1040.7 2.665 10 1   T 2.3 10 3   T 2 ) X f a t ( 1.011 + 9.76 10 3   T 4.81 10 5   T 2 )   [ k g m 3 ]
C p , m ( T , X f a t ) = 4.0283 + [ 0.0304 exp ( 0.0069 ( T 19.2424 ) 2 ) + 0.0071 exp ( 0.0686 ( T 34.1363 ) 2 ) ] X f a t 10 2 0.0182 X f a t 10 2   [ J g   K ]
where T is in °C and X fat is the milk fat fraction, which was set equal to 0.015.
The properties were evaluated for each heat exchanger by taking the average of the property evaluated at the inlet and outlet port temperature. This estimates the property change throughout the heat exchanger. Parameter values and nominal condition used in the simulation of the process are reported in Table 1.

3. Impact of Process Disturbances and Control Strategy

The effect of disturbances on the pasteurization process is important to evaluate, in order to find the best control strategy that guarantees adherence to the desired conditions. Following the analysis reported in Morison [23], and using the model proposed by Niamsuwan et al. [14], two possible disturbances were introduced in the system using simple step-variations. The analysis considers the effect of the variations of inlet milk temperature and ambient temperature variations on the milk temperature at the exit of the holding tube, and on the cold milk temperature at the end of the fourth heat exchanger. Conditions and results are reported in Table 2 in deviations variables, where the initial values are calculated at the nominal conditions.
In order to avoid unnecessary heat exposure of proteins and enzymes, which can lead to a product of poorer quality, variations of no more than about +0.5 °C are allowed for the milk in the holding tube [14]. On the contrary, if the temperature does not reach 72 °C, the milk will be rejected and must be treated once more. The results reported in Table 2 show that such conditions are always exceeded when the inlet milk temperature or the ambient temperature are increased. Furthermore, the effects on cold milk leaving the pasteurizer are not negligible and a controller to reject them is essential for the product quality and the process economy.
Considering the effect of disturbances, two feedback controllers were designed. One was used to control the temperature at the exit of the holding tube, as it is related to the safety of the milk and the product quality. The manipulated input is the fuel flow rate at the boiler. The second controlled variable is the temperature of cold milk that exits the pasteurizer, using the rotating speed at the ripple plate. A low outlet temperature is required to diminish the growth rate of microorganisms still present after the pasteurization, avoiding recontamination of the product downstream in the pasteurizer. The controller is also useful to prevent the outlet temperature becoming lower than the required value, thus reducing energy waste.

4. Control Design

As evidenced by Niamsuwan et al. [14], the sluggish dynamics and the nonlinear behavior make the pasteurization process difficult to control by means of a simple feedback action. Nevertheless, operators in the industry are most familiar with traditional PID controllers and such a scheme is usually preferred. Model-based feedback control can be a valid alternative to the use of the linear model-free PI algorithm, and the use of adaptive parameters along the process motions could represent a possible solution for controlling nonlinear time-variant processes [24]. The algorithm proposed in the present paper uses the observer-based controller method [21], but it differs from the original approach, because the parameters of the control action are adjusted as a function of the process conditions.

4.1. Observer-Based Control

The adaptive algorithm proposed here is derived from the control law proposed by Castellanos-Sahagun et al. [21]. The original algorithm is based on the nonlinear constructive theory and passivity concepts and it has the following characteristics: (i) systematic tuning design and (ii) antiwindup capability, which is obtained by the observer-based structure with linear decentralized PI components.
The linearized form of the equation describing the generic controlled output (y) and the manipulated input (u) is considered as reported in Equations (11) and (12):
d y d t = a u + b
a = f d u | o
where the coefficient a is the derivative of the vector field f describing the dynamic of y, and in the original form of this approach, it is calculated at a reference condition.
The control problem can be solved for u if a reference trajectory is set for the output y, as reported in Equation (12):
dy dt = K c y
where Kc is a constant gain.
Using Equation (11) in Equation (13), u can be obtained as reported in Equation (14):
u = ( K c y + b ) a
The variable b in Equation (14) is unknown and it represents the modeling errors, or more specifically, the nonlinear part of the model and the effect of disturbances [25]. The variable b can be estimated using the measured input u and the measured output y as reported in Equation (15):
d b ^ dt = K o ( b b ^ ) = K o ( dy dt a u b ^ )
where K o is the observer gain. It is better to avoid the use of derivative in the algorithm, therefore a coordinate change is introduced in Equations (16) and (17):
χ = b ^ K o y
d χ dt = K o χ K o ( a u + K o y )
The controller can be easily tuned using the rules proposed by Castellanos-Sahagun et al. (2005), reported in Equations (18) and (19):
K c = 1 ÷ 2 ω
K o = 5 ÷ 10 K c
The controller gain in Equation (18) represents the characteristic frequency of the controller, which should be equal or twice as large as the characteristic frequency ω of the system. In order to obtain high efficiency, the estimator should converge more quickly than the controller.

4.2. Adaptive Scheme

In this work, the algorithm proposed by Castellanos-Sahagun et al. [21] is modified by updating the value of the parameter a in Equation (12) at each sampling time. Indicating with a 1 the parameter used in the adaptive PI controller of the holding tube temperature, the approximated equation reported in Equation (20) is used to adjust the gain of the controller at the varying conditions.
a 1 = f m o , 2 m f f m o , 2 T h , o T h , o m f
where f m o , 2 indicates the function describing the dynamics of the milk at the exit of the second heat exchanger, T h , o is the hot water temperature at the exit of the boiler, and m f is the fuel rate.
The same dependence can be obtained for the cold milk temperature at the outlet of the fourth heat exchanger. Such output is controlled by manipulating the rotating speed at the ripple plate and the adjustment of the controller gain is calculated with Equation (21):
a 2 = f m o , 4 n s f m o , 4 T i , o T i , o n s
where f m o , 4 is the function describing the dynamics of the cold milk, T i , o is the iced water at the outlet of the ripple plate, and n s indicates the rotating speed.
The two coefficients a 1 and a 2 are updated at each sampling time of the controller, in order to take into account the changes in the process characteristics due to the nonlinear behavior.
The controller scheme is reported in Figure 2, where y = ( T m o 2 d ,   T m o 4 ) represents the controlled output vector, ( T m o 2 d ) is the milk temperature at the exit of the holding tube,   T m o 4 is the outlet milk temperature, y r = ( T m o 2 d r ,   T m o 4 r ) is the vector of the reference trajectories, u = ( m f , n s ) is the manipulated input vector, x = ( T m o 1 ,   T m o 2 , T h i , T h 0 , T i i , T i o , T m o 3 ) is the auxiliary measurements vector, and e is the mismatch between set-points and controlled outputs. The block “Observer” indicates the integration of Equation (15).

5. Results

The effect of disturbances in the plant can impact the properties of the product (pasteurization temperature) and the consumption of energy, both for heating and cooling the milk. For this reason, the presence of controllers in the plant can improve the efficiency of the process and guarantee the desired pasteurization conditions. As reported in the previous section, two control loops were designed. A one-way interaction was present between the two control loops, because changes in the ripple plate do not affect the pasteurization temperature. Nevertheless, when the temperature of the hot milk was changed by the controller, all downstream temperatures were affected.
A decentralized MIMO (multi input-multi output) system was used to control the pasteurization process, using the observer-based approach with adaptive gain described in the previous section. Different simulations were performed in order to evaluate the controller performances.
First, set-point tracking in the absence of input disturbances is reported in Figure 3. The temperature at the exit of the holding tube was first set to 72 °C, then to 74 °C, and again to 72 °C. The temperature of the outlet milk was set first to 8 °C and then to 4 °C. The two controllers manipulated the fuel flow rate ( m f ) and the rotating speed of the riple plate ( n s ) and were able to track the system to the desired conditions. In this simulation, the ambient temperature was set equal to 30 °C, according to the nominal conditions in Niamsuwan et al. [14].
In order to evaluate the system in a more demanding situation, a colder scenario with respect to the previous situation was considered, applying the ambient temperature variation reported in the top panel of Figure 4. Set-point changes were also applied to make the test more challenging for the controllers: milk temperature at the exit of the holding tube was first set to 74 °C, then to 72 °C. Outlet milk temperature was set first to 4 °C and then to 8 °C. The performance was still satisfactory and, even if a small offset was present due to the impossibility to completely reject the disturbances, the mismatch was always within +/−0.5 °C. A more aggressive controller could eliminate the offset, but it led to large overshoots in the holding tube temperature response and this behavior can be detrimental to the quality of the product.
A further test was conducted by varying the input milk temperature, using step changes, as reported in the top panel of Figure 5. The responses for the holding tube and exit milk temperature are shown in Figure 5, showing a good performance in terms of response rate and very small offsets. The impact on the cold milk temperature is very low, and a small amount of work of the controller is able to guarantee the desired conditions, which are 72 °C for the milk temperature leaving the holding tube and 8 °C for the outlet milk temperature.

6. Conclusions

The control problem of a milk pasteurization process was addressed by means of an adaptive PI control algorithm. In particular, the proposed control strategy was applied to a simulation model proposed in the literature, where the dependence on temperature of the milk’s physical parameters was introduced, in order to obtain a more realistic representation of the process. Two feedback control loops were designed in order to reduce disturbance effects on pasteurization temperature and milk temperature leaving the pasteurizer. The main purpose was to avoid process conditions that can damage the product, reduce microbial risks, and avoid energy consumption in the cooling phase of the product. Since the process has a nonlinear behavior, the gains of the two feedback controllers were adapted online, such that the performances of the systems were also maintained when the distance from the nominal conditions was large. The proposed approach also has the advantage of being characterized by a systematic tuning procedure and antiwindup properties. The obtained control system was evaluated in the presence of large variations of disturbances and set-point changes, and it is important to note that it was possible to reject them with very small oscillations of the controlled variables.

Author Contributions

Conceptualization M.E. and S.T.; software and formal analysis, S.T. and E.D.J.R.; methodology and resources S.T.; writing—review and editing S.T. and M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A p p , A b , A f , A r p Heat exchanger area (plate pasteurizer), boiler, peripheral flow area of ripple plate compressor and ripple plate [m2]
C p , m , C p , h , C p , t , C p , i , C p , w Heat capacity at constant pressure of milk, hot water, tap water and ice water [J kg−1 °C−1]
dfImpeller diameter at ripple plate compressor [m]
dpDiameter of holding tube [m]
EEvaporation rate for cooling tower [m3 s−1]
Fm, Fh, Ft, Fi, Fw, FcVolumetric flow of milk, hot water, tap water, ice water, make-up water and circulation water [m3 s−1]
HfHeating value of fuel [kJ kg−1]
hvLatent heat of vaporization of water [kJ kg−1]
hAHeat transfer coefficient at water surface [W K−1]
hfgLatent heat of vaporization of refrigerant at ripple plate [kJ kg−1]
LdMechanical drift loss at cooling tower [%]
LpLength of holding tube [m]
m ˙ f Fueling rate at boiler [kg s−1]
nsImpeller rotational speed at ripple plate compressor [rpm]
T m , i , T m , o Temperature of milk at inlet or outlet of heat exchanger [°C]
T h , i , T h , o Temperature of hot water at inlet or outlet of heat exchanger [°C]
T t , i , T t , o Temperature of hot water at inlet or outlet of heat exchanger [°C]
T i , i , T i , o Temperature of hot water at inlet or outlet of heat exchanger [°C]
T w Temperature of make-up water [°C]
TaAmbient temperature [°C]
tTime [s]
Δ T l m Logarithmic mean temperature difference [°C]
Upp, Ub, Urp, UpOverall heat transfer coefficient of heat exchanger, boiler, ripple plate and holding tube [W m−2 K−1]
V p p ,   V b ,   V c t ,   V r p Water volume for heat exchanger, boiler, cooling tower and ripple plate [m3]
ρ m , ρ h , ρ t , ρ i , ρ w Density of milk, hot water, tap water, ice water and make-up water [kg m−3]

References

  1. Early, R. Dairy products and milk-based food ingredients. In Natural Food Additives, Ingredients and Flavourings; Woodhead Publishing Series in Food Science; Technology and Nutrition: Sawston, UK, 2012; pp. 417–445. [Google Scholar]
  2. Rankin, S.A.; Bradley, R.L.; Miller, G.; Mildenhall, K.B. A 100-year review: A century of dairy processing advancements-pasteurization, cleaning and sanitation, and sanitary equipment design. J. Dairy Sci. 2017, 100, 9903–9915. [Google Scholar] [CrossRef]
  3. Westhoff, D.C. Heating milk for microbial destruction: A historical outline and update. J. Food Prot. 1978, 41, 122–130. [Google Scholar] [CrossRef] [PubMed]
  4. Stancl, J.; Zitny, R. Milk fouling at direct ohmic heating. J. Food Eng. 2010, 99, 437–444. [Google Scholar] [CrossRef]
  5. Valero, E.; Villamiel, M.; Sanz, J.; Martinez-Castro, I. Chemical and sensorial changes in milk pasteurized by microwave and conventional systems during cold storage. Food Chem. 2000, 70, 77–81. [Google Scholar] [CrossRef]
  6. McAuley, C.M.; Singh, T.K.; Haro-Maza, J.F.; Williams, R.; Buckow, R. Microbiological and physicochemical stability of raw, pasteurized or pulsed electric field-treated milk. Innov. Food Sci. Emerg. Technol. 2016, 38, 365–373. [Google Scholar] [CrossRef]
  7. Monteiro, S.H.M.C.; Silva, E.K.; Alvarenga, V.O.; Morares, J.; Freitas, M.Q.; Silva, M.C.; Raices, R.S.L.; Sant’Anna, A.S.; Meireles, M.A.A.; Cruz, A.G. Effects of ultrasounds energy density on the non-thermal pasteurization of chocolate milk beverage. Ultrason. Sonochem. 2018, 42, 1–10. [Google Scholar] [CrossRef] [PubMed]
  8. Lim, S.-Y.; Benner, L.C.; Clark, S. Neither thermosonication nor cold sonication is better than pasteurization for milk shelf life. J. Diary Sci. 2019, 102, 3965–3977. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, L.; Sunden, B.; Manglik, R.M. Plate Heat Exchangers: Design, Applications and Performance, 1st ed.; WIT Press: Southampton, UK, 2007. [Google Scholar]
  10. Kanade, P.; Subramani, A. Hygienic design aspects of pasteurizer to assure effective pasteurization of milk. J. Hyg. Eng. Des. 2014, 8, 19–29. [Google Scholar]
  11. Bon, J.; Clemente, G.; Vaquiro, H.; Mulet, A. Simulation and optimization of milk pasteurization processes using a general process simulator (ProSimPlus). Comput. Chem. Eng. 2010, 34, 414–420. [Google Scholar] [CrossRef]
  12. Arteaga, G.E.; Vazquez-Arteaga, M.C.; Nakai, S. Dynamic optimization of the heat treatment of milk. Food Res. Int. 1994, 27, 77–82. [Google Scholar] [CrossRef]
  13. EUROSTAT. 2020. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php/Milk_and_milk_product_statistics (accessed on 15 June 2020).
  14. Niamsuwan, S.; Kittisupakorn, P.; Mujtaba, I.M. Control of milk pasteurization process using model predictive approach. Comput. Chem. Eng. 2014, 66, 2–11. [Google Scholar] [CrossRef] [Green Version]
  15. Khadir, M.T.; Ringwood, J. Linear and nonlinear model predictive control design for a milk pasteurization plant. Control Intell. Syst. 2003, 31, 1–8. [Google Scholar]
  16. Niamsuwan, S.; Kittisupakorn, P.; Mujtaba, I.M. Optimization approach to minimize energy consumption in pasteurized milk process. In International Proceedings of Chemical, Biological & Environmental Engineering, Bangkok, Thailand; IACSIT Press: Singapore, 2011; pp. 35–39. [Google Scholar]
  17. Kayalvizhi, M.; Manamalli, D.; Bhuvanithaa, K. Cascade control of HTST milk pasteurization process with and without pre-cooling stage. In Proceedings of the International Conference on Energy, Communication, Data Analytics and Soft Computing (ICECDS), Chennai, India, 1–2 August 2017; pp. 1807–1812. [Google Scholar] [CrossRef]
  18. Pour, F.K.; Puig, V.; Ocampo-Martinez, C. Comparative assessment of LPV-based predictive control strategies for a pasteurization plant. In Proceedings of the 4th International Conference on Control, Decision and Information Technologies (CoDIT), Barcelona, Spain, 5–7 April 2017; pp. 0821–0826. [Google Scholar] [CrossRef]
  19. Pour, F.K.; Puig, V.; Ocampo-Martinez, C. Multi-layer health-aware economic predictive control of a pasteurization pilot plant. Int. J. Appl. Math. Comput. Sci. 2018, 28, 97–110. [Google Scholar] [CrossRef] [Green Version]
  20. Munir, M.T.; Zhang, Y.; Yu, W.; Wilson, D.I.; Young, B.R. Virtual milk for modelling and simulation of dairy processes. J. Dairy Sci. 2015, 99, 3380–3395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Castellanos-Sahagún, E.; Alvarez-Ramírez, J.; Alvarez, J. Two-point temperature control structure and algorithm design for binary distillation columns. Ind. Eng. Chem. Res. 2005, 44, 142–152. [Google Scholar] [CrossRef]
  22. Siguemoto, E.S.; Purgatto, E.; Hassimotto, N.M.A.; Gut, J.A.W. Comparative evaluation of flavor and nutritional quality after conventional and microwave assisted pasteurization of cloudy apple juice. LWT Food Sci. Technol. 2019, 11, 853–860. [Google Scholar] [CrossRef]
  23. Morison, K.R. Steady-state control of plate pasteurisers. Food Control 2005, 16, 23–30. [Google Scholar] [CrossRef]
  24. Tronci, S.; Baratti, R. A Gain-Scheduling PI Control Based on Neural Networks. Complexity 2017. [Google Scholar] [CrossRef] [Green Version]
  25. Tronci, S.; Grosso, M.; Alvarez, J.; Baratti, R. On the global nonlinear stochastic dynamical behavior of a class of exothermic CSTRs. J. Process Control 2011, 21, 1250–1264. [Google Scholar] [CrossRef]
Figure 1. Schematic flow diagram for the pasteurization process considered.
Figure 1. Schematic flow diagram for the pasteurization process considered.
Processes 08 00930 g001
Figure 2. Block diagram of the adaptive control scheme.
Figure 2. Block diagram of the adaptive control scheme.
Processes 08 00930 g002
Figure 3. Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (top left panel), and rotating speed at the refrigeration rate, ns (top right panel).
Figure 3. Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (top left panel), and rotating speed at the refrigeration rate, ns (top right panel).
Processes 08 00930 g003
Figure 4. Process disturbance is the ambient temperature, Ta (top panel). Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (middle left panel), and rotating speed at the refrigeration rate, ns (middle right panel).
Figure 4. Process disturbance is the ambient temperature, Ta (top panel). Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (middle left panel), and rotating speed at the refrigeration rate, ns (middle right panel).
Processes 08 00930 g004
Figure 5. Process disturbance is the inlet milk temperature, Tm1i (top panel). Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (middle left panel), and rotating speed at the refrigeration rate, ns (middle right panel).
Figure 5. Process disturbance is the inlet milk temperature, Tm1i (top panel). Controlled outputs are the temperature at the outlet of the holding tube, Tmo2d (bottom left panel) and milk temperature at the outlet of the fourth heat exchanger, Tmo4 (bottom right panel). Set-points for both Tmo2d and Tmo4 are reported with a grey dash-dotted line; manipulated inputs are fuel flow rate at boiler, mf (middle left panel), and rotating speed at the refrigeration rate, ns (middle right panel).
Processes 08 00930 g005
Table 1. Physical properties and process data for the simulation.
Table 1. Physical properties and process data for the simulation.
Heat   Transfer   Areas   [ m 2 ] Boiler Characteristics
Regenerative stage; A 1 2Fuel rate [ kg / s ] ; m f 1.8 × 10 3
Heating stage; A 2 2.2Fuel heating value [ kJ / kg ] ; H f 49888
Precooling stage; A 3 1Heat transfer area [ m 2 ] ; A b 2.3
Cooling stage; A 4 5.0
Overall   Heat   Transfer   Coefficients   [ W / m 2 K ]
Regenerative stage; U 1 940Cooling tower characteristics
Heating stage; U 2 940Heat transfer coefficient [ W / K ] ; h A 2000
Precooling stage; U 3 950Heat of vaporization [ kJ / kg ] ; h v 2410
Cooling stage; U 4 1000Circulation water [ m 3 / s ] ; F c 4.8 × 10 3
Boiler; U b 227Drift loss [%]; L d 2
Flowrates   [ m 3 / s ] Make - up   Water   Temperature   [ ° C ] ;   T w 27
Milk; F m 4 × 10 4 Evaporation rate [ m 3 / s ] ; E 1.7 × 10 5
Hot water; F h 1.6 × 10 3
Tap water; F t 4.8 × 10 3 Specifications for ripple plate compressor
Ice water; F i 1.92 × 10 3 Rotating speed [rpm]; n s 5000
Volumes   [ m 3 ] Diameter   of   Impeller   [ m ] ;   d f 0.4
Regenerative stage; V 1 0.3803Peripheral flow area [ m 2 ] ; A f 0.002
Heating stage; V 2 0.4079Heat transfer area [ m 2 ] ; A r p 2.5
Precooling stage; V 3 0.2873
Cooling stage; V 4 0.4565Specifications for refrigerant
Boiler; V b 1.2Heat of vaporization [ kJ / kg ] ; h f 217
Cooling tower; V c t 0.05Specific volume at the exit; v 0.5
Ripple plate; V r p 0.5Overall heat transfer [ W / m 2 K ] ; U r p 1120
Ambient air temperature [°C]; T a 30
Table 2. Effects of input changes on steady-state conditions of the pasteurization process.
Table 2. Effects of input changes on steady-state conditions of the pasteurization process.
Disturbance VariableDisturbance AmountChange in Temperature at the Exit of Holding TubeChange in Temperature at the Exit of Fourth Heat Exchanger
Milk feed temperature+4 °C+2.6 °C+0.94 °C
−4 °C−2.6 °C−0.92 °C
Environment temperature+4 °C+1.48 °C3.05 °C
−4 °C−1.47 °C−3.04 °C

Share and Cite

MDPI and ACS Style

Rasmussen, E.D.J.; Errico, M.; Tronci, S. Adaptive Feedback Control for a Pasteurization Process. Processes 2020, 8, 930. https://doi.org/10.3390/pr8080930

AMA Style

Rasmussen EDJ, Errico M, Tronci S. Adaptive Feedback Control for a Pasteurization Process. Processes. 2020; 8(8):930. https://doi.org/10.3390/pr8080930

Chicago/Turabian Style

Rasmussen, Emil D.J., Massimiliano Errico, and Stefania Tronci. 2020. "Adaptive Feedback Control for a Pasteurization Process" Processes 8, no. 8: 930. https://doi.org/10.3390/pr8080930

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop