Разработка математической модели процесса сушки плодов черной смородины в вакуум-аппарате с СВЧ - энергоподводом Текст научной статьи по специальности «Физика»

Процессы и аппараты пищевых производств

УДК 51-74

Профессор С.Т. Антипов, доцент Д.А. Казарцев, доцент А.В. Журавлев, аспирант С.А. Виниченко

(Воронеж. гос. ун-т. инж. технол.) кафедра машин и аппаратов пищевых производств. тел. (473) 255-38-96 E-mail: vini-46@yandex.ru

Professor S.T. Antipov, associate Professor D.A. Kazartsev, associate Professor A.V. Zhuravlyov, post-graduate SA. Vinichenko

(Voronezh state university of engineering technologies) Department of machines and apparatuses of the food industry. Tel.: (473) 255-38-96 E-mail: vini-46@yandex.ru

Разработка математической модели процесса сушки плодов черной смородины в вакуум-аппарате с СВЧ-энергоподводом

Development of mathematical model of black currant berries drying process in vacuum device with the microwave power supply

Реферат. Разработана математическая модель процесса сушки плодов черной смородины в вакууме с СВЧ-энергоподводом, отличающаяся высокой пространственной и временной детализацией, учитывается структура отдельных плодов и слоя плодов, теплофизические параметры зависят от координаты и времени, а также с учетом изменения формы плодов и структуры слоя плодов по мере обезвоживания. Для моделирования механического поведения плодов используется метод динамики частиц, который в последние десятилетия все чаще используется в различных отраслях науки и техники. Для того чтобы модель обладала высоким пространственным разрешением, моделируемый слой плодов состоит из 20-50 отдельных плодов, при этом каждый плод, в зависимости от диаметра, разбивается приблизительно на 100 отдельных элементов. Элементы плодов по своим физическим свойствам делятся на 3 типа (кожура, мякоть, семена). Поэтому в целом, с точки пространственной детализации, слой плодов состоит из 2000-5000 элементов. Моделирование производится в двухмерном декартовом пространстве X-Z. Состояние каждого элемента-круга задается четырьмя переменными: декартовыми координатами его центра и двумя составляющими скорости. Механическое взаимодействие элементов между собой принято вязкоупругим, что позволяет учитывать основные механические свойства плодов - модуль упругости, коэффициент внутреннего трения. В рамках разработанной модели учитывается, что между соседними элементами могут возникать силы отталкивания (при внедрении элементов друг в друга), или притяжения (при отдалении сцепленных элементов друг от друга). Описание процессов тепло- и влагообмена между элементами, а также между элементами и окружающей средой базируется на основе общепринятых уравнений сушки. В модели считается, что подводимая мощность СВЧ - излучения перераспределяется в нагреваемом объеме пропорционально влажности элементов.

Summary. The mathematical model of black currant berries drying process in vacuum device with the microwave power supply, differing in high space and temporary specification is developed, the structure of separate berries and a layer of berries is considered, heat- and physical parameters depend on coordinate and time as well as in the accounting of berries form change and a berries layer structure in the dehydration process. We used the particles dynamics method for modeling of mechanical behavior of berries which is increasingly being used now in different branches of science and technology. To give the model the high space specification the modeled berries layer consists of 20-50 separate berries, each berry depending on diameter being broken into approximately 100 separate elements. Berries elements are divided into 3 types depending on the physical properties (peel, pulp, seeds). Therefore, in general, from the point of space detail, a layer of berries consists of2000-5000 elements. Modeling is carried out in two-dimensional Cartesian space X-Z. The condition of each element of circle is set by four variables: Cartesian coordinates of its center and two components of speed. Mechanical interaction of elements among themselves is accepted to be viscous and elastic that allows to consider the main mechanical properties of berries - the elasticity module, internal friction index. Within the developed model it is considered that between the neighbor elements there can be pushing away forces (at the introduction of elements into each other), or attractions (at a distance of the linked elements from each other). The description of the processes of warm and moisture exchange between the elements as well as between the elements and environment is based upon the standard equations of drying. In model it is considered that the microwave power brought is redistributed in the heated volume in proportion to elements moisture.

Ключевые слова: сушка, черная смородина, математическая модель, вакуум, тепло- массообмен, СВЧ-энергоподвод.

Keywords: drying, black currant, mathematical model, vacuum, heat and mass exchange, the microwave power supply.

© Антипов С.Т., Казарцев Д.А., Журавлев А.В., Виниченко С.А, 2014

БД Agris

7

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

Многие известные модели процесса сушки используют приближения постоянства геометрических параметров объекта сушки и постоянства его физических свойств [1-5]. Однако при описании СВЧ-сушки плодов черной смородины такие модели имели бы низкую адекватность. Поэтому в настоящей работе поставлена задача разработать математическую модель процесса сушки, обладающую высокой детализацией и высокой адекватностью, базирующейся на общепринятых уравнениях сушки, но использующей в полной мере вычислительные возможности современных компьютеров [6].

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

- изменение формы отдельного плода;

- изменение структуры слоя плодов;

- внутренний тепло- и влагообмен внутри плода;

- тепло- и влагообмен в слое плодов;

- тепло- и влагообмен с сушильным агентом.

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

Для моделирования механического поведения плодов используется метод динамики частиц, который в последние десятилетия все чаще используется в различных отраслях науки и техники [7-8].

Для того чтобы модель обладала высоким пространственным разрешением, моделируемый слой плодов состоит из 20-50 отдельных плодов, при этом каждый плод, в зависимости от диаметра, разбивается приблизительно на 100 отдельных элементов (рисунок 1). Элементы плодов по своим физическим свойствам делятся на три типа (кожура, мякоть, семена). Поэтому в целом, с точки пространственной детализации, слой плодов состоит из 2000-5000 элементов.

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

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

Рисунок 1. Плоды черной смородины (а, в) и их представление в модели (б, г): а, б - до сушки; в, г -после сушки

Рисунок 2. Слой плодов черной смородины (а, в) и его представление в модели (б, г): а, б - до сушки; в, г - после сушки

Моделирование производится в двухмерном декартовом пространстве Х-Х. Состояние каждого элемента-круга Е задается четырьмя переменными: декартовыми координатами его центра (х 1, и двумя составляющими скорости (Ух1, V21). Механическое взаимодействие элементов между собой принято вяз-коупругим, что позволяет учитывать основные механические свойства плодов - модуль упругости, коэффициент внутреннего трения.

В рамках разработанной модели учитывается, что между соседними элементами мо-

б

а

в

г

в

г

гут возникать силы отталкивания (при внедрении элементов друг в друга), или притяжения (при отдалении сцепленных элементов друг от друга) (рисунок 3).

Уравнения движения элементов составляются на основе второго закона Ньютона:

m,.

d x¡ "с/2

■ • Э

= Е

1=1 1 *

п л+х -х1)

сч- г)—-

1 2 1 к

е. + е. о, к >——1 + де„;

е + е.

+ кч(к —2Г~){У*1 - ), Г1 <"

е + е. 2

+ де

т..

■*' Э

= -т& + Е

^ + е (21 - х) .

— - г..)-— + к., (г..---)(у . - V ), г.. < —

п/ 1х- 1 г\ ■'Ун и

1=1 1 *'

1 2

е. + е.

о, к > --1

1 2

е + е..

+ де

2

е + е..

2

+ де

(1)

2

г

где . - номер элемента; и е - масса и диаметра ¡-го элемента; хь х 1 - декартовы координаты элемента; / - время; Ыэ - количество элементов; 1 - номер элемента, возможно контактирующего с .-м элементом; с ц и - коэффи-

циенты жесткости и вязкости взаимодействия элементов . и 1; гц - расстояние между центрами элементов . и 1; Vх1, V^ - декартовы составляющие скорости ¡-го элемента; е в - расстояние ограничения взаимодействия между элементами; g - ускорение свободного падения.

г и

..

а б

Рисунок 3. Силы упругого взаимодействия элементов плодов смородины при контакте: а - силы отталкивания при внедрении элементов друг в друга; б - силы притяжения при удалении элементов друг от друга

Расстояние г и между центрами элементов рассчитывается через координаты центров

по теореме Пифагора: ^ = -^(х1 - х1 )2 + (у - у1 )2.

Уравнения (1) представляют собой дифференциальные уравнения второго порядка и решаются в процессе моделирования численным методом - методом Рунге-Кутта второго порядка [9]:

хГ1 = хТ+ VI, •Д/ + аТ •(д)2 /2;

V Т

(2)

•Т+1 = VТ+ аТ • Д/;

х Т+1 = х;+ ^-Д + аТ .(Д)2 /2; VI

■Т+1 = VI + аТ • Д/,

где - номер элемента, т и т+1 - индексы текущего и следующего временного шага; А/ -шаг интегрирования по времени; хь VI, а1 - ко-

ордината, скорость, ускорение элемента.

Создание слоя плодов в модели производится следующим образом. В зависимости от заданной толщины слоя продукта к определяется количество плодов Ып, которое необходимо поместить в модельное пространство:

N

П

4 Ьк

к уп >5П

к пОП

уп П .ср

(3)

где 5с и 5п - площадь слоя плодов (в двухмерной модели рассматривается не объем продукта, а площадь) и отдельного плода; куп - коэффициент упаковки плодов для учета пористости слоя (составляет около 1,4); Ьх - длина модельного пространства; ^п.ср - средний диаметр плодов.

Плоды размещаются по всей области моделирования случайным образом по равно-

х

х

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

После того, как плоды размещены в пространстве (в виде кругов заданного диаметра), начинается заполнение кругов отдельными элементами. Сначала размещаются элементы кожуры по окружности с равным шагом, таким образом, чтобы они касались друг друга. Затем размещаются элементы мякоти внутри круга случайным образом по равномерному закону с исключением их перекрытия. Количество элементов Кэш, на которые разбивается плод 7, зависит от диаметра плода Бш:

КэШ =

£

т

Б

т

к уп £ э

Л и ^

уп Э. нач

(4)

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

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

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

Описание процессов тепло- и влагооб-мена между элементами, а также между элементами и окружающей средой базируется на основе общепринятых уравнений сушки А.В. Лыкова. Уравнения сушки капиллярно-пористого тела в векторной форме может быть записано следующим образом [5]:

с § ~е1й=?(т)-с (вш+Вт ЧТ +Qcвч ж ,т) (5)

(1 - е)— = у(рш)+ У(БТ УТ),

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

ния посредством испарения и конденсации к изменению влагосодержания за счет переноса

жидкости; У = — 7 +— к - дифференциаль-

дх дг

ный оператор набла; X - локальный коэффициент теплопроводности; С - теплоемкость жидкости; Б, Б т - коэффициенты, характеризующие данное пористое тело и определяемые эмпирически; Qсвч - теплота, выделяющаяся в среде под действием СВЧ-излучения.

Уравнения (5) для процесса СВЧ-сушки плодов черной смородины являются чрезвычайно сложными, и не допускают аналитического решения путем введения сколько-нибудь обоснованных допущений, так как большинство переменных в уравнениях являются функциями как координат, так и времени: Т(х, г, 0; Ж(х, г, 0; Сэ(х, г, 0; С(х, г, 0; е(х, г, 0; Цх, г, 0; Б(х, г, 0 ; Бт(х, г, t); Qcвч (х, г, 0. Поэтому данные уравнения целесообразно решать только численными конечно-разностными методами, с помощью предварительной дискретизации модельного пространства. Для этого плоды черной смородины разбиваются на отдельные узлы, совпадающие с центрами элементов круговой формы для механического моделирования. В этом случае каждый 7-й элемент плода характеризуется температурой Т и содержанием влаги Ж1.

Получившаяся сетка для конечно-разностных методов является нерегулярной (рисунок 4). Считается, что между двумя узлами сетки присутствует связь (возможен обмен теплом и влагой), если соответствующие элементы механически контактируют друг с другом: Гц < кк(и + ^|)/2, где кК - пороговый коэффициент контакта (равен 1,2 для изображенной на рисунке 4 сетки).

Рисунок 4. Представление слоя плодов в модели (а)

2

а

б

ВестникВТУИЖ №1, 2014_

и соответствующая сетка для решения уравнений сушки (б)

С течением времени температура и вла-госодержание каждого элемента изменяются. При этом расчет Т и Ж1 производится с использованием уравнений математической физики в конечных разностях. На каждом шаге интегрирования т рассчитывается новое значение содержания влаги и температуры элемента 7, находящегося в контакте с элементом ], по следующим итерационным формулам:

W = Wт + Dne

ko (T,Z-TK)

w :- wт

TX+1 = T *+x-

T T T T

r

J

-At,

Ai; (6)

(7)

где At - шаг интегрирования по времени; Б0 -коэффициент диффузии влаги при комнатной температуре; Лб - коэффициент экспоненциального роста коэффициента диффузии; Тк -комнатная температура (20ОС). х - коэффициент температуропроводности.

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

(

w;~ Wolp P^^At; (8)

ЖХ+1 = Жг- к

1 1 I

тг1 = т;-Хг (т-- т;кр) At, (9)

где Лв - коэффициент передачи влаги при испарении с границы; %т - коэффициент передачи температуры от окружающей среды к элементу; Жокр и Т^кр - влажность и температура газа вблизи слоя плодов (принято допущение, что данные параметры не зависят от времени); ратм - атмосферное давление; р - давление

окружающей среды.

Так как влажность Wi элементов мякоти постепенно уменьшается, в модели регулярно производится перерасчет массы и диаметра элементов:

m (W ) = mco + kW; (10)

dt (Wi ) = V kd i + kd W

(11)

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

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

(12)

P At

Tx+1 = t. х + ' = р.

i i X

cimi

WT At

IW *

где Р1 - мощность, подводимая к элементу ; с 1 - удельная теплоемкость плодов черной смородины.

Таким образом, общую систему уравнений, описывающую процесс СВЧ-сушки плодов черной смородины можно записать следующим образом:

где Г - радиус-вектор, задающий положение рассматриваемой точки.

В отличие от существующих научных работ по сушке плодов, в настоящей работе:

- уравнения сушки решаются с высокой пространственной и временной детализацией (все функции и коэффициенты зависят от положения в пространстве Г и от времени 1);

- добавлены уравнения деформации и движения плодов.

r

X

'=1

cs (r, t У-ЩА-s^T = v(a(F' t V (?, t))-c(f, t №r, t )VW (r, t) +DT (r, t )VT (r, t ))VT (r, t)+q^ (f, t);

(1 -s(r, t ))?КЫ = v(D(r, t )VW (r, t)) + V{DT (r, t )VT (r, t)); dt

(

d2 x & m-L = I

' Л! ^

dt2

d2 z.

m.-'- = -mg +1

' dt2 lS

d + d ( x - x ) d + d d + d с..(--^ - r..)—-J- + k.(r —'-J-)(v . - v.), r.. <--^ + Ad„

УУ rs у/ J V IJ rs -/v x' xx^ IJ rs В

2 rJ 2 2

d +d

0, r >--J- + AdB;

J 2 В

d + d. (z. - z ) d + d. d + d. ^

C'J (^T^ - r. У-^ + kJ (r.- v. ) rJ <-L-. + AdB 2 r 2 2

(13)

0, r >

d +d

+ Ad

11

ВестникВТУИЖ №1, 2014_

ЛИТЕРАТУРА

REFERENCES

1 Лыков А. В. Теория сушки. М.: Энергия, 1968. 472 с.

2 Лыков А.В. Тепломассообмен. М.: Энергия, 1978. 480 с.

3 Гинзбург А.С. Технология сушки пищевых продуктов. М.: Пищевая промышленность, 1976. 248с.

4 Гинзбург А.С. Расчет и проектирование сушильных установок пищевой промышленности. М.: Агропромиздат, 1985. 336 с.

5 Решетин О.Л., Орлов С.Ю. Теория переноса тепла и влаги в капиллярно-пористом теле // Журнал технической физики. 1998. Т. 68. № 2. С. 140-142.

6 Советов Б. Я. Моделирование систем: учебное пособие. М.: Высшая школа, 1998. 319 с.

7 Allen M.P., Tildesley D.J. Computer Simulation of Liquids. Oxford: Clarendon Press, 1987. 408 p.

8 Кривцов А.М. Деформирование и разрушение тел с микроструктурой. М.: ФИЗМАТЛИТ, 2007. 304 с.

9 Инженерные расчеты на ЭВМ: справочное пособие под ред. В.А. Троицкого. Л.: Машиностроение, 1979. 288 с.

1 Lykov A. V. Teoriia sushki [Drying theory]. Moscow, Energiia, 1968. 472 p. (In Russ.).

2 Lykov A.V. Teplomassoobmen [Heat and mass exchange]. Moscow, Energiia, 1978. 480 p. (In Russ.).

3 Ginzburg A.S. Tekhnologiya sushki pishche-vykh produktov [Technology of drying foodstuff]. Moscow, Agropromizdat, 1976. 248 p. (In Russ.).

4 Ginzburg A.S. Raschet i proektirovanie sushi-l'nykh ustanovok [Calculation and design of drying installations of the food industry]. Moscow, Agropromizdat, 1985. 336 p. (In Russ.).

5 Reshetin O. L., Orlov S. Iu. The theory of transfer of heat and moisture in a capillary and porous body. Zhurnal tekhnicheskoi fiziki. [Magazine of technical physics], 1998. vol. 68, no. 2, pp. 140142. (In Russ.).

6 Sovetov, B. Ia. Modelirovanie system [Modeling of systems]. Moscow, Vysshaia shkola, 1998. 319 p. (In Russ.).

7 Allen M.P. Tildesley D.J. Computer Simulation of Liquids. Oxford, Clarendon Press, 1987. 408 p.

8 Krivtsov A.M. Deformirovanie i razrushenie tel s mikrostrukturoi [Deformation and destruction of bodies with a microstructure]. Moscow, Fizmatlit, 2007. 304 p. (In Russ.).

9 Troitskii V.A. Inzhenernye raschety na EVM [Engineering calculations on the computer], Leningrad, Mashinostroenie, 1979. 288 p. (In Russ.).