Максимізація продуктивності ємнісних пастеризаторів шляхом математичної оптимізації температурних режимів

Ці розрахунки — математична модель до статті Оптимізація випоювання телят: ємнісні пастеризатори УЗМ-Е для молочних ферм

Ціль цього дослідження — знаходження абсолютного математичного мінімуму тривалості повного циклу пастеризації для ємнісних апаратів.

Класичний режим пастеризації (65 °C із витримкою 30 хвилин) є технологічним рудиментом, який суттєво обмежує продуктивність обладнання. Загальний час циклу пастеризації формується з трьох фізичних етапів:

• Часу нагрівання молока до цільової температури.
• Часу стаціонарної витримки при цій температурі.
• Часу подальшого охолодження до температури випоювання.

Мінімізація загального часу вимагає побудови комплексної математичної моделі, яка враховує термодинаміку теплообмінника та кінетику відмирання мікрофлори.

1. Математична модель складових циклу

1.1. Час стаціонарної витримки Основою для визначення необхідного часу експозиції є рівняння Дальберга-Кука, яке описує кінетику знезараження туберкульозної палички в молоці. Експоненційна залежність часу витримки від температури пастеризації на основі базових режимів (65 °C / 30 хв; 74 °C / 20 с; 85 °C / 0 с) має такий вигляд:

$$t_{витримки}(T_{пастер}) = e^{36.7459 - 0.45 \cdot T_{пастер}}$$

1.2. Час нагрівання Тривалість розгону системи залежить від підведеної теплової потужності та маси продукту. Для ємнісного пастеризатора УЗМ (потужність P = 18000 Вт, маса молока Mмол = 200 кг, середня питома теплоємність Cмол~3900 Дж/(кг·°C) функція часу нагрівання (в секундах) описується рівнянням:

$$t_{нагріву}(T_{пастер}) = \frac{(T_{пастер} - T_{мол0}) \cdot C_{мол} \cdot M_{мол}}{P}$$

1.3. Час охолодження Охолодження ємнісного апарата підпорядковується закону Ньютона-Ріхмана і має логарифмічну залежність. За результатами чисельного моделювання температурної кривої, час охолодження від температури пастеризації до температури випоювання (40 °C) крижаною водою (10 °C) описується функцією:

$$t_{охолодження}(T_{пастер}) = -\frac{1250}{31} \cdot \ln\left(\frac{T_{випоювання} - T_{води}}{T_{пастер} - T_{води}}\right)$$

2. Оптимізація та знаходження екстремуму функції

Щоб знайти найшвидший режим роботи, необхідно знайти екстремум (мінімум) сумарної функції часу, яка об'єднує всі три складові (з приведенням секунд до хвилин):

$$t_{повне}(T_{пастер}) = \frac{t_{нагріву}(T_{пастер}) + t_{витримки}(T_{пастер})}{60} + t_{охолодження}(T_{пастер})$$

Аналіз поведінки функції демонструє конфлікт двох фізичних процесів:

1. Підвищення температури пастеризації експоненційно скорочує час витримки.
2. Підвищення температури пастеризації лінійно збільшує час нагрівання та логарифмічно розтягує час подальшого охолодження.

Розрахунок у середовищі Mathcad виявляє дві критичні точки:

• Якщо оптимізувати лише час експозиції (нагрівання + витримка), ігноруючи охолодження, математичний мінімум лежить на позначці 71.508 °C.
• При врахуванні масивного "хвоста" охолодження, екстремум загальної функції зміщується вниз. Оптимальна температура пастеризації, що забезпечує найменший загальний час циклу, становить 70.047 °C.

3. Економічна доцільність: Порівняння режимів

Математично доведений оптимум (70.047 °C) кардинально перевершує застарілий класичний режим (65 °C). Зведемо розрахункові дані у порівняльну таблицю.

Параметр циклу	Режим 1: Класичний	Режим 2: Оптимізований (УЗМ)
Температура пастеризації	65 °C	70.047 °C
Час нагрівання	~44 хв	~47.6 хв
Час витримки	30 хв	~3 хв (динамічна)
Час охолодження	24.4 хв	28.3 хв
Загальний час циклу	~98.4 хв	~78.9 хв

Відмова від стаціонарної 30-хвилинної витримки та перехід на розрахований температурний оптимум заощаджує понад 19 хвилин чистого часу на кожному циклі. Враховуючи скорочення допоміжних операцій завдяки автоматизації, реальний виграш перевищує 25 хвилин. Для середнього молочного підприємства це означає можливість проведення однієї додаткової варки за зміну без жодних інвестицій у розширення парку обладнання.

Наступний інженерний крок

Знайдений оптимум температури (70.047 °C) без стаціонарної витримки беззаперечно максимізує продуктивність апарата. Проте, з точки зору мікробіології та фізичної хімії білків, цей розрахунок є неповним, оскільки молоко акумулює летальність та теплове навантаження також під час розгону та тривалого охолодження.

Чи є цей форсований швидкісний режим безпечним для продукту? Щоб довести це, необхідно перейти від простого балансу часу до інтегрального числення теплового впливу на всьому криволінійному графіку циклу. Цьому присвячена наступна стаття.