CC BY
21
Биотехнология, бионанотехнология и технология сахаристых продуктов
Оригинальная статья/Original article_
УДК 66.061 ~ DOI: http://doi.org/10.20914/2310-1202-2016-4-233-237_
Расчет периода раскачки утфеля раствором сахарозы _в вакуум-аппарате_
Анатолий А. Славянский 1 anatoliy4455@yandex.ru Евгений В. Семенов 1 sem-post@mail.ru Наталья Н. Лебедева 1 mgutu-sahar@mail.ru _Сергей Т. Антипов_2 ast@vsuet.ru_
1 кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Москов. гос. ун-т технол. и управл. им. К.Г. Разумовского (ПКУ), ул. Земляной вал, 73, г. Москва, 109004, Россия
2 кафедра машин и аппаратов пищевых производств, Воронеж. гос. ун-т. инж. техн., пр-т Революции, 19, г. Воронеж, 394066, Россия Реферат. Важным и ответственным элементом в управлении технологическим потоком производства кристаллического сахара товарных кондиций является период времени обогащения утфеля в виду его обессахаривания в процессе кристаллообразования в вакуум-аппарате. Хотя в настоящее время предложен и внедрен в промышленность, в том числе и в качестве интерфейсного сопровождения процесса, ряд способов, в реальных условиях сахарного производства во многих случаях приходится прибегать к услугам опытного варщика. Очевидно, что в любом случае необходимо располагать зависящими от объемного содержания сахара в растворе данными по времени подачи сахарсодержащего раствора в вакуум-аппарат. Касательно вопроса о периоде протекания процесса обогащения обедненного сахарозой раствора интродуцируемым в эту субстанцию раствором сахарозы следует отметить, что данная проблема теоретически проработана все еще недостаточно. При этом очевидно, что в практических целях желательно иметь простые и удобные для инженерного расчета зависимости времени обработки сахаром заданного объема воды от режимных параметров процесса (требуемой концентрации раствора по сахарозе, температуре раствора, условий перемешивания и др.). Проблема количественного анализа процесса кристаллизации сахарозы в вакуум-аппарате, включая сроки обогащения межкристального раствора сиропом, периода времени обработки утфеля в целом в данном аппарате, исследовалась во многих работах. Однако вследствие своей важности для решения задачи получения товарного сахара высоких кондиций эта проблема нуждается в дальнейшем углубленном изучении. В статье с целью обоснования процесса обогащения раствора сахарозой в результате обессахаривания этой жидкостной системы в вакуум-аппарате, с позиций диффузионной теории кристаллообразования в растворе, предлагается более обоснованный, по сравнению с известными до сих пор аналогами, количественный анализ данного процесса, где в качестве коллектива кристаллов сахарозы рассматривают систему шаров, равномерно распределенных в моделирующем вакуум-аппарат объеме. На базе решения поставленной задачи получены расчетные формулы, результаты численного эксперимента отражены графическим путем.
^лючевые^лова^ваку^м-^ппиЕии^^иииииЕиихиЕизы^^^ийеии
Calculation of period processing solution syrup in vacuum apparatus
Anatolii A. Slavyanskii 1 anatoliy4455@yandex.ru Evgenii V. Semenov 1 sem-post@mail.ru Nataliya N. Lebedeva 1 mgutu-sahar@mail.ru _Sergei T. Antipov_2 ast@vsuet.ru_
j technologies of products from vegetable raw materials and prafyumery-cosmetic products department, K.G. Razumovsky Moscow State
University of technologies and management (the First Cossack University), Zemlyanoi val st., 73, Moscow, 109004, Russia 2 machinery and equipment for food production department, Voronezh state university of engineering technologies, Revolution Av., 19, Voronezh,
394066, Russia_
Summary. Important and crucial element in the management of the technological flow of production of sugar product standards is the period of time the enrichment of massecuite, since its neutralization in the process of crystal formation in vacuum apparatus, excess sugar solution. Although currently proposed and implemented in the industry, including as a front-end accompany the process, a number of ways in the real world sugar production in many cases have to resort to the services of an experienced operator. It is obvious that in any case it is necessary to have a surround-dependent glucose solution data on time for the excess sugar solution into the vacuum apparatus. With regard to the period of the enrichment of depleted sucrose solution are entered into this substance excess sucrose solution, it should be noted that this problem is theoretically still insufficiently developed. It is obvious that for practical purposes it is desirable to have a simple and convenient for engineering calculation of sugar processing time dependencies of the specified volume of water from the operating parameters of the process (the required concentration of sucrose, temperature of the solution stirring). The problem is the quantitative analysis of sucrose crystallization in vacuum apparatus, including the timing of enrichment solution to the excess syrup, period of time processing massecuite total this apparatus has been investigated in many works. However, due to its importance to the task of obtaining commercial sugar high standards this issue required further in-depth examination. In the article to support the enrichment process solution sucrose due to neutralize this solvent system in vacuum apparatus, from the standpoint of diffusion theory provides a more reasonable compared to known so far, quantitative analysis of this process. Where as sucrose crystals team are considering a system of balls, uniformly distributed in vacuum apparatus. On the basis of the solution obtained by calculation formulas, results of numerical experiment reflects the graphic way. Keywords: vacuum apparatus, solution of sucrose, diffusion process, processing period
Для цитирования Славянский А.А., Семенов Е.В., Лебедева Н.Н., Антипов С.Т. Расчет периода раскачки утфеля раствором сахарозы в вакуум-аппарате // Вестник ВГУИТ. 2016. № 4. С. 233-237. doi:10.20914/2310-1202-2016-4-233-237
For citation
Slavyanskiy A.A., Semenov E.V., Lebedeva N.N., Antipov S.T. Calculation of period processing solution syrup in vacuum apparatus. Vestnik VSUET [Proceedings of VSUET]. 2016. no. 4. pp. 233-237. (in Russian). doi:10.20914/2310-1202-2016-4-233-237
BecmHUKßry^T/Proceeätngs of VSUET, № 4, 2016L
Введение
Одним из важных и ответственных элементов в управлении технологическим потоком производства кристаллического сахара товарных кондиций является период времени обогащения утфеля, в виду его обессахаривания в процессе кристаллообразования в вакуум-аппарате. Несмотря на то, что в настоящее время предложен и внедрен в промышленность в том числе, и в качестве интерфейсного сопровождения процесса, ряд способов, в реальных условиях сахарного производства во многих случаях приходится прибегать к услугам опытного варщика. Очевидно, что в любом случае необходимо располагать зависящими от объемного содержания сахара в растворе данными по времени подачи сахарсодержа-щего раствора в вакуум-аппарат.
Касательно вопроса о периоде протекания процесса обогащения обедненного сахарозой раствора интродуцируемым в эту субстанцию раствором сахарозы следует отметить, что данная проблема теоретически проработана все еще недостаточно. При этом очевидно, что в практических целях желательно иметь простые и удобные для инженерного расчета зависимости времени обработки сахаром заданного объема воды от режимных параметров процесса (требуемой концентрации раствора по сахарозе, температуре раствора, условий перемешивания и др.).
Проблема количественного анализа процесса кристаллизации сахарозы в вакуум-аппарате, включая сроки обогащения межкристального раствора сиропом, периода времени обработки утфеля в целом в данном аппарате, исследовалась в работах [1-6]. Однако вследствие своей важности для решения задачи получения товарного сахара высоких кондиций эта проблема, по нашему мнению, нуждается в дальнейшем углубленном изучении.
Ниже, с целью обоснования процесса обогащения раствора сахарозой в результате обессахаривания этой жидкостной системы в вакуум-аппарате, с позиций диффузионной теории кристаллообразования в растворе, предлагается более обоснованный, по сравнению с известными до сих пор аналогами, количественный анализ данного процесса.
Постановка и решение задачи
Пусть в некотором рабочем объеме имеют гетерогенную жидкостную систему, включающую жидкую среду и коллектив моделируемых как шары радиусом Я кристаллов сахарозы. Предполагая, что шары статистически равномерно распределены в жидкости, в качестве модели процесса массопереноса от жидкой
среды к шарам выбирают изолированный шар, к которому из прилегающей к нему области (Я < г < Ь) (где г - радиальная координата, 21 - среднее расстояние между центрами шаров) подаются молекулы сахарозы (рисунок 1).
Кристалл сахарозы Sucrose crystal
Сахарсодержащий раствор (утфель) Solution of sucrose
Поток концентрации Flow concentration
Рисунок 1. Схема к расчету процесса кристаллообразования сахарозы в утфеле
Figure 1. The scheme to calculation of process of grain formation of sucrose in the massecuite
При этом, если Уш и V- соответственно, суммарный объем шаров и жидкости в рабочем объеме и £ - объемная концентрация шаров в жидкостной системе, то в соответствии с принятой геометрической моделью жидкостной системы "жидкость + система шаров" половина расстояния l между двумя соседними шарами составит (рисунок 2)
I = RÇ-1'3, 0 <£< 1.
(1)
* У
А
Сахарсодержащий раствор
l
Рисунок 2. Схема к расчету процесса массопереноса сахарозы в жидкостной системе
Figure 2. The scheme to calculation process of a mass transfer of sucrose in liquid system
Полагают, что массоперенос к шарам от жидкой среды реализуется симметричным образом, и, кроме того, для выбранной расчетной модели, с удовлетворительной точностью выполнены все, обычно принимаемые при анализе процесса массопереноса в сплошной среде, допущения по физико-механическим свойствам изучаемого объекта.
ФестникВГУИТ/Огоаве&чр of ЪЗ&ЕТ, № 4, 2016,
Поток концентрации q сахарозы в направлении оси г в принятой одномерной модели массопереноса в сферической системе отсчета рассчитывают по формуле
дс
Я = ~В —, (2)
дг
где В - коэффициент диффузии, с - объемная концентрация сахарозы, г - радиальная координата.
В качестве исходного соотношения, описывающего кинетику распределения концентрации в растворе, выбирают, как обычно, отнесенное к сферическим координатам, с началом координат в центре О1 шара, уравнение нестационарной диффузии (второй закон Фика) (рисунок 2).
^ = В
дл
Г Я 2
д2с дг2
2 дс г дг
\
(Я<г<Ь ),
(3)
V /
где л - время.
Начальное условие по исходному распределению концентрации принимают в виде
с(г, 0) = сп (Я < г < I), (4)
где сп - концентрация пересыщенного раствора. Граничное условие на поверхности шара
с(Я, г )= с (0< г«ю), (5)
где с н - концентрация сахарозы на поверхности шара (кристалла сахарозы).
Граничное условие симметричности поля концентрации между соседними шарами
дс(Ц) дг
= 0 (0<л<х),
(6)
где I связана с Я зависимостью (1).
Хотя решение задачи (2)-(6) известно, оно малопригодно для практических расчетов, так как имеет вид ряда. Поэтому, как поступают в некоторых случаях, ниже применяют приближенный метод решения данной задачи, аналогичный методу осреднения.
С этой целью, как обычно, предварительно проводят замену переменной
с = и + сн, (7)
в результате чего граничные условия (5), (6) становятся однородными по переменной и - приведенной концентрации
и(Я, 0 = 0, (8)
ди(1, Л)
дг
= 0, 0<?<о>.
(9)
Начальное условие вследствие (4) принимает форму
и(г,0)=Ас(Я< г< I), (10)
где Ас = сп - сн > 0 - перепад концентрации в начальный момент времени ^ = 0.
При этом уравнение (3) по переменной и вследствие (7) сохраняет свой вид
^ = В дл
дг2
2 ди
г дг
(11)
В целях упрощения решения задачи начальное условие (10) заменяют его усредненным значением по прилегающей к шару области г е (Я, 1)
1
13 - Я:
-|и(г,0)г Л = Ас. (12)
В свою очередь, усредняя левую часть уравнения (11) по тому же объему, принимают
1 .Г ^ г Л = Бф(0,
13 - Я3 [ дЛ
(13)
где ф(т) - функция, подлежащая определению.
В результате уравнение (11) заменяют приближенным
^д2и 2 ди^
Vдг2 г дг у
= ф(л).
(14)
Так как правая часть уравнения (14) зависит лишь от 1 то, интегрируя его дважды по г, находят
и(г,Л)= В + В /г+г2ф(г)/6,
(15)
где £1, £ - произвольные величины, зависящие от I
Согласуя (15) с граничными условиями (8), (9), получают
^ Я3+2/3 , ^ 13 ,,ч
В =--ф(0, в2 = -ф(0
1 6Я 2 3
и поэтому вместо (15) имеют
и(г, 0=ИТ)ф(0, (16)
где
Я3 + 213 I3 г2
^(г) = - -6+Я2- + >+т (17)
Подставляя (16) в (13), приходят к дифференциальному уравнению относительно искомой функции ф
Т-Ъё I -<г) Л
[ у/(г) ф 2с1г = Щ(Л\
или
где
А = Вф Л
(18)
А=-
1
13 - Я3
[ у/(г~)г2 Лг, (19)
I связано с Я и концентрацией с зависимостью (1). Интегрируя (18), получают
ф(Л) = В ехрфг / А), (20)
где £з - произвольная постоянная.
BecmnvKjBTy^T/Proceeémgs of VSUET, № 4, 2016L
Поэтому, в соответствии с (16), (20) имеют
u(r, t )= D3^(r)exp(D )
(21)
где iy( r ) задана по (17).
С целью определить D3, подставляют (21) в (12) 1
l3 - R3
J D3y(r )r2 dr = Ac,
откуда приходят к соотношению D А = Ac,
(22)
где А определяют по (20).
В результате чего, вследствие (21), (22) окончательно имеют
. . Ac . . .Dt. u(r, t )= —^(r)exp(—). А А
(23)
Как видно, что найденное по (23) выражение приведенной концентрации и(г, ^ имеет простую и удобную для количественного и качественного анализа мультипликативную форму. А именно, приведенная концентрация пропорциональна перепаду концентрации Ас, и, поскольку, согласно расчету, А < 0, р(г) < 0 при R < г < \, в соответствии с (7), (23) концентрация жидкостной среды положительна, убывает по экспоненте с течением времени t и, тем быстрее, чем больше коэффициент диффузии Б.
Если входящее в (23) выражение р(г) заменить его усредненным значением р = А согласно (19), то (23) упрощается и принимает вид
u = Ac exp( Dt ),
(24)
где Ас = Сп - сн.
На базе количественного анализа (24) судят о распределении концентрации в обессахариваемом утфеле в заданный период времени.
Поток концентрации сахарозы от жидкостной среды (утфеля) к поверхности шара (кристалла сахара) согласно (2), (23) (с точностью до знака) вычисляют по зависимости (м/с)
q(t ) = D
R3 -13
х
du( R, t) D
dr
3A X
л I Dt ■Ac ■exp I —
(25)
Полный удельный объем Q сахарозы, подводимый к единице поверхности шара за время Т, получают, интегрируя (25) по данному временному интервалу Т (м3/м2)
Q = J q(z)dz =
( l3 - R3) Aс
3R
1 - exp
DT
(26)
В соответствии с (26) предельное значение подводимого к единице поверхности шара объема сахарозы составляет
Qnped = limrQ =
(l3 - R3) Ac
(27)
пред ^да ^ З./^
Тогда если значение периода времени т между двумя раскачками утфеля пересыщенным раствором определяют исходя из того, чтобы выраженный в долях бпред коэффициент пересыщения к рассчитывался по отношению
а'апРед=к , (28)
откуда, исходя из (26)-(28), получают расчетное значение периода (мин) между двумя раскачками утфеля
т = —А— 1п(1 -к), (29)
600
где А вычисляют по (19).
Численный эксперимент
Рассматривается пример расчета периода времени между двумя раскачками утфеля са-харсодержащим раствором.
Коэффициент диффузии принимаем О = 2-10 - 10 м2/с, что примерно соответствует значению данного коэффициента при температуре раствора 70 °С. Полагали, что кристаллы сахарозы в растворе монодисперсны, диаметр частиц d = 20...40 мкм, объемная концентрация частиц £ = 40.50%.
По выбранным таким образом исходным параметрам, на базе (29), параметрические по размеру d частиц и их объемной концентрации £
рассчитаны графические зависимости периода т между двумя раскачками утфеля от коэффициента пересыщения к (рисунок 3).
20г
г1( к)
15
г2( к) г3(к)10 г4( к)
1 1 1 4
- 2 -
3 ~
4
— _ . .-»■""" , . ' _
1 1 I 1
0
0.9 0.92 0.94 0.96 0.98 1 к
Рисунок 3. Зависимости периода г (мин) между двумя раскачками утфеля от коэффициента пересыщения к для различных значений диаметра d (мкм) и объемной концентрации £ (%) частиц (d = 20 мкм: 1 - £ = 40, 2 - £ = 50%; d = 40 мкм: 3 - £ = 40, 4 - £ = 50%)
Figure 3. Dependences of the period г (min.) between two swings massecuite from supersaturation coefficient rfor various values of diameter d (mcm) and volume concentration £ (%) of particles (d = 20 мкм: 1 - £ = 40, 2 - £ = 50%; d = 40 мкм: 3 - £ = 40, 4 - £ = 50%)
0
BecmnvKjBTy^T/Proceedings of VSUET, № 4, 2016
Из вида графиков рисунка 3 выявляется согласие полученных результатов количественного моделирования рассматриваемого процесса с физическим смыслом исследуемого явления. А именно, естественное увеличение периода т раскачки в вакуум-аппарате когда возрастает коэффициент пересыщения к (все кривые рисунок 3). Кроме того, отмечается уменьшение периода раскачки для частиц одного размера вместе с ростом их концентрации (кривые 2 и 4 расположены, соответственно, ниже кривых 1 и 3), а также в случае, когда при одинаковой концентрации £ увеличивается размер частиц.
ЛИТЕРАТУРА
1 Saska M. Modelling crystallization and inclusion formation in sucrose crystals // Zuckerindustrie. 2011. V. 113. № 3. Р. 224-229.
2 Tuzhilkin V., Urazbayeva K., Alibekov R. Semicontinuos sugar crystallization // Jour. of industrial technology and engineering. 2012. № 2 (3). P. 11-18.
3 Beckmann W. Crystallization: basic concepts and industrial applications. Weincheim: Wiley-VCH Verlag GmbH & Co, 2013. 346 p.
4 Гнездилова А.И, Червецов В.В. Интенсификация процессов кристаллизации при производстве молочных продуктов. Воронеж: "Издательство Ритм", 2011. 196 с.
5 Семенов Е.В., Славянский А.А., Сергеева Е.А., Щитова Т.А. Особенности диффузионного процесса кристаллизации сахарозы // Сахар. 2013. № 3. С. 46-50.
6 Тужилкин В.И., Прокофьев Е.Н., Коваленок В.А., Клемешов Д.А. и др. Технология полунепрерывной кристаллизации сахара // Сахар. 2014. № 2. С. 30-33.
СВЕДЕНИЯ ОБ АВТОРАХ Анатолий А. Славянский д. т. н., заведующий кафедрой, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Земляной вал, 73, г. Москва, 109004, Россия, anatoliy4455@yandex.ru Евгений В. Семенов д. т. н., профессор, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Земляной вал, 73, г. Москва, 109004, Россия, sem-post@mail.ru Наталья Н. Лебедева к. т. н., старший преподаватель, кафедра технологии продуктов из растительного сырья и парфюмерно-косметических изделий, Московский государственный университет технологий и управления им. К.Г. Разумовского (ПКУ), ул. Земляной вал, 73, г. Москва, 109004, Россия, mgutu-sahar@mail.ru Сергей Т. Антипов д. т. н., профессор, проректор по научной и инновационной деятельности, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия, ast@vsuet.ru
КРИТЕРИЙ АВТОРСТВА Анатолий А. Славянский предложил научную концепцию работы, корректировка рукописи статьи
Евгений В. Семенов обоснование физико-математической модели процесса, несет ответственность за плагиат
Наталья Н. Лебедева обзор литературных источников по исследуемой проблеме, выполнила численный эксперимент и его графическую иллюстрацию Сергей Т. Антипов консультация в ходе исследования КОНФЛИКТ ИНТЕРЕСОВ Авторы заявляют об отсутствии конфликта интересов.
ПОСТУПИЛА 14.11.2016 ПРИНЯТА В ПЕЧАТЬ 01.12.2016
Заключение
На основе модели молекулярного массопереноса сахарозы от межкристального раствора к коллективу шарообразных частиц, имитирующих кристаллы сахара в вакуум-аппарате, предлагается аналитический аппарат по расчету периода времени между раскачками утфеля как функции коэффициента пересыщения и основных физико-механических параметров процесса кристаллизации сахарозы.
REFERENCES
1 Saska M. Modelling crystallization and inclusion formation in sucrose crystals. Zuckerindustrie, 2011, vol. 113, no. 3, pp. 224-229.
2 Tuzhilkin V., Urazbayeva K., Alibekov R. Semicontinuos sugar crystallization. Jour. of industrial technology and engineering, 2012, no. 2 (3), pp. 11-18.
3 Beckmann W. Crystallization: basic concepts and industrial applications. Weincheim, Wiley-VCH Verlag GmbH & Co, 2013. 346 p.
4 Gnezdilova A.I, Chervetsov V.V. Intensifikaciya processov kristallizacii pri proizvodstve molochnyh produktov [An intensification of processes of crystallization by production of dairy products]. Voronezh, "Publishing House Rhythm". 2011, 196 p. (in Russian)
5 Semenov E.V., Slavyanskii A.A., Sergeeva E.A., Shchitova T.A Features ofthe diffusionprocess of crystallization of the sucrose. Sakhar [Sugar] 2013, no. 3, pp. 46-50. (in Russian).
6 Tuzhilkin V.I., Prokofev E.N., Kovalenok V.A., Klemeshov D.A. et al. Semi-continuous sugar crystallization technology. Sakhar [Sugar] 2014, no. 2, pp. 30-33. (in Russian).
INFORMATION ABOUT AUTHORS
Anatolii A. Slavyanskii doctor of chemical sciences,, head of the department, technologies of products from vegetable raw materials and perfume-cosmetic products department, K.G. Razumovsky Moscow State University of technologies and management (the First Cossack University), Zemlyanoi val str, 73, Moscow, 109004, Russia, anatoliy4455@yandex.ru Evgenii V. Semenov doctor of chemical sciences, professor, technologies of products from vegetable raw materials and perfume-cosmetic products department, K.G. Razumovsky Moscow State University of technologies and management (the First Cossack University), Zem-lyanoi val str, 73, Moscow, 109004, Russia, sem-post@mail.ru Nataliya N. Lebedeva candidate of technical sciences, chief lecturer, technologies of products from vegetable raw materials and perfume-cosmetic products department, K.G. Razumovsky Moscow State University of technologies and management (the First Cossack University), Zem-lyanoi val str, 73, Moscow, 109004, Russia, mgutu-sahar@mail.ru Sergei T. Antipov doctor of technical sciences, professor, vice-rector for research and innovation, Voronezh state university of engineering technologies, Revolution Av., 19, Voronezh, 394036, Russia, ast@vsuet.ru
CONTRIBUTION Anatolii A. Slavyanskii proposed a scientific concept, correction of the manuscript
Evgenii V. Semenov rationale physic-mathematical model of the process, is responsible for plagiarism
Nataliya N. Lebedeva review of the literature on an investigated problem, performed numerical experiment and graphical illustration Sergei T. Antipov consultation during the study
CONFLICT OF INTEREST The authors declare no conflict of interest. RECEIVED 11.14.2016 ACCEPTED 1.12.2016
