ВведениеОбработка импульсным электрическим полем (ИЭП) представляет собой новый вид электрофизического воздействия на клеточном уровне в пищевой промышленности, которая требует глубокого исследования ее механизмов и эффектов. Одним из направлений использования данного вида воздействия в РФ является совершенствование процессов переработки масличных материалов [1]. В процессе обработки ИЭП происходит разрушение целостности мембран масличных клеток за счет прохождения импульсов высокого напряжения через структуру материала, что приводит к облегченному выходу внутриклеточных компонентов. Для расширения спектра применения ИЭП к различным пищевым продуктам и процессам в промышленном масштабе необходимо четкое понимание формируемых изменений внутренней структуры материалов и их реологических свойств [2]. Особый интерес представляют данные течения материала в канале шнекового транспортера, для осуществления непрерывного режима обработки импульсным электрическим полем. Процессы переработки сырья в агропромышленном комплексе имеют ряд особенностей, связанных с превалирующей долей стоимости сырья в себестоимости готовой продукции. Поэтому разработка перспективных направлений эффективного производства пищевой продукции связана с последовательными процессами многоэтапного фракционирования с выделением пищевых компонентов. Наиболее интересными в этом случае являются такие процессы фракционирования, как теплоперенос, массоперенос и фильтрация. Объединяет эти процессы наличие диффузионного и реологического механизмов переноса компонентов в рабочей зоне аппарата [3].Численный анализ этих процессов показывает значительное влияние консистенции материала на скорость процессов переноса целевых компонентов [4]. Учитывая ограничения математического моделирования, связанные с параметризацией уравнений переноса, особую актуальность приобретают физико-химические методы анализа компонентов сырья, основных и вспомогательных материалов, участвующих в этих процессах. В ряде случаев математическое моделирование позволяет получить косвенные характеристики процесса [5, 6], не наблюдаемые непосредственно. В этом случае решающим фактором, определяющим адекватность модели реальному процессу, является экспериментальное определение этих характеристик [7, 8].Целью данной работы является исследование параметров течения пластичной массы масличного материала в канале шнекового транспортера с оценкой влияния предварительной обработки импульсным электрическим полем. Объекты и методы исследованийИдентификация течения вязко-пластичного материала в канале шнекового транспортера основана на определении эффективной вязкости. Учитывая влияние давления, развиваемого в зоне межэлектродного пространства [9], ячейка вискозиметра была оборудована гидравлической системой регулируемого давления на ротор вискозиметра (рис. 1). Учитывая пластичность исследуемого материала, зазор контакта «ротор – поршень» составлял 2 мм.   поршень ротор материал ячейка    Рисунок 1 – Ротационный вискозиметрс гидравлической системой нагрузкиFigure 1 – Rotational rheometer with hydraulic loading system Выбор ротационного способа измерения основан на возможности измерения вязкости, как ньютоновских, так и неньютоновских [9], т. е. истинных и структурированных сред. Структура реологических потоков в шнековом транспортере в значительной степени определяется выбором реологического уравнения течения, которое влияет на объемную производительность экструдера. Поэтому экспериментальное исследование было направлено на определение зависимости вязкости от напряжения и скорости сдвиговой деформации, а также гидравлического давления на исследуемый материал до и после обработки импульсным электрическим полем.Учитывая, что измеряемая вязкость связана с процессами структурообразования при увеличении скорости сдвиговой деформации, то измеряемая эффективная вязкость стремится к ньютоновской.В качестве экспериментального материала использовали:мезгу, выходящую из жаровни и поступающую в экструдер;предварительно обработанную мезгу при параметрах напряженности поля E = 8 кВ/см и количестве импульсов n = 100, 200 и 300 шт., при ширине импульса 100 мкс и паузой между импульсами 900 мкс.   Режим обработки был выбран на основе данных электрических и диэлектрических параметров обрабатываемого материала с показателем удельной проводимости 3,93 ∙ 10-4 см/м [10]. Учитывая характер обработки и незначительный рост температуры (0,5–1 °С), температура принималась равной для всех экспериментов. Обработку проводили на базе ЦКП «Исследовательский центр пищевых и химических технологий» ФГБОУ ВО «КубГТУ» с использованием энергетической установки Matsusada в разработанной камере обработки с плоскопараллельными электродами [11]. Вязкость измеряли на ротационном вискозиметре Fungilab One S при избыточном гидростатическом давлении 981 Па. Высота слоя материала в измеряемой ячейке составила 38 мм. Скорость сдвиговой деформации в процессе измерения варьировалась от 1 до 10 с-1. Полученные экспериментальные данные подвергались статистической обработке в программной среде Excel и ANOVA. Результаты и их обсуждениеТипичный график зависимости коэффициента эффективной вязкости представлен в осях μэф  (эффективная вязкость) –  γ  (скорость сдвиговой деформации) (рис. 2).Как видно из представленного графика, при увеличении скорости сдвига вязкость мезги монотонно снижается, что подтверждает выдвинутое предположение о структуро-образовании мезги, которое сохраняется и для образцов после обработки ИЭП. При увеличении количества посылаемых импульсов от n = 100 до 300 шт. вязкость мезги возрастает. Далее для сравнения будут использованы данные для материала с обработкой ИЭП при n = 300 шт. ·Учитывая, что коэффициент эффективной вязкости, регистрируемый на ротационном вискозиметре, представляет собой отношение напряжения сдвига к скорости сдвиговой деформации , произвели обработку первичных реологических данных (табл. 1).  ·Для аппроксимации данных таблицы использовали линеаризацию исходных реологических показателей в обратных величинах: ; . В этом случае модельная зависимость кривой течения материала в шнековом транспортере может быть представлена в следующем виде:        (1)   где b0, b1 – коэффициенты линейной аппроксимации обратных величин кривой течения материала в канале шнекового транспортера (b0 = 0,03826 Па-1) и (b1 = 0,01009 Гц/Па). Для оценки вида этого структурообразования построили график консистенции исследованного вязкопластичного материала в осях «напряжение сдвига – скорость сдвиговой деформации» (рис. 3).  Рисунок 2 – Зависимость эффективной вязкостиот скорости сдвиговой деформации для материала дои после обработки при напряженности поля E = 8 кВ/смFigure 2 – Dependence between effective viscosity and shear strain rate for the material before and after treatmentat field intensity 8 кV/сm  τ = 0,0917 𝛾 + 24,522   Рисунок 3 – График консистенции вязкопластичного материала до обработки ИЭПFigure 3 – Viscoplastic materialthickness before treatment with pulse electric field Таблица 1 – Реология структурообразования мезги, поступающей в шнековый транспортерTable 1 – Rheology of the pulp structure formation when it enters a conveyor auger feeder Скорость сдвиговой деформации γ , с-1Напряжение сдвига τ, ПаИнженерная модель τR, ПаБингамовская модель τлин, ПаНевязка инженерной модели (τ-τR)/τ, %Невязка бингамовской модели (τ-τлин)/τ, %Доверительный интервал δτ/τ, %10,4725,4825,4925,480,0 0,0 1,8 5,2425,0024,8825,000,5 0,0 2,43,1424,0924,1124,810,1 3,0 2,2 2,0923,0923,2124,710,5 7,0 3,3 1,0520,9220,8824,620,2 17,7 9,6        τ = 0,1143 γ  + 22,889    Рисунок 4 – График консистенции вязкопластичного материала после обработки ИЭП при n = 300Figure 4 – Viscoplastic material thickness after treatmentwith pulse electric field when n = 300 Для определения эффекта влияния ИЭП на реологические параметры масличного материала построили график консистенции вязкопластичного материала после обработки импульсным электрическим полем (рис. 4).Инженерная модель реологии мезги, представленная уравнением (1) позволяет определить предельные значения реологических свойств за счет построения асимптоты напряжения сдвига τ∞, определяемой следующим соотношением:        (2)   Для уточнения параметров инженерной реологической модели (2) необходима гладкая функциональная зависимость в виде сплайн-аппроксимации, заданной в точках графика (рис. 3, 4) зависимости на отрезке [a = 1,05; b = 10,47] Гц, разбитом на части γi (табл. 1). Для аппроксимации использовали кубический сплайн дефекта 1, представляющий собой функцию, которая:•     на каждом отрезке является многочленом степени не выше третьей; •     имеет непрерывные первую и вторую производные на всем отрезке [a, b]; •     в экспериментальных точках выполняется равенство сплайна интерполирующей функции. Для однозначного задания сплайна накладываем дополнительные требования на границах сплайна: τ''(a) = τ''(b) = 0. В этом случае, согласно теореме Шенберга – Уитни об условиях существования интерполяционного сплайна, существует только один сплайн τs(γ), удовлетворяющий перечислен-ным выше условиям. В этом случае интегральная относительная невязка пробной инженерной реологической функции может быть представлена в виде целевой функции Z(b0, b1):       (3)  Минимизация функционала (3) позволила уточнить параметры инженерной модели по сравнению с их квазилинейной аппроксимацией (b0 = 0,03817 Па-1) и (b1 = 0,01052 Гц/Па) и (b0 = 0,0399 Па-1) и (b1 = 0,0156 Гц/Па) для образцов после обработки. Реологическое уравнение течения, характерное для обработки материала ИЭП в процессе извлечения масла из мезги, определяется интервалом скоростей сдвиговой деформации от 5 до 11 c-1. В этом случае для определения реологических параметров течения мезги в канале шнекового транспортера мезги наиболее реалистичным уравнением течения является идеально-пластическая модель Бингама:         (4)   где τ0 – предел текучести; μпл – коэффициент пластической вязкости. Параметры уравнения (4) могут быть определены на основе линейной аппроксимации в указанном интервале скоростей сдвиговой деформации с учетом найденной асимптоты (2) из параметров инженерноймодели (1). Из графиков линейных аппрокси-маций (рис. 4) следует, что τ0 = 24,522 Па;μпл = 0,0917 Па∙с и τ0 = 22,889 Па;μпл = 0,1143 Па∙с для образцов после обработки ИЭП соответственно. Уточняя начальное прибли-жение этих параметров для скоростей сдвиговой деформации, используемых в процессах транспортировки масличного материала в шнековом транспортере, применяли интегральную относительную невязку идеально-пластической модели Бингама (4) относительно инженерной реологической функции:       (5)   Минимизация функционала (5) позволила уточнить параметры идеально-пластической модели Бингама относительно инженерной реологической функции (τ0 = 24,3617 Па; μпл = 0,1168 Па∙с).  ВыводыУчитывая тот факт, что пластическая вязкость мезги соответствует вязкости растительногомасла [12], входящего в состав этого вязко-пластичного материала, получили подтверждение о выдвинутом предположении бингамовской реологии этого материала и сохранении ее после обработки импульсным электрическим полем. Как видно из представленных данных, график консистенции представляет собой бингамовскую вязкопластичную жидкость.Проведенные экспериментальные исследования показали наличие структурообразования материала, поступающего в шнековый транспортер. При увеличении скорости сдвиговой деформации эффективная вязкость снижается. Консистенция материала соответствует бингамовской реологии. Влияние обработки ИЭП на реологические параметры масличного материала присутствуют и подтверждают эффект изменения структуры материала после предварительной обработки со снижением предельного напряжения сдвига от 24,36 до 22,89 Па для исследуемого материала после обработки напряженностью поля E = 8 кВ/см и количеством импульсов n = 300. Снижение предель-ного напряжения сдвига за счет обработки ИЭП позволяет затрачивать меньшее количество энергии на процесс транспортировки или отжима масличного материала в случае процесса прессования. Данный фактор является положи-тельным дополнением к существующему эффекту повышения выхода масла после предварительной обработки ИЭП. Полученные данные о слоистом представлении масличного материала формата «масло – твердые частицы – масло» и его реологические характеристики являются важными в процессах моделирования современного процесса обработки пищевых продуктов импульсным электрическим полем. Полученные параметры инженерной модели позволяют прогнозировать реологию вязко-пластичного течения в широком диапазоне скоростей сдвиговой деформации в каналах шнекового транспортера.