Russian Federation
Russian Federation
Russian Federation
Russian Federation
Russian Federation
During cold drying of fish, the surface layers, which have lost a significant part of the water, become compacted and prevent the diffusion of water from the inner layers of the fish to the surface. In order to restore the diffusion properties of the surface layer of fish, it is necessary to create conditions for the redistribution of moisture along the thickness of the object, for example, temporarily limit external mass transfer. To do this, reduce the temperature, the speed of air movement in the drying chamber and increase its relative humidity. Due to diffusion, water from the central layers of the fish will flow to the dehydrated surface layer and moisten it. As a result, the mass-conducting properties of the surface layer will increase. The purpose of this work is to study the effect of relaxation processes on changes in the thickness of the dehydrated surface layer of fish for the subsequent development of scientifically sound rational drying regimes. The calculation of moisture redistribution in this work was carried out under the condition (assumption) that there is no external mass transfer of moisture from the surface of the processing object. The kinetics curve and water distribution curves in the processing facility were constructed during the production of dried fish products in 2; 2.5; 13.5; 55 hours respectively. The average humidity of the whiting after dehydration lasting 2 hours was 333% in terms of dry matter. The thickness of the dehydrated layer S = 1.6 mm. After 1.5 hours of relaxation, the moisture on the surface of the fish increased from 272 to 330% and approached the average humidity of the sample. During dehydration of fish lasting 13.5 hours, the average moisture content of the sample per dry mass was 217%, the thickness of the dehydrated surface layer S = 0.9 mm. After 2.5 hours of relaxation, the moisture content on the surface of the layer was 174%. It was found that at the beginning of the drying process, the dry dehydrated layer has a positive tendency to restore the lost moisture-conducting (diffusion) properties, but at the end of the drying process it takes more time to restore them. The paper proposes a method for calculating relaxation curves, which can be used to study the processes of moisture redistribution in the processing object in order to create rational drying modes.
drying process, humidity, dehydration, fish, whiting, relaxation of the dehydrated layer, dry matter, surface layer, moisture redistribution
Введение
Согласно «Стратегии развития рыбохозяйственного комплекса Российской Федерации на период до 2030 г. и плана мероприятий по ее реализации» [1], тресковые виды рыбы являются одним из основных видов рыбных ресурсов и обладают высоким потенциалом глобального рыночного спроса. В данном документе изложены плановые показатели, среди которых особый интерес представляет увеличение инвестирования средств в рыбодобывающую и рыбоперерабатывающую отрасли, а также повышение доли рынка в сегменте переработки рыб семейства тресковых. При этом одним из перспективных направлений производства является изготовление сушеной, вяленой и соленой рыбной продукции (клипфиска) [1]. Также в данном документе указано, что одним из приоритетов рыбохозяйственной отрасли должно быть повышение энергоэффективности технологических процессов в рамках проекта комплексного развития данной сферы производства [1, 2].
Вяление – традиционный способ получения продукции, пользующейся повышенным спросом. Однако изменение сырьевой базы, достаточно высокая энергоемкость и длительность процесса вяления, высокие требования к качеству выпускаемой продукции вынуждают искать пути совершенствования технологии и контроля производства [3]. Особый интерес, по мнению исследователей, представляет конвективная сушка при низких температурах, т. к. при этом можно получить продукцию наиболее высокого качества, поскольку при данных температурах нежелательные реакции, такие как окислительные процессы липидов, деструкция белков, происходят наименее активно [4–6].
Для внедрения в производство новых технологий обезвоживания необходимы технические решения, позволяющие интенсифицировать процесс сушки рыбы, сократить энергетические затраты на производство и увеличить выпуск при улучшении качества вяленой продукции. Данных результатов невозможно достичь без теоретического и экспериментального исследования тепло-массообменных процессов, происходящих в объектах, которые обладают лабильными свойствами [7, 8].
Известно, что при конвективной сушке рыбы по достижении содержания общей влаги в ней 55–60 % процесс обезвоживания значительно замедляется. Длительность первого периода сушки, протекающего с постоянной скоростью, относительно времени всего процесса обезвоживания достаточно мала [9, 10]. Изучение кинетики сушки имеет важное значение при решении вопросов интенсификации данных процессов, т. к. это позволяет оценивать скорость обезвоживания в первый и второй периоды сушки при различных режимных параметрах процесса [10, 11]. Установлено, что в процессе сушки на поверхности рыбы формируется уплотненный слой, препятствующий диффузии влаги из продукта [9], при этом оставшаяся влага в основном сосредотачивается в центре образца, и от центра к поверхности формируется градиент влажности. В итоге внешний обезвоженный слой образует изоляционный барьер, препятствующий отводу влаги из продукта. Замедление скорости сушки также связано с тем, что вода, оставшаяся в центре объекта, проходит более длительный путь к поверхности, чем приповерхностная вода в начале сушки [6, 8, 12, 13].
Для обеспечения эффективности обезвоживания сырья животного и растительного происхождения технологические параметры необходимо менять в зависимости от стадии процесса обезвоживания и, кроме этого, требуется учитывать внутренние свойства обрабатываемого сырья [14]. Например, в работе А. Ключникова [7] рассмотрена кинетика сушки продукции растительного происхождения при конвективном низкотемпературном многостадийном способе, определены этапы сушки и оптимальные параметры процесса, отмечено положительное влияние низких температур на состояние конечного продукта.
Для исследования процессов обезвоживания пищевых продуктов создаются имитационные математические модели, которые могут использоваться для изучения распределения температуры и влажности в объекте обработки и других параметров, например усадки [12]. Применение математических моделей также позволяет проектировать, оценивать производительность сушильного оборудования, а также управлять процессом и оптимизировать его [15–17]. При моделировании изменения влагосодержания по толщине рыбы в процессе сушки поиск аналитических решений крайне затруднен, т. к. коэффициент массопереноса является, в свою очередь, функцией влагосодержания [18, 19]. Установлено, что при обезвоживании коэффициент массопереноса (диффузии влаги) может уменьшаться в несколько раз, т. к. функционально связан с содержанием воды в рыбе [18]. Этот факт необходимо учитывать при моделировании процессов холодной сушки рыбы. По нашему мнению, моделирование конвекционных процессов обезвоживания может быть выполнено на основе обобщенных зависимостей.
Обобщение экспериментальных зависимостей кинетики обезвоживания провесной и вяленой рыбы проводилось на основе безразмерных чисел подобия [20]. Применение обобщенных зависимостей дает возможность использовать математические выражения для расчета коэффициентов массопереноса в зависимости от текущей влажности рыбы. Это позволяет моделировать изменение влажности по толщине объекта сушки для всего процесса. На основании литературного обзора по тематике исследования можно сделать вывод, что при моделировании процессов обезвоживания пищевых продуктов необходимо учитывать изменчивость диффузионных свойств и наличие усадки обрабатываемого сырья [9]. Моделирование процессов обезвоживания с учетом усадки обрабатываемого сырья дает более точные результаты [12]. Для построения расчетной кривой кинетики обезвоживания и последующего построения кривых распределения влаги по толщине объекта с целью изучения релаксационных процессов принято решение воспользоваться методикой, представленной ранее в работе [20]. Релаксация – это процесс перераспределения влаги в объекте обработки, который применяется с целью повышения влагопроводных свойств обезвоженного поверхностного слоя рыбы. При этом необходимо создать условия для перераспределения влаги по толщине объекта, например временно ограничить внешний массоперенос. Для этого снижают температуру, скорость движения воздуха в сушильной камере и повышают его относительную влажность. За счет диффузии вода из центральных слоев рыбы будет поступать к обезвоженному поверхностному слою и увлажнять его. В итоге массопроводные свойства поверхностного слоя будут повышаться. Релаксация способствует более эффективному даль-
нейшему процессу обезвоживания за счет увлажнения поверхностного слоя рыбы влагой центральных слоев. Представленная в работе [20] методика позволяет строить кривые кинетики обезвоживания и находить расчетные кривые распределения воды в рыбе с учетом изменения диффузионных свойств объекта, а также учитывает начальную влажность рыбы, геометрические размеры и параметры (температура, относительная влажность воздуха) технологического процесса при скорости потока теплоносителя от 2 до 3 м/с.
Цель данной работы – изучение влияния процессов перераспределения влаги (релаксации) на изменение толщины обезвоженного поверхностного слоя рыбы для последующей разработки научно обоснованных рациональных режимов сушки.
Задачи исследования:
– произвести расчет распределения влаги по толщине рыбы в процессе холодной сушки;
– осуществить расчет процессов перераспределения воды внутри рыбы при релаксации в различные моменты времени;
– определить расчетным и экспериментальным путем значения толщины сухого поверхностного слоя рыбы.
Материалы и методы
Объект исследования – путассу северная Micromesistius poutassou (A. Risso, 1827).
Начальная влажность путассу w0 = 79,4 %. Удельная поверхность путассу s / m = 0,17 м2/кг, где s – площадь рыбы, м2 ; m – масса рыбы, кг.
Параметры воздуха в сушильной камере при обезвоживании: температура tв = 27,5 °С; относительная влажность φ = 20,3 %; скорость движения v = 2 м/с.
Определение массовой доли воды выполнялось согласно ГОСТ 7636-85 [21].
Плотность структуры (твердость) обезвоженного поверхностного слоя определялась по показателю «Усилие резания» на прочностномере Food Checker (Япония), инвертор – нож.
Для построения экспериментальных кривых кинетики обезвоживания определялась влажность рыбы в расчете на сухую массу , %, в соответствующие моменты времени по формуле
, (1)
где mi - масса рыбы в определенный момент времени, кг; mc - масса сухого вещества в рыбе, кг.
Расчет массы сухого вещества в рыбе:
(2)
где m0 – масса рыбы в начальный момент времени, кг; - начальная влажность рыбы на общую массу, %.
Построение кривых распределения влаги по толщине рыбы осуществлялось с помощью численных методов расчета [20], которые базируются на решении дифференциального уравнения в частных производных диффузионного типа, например для одномерного (по пространству) дифференциального уравнения в частных производных [14]:
, (3)
где – изменение влажности в пересчете на сухое вещество wc в единицу времени t; am – коэффициент диффузии влаги;
– изменение градиента влажности; х – пространственная координата.
Построение релаксационных кривых производилось на основе решения дифференциального уравнения влагопереноса с граничными условиями третьего рода [22]:
(4)
где q(x, τ) – безразмерная влажность рыбы на расстоянии от центра пластины x для момента времени τ; n – количество равноотстоящих точек на той или иной части кривой динамики релаксации (n = 4–7); Ai – постоянный коэффициент, свой для каждого члена ряда (не зависящий ни от координат, ни от времени τ); μi – корни характеристического уравнения; Fo – критерий Фурье, Fo = Dt / R2, где D – коэффициент диффузии влаги; R – полутолщина рыбы.
Расчет влажности левой части релаксационной кривой относительно точки А производился по формуле [22]
(5)
где – значение влажности в точках поворота релаксационных кривых, j = 0 – для начальной, j = 1 – первой, …, j = n – последней релаксационной кривой; w(x, τ) – искомая влажность рыбы на расстоянии от центра пластины x для момента времени τ;
– влажность рыбы в соответствующей точке предыдущей релаксационной кривой.
Расчет влажности правой части релаксационной кривой относительно точки А производился по формуле [22]
(6)
При холодной сушке рыбы, согласно технологическим инструкциям, необходимо делать перерывы в процессе обезвоживания. Это нужно для перераспределения воды по толщине объекта тепловой обработки. Дело в том, что при температурах сушильного агента менее 40 °С интенсивность внешнего массообмена превышает интенсивность внутреннего массопереноса. При обезвоживании поверхностные слои, потерявшие значительную часть воды, уплотняются и препятствуют диффузии воды из внутренних слоев рыбы к поверхности. Установлено, что коэффициенты диффузии воды в значительной степени зависят от содержания воды в рыбе. Данные коэффициенты могут уменьшаться в несколько раз по мере сушки сырья. Для восстановления диффузионных свойств поверхностного слоя рыбы применяют специальные технологические режимы для перераспределения влаги по толщине объекта. Это позволяет более рационально использовать теплоноситель (подготовленный воздух, дымовоздушную смесь) и снижать затраты электрической энергии на процесс сушки. Для создания условий для перераспределения влаги в толще объекта необходимо на время снизить внешний массоперенос, например прекратить подачу сушильного агента в камеру или снизить температуру, скорость движения теплоносителя и повысить его относительную влажность. Однако изменение влагосодержания по толщине рыбы в процессах релаксации влаги без воздействия сушильного агента определить расчетным путем достаточно сложно. Хотя задача решается, если использовать кинетические зависимости, определенные ранее [22], а затем рассчитать и построить кривые динамики распределения влаги (рис. 1).
Рис. 1. Кривые распределения влаги по толщине рыбы: 1 – начало процесса релаксации – расчет методом сеток
на основе решения уравнения (3) [20]; 2 – продолжительность релаксации 0,5 ч – расчет методом сеток
на основе решения уравнения (3); 3 – продолжительность релаксации 0,5 ч – расчет с помощью уравнений (4)–(6);
А – точка пересечения кривых 1 и 2; B – точка пересечения касательных к кривым 1 и 2;
точки C и D характеризуют различие кривых 2 и 3; S – толщина зоны релаксации
Fig. 1. Curves of the dynamics of moisture relaxation in fish: 1 – the beginning of the relaxation process – calculation
by the grid method based on the solution of equation (3) [20]; 2 – relaxation time 0.5 hours – calculation
by the grid method based on the solution of equation (3); 3 – relaxation time 0.5 hours – calculation using equations (4)–(6); Point А is the point of intersection of curves 1 and 2; Point B is the point of intersection of the tangents to curves 1 and 2; Points C and D characterize the difference between curves 2 and 3; S is the thickness of the relaxation zone
Кривая 1 отображает распределение влаги по толщине рыбы в момент прекращения подачи сушильного агента. Кривые 1–3 пересекаются в точке А, которую назвали точкой поворота левой и правой части кривых. Релаксационные кривые 2, 3 в левой и правой части от точки А практически совпали. Отклонение одной кривой от другой в правой части от точки А составило менее 1 %. Толщина зоны релаксации S характеризует толщину слоя с низкими влагопроводными свойствами.
Таким образом, можно заключить, что расчет кривых динамики релаксации следует вести двумя методами: в начале релаксации расчет ведут численным методом, кривая 1 (см. рис. 1). После построения методом сеток кривой 1 расчет продолжают, используя уравнения (4)–(6). Если расчет процесса релаксации вести, используя только численные методы, то в дальнейшем при продолжительности релаксации более получаса расчетные релаксационные кривые будут существенно отличаться от таковых, полученных путем эксперимента.
Рассмотрим на конкретном примере изменение влагопроводных свойств при сушке рыбы.
На рис. 2 представлены экспериментальная и расчетная кривые кинетики обезвоживания путассу.
Рис. 2. Кривые кинетики обезвоживания путассу: K1 – критическая точка, характеризующая завершение первого
периода обезвоживания, когда скорость сушки максимальна, данной точке соответствует влажность рыбы wK1;
K2 – критическая точка, характеризующая существенное замедление скорости обезвоживания, которое связано
с изменением диффузионных свойств объекта сушки, данной точке соответствует влажность рыбы wK2
Fig. 2. Curves of blue whiting dehydration kinetics: K1 is a critical point characterizing the completion of the first dehydration period, when the drying rate is maximum, this point corresponds to the fish moisture content wK1;
K2 is a critical point characterizing a significant slowdown in the dehydration rate,
which is connected with a change in the diffusion properties of the drying object,
this point corresponds to the moisture content of the fish wK2
Для построения экспериментальной кривой использовались формулы (1) и (2). Для построения расчетной кривой применялась методика, представленная в работе [23]. На кривой кинетики обезвоживания путассу отмечены две критических точки, K1 и K2.
Критические точки (K1 и K2) характеризуют переход к удалению влаги, имеющей, соответственно, меньшую и большую энергии связи влаги с материалом [14]. При достижении на кривой кинетики обезвоживания второй критической влажности wK2
скорость удаления влаги значительно снижается. Это происходит из-за ухудшения влагопроводных свойств рыбы. По мере высушивания происходит усадка поверхностных слоев рыбы, в результате чего уменьшаются размеры капилляров. Это является основной причиной замедления скорости обезвоживания. Если влагопроводные свойства восстановить, то можно снизить затраты энергии на процессы копчения, вяления и сушки рыбы. На экспериментальной кривой сушки точка K2 возникла через 5 ч обезвоживания.
На рис. 3 представлены графики кривых распределения влаги по толщине (релаксации) после двухчасового обезвоживания путассу.
Рис. 3. Кривые распределения влаги по толщине рыбы после двухчасового обезвоживания:
1 – начало процесса релаксации – расчет методом сеток на основе решения уравнения (3) [20];
2 – релаксация 0,5 ч – расчет методом сеток на основе решения уравнения (3) [20];
3, 4, 5 – релаксация 0,5; 1,0; 1,5 ч соответственно, расчет с помощью уравнений (4)–(6);
6 – средняя влажность рыбы; А – точка поворота кривых динамики релаксации;
А1 – смещение точки поворота кривых динамики релаксации; В – точка пересечения
касательных к релаксационным кривым; О – значение средней влажности рыбы;
l0 – толщина слоя после релаксации; – влажность рыбы в точке поворота начальной релаксационной кривой;
– влажность рыбы в точке поворота следующей релаксационной кривой
Fig. 3. Moisture distribution curves across fish thickness after two-hour dehydration:
1 – the beginning of the relaxation process – calculation by the grid method based on the solution of equation (3) [20];
2 – relaxation for 0.5 hours – calculation by the grid method based on the solution of equation (3) [20];
3, 4, 5 – relaxation for 0.5, 1.0, 1.5 hours, respectively, calculation using equations (4)–(6);
6 – average moisture content of fish; А is the turning point of the relaxation dynamics curves;
А1 – shift of the turning point of the relaxation dynamics curves; В is the point of intersection of the tangents
to the relaxation curves; О – the value of the average moisture content of the fish;
l0 – layer thickness after relaxation; – fish moisture at the turning point of the initial relaxation curve;
– the moisture content of the fish at the turning point of the next relaxation curve
Средняя влажность путассу после двухчасового обезвоживания была достаточно высокой и составляла 333 %. Однако, несмотря на высокое среднее влагосодержание, на поверхности рыбы образовался достаточно обезвоженный слой толщиной S = 1,6 мм. После полуторачасовой релаксации влажность на поверхности рыбы увеличилась с 272 до 330 %.
В левой и правой частях образца релаксационная кривая достаточно близко расположилась к линии, отражающей среднюю влажность образца в момент завершения процесса релаксации влаги в нем.
Точка А характеризует положение точки поворота кривых динамики релаксации (см. рис. 3). По мере дальнейшей релаксации точка поворота перемещается вправо и вверх (точка А1). Влажность рыбы в точке А1 через 1,5 ч релаксации выше, чем среднеобъемная влажность образца в момент начала релаксации (линия среднеобъемной влажности), что совершенно логично, т. к. повышение влажности в левой части образца происходит за счет притока влаги из его правой части, в которой влажность выше. Значения и в уравнениях (5), (6) находят по оси ординат в точках А и А1 (см. рис. 3).
В процессе релаксации влаги сухой слой увлажняется и его граница перемещается вправо до точки А1 или далее, пока влажность в образце не выровняется до среднеобъемной влажности (см. рис. 3).
В начале релаксации толщина сухого слоя составляла S = 1,60 мм. По мере развития релаксационного процесса толщина слоя возросла до l0 = 2,23 мм. Увеличение толщины обезвоженного поверхностного слоя сопровождается и частичным восстановлением влагопроводных свойств рыбы.
Увеличение продолжительности обезвоживания перед релаксацией на 0,5 ч мало повлияло на интенсивность увлажнения сухого слоя в период релаксации. Его толщина S составила 1,45 мм, а l0 = 2,20 мм (рис. 4).
Рис. 4. Кривые распределения влаги по толщине рыбы после обезвоживания 2,5 ч:
1 – начало процесса релаксации – расчет методом сеток на основе решения уравнения (3) [20];
2 – релаксация 0,5 ч – расчет методом сеток на основе решения уравнения (3) [20];
3, 4, 5 – релаксация 0,5; 1,0; 1,5 ч соответственно – расчет с помощью уравнений (4)–(6);
6 – средняя влажность рыбы; А – точка поворота кривых динамики релаксации;
А1 – смещение точки поворота кривых динамики релаксации;
В – точка пересечения касательных к релаксационным кривым
Fig. 4. Curves of moisture distribution across the thickness of fish after dehydration for 2.5 hours:
1 – the beginning of the relaxation process – calculation by the grid method based on the solution of equation (3) [20];
2 – relaxation for 0.5 hours – calculation by the grid method based on the solution of equation (3) [20];
3, 4, 5 – relaxation for 0.5, 1.0, 1.5 hours, respectively – calculation using equations (4)–(6);
6 – average moisture content of fish; A is the turning point of the relaxation dynamics curves;
А1 – shift of the turning point of the relaxation dynamics curves;
B is the point of intersection of the tangents to the relaxation curves
Однако при обезвоживании рыбы перед релаксацией в течение 13,5 ч точка поворота кривых динамики релаксации постепенно перемещается влево (рис. 5).
Рис. 5. Кривые распределения влаги по толщине после обезвоживания 13,5 ч: 1 – начало процесса релаксации –
расчет методом сеток по уравнению (3); 2 – релаксация 0,5 ч – расчет методом сеток по уравнению (3);
3, 4, 5, 6, 7 – релаксация 0,5; 1,0; 1,5; 2,0; 2,5 ч соответственно, расчет с помощью уравнений (4)–(6);
8 – средняя влажность рыбы; А – точка поворота кривых динамики релаксации;
А1, А2 – смещение точек поворота кривых динамики релаксации (точка А2 смещается далее вправо
на релаксационной кривой при увеличении времени релаксации); В – точка пересечения касательных
к релаксационным кривым; B1 – точка пересечения касательных
к релаксационным кривым при дальнейшем увеличении времени релаксации
Fig. 5. Moisture distribution curves by thickness after dehydration for 13.5 hours:
1 – the beginning of the relaxation process – calculation by the grid method according to equation (3);
2 – relaxation for 0.5 hours - calculation by the grid method according to equation (3);
3, 4, 5, 6, 7 – relaxation for 0.5, 1.0, 1.5, 2.0, 2.5 hours, respectively, calculation using
equations (4)–(6); 8 – average moisture content of fish; A is the turning point of the relaxation dynamics curves;
A1, A2 – displacement of the turning point of the relaxation dynamics curves (point А2 shifts further to the right on the relaxation curve as the relaxation time increases); B is the point of intersection of the tangents to the relaxation curves;
B1 – the point of intersection of tangents to relaxation curves with a further increase in relaxation time
Средняя влажность образца составляла 217 %. В начале релаксации толщина обезвоженного поверхностного слоя S составляла 0,925 мм, через 2,5 ч релаксации толщина обезвоженного поверхностного слоя увеличилась до l0 = 1,18 мм. Влажность на поверхности слоя составила 174 %.
Если путассу непрерывно обезвоживать перед релаксацией более 55 ч, то толщина сухого слоя составит S = 0,89 мм, а точка поворота совпадет с толщиной сухого слоя даже при непрерывной последующей релаксации в течение 4 ч. Поэтому длительное обезвоживание перед релаксацией не способствует снижению энергетических затрат на сушку рыбы. Относительная погрешность расчетного значения толщины сухого слоя и найденного экспериментальным путем составила 14 %.
В таблице представлены экспериментальные данные по изменению показателей обезвоженного поверхностного слоя путассу (толщины и твердости структуры) при температуре сушильного агента 25 °С при конечной влажности рыбы на общую массу от 62 до 64 %.
Показатели обезвоженного слоя путассу (толщина, твердость)
в зависимости от относительной влажности сушильного агента при температуре 25 °С
Indicators of the dehydrated layer of blue whiting (thickness, hardness)
depending on the relative humidity of the drying agent at a temperature of 25°C
Относительная влажность |
Толщина |
Усилие резания |
Темп обезвоживания, |
20 |
1,20 ± 0,15 |
440 ± 151 |
29,0 |
40 |
1,10 ± 0,12 |
264 ± 62 |
21,1 |
50 |
1,03 ± 0,12 |
166 ± 53 |
14,7 |
Толщина слоя незначительно уменьшается при сушке теплоносителем с достаточно высокой относительной влажностью. Следует отметить, что усилие резания возрастает при сушке теплоносителем с низкой относительной влажностью. Это связанос более интенсивным обезвоживанием. Для темпа обезвоживания 29,0 %/ч усилие резания принимало наибольшее значение из рассмотренных вариантов. Данный факт свидетельствует об увеличении твердости обезвоженного поверхностного слоя при использовании интенсивных режимов сушки. Для данных процессов конвективной холодной сушки рыбы целесообразно использовать режимы релаксации с целью рационального использования теплоносителя и экономии энергоресурсов. Так, во время релаксации нагревательные элементы сушильной установки выключаются и в установку подается теплоноситель температурой на 4–6 °С ниже температуры непосредственного обезвоживания, относительная влажность теплоносителя от 60 до 80 % (экспериментально установлено, что при температуре 20 °С и относительной влажности более 70 % выраженный сухой поверхностный слой не образуется). Для осуществления процесса релаксации возможно использование смешанного свежего уличного и цехового воздуха.
Время включения в процесс сушки, по нашему мнению, должно совпадать с образованием на кривой кинетики сушки критической точки K2 (см. рис. 2).
В дальнейшей работе планируется расчетным и экспериментальным путем обосновать периодичность и длительность применения фаз релаксации в процессе сушки рыбы.
Расчетным путем получены графические зависимости распределения влажности в рыбе при температуре процесса сушки 27,5 °С. При продолжительности процесса сушки 2 ч средняя влажность рыбы в пересчете на сухое вещество принимала значение 333 %. Для этого момента времени определена толщина поверхностного обезвоженного слоя – 1,6 мм. Влажность на поверхности этого слоя составила 272 % в расчете на сухое вещество. После релаксации в течение 1,5 ч влажность на поверхности рыбы составила 330 % в пересчете на сухое вещество, при этом толщина обезвоженного поверхностного слоя увеличилась до 2,23 мм. Об эффективности процесса перераспределения влаги можно судить, если сравнить значения влажности на поверхности и средней влажности образца. Эффективность перераспределения влаги для длительности обезвоживания в течение 2 ч имела максимальное значение.
Средняя влажность путассу после обезвоживания продолжительностью 2,5 ч составила 323 % в расчете на сухое вещество. Толщина обезвоженного слоя – 1,45 мм. Через 1,5 ч релаксации влажность на поверхности рыбы увеличилась с 243 до 315 % в расчете на сухое вещество, а толщина обезвоженного слоя возросла до 2,2 мм. Данные результаты релаксации сопоставимы с результатами при обезвоживании в течение 2 ч. В этих случаях наблюдается достаточно быстрое повышение влажности обезвоженного слоя. Также следует отметить существенное увеличение толщины обезвоженного поверхностного слоя после релаксации примерно на 0,7 мм. Это также свидетельствует о том, что слой на данном этапе сушки способен к насыщению водой из центральных слоев рыбы, он не имеет твердой структуры.
При продолжительности процесса сушки в течение 13,5 ч средняя влажность рыбы имеет значение 217 % в расчете на сухое вещество. Для этого момента времени обезвоженный поверхностный слой имеет толщину 0,9 мм. При достаточно длительной релаксации (2,5 ч) толщина поверхностного слоя увеличилась на 0,3 мм, а влажность на поверхности этого слоя достигла 174 % в расчете на сухое вещество.
При непрерывном обезвоживании путассу перед релаксацией более 55 ч толщина поверхностного обезвоженного слоя составила 0,89 мм. Толщина сухого слоя остается неизменной при перераспределении влаги в течение 4 ч. Результаты моделирования толщины обезвоженного слоя сопоставимы со значениями, найденными экспериментальным путем.
Обезвоженный слой имеет низкие влагопроводные свойства и сдерживает процесс сушки в целом, поэтому актуальным остается поиск технологических режимов обезвоживания, снижающих возникновение на поверхности рыбы слоя с низкими диффузионными свойствами. Для повышения эффективности процесса вяления рыб необходимо рационально применять режимы перераспределения воды в процессе сушки и не допускать образования плотного обезвоженного поверхностного слоя. Длительное обезвоживание перед релаксацией не способствует снижению энергетических затрат на сушку рыбы. В последующих работах планируется провести дополнительные исследования по влиянию длительности релаксации на энергоэффективность процесса сушки.
В работе В. А. Наумова, А. Э. Суслова, Ю. А. Фатыхова [18] найдены схожие зависимости распределения влажности по толщине рыбы путем моделирования. Эти зависимости получены для безразмерного влагосодержания при различных значениях числа подобия Био [18]. Однако в настоящей работе данные кривые имеют размерность для толщины
и влажности рыбы. В исследовании L. Mayor, A. M. Sereno [12] установлено, что пищевые продукты могут подвергаться усадке в процессе сушки. Усадка пищевых материалов увеличивается с увеличением объема удаленной воды. Чем больше удалено воды, тем больше напряжений сжатия возникает в материале [12]. В данной работе моделирование обезвоженного слоя осуществлялось без учета влияния усадки всего высушиваемого образца. Однако уменьшение толщины и влажности обезвоженного поверхностного слоя, найденные в настоящей работе, свидетельствует о его уплотнении и усадке по мере процесса сушки. В работе N. N. Potter, J. Hotchkiss [13] сообщается, что усадка материала происходит неравномерно по его объему, т. к. влага из продукции не удаляется равномерно, а также пищевые продукты редко обладают идеальной эластичностью. Исследователями А. Э. Сусло-
вым, Ю. А. Фатыховым [9] установлено, что повышение температуры воздуха до 30 °С при обезвоживании окуня обеспечивает интенсификацию процесса, но также приводит к ухудшению условий сушки, т. к. пересушивание поверхностных слоев способствует образованию плотного поверхностного слоя, препятствующего термовлагодиффузии. В настоящей работе найдена толщина сухого обезвоженного слоя расчетным методом. Получены экспериментальные значения толщины обезвоженного слоя при температуре теплоносителя 25 °С и относительной влажности от 20 до 50 %.
Предложен метод расчета релаксационных кривых, который может применяться для исследования процессов перераспределения влаги в объекте обработки с целью создания рациональных режимов сушки. Определена толщина сухого слоя расчетным и экспериментальным путем. В результате анализа полученных расчетных релаксационных кривых распределения влаги по толщине рыбы можно сделать вывод о том, что применение процессов релаксации целесообразно после завершения периода постоянной скорости сушки. В начальный период сушки обезвоженный слой достаточно быстро насыщается влагой из центральных слоев, при этом толщина его увеличивается. К концу процесса сушки обезвоженный слой уплотняется, а эффективность релаксации снижается. Также установлено, что на формирование обезвоженного слоя рыбы
в значительной степени влияют интенсивные режимы сушки.
1. Ob utverzhdenii strategii razvitiia rybokhoziaistven-nogo kompleksa Rossiiskoi Federatsii na period do 2030 goda i plan deistvii po ee realizatsii: Postanovlenie Pravi-tel'stva RF ot 26 noiabria 2019 g. № 2798-r [On approval of the strategy for the development of the fisheries complex of the Russian Federation for the period up to 2030 and the action plan for its implementation: Decree of the Government of the Russian Federation dated November 26, 2019 No. 2798-r]. Available at: //www.garant.ru/products/ipo/prime/doc/72972854/#1000.%2010%20oct.%202022 (accessed: 04.11.2023).
2. Skul'skaia L. V., Shirokova T. K. O problemakh ry-bokhoziaistvennogo kompleksa Rossii [About the problems of the Russian fisheries complex]. Ekonomika i biznes: teoriia i praktika, 2022, no. 1-1 (83), pp. 193-203. DOI:https://doi.org/10.24412/2411-0450-2022-1-183-193-203.
3. Bessmertnaia I. A., Ermakova Iu. A. Sovershenstvo-vanie tekhnologii vialenoi ryboproduktsii [Improving the technology of dried fish products]. Innovatsii v tekhnologii produktov zdorovogo pitaniia: materialy Mezhdunarodnoy nauchnoy konferentsii (Kaliningrad, 26 maia 2016 g.):. Kaliningrad, Izd-vo Kaliningrad. GTU, 2016. Pp. 54-60.
4. Davis L. et al. Lipid oxidation in salted dried fish: The effect of temperature and light on the rate of oxidation of a fish oil. Journal of the Science of Food and Agriculture, 1993, vol. 62, no. 4, pp. 355-359.
5. Raghunath M. R., Sankar Th. V., Ammu K., Devadasan K. Biochemical and nutritional changes in fish proteins during drying. Journal of the Science of Food and Agriculture, 1995, vol. 67, no. 2, pp. 197-204.
6. Siddhnath Kumar, Amit Ranjan, Bimal Prasanna Mohanty, Pooja Saklani. Dry fish and its contribution towards food and nutritional security. Food Reviews International, 2022, vol. 38, no. 4, pp. 508-536.
7. Kliuchnikov A. Development of new method of drying at energy-saving universal dryer to improve quality of crops used in fodder production. Proceedings of the 18th International Scientific Conference “Engineering for Rural Development” (Jelgava, Latvia, 22–24 May 2019). Latvia University of Life Sciences and Technologies, 2019. Pp. 105-111.
8. Boziaris I. S. Seafood processing: Technology, quality and safety. John Wiley & Sons, 2014. 512 p.
9. Suslov A. E., Fatykhov Iu. A. Rezul'taty issledovaniia protsessa sushki okunia [The results of the study of the dry-ing process of perch]. Izvestiia KGTU, 2011, no. 23, pp. 133-140.
10. Suslov A. E., Fatykhov Iu. A. Rezul'taty issledovaniia protsessa sushki ryby (plotva i chekhon') [The results of the study of the drying process of fish (roach and chehon)]. Protsessy i apparaty pishchevykh proizvodstv, 2012, no. 1, p. 51.
11. Jain D., Pathare P. B. Study the drying kinetics of open sun drying of fish. Journal of food Engineering, 2007, vol. 78, no. 4, pp. 1315-1319.
12. Mayor L., Sereno A. M. Modelling shrinkage during convective drying of food materials: a review. Journal of food engineering, 2004, vol. 61, no. 3, pp. 373-386.
13. Potter N. N., Hotchkiss J. H. Food science. Springer Science & Business Media, 2012. Pp. 200-211.
14. Lykov A. V. Teoriia sushki [Drying theory]. Mos-cow, Energiia Publ., 1968. 372 p.
15. Castro A. M., Mayorga E. Y., Moreno F. L. Mathe-matical modelling of convective drying of fruits: A review. Journal of food engineering, 2018, vol. 223, pp. 152-167.
16. Elmesery H. S., Tolba N. M., Kamel R. M. Mathe-matical modelling and performance analysis of airflow distribution systems inside convection hot-air dryers. Alexandria Engineering Journal, 2023, vol. 62, pp. 237-256.
17. Heilporn C., Debaste F., Boey C., Haut B., van der Pol F., Nonclercq A. Implementation of a rational drying process for fish conservation. Food Security, 2010, vol. 2, pp. 71-80.
18. Naumov V. A., Suslov A. E., Fatykhov Iu. A. Gra-nichnye usloviia pri konvektivnoi sushke ryby [Boundary conditions for convective drying of fish]. Vestnik MGTU, 2015, vol. 18, no. 4, pp. 647-653.
19. Naumov V. A., Suslov A. E., Fatykhov Iu. A. O vybore ratsional'nykh parametrov protsessa konvektivnoi sushki presnovodnykh vidov ryby [On the choice of rational parameters of the process of convective drying of freshwater fish species]. Tekhnologii pishchevoi i pererabatyvaiushchei promyshlennosti APK – produkty zdorovogo pitaniia, 2017, no. 2 (16), pp. 59-64.
20. Ershov M. A., Ershov A. M., Votinov M. V. Mass transfer processes modeling in the production of dried and air-dried fish. IOP Conference Series: Earth and Environmental Science. IOP Publishing, 2022. Vol. 1052. No. 1. P. 012067.
21. GOST 7636-85. Ryba, morskie mlekopitaiushchie, morskie bespozvonochnye i produkty ikh pererabotki. Metody analiza [ISS 7636-85. Fish, marine mammals, marine invertebrates and products of their processing. Methods of analysis]. Moscow, Standartinform Publ., 2010. 121 p.
22. Ershov M. A. The relaxation processes calculation of fish dehydrated surface layer during drying and smoking. IOP Conference Series: Earth and Environmental Science. IOP Publishing, 2019. Vol. 302. No. 1. P. 012021.
23. Ershov A. M., Ershov M. A., Nikolaenko O. A. Metodika rascheta krivykh kinetiki i dinamiki obez-vozhivaniia v protsessakh vialeniia i kholodnogo kopcheniia ryby [Method of calculation of kinetics curves and dynamics of dehydration in the processes of drying and cold smoking of fish]. Vestnik MGTU, 2010, vol. 13, no. 4-2, pp. 947-950.