CC BY
20
ВестнщФТУМт:/Proceedings of VSUET ISSN 2226-910X E-ISSN 2310-1202
DOI: http://doi.org/1Q.2Q914/231Q-1202-2Q21-1-62-7Q_Оригинальная статья/Research article
УДК 664.1.054_Open Access Available online at vestnik-vsuet.ru
Моделирование процесса роста кристаллов сахарозы в сахарсодержащем растворе
Евгений В. Семенов 1 Анатолий А. Славянский 1 Вера А. Грибкова 1 Дарья П. Митрошина 1 Сергей Т. Антипов_2
[email protected] 0000-0001-6203-1783 anatoliy445 [email protected] 0000-0002-0262-8841 [email protected] 0000-0002-5079-8132 [email protected] 0000-0002-8724-3368 [email protected]_0000-0003-4299-1538
1 Московский государственный университет технологий и управления им. К.Г. Разумовского, ул. Талалихина, 31, г. Москва, 109004, Россия
2 Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия Аннотация. Кристаллизация сахарозы представляет собой один из наиболее важных этапов в технологическом потоке получения товарного сахара. По своей природе - это развивающийся в достаточно широком температурном диапазоне, условиях изменяющейся концентрации раствора, его пептизации и другого массообменный процесс образования и роста кристаллов сахара в сахарном растворе путем диффузионного переноса в нем молекул сахарозы. К этому нужно добавить еще и такие осложняющие физико-математическую формализацию явления, влияющие на кристаллизацию, как форма кристалла, чистота, вязкостные свойства и поверхностное натяжение раствора, наличие в нем несахаров и т.п. Поэтому среди обусловливающих в целом протекание процесса факторов обычно проводят, сохраняя среди них лишь приоритетные и наиболее важные с позиций исследования. С точки зрения классической диффузионной теории явление кристаллизации трактуется как молекулярный с бесконечной скоростью распространения возмущения от источника диффузионный перенос в растворе молекул сахарозы, вследствие чего концентрация в сахарсодержащей жидкостной системе также полагается изменяющейся мгновенно в каждой точке. Однако, поскольку в реальных условиях этот эффект не наблюдается, то, с целью разрешить это противоречие, вводят понятия возмущенного и невозмущенного, разделенных называемым диффузионным фронтом концентрации областей этих фазовых состояний. Используя кинематические характеристики фронта, определяют время протекания процесса кристаллизации в целом по обеим фазам. Что позволяет с большей, по сравнению с полученными
на базе классической теорией результатами, точностью рассчитать период обработки раствора._
Ключевые слова: вакуум-аппарат, концентрация раствора, затравка, кинетика кристаллизации, сахароза
Simulation of the growth of sucrose crystals in a sugar-containing solution
Evgeny V. Semenov 1 Anatoly A. Slavyansky 1 Vera A. Gribkova 1 Daria P. Mitroshina 1 Sergey T. Antipov 2
[email protected] 0000-0001-6203-1783 anatoliy445 [email protected] 0000-0002-0262-8841 [email protected] 0000-0002-5079-8132 [email protected] 0000-0002-8724-3368 [email protected]_0000-0003-4299-1538
1 Moscow State University of Technology and Management. K.G. Razumovsky, st. Talalikhina, 31, Moscow, 109004, Russia
2 Voronezh State University of Engineering Technologies, Revolution Av., 19 Voronezh, 394036, Russia_
Abstract. Crystallization is one of the most important steps in the technological flow of commodity sugar. By its nature, it develops in a fairly wide temperature range, conditions of varying concentration of the solution, its peptization and other mass-exchange process of formation and growth of sugar crystals in a sugar solution by diffusion transfer in it sucrose molecules. To this it is necessary to add such complicating physical and mathematical formalization phenomena, affecting crystallization, such as the shape of the crystal, purity, viscous properties and surface tension of the solution, the presence of non-sugars, etc. Therefore, among the factors that determine the whole process, the factors usually take place, keeping among them only priority and most important from the point of view of the study. The article, unlike known analogues, provides justification and quantitative analysis of the kinetics of the sucrose crystallization process, taking into account the features of the crystal shape and diffusion properties of the solution. From the point of view of classical diffusion theory, the phenomenon of crystallization is interpreted as molecular with infinite rate of perilation from the source of diffusion transfer in the solution of sucrose molecules, so that the concentration in the sugar-containing liquid system is also supposed to change instantly at each point. However, since in the real world this effect is not observed, in order to resolve this contradiction, introduce the concepts of the indignant and unperturbed, separated by the diffusion front of the concentration of the regions of these phase states. Using the kinematic characteristics of the front, determine the time of the crystallization process in general in both phases. This allows us to calculate the processing period of the solution with
greater accuracy than the results obtained on the basis of classical theory._
Keywords: vacuum apparatus, solution concentration, seed, crystallization kinetics, sucrose
Для цитирования
Семенов Е.В., Славянский А.А., Грибкова В.А., Митрошина Д.П., Антипов С.Т. Моделирование процесса роста кристаллов сахарозы в сахарсодержащем растворе // Вестник ВГУИТ. 2021. Т. 83. № 1. С. 62-70. doi: 10.20914/2310-1202-2021-1 -62-70
© 2Q21, Семенов Е.В. и др. / Semenov E.V. et al.
For citation
Semenov E.V., Slavyansky A.A., GribkovaV.A., Mitroshina D.P. Antipov S.T. Simulation of the growth of sucrose crystals in a sugar-1 containing solution. Vestnik VGUIT [Proceedings of VSUET]. 2021. vol. 83. no. 1. pp. 62-70. (in Russian). doi:10.20914/2310-1202-2021-
1-62-70_
This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License
Семенов Е.В. и др.Вестник.ФГУИШ, 2021, Т. 83, №. 1, С- 62-Введение
Хотя в настоящее время известно четыре научные гипотезы кристаллизации сахарозы, в число которых входит термодинамическая, диффузионная, молекулярно-кинетическая и дислокационная, по нашему мнению, наиболее простой и эффективной для прогнозирования протекания процесса кристаллизации является диффузионная гипотеза П.М. Силина [1-6]. В ней кристаллизация определяется как двухстадийное физическое явление: сначала инициация в чистом сахарсодержащем растворе в виде центров кристаллизации (первая стадия процесса), затем их эволюция - рост кристаллов (вторая стадия). По завершении интродукции в пересыщенный раствор затравочного материала, представляющего собой первую стадию процесса, сам процесс кристаллизации (вторая стадия) осуществляется при температуре Т > 50 °С в условиях градиента концентрации сахарозы в растворе на границе раздела фаз: "жидкость (раствор сахарозы) + твердое (кристалл)", когда молекулы сахарозы адсорбируются на поверхности кристалла и встраиваются в кристаллическую решетку.
Следует отметить, что современная диффузионная теория процесса кристаллизации не учитывает влияния формы кристаллов на динамику кристаллообразования. Авторы статьи считают наиболее близкой к истине концепцию П.М. Силина [1], полагавшего, что на 1 м2 поверхности кристаллов конденсируется одинаковое количество сахара независимо от их размера.
Стоит обратить внимание на то, что в процессе кристаллизации условия роста кристаллов изменяются беспрерывно. Так, уменьшение объема и толщины слоя межкристального пространства, а также его концентрации наблюдается при росте размеров кристаллов. Помимо этого, вследствие имеющего место возрастания температуры в вакуум-аппарате в результате конденсации, интенсифицируется и собственно диффузионный перенос молекул сахарозы к поверхности кристалла.
Обзор посвященных научно-техническому обоснованию проблемы кристаллообразования сахарозы в вакуум-аппарате зарубежных литературных источников, с одной стороны, показывает, что по-прежнему, усилия специалистов-сахарников в большой степени направлены на изучение традиционной для сахарного производства задачи количественного анализа, в условиях изменения размера кристалла, процесса кристаллообразования в вакуум-аппарате. Так, работы [8, 9] посвящены оценке кинетики кристаллизации сахарозы на основе полученных в кристаллизаторах опытных данных при равномерно-монотонном
убывании температуры раствора. Для чего при измерении размера кристалла был применен оптический микроскоп совместно со сканирующим электронным микроскопом типа SEM. При этом рассчитанные результаты сравниваются с полученными в других работах данными.
С другой стороны, особенность литературных источников, с очевидностью, указывает на ускорение перехода научного сообщества в области сахарного производства к использованию в исследованиях углубленного научного анализа, основанного на проведении опытов с использованием новейших образцов в области экспериментального оборудования [10-20]. Это связано с тем, что на практике в растворах, в том числе, и сахарсодержащих, явление кристаллообразования часто протекает в условиях, при которых, увеличивая свой размер, кристалл по габитусу отклоняется от своей классической формы в чистом растворе, т. е. становится дефектообразным. Данный фактор во многих исследованиях служит предметом пристального внимания ученых. Причем, поскольку этот процесс может инициироваться и выступающими в качестве центров кристаллизации имеющими размеры порядка одного нанометра молекулами сахарозы, то в статьях [10, 11] были изучены темпы роста кристаллизации сахарозы в чистых растворах при температуре в кристаллизаторе 313 К и пересыщении 0,094-0,181. При этом полученный кинетический коэффициент (коэффициент диффузии) изменился при начальном перенасыщении с D = 2,79* 10-7 до конечного D = 9,45 х10-5 м2/с, соответствующий размер кристалла варьировался от 1,46 *10-9до7,47* 10-9 м.
С теоретической позиции концепция трактовки явления кристаллизации, в том числе, и сахарозы как молекулярного диффузионного переноса массы (или, аналогично, явления переноса тепловой энергии в классической теории теплопроводности) скорость распространения возмущения от источника массы или энергии принимается бесконечной, вследствие чего концентрация раствора или температура среды изменяются мгновенно в каждой его точке. Чтобы учесть этот фактор вводят понятие «диффузионного фронта» x =р(0 (рисунок 1). При этом для раствора на фронте осуществляются условия:
1) концентрация с равна концентрации невозмущенной области;
2) производная дс/дк на фронте полагается равной нулю, что влечет за собой плавное непрерывное распределение концентрации, а также потока концентрации в возмущенной области и во всей исследуемой области.
Semenov E..V. et aC.Proceedings of VSUET, 2021, voC. 83, no. 1, pp. Принятые допущения означают, что процесс распространения концентрации осуществляется в две фазы.
Сахарсодержащий р-р Sugar-containing solution
X, £
Кристалл сахарозы | Sucroseciystal
Рисунок 1. Схема к расчету процесса кристаллизации сахарозы в вакуум-аппарате
Figurel.Scheme for calculating the process of crystallization of sucrose in a vacuum apparatus
Ниже, согласно предложенной П.М. Силиным диффузионной модели, с учетом перечисленных особенностей кинетики кристаллообразования в условиях молекулярной конденсации, проводится аналитический и численный анализ процесса кристаллизации сахарозы в вакуум-аппарате.
Постановка задачи
Ставится задача: при заданных геометрических и физико-механических параметрах процесса кристаллообразования получить, с учетом фактора наличия фронта концентрации, количественную оценку по времени эволюции роста размера кристаллов сахарозы, их объемного содержания при обессахаривании межкристального раствора в вакуум-аппарате.
Предварительно, с целью обосновать среднее значение h расстояния между частицами сахарозы исходя из величины объемной концентрации их в растворе принимается, что если V4c и Vym - соответственно, общий объем кристаллов и объем утфеля в вакуум-аппарате, ^ - объемная концентрация кристаллов в томже объеме, то в соответствии с принятой геометрической моделью жидкостной системы половина расстояния h между двумя соседними частицами составит [6]
h = 5/(2^1/3), 0 1, (1)
где Учс / Уут, 5 - приведенный диаметр частицы сахарозы.
Считают, что массоперенос к частице из жидкой среды осуществляется симметрично и, кроме того, для выбранной расчетной модели с удовлетворительной точностью выполняются все предположения относительно физико-механических свойств исследуемого объекта.
62-70 [email protected]
Поток концентрации q (объем вещества, отводимого с единицы поверхности в единицу времени) целевой субстанции (сахарозы) в направлении оси х в принятой одномерной модели массопереноса рассчитывают по формуле согласно первому закону Фика
= -D 8-С,
дх
(2)
где D - коэффициент диффузии, с - объемная концентрация сахарозы.
При формализации физической модели предполагают, что концентрация сахарозы в исходном растворе одинакова по всему объему, центры кристаллизации распределены равномерно, кристаллы в процессе своего роста имеют форму, близкую прямоугольному параллелепипеду и сохраняют приблизительно схожие размеры. Кроме того, дополнительно считают, что начальный пересыщенный объем раствора включает как центры кристаллизации, достаточно большое количество молекул и частиц этого продукта настолько малым размером, что эти частицы диффундировали к поверхности кристалла и впоследствии встраивались в кристаллическую решетку.
В дальнейшем полагают, что процесс кристаллообразования сахарозы обусловлен ростом кристалла лишь в направлении большей из его граней, то есть, в направлении оси х. Тогда в качестве приближенной геометрической модели кристалла принимается полупространство х > h, ограниченное поверхностью х = h, к которому из области 0 <х < h раствора диффундируют молекулы данного вещества (рисунок 1).
Более того, считают, что исследуемый процесс развивается, если концентрация сахарозы превышает растворимость ее или так называемую равновесную концентрацию, т. е. при наличии разницы в концентрации растворенного вещества в пересыщенном растворе.
Тогда с целью количественного анализа задачи о росте кристалла в соответствии с законом сохранения массы приходят к кинетическому уравнению [1]
8с d 8 2c
8t ~ 8х2''
(3)
иде t - время, D - полагаемый постоянной величиной коэффициент диффузии.
В таком случае, если Сп - исходная концентрация сахарозы в пересыщенном растворе, то решение уравнения (3) согласуют с начальным условием
с(х, 0) =Сп 0 <х (4)
Семенов Е.В. и др.Вестник.ФТУИТ, 2021, Т. 83, №. 1, С- 62-70
В свою очередь, поскольку посередине расстояния х = 0 между двумя соседними кристаллами концентрация раствора сахарозы достигает максимальное значение, то граничное условие имеет форму 8с (0, г)
8х
= 0,0 <г <<х>.
(5)
В соответствии с концепцией Л.П. Хворо-вой и В.А. Коваленко [4, 5] через определенный промежуток времени по мере роста кристалла в растворе в контактирующем с поверхностью кристалла слое раствора скорость осаждения вещества определяется законом (2), что согласуется с представлением о наличии в растворе фронта концентрации.
Тогда граничное условие на границе раздела раствор-кристалл х = 0 имеет следующий вид 8с(0, г)
8х
= у[с(0, г) - сн ]
или
(6)
8с(0, г) 8х
= а[с(0, г) - Сн ]
где а =у/В, у - так называемый коэффициент кристаллохимической реакции, характеризующий скорость поверхностных процессов, м/с, а - удельный коэффициент диффузии, м-1, Сн -концентрация насыщенного раствора сахарозы вблизи кристалла.
Решение задачи
В дальнейшем, с целью упрощения проводимого исследования вводят безразмерные величины по зависимостям
с - сп х е Вг
-- = и,— = С,,—Г = х,
Ас Ь Ь2
уh
аЬ = — = а, Ас = с - с < 0,
В н п
(7)
где х - аналог числа Фурье (диффузионный), а - аналог числа Био(диффузионный).
В таком случае задача (3)-(6) примет форму
8и 8 2и
8х 8С2 ' и(С, 0) = 0,
8(0, г)
= 0,
8и(1, г)
8С
аи (1, г) = а.
(8) (9) (10)
(11)
роН@уЫп1к-угиИ. ги
Прибегая к так называемому методу осреднения, при исследовании задачи с учетом наличия фронта концентрации, решают ее по двум фазам, заменяя левую часть уравнения (8) ее осредненным значением по интервалу (р, 1), а именно, вводят в рассмотрение функцию [7]
2ф(Р)(х) = -±- Г^С, (12)
1 - рJ 8х
где р - индекс фазы (Р = 1,2);р - безразмерная координата фронта концентрации.
Тогда вместо уравнения (8) используют
приближенное соотношение
82и(1)
8иг = 2ф(1)(х), (13)
удовлетворяющее, согласно (8)-( 11) условиям ы0)(р,х) = 0, (14)
[8м(1)(2,,х) /8^ = 0, (15)
[8м(1)(^,х) /84 + аи(1)(^,х)]?=1 =а, (16) и, кроме того, вытекающему из (12) равенству
1 г 8и °
2ф(1)(х) = Г8-—
1 - рГ 8х
(17)
Проводя последовательно интегрирование в (13), имеют
8и(1)(4,х) /84 = 2ф(1)(х)4 +АЬ (18) м(1)(4,х) =ф(1)(х)42 +А:4 +ВХ. (19) Согласуя (15)-( 17) с условиями (14)-(16), для определения неизвестных ф(1)(х), А\, В\, с использованием (18), (19), получают систему алгебраических уравнений
(а + 2)ф(1) + (а + 1) А1 + аВ1 =а, (20) р2ф(1) + рА1 + В1 = 0,2рф(1) + А1 = 0, (21) решением (20) является
ф(1)(х) =а/Вь А1 =-2ар/А, В1 =ар2/Вь (22) где обозначено
В1 =а(1 - р)2 + 2 (1- р).
(23)
Подставляя (21) в (19), находят
и(1)(4,х) =а(4 - р)2/В1.
Функция р =р(х) определяется с помощью (17)
Поскольку, согласно (23)
8и
а йр
В
= --т-ЯВС-р) + (С - р)2]
8х В йх
йр
то после подстановки в уравнение (17) значения ф(1)(х) и значения 8и(1)/8х, получают
Бетепоу Е.М et аСФгосеебтр о/МЗЦЕТ, 2021, 83, по. 1, рр.
А(1 ~Р) =
= - ^[2ДГ'(<-р) С< + ^ 1\<-р)2 С<], ат ар
откуда, после интегрирования
6 £1(1 - р)= -(1 - р)2 ^[Щ + (1 - р).
dг
Поэтому с учетом (22) имеют СЮ1
Ср
С р
= - 2[а(1 -р) + 1],
сЮ,
3А + (1 - р)-^ = а(1 - р)2 + 4(1 - р). а р
В результаты чего приходят к дифференциальному уравнению относительно р =р(г)
6[а(1 - р) + 2] = - (1 - р)Ср[ет(1 - р) + 4],
Сг
или, после разделения переменных к соотношению
л _ _ а -р)[а(1 -р) + 4] ср. (24) 6[а(1 -р) + 2]
Поскольку выполняется начальное условие г = 0 при р = 1, то интегрируя это выражение слева по г от нуля до г, справа - по р от р до единицы, исходя из (24), т. е., выполняя процедуру
1(1 -р)[а(1 -р) + 4]
I Сг = -!(
6[а(1 -р) + 2]
-Ср,
получают частное решение
г = £-рр1 + -Лшí 1 + 0!^) I. (25)
12
3а 3а2
Продолжительность первой фазы находят путем подстановки р =0 в (25)
1 1
Г1 =— +
2 а + 2
-1п(—:—).
2 -V (26)
12 3а 3а2 2 Если (численным путем) разрешить (25) относительно р как неявную функцию, то тем самым определяют и фронт концентрации
р =р(г).
Вторая фаза. Функцию ф(2)(г) вводят по соотношению
1 ды ()
2ф(2)(г) =
(27)
Задача ставится таким образом, чтобы получить решение уравнения
д 2ы (2)
д-Ыг = 2ф(2)(г), (28)
удовлетворяющее условиям:
и(2)(0,п) = 0, (29)
ди(2)(0,г) /д£ = 0, (30)
ди(2)( 1,г) /д£ + аи(2)( 1,г) =а, (31)
62-70 [email protected]
Интегрируя, последовательно, по £ уравнение (28), имеют:
ди(2)(£,г) /д£ =2ф(2)(х)£+Л2(х), (32) и(2)(£,г) =ф(2)(г)£2 +^2. (33)
Поскольку из выражений (30), (32) вытекает, что ^2(г) = 0, то в таком случае согласно (33) получают
и(2)(£,г) =ф(2)(т)£2+£2. (34)
Удовлетворяя условию (31), имеют 2ф(2)(г) + а[ф(2)(г) + В2] =а, откуда находят
В2 = 1 - (а + 2)ф(2)(г)/а.
В таком случае (34) принимает вид и(2)(£,г) = 1 - ф(2)(г) [(а + 2) /а - £2]. (35)
Функция ф(2)(г) определяется на основе условий (27), (29), с учетом того, что, с одной стороны,ди(2)(1,г) /дг = -дф(2)(г)/дг-[(а + 2)/а -£2], с другой стороны
1(2)
2а + 6
ф) (,).. 2 с (£±+ _<)<- 1---
2 Сг „ а 6 Сг а
и поэтому, разделяя далее переменные в этом уравнении и интегрируя полученное в результате соотношение, находят общее решение ф(2)(г) =Ю2ехр [-3аг/(а +3)], где Ю2 - произвольная постоянная.
Тогда в соответствии с (35)
(<, ' )= 1- Ю2
а + 2
а
- < \ехр
-3s(t - /1)
(36)
где Г1определяют по (26).
Далее, принимая во внимание (29), находят
Ю2 =а/(а + 2) ехр [3аг1 /(а +3)].
После чего, подставляя Ю2 в (36), окончательно получают
Г-3я(/-/1)
и(2) <) = 1-
(* + 2)/
а + 2 2 | --< | х ехр
(* + 3)
. (37)
Таким образом, в рамках принятой физико-математической модели задача об определении поля концентрации сахарозы в области между двумя смежными кристаллами сахарозы в вакуум-аппарате решена полностью.
В дальнейшем, возвращаясь согласно (7) к исходному значению с концентрации сахарозы в сахарсодержащем растворе, имеют
с (£,г) =Сп+ и(2)(£,г)Дс, (38)
где и(2)(£,г) определяется по соотношению (37).
Для того чтобы рассчитать эффективность процесса кристаллизации сахарозы, предварительно, на базе зависимости (2), (37) определяется поток q концентрации.
0
р
Семенов Е.В. и др.Вестник.ФТУИШ, 2021, Т. 83, №. 1, С- 62-70
При этом, поскольку в соответствии с (37), (38) имеют
дс ди (2\ 2аЕ-Дс
— =-Дс =-ехр
д£ д£ а + 2
-3s (t - tj)
(39)
, м-1, (40) то в соответствии с (2), (40)
(^ + 3)
или, согласно (7), в размерных переменных дс _ дс дх hд% справедлива формула q = D дс/ дх.
Поэтому поток q концентрации на поверхности х =^4 =1), в явной форме выражается в виде
2аДcD
q =-ехр
(а + 2)h
-3s (t - t1)
. (s +3) .
В свою очередь, так как за время dx на той же поверхности осаждается объем сахарозы
dQ =
qh2dt 2аДск
D (а + 2)
ехр
-3а (т -т1)
. (а +3) .
dt
(4J)
то количество АУ сахарозы, оседающей за период времени т получают, интегрируя выражение (41) в пределах от 0 до т. В результате имеют
Q (1, А, И) = АУ
2аДск
Т ехр
(а + 2) о 2(а + 3)НДс
-3а (т -т1)
(а +3) . (а + 3)
dT =
ехр
exp
(а + 3).
3(а + 2) или по массе
Q(t, D, h) =Q(t, D, h)xpj, кг/м2 где pj - плотность сахарозы.
,м3/м2
(42)
роНФуЫпк-угиИ. ги Численный эксперимент
Пусть в одном кубометре утфеля содержится, в равных объемных долях кристаллы сахарозы размером 5 =10-4 м и межкристального раствора коэффициентами, соответственно, пересыщения Сп = 0,8 и насыщения сн = 0,7. И пусть, по соглашению, коэффициент кри-сталлохимической реакции у = 10-6 м2/с [4, 5], коэффициент диффузии раствора De [10-10, 5х10_п], м2/с [2]. Требуется на базе выбранных геометрических и физических параметров процесса определить по массе отнесенный к 1 м3утфеля полный молекулярный расход межкристального раствора в зависимости от времени проведения процесса кристаллизации.
Предварительно, поскольку, согласно формуле (1) при ^ = 0,5 половина расстояния между частицами сахарозы составляет h = 6,3 х10-5 м, то в расчетах принимали he[5х10"5; 10-4], м. В свою очередь, так как количество частиц сахарозы равно ^х6/(л53), то их общая поверхность составляет 5 = 3/5 = 3х104 1/м3.
С учетом (38) на рисунок 2 представлены графики зависимости концентрации сахарозы на середине расстояния между двумя смежными частицами (т. е. при х = 0)
Как следует из анализа поведения кривых (рисунок 2), к моменту т = 5 концентрация сахарозы посередине между двумя кристаллами в результате конденсации молекул сахарозы на кристаллах примерно выравнивается, приближаясь к значению сн = 1,1 насыщения.
с
0.S5
i i i i
clCO
c2(r) 03 \
c3(x) i i T
4 it " . ■ 3 ,
C.7 i i ——
г 5
Рисунок 2. Зависимости концентрации сахарозы между двумя смежными частицами от безразмерного времени т процесса обработки центров кристаллизации в сахарсодержащем растворе при различных значениях коэффициента D диффузии и половины расстояния h между кристаллами (h = 5x10-5 м: 1 - D = 5x10-11, 2 - D =10-10 м2/с; h = 10-4 м: 3 - D = 5x10"", 4 - D = 10-10 м2/с)
Figure 2. Dependences of the concentration of saccharose between two adjacent particles on the dimensionless time treatments of the processing of crystallization centers in a sugar-containing solution at different values of the diffusion coefficient D and half the distance h between crystals (h = 5x 10-5 m:1 - D = 5x 10-11, 2 - D = 5x 10-10 m2/s; h =10-4 m: 3 - D = 5x 10-11, 4 - D = 5x 10-10 m2/s)
Semenov E..V. et aC.Proceedings of VSUET, 2021, voC. 83, no. 1, pp. 62-70
тжт
post@vestnik-vsuet. ru
500 М1(т)лоо-
МВД .
I-г
3 4
~i-1—
. - - - - -.M '
\bZ
3M1(t>
змздш
ЗЩ) 120 100
%
3 -■*' 4
1 1
2
)* Т II 1 1
т,кг -
1 1 1
400 - -
300 М(т) — :оо - 1 -
:оо -
0 1 1 I
0 : 10 1Г г,мин
С D 60
Рисунок 3. Зависимости от временит массы m конденсированной на поверхности кристалла сахарозы при различных значениях коэффициента D диффузиии половины расстояния h между кристаллами (h = 5х10-5 м: 1 - D = 5х10-11, 2 - D = 10-10 м2/с; h = 10-4 м: 3 - D = 5х10-11, 4 - D = 10-10 м2/с) Figure 3. Time т dependences of the mass m of sucrose condensed on the surface of the crystal at different values of the coefficient D diffusion and half distance h between crystals (h = 5х10-5 m: 1 - D = 5х10-11 m2/s, 2 - D = 10-10 /с; h = 10-4 м: 3 - D = 5х10-11, 4 - D = 10-10 m2/s)
Рисунок 4. Зависимости от времени т относительной массы 5m/m0, % конденсированной на поверхности кристалла сахарозы при различных значениях коэффициента D диффузии и половины расстояния h между кристаллами (h = 5х10-5 м: 1 - D = 5х10-11, 2-D= Ю"10 м2/с; h = 10"4 м: 3 - D = 5х 10"11, 4-D= 10"10 м2/с) Figure 4. Time т dependences of the relative mass 6m niu. % of sucrose condensed on the surface of the crystal for different values of the diffusion coefficient D and half the distance h between the crystals (h = 5х 10-5 m: 1 - D = 5х10-11, 2 - D = 10-10 h = 10-4 м: 3 - D = 5х10-11, 4 - D = 10-10 m2/s)
Рисунок 5. Зависимость от времени массы осадившейся на поверхности взвеси из кристаллов сахарозы Figure 5. Time dependence of the mass of a suspension of sucrose crystals deposited on the surface
Отраженные графиками рисунок 3, 4 результаты вычислений на базе (42) обнаруживают согласованность полученных результатов количественного моделирования исследуемого процесса с физическим смыслом задачи.
А именно, с одной стороны, обессахарива-ние раствора по массе в связи с кристаллизацией сахарозы, как по абсолютному (рисунок 3), так и относительному к исходному (рисунок 4) значениям, по прошествии времени экспоненциально возрастает, и тем больше, чем выше значение коэффициента диффузии (например, на рисунок 3, кривая 2 расположена выше кривой 1).
С другой стороны, аналогичным образом, графически, например, рисунком 4 отражена зависимость прироста массы от времени: при коэффициенте диффузии D = 5 х 10-11 м2/с, с увеличением расстояния h (и, следовательно, с уменьшением счетной концентрация частиц в растворе), интенсивность процесса кристаллизации снижается - кривая 3 ниже кривой 1.
Если, при тех же исходных данных по параметрам процесса расчет проводится в реальном временном режиме, то в таком случае, например, при значениях параметров h = 10-4 м и D = 5х10-11 м2/с графическая зависимость от времени t массы т осадившейся на поверхности взвеси из кристаллов сахарозы в вакуум-аппарате имеет показанный на рисунок 5 вид.
Как видно по рисунок 5, процесс кристаллизации завершается практически к пятнадцатой минуте и для возобновления его необходимо раствор "раскачивать", обогащая его пресыщенной сахарсодержащей субстанцией.
Однако следует отметить, что отраженные рисунками 2-5 результаты расчетов, из-за того в них не учитывается эффект истощения раствора сахарозой (в частности, коэффициент Д = Сн - Сп полагается постоянной величиной), приводят к завышенным, по сравнению с наблюдаемыми на практике, данными.
т
т
Семенов Е.В. и др.Вестник,ФГУИТ, 2021, Т. 83, №.. 1, С- 62-70 [email protected]
Заключение кристаллизации сахарозы в данном растворе при
На основе модели диффузионного массо- учете кинетики #онта концентрации. переноса целевого продукта, такого как сахароза, Предлагается даал™^™ аппарат для
из пересыщенного межкристального раствора расчета зависимости массы целевого продукта
к кристаллу сахарозы, основанного на уравнении от времени в результате процесса крисгаллизгщда. нестационарной диффузии (второго закона На основе разработанного аналигачеаюго
Фика), обосновывается физико-математическая аппарата проводится количественное модели-
модель по прогнозированию протекания процесса рование эффективности данного процесса.
Литература
1 Силин П. М. Технология сахара. Рипол Классик, 2013.
2 Магомедов М. Г. Технология получения пасты из сахарной свеклы //Вестник ВГУИТ. 2015. № 3. С. 138-141. doi: 10.20914/2310-1202-2014-3-138-141
3 Славянский А.А., Семенов Е., Лебедева Н., Неделькин В. и др. Математическое моделирование кинетики кристаллизации гидратной глюкозы // Сахар. 2014. № 5. С. 45-49.
4 Хворова JI.C. Трехпродуктовая технологическая схема получения глюкозы с кристаллизацией двух продуктов в ангидридной форме // Пищевая промышленность. 2017. № 9.
5 Хворова JI.C., Баранова JI.B., ГоменюкВ.А. Способ получения кристаллической ангидридной глюкозы. 2018.
6 Семенов Е.В., Славянский А.А. Моделирование процесса кристаллизации и центрифугирования. М.: Спутник+, 2015. 217 с.
7 Кулинченко В.Р., Мирончук В.Г. Промышленная кристаллизация сахаристых веществ. 2012.
8 Berzi К., Suryanarayanan R. Compression-induced crystallization in sucrose-polyvinylpyrrolidone amorphous solid dispersions//Crystal Growth & Design. 2018. V. 18. №2. P. 839-848. doi: 10.1021/acs.cgd.7b01305
9 Hubbes S.S., Braun A., Foerst P. Sugar particles and their role in crystallization kinetics and structural properties in fats used for nougat creme production // Journal of Food Engineering. 2020. V. 287. P. 110130. doi: 10.1016/j.jfoodeng.2020.110130
10 Velazquez-Camilo O., Bolanos-Reynoso E., Rodriguez E., Alvarez-Ramirez J. Characterization of cane sugar crystallization using image fractal analysis // Journal of Food Engineering. 2010. P. 77-84. doi: 10.1016/j.jfoodeng.2010.03.030
11 Khaddour I., Bento L., Ferreira A., Rocha F. Kinetics and thermodynamics of sucrose crystallization from pure solution at different initial supersaturations // Surface Science. 2010. P. 1208-1214. doi: 10.1016/j.susc"2010.04.005
12Frenzel S. Crystallization Schemes in the Sugar Industry// ChemBioEng Rev. 2020. №5. P. 1-9. doi: 10.1002/cben.202000010
13 Thakral S., Sonje J., SuryanarayananR. Anomalous behavior of mannitol hemihydrate: Implications on sucrose crystallization in colyophilized systems // International journal of pharmaceutics.2020. V. 587. P. 119629. doi: 10.1016/j.ijpharm.2020.119629
14 de Castro B. J. C. et al. Sucrose crystallization: modeling and evaluation of production responses to typical process fluctuations //Brazilian Journal of Chemical Engineering" 2019. V. 36. № 3. P. 1237-1253. doi: 10.1590/0104-6632.20190363s20180240
15 SinghK, Gupta S.P., Kumar A., Kumar A. The effect of high intensity ultrasound (HIU) on the kinetics of ciystallizationof sucrose: Elimination of latent period//Ultrasonics sonochemistiy. 2019. V. 52. P. 19-24. doi: 10.1016/j.ultsonch.2018.05.030
16Chakir A., Romane A., Marcazzan G.L., Ferrazzi P. Physicochemical properties of some honeys produced from different plants in Morocco // Arabian Journal of Chemistry. 2016. V. 9. P. S946-S954. doi: 10.1016/j.arabjc.2011.10.013
17Nayik G.A., Dar B.N., Nanda V. Physico-chemical, rheological and sugar profile of different unifloial honeys from Kashmir valley of India// Arabian Journal of Chemistry. 2019. V. 12. № 8. P. 3151-3162. doi: 10.1016/j.arabjc.2015.08.017
1 8E1 Batal H., Hasib A., Ouatmane A., Dehbi F. et al. Sugar composition and yield of syrup production from the pulp of Moroccan carob pods (Ceratonia siliqua L.) // Arabian Journal of Chemistry. 2016. V. 9. P. S955-S959. doi: 10.1016/j.arabjc.2011.10.012
19Boussaid A., Chouaibi M., Rezig L., Hellal R. et al. Physicochemical and bioactive properties of six honey samples from various floral origins from Tunisia // Arabian journal of chemistry. 2018. V. 11. № 2. P. 265-274. doi: 10.1016/j.arabjc.2014.08.011
20Torr6nen R., Hellstrom J., Mattila P., Kilpi K. Postprandial glycaemic response to berry nectars containing inverted sucrose // Journal of nutritional science. 2017. V. 6. doi: 10.1017/jns.2016.44
References
1 Silin P.M. Technology of sugar. Ripol Classic, 2013. (in Russian).
2 Magomedov M.G. Technology of obtaining paste from sugar beet. Proceedings of VSUET. 2015. no. 3. pp. 138-141. doi: 10.20914/2310-1202-2014-3-138-141 (in Russian).
3 Slavyansky AA, Semenov E., Lebedeva N., Nedelkin V. et al. Mathematical modeling of the kinetics of crystallization of hydrated glucose. Sugar. 2014. no. 5. pp. 45-49. (in Russian).
4 Khvorova L.S. Three-product technological scheme for producing glucose with crystallization of two products in anhydride form. Food Industry. 2017. no. 9. (in Russian).
5 Khvorova L.S., Baranova L.V., Gomenyuk V.A. A method of obtaining crystalline anhydride glucose. 2018. (in Russian).
6 Semenov E.V., Slavyansky A.A. Simulation of the crystallization and centrifugation process. Moscow, Sputnik +, 2015. 217 p. (in Russian).
7 Kulinchenko V.R., Mironchuk V.G. Industrial crystallization of sugary substances. 2012. (in Russian).
8 Berzi K., Suryanarayanan R. Compression-induced crystallization in sucrose-polyvinylpyrrolidone amorphous solid dispersions. Crystal Growth & Design. 2018. vol. 18. no. 2. pp. 839-848. doi: 10.1021/acs.cgd.7b01305
9 Hubbes S.S., Braun A., Foerst P. Sugar particles and their role in crystallization kinetics and structural properties in fats used for nougat creme production. Journal of Food Engineering. 2020. vol. 287. pp. 110130. doi: 10.1016/j.jfoodeng.2020.110130
Semenov E..V. et aC.Proceedings of VSUET, 2021, voC. 83, no. 1, pp. 62-70
postevesinik-vsuet. ru
10 Velazquez-Camilo O., Bolanos-Reynoso E., Rodriguez E., Alvarez-Ramirez J. Characterization of cane sugar crystallization using image fractal analysis. Journal of Food Engineering. 2010. pp. 77-84. doi: 10.1016/j.jfoodeng.2010.03.030
11 Khaddour I., Bento L., Ferreira A., Rocha F. Kinetics and thermodynamics of sucrose crystallization from pure solution at different initial supersaturations. Surface Science. 2010. pp. 1208-1214. doi: 10.1016/j.susc.2010.04.005
12 Frenzel S. Crystallization Schemes in the Sugar Industry. ChemBioEng Rev. 2020. no. 5. pp. 1-9. doi: 10.1002/cben.202000010
13 Thakral S., Sonje J., Suryanarayanan R. Anomalous behavior of mannitol hemihydrate: Implications on sucrose crystallization in colyophilized systems. International journal of pharmaceutics. 2020. vol. 587. pp. 119629. doi: 10.1016/j.ij pharm.2020.119629
14 de Castro B.J.C. et al. Sucrose crystallization: modeling and evaluation of production responses to typical process fluctuations. Brazilian Journal of Chemical Engineering. 2019. vol. 36. no. 3. pp. 1237-1253. doi: 10.1590/0104-6632.20190363s20180240
15 Singh K., Gupta S.P., Kumar A., Kumar A. The effect of high intensity ultrasound (HIU) on the kinetics of crystallization of sucrose: Elimination of latent period. Ultrasonics sonochemistry. 2019. VOL. 52. PP. 19-24. doi: 10.1016/j.ultsonch.2018.05.030
16 Chakir A., Romane A., Marcazzan G.L., Ferrazzi P. Physicochemical properties of some honeys produced from different plants in Morocco. Arabian Journal of Chemistry. 2016. vol. 9. pp. S946-S954. doi: 10.1016/j.arabjc.2011.10.013
17 Nayik G.A., Dar B.N., Nanda V. Physico-chemical, rheological and sugar profile of different unifloral honeys from Kashmir valley of India. Arabian Journal of Chemistry. 2019. vol. 12. no. 8. pp. 3151-3162. doi: 10.1016/j.arabjc.2015.08.017
18 El Batal H., Hasib A., Ouatmane A., Dehbi F. et al. Sugar composition and yield of syrup production from the pulp of Moroccan carob pods (Ceratonia siliqua L.). Arabian Journal of Chemistry. 2016. vol. 9. pp. S955-S959. doi: 10.1016/j.arabjc.2011.10.012
19 Boussaid A., Chouaibi M., Rezig L., Hellal R. et al. Physicochemical and bioactive properties of six honey samples from various floral origins from Tunisia. Arabian journal of chemistry. 2018. vol. 11. no. 2. pp. 265-274. doi: 10.1016/j.arabjc.2014.08.011
20 Torronen R., Hellstrom J., Mattila P., Kilpi K. Postprandial glycaemic response to berry nectars containing inverted sucrose. Journal of nutritional science. 2017. vol. 6. doi: 10.1017/jns.2016.44
Сведения об авторах Евгений В. Семенов д.т.н., профессор, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Талалихина, 31, г. Москва, 109004, Россия, зет-ровШтаП.ги
https://orcid.org/0000-0001-6203-1783 Анатолий А. Славянский д.т.н., профессор, зав. кафедры технология продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Талалихина, 31, г. Москва, 109004, Россия, апа1оЬу4455(й!уа.ги
https://orcid.org/0000-0002-0262-8841 Вера А. Грибкова к.т.н., доцент, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Талалихина, 31, г. Москва, 109004, Россия, твййлйррйуа.ги
https://orcid.org/0000-0002-5079-8132 Дарья П. Митрошина ассистент, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Талалихина, 31, г. Москва, 109004, Россия, d_mitr96(й)mail.ru
https://orcid.org/0000-0002-8724-3368 Сергей Т. Антипов д.т.н., профессор, зав. кафедры машины и аппараты пищевых производств, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия, авШувие!™ https://orcid.org/0000-0003-4299-1538
Вклад авторов
Все авторы в равной степени принимали участие в написании рукописи и несут ответственность за плагиат
Конфликт интересов
Авторы заявляют об отсутствии конфликта интересов.
Information about authors Evgeny V. Semenov Dr. Sci. (Engin.), professor, technology of products from plant raw materials and perfumery and cosmetic products department, Moscow State University of Technology and Management. K.G. Razumovsky (PKU), st. Talalikhina, 31, Moscow, 109004, Russia, sem-post(S!mail.ru
https://orcid.org/0000-0001-6203-1783 Anatoly A. Slavyansky Dr. Sci. (Engin.), professor, head of department, technology of products from plant raw materials and perfumery and cosmetic products department, Moscow State University of Technology and Management. K.G. Razumovsky (PKU), st. Talalikhina, 31, Moscow, 109004, Russia, anatoliy4455(S)ya.ru https://orcid.Org/0000-0002-0262-8841
Vera A. Gribkova Cand. Sci. (Engin.), associate professor, technology of products from plant raw materials and perfumery and cosmetic products department, Moscow State University of Technology and Management. K.G. Razumovsky (PKU), st. Talalikhina, 31, Moscow, 109004, Russia, instituttpp(S)ya.ru https://orcid.org/0000-0002-5079-8132
Daria P. Мitroshina assistant, technology of products from plant raw materials and perfumery and cosmetic products department, Moscow State University of Technology and Management. K.G. Razumovsky (PKU), st. Talalikhina, 31, Moscow, 109004, Russia, d_mitr96(S!mail.ru https://orcid.org/0000-0002-8724-3368
Sergey T. Antipov Dr. Sci. (Engin.), professor, head of department machines and apparatus for food production, Voronezh State University of Engineering Technologies, Revolyutsii Ave., 19, Voronezh, 394036, Russia, ast(S)vsuet.ru https://orcid.org/0000-0003-4299-1538 Contribution
All authors are equally involved in the writing of the manuscript and are responsible for plagiarism
Conflict of interest
The authors declare no conflict of interest.
Поступила 11/01/2021_После редакции 22/02/2021_Принята в печать 01/03/2021
Received 11/01/2021 Accepted in revised 22/02/2021 Accepted 01/03/2021
