Расчет процессов релаксации при обезвоживании рыбы Текст научной статьи по специальности «Механика и машиностроение»

ВестпикВТУИТ/Proceedings of VSUET, № 3, 2016L

Оригинальная статья/Original article_

УДК 664.956

DOI: http://doi .org/10.20914/2310-1202-2016-3-37-43

Расчет процессов релаксации при обезвоживании рыбы

Михаил А. Ершов, Александр М. Ершов. Владимир А. Гроховский,

[email protected] [email protected] [email protected]

1 кафедра технологий пищевых производств, Мурманск. гос. техн. ун-т, ул. Спортивная, 13, г. Мурманск, 183010, Россия Реферат. В данной работе разработаны методы изучения релаксационных процессов, применяемых при обезвоживании рыбы, с учетом влияния химического состава, геометрических размеров объектов обработки и режимных параметров процессов сушки. Использование численных методов расчета на основе решений дифференциального уравнения второго порядка для пластины с граничными условиями третьего рода позволяет с достаточной для инженерной практики точностью моделировать различные условия ведения процесса обезвоживания с периодическим восстановлением влагопроводных свойств объектов сушки. Предлагаемые расчетные методы применимы для процессов обезвоживания, состоящих из непрерывного начального периода и последующих комбинированных периодов сушки рыбы и релаксации поверхностного слоя объекта обезвоживания. Необходимость применения релаксации обусловлена тем, что по мере обезвоживания поверхностные слои, потерявшие часть влаги, уплотняются. Уменьшаются размеры капилляров для прохода влаги через поверхностные слои. Вблизи поверхности появляется зона, свободная от подавляющей массы влаги, а, следовательно, имеющая низкие диффузионные свойства. В результате замедляется процесс обезвоживания всего образца. Применение релаксации позволяет восстановить влагопроводные свойства поверхностного слоя рыбы. Во время релаксации прекращается подача электрической энергии на нагревательные элементы. Снижается скорость циркуляции сушильного агента. В сушильную установку подается воздух более низкой температуры и более высокой относительной влажности, чем сушильный агент. В сушильной установке создаются условия, сдерживающие внешний массообмен и способствующие к перераспределению влаги по толщине рыбы. Во время релаксации влага постепенно перемещается от центральных слоев, где обезвоживание еще не наступило, к обезвоженным поверхностным слоям.

Ключевые слова: обезвоживание, п

[роцещрелаксациЯ^КоЭФФииииниЫииЫУЗ

ии влаги, влагопр

роводностьмассоотдач^

The calculation of the relaxation processes during dehydration of fish

Mikhail A. Ershov, Aleksandr M. Ershov, Vladimir A. Grokhovskii,

[email protected] [email protected] [email protected]

1 department of technology offood manufactures Russia, Murmansk state technical university, Sportivnaya st., 13, Murmansk, 183010, Russia Summary. The methods of study of relaxation processes that are applicable during fish dehydration considering the influence of the chemical composition, geometric size of the processing object and regime parameters of the dehydration processes are developed in this work. The usage of the numerical methods of calculation on the basis of solution of the second-order differential equation with boundary conditions of the third kind allow with sufficient accuracy for engineering practice modeling the different conditions of the dehydration process conducting with periodic restoration of hydraulic conductivity properties of the objects of dehydration. The proposed calculation methods applicable to dehydration processes. Which consists of a continuous initial period and the combined periods of drying fish and relaxation dehydration facility. During the relaxation is provided an exposure to a drying agent for dewatering for a certain time. During relaxation stops the supply of electric power to the heating elements. Reduces the rate of circulation of the drying agent. In the drying installation is supplied with air of lower temperature and higher relative humidity than the drying agent. In drier conditions are created that constrain external mass transfer and promotes relaxation of moisture, that is, to its redistribution in the thickness of the fish. During the relaxation of the moisture is gradually shifting from the central layers where dehydration has not yet come to the dehydrated surface layers. The appearance of moisture inside the dehydrated surface area leads to putting up and expansion of capillaries. At the next interval change product moisture re-enters the dehydration process, the high conductive properties throughout its volume. Keywords: dehydration, process, relaxation, moisture diffusion coefficient, hydraulic conductivity

Введение

На интенсивность процесса массопере-носа оказывают влияние режимные параметры, геометрические размеры тела, его химический состав. Исследователями установлено [1], что скорость движения сушильного агента выше двух метров в секунду не оказывает влияния на интенсивность обезвоживания после достижения первой критической влажности на кривых кинетики сушки. Это явление можно объяснить тем, что процесс испарения влаги с поверхности в сушильную среду протекает с большей скоростью, чем перенос влаги из центра к поверхности

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

Для цитирования

Ершов М. А., Ершов А. М., Гроховский В.А. Расчет процессов релаксации при обезвоживании рыбы // Вестник ВГУИТ. 2016. № 3. С. 37-43. doi:10.20914/2310-1202-2016-3-37-43

For citation

Ershov M.A., Ershov A.M., Grohovskii V. A. The calculation of the relaxation processes during dehydration of fish. Vestnik VSUET [Proceedings of VSUET]. 2016. no. 3. pp. 37-43. (in Russian). doi:10.20914/2310-1202-2016-3-37-43

Весттк&ТУМТ/Фгосе^^ о/'ТЯНЕТ, № 3, 2016

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

1.1Цель и задачи исследований

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

1. Разработать методику изучения процесса релаксации.

2. Выявить основные закономерности, протекающие при перераспределении влаги внутри объекта обезвоживания.

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

1.2 Описание процесса релаксации

Количество влаги йМ\, перемещающейся из внутренних слоёв тела к поверхности за время йх, можно записать, воспользовавшись законом Фурье. Перенос влаги из центра тела к его поверхности равен:

ЛМ1 = —ЛтгЛ¥ (дм / дп\ сИ, (1)

где Хш\ - удельный коэффициент массопровод-ности, кг / (м-%-с); ¥ - поверхность тела, м2; w - содержание влаги в теле, %; т - продолжительность, с; (дw/дh)\ - производная по направлению, перпендикулярному к поверхности ¥ (градиент влажности), % / м.

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

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

ЛМ2 = —Лт2Л¥ (дм / дп)2 Л, (2)

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

Определенное количество влаги предается от поверхности тела к окружающей среде:

СМ3 = рт ( - )Л) (3)

где рт - удельный коэффициент массоотдачи от поверхности тела в окружающую среду, кг / (м2-%-с); Wп - содержание воды на поверхности тела, %; wв - содержание воды в сушильном агенте, %.

При отсутствии движения сушильного агента в период релаксации коэффициент рт существенно уменьшается, стремясь к нулю. Тогда йМз = 0. Принимая во внимание, что йМ\ = йМ2, получим следующее выражение:

—Лтс1¥(дм / дп) && —Лт2Л¥(дм / дп)2 Л (4)

Разделим обе части уравнения на dFdт, —Лт (( /дп))—Лт2 (дм / дп) (5)

Для удобства анализа перейдем к одной пространственной координате и будем рассматривать случай неограниченной пластины толщиной 21. Из уравнения (5) следует, что в процессе перераспределения влаги релаксационные кривые имеют общую точку с одинаковой влажностью, в которой они пересекаются. Общая точка с постоянной влажностью после достаточно глубокого обезвоживания перед релаксацией располагается вблизи внутренней границы «сухого слоя». Коэффициент удельной влагопроводности 1т2 по мере насыщения поверхностного слоя влагой увеличивается более значимо, чем уменьшается \т\. В таком случае общая точка будет постепенно перемещаться вверх и вправо на кривых динамики обезвоживания от точки А до точки А1 (рисунок 1). Разницу в динамике изменения Хт\ и 1т2 можно объяснить тем, что объем приповерхностного сухого слоя значительно меньше объема слоя с высокими влагопроводными свойствами. Поворачиваясь вокруг общей точки одинаковой влажности, левая часть кривых динамики релаксации поднимается вверх, а правая часть опускается вниз. Точку одинаковой влажности на кривых динамики релаксации назовем точкой поворота. Через довольно продолжительное время обе части достигнут среднеобъемной влажности.

1.3 Особенности релаксации

расчета процессов

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

и-1

Zwi —

i=1 n

(6)

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

образца; w - среднеобъёмная влажность, на значение которой оказывают влияние следующие параметры:

w = iy(r,N,Xp, w0, s / m), (7)

где N - скорость обезвоживания в первый период сушки; Xp - жесткость режима процесса обезвоживания; W0 - начальная влажность объекта обработки; s/m - удельная поверхность объекта обработки.

Условием расчёта процесса релаксации является наличие кривой динамики сушки в момент прекращения подачи сушильного агента. Однако кривые динамики релаксации объекта сушки, полученные расчетом методом сеток, в приповерхностой зоне имеют выпуклость вниз. Это показывает, что влажность на поверхности растёт несколько быстрей, чем в приповерхностной зоне, т. е. в данном случае допущение dM3 = 0 приводит к несколько искажённым результатам при расчетах процесса релаксации с использованием численных методов.

Рисунок 1. Кривые динамики релаксации рыбы: 1 - кривая на момент начала релаксации; 2 - кривая на момент релаксации 0,5 часа; 3 - кривая на момент релаксации 1 час; 4 - кривая на момент релаксации 1,5 часа; 5 - кривая на момент релаксации 2 часа; 6 - среднее влагосодержание в образце рыбы

Figure 1. The curves of the dynamics relaxation fish: 1 - the curve at the beginning of the relaxation; 2 - relaxation curve at the time of 0,5 hours, 3 - relaxation curve at the time of 1 hour; 4 - curve at the time of relaxation of 1,5 hours; 5 - curve at the time of relaxation 2 hours; 6 - average moisture level in a sample of fish

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

Для уменьшения погрешности определения влажности на поверхности объекта релаксации предпочтительно использовать данные расчёта численными методами только для нахождения двух кривых динамики распределения влаги непосредственно после обезвоживания, т. е. в начальный момент релаксации при Тр = 0 и в заключительный момент этого периода релаксации при Тр = Т1, например Тр = 0,5 часа. В дальнейших расчетах применяем решение дифференциального уравнения второго порядка с граничными условиями третьего рода для неограниченной пластины, которое выглядит следующим образом:

0

(x,)

= Z4

e-M'Fo -cos(щx/R)

(8)

здесь п - количество равноотстоящих точек на той или иной части кривой динамики релаксации (п = 4 - 7); Л, - постоянный коэффициент, свой для каждого члена ряда (не зависящий ни от координат, ни от времени т); Ео - критерий Фурье, Ео = Вт / Я2, где В - коэффициент диффузии влаги; ц - корни характеристического уравнения:

ctg^ = — Bi

(9)

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

При расчете левой половины кривой динамики релаксации Я = S, где S толщина зоны релаксации. Для правой части кривой динамики релаксации от точек А и А1 значение Я = I, здесь I половина толщины образца (рисунок 1). При Ео > 0,5 достаточно воспользоваться только одним корнем ц выражения (8). Использование выражения (8) для расчёта левой половины от точки А кривой динамики в процессе релаксации:

0( x ,т) W°-w0 - w (x T) w (Hl)

- w (x T) 0

w, - Wi Л (Hl)

w - w

(10)

w - w,

где w0, w1, wn - значение влажности в точках поворота для начальной, первой и последующих групп релаксационных кривых (рисунок 1), пересекающихся в одной общей для группы кривых точке: В, В1, Вп1; w(x,т) - искомая влажность образца в какой-либо его точке через промежуток времени, равный т = Тр; w(ш) - влажность образца в соответствующей точке предыдущей релаксационной кривой.

Для правой части от точки А (рисунок 1) кривой динамики релаксации искомое значение компоненты в безразмерном виде выглядит следующим образом:

0( X,) =

w - wn

(x,) 0

w, .> - w'

(rn) 0

\x,,)

(11)

w - w,

. (x,) 1

w, ., - w0

(Hl) 1

Л x,)

w - wn

(x,) n

w, . > - w'

( Hl ) n

Из уравнения (8) следует, что для релаксационных процессов распределения влажности по толщине «сухого слоя» для разных моментов времени т свойственны симметричные кривые, монотонно возрастающие от наружной поверхности «сухого слоя» образца к его внутренней поверхности (рисунок 1). Касательные, проведенные к кривым динамики релаксации, проходят через одну точку В. Если касательные к кривым пересекаются, то эти кривые должны сойтись в какой-то точке А.

Исходными данными для расчета служили левые и правые части кривых 1 и 2 (рисунок 1). Эти кривые получены расчётным путем с использованием метода сеток. Точки на кривой 1 являлись исходными для нахождения Л, и fit в уравнении (8).

Для расчёта релаксации рыбы в приповерхностном «сухом слое» (слева от точки А) использовалась величина S в качестве определяющего размера. Критерий Fo находили для левой и правой частей релаксационных кривых из выражения: Fo = Dt/S2. Определение коэффициентов по-тенциалопроводности D (коэффициентов диффузии влаги) производили по методике, разработанной авторами ранее [4].

i =1

1 Общие точки пересечения релаксационных кривых В2, Вз,..., Вп

на рис.1 не показаны. Они характерны для релаксационных процессов после достаточно глубокого обезвоживания рыбы.

Значение D определяли для приповерхностного «сухого слоя». Подбирая значения Ai и fit, добивались близкого совпадения значений, полученых из уравнения (8), с таковыми на первой релаксационной кривой 2 (рисунок 1) при Тр = 0,5 часа. В этом случае сохраняются в дальнейших расчётах граничные условия третьего рода, которые были учтены в расчетах процессов обезвоживания численными методами. Продолжительность релаксации Тр = 0,5 часа является минимально-возможной в реальных промышленных установках с точки зрения рационального управления системой автоматического поддержания заданных технологических режимов. При достаточно длительном обезвоживании пересечение кривых 1 и 2 (точка А) находится вблизи внутренней границы приповерхностного «сухого слоя». Далее строили касательную к левой части кривой 1. На касательной отмечали точку В в месте пересечения касательной и вертикальной линии, проходящей через точку А. Касательные, проведенные к двум - трем релаксационным кривым, пересекаются в точке В, а кривые 1, 2 и 3 проходят через точку А (точка поворота релаксационных кривых). Касательные, проведенные к кривым 4 и 5, имеют уже другую точку пересечения В1 и, как следствие, сами кривые пересеклись в точке А1, которая находится правее и выше точки А.

Значение влажностей w'o и w\ в точках В и В1 (рисунок 1 ) являются начальными значениями при определении безразмерной величины 0. При расчете процесса релаксации длительностью 2,5 - 3 часа возможно, что количество точек «В» возрастет до трех или четырех. В этом случае в расчетах 0 меньшее значения влажности w'o в точке В заменяют на большое в точках В1 или В2, а, возможно, и В3. Следовательно, точка поворота перемещается вправо и вверх к точкам А2, А3 (на рисунке 1 не показаны).

В реальности точка поворота А по мере релаксации плавно перемещается вправо и вверх (от точки А до точки А1). Чем меньше отношение Wc/Wx2 (отношение средней влажности к влажности второй критической точки) или wq/wo (отношение средней влажности к начальной влажности рыбы), тем меньше значение ординаты точки А1 отличается от таковой в точке А и наоборот. При больших значениях wcp/wK2 или w^/wo может наблюдаться существенное удаление точки А1 от точки А. Это явление объясняется тем, что при сравнительно высоких влажностях по мере релаксации быстрее восстанавливаются влагопроводные свойства так называемого приповерхностного сухого слоя, и процесс релаксации происходит при меньших

значениях Дw = w'o - w/, постепенно достигая равновесного состояния как в правой, так и в левой от точки А частях кривой динамики релаксации.

При бесконечно большом времени релаксации кривые динамики релаксации достигнут средней влажности образца. Тангенс угла а наклона касательных (рисунок 1) равен:

B1C да tgа =-=--

* FC дх Х=1

С другой стороны, продолжив касательную влево от точки ¥, очевидно, что

tgа = tgq>, tgф= Оу - мв) / Я,

где w/ - содержание воды на поверхности рыбы (рисунок 1).

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

-Лт (дм / 8х) = Ът - К )

или

дм / дх = tgа = tgф = - мв) /Я = Ът / Лт (м} — мв). Отсюда следует, что Я = Рт/^т. Таким образом:

dw дх

= (Wf -we)/(Хт/вт) = (Wf -we)• Bi/l (12)

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

Рассмотрим вариант, в котором большая часть влаги остается в левой части образца для восстановления влагопроводных свойств последнего, а часть - испаряется в окружающую среду. Обозначим долю влаги, испарившейся с поверхности: а = йМз / йМ2. В этом случае основное уравнение массоотдачи в дифференциальном виде примет вид:

аЛт2(дм / 8х) = Ът - ^в ) (13)

После небольших преобразований получим следующее выражение:

dw дх

= ( Wf - We )/(a • ^m2 jвт ) = --(wf - we )•Bi/(a •l)

(14)

Последнее выражение можно использовать для расчета критерия Ы, когда значение коэффициента влагообмена Рт > 0.

На рисунке 2 представлен расчет процесса релаксации в соответствии с методическим подходом, изложенным выше. Результаты

х=1

расчета с использованием численных методов практически совпадают с аналитическим решением дифференциального уравнения второго порядка с граничными условиями третьего рода.

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

Рисунок 2. Кривые динамики обезвоживания и релаксации камбалы-ерша: 1 - кривая на момент начала релаксации; 2 - кривая на момент релаксации 0,5 часа; 3 - кривая на момент релаксации 1 час; 4 - кривая на момент релаксации 1,5 часа; 5 - кривая на момент релаксации 2 часа; 6 - кривая на момент релаксации 0,5 часа, построенная методом сеток; 7 - среднее влагосодержание в образце рыбы

Figure 2. The curves of the dynamics of dehydration and relaxation rough dab: 1 - the curve at the beginning of the relaxation; 2 - relaxation curve at the time of 0,5 hours, 3 - relaxation curve at the time of 1 hr; 4 - curve at the time of relaxation of 1,5 hours; 5 - curve at the time of relaxation 2 hours; 6 is a curve at the time of 0,5 hours of relaxation constructed by numerical method; 7 - average moisture level in a sample of fish

Заключение

Заложены теоретические основы изучения механизма релаксации приповерхностного сухого слоя рыбы посредством анализа кривых динамики распределения влаги, полученных

ЛИТЕРАТУРА

1 Глазунов Ю.Т. и др. Процессы сушки, копчения и вяления рыбы: учеб. пособие для студентов образоват. орг. высш. образования, обучающихся по направлению подгот. 15.03.02 (151000) «Технологические машины и оборудование» уровня бакалавриата. М.: Моркнига, 2015. 261 с.

2 Ortiz J., Lemus-Mondaca R., Vega-Galvez A., Ah-Hen K. et al. Influence of air-drying temperature on drying kinetics, color, firmness and biochemical characteristics of Atlantic salmon (Salmosalar L.) fillets // Food Chemistry. 2013. V. 139. P. 162-169.

расчетным путем с использованием численных методов и аналитических решений дифференциального уравнения второго порядка для неограниченной пластины с граничными условиями третьего рода.

3 Martins M.G., Martins D.E.G, Pena R.S. Drying kinetics and hygroscopic behavior of pirarucu (Ara-paima gigas) fillet with different salt contents // LWT -Food Science and Technology. 2015. V. 62. P. 144-151.

4 Ершов А.М., Ершов М.А., Никола-енко О.А. Методика расчета кривых кинетики и динамики обезвоживания в процессах вяления и холодного копчения рыбы // Вестник МГТУ. 2010. Т. 13. № 4. Ч. 2. С. 947-950.

5 Похольченко В.А., Ершов А.М., Ершов М.А. Особенности обобщения подобных явлений в процессах тепловой обработки рыбы // Вестник ВГУИТ. 2015. № 2 (64). С. 44-50.

REFERENCES

1 Yu.T. Glazunov Yu.T. et al. Protsessy sushki, kop-cheniya i vyaleniya ryby [Drying processes, smoked and dried fish: the textbook for students of educational institutions of higher education, in the direction of 03.15.02 (151000) "Technological machines and equipment" level undergraduate] Moscow, MORKNIGA, 2015. 261 p. (in Russian).

2 Ortiz J., Lemus-Mondaca R., Vega-Galvez A., Ah-Hen K. et al. Influence of air-drying temperature on drying kinetics, colour, firmness and biochemical characteristics of Atlantic salmon (Salmosalar L.) fillets. Food Chemistry, 2013, vol. 139, pp. 162-169.

СВЕДЕНИЯ ОБ АВТОРАХ

Михаил А. Ершов канд. техн. наук, ст. науч. сотр., кафедра технологий пищевых производств, Мурман. гос. техн. ун-т, ул. Спортивная, 13, г. Мурманск, 183010, Россия, [email protected] Александр М. Ершов д-р. техн. наук, профессор, гл. науч. сотр., кафедра технологий пищевых производств, Мурман. гос. техн. ун-т, ул. Спортивная, 13, г. Мурманск, 183010, Россия, [email protected] Владимир А. Гроховский д-р. техн. наук, профессор, заведующий кафедрой технологий пищевых производств, Мурман. гос. техн. ун-т, ул. Спортивная, 13, г. Мурманск, 183010, Россия, [email protected]

КРИТЕРИЙ АВТОРСТВА

Михаил А. Ершов написал рукопись, корректировал рукопись до подачи в редакцию, несет ответственность за плагиат

Александр М. Ершов консультация по расчетной части исследования

Владимир А. Гроховский консультация по тематике исследования

КОНФЛИКТ ИНТЕРЕСОВ

3 Martins M.G., Martins D.E.G, Pena R.S. Drying kinetics and hygroscopic behavior of pirarucu (Arapaima gigas) fillet with different salt contents. LWT - Food Science and Technology. 2015, vol. 62, pp. 144-151.

4 Ershov A.M., Ershov M.A., Nikolaenko O.A. Methods of calculating the curves of the kinetics and dynamics of dehydration in the process of drying and cold smoked fish. Vestnik MGTU [Proceedings of MSTU] 2010, vol. 13, no. 4, part 2, pp. 947-950. (in Russian).

5 Pokholchenko V.A., Ershov A.M., Ershov M.A. Features generalization of such phenomena in the heat treatment process fish. Vestnik VGUIT [Proceedings of VSUET] 2015, no. 2 (64), pp. 44-50. (in Russian).

INFORMATION ABOUT AUTHORS

Michael A. Ershov Ph.D., senior researcher, Department of Food Production Technology, Murman. state. tehn. University Press, str. Sport, 13, Murmansk, 183010, Russia, [email protected]

Aleksandr M. Ershov Doctor of Technical Sciences, Professor, Chief Researcher, Department of Food Production Technology, Murman. state. tehn. University Press, str. Sport, 13, Murmansk, 183010, Russia, [email protected] Vladimir A. Grokhovskii Dr. tehn. Sciences, Professor, Head of the Department of Food Production Technology, Murman. state. tehn. University Press, str. Sport, 13, Murmansk, 183010, Russia, [email protected]

CONTRIBUTION

Michael A. Ershov wrote the manuscript, corrected manuscript before the submission to the editorial board and is responsible for plagiarism

Aleksandr M. Ershov consultation calculated of the study

Vladimir A. Grokhovskii consultation on the subject of research

CONFLICT OF INTEREST

Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.

ПОСТУПИЛА 03.08.2016 RECEIVED 8.3.2016

ПРИНЯТА В ПЕЧАТЬ 22.08.2016 ACCEPTED 8.22.2016