автореферат диссертации по технологии продовольственных продуктов, 05.18.14, диссертация на тему:Холодильная технология рыб и рыбопробуктов

кандидата технических наук
Грайзи, Жамиль Афиф
город
Одесса
год
1997
специальность ВАК РФ
05.18.14
Автореферат по технологии продовольственных продуктов на тему «Холодильная технология рыб и рыбопробуктов»

Автореферат диссертации по теме "Холодильная технология рыб и рыбопробуктов"

МИНИСТЕРСТВО ОБРАЗОВАНИЯ УКРАИНЫ ОДЕССКАЯ ГОСУДАРСТВЕННАЯ АКАДЕМИЯ ХОЛОДА

На правах рукописи

ГРАЙЗИ ЖАМИЛЬ АФИФ

ХОЛОДИЛЬНАЯ ТЕХНОЛОГИЯ РЫБ И РЫБОПРОДУКТОВ

Специальность 05Л8.14-Хранение и холодильная технология

пищевых продуктов

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Одесса - 1997

Работа выполнена в Одесской Государственной Академии Холода

Научный руководитель - кандидат технических наук, профессор В.П. Оншценко

Официальные оппоненты - доктор технических наук, профессор В.А. Загоруйко - кандидат технических наук Д.Н. Ильинский

Защита состоится " апреля 1997 г. в {{ часов на заседании специализированного Совета Д.05.20.01 при Одесской Государственной Академии Холода по адресу. 270100, г. Одесса, ул. Дворянская, 1/3.

С диссертацией можно ознакомиться в библиотеке ОГАХ. Автореферат разослан марта 1997 г.

Ученый секретарь

Ведущая организация - Одесская Государственная Академия Пищевых Технологии

специализированного Совета д.т.н., профессор

В.А. Календерьян

Исх.№

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Рыба и рыбопродукты являются важными компонентами питания человека, источником биологически ценных веществ, наименее опасны с точки зрения развития онкологических и сердечно-сосудистых заболеваний. Охлаждение, замораживание, последующее холодильное хранение и транспортировка являются главными составными процессами первичной переработки рыб с целью сохранения их исходного качества, биологической ценности. Однако, наряду с разнообразием способов холодильной обработки и хранения, в настоящее время холодильная обработка и хранение недостаточно стандартизированы как по видам рыб и рыбопродуктов, так и по конечным характеристикам достигаемых температурных полей, интенсивности охлаждения или заморажяваншг, протекающих биохимических реакций в процессах холодильной обработки, хранения, последующего размораживания. Последние существенным образом определяют качество рыбопродуктов, готовых к употреблению, их роль в питании и здоровье человека.

С другой стороны, кинетика биохимических реакций определяется достигаемыми уровнями кристаллизации воды, температурными полями и тепловыми потоками на тех или иных временных интервалах. Поэтому поиск оптимальных холодильных технологий для сохранения исходного качества рыбного сырья не может быть проведен без предварительного математического моделирования процессов холодильной технологии, охватывающего все разнообразие способов и технических средств холодильной обработки, хранения, размораживания.

Имеющиеся в мировой научно-технической литературе результаты исследований Венгер К., Гинзбурга А., Ионова А., Конокотина Г., Лобзина П., Фикиина А., Чижова Г., Dickerson R., Cleland D., Kamhar В., Katayama К., Levin D., Lorentzen G. и их научных школ, других авторов отражают чаще температурные поля в теле рыб и рыбопродуктов для конкретных условий экспериментов. Обобщение экспериментальных результатов проведено только по временной длительности процессов охлаждения, замораживания, размораживания в рамках известной формулы Планка для процессов намораживания льда из чистой воды. Изменение температурных полей и тепловых потоков в процессах замораживания и размораживания по сути почти не исследовано. Необходимое в решении этих задач знание температурной зависимости теплофизических свойств (ТФС) различных рыб, динамики кристаллизации и таяния льда в широком диапазоне температур не достигнуто, а ТФС рыб и рыбопродуктов в процессах их размораживания в литературе вовсе отсутствуют. Поэтому математическое моделирование ТФС, температурных полей, тепловых потоков, доли вымороженной воды в различных рыбах и различных процессах охлаждения, замораживания и размораживания представляется необходимым на пути дальнейшего расчетно-теоре-тического анализа изменений качества рыбного сырья.

Целями и задачами настоящей работы являются: - разработка методики прогнозирования и таблиц ТФС (эффективная теплоемкость, плотность, теплопроводность, доля вымороженной воды) различных видов рыб в широком диапазоне температур и начальных влагосодержаний. При этом необходимо учесть гистерезис формирования значений доли

- вымороженной воды при осуществлении процессов замораживания-размораживания;

- разработать математическую модель процессов охлаждения, замораживания и размораживания рыб как тел неправильной геометрической формы в различных охлаждающих средах, с учетом сильной зависимости ТФС от температуры;

- разработать численный алгоритм решения уравнений математической модели, обеспечивающий вычисление температурных полей, тепловых потоков, сред-неэнтальпийных температур в различные моменты времени протекания процессов охлаждения, замораживания и размораживания;

- провести апробацию методики прогнозирования ТФС, математической модели путем сравнения расчетных величин с экспериментальными данными для отдельных видов рыб и рыбопродуктов, в различных процессах их холодильной обработки и размораживания.

Научную новизну н научные результаты работы составляют:

- для двадцати шести видов рыб и рыбопродуктов получены впервые зависимости для ТФС "сухой" компоненты в диапазоне температур от минус 40 до 30 °С, которые в сочетании с аналогичными зависимостями для воды, льда, воздуха позволили описать ТФС рыб и рыбопродуктов с погрешностью, не выше погрешностей экспериментальных данных различных авторов, как в процессах замораживания, так и в процессах размораживания;

- разработан алгоритм расчета и представлены подробные таблицы для ТФС рыб и рыбопродуктов как в процессах замораживания, так и в процессах размораживания. Они рекомендуются для практического использования в отрасли переработки рыбного сырья;

- разработана математическая модель и реализующий ее численный алгоритм для вычисления всех необходимых характеристик процессов замораживания и размораживания рыб и рыбопродуктов в различных средах. Математическая модель представляет собой нелинейную нестационарную краевую задачу теплопроводности с граничными условиями второго и третьего рода, а также с интерполяционным параметром, который отражает соотношение между теплообменной поверхностью, объемом и характерным размером рыб и рыбопродуктов;

- проведена апробация (сравнение расчетных и экспериментальных величин) математической модели по температурным полям, тепловым потокам, длительности процессов замораживания и размораживания для девяти видов рыб и рыбопродуктов, а также по длительности процессов замораживания и размораживания модельных тел (шар, конечный цилиндр, конечная пластина, яйцо, пирамида, рыба) из синтетического вещества - тилозы, строго заданного состава. Это позволяет рекомендовать разработанные математическую модель и ее алгоритм для проведения дальнейших исследований и разработки систем контроля процессов холодильной обработки и размораживания рыб и рыбопродуктов.

На основе полученных научных результатов сформулировано и обосновано научное положение:

гистерезис значений доли вымороженной воды в процессах замораживания и размораживания пищевого сырья существенным образом определяет различие в закономерностях формирования температурных полей, тепловых потоков. Дальнейшее развитие процессов холодильной технологии и систем ее контроля требует проведения прецизионных экспериментальных исследований деталей процессов кристаллизации и таяния льда в пищевом сырье. Достоверность полученных научных результатов, научного положения и рекомендаций обеспечены подробным сравнением расчетных величин и экспериментальных данных различных авторов, которые демонстрируют отклонения не выше погрешностей опытных данных. Достоверность полученных данных по ТФС подтверждена не только явным сравнением расчетных и экспериментальных значений самих ТФС, но и неявно, через их использование для расчета температурных полей, тепловых потоков, длительности процессов холодильной обработки и размораживания.

Практическая ценность результатов работы состоит в том, что получен обширный обобщающий материал для решения задач анализа протекающих биохимических реакций и формирования качественных изменений в различных видах рыб и рыбопродуктов, последующей оптимизации процессов холодильной технологии, проектирования и эксплуатации техники холодильной обработки и размораживания рыб и рыбного сырья.

Апробация работы: основные результаты исследований докладывались и обсуждались на 60-й учебно-методической и научно-технической конференции профессорско-преподавательского состава, научных сотрудников и студентов ОГАХ "Теория и практика вузовской науки" (г. Одесса, 1995 г.), на

IV международной конференции по экологии ОГАХ "Экология, продукты питания, здоровье" (г. Одесса, 1995 г.), на международной научно-технической конференции "Холод и пищевые производства" (г. Санкт-Петербург, 1996 г.), на

V научно-методической конференции ОГАХ "Людина та навколишне середовище-проблеми безперервно'1 еколопчно'1 освгги в вузах" (г. Одесса, 1996 г.), на международной конференции " Refrigeration and Aquaculture " Международного института холода (г. Париж, 1996 г.).

Публикации: по теме диссертации опубликовано две статьи в трудах международной конференции "Refrigeration and Aquaculture", проведенной Международным Институтом Холода (г. Париж, 1996 г.), три доклада в трудах Одесской Государственной Академии Холода, один доклад в трудах международной конференции Саньсг-Петербурской Государственной Академии Холода и Пищевых Технологий.

Структура и объем диссертации: диссертация состоит из введения, четырех глав, заключения, списка использованной литературы и содержит 95 стр. текста, 56 рисунков, 13 таблиц и 3 приложения на 26 страницах. Список литературы состоит из 161 источников, из них 81 иностранные.

СОДЕРЖАНИЕ РАБОТЫ

В первой главе представлен анализ современного состояния теории и практики в реализации процессов охлаждения, замораживания и размораживания рыбы, в частности, данных по ТФС рыб и рыбопродуктов, способов и режимов холодильной обработки и размораживания, их математического моделирования. На основе проведенного анализа формулируются цели и задачи диссертационного исследования.

Вторая глава посвящена восстановлению значений ТФС при охлаждении и замораживании рыб в широком диапазоне температур (-40+30 °С) по имеющимся ограниченным экспериментальным литературным данным различных исследователей. В анализе использовались данные для двадцати шести видов рыб (треска, лосось, пикша, морской окунь, тунец, скумбрия, лещ, шпрота, ставрида, килька, форель, судак, сом, сельдь и др.). Поскольку в литературе отсутствуют данные о свойствах жировой компоненты рыб, они рассматривались как объекты четырехкомпонентного состава - сухая компонента, воздух, вода, лед. При заданных зависимостях от температуры ТФС воды (включая переохлажденное состояние), льда и воздуха .получено, что для сухой компоненты всех видов рыб и рыбопродуктов имеют место зависимости:

- для удельной теплоемкости: Сс(Т) = 450 + 9.9г,

- для плотности: рс(Т) = 1330-02557*, (1)

- для теплопроводности: ^с/П = 0.131 + 2- 10~3Г.

Суммарные ТФС рассчитывались по аддитивным соотношениям вкладов соответствующих свойств компонент (массовым или по удельньм объемам), апробированным для мясного сырья, овощей, молочных продуктов: С(Т) = СсШс + СВЖВ (1 - ш) + СлГва>,

1 ^-со^иуо [ Д^ (2)

Р<т) Рс Р в Рп Рг

-1 ^М) ,, Гйа> ДГГ

Р(т) Рс Ре Рл Рг

Эффективная удельная теплоемкость рассчитывалась (соответственно источни-ковому члену уравнения теплопроводности) по формуле:

¿т""

а для равновесной доли вымороженной воды ш(Т) и ее производной по температуре наиболее удачной оказалась аппроксимация Риделя:

Для рыб и рыбопродуктов получено, что Тэ=77 К, а=1.10, у=1.12 в различных процессах замораживания. На основании полученных соотношений рассчитывалась также удельная энтальпия, объемная эффективная теплоемкость, теплота фазового перехода вода-лед, необходимые в дальнейшем анализе процессов холодильной обработки рыб.

СЭФФ(Г) = С(Г)-1Г,ЦГ)-г£г1, (3)

Сравнение расчетных н экспериментальных значений показало, что отклонения достигают по теплоемкости 6...7%, по плотности не более 4%, по теплопроводности 8...10%, но в области температур начала кристаллизации воды возрастают до 20%. Такое описание принято как удовлетворительное, поскольку экспериментальные данные различных авторов различаются с аналогичными отклонениями. Типичное согласование экспериментальных и расчетных данных по ТФС представлено на рис. 1,2.

Гит — «.м 'О О Д «.-«.ямвымися У «.М 11. ({«»•) О *»• 4» Гю*г9г А.С.

« УК

«.Я *.< 1*«»

й^.иима««!!

О*.-Ш Лм» Л.С Д^- 1Я Чттт ПК. (1ГЩ

зЁзЁ^,

10 20 30 «Г

5,%30

ь> 0ТиоОа1&°и[

а ОПРПООРРООЯ

- «.'С

-40 -33 -30 -23 -20 -13 -10 -3 0 3 10 15 20

Рис. I - Сравнение расчетных и экспериментальных данных по теплоемкости при замораживании пикши (а-сравнение абсолютных величин, Ь- относительные отклонения)

Й

<Ъ □

□ »,-«3

О

О »,-ю

А

Ж

о о

1*С

-40 -30 -20 -10 0 10 29 30 40

Рис, 2 - Сравнение расчетных и экспериментальных данных по теплопроводности при замораживании трески (а- сравнение абсолютных величин, Ь- относительные отклонения)

Третья глава посвящена математическому моделированию процессов охлаждения и замораживания рыб и рыбного сырья с целью описания длительности процессов, изменений температурных полей и тепловых потоков во времени, достигаемых средиеэитальпийяых температур. Математическая модель представляет собой нелинейную одномерную краевую задачу теплопроводности с интерполяционным параметром Г, характеризующего форму объектов холодильной обработки как тел неправильной геометрической формы:

Щ;с, г)

хт ас

ск х ¿к ди

Г(х,0) = /(*), л б [0, Л],

- г+а(т)[Гс(г)_ Г(Я) т)] = О,

<ЯГ(0, г)

&

Одномерность краевой задачи обусловлена отсутствием достаточно полных данных по локальным коэффициентам теплоотдачи, фактическая усредненность а(т) по угловым координатам поверхности объектов холодильной обработки. Зависимости ТФС от температуры соответствуют результатам второй главы.

Конечно-разностная аппроксимация краевой задачи (5) системой алгебраических уравнений проведено путем усреднения значений ТФС в окрестности е каждого узла интегрирования (общим числом до 80... 100) в рамках соотношения

F^WH^UA- (6)

Для граничных условий краевой задачи использовались соотношения для односторонней аппроксимации производных. Полученная система алгебраических уравнений решалась методом прямой и обратной прогонки, а численные значения температур в узлах конечно-разностной сетки устойчиво уточнялись в режиме "прогноз-коррекция" с заданного временного слоя на два последующих.

Интерполяционный параметр Г геометрической формы объектов холодильной обработки определялся как соотношение между объемом V, величиной теплообменной поверхности S и характерным размером R объектов

г 6f0,2b (7) Значение коэффициентов теплоотдачи а задавались численно, соответственно условиям конкретных экспериментальных режимов холодильной обработки, либо вычислялись по интерполяционным соотношениям Гнилинского и др. для условий обтекания тел неправильной геометрической формы газами или жидкостями с заданной скоростью потока

Nu = Numin+VNuL+Nufu,b; Nu)anl = 0.664 Re"2Pr"3;

Nu.^O.OSTRe11"-P[ , ,„ , ; (8)

,urb 1 +2443 Re ^Pr - ij

Numin=0.3...2

Апробация математической модели и реализующего ее алгоритма проводилась, прежде всего, при описании длительности процессов замораживания (достижения в опытах температуры "центра" объекта, равной минус 10 °С) объектов из модельного, строго заданного состава, вещества - тилозы. Результаты экспериментов Gelant! et, al, и наших расчетов приставлены в тзбл, 1. Представленные здесь значения средиеэнтальпийной температуры в конце процесса указывают насколько некорректно использовать температуру центра в качестве конечной цели процесса. Сравнение расчетных и экспериментальных профилей температурных полей, тепловых потоков, их изменение во времени проведено для девяти видов рыб. Расхождение составляет несколько градусов по температурам, Ю...15% по тепловым потокам с поверхности. Учитывая известную неопределенность положения термопар (как и расчетных координат) в объектах неправильной геометрической формы полученное согласование расчетных и

Таблнца 1 - Сравнение расчетных мени замораживания модельного вещества

и экспериментальных данных по вре-различных геометрических форм из - тилозы

№ Форма и значение Г 'ч а , яр •я ^ЗКСЯ Трасч 5 Чн Чк и"

1 Шаг». 2.0 -10.0 51 4 ?0 5 -21.| 1 оз 3.45 13.86 7051 707 -16 1

2 Шар. 2.0 -10.0 51.4 36.7 -22.5 3.55 3.53 -0.42 2931 227 -17.2

3 Шар. 2.0 -10.0 51.4 4.00 -37.9 1.83 1.87 2.57 2056 392 -27.0

4 Шар. 2.0 -10.0 51.4 30.7 -38.1 1.99 2.22 11.4 3401 394 -27.3

5 Цилиндр. 1.0 -10.0 28.0 18.1 -19.7 9.83 9.54 -2.95 1030 143 -12.6

6 Цилиндр. 1.0 -10.0 28.0 2.80 -23.6 6.77 7.00 3.4 717 187 -14.0

7 Цилиндо.1.0 -10.0 28.0 33.6 -25.0 7.88 8.36 6.1 1601 203 -14.6

8 Цилиндр. 1.0 -10.0 28.0 19.4 -27.0 7.34 7.12 3.0 1265 224 -15.4

9 Цилиндр. 1.0 -10.0 28.0 2.70 -30.6 5.82 5.66 2.75 905 279 -16.1

Ю Пилиндо. 1.0 -10.0 28.0 19,4 -33.9 5.45 6.04 -10.83 1452 318 -17.6

11 Цилиндр. 1.0 -10.0 28.0 23.9 -38.5 5.18 5.56 -7.34 1702 373 -19.2

12 Цилиндр. 1.0 -10.0 34.2 3.30 -19.8 4.51 4.78 -5.99 769 225 -11.3

13 Цилиндр. 1.0 -10.0 34.2 19.3 -27.2 3.78 4.10 -8.47 1553 360 -12.9

14 Цилиндр. 1.0 -10.0 34.2 3.10 -30.6 3.01 3.21 -6.64 1122 420 -14.9

15 Цилиндо.1.0 -10.0 34.2 32.4 -39.0 2.71 2.80 -3.32 2388 694 -17.8

16 Цилиндр.1.0 -10.0 20.3 14.9 -22.2 6.27 6.85 -9.25 742 200 -11.2

17 Цилиндо.1.0 -10.0 20.3 19.7 -34.1 4.12 4.05 1.7 1077 365 -12.6

18 Пирамида,0.5 -10.0 41.0 35.8 -25.6 3.62 3.92 8.43 2463 412 -12.8

19 Пирамида.0.5 -10.0 41.0 5.10 -29.2 2.70 2.72 0.93 1370 492 -13.5

20 Пирамида,0.5 -10.0 41.0 19.8 -38.0 2.48 2.47 -0.40 2314 686 -16.0

21 Яйцо. 1.6 -10.0 51.4 20.6 -22.6 5.76 5.67 -1.48 2105 147 -17.3

22 Яйцо. 1.6 -10.0 51.4 34.9 -29.1 4.55 4.90 7.7 3133 206 -21.8

23 Яйцо. 1.6 -10.0 51.4 19.6 -38.6 4.11 3.51 -14.48 2835 301 -28.2

24 Рыба. 1.1 -10.0 51.4 5.90 -21.9 4.13 4.15 0.48 1368 238 -15.0

2? Рыба. 1.1 -10.0 51.4 34.4 -29.1 3.65 3.90 6.99 3147 350 -18.4

26 Рыба. 1.1 -10.0 51.4 16.3 -38.5 2.81 2.72 -3.2 2704 503 -23.0

экспериментальных величин можно считать удовлетворительным. В свою очередь, полученные результаты свидетельствуют об апробации данных, полученных нами по ТФС рыб и рыбопродуктов.

Типичные примеры согласования расчетных и экспериментальных величин иллюстрируются на рис. 3-6.

Рис. 3 - Сравнение расчетных и экспериментальных Рас. 4 - Сравнение расчетных н экспериментальных ланных по температурным полям во времени данных по температурным полям во времени

для процесса замораживания тунца в рассоле для процесса замораживания форели в рассоле

Эмщмип: СсфмвкмМ т-И кг, ^ -4 м^ ( *-09*С ^ мгл)<!мле(1013м -

5!

1 ■

м 4 •*

< •<

44 ••

Рис. 5 - Сравнение расчетных и экспериментальных данных по температурным полям во времени для процесса замораживания трески в воздухе

100 ИЗ 130 175

Рис. 6 - Сравнение расчетных и экспериментальных данных по температурным полям во времени для процесса замораживания сома в воздухе

Четвертая глава посвящена исследованию процесса размораживания рыб, на протекании которого сказываются необратимые эффекты при замораживании в связях воды с "сухой" частью продукта. Интегрально они сказываются на долях вымороженной воды, существенно отличных от значений, формирующихся при замораживании (гистерезис). В итоге существенно другими становятся и ТФС водосодержагцих продуктов, начиная с эффективной теплоемкости, которая в отличие от (3), определяется соотношением

СЭФФ(Т) =С(Т)+ ЖВЦТ) . (9)

Равновесная доля вымороженной воды для рыб при их размораживании экспериментально не исследована. Поэтому в основу анализа положены экспериментальные данные по ю(Т) для говядины, полученные Роговым И.А., Кулагиным В.Н. методом ЯМР в режиме нагрева образца, замороженного в жидком азоте. Показано, что выбором значений параметров а, у формулы Риде-ля (4) такие данные удовлетворительно описываются и составляют величины: а=1.55, 7=1.13 - для говядины, а=1.65, у=1.48 - длятилозы, а=1.52, у=1.45 - для рыб. На рис. 7, 8 представлены типичные кривые для доли вымороженной воды в процессах замораживания и размораживания, а на рис. 9-12 соответствующие кривые для эффективной теплоемкости и теплопроводности. Столь сильные различия в значениях ТФС при замораживании и (прогнозируемых) размораживании требуют их тщательной апробации для описания температурных полей, длительности процессов размораживания.

Такая апробация проведена, прежде всего, на объектах различной геометрической формы (шар, пирамида, пластина, цилиндр, яйцо, рыбы) из тилозы. Использованы данные 130 опытов различных исследователей. Согласование расчетных и экспериментальных значений длительности процесса размораживания, как времени достижения "центром" температуры О °С, состав-

I— Г1 ...1 1 I-.. \

1-. 1, \

Заотчмжямгае -1175 Ргмюрапваиак ---ЩКгмгтсИН.

*

-80 -Ю -30

■29 -Ю 0

1.С

.... ... — N.

\

-«,-о.м Гтщмпм вяк - -

■ТО -60 -Я -« -Ю -20 -10

ГС

Рис. 7 - Гистерезисная зависимость доли вымороженной воды от температуры в процессе замораживания и обратном процесс« размораживания говядины

Рис. 8 - Гистерезис на* зависимость доли вымороженной воды от температуры в процессе замораживания и обратном процессе размораживания рыбы

(о 1

С^.83 «¡Ми-К

60

-w.-i.7S

-43

■Я

-20

-10

1,'С

Рис. 9 - Гистерезисная зависимость эффективной теплоемкости от температуры а процессе замораживания и обратном процессе размораживания говядины

с*»-

кЕМя-Ю 73

-а»

-20

-10

СС

Рис. 10 - Гистерезисная зависимость эффективной теплоемкости от температуры в процессе замораживания и обратном процессе размораживания рыбы

К I»

ВДмК)

1.6

М П 1.

М

116 <М

...

N

\

Заморааогмгас: -»,-«.75

• -

-Ю -ТО -60 -50 -» -Ю -29 -10 0 Ю 20

Заморгай««*?

Гимо ■ПС 1.80

\

\

-80 -ТО -60 -Я -Л) -30 -20 -10 О 10 20

Рис. 11 - Гистерезисная зависимость теплопроводности от Рис. 12 - Гистерезисная зависимость теплопроводности от температуры в процессе замораживания и обрат- температуры в процессе замораживания и обрат-

ном процессе размораживания говядины ном процессе размораживания рыбы

ляет не более 5...7%, результаты локазывшот о столь же некорректном использовании температуры "центра". Аналогичные расчеты проведены для размораживания говяжьих пластин (18 опытов) и дали такие же результаты. Проведено также сопоставление результатов наших расчетов с результатами расчетов других авторов и другими методами (включая зависимости, аппроксимирующие опытные данные). Результаты такого сопоставления представлены в табл. 2-3.

Таблица 2 - Сравнение расчетных и экспериментальных данных по времени размораживания различных геометрических форм из модельного вещества - тилозы

№ Форма и значение Г а , *р •н ^ЭКСЛ Трасч Чн

1 Шар, 2.0 51.4 -11.9 12.8 5.94 6.01 -1.18 1238 45.2 10.3

2 Шар, 2.0 51.4 -20.9 5.6 11.14 10.61 4.76 1332 16.5 4.64

3 Цилиндр, 1.0 28.0 -15.2 21.1 9.95 9.78 1.71 1002 95.5 12.8

4 Цилиндр, 1.0 28.0 -16.6 20.8 10.6 9.92 6.42 1032 94.0 12.6

5 Цилиндр, 1.0 28.0 -13.7 13.7 13.79 13.21 4.21 756 57.8 8.5

6 Цилиндр, 1.0 28.0 -17.2 13.6 16.12 13.41 16.81 850 57.3 8.4

7 Цилиндр, 1.0 28.0 -16.1 5.6 26.23 24.84 5.0 599 21.9 3.52

8 Цилиндр, 1.0 34.2 -15.7 21.1 5.13 5.19 -1.17 1242 142 12.2

9 Цилиндр, 1.0 34.2 -12.7 13.6 7.19 6.94 3.48 887 84.7 8.1

10 Цилиндр, 1.0 34.2 -12.6 8.5 9.35 9.67 -3.42 712 50.0 5.17

11 Цилиндр, 1.0 34.2 -32.5 5.6 15.17 15.47 -1.98 1288 34.0 3.38

12 Цилиндр, 1.0 20.3 -28.1 5.3 18.83 17.98 4.51 670 30.8 2.73

13 ПластинаДО 281 -22.7 14.6 6.44 5.96 7.45 9730 194 6.3

14 Пластина,0.0 281 -11.5 5.6 11.93 12.01 -0.67 4434 75.1 2.3

15 ПластинаДО 79.1 -10.3 41.3 3.45 3.21 6.96 3986 550 15.8

16 ПластинаДО 79.1 -23.5 13.2 8.43 8.05 4.51 2842 158 5.1

17 ПластинаДО 79.1 -27.4 5.4 15.99 15.22 4.82 2540 65.7 1.95

18 Пирамида,0.5 41.0 -13.5 21.0 4.52 4.61 -1.99 1397 206 9.50

19 Пирамида,0.5 41.0 -15.6 9.8 7.65 7.32 4.25 1028 87.7 4.5

20 Пирамида,0.5 41.0 -28.2 5.6 12.28 12.47 -1.55 1370 49.7 2.48

21 Яйцо, 1.6 51.4 -14.9 30.1 5.62 6.35 -13.0 2247 112 22.5

22 Яйцо, 1.6 51.4 -28.6 12.8 11.46 11.54 -0.7 2072 37.6 10.0

23 Яйцо, 1.6 51.4 -12.0 5.7 21.22 19.32 8.95 883 14.7 4.54

24 Рыба, 1.1 51.4 -13.0 14.7 8.29 8.11 2.17 1391 92.4 8.41

25 Рыба, 1.1 51.4 -13.7 10.8 10.48 10.16 3.1 1231 65.8 6.3

26 Рыба, 1.1 51.4 -28.4 5.6 18.76 16.56 12.1 1711 33.2 3.23

-И -

Таблица 3 - Суммарные отклонения между расчетными и экспериментальными данными по времени размораживания пластин, цилиндров, шаров из модельного вещества - тилозы полученные различными методами и авторами

№ Источник Форма 5min 4.Í

Все 3 формы -4.8 7.6 -0.3 2.4

1 Автор (используя экспери- Пластина -4.6 2.7 -1.1 2.1

ментальные данные Cleland Цилиндр -5.5 5.2 0.0 2.4

et.al.) Шар -4.8 7.6 0.3 2.7

Все 3 формы -9.5 9.1 -1.3 4.0

2 Cleland и др. Пластина -6.3 6.8 -0.8 3.5

(метод конечных разностей) Цилиндр -6.8 3.6 -2.2 2.4

Шар -9.5 9.1 -0.9 5.5

Все 3 формы -9.2 11.4 -0.9 4.4

3 Cleland и др. Пластина -7.0 6.5 -2.0 3.6

(метод конечных элементов) Цилиндр -6.0 4.6 -1.7 2.5

Шар -9.2 11.4 1.2 5.8

Все 3 формы -10.2 15.2 0.1 5.7

4 Cleland и др. Пластина -7.2 15.2 0.4 6.0

(с использованием формулы Цилиндр -7.7 5.2 -1.6 3.2

A.Calvelo) Шар -10.2 14.6 1.5 6.9

Все 3 формы -9.2 12.0 -0.2 4.9

5 Cleland и др. Пластина -7.0 10.2 -0.2 4.5

(с использованием формулы Цилиндр -9.2 5.9 -1.8 3.3

DJ. Cleland et. al.) Шар -9.0 12.0 1.3 6.0

Все 3 формы -12.6 14.3 0.0 5.4

6 Cleland и др. Пластина -7.0 14.3 1.7 6.0

(с использованием формулы Цилиндр -7.1 4.4 -1.1 3.1

К.Т. Pham) Шар -12.6 10.7 -0.5 6.2

Все 3 формы -10.2 10.3 0.3 5.6

7 Cleland и др. Пластина -8.8 11.1 0.1 5.2

(с использованием формулы Цилиндр -8.8 5.9 -1.3 3.9

К.Т. Phara) Шар -10.2 14.8 2.1 6.9

8 Plank Все 3 формы -28.0 50.5 6.0 21.4

9 Goodman Все 3 формы -15.3 53.7 13.3 17.6

10 Nagaoka и др. Все 3 формы 45.4 122.6 79.1 15.2

11 Mellor and Seppings Все 3 формы -36.4 30.7 -4.0 14.4

12 Pham Все 3 формы 0.6 35.0 14.7 7.2

13 Hung and Thompson Все 3 формы 14.7 134.2 92.2 26.6

14 Cálvelo Все 3 формы -12.1 17.1 0.0 6.0

15 Creed and James Все 3 формы -20.0 21.9 0.2 9.1

Температурные поля блоков рыбы при размораживании рассчитывались по методике третьей главы, но со значениями ТФС, характерных для размораживания. Сопоставление с экспериментальными данными проводилось для трески и кильки, для которых имелись экспериментальные данные. Расхождение расчетных и экспериментальных значений температур такое же, как и в случае замораживания. Это подтверждает "работоспособность" математической модели, а также подтверждает правомочность сделанного прогноза значений ТФС при размораживании. Типичное согласование расчетных и экспериментальных значений температур в размораживаемом рыбном блоке иллюстрируется на рис. 13,14.

ш

1 /

Г

Ом >и 1 В 417 ■ Д тл)> ■■ ЧП ■

1"С»

О 25 30 73 Ю0 125 150 175 200 225 250

Ряс. 13 - Сравнение расчетных и экспериментальных данных по температурным полям во времени для процесса размораживания блока кильки в условиях свободного движения воды

уа

А Л

1 г Н А

■гЫ. н № й± 1ш, пи. г31

ьт <пт

1 2 еепуе—еян*

1

20

60

100 Ш 140 160 1Ю

Рис. 14 - Сравнение расчетных и экспериментальных данных по температурным полям во времени для процесса размораживания трески в воздухе

ОБЩИЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ

Представленные научные результаты во второй, третьей и четвертой главах диссертации позволяют сделать вывод о том, что поставленные цели и задачи работы выполнены, данная работа в этом смысле является законченной научно-исследовательской работой. С другой стороны, обобщение полученных научных результатов позволяет сделать следующие выводы, которые явно не отражены в конце каждой из глав диссертации:

1. Аддитивные по массовым долям и по удельньм объемам компонент рыб и рыбных продуктов соотношения для ТФС обеспечивают приемлемое, с погрешностью не выше погрешности существующих экспериментальных данных, описание зависимостей ТФС, несмотря на включение жировой компоненты в общую "сухую" часть.

2. Газовая (воздушная) компонента рыб н рыбопродуктов существенно влияет на значения плотности и особенно теплопроводности в отрицательной области температур.

3. Гистерезисный характер формирования значений доли вымороженной воды для процессов замораживания и размораживания рыбного сырья, говядины, тилозы описан корреляционной зависимостью достаточно приемлемо в смысле согласования расчетных и экспериментальных значений ТФС, а затем и

данных по температурным полям, тепловым потокам, продолжительности процессов. Однако, в количественном отношении по доли вымороженной воды необходимы экспериментальные исследования, которые позволят уточнить специфику необратимых изменений в продуктах при их замораживании с последующим размораживанием.

4. Температурные исследования тел неправильной геометрической формы в экспериментах различных авторов содержат существенную некорректность в фиксации положения термопар (где в центре объекта?, на глубине от поверхности, но с какой стороны? и т.п.). Несмотря на это, в рамках столь же некорректной постановки, одномерная краевая задача теплопроводности с интерполяционным параметром формы позволяет достичь описания экспериментальных температурных полей объектов неправильной геометрической формы с погрешностью на уровне отклонений в 1...2 °С как для быстротекущих процессов холодильной обработки (200...300 с), так и для медленнотекущих (6...7 часов и более). Учитывая отсутствие в литературе каких-либо обобщенных зависимостей для значений локальных коэффициентов теплоотдачи при обтекании газами или жидкостями тел различной геометрической формы, представляется, что развитие разработанной здесь методики в направлении двух- и трехмерных задач теплопроводности не обеспечит существенного уточнения результатов и, по-видимому, не нужно.

5. Использование температуры "псевдоцентра" тел неправильной геометрической формы в качестве характерной конечной температуры в процессах холодильной обработки и размораживания неправомочно даже в полуэмпирических соотношениях ученых новозеландской, австралийской, аргентинской научных школ для продолжительности процессов. Полученные нами результаты показывают, что при одинаковой для различных тел конечной температуре центра их среднеэнтальпийная температура существенно различна (до 10 °С), т.е. от таких тел отводится (или подводится) существенно различное количество теплоты. Поэтому в качестве конечной температуры следует использовать значения среднеэнтальпийной температуры, а ее измерение (в системах контроля) проводить в характерных для нее точках объектов холодильной обработки или размораживания.

6. В ситуации, когда нами достигнуто полномасштабное описание процессов холодильной обработки и размораживания, возникает впечатление, что на этой основе можно приступить к решению задач поиска оптимальных режимов холодильной обработки или размораживания при использовании различных охлаждающих или нагревающих сред (газовых, жидких, твердых поверхностей). При этом оптимизация возможна по критериям энергетических затрат, в рамках технико-экономического или термо-экономического (эксергетического) критериев, в условиях "плавающих" цен на рыбу и рыбопродукты. Тем не менее ясно, что изменения качества, формирующиеся протекающими биохимическими процессами в рыбном сырье при его холодильной обработке, хранении, размораживании, являются гораздо более весомым фактором в формировании оптимизационного функционала, чем энергетические или другие. В этом плане достигнутые нами описания температурных полей при различных способах холодильной обработки и размораживания следует рассматривать как создание основы для построения

8. математических моделей протекающих биохимических процессов, основой для проведения дальнейших актуальных исследований.

9. Разработанные нами алгоритмы расчетов, с одной стороны, благодаря интерполяционному параметру формы объектов, обеспечивают получение температурных полей, тепловых потоков, среднеэнтальпийных температур в режиме реального времени, с другой стороны, реализуют расчеты на современных компьютерах достаточно быстро. Поэтому их можно рекомендовать не только для проведения технологических расчетов, использования в решении задач проектирования технических средств холодильной обработки и размораживания, но и в системах контроля работы технических средств и технологических процессов.

Ю.Конкретные виды рыб и рыбопродуктов, как и другое пищевое сырье, включая растительное, выступают в качестве неких абстрактных объектов, поскольку существенно отличаются влагосодержанием и другими характеристиками в зависимости от места лова, особенностей климата и т.п. В связи с этим не представляется возможным стандартизировать ТФС рыб по их конкретному названию. Поэтому целесообразнее стандартизировать алгоритм расчета или таблицы свойств в широком (возможном) диапазоне влагосо-держаний.

Основное содержание диссертации изложено в работах:

1. Graizi J.A., Chumak I.G., Onistchenko V.P., Terziev G.S. Mathematical Model of Various Freezing Processes of Fish and Fish-Products II Proc. Of the Conf. Of IIR, Commission C2, Bordeaux, France, March 1996. - Paris, IIR, 1996. - Vol. 2. - p. 379-386.

2. Graizi J.A., Onistchenko V.P., Chepurnenko V.P. Thermophysical Properties of Different Kinds of Fish // Proc. Of the Conf. Of IIR, Commission C2, Bordeaux, France, March 1996. - Paris, IIR, 1996. - Vol. 2. - p. 283-290.

3. Грайзи Ж.А., Кутыркин О.Ф., Оншценко В.П. Теплофизические свойства рыб и рыбного сырья,- В кн.:Тр. 60-й учебно-методической и научно-технической конференции профессорско-преподавательского состава, научных сотрудников и студентов "Теория и практика вузовской науки", (г. Одесса, 1995 г.): Тез. докл.- Одесса: ОГАХ.- 1995.-е. 14.

4. Онищенко В.П., Желиба Ю.А., Головский С.Е., Грайзи Ж.А. Научно-методические аспекты задач замораживания биологических объектов растительного и животного происхождения.- В кн.: IV международная конференция по экологии "Экология, продукты питания, здоровье", (г. Одесса, 1995 г.): Тез. докл.- Одесса: ОГАХ,- 1995,- с. 52.

5. Грайзи Ж.А., Онищенко В.П. Математическое моделирование процессов охлаждения и замораживания рыбы.- В кн.: Международная научно- техническая конференция "Холод и пищевые производства", (г. Санкт-Петербург, 1996 г.): Тез. докл.- Санкт-Петербург: СПАХиПТ,- 1996,- с. 129.

6. Грайзи Ж.А. Экологические аспекты холодильной технологии рыбы и рыбопродуктов.- В кн.: V науково-методична конференщя "Людина та навколишне середовище-проблеми безперервноУ еколопчноТ ocbíth в вузах" (г. Одесса, 1996 г.): 36. наук, праць - Одесса: ОГАХ.- 1996,- с. 129.

условные обозначения

№■ Условные обозначения Наименование величин Единица измерений

I С удельная теплоемкость кЛж/ГкгЮ

2 X теплопроводность Вт/ (м К)

3 а температуропроводность м7с

4 о плотность кг/м3

5 (В доля вымороженной воды .

6 1 удельная энтальпия Дж/кг

7 Ь скрытая теплота фазового перехода Дж/кг

8 XV относительная доля вещества в продукте -

9 тт температура К(°С)

10 г интерполяционный параметр -

И ф коэффициент формы (гомохронности) -

12 я характерный размер продукта м

13 8 площадь м2

14 V объем м3

15 ш масса кг

16 V скорость м/с

17 Р удельный тепловой поток Вт/м2

18 а коэффициент теплоотдачи Вт/ (м2 К)

19 т время ч (мин., с)

20 X текущая координата -

21 N11 число Нуссельта -

22 Яе число Рейнольдса

23 Рг число Прандтля

27 8 погрешность, (толщина продукта) %, (м)

Подстрочные индексы

1 в вода 13 н начальная

2 л лед 14 к конечная

3 г газ 15 с среда

4 ж жир 16 ср средный

5 с сухие вещества 17 п поверхность

б кр криоскопическая разность 18 Ц центр

7 эфф эффективная 19 тт минимальный

8 1 энталытийная 20 шах максим альный

9 э эвтектика 21 1ат ламинарный

10 пр продукт 22 СОПУ конвективный

11 эксп. Эксперимент 23 об общий

12 расч. Расчет 24 СТ стандартный

АНОТАЦ1Я

ГРАЙ31 ЖАМ1ЛЬ АФ1Ф. Дисертацк на здобуггя вченого ступеня кандидата техшчннх наук за спещальшстю 05.18.14 "Збереження та холодильна технолопя харчових продукпв". Одеська державна академгя холоду. Одеса, 1997 р.

- У днсертацн наведеш узагальнеш сшвв1дношекня для розрахунку теплофви-чних властивостей (теплоемшсть, густина, теплопровщшсть, ентальп1я) для двадцати шести вщцв риб (рыбопродукт), детально подаш таблищ значень цих величин у д1апазот температур вщ 233 до 303 К при значениях вологовм1сту риб вщ 0.60 до 0.85;

- розроблеш математична модель процесш заморожування риб i рибопродукпв, як нелшшна штерполяцшна одном1рна гранична задача та реамзуючий и алгоритм розрахунклв по тривалосги процеав, змшам температурних пол1в та теплових потоыв у час-1, середньоентальшйних температур. Подано результата апробацп математичноУ модел1 шляхом с ni вставления розрахункових та експериментальних даних як для заморожування рЬних вид ¡в риб, так i тш pi3Ho'i геометр ичноТ форми з модельног речовини - тшози;

- зпрогнозоваш значения теплоф!зичних властивостей риб та рибопродуюзв при ix розморожуванш, що враховують незворотш змши в продуктах пщ час заморожування;

- розроблена математична модель та приведен! результата и апробацн для опису процес^в розморожування риб та рибопродукпв.

Захищаеться наукове положения, що конкретизуе вплив г1стерезису по дол1 виморожено1 води пщ час заморожування i розморожування харчово! сировини на енергетику цих nponeciB та яисть продукт!в.

Юпочов! слова: риба та рибопродуктм, теплоф'|зич!н властивост1, охолодження, заморожування, розморожування, математичне моделювання.

summary

Graizi Jamil Afif. The dissertation for the scientific degree of the candidate of technical sciences on speciality 05.18.14 - storage and refrigerating technology of food products. Odessa State Academy of Refrigeration. Odessa, 1997.

• In the dissertation it was proposed additive correlations for calculations of thrmo-physical properties (heat capacity, density, heat conductivity, enthalpy) for 26 species of fishes (fish-products). Results of the work are accompanied by tables in a wide range of temperatures (233 - 303 K) and weight fractions (0.6 - 0.85);

• Mathematical modelling of fish and -fish produce freezing process is based on one-dimensional non-linear interpolating boundary problem of heat conductivity. The algorithm gives the opportunity to calculate temperature profile, heat flow, mean enthalpy temperature, thickness of several frozen layers in any point of time and finite time of freezing for objects with arbitrary shapes. Results calculated are compared with a big variety of experimental data on processes of freezing for different kinds of fish and tylose objects;

• Forecasted the values of thermophysical properties of fish and fish produce on processes of thawing, which take into account irreversible changes in products during the freezing processes;

• Worked out the Mathematical model and is developed results of its approbation for description of thawing processes of fish and fish produce.

Maintained is the scientific statement concretizing the influence of hysteresis in equilibrium fraction of frozen-out water o(T) during the freezing and thawing processes of food raw material on maintaining product quality and saving energy.

Key words: fish and fish produce, thermophysical properties, cooling,freezing, thawing, mathematical modeling.