Startseite Quadratic dynamic predictive control framework for continuous high-temperature-short-time pasteurization process
Artikel
Lizenziert
Nicht lizenziert Erfordert eine Authentifizierung

Quadratic dynamic predictive control framework for continuous high-temperature-short-time pasteurization process

  • Marimuthu Indumathy , Sobana Subramani , Rames C. Panda ORCID logo EMAIL logo und Atanu Panda
Veröffentlicht/Copyright: 28. Mai 2025
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

The present study focuses on developing advanced model-based controllers to enhance pasteurization temperature control in the dairy industry. The analysis and discussions of various multivariable control systems, including the cascade algorithm, model-predictive-control (MPC), and minimum-variance (MV) controller, are presented. A mathematical-model of the plant has been designed using mass and energy-balance equations for the integrated three-stage sections of the system, including the regeneration, heating, and cooling of counter-current gasketed PHE (plate heat exchanger), holding tube, holding tank, boiler, and chiller unit. The plant is modeled as a two-input and single-output system and used for the synthesis of controllers. This paper also uses a model reference (MR) framework for multivariable systems to construct the data-driven quadratic dynamic matrix controller (QDMC). Without building a mathematical model for the process, the data-driven controller is developed using one or more batches of plant data. The efficiency of the controller is evaluated by minimizing the objective function that contains the closed-loop error, and its ability to maintain control despite measurement noise is also examined. The closed-loop results (performances) are compared using the integral absolute error (IAE) performance criteria.


Corresponding author: Rames C. Panda, Chemical Engineering Department, RajaLakshmi Engg College, Chennai, India, E-mail:

Acknowledgements

The authors thank their institutes for providing facilities to carry out this research.

  1. Research ethics: Not applicable.

  2. Informed consent: Not applicable.

  3. Author contributions: The author ‘M.I.’ has carried out data curation, investigation, simulation, observation and writing the rough draft of the paper. The second author ‘S.S.’ has supervised and helped in collecting data. The third author “R.C.P.” has supervised the work, did administration, provided concepts and has corrected the final draft of the paper. The fourth author ‘A.P.’ has contributed through estimation and control related simulations, and drafting for the paper.

  4. Use of Large Language Models, AI and Machine Learning Tools: None declared.

  5. Conflict of interest: The authors state that there is no conflict of interests in this publication.

  6. Research funding: None declared.

  7. Data availability: Not applicable.

Appendix A (Augmented state and parameter estimation using adaptive square root unscented Kalman filter) [48, 50, 51]

Step 1:

Augmented state initialization

(A.1) x ˆ 0 = E x ˆ 0 P 0 = E x 0 x ˆ 0 x 0 x ˆ 0 T S 0 = chol P 0

Time updated equations

Step 2:

Computation of sigma points, weights

(A.2) χ i k 1 = x ˆ k 1 ; w i , m = λ n + λ ; w i , c = λ n + λ + 1 α 2 + β ; i = 0 χ i k 1 = x ˆ k 1 + n + λ P k 1 ; w i , m = w i , c = λ 2 n + λ ; i = 1 , . . . , n χ i k 1 = x ˆ k 1 n + λ P k 1 ; w i , m = w i , c = λ 2 n + λ ; i = n + 1 , . . . , 2 n

Step 3:

Priori estimation with updation of error co-variance

(A.3) χ i k | k 1 = f χ i k 1 , u k 1 x ˆ k | k 1 = i = 0 2 n w i , m χ i k | k 1 ; i = 0 , . . . , 2 n s x k | k 1 = q r w 1 : 2 n , c χ 1 : 2 n k | k 1 x ˆ k | k 1 Q k 1 s x k | k 1 = Cholupdate s x k | k 1 , w 0 , c χ 0 k | k 1 x ˆ k | k 1

Measurement updated equations:

Step 4:

Determining predicted observation and the covariance matrixes

(A.4) y i k | k 1 = h χ i k | k 1 , u k 1 y ˆ k | k 1 = i = 0 2 n w i , m y i k | k 1 ; i = 0 , . . . , 2 n s y k = qr w 1 : 2 n , c y 1 : 2 n k | k 1 y ˆ k | k 1 R k 1 s y k | k 1 = Cholupdate s y k , w 0 , c y 0 k | k 1 y ˆ k | k 1 P xy = i = 0 2 n w i , c χ i k | k 1 x ˆ k | k 1 y i k | k 1 y ˆ k | k 1 T P yy = i = 0 2 n w i , c y k y ˆ k | k 1 y k y ˆ k | k 1 T

Step 5:

Deriving filter gain, residual, multi step innovation vector, noise covariances

(A.5) κ k = P x y k s y k s y T k 1 e ̲ k = y k y ˆ k | k 1 ε k = 1 m i = k m + 1 k e ̲ i e ̲ T i Q k = κ k ε k κ T k R k = ε k + P y y k

Step 6:

Updation of augmented states and covariance

(A.6) x ˆ k = x ˆ k | k 1 + κ k e ̲ k P k = P k | k 1 κ k P yy k κ T k s k = cholupdate s x k | k 1 , κ k P xy , 1

Step 7:

Updation of augmented states and covariance

(A.7) k = k + 1

Appendix B:

Proof: Assuming Lyapunov candidate function as

(B.1) V k = x T k 1 P 1 k 1 x T k 1

Plunging matrix inversion lemma in 2nd equation of Eq. (A.6), it yields

(B.2) P 1 k = P 1 k | k 1 + P 1 k | k 1 P xy k | k 1 P xy k | k 1 P xy T k | k 1 P 1 k | k 1 P xy k | k 1 1 P xy T k | k 1 P 1 k | k 1 = P 1 k | k 1 + C T k R 1 k C k

Thus, 1st equation of Eq. (A.5) reduces to

(B.3) κ k = P k | k 1 C T k 1 C k 1 P k | k 1 C T k 1 + R k 1 = P k C T k R 1 k

Putting Eq. (A.6) into Eq. (67), it becomes

(B.4) x k = x k x ˆ k | k 1 + κ k y k y ˆ k = x k | k 1 κ k y k

Combining Eq. (A.5, B.1, B.4), Lyapunov function reduces to

(B.5) V k = x k | k 1 κ k e k 1 P 1 k x k | k 1 κ i e k = x T k | k 1 P 1 k x k | k 1 x T k | k 1 C T k κ T k P 1 i x k | k 1 x T k | k 1 P 1 k κ k C k x k | k 1 + x T k | k 1 C T k κ T k P 1 k κ k C i x k | k 1

Merging Eq. (B.3, B.5), it yields

(B.6) V k = x T k | k 1 P 1 k C T k R 1 k C k + C T k R 1 k C k P k C T k R 1 k C k x k | k 1 = x T k | k 1 P 1 k C T i R 1 k R 1 k C k P k C T k R 1 k C k x k | k 1

However, after simplification, the below terms reduced as

(B.7) R 1 k R 1 k C k P k C T k R 1 k = R 1 k R 1 k C k P k | k 1 C T k C i P k | k 1 C T k + R k 1 = R 1 k C i P k | k 1 C T k + R k R 1 k C k P k | k 1 C T k C k P k | k 1 C T k + R k 1 = C k P k | k 1 C T k + R k 1 = P yy 1 k

Combining Eq. (72, B.6, B.7), the Lyapunov function can be derived as

(B.8) V k = x T k | k 1 F T k η k P 1 k | k 1 C T k P yy 1 k C k η k F k x k | k 1

Since, V(k) would be a decreasing sequence if there exists a positive scalar (1>ρ>0) such that

(B.9) V k 1 ρ V k 1

Eq. (B.9) can be extended into the below LMI

(B.10) x T k | k 1 F T k η k P 1 k | k 1 C T k P yy 1 k C k η k F k 1 ρ P 1 k 1 x k | k 1 0

To obey Eq. (B.10), below necessary condition is introduced

(B.11) μ l F T k η k C T k P yy 1 C k η k F k + μ l 1 ρ P 1 k 1 μ u F T k η k P 1 k | k 1 η k F k

To corroborate Eq. (B.11), below inequalities are introduced

(B.12) μ u F 2 k μ u η k 2 μ u P 1 k | k 1 μ u F T k η k P 1 k | k 1 η k F k μ l F 2 k μ l η k 2 μ l C 2 k μ l P yy 1 μ l F T k η k C T k P yy 1 C k η k F k

Eq. (B.11) again transforms as

(B.13) μ l F 2 k μ l η k 2 μ l C 2 k μ l P y y 1 + 1 ρ μ l P 1 k 1 μ u F 2 k μ u ς 2 k μ u P 1 k | k 1

Thus, Eq. (B.13) yields as

(B.14) μ l η k 2 η k F 2 k η k C 2 k μ u F 2 k μ u P 1 k | k 1 μ u P yy 1 + 1 ρ μ l P 1 k 1 μ u F 2 k μ u P 1 k | k 1 μ u η k 2 i . e . μ l η k 2 μ l C 2 k μ l P k | k 1 μ u P yy 1 + 1 ρ μ l P 1 k 1 μ l P k | k 1 μ u F 2 k μ u η k 2

Eq. (B.14) shows that the Lyapunov function (as defined in Eq. (B.1)) is a monotonically decreasing function. To bring forward, below criteria must be met.

(B.15) lim K x T k = 0

Thus, Eq. (B.15) shows that error due to plant-model mismatch tends to zero asymptotically.

References

1. Vallery-Radot, R. Life of Pasteur. Whitefish, Montana, USA: Kessinger Publishing Co.; 1928:113–4 pp. 978-0-7661-4352-4Suche in Google Scholar

2. Goff, H, Douglas, H, Richard, W. Ice cream, 7th ed. New York, NY: Springer; 2013.10.1007/978-1-4614-6096-1Suche in Google Scholar

3. Shieh, JS, Chen, HC, Ferng, LH. Application of a fuzzy logic controller in temperature control of a pilot high-temperature short-time heat exchanger. J Food Control 1992;3:91–6. https://doi.org/10.1016/0956-7135(92)90038-c.Suche in Google Scholar

4. Negiz, A, Ramanauskas, P, Ali, C, Joseph, S, E, Armstrong, J. Automated control of high temperature short time pasteurization. J Food Control 1996;7:309–15.David10.1016/S0956-7135(96)00048-5Suche in Google Scholar

5. Schlesser, JE, Armstrong, DJ, Cinar, A, Ramanauskas, P, Negiz, A. Automated control and monitoring of thermal processing using high temperature, short time pasteurization. J Dairy Sci 1997;80:2291–6. https://doi.org/10.3168/jds.s0022-0302(97)76178-2.Suche in Google Scholar

6. Khadir, MT, Ringwood, JV. First principles modelling of a pasteurization plant for model predictive control. Math Comput Model Dyn Syst 2003;9:281–301. https://doi.org/10.1076/mcmd.9.3.281.24150.Suche in Google Scholar

7. Khadir, MT, Ringwood, VJ. Higher Order Predictive Functional Control versus dynamical matrix control for a milk pasteurization. Int J Food Eng 2013;10:17–27.10.1515/ijfe-2012-0006Suche in Google Scholar

8. Vidhyaa, CB, Vishnu Anirudh, PD, Panda, B, Sobana, S, Panda, RC, Laskar, JH. Modeling and control of vascular dementia disease by exact dosing of medicines. Int J Robust Nonlinear Control 2023:1–19. https://doi.org/10.1002/rnc.7067.Suche in Google Scholar

9. Morison, KR. Steady-state control of plate pasteurizers. J Food Control 2005;16:23–30. https://doi.org/10.1016/j.foodcont.2003.11.001.Suche in Google Scholar

10. Riverol, C, Ricart, G, Carosi, C, Di Santis, C. Application of advanced soft control strategies into the dairy industry. Innov Food Sci Emerg Technol 2008;9:298–305. https://doi.org/10.1016/j.ifset.2007.07.002.Suche in Google Scholar

11. Carola, GCCG, Gabriel, ND, Jorge, AWG. Dynamic simulation of a plate pasteurizer unit: mathematical modeling and experimental validation. J Food Eng 2014;131:124–34. https://doi.org/10.1016/j.jfoodeng.2014.01.029.Suche in Google Scholar

12. Munir, MT, Zhang, Y, Yu, W, Wilson, DI, Young, BR. Virtual milk for modeling and simulation of dairy processes. J Dairy Sci 2015;99:3380–95. https://doi.org/10.3168/jds.2015-10449.Suche in Google Scholar PubMed

13. Kayalvizhi, M, Manamalli, D, Bhuvanithaa, K. Cascade control of HTST milk pasteurization process with and without pre-cooling stage. In: International Conference on Energy, Communication, Data Analytics and Soft Computing (ICECDS-2017), Chennai, India; 2017:1807–12 pp.10.1109/ICECDS.2017.8389760Suche in Google Scholar

14. Widiatmo, JS, Hendrarsakti, J. Process control of milk pasteurization using Geo-thermal brine under Geothermal brine temperature and flow rate disturbance. IOP Conf Ser Earth Environ Sci 2019;254:1–8. https://doi.org/10.1088/1755-1315/254/1/012018.Suche in Google Scholar

15. Andrade, LJM, Rosa, LG. Temperature control of a solar furnace for material testing. Int J Syst Sci 2011;42:1253–64, https://doi.org/10.1080/00207721.2011.588894.Suche in Google Scholar

16. Rasmussen, EDJ, Errico, M, Tronci, S. Adaptive feedback control for a pasteurization process. Processes 2020;8:1–13. https://doi.org/10.3390/pr8080930.Suche in Google Scholar

17. Lalwani, S, Lewerentz, F, Håkansson, A, Babic, M, Madina, P, Edén, J, et al.. A pilot scale study on the effects of various temperature and time conditions on milk quality using direct and indirect heating systems. Int Dairy J 2024;155:105951. https://doi.org/10.1016/j.idairyj.2024.105951.Suche in Google Scholar

18. Astrom, KJ. Introduction to stochastic control theory. New York: Academic Press; 1970.Suche in Google Scholar

19. Borisson, U, Wittenmark, B. Moisture content control on a paper machine: an application of a self tuning regulator. Lund, Sweden: Department of Automatic Control, Lund Institute of Technology; 1973. Technical report (TFRT-7038).Suche in Google Scholar

20. Harris, T. Assessment of closed-loop performance. J Chem Eng 1989;67:856–61. https://doi.org/10.1002/cjce.5450670519.Suche in Google Scholar

21. Harris, T, Seppala, C, Desborough, L. A Review of performance monitoring and assessment techniques for univariate and multivariate control Systems. J Process Control 1999;9:1–17. https://doi.org/10.1016/s0959-1524(98)00031-6.Suche in Google Scholar

22. Desborough, L, Harris, T. Performance assessment measures for univariate feedforward/feedback Control. J Chem Eng 1993;71:605–16. https://doi.org/10.1002/cjce.5450710413.Suche in Google Scholar

23. Modares, A, Sadati, N, Modares, H. Data-driven safe gain-scheduling control. Asian J Control 2023;25:4171–82H. https://doi.org/10.1002/asjc.3169.Suche in Google Scholar

24. Başaran, A, Yılmaz, T, Çivi, C. Experimental investigation and evaluation of the thermodynamical performance of a novel hybrid design for milk pasteurization. Therm Sci Eng Prog 2023;44:102028. https://doi.org/10.1016/j.tsep.2023.102028.Suche in Google Scholar

25. Wang, J. Minimum variance control strategy for closed loop linear time invariant system. Int J Syst Syst Eng 2019;9:62–74. https://doi.org/10.1504/IJSSE.2019.097908.10.1504/IJSSE.2019.097908Suche in Google Scholar

26. Garcia-Sanz, M, Guillen, JC, Ibarrola, JJ. Robust Controller design for uncertain systems with variable time delay. J Control Eng Pract 2001;9:961–72. https://doi.org/10.1016/s0967-0661(01)00040-5.Suche in Google Scholar

27. Ibarrola, JJ, Sandoval, JM, Garcia-Sanz, M, Pinzolas, M. Predictive control of a high temperature-short time pasteurization process. Control Eng Pract 2002;10:713–25.10.1016/S0967-0661(02)00029-1Suche in Google Scholar

28. Niamusuwan, S, Kittisupakorn, P, Mujtaba, IM. Control of milk pasteurization process using model predictive approach. Comput Chem Eng 2014;66:2–11. https://doi.org/10.1016/j.compchemeng.2014.01.018.Suche in Google Scholar

29. Hjalmarsson, H, Gevers, M, Gunnarsson, S, Lequin, O. Iterative feedback tuning: theory and applications. IEEE Control Syst Mag 1998;18:26–41.10.1109/37.710876Suche in Google Scholar

30. Karimi, A, Miskovic, L, Bonvin, D. Iterative correlation-based controller tuning. Int J Adapt Control Signal Process 2004;18:645–64.10.1002/acs.825Suche in Google Scholar

31. Miskovic, L, Karimi, A, Bonvin, D, Gevers, M. Correlation-based tuning of decoupling multivariable controllers. Automatica 2007;43:1481–94. https://doi.org/10.1016/j.automatica.2007.02.006.Suche in Google Scholar

32. Daniel, DH, da Silva, GRG, Lucíola, C. Data-driven control design by prediction error identification for a refrigeration system based on vapor compression. IFAC-Papers On Line 2018;51:704–9. https://doi.org/10.1016/j.ifacol.2018.06.186.Suche in Google Scholar

33. Bazanella, AS, Campestrini, L, Eckhard, D. Data-driven controller design: the H2 approach. Netherlands: Springer; 2012.Suche in Google Scholar

34. Kadakal, C, Duman, T, Ekinci, R. Thermal degradation kinetics of ascorbic acid, thiamine and riboflavin in rosehip (Rosa canina L) nectar. Food Sci Technol 2018;38:667–73. https://doi.org/10.1590/1678-457x.11417.Suche in Google Scholar

35. Indumathy, M, Sobana, S, Ramesh, CP. Identification of high temperature short time milk pasteurization unit. J Food Process Eng 2020. https://doi.org/10.1111/jfpe.13410.Suche in Google Scholar

36. Indumathy, M, Sobana, S, Ramesh, CP. Modelling of fouling in a plate heat exchanger with high temperature pasteurization process. Appl Therm Eng 2021;189:116674.10.1016/j.applthermaleng.2021.116674Suche in Google Scholar

37. Stancl, J, Zitny, R. Milk fouling at direct ohmic heating. J Food Eng 2010;99:437–44. https://doi.org/10.1016/j.jfoodeng.2009.11.019.Suche in Google Scholar

38. U.S. Food and Drug Administration. Grade “A” Pasteurized Milk Ordinance, 1989 Revision; 1989 [Retrieved 09 Dec 2024].Suche in Google Scholar

39. Qin, SJ, Badgwell, TA. A survey of industrial model predictive technologies. Control Eng Pract 2003;11:733–64. https://doi.org/10.1016/s0967-0661(02)00186-7.Suche in Google Scholar

40. Helena, FA, Jorge, AWG. Continuous HTST pasteurization of liquid foods with plate heat exchangers: mathematical modeling and experimental validation using a time-temperature integrator. J Food Eng 2014;123:78–86. https://doi.org/10.1016/j.jfoodeng.2013.09.022.Suche in Google Scholar

41. Afel, B, Mohamed, NM, Mohamed, HB. Cascade control of grid connected PV systems using TLBO based fractional order PID. Int J Photoenergy 2019:1–17:4325648.10.1155/2019/4325648Suche in Google Scholar

42. Keviczky, L, Hetthessy, J. Minimum variance control-A review and outlook. Budapest: Department of Automation, dissertation, Technical University; 1976.Suche in Google Scholar

43. Panda, A, Panda, RC. Adaptive nonlinear model-based control scheme implemented on the nonlinear processes. Nonlinear Dyn 2018;91:2735–53. https://doi.org/10.1007/s11071-017-4043-7.Suche in Google Scholar

44. Panda, RC, Subramaniam, S. Design of a PID controller from a predictive control algorithm. Chem Prod Process Model 2009;4. https://doi.org/10.2202/1934-2659.1362.Suche in Google Scholar

45. De Callafon, RA, Van den Hof, PMJ. Suboptimal feedback control by a scheme of iterative identification and control design. Math Model Syst 1997;3:77–101. https://doi.org/10.1080/13873959708837050.Suche in Google Scholar

46. Wei, XZ. On least square identification of ARMAX models. Spain: 15th terminal congress; 2002.Suche in Google Scholar

47. Wang, J, Song, J, Li, Y, Ren, T, Yang, Z. State of charge estimation for lithium-ion battery based on improved online parameters identification and adaptive square root unscented Kalman filter. J Energy Storage 2024;77:109977. https://doi.org/10.1016/j.est.2023.109977.Suche in Google Scholar

48. Aguilar, JA, Cetto, JA, Husar, A. Control-oriented estimation of the exchange current density in PEM fuel cells via stochastic filtering. Int J Energy Res 2022;46:2516–22529.10.1002/er.8555Suche in Google Scholar

49. Niva, L, Yli-Korpela, A. Control of a benchmark boiler process model with DMC and QDMC. IFAC Proc 2012;45:190–5. https://doi.org/10.3182/20120902-4-fr-2032.00035.Suche in Google Scholar

50. Desborough, L, Harris, T. Performance Assessment measures for univariate feedback control. J Chem Eng 1992;70:1186–97. https://doi.org/10.1002/cjce.5450700620.Suche in Google Scholar

51. Li, L, Xia, Y. Stochastic stability of the unscented Kalman filter with intermittent observations. Automatica 2012;48:978–81. https://doi.org/10.1016/j.automatica.2012.02.014.Suche in Google Scholar

52. Li, L, Yu, D, Xia, Y, Yang, H. Stochastic stability of a modified unscented Kalman filter with stochastic nonlinearities and multiple fading measurements. J Franklin Inst 2017;354:650–67. https://doi.org/10.1016/j.jfranklin.2016.10.028.Suche in Google Scholar

53. Ghaffari, A, Tabatabaei Oreh, SH, Kazemi, R, Karbalaei, MA, R. An intelligent approach to the lateral forces usage in controlling the vehicle yaw rate. Asian J Control 2011;13:213–31. https://doi.org/10.1002/asjc.312.Suche in Google Scholar

54. Shamsuzzoha, M, Skogestad, S. The setpoint overshoot method: a simple and fast closed-loop approach for PID tuning. J Process Control 2010;20:1220–34. https://doi.org/10.1016/j.jprocont.2010.08.003.Suche in Google Scholar

55. Panda, RC. Estimation of parameters of under-damped second order plus dead time processes using relay feedback. Comput Chem Eng 2006;30:832–7. https://doi.org/10.1016/j.compchemeng.2005.12.010.Suche in Google Scholar

56. Nguyen, NH, Nguyen, PD. Overshoot and settling time assignment with PID for first-order and second-order systems. IET Control Theory Appl 2018;12:2407–16. https://doi.org/10.1049/iet-cta.2018.5076.Suche in Google Scholar


Supplementary Material

This article contains supplementary material (https://doi.org/10.1515/ijfe-2025-0038).


Received: 2024-11-14
Accepted: 2025-05-09
Published Online: 2025-05-28

© 2025 Walter de Gruyter GmbH, Berlin/Boston

Heruntergeladen am 5.10.2025 von https://www.degruyterbrill.com/document/doi/10.1515/ijfe-2025-0038/html
Button zum nach oben scrollen