№ 10 (67) октябрь, 2019 г. Потери давления ΔРж для преодоления гидравли- В уравнение (14) входит 2 неизвестных пара- ческого сопротивления отверстий при входе жидко- метра: w0 и w1. Для их нахождения необходимо обес- сти в газораспределительную насадку можно опреде- печивать работу экстрактора с заранее известной ве- лить, как личиной газовой подушки «h». Согласно данным [5], для равномерного поступления газа во все смеси- Pc вх wc2 c (13) тельные элементы экстрактора величина газовой по- 2 душки должна быть равна h=5dг, где dг – диаметр га- зораспределительной насадки. где е› – коэффициент сопротивления при входе жид- Проведенные нами экспериментальные исследо- кости в газораспределительную насадку. вания [1; 2] позволили установить основные пара- Выполнив необходимые математические преоб- метры, входящие в уравнение (14). разования, получаем следующее уравнение для рас- Коэффициенты сопротивления отверстий для по- чета величины газовой подушки: дачи газа в смесительные зоны экстрактора можно определить от зависимости [4; 5]: h o w02 Г 2c g 1 1 H 1w12 Г вх wc2 c (м) 2 2g c Г 1 0,5 1 0,37 1 0,37 d0 , ж 0 ж 0 0 (14) (16) где ρг – плотность газа, кг/м3; ρж – плотность жидко- где σ – поверхностное натяжение экстрагируемой сти, кг/м3; g – ускорение силы тяжести, 9,8 м/с2; Н – длина зоны смешения кольцевого канала, м; φ1 – га- жидкости, н/м; σв –поверхностное натяжение воды, зосодержание в кольцевом канале; ξ0 – коэффициент н/м; τ0 – коэффициент, зависящий от формы отвер- сопротивления отверстий газораспределительной стия; λ – коэффициент трения; δ – толщина стенки га- насадки; w0 – скорость газа в отверстии газораспре- делительной насадки, м/с; w1 – скорость газа в отвер- зораспределительного элемента, м; d0 – диаметр от- стии газопроводящего патрубка, м/с; wж – скорость верстия для прохода газа, м. жидкости в газораспределительной насадке, м/с; е› Значения τ0, ξ0 ,ξ1 в зависимости от δ/d0 можно определить по полученным экспериментальным гра- – коэффициент сопротивления входу жидкости в га- фикам, представленным на рис. 5 и рис. 6. зораспределительной насадке; 1 – коэффициент со- График влияния отношения площадей отверстий противления отверстия газопроводящего патрубка. к площади барботажной трубы в зависимости от по- верхностного натяжения представлен на рис. 7. Рисунок 5. Зависимость 0 от / d0 50
№ 10 (67) октябрь, 2019 г. Рисунок 6. Зависимость коэффициентов сопротивления на выходе газа из отверстий от соотношения δ/d Рисунок 7. Влияние поверхностного натяжения на коэффициент сопротивления входу жидкости в барботажный патрубок Полученные результаты. изменение газовой подушки. Эксперименты проводились отдельно для каждого отверстия во Для проверки уравнения 14 газовой подушки «h» внутренней зоне перемешивания d0 = 1,5 и d1=2,3 мм. в результате теоретических исследований было При различных значениях расхода газа Qг = 0,15; выполнено исследование экспериментальной 0,25; 0,35; 0,45; 0,55; 0,65; 0,75, подаваемого через установки экстрактора (рис. 4). отверстие d0 = 1,5 мм во внутренней зоне смешения, соответственно изменялось значение газовой Экспериментальные исследования подушки h = 6; 14; 31; 52; 74; 101; 122 мм. проводились в следующем порядке: При этом не удалось определить значение а) Отверстия для подачи газа внешней зоны газовой подушки со значением Qг = 0,85 м3/час расхода газа. Причина этого – длина смешения аппарата закрыты. Газ подавался только к газораспределительной насадки была недостаточной. внутренней зоне смешения через отверстия d = 1,5; Значение газовой подушки в соответствии с 2,3 мм. Значение расхода газа составляет до Qг = вышеупомянутыми значениями расхода газа через 0,15÷0,850 м3/час, передается с шагом 0,1 м3/час, а отверстие d1=2,3 мм составляет, h = 2; 3; 5; 8; 13; 17; значение расхода жидкости составляет до 28; 39 мм. Проверка показала, что увеличение QС=0,070,39 м3/час 0,08 м3/час, передается с шагом расхода жидкости на аппарате не повлияло на 0,08 м3/ час. В качестве легкой жидкости была взята газовую подушку. вода. Для каждого значения расхода жидкости Qж = 0,07; 0,23; 0,31; 0,39 м3/час, подавался отдельный расход газа со значением Qг = 0,15; 0,25; 0,35; 0,45; 0,55; 0,65; 0,75; 0,85 м3/ч, были проверены на 51
№ 10 (67) октябрь, 2019 г. Как видно из приведенных выше экспериментов, диаметром d0=1,5мм. Соответствующее высоту при подаче газа через отверстие d0 = 2,3 мм во газопроводящего патрубка во внешней зоны внутренней зоне смешения аппарата максимальное значение газовой подушки составляло h = 39 мм. смешения Н1=315; 240; 165мм, было открыто по порядку отверстии d1=0,6; 0,8; 1,0 мм. Если открыть отверстия для подачи газа в кольцевом канале, значение газовой подушки будет Отдельные эксперименты были проведены для недостаточным. каждого отверстия. Расход жидкости зоны Поэтому мы не будем использовать это отверстие смешивания передан со значениями Qж = 0,07; 0,15; d0=2,3 мм на следующих этапах эксперимента. 0,23; 0,31; 0,39 м3/час (с шагом 0,08 м3/час). Расход Результаты вышеупомянутых экспериментов газа 4 разных значений для каждого из значений были обработаны на ЭВМ, и было получено уравне- ние регрессии для расчета газовой подушки с расхо- расхода жидкости Qг = 0,55; 0,65; 0,75; 0,85 м3/час, дом газа через отверстие d0 = 1,5 мм [2]. определены экспериментальные значения газовой Полученное уравнение регрессии выглядит сле- подушки. дующим образом: Результаты экспериментов были обработаны на y0= -11,576+79,32x+136,5x2 (18) основе компьютерной программы, и получено где х – расход газа, м3/час; у0 – соответствующая рас- уравнение регрессии. Ниже приведены уравнения ходу газа величина газовой подушки, м. регрессии. Экспериментальные результаты, полученные на Размер отверстия для подачи газа во внутреннюю основе этого уравнения, были теоретически проверены с использованием уравнения (1), и был зону смешивания составляет d0 = 1,5 мм, в то время построен сравнительный график (рис. 8). как отверстие для подачи газа во внешнюю зону Экспериментальные исследования полностью подтвердили это уравнение; смешивания составляет d1 = 0,6 мм; б) Было открыто газопроводящее отверстие во y1*= -36,89+137,96х+28,92х2 (19) внутреннюю зону перемешивания аппарата При d0=1,5 мм и d1=0,8мм; (20) (21) y2*= 130,2 - 368,15х+346,28х2 При d0=1,5мм и d1=0,8мм; y3*= 89,75 - 316,46х+298,2х2 Рисунок 8. Зависимость величины газовой подушки от расхода барботирующего газа 52
№ 10 (67) октябрь, 2019 г. Рисунок 9. Зависимость величины газовой подушки от расхода газа (d0=1,5 мм – cоnst) 1; 1’ - d1=0,6 мм; Н1=315 мм В зависимости от заданного расхода газа были 2; 2’ - d1=0,8 мм; Н1=240 мм рассчитаны значения скоростей газа, проходящего 3; 3’ - d1=1,0 мм; Н1=165 мм через отверстия для подачи газа в зоны смешения. Теоретическое значение газовой подушки На основании проведенных экспериментальных находилось по уравнению (14) газовой подушки. исследований был определен расход распределяемого газа к внутренним и внешним зонам Теоретически найденные значения сравнивались смешения аппарата. с экспериментальными значениями (рис. 9). Эксперименты полностью подтвердили уравнение Для этого мы вычтем значения газовых подушек (14). Погрешность между теоретическими и на графике с рисунка 8, которые соответствуют экспериментальными значениями высоты газовой расходу газа на графике на рисунке 9, то есть h1 = h0 подушки не превышает ± 10%. – h1* мм; h2 = h0 – h2* мм; h3 = h0 – h3* мм. Эта величина представляет собой количество газа, Вывод распределенного во внешней зоне смешивания. Теоретические и экспериментальные Расход газа, распределяемого во внешнюю зону исследования полностью подтвердили уравнения для смешения, равен расчета газовой подушки (14). В результате оказалось возможным определить величину расхода QТ Qобш Q0 , м3/ час (22) газа при равномерной и устойчивой подаче инертного газа в зоны смешения барботажного Для каждого случая экспериментов, проводимых экстрактора при постоянных размерах газовых одним и тем же способом, распределяемое значение подушек. Уравнение 14 можно использовать для расхода газа рассчитано на основе общего значения проектирования промышленного устройства расхода газа в зонах смешения аппарата. барботажного экстрактора. Список литературы: 1. Каримов И.Т., Алиматов Б.А., Тожиев Р.Ж. Гидродинамика барботажных экстракторов // Научно-техниче- ский журнал Ферганского политехнического института. – 2000. – № 3. – С. 109-112. 2. Каримов И.Т. Исследование гидродинамических процессов в смесительных зонах барботажного экстрактора: Дис. канд. техн. наук. – Ташкент: Ташкентский гос. техн. ун-т, 2001. – 131 с. 3. Многоступенчатый барботажный экстрактор / Б.А. Алиматов, В.Н. Соколов, Х.М. Саъдуллаев, И.Т. Каримов // А.С. № 1607859 (СССР), БИ № 43, 1990. 4. Соколов В.Н., Геллис Ю.К. Гидродинамика барботажного кожухотрубного реактора // Химическая промыш- ленность. – 1962. – № 10. – С. 757-761. 5. Соколов В.Н., Доманский И.В. Газожидкостные реакторы. – Л.: Машиностроение, 1976. – 216 с. 53
№ 10 (67) октябрь, 2019 г. ИССЛЕДОВАНИЕ ГИДРАВЛИЧЕСКОГО СОПРОТИВЛЕНИЯ РОТОРНО-ФИЛЬТРУЮЩЕГО АППАРАТА Исомидинов Азизжон Саломидинович докторант (PhD), Ферганский политехнический институт, Республика Узбекистан, г. Фергана E-mail: [email protected] RESEARCH OF HYDRAULIC RESISTANCE OF ROTOR-FILTER APPARATUS Azizjon Isomidinov PhD, Fergana polytechnic institute, Republic of Uzbekistan, Fergana АННОТАЦИЯ В статье приведены результаты экспериментов по определению гидравлического сопротивления пылеулавливающего роторно-фильтрующего аппарата с мокрым способом. Предложена формула для определения коэффициента гидравлического сопротивления в зависимости от соотношения активной и пассивной поверхностей фильтрующего материала сетки на роторно-фильтрующем аппарате. На основе опытных данных общее потерянное давление на аппарате определено в зависимости от диаметра фильтрующей сетки. ABSTRACT In the article, the results of experiments to determine the hydraulic resistance of a dust-collecting rotor-filter apparatus with a wet method were presented. A formula is proposed for determining the hydraulic resistance coefficient depending on the ratio of the active and passive surfaces of the filter material of the mesh on the rotor-filter apparatus. Based on the experimental data, the total lost pressure on the apparatus is determined depending on the diameter of the filter section. Ключевые слова: активная и пассивная поверхность, ротор-фильтр, гидравлическое сопротивление, скорость потока, запыленный воздух, вторичный газ, коэффициент сопротивления, диффузор, конфузор, фильтрующая сетка. Keywords: active and passive surface, rotor filter, hydraulic resistance, flow rate, dusty air, secondary gas, drag coefficient, diffuser, confuser, filter mesh. _________________________________________________________________________________ _______________ Введение. Выбор оптимальной нагрузки, созданной на рабочих органах аппарата, влияет на скорость запыленного воз- духа, движущегося на пылеулавливающих аппаратах с мокрым способом, считается одной из основных характе- ристик. Эта характеристика определяет гидравлическое сопротивление аппарата и допустимое количество рабочей жидкости, выходящей из него [1, 2, 3]. Кроме этого, повышение гидравлического сопротивления на рабочем органе аппарата влияет на улучшение эффективности улавливания, уменьшение производительности аппарата и затекаемости частиц пылей в трубах. В свою очередь, это приводит к увеличению расхода энергии. Исследования были проведены с целью определения оптимальных значений этих факторов [6, 8, 9 и др.]. Однако полное решение для достижения максимальной эффективности очистки и производительности при минимальных значениях гидравлического сопротивления не было разработано. Общее исследование проблем устройства с роторным фильтром для очистки газов от пыли и запыленного воздуха мокрым способом [10] проведено для определения влияния на гидравлическое сопротивление устройства, на эффективность и результативность очистки. На рисунке 1 дан обзор экспериментального устройства. __________________________ Библиографическое описание: Исомидинов А.С. Исследование гидравлического сопротивления роторно-филь- трующего аппарата // Universum: Технические науки : электрон. научн. журн. 2019. № 10(67). URL: http://7universum.com/ru/tech/archive/item/7926
№ 10 (67) октябрь, 2019 г. Рисунок 1. Общий вид экспериментального аппарата: 1 – корпус аппарата; 2 – диффузор; 3 – ротор-фильтр; 4 – конфузор; 5 – ванна для рабочей жидкости; 6 – отвод для шлама; 7 – штуцер для выхода жидкости; 8 – зонт; 9 – патрубок Прандтля; 10 – электромотор; 11 – насос; 12 – вентилятор; 13 – управляющий аппарат скорости электромотора; 14 – регулятор рабочей жидкости Эксперименты проводились на приведенных где ∆Pоб – общее гидравлическое сопротивление ниже параметрах в пределах: диаметры отверстия аппарата, Па; Pдиф –гидравлическое сопротивление в фильтрующего материала покрытия на ротор dф= 2; диффузоре аппарата, Па; Pконф –гидравлическое 3; 4 мм; частота вращения ротора – 15; 25; 35 об/мин; сопротивление в конфузоре аппарата, Па; Pрп – скорость газа ʋг= 7,67÷34,4 м/с (указанный диапазон гидравлическое сопротивление в рабочих скоростей выбран вблизи с диапазоном скоростей, применяемых в промышленности). Эксперименты поверхностях аппарата, Па. проводились на жидкой и газовой системах при температуре 20±20С . Pрп определяется по следующей формуле: Объект и метод исследования. Pрп P1 P2 , Па (2) Потеря давления наблюдается при перемещении запыленного воздуха на всех типах где P1 – гидравлическое сопротивление в рабочем воздухоочистительных аппаратов. Это объясняется органе А, Па; P2 – гидравлическое сопротивление в конструктивной структурой аппарата и количеством рабочих органов. рабочем органе Б, Па. Эксперименты для Разработанный аппарат состоит из диффузора (2), рабочих поверхностей (А и Б), фильтра (3), определения общего гидравлического конфузора (4). Движущийся запыленный поток воздуха проходит через диффузор (2), поток, сопротивления роторно-фильтрующего аппарата расширяясь на рабочих поверхностях (А и Б) фильтра, переходит на линейное движение. проводились в неоросительном состоянии очистки. Прямолинейное движение потока осуществляется распределительным зонтом (8) и жидкостью, Общее гидравлическое сопротивление не в очищающейся в ванне (5). Очищенный воздух на поверхности фильтра (3) выходит через конфузор (4). оросительном состоянии ротор-фильтрующего Потеря давления происходит на трех рабочих органах аппарата. аппарата определяется по формуле Дарси – Тогда общее гидравлическое сопротивление аппарата можно записать следующим образом [4]: Вейсбаха: Pсух об г 2 г , Па (3) 2 Pоб Pдиф Pрп Pконф , Па (1) где: газ – плотность газа, кг/м3; газ – потерянная скорость газа в рабочих органах аппарата, сух – коэффициент общего сопротивления аппарата. Определение и расчет коэффициентов сопротивления имеет сложную природу и требует различных отклонений, которые могут быть определены только путем экспериментов. 55
№ 10 (67) октябрь, 2019 г. Коэффициент гидравлического сопротивления расширения диффузора; конф – коэффициент роторно-фильтрующего аппарата определяется вышеуказанными факторами и вычисляется по сопротивления конфузора, и это определяется по следующей формуле [4]: следующей формуле: об диф рп конф (3) конф конф 1 (5) 8sin тб 1 n2 где диф – коэффициент сопротивления определяется приблизително по формуле И.Е. Идельчика и конф Флингера [7]. 2 диф диф 1 sin кб 1 2 (4) где конф – коэффициент трения конфузора; ус – 8sin кб 1 n2 2 1 nдиф угол сжатия конфузора; nконф – степень сжатия диф конфузора. 2 Коэффициенты сопротивления в диффузоре и где диф – коэффициент трения в диффузоре; nдиф – конфузоре роторно-фильтрующего устройства были степень расширения диффузора; ур – угол определены экспериментально в соответствии с приведенными выше формулами (4) и (5). Значения коэффициентов гидравлического сопротивления приведены в таблице 1. Таблица 1. Значения коэффициентов гидравлического сопротивления υгаз ,м/с 34,4 28 23,8 18,79 7,67 диф 0,3223 0,3085 0,3127 0,2978 0,3117 0,2154 0,2081 0,2113 0,2205 0,2178 конф Из значений, приведенных в таблице 1, видно, Sакт 2,38 показатель равен k 0,91, при что коэффициенты местного сопротивления в Sпас диффузорах и конфузорах в незначительной степени влияют на движущийся воздушный поток. Одной из значениях Sакт 4,7 показатель равен основных причин этого, можно сказать, является Sпас снижение силы трения за счет резкого расширения диффузора и резкого сужения конфузора. Тогда k 0, 48 ). можно считать средним коэффициент локального сопротивления в диффузоре и конфузоре равным 0,5 Это связано с соотношением общего с достаточной точностью. коэффициента сопротивления к активной поверхности (рис. 2). рп – коэффициент сопротивления рабочей поверхности, определение зависит от отношения между активными и пассивными поверхностями фильтрующего материала [5]. рп k Sакт , (5) Sпас где Sакт – поверхность фильтрующего материала аппарата без отверстий, Sпас – поверхность фильтрующего материала аппарата с отверстиями, k – поправочный коэффицент, определяется экспериментально. k значения были определены экспери- Рисунок 2. График зависимости k от Sакт Sпас ментально в зависимости от отношения активной поверхности к пассивной поверхности (рис. 2). Как показано на рисунке 2, увеличение активной поверхности приводит к уменьшению поправочного коэффициента (при значениях Sакт 1,65 В этом случае формула (3) записывается в Sпас следующем виде показатель равен k 1,1, при значениях 56
№ 10 (67) октябрь, 2019 г. об 0,5 k Sакт (6) Sпас Общий коэффициент гидравлического сопротив- ления для неорошаемого аппарата определяется по формуле (6). Общий коэффициент сопротивления опреде- лялся как соотношение активной поверхности сетки к пассивной поверхности сетки фильтра (рис. 3). Рисунок 3. График зависимости об от Sакт Рисунок 4. График зависимости Pсух от г на Sпас сухом аппарате Как видно из рисунка 3, увеличение активной поверхности приводит к увеличению сопротивления 1- Sакт 0.268м2 ; 2- Sакт 0.229м2 ; устройства (при Sакт 0, 268м2 об 2,1; при 3- Sакт 0.202м2 . Sакт 0, 229м2 об 2, 2 ; при Sакт 0, 202м2 В роторном фильтрующем устройстве распыли- ваемая жидкость из отверстия щтуцера (7) через зонт об 2, 6 ). (8) на активные и пассивные поверхности фильтра образует пленку на поверхности фильтра (3). Используя полученные экспериментальные Расстояние близости ротора-фильтра (3) и зонта (8) позволяет проходящему запыленному воздуху быть данные, гидравлическое сопротивление перпендикулярно направленым через пленки жидкости на активной и пассивной поверхностях. неорошаемого аппарата можно определить по Пленка жидкости контактирует с запыленным воздухом, вследствие этого улавливаются частицы формуле (3). Небольшая частота вращения ротора не пыли. Рабочая жидкость, собранная в жидкостной ванне (5) под устройством, будет вымывать частицы создает эффект вентиляции, однако вследствие этого пыли, удерживаемые в поверхностном слое роторного фильтра, и очищать фильтр до рабочего создается небольшое гидравлическое сопротивление состояния. потока воздуха. Следовательно, влияние частоты Толщина слоя пленки зависит от количества подаваемой жидкости. Чем толще слой, тем выше вращения ротора не учитывается при расчете эффективность очистки. Но при этом увеличивается гидравлическое сопротивление. Исследование слоя гидравлического сопротивления. Вместе с тем пленки представляет собой сложный процесс, который зависит от многих факторов. Поскольку частота вращения ротора обеспечивает быстрое исследование толщины пленки не было главной обновление рабочей поверхности и увеличение целью работы, то гидравлическое сопротивление поверхности контакта. устройства определялось количеством потребляемой жидкости (рис. 5, 6, 7). Эксперименты выявили Данные, представленные на рисунке 4, высокие гидравлические нагрузки в аппарате в показывают, что при небольших значениях скорости зависимости от максимальной работоспособности воздуха, передаваемой в устройство, расход невелик, жидкости, подаваемой в аппарат. а потеря давления в устройстве близка друг к другу. При увеличении скорости воздуха поток Для экспериментов были выбраны щтуцеры с уплотняется, и вследствие этого увеличивается диаметром отверстия dш=1; 2; 3 мм, которые гидравлическое сопротивление. установлены для распыления жидкости фильтрующега ротора. Гидравлическое сопротивление при орошении определялось в диапазоне скоростей газа 7,67- Результаты проведенных экспериментальных 34,4 м/с. исследований приведены в следующих графиках: 57
№ 10 (67) октябрь, 2019 г. Рисунок 5. График зависимости Pорш от г c Рисунок 7. График зависимости Pорш от г c орошением аппарата при диаметре штуцера орошением аппарата при диаметре штуцера dш=1мм-const dш=3 мм-const Среднеквадратическое отклонение полученных экспериментальных данных составляет 1. Отклонение по газу – 12,7%; 2. Отклонение по изменению давления – 10%. Экспериментальные исследования (рис. 5, 6, 7) показывают, что гидравлическое сопротивление увеличивается с увеличением потока жидкости к роторно-фильтрующему аппарату. Это позволяет повысить эффективность очистки аппарата. По результатам экспериментальных исследований удалось определить основные рабочие параметры роторно-фильтрующего аппарата. Рисунок 6. График зависимости Pорш от г c орошением аппарата при диаметре штуцера dш=2 мм-const Список литературы: 1. Алиев Г.М. Техника пылеулавливания и очистки промышленных газов. – М.: Металлургия, 1986. – 544 с. 2. Andersen B.O., Nielsen N.F., Walthe J.H. Численное и экспериментальное исследование импульсно-струйной очистки в тканевых фильтрах: Дисс. канд. филос. наук. – Где издано, 2017. – 203 с. 3. Валдберг А.Ю., Николайкина Н.Е. Процессы и аппараты защиты окружающей среды. – М.: Дрофа, 2008. – 239 с. 4. Исомидинов А.С., Мадалиев А.Н. Гидродинамика и аэродинамика ротор- фильтрного аппарата очищающего запыленных газов // Междунар. науч.-практ. конф. «Международный научный обзор проблем и перспектив современной науки и образования» (22-23 ноября 2018 г., Бостон, США). – Бостон, 2018. – 32 с. 5. Исомиддинов А.С., Каримов И.Т., Тожиев Р.Ж. Определение активной и пассивной поверхностей ротор- фильтрного аппарата для фильтрации запыленного воздуха // I Междунар. науч.-практ. конф. «Актуальные проблемы внедрения инновационной техники и технологий на предприятиях по производству строительных материалов, в химической промышленности и в смежных отраслях» (24-25 мая 2019 г., Фергана). – Фергана, 2019. – Т. 3. – С. 429-431. 6. Кабалдин Г.С. Модернизация распылительных и барабанных сушильных установок / Г.С. Кабалдин. – М.: Металлургия, 1991. – 205 с. 7. Латипов К.Ш. Гидравлики, гидромашины и гидроприводы. – Ташкент: Учитель, 1992. – С. 75-80. 8. Нечаева Е.С. Исследование основных характеристик роторного распылительного пылеуловителя: Дисс канд. техн. наук. – Кемерово, 2014. – 149 с. 9. Сорокопуд А.Ф. Совершенствование мокрого пылеулавливания в технологии сыпучих пищевых продуктов / А.Ф. Сорокопуд, М.И. Даниленко // Хранение и переработка сельхозсырья. – 1997. – № 4. – С. 3-8. 10. Тожиев Р.Ж., Каримов И.Т., Исомидинов А.С. Роторный мокрый пылеуловитель // ФерПИ: научно- технический журнал. – 2018. – Вып. 1. – С. 195-198. 58
№ 10 (67) октябрь, 2019 г. ТЕХНОЛОГИЯ ПРОДОВОЛЬСТВЕННЫХ ПРОДУКТОВ БЕНТОНИТОВЫЕ АДСОРБЕНТЫ МОДИФИЦИРОВАННЫЕ РАСТВОРОМ КАРБАМИДА Ахмедов Азимжон Нормуминович канд. техн. наук, доцент Каршинского инженерно-экономического института, Узбекистан, г.Карши E-mail:[email protected] Абдурахимов Саидакбар Абдурахмонович д-р техн.наук, профессор Ташкентского химико-технологического института, Узбекистан, г.Ташкент E-mail:[email protected] Азимов Юсуф Хидирович ассистент Каршинского инженерно-экономического института, Узбекистан, г.Карши E-mail: [email protected] BENTONITE ADSORBENTS MODIFIED BY THE CARBAMIDE SOLUTION Azimjon Akhmedov Candidate of Technical Sciences, Associate Professor of Karshy Engineering-Economic Institute, Uzbekistan, Karshy Saidakbar Abdurakhimov Doctor of Technical Science, Professor of Tashkent Chemical-Technological Institute, Uzbekistan, Tashkent Yusuf Azimov Assistant of Karshy Engineering-Economic Institute, Uzbekistan, Karshy АННОТАЦИЯ В данной статье приведена информация о сочетании мочевины с госсиполом, различных местных видов бен- тонита и методах получения новых модифицированных бентонитов и основные показатели модифицированных адсорбентов и результаты использования их при очистке сырого форпрессового хлопкового масла. ABSTRACT The article deals with information on the combination of urea with gossypol, various local types of bentonite and methods to produce new modified bentonites and main indicators of modified adsorbents, and results of their use in the purification of crude forepress cottonseed oil. Ключевые слова: I-II и III-IV сортов семян хлопчатника. хлопковое масло, полученное из низкосортных семян хлопчатника, госсипол, показателей сырых форпрессовых масел, каратиноид, хлорофилл, бентонит. Keywords: I-II and III-IV varieties of cotton seeds; cottonseed oil, derived from low grade cotton seeds, gossypol, indicators of crude forepress oils, carotenoid, chlorophyll, bentonite. ________________________________________________________________________________________________ Введение: Сырые масла, получаемые из низко- цветность хлопковых масел, полученных из семян сортных (III и IV сорта) и нестандартных семян хлоп- III-IV сортов существенное влияние оказывает пиг- чатника имеют темный (или непросматриваемый) мент протоплазмы, локализованный в пластидах - цвет, который обусловлен высоким содержанием хлорофилл, который при омылении щелочью, дает госсипола, хлорофилла и их производных [1]. На __________________________ Библиографическое описание: Ахмедов А.Н., Абдурахимов С.А., Азимов Ю.Х. Бентонитовые адсорбенты моди- фицированные раствором карбамида // Universum: Технические науки : электрон. научн. журн. 2019. № 10(67). URL: http://7universum.com/ru/tech/archive/item/7951
№ 10 (67) октябрь, 2019 г. продукты реакции, окрашенные в бурый или желтый при хранении в условиях повышенной влажности и цвет. госсикаэрулин - продукт окисления госсипола, в ко- тором гидроксильная и хинонная группы находятся в Отличительной особенностью хлопковых масел смежных положениях. от других растительных жиров является наличие кра- сящих пигментов: госсипола и его производных В процессе хранения хлопковых семян увеличи- форм, которые локализуются в ядре семени в госси- вается их засоренность, возрастает кислотность и со- половых железках, пропитанных пектином, гемицел- держание связанных форм госсипола, снижается мас- люлозой и другими неидентифицированными веще- личность, доля водо- и солерастворимой части бел- ствами. ков, моносахаридов, количество липопротеиновых комплексов. С ухудшением сортности семян снижа- В железках, кроме госсипола содержатся еще ме- ется качество и потребительские свойства (рафини- нее изученные пигменты, построенные на основе руемость) хлопковых масел из-за значительного со- сходного с ним углеводородного скелета. Исследова- держания в них свободных жирных кислот, связан- ниями установлено, что в госсиполовых железках в ных и измененных форм госсипола и хлорофилла. зависимости от степени зрелости семян содержатся: свободный госсипол - 14,5-45,0%, связанный госси- В работах [1, 2] показана возможность рафина- пол - 0,57-2,0%, госсипурпурин - 0,15-6,3%, госси- ции хлопковых масел с высоким содержанием крася- вердурин - до 2,0 %. щих пигментов с использованием карбамида, где гос- сипол, взаимодействуя с ним, образует соединения Кроме перечисленных пигментов в хлопковых указанные в рис. 1. маслах содержатся госсифульфин - аминопроизвод- ное госсипола, образующийся в хлопковых семенах H-C=O OH N-CO-NH2 H-C OH HO NH2 HO 2H2O 2C O HO CH3 HO CH3 NH2 CH3-CH-CH3 2 CH3-CH-CH3 2 госсипол карбамид карбамид-госсипол Рисунок 1. Схема взаимодействия госсипола с карбамидом При этом, в масляной смеси происходит образо- концентрации). В обоих случаях трудность заключа- вание как монокарбамид-госсипола, так и дикарба- ется в разделении полученных продуктов и неудовле- мид-госсипола [2, 3]. творительном качестве получаемого соапстока, кото- рый трудно поддаётся разложению и т.д. С практической точки зрения ценным является возможность взаимодействия карбамида с некото- С другой стороны, карбамид затрудняет про- рыми темноокрашенными измененными формами мывку масла и увеличивает его потери с водой. госсипола, которые весьма трудно выводятся из со- става масла (или вовсе не удаляются из него) при Учитывая выше отмеченные недостатки нами, обычной щелочной рафинации. для предварительной рафинации масел, получаемых из низкосортных и нестандартных семян хлопчат- Выход и качество рафинированного хлопкового ника предложены бентонитовые адсорбенты, моди- масла зависит ряда условний: температуры, концен- фицированные карбамидом. При этом в качестве гли- трации, условий ведения реакции, от способа взаимо- нистых адсорбентов использованы щелочные и ще- действия с щелочным реагентом в т.ч. и карбамидом. лочноземельные бентониты Навбахорского место- Карбамид может быть введен в масло в виде сухого рождения Навоийской области [4]. Химический со- порошка или его водного раствора (при оптимальной став данных бентонитов представлен в таблице 1. Таблица 1. Химический состав бентонитов Навбахорского месторождения Наименование компонентов Содержание, % в пересч. на абс. сух. веществ SiO2 щелочной бентонит щелочно-земельный бентонит Al2O3 Fe2O3 53,35 52,39 CaO TiO2 11,78 13,76 MgO 9,29 9,74 0,24 2,76 1,32 1,51 1,77 1,82 60
№ 10 (67) октябрь, 2019 г. Na2O 2,23 1,03 2,53 2,06 K2O 0,06 0,06 0,15 0,13 MnO 14,31 13,21 97,03 99,47 P2O п.п.п Сумма Из табл. 1 видно, что щелочно-земельном бенто- Результаты и обсуждения: В табл.2 представ- ните содержание СаО составляет 2,76%, а в щелоч- лены основные физико-химические и технологиче- ном-0,24%. ские показатели модифицированных бентонитовых адсорбентов раствором карбамида. Степень очистки Материалы и методы: Модификацию данных хлопкового масла от госсипола, каротиноидов и их бентонитов карбамидом осуществляли согласно производных, а также от хлорофилла, феофитина и блок-схеме, приведенной в рис. 2: их производных определялась по спектральному ме- тоду [5, 6]. Рисунок 2. Блок-схема получения бентонитового адсорбента модифицированного раствором карбамида Таблица 2. Основные показатели модифицированных адсорбентов и результаты использования их при очистке сырого форпрессового хлопкового масла Степень очистки масла Показатели адсорбента Вид от госсипола, от хлоро- масло- фильтруе- , удельная содержание адсорбента, каротиноидов и филла, и ёмкость, % мость масла поверх- структурной бентонит: их производных его произ- ность, 5х103 воды, моль/кг водных мл/5 сек Щелочной м2/кг Щелочно- земельный До обработки адсорбента (контроль) Щелочной 1,45 1,33 48,5 14,3 51,0 4,18 Щелочно- земельный 1,57 1,41 49,8 12,7 49,0 4,05 После обработки адсорбента 20% ным раствором карбамида (опыт) 5,05 1,98 1,54 46,3 14,0 53,1 2,15 1,66 47,5 12,3 51,2 5,13 Из табл.2 видно, что модифицированные карба- полученных из низкосортных и нестандартных семян мидом щелочные и щелочно-земельные бентонито- хлопчатника. Повышение степени очистка масла вые адсорбенты показывают более высокую степень наблюдается и по хлорофиллу, феофитинам и их про- очистки масла от госсипола, каротиноидов и их про- изводным (1,54 и 1,66), что также важно для предва- изводных (1,98 и 2,15). рительной рафинации сырых хлопковых масел. Дру- гие показатели бентонитовых адсорбентов до и после Это ещё раз подтверждает эффективность пред- лагаемого способа карбамидной рафинации масел, 61
№ 10 (67) октябрь, 2019 г. их обработки 20% ным раствором карбамида меня- натриевыми или другими солями, вступать во взаи- ются незначительно. модействие со свободным госсиполом с образова- нием оснований Шиффа, взаимодействовать с неко- Следовательно, бентонитовые адсорбенты, мо- торыми темноокрашенными измененными формами дифицированные карбамидом, могут быть использо- госсипола и таким путем удалить их из состава масла. ваны как эффективные реагенты при рафинации ма- сел, получаемых из низкосортных и нестандартных Таким образом, получение модифицированных семян хлопчатника. карбамидом бентонитовых адсорбентов и их приме- нение при рафинации сырых масел, получаемых из Выводы: Технологические свойства карбамида, с низкосортных и нестандартных семян хлопчатника точки его применения в масложировой промышлен- позволяет значительно повысить технико-экономи- ности, изучены авторами [1,2,7,8] и выявлены поло- ческие показатели рафинационного производства и жительные эффекты при очистки от сопутствующих качества выпускаемых продуктов. веществ. По данным, при рафинации хлопкового масла, карбамид проявляет способность образовывать со- единения со свободными жирными кислотами, их Список литературы: 1. Ильясов А.Т., Ураков Р.М. Решение актуальных проблем рафинации и демаргаринизации хлопкового масла: Ташкент, ТГТУ, 1996-85 с. 2. Ильясов А.Т., Пак В. Двухступенчатая рафинации высококислотного хлопкового масел // Пищевая промыш- ленность, 1991, №3, с. 66-68. 3. Пулатова С.Р. Повышение качестве экстракционного хлопкового масла и шрота. Автореф.канд.дисс. Таш- кент, 1989 – 16 с. 4. Мирзаев А.У., Черненко Т.В., Глушенкова А.И., Чинникулов Х. Сорбционные свойства бентонитовых глин Навбахарского месторождения //Узб.хим.журнал, 1999, №5-6, с. 34-36. 5. Абдурахимов С.А., Эргашева Д.К., Назиркулов А.Н. Спектральная оценка качества состава рафинированных масел //Изв.вузов, Пищевая технология, 1990, №6, с.77-78. 6. Абдурахимов С.А., Бахтияров С.Б., Салимов З.С., Латипов Б.Х., Тиллаева Г.У. Оценка степени осветления масла и избирательности адсорбента к сопутствующим триацилглицеринам веществ //Узб.хим.журнал, 1998, №1, с.66-69. 7. Ахмедов А.Н. Совершенстовование технологии комплексной рафинации масел, полученных из низкосорт- ных семян хлопчатника. Дисс. … канд. тенх. наук. - Карши, 2012. - 110 с. 8. Ахмедов А.Н. Исследование показателей хлопкого масла, полученного методом форпрессования из низкосортных семян хлопчатника //Universum журнал: Технический науки, 2019, Выпуск:№4(61).-М.: -23- 26 с. 62
№ 10 (67) октябрь, 2019 г. РАСЧЕТ КОНВЕКТИВНОЙ СУШИЛЬНОЙ УСТАНОВКИ С ПРИМЕНЕНИЕМ ТЕПЛОВЫХ ТРУБ Холиков Алижон Абдираупович ст. препод. Бухарского инженерно-технологического института, Республика Узбекистан, г. Бухара E-mail: [email protected] Мирзаева Шохиста Усмоновна докторант Бухарского инженерно-технологического института, Республика Узбекистан, г. Бухара E-mail: [email protected] Рахматов Умид Ражабович ассис. Бухарского инженерно-технологического института, Республика Узбекистан, г. Бухара E-mail: [email protected] CALCULATION OF CONVECTIVE DRYING UNIT WITH APPLICATION OF HEAT PIPES Alijon Kholikov senior lecturer, Bukhara Engineering and Technology Institute, Republic of Uzbekistan, Bukhara Shokhista Mirzaeva PhD student, Bukhara Engineering and Technology Institute, Republic of Uzbekistan, Bukhara Umid Rakhmatov assis. Bukhara Engineering and Technology Institute Republic of Uzbekistan, Bukhara АННОТАЦИЯ В статье изложены пути решения энергетических затрат с применением новых методов энергетического под- вода, обеспечивающих экономию энергии в процессе сушки, и разработана математическая модель расчета кон- вективной сушильной установки с применением тепловых труб. ABSTRACT The article describes ways to solve energy costs using new methods of energy supply, which provide energy savings during the drying process and a mathematical model for calculating a convective drying plant using heat pipes is developed. Ключевые слова: тепловая труба, тепловой насос, сушка, энергия, метод, переработка, масса, влага, температура, энтальпия, удельное. Keywords: heat pipe, heat pump, drying, energy, method, processing, mass, moisture, temperature, enthalpy, specific. ________________________________________________________________________________________________ Задача повышения энергетической эффективно- Плодоовощное сырье представляет собой неза- сти тепломассообменных установок может быть менимый источник важнейших физиологически ак- успешно решена, в частности, путем применения тивных веществ – витаминов, углеводов и минераль- тепловых насосов (ТН), тепловых труб (ТТ), позволя- ных веществ, необходимых для нормальной жизне- ющих создавать рациональные схемы использования деятельности человека. Вместе с тем в условиях жар- энергии, утилизировать вторичные энергоресурсы, кого летнего периода массового созревания и уборки применять нетрадиционные возобновляемые источ- многие виды плодов и овощей могут сохраняться не- ники энергии. долго. Длительное их хранение возможно только по- сле сушки. __________________________ Библиографическое описание: Холиков А.А., Мирзаева Ш.У., Рахматов У.Р. Расчет конвективной сушильной установки с применением тепловых труб // Universum: Технические науки : электрон. научн. журн. 2019. № 10(67). URL: http://7universum.com/ru/tech/archive/item/7989
№ 10 (67) октябрь, 2019 г. Для этого необходимо разработать высокоэффек- L W (7) тивный способ обработки, основанный на примене- d2 do нии новых методов энергетического подвода, обеспе- чивающих экономию энергии в процессе сушки и по- Удельный расход воздуха на 1 кг влаги: вышение качества продукта. Это могут быть тепло- массообменные установки с применением ТТ одно- L l или L 1 (8) целевого и комплексного назначения, особенностью W d2 do которых является максимальная утилизация вторич- ных энергоресурсов, использование нетрадиционных Влагосодержание воздуха do и d2 определяется возобновляемых источников энергии, применение ТТ для предварительного нагревания продукта. Для по I-d диаграмме Рамзина. этого нужно разработать математическую модель Для расчета сушилки необходимо знать пара- расчета конвективной сушильной установки с прменением ТТ [2; 7; 9; 10]. метры атмосферного воздуха: температуру to и Cперва нужно рассчитать материальный баланс. влажность o ; температуру воздуха на входе в су- Для определения количества влаги, удаляемой из ма- териала в процессе сушки, составляют материальный шилку t1 и один из параметров воздуха на выходе из баланс: сушилки: t2 или 2 . G1 G2 W (1) По заданным параметрам to и o находим на их пересечении точку А и определяем соответствующие где G1 u G2 – масса влажного и высушенного мате- ей значения влагосодержания do и энтальпии Io . В калорифере процесс нагревания воздуха проходит риала, кг; W – влага, удаленная из материала в процессе при постоянном влагосодержании do d1 const . Проведя через точку А вертикальную линию посто- сушки, кг/с. Материальный баланс по абсолютно сухому ве- янных влагосодержаний до пересечения с линией ществу, количество которого не меняется в процессе температуры t1 , получим точку В и определим соот- сушки: ветствующее ей значение I1 . Линия АВ характери- G1 100 Wн G2 100 Wk зует процесс нагревания воздуха в калорифере. 100 100 (2) В действительной сушилке из-за потерь тепла и изменения энтальпии воздуха в процессе сушки где Wн и Wk – содержание влаги во влажном и высу- определяется величиной и знаком поправки к рас- шенном материалах. чету действительной сушилки: Решая уравнения (1) и (2), определяют количе- ство испаренной влаги W и высушенного материала Коэффициент поправки определяется G2 : следующим формулой. W G1 G2 l J2g J1 (9) (10) или J2g J1 l (3) G2 G1 1 Wн 1 Wk (4) Если 0 то J2 J1 100 100 Если 0 то J2 J1 W G1 WнWk (5) 100 Wk Если 0 то J2 J1 Из уравнения материального баланса влаги Построение процесса сушки в I-d-диаграмме для действительной сушилки проводится в следующем можно определить удельный расход воздуха L на порядке: испарение 1 кг влаги из материала. При установив- 1. На линии берем произвольную точку и через нее проводим линию до пересечения с линией точки. шемся процессе сушки количество влаги, внесенной 2. На полученной линии откладываем в милли- в сушилку с материалом и воздухом, должно быть метрах отрезок, где m=1000 – отношение масштабов осей диаграммы. равно количеству влаги, унесенному из сушилки. 3. Точку В соединяем с точкой А, продолжая ее Уравнение баланса влаги имеет вид: до пересечения с заданной линией температуры или влажности, получаем конечную точку процесса [1; 3; G1Wн Ldo G2Wk Ld2 (6) 100 100 4-6; 8; 10]. где L – расход сухого воздуха, кг/с; dou d2 – влаго- содержание воздуха на входе и выходе из сушилки. Из уравнения (6) определим расход сухого воз- духа: 64
№ 10 (67) октябрь, 2019 г. Q5 L J2 (15) где J2 – энтальпия воздуха на выходе из сушилки (определяется по I-d- диаграмме). Тепло, унесенное из сушилки парами влаги, ис- парившейся из материала, Вт: Q6 W Jn (16) где J n – энтальпия пара, Дж/кг. Тепло, ушедшее с высушенным материалом, Вт: Q7 G2 C2 tM (17) Рисунок 1. I-d-диаграмма где G2 – производительность сушилки по высушен- Тепловой баланс конвективной сушилки имеет ному материалу, кг/с; вид: tM – теплоемкость сухого материала, Дж/кг К, оC Тепло, поступающее в сушилку с воздухом, Вт: Тепло, унесенное транспортными устройствами, Вт: Q8 GTРCtТР (18) Q1 L J0 (11) tTP t2 – температура транспортных устройств где J0 – энтальпия воздуха на входе в калорифер, на выходе из сушилки. Тепловые потери, Вт: определяется по I-d- диаграмме. Тепло, поступающее в сушилку с материалом, QПОТ Q1 Q6 (19) Вт: Q2 G2 C2 tm (12) Из уравнения теплового баланса определим ко- личество подводимого тепла: где C2 – удельная теплоемкость продукта, Дж.кг К; Q L J2 J0 W Jn Cвл tвл G2C2 t1M G2C2tM tm = 15-20оС – начальная температура материала GTPCTP ( tT1P tTP) QПОТ (20) Тепло, поступающее с испаряемой влагой мате- Действительный удельный расход тепла на 1 кг риала, Вт: испаренной влаги: Q3 W Cвл евл (13) где Cвл – удельная теплоемкость влаги, Дж.(кг К); tвл q Qкал 1000 J1 J0 l J1 J0 (21) W – температура испарения влаги, оС. d2 d0 Тепло транспортных устройств, Вт (определя- Расход тепла в калорифере, Вт: ется для конвейерных сушилок): Q4 GТР СТР tТР (14) Qкал L J1 J0 (22) где GТР – масса транспортных устройств, кг; CTP – где J0 и J1 – энтальпии воздуха на входе в калорифер удельная теплоемкость транспортных устройств, и на выходе из него (определяются по I-d-диа- Дж.кг К; tTP =15-20оС – температура транспортных грамме). устройств на входе в сушильную установку. Для определения конструктивных размеров су- Qкал – тепло, сообщенное воздуху в калорифере, шилок нужно учитывать и знать некоторые пара- Вт; метры. Габаритные размеры сушильной камеры (длина, диаметр) зависят от заданной производитель- Qдоп – дополнительно подведенное тепло, Вт. ности сушилки и продолжительности процесса сушки. С целью увеличения производительности су- Значения Qкал и Qgop определяются по равности шилки, обеспечения равномерности сушки и эконо- мичности установки нужно максимально заполнить между приходной и расходной частями теплового ба- объем камеры материалом. При сушке дисперсных ланса. материалов должно быть обеспечено обмывание ча- стиц материала воздухом, что определяет степень за- Тепло, унесенное сушильным агентом, Вт: грузки сушильной камеры. 65
№ 10 (67) октябрь, 2019 г. Для сушки кусковых (штучных) пищевых про- LТ n Lвaг Lo (25) дуктов применяют туннельные и камерные сушилки. Камерные сушилки периодического действия бы- Lвaг – длина вагонетки, м вают шкафного типа и со стационарными вагонет- ками. Величина Lo определяется по конструктивным Более совершенными являются установки тун- соображениям в зависимости от вместимости тун- нельного типа (конвейерные, ленточные, вагонеточ- неля и способа подачи воздуха. ные), которые работают как цикличные или непре- рывно действующие аппараты. Ширина туннеля: Особенностью расчета указанных установок яв- В Вваг 40 70мм ляется методика определения габаритных размеров сушильной камеры. Bваг – ширина вагонетки, м. Для расчета туннельной сушилки с вагонетками Выводы. нужно рассчитать следующие параметры. Успешно решена задача расчета повышения энергетической эффективности тепломассообмен- Вместимость одного туннеля по высушенному ных установок путем применения ТТ, это доказано в материалу (кг): работе путем применения математической модели в процессе конвективной сушки, позволяющей утили- GT G2 r (23) зировать вторичные энергоресурсы, применять не- N традиционные возобновляемые источники энергии. Математическая модель позволит разработать где G2 – производительность сушилки по высушен- тепломассообменные установки с применением ТТ одноцелевого и комплексного назначения, особенно- ному материалу, кг/ч; r – продолжительность стью которых является максимальная утилизация вторичных энергоресурсов, использование нетради- сушки, ч; N – количество туннелей (обычно кратно ционных возобновляемых источников энергии, при- двум). менение ТТ для предварительного нагревания про- дукта. Кроме того, ТТ создает возможность осу- Число вагонеток, одновременно находящихся в ществлять в одной и той же тепломассообменной туннеле: установке нагревание и охлаждение различных тех- нологических потоков, что особенно важно для со- n Gr G2 r (24) временных фермерских хозяйств, которые являются Lвaг 24 N Lвaг потребителями большого количества тепла при пере- работке сырья. gвaг – масса вагонетки, кг. Для широкого внедрения этих установок обосно- вана целесообразность их использования, что позво- Полученную по расчету величину n округляют лит провести научную разработку энергосберегаю- до целого числа вагонеток. Размеры туннеля опреде- щих установок при переработке некоторых пищевых ляются габаритами вагонеток и числом вагонеток, и сельскохозяйственных продуктов, методы расчета находящихся в туннеле. и проектирования установок. Список литературы: 1. Гинзбург А.С. Расчет и проектирование сушильных установок пищевой промышленности. – М.: Агропром- издат, 1985. – 336 с. 2. Дан П.Д., Рей Д.А. Тепловые трубы / Пер. с англ. – М.: Энергия, 1979. – 272 с., ил. 3. Кавецкий Г.Д., Королев А.В. Процессы и аппараты пищевых производств. – М.: Агропромиздат, 1991. – 432 с. 4. Касаткин А.Г. Основные процессы и аппараты химической технологии. – М.: Химия, 1973. – 754 с. 5. Кувшинский М.Н., Соболева А.П. Курсовое проектирование по предмету «Процессы и аппараты химической промышленности». – М.: Высшая школа, 1980. – 223 с. 6. Процессы и аппараты пищевых производств / В.Н. Стабников, В.Д. Попов, Ф.А. Редько и др. – М.: Пищевая промышленность, 1976. – 663 с. 7. Саламатов Ю.П. Идеализация технических систем [Электронный ресурс] – Режим доступа: http://www.triz- minsk.org/e/21102100.htm (дата обращения: 22.05.2013). 8. Хикматов Д.Н. Совершенствование процесса комбинированной сушки абрикоса: Автореф. дисс. на соиск. уч. ст. канд. техн. наук. – Ташкент, 2011. – 118 с. 9. Чи С. Тепловые трубы: Теория и практика / Пер. с англ. – М.: Машиностроение, 1981. – 207 с. 10. Энергоресурсосберегающие технологии при переработке плодов и овощей / А.Ф. Сафаров, К.Х. Гафуров, Д.Н. Хикматов, А.А. Холиков. – Бухара: Дурдона, 2013. – 248 с. 66
№ 10 (67) октябрь, 2019 г. ОПТИМИЗАЦИЯ ПРОЦЕССА ОТБЕЛКИ СОЕВОГО МАСЛА Хамракулова Муборак Хакимовна д-р техн. наук (PhD), Ферганский политехнический институт, Узбекистан, г. Фергана Абдуллаева Масохат Абдулбориевна ст. преп., Ферганский политехнический институт, Узбекистан, г. Фергана Хошимов Илхомжон Эркин угли ст. преп., Ферганский политехнический институт, Узбекистан, г. Фергана Турдибоев Илхомжон Хаётжон угли ассистент, Ферганский политехнический институт, Узбекистан, г. Фергана E-mail: [email protected] BLEACHING PROCESS OPTIMIZATION OF SOYBEAN OIL Muborak Hamrakulova Doctor of Technical Science, PhD, Ferghana Polytechnic Institute, Uzbekistan, Ferghana Masokhat Abdullayeva Senior Lecturer, Ferghana Polytechnic Institute, Uzbekistan, Ferghana Ilkhomjon Hoshimov Senior Lecturer, Ferghana Polytechnic Institute, Uzbekistan, Ferghana Ilkhomjon Turdiboyev Assistant, Ferghana Polytechnic Institute, Uzbekistan, Ferghana АННОТАЦИЯ В статье приводятся результаты опытов с применением наиболее эффективного способа определения опти- мальных значений технологических параметров процесса отбелки соевого масла на активированной глине «Уз- бекистан». ABSTRACT The article presents the results of experiments of the most effective way to determine the optimal values of the tech- nological parameters of the process of bleaching soybean oil on activated clay \"Uzbekistan\". Ключевые слова: глина, адсорбционная очистка, адсорбент, отбелка, нейтрализация, госсипол, критерий оптимизации, осветление, критерий Кохрена, дисперсия. Keywords: clay, adsorption treatment, adsorbent, bleaching, neutralization, gossypol, optimization criterion, clarifi- cation, Cokhren criterion, dispersion. ________________________________________________________________________________________________ Исследование глин «Навбахор», «Зарафшон» и К сожалению, работ по отбелке соевого масла «Узбекистан» при отбелке нейтрализованного сое- мало, а вопросы оптимизации данного процесса вого масла показало, что наиболее эффективным ад- практически не рассмотрены. сорбентом оказалась последняя [4, с. 17-19]. Как показали результаты предварительных ис- Ранее [10; 14] проведены исследования местных следований местных глин в процессе отбелки сое- глин в процессе адсорбционной очистки хлопкового вого масла [4, с. 17-19], для активированной глины масла, содержащего остатки госсипола, лорофилла и «Узбекистан» необходимо определить оптимальные их производных. технологические режимы, обеспечивающие макси- мальное осветление целого продукта. __________________________ Библиографическое описание: Оптимизация процесса отбелки соевого масла // Universum: Технические науки : электрон. научн. журн. Хамракулова М.Х. [и др.]. 2019. № 10(67). URL: http://7universum.com/ru/ tech/ar- chive/item/7993
№ 10 (67) октябрь, 2019 г. Планирование эксперимента [14] является наибо- температура процесса, оС; х2 – количество активи- лее эффективным способом определения оптималь- ных значений технологических параметров процесса рованной глины «Узбекистан»; х3 – время отбелки, отбелки соевого масла на активированной глине «Уз- час. бекистан». Из качественных показателей в качестве крите- При этом опыты проводятся по полному фактор- ному эксперименту (ПФА) по глине ������ = 23[6]. рия оптимизации (У) выбрана цветность соевого Из анализа процесса отбелки соевого масла вы- масла, измеряемая в мг J2. По методике [14] составлены интервалы варьи- явлены следующие переменные факторы: х1 − рования и уровни переменных факторов , ко- торые представлены в табл. 1. Таблица 1. Интервалы варьирования и уровни переменных факторов Переменные факторы Основной Интервал Уровень факторов уровень варьирования Верхний (+) Нижний (-) Температура процесса, х1, оС Количество активированной глины, х2, г 65 10 75 55 Время отбелки, х3, час 1,5 0,5 2,0 1,0 0,75 0,25 1,0 0,5 Экспериментальные исследования проводились вергали сушке под вакуумом. Опыты проводили со- на лабораторной установке для отбелки раститель- гласно рандомизированной матрице по двум парал- ных масел при интенсивном перемешивании – 60 лельным пробам [8]. об/мин. Для отбелки использовали нейтрализованное соевое масло с цветностью 55 мг/г и кислотным чис- Результаты экспериментов представлены в лом 0,29 мг кон/г. Перед отбелкой соевое масло под- табл. 2. По данным табл. 2, проведены расчеты коэффи- циентов следующего регрессионного уравнения (мо- дели): У = ������0 + ������1������1 + ������2������2 + ������3������3 + ������1,2������1������2 + ������2,3������2������3 + ������1,3������1������3 … (1) где У – функция отклика (критерий оптимизации), ������2{У} = ∑���������=��� 1(У̅������−У������)2 (3) например цветность соевого масла, мг J2; В0 и В������ – свободный член и коэффициенты регрессионного ������−1 уравнения У, которые вычисляются по формуле: где – среднее арифметическое значение параметра Вi = ∑iN=1 XijУ̅ (2) N оптимизации; n – число параллельных опытов (на нашем примере n=2). где i, j – номер опыта и фактора; N – число опытов. Используя уравнение (2), мы рассчитали коэффи- Результаты расчета значений дисперсий для уравнения У представлены в табл. 3. циенты Bi уравнения регрессии У. Результаты компьютерного расчета представ- Критерий Кохрена ������, оценивающий воспроизво- димость эксперимента, рассчитывали по следующей лены в табл. 2. формуле [3]: Дисперсию для каждого опыта определяли по ������ = ���������2���������������{У} (4) следующей формуле [12]: ∑���������=��� 1 ���������2��� {У} Таблица 2. Порядок реализации эксперимента и результаты наблюдений и расчета У для глины «Узбекистан» Номер опыта Х0 Х1 Х2 Х3 Х1Х2 Х1Х3 Х2Х3 Уl Уl У̅ср У̂выч 1 + + + - + - - 22,0 22,0 22,0 17,25 2 + - - - + + + 23,0 25,0 24,0 23,75 3 + + - - - - + 27,0 25,0 26,0 30,75 4 + - + - - + - 19,0 21,0 20,0 20,25 5 + + + + + + + 17,0 19,0 18,0 18,25 6 + - - + + - - 13,0 15,0 14,0 9,75 7 + + - + - + - 22,0 22,0 22,0 21,75 8 + - + + - - + 11,0 13,0 12,0 16,25 19,75 2,25 -1,75 -3,25 -2,5 1,25 2,5 68
№ 10 (67) октябрь, 2019 г. Таблица 3. Расчетные значения дисперсий опыта и критерия Стьюдента для уравнения У Значения уравнения У Показатели уравнения У 0 2,0 2,0 2,0 2,0 2,0 0 2,0 12,0 S{B} 1,666 B 0,15 0,0188 0,124 0,286 где ���������2���������������{У} – максимальное значение дисперсии. У = 19,75 + 2,25X1 − 1,75X2 − 3,25X3 − Значение критерия Кохрена для уравнения У 2,5X1X2 + 1,25X1X3 + 2,5X2X3 (9) представлено в табл. 3. Из табл. 3 видно, что экспериментальное значе- Анализ значений коэффициентов уравнения У показал, что по количественному вкладу исследован- ние критерия Кохрена в уравнении У оказалось ма- ные факторы ������1 ÷ ������2 располагаются в следующем лым по сравнению с табличным [3]. Это говорит о порядке убывания: том, что для дисперсии опыта можно использовать среднее значение из дисперсий, которое определя- ется по формуле [2]: ���������2���������{У} = ∑���������=��� 1 ���������2��� {У} (5) ������3 > ������2������3 > ������1������2 > ������1 > ������2 > ������1������3 ������ Далее нами изучена адекватность уравнения У Результаты расчета ���������2���������{У} представлены в табл. истинной поверхности отклика. Для этого рассчиты- ваем остаточную дисперсию (������������), критерий Фишера 3. (������) и сравниваем с их табличными значениями [7; 9]. Дисперсию коэффициентов в уравнении У опре- Остаточную дисперсию рассчитываем по следу- деляли по следующей формуле [13]: ющей формуле [9]: Отсюда: ������2{������} = ���������2���������{У} (6) ���������2��� = ∑������������=1(У̅������ − У̂������)2 (10) (7) ������ Из этого: ������{������} = √������2{������} Значения ������2{������} и для уравнения У приве- ������ = ∑���������=��� 1 ���������2��� (11) ������������ дены в табл. 3. Отсюда: (12) При этом оценка значимости коэффициентов ������������ = ������ − ������ − 1 осуществляется по критерию Стьюдента-Фишера где: Kи – число факторов и степеней свободы урав- при уровне значимости ������ = 0,95 [5]: нения У. Расчетное значение критерия Фишера определя- ∆������������ = ±������ ∙ ������{������} (8) ется по следующей формуле [7]: где – значение критерия Стьюдента при выбранной ������расч = ������ (13) доверительной вероятности или уровне значимости ������2(У������) (в нашем ������ = 2,31). Значение для уравнения У приведено в табл. 3. Полученное уравнение У имеет следующий вид: 69
№ 10 (67) октябрь, 2019 г. Показатели адекватности уравнения У истинной проводить опыты более высокого порядка, т. е. «кру- поверхности отклика представлены в табл. 4. тое восхождение» [7]. Как видно из данных табл. 4, расчетные значения Таким образом, следующие оптимальные техно- критерия Фишера меньше, чем их табличные [7]. Это говорит о том, что уравнение У отвечает истинной логические режимы: ������1 = 55℃; ������2 = 2.0 г ; Х3 = поверхности отклика. 1,0 час – могут быть использованы при отбелке сое- вых масел на активированной глине «Узбекистан». Из табл. 2 видно, что наибольшая глубина освет- ления соевого масла на активированной глине «Узбе- Опытно-производственные испытания предлага- кистан» получена в опыте 8, где цветность масла со- емых оптимальных технологических режимов от- ставляет 12,0 ед. мг J2. Поэтому нет необходимости белки соевого масла на глине «Узбекистан» подтвер- дили правильность полученных ранее результатов. Список литературы: 1. Бахтияров С.Б. Совершенствование технологии получения активированных адсорбентов из местных глин: Автор. дисс. канд. техн. наук. – Ташкент, 1998. – С. 22. 2. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. – М.: Статистика, 1973. – 392 с. 3. Иванов А.З., Круг Г.К., Филоретов Г.Ф. Статистические методы в инженерных исследованиях. Регрессион- ный анализ. – М.: Химия, 1977. – 78 с. 4. Исследование процесса отбелки нейтрализованного соевого масла с применением местных отбельных глин / М.Х. Хамракулова, Ю. Кадиров, Д.Н. Далимов, Р.М. Турсунова // Узбекский химический журнал. – 2002. – № 5. – С. 17-19. 5. Кафаров В.В. Методы кибернетики в химии и химической технологии. – М.: Химия, 1970. – 418 с. 6. Налимов В.В. Применение математической статистики при анализе веществ. – М.: Физматиз, 1960. – 27 с. 7. Налимов В.В. Теория эксперимента. – М.: Наука, 1971. – 208 с. 8. Налимов В.В., Чернова И.М. Статистические методы планирования экстремальных экспериментов». – М.: Наука, 1965. – 125 с. 9. Панфилов В.А. Научные основы развития технологических линий пищевых производств. – М.: Агропромиз- дат, 1986. – 245 с. 10. Повышение качества пищевых хлопковых масел методом их адсорбционной очистки на активированных местных глинах / С.А. Абдурахимов, С.Б. Бахтияров, З.С. Салимов и др. // Сб. науч. тр. Междунар. науч.- практ. конф., посвящ. 1200-летию со дня рождения Ахмад Аль-Фаргоний. – Фергана, 1998. – С. 21-22. 11. Подбор адсорбентов для облагораживания растительных масел и жиров / С.Б. Бахтияров, С.А. Абдурахимов, З.С. Салимов, А.Р. Сарыксакходжаев // Тез. докл. Узбекско-Малайзийского симпозиума. – Ташкент, 1997. – С.160-161 12. Рузиков Л.М., Слободгикова Р.И. Планирование эксперимента в химии и химической технологии. – М.: Хи- мия, 1980. – 280 с. 13. Рузимов Л.П. Статистические методы оптимизации химических процессов. – М.: Химия, 1972. – 199 с. 14. Хартман К., Лецкий Э., Шифер В. Планирование эксперимента в исследовании технологических процессов. – М.: Мир, 1977. – 552 с. 15. Дадаходжаев А. Т., Маматалиев Н. Н. СПОСОБЫ ИЗВЛЕЧЕНИЯ НИКЕЛЯ ИЗ ОТХОДОВ ПРОИЗВОДСТВ И ЕГО ПРИМЕНЕНИЕ //Universum: технические науки. – 2019. – №. 4 (61). 70
№ 10 (67) октябрь, 2019 г. ИССЛЕДОВАНИЕ И ОПИСАНИЕ ЭКСПЕРИМЕНТАЛЬНЫХ ЛАБОРАТОРНЫХ УСТАНОВОК Жумаев Ботир Мелибаевич соискатель, Ташкентский государственный технический университет, Узбекистан, г. Ташкент Султанова Шахноза Абдувахитовна доцент, PhD, Ташкентский государственный технический университет, Узбекистан, г. Ташкент E-mail: [email protected] Сафаров Жасур Эсиргапович д-р техн. наук, Ташкентский государственный технический университет, Узбекистан, г. Ташкент E-mail: [email protected] RESEARCH AND DESCRIPTION OF EXPERIMENTAL LABORATORY MACHINES Botir Jumayev Degree-seeking student, Tashkent State Technical University, Uzbekistan, Tashkent Shakhnoza Sultanova Associate Professor, PhD, Tashkent State Technical University, Uzbekistan, Tashkent Jasur Safarov Doctor of Technical Science, Tashkent State Technical University, Uzbekistan, Tashkent АННОТАЦИЯ В статье рассмотрено исследование, включая описание экспериментальных лабораторных установок. Разра- ботанная нами технология конвективной сушки с предварительной обработкой включает в себя процессы кон- вективного энергоподвода на стадии постоянной скорости сушки и ИК-энергоподвода на стадии падающей ско- рости сушки. Интенсификация процесса сушки способствует улучшению качества сушеного продукта и сниже- нию потерь витаминов и других ценных питательных веществ. Однако толщина кусочка может быть уменьшена до определенной величины (2 мм), так как резка на более тонкие кусочки приводит к образованию большого количества крошек. ABSTRACT The article deals with research including a description of experimental laboratory machines. Our developed convec- tive drying technology with pre-treatment includes processes of convective energy supply at the stage of constant drying speed and IR energy supply at the stage of decreasing drying speed. Intensification of the drying process improves quality of the dried product and reduces loss of vitamins and other valuable nutrients. However, the piece thickness can be reduced to a certain size (2 mm), since cutting into thinner pieces leads to the formation of a large number of crumbs. Ключевые слова: сушка, установка, тепловой процесс, качество. Keywords: drying; machine; heating process; quality. ________________________________________________________________________________________________ Сушка – тепловой процесс обезвоживания про- На качество сушеных объектов влияние оказы- дуктов путем испарения влаги и отвода образую- вают подготовительные технологические операции: щихся паров. При этом в веществе происходит пере- форма нарезки и время предварительной тепловой нос тепла и диффузионное перемещение влаги. Про- обработки. цесс сушки используется во многих отраслях произ- водства [1-2]. При подготовке к сушке инжир и яблоки режут на кусочки различных размеров и формы: столбики, __________________________ Библиографическое описание: Жумаев Б.М., Султанова Ш.А., Сафаров Ж.Э. Исследование и описание экспери- ментальных лабораторных установок // Universum: Технические науки : электрон. научн. журн. 2019. № 10(67). URL: http://7universum.com/ru/tech/archive/item/8030
№ 10 (67) октябрь, 2019 г. кружки, дольки, стружка, кубики и пластинки. При бланшировании сахара клейстеризованный Форма и размеры кусочков оказывают большое вли- крахмал и желирующие пектиновые вещества прони- яние на скорость сушки, а, следовательно, на произ- кают в межклеточное пространство и закупоривают водительность сушильной установки. С уменьше- поры, тем самым затрудняется удаление влаги при нием толщины кусочков продукта сокращается про- последующей сушке. Также в процессе предвари- должительность обезвоживания и ускоряется время тельной тепловой обработки наблюдаются ощути- восстановления сушеного продукта при его кулинар- мые потери красящих веществ, что отрицательно ска- ной обработке. зывается на товарных свойствах готовых продуктов. Если продукты нарезаны на мелкие кусочки, то Тепловая обработка позволяет снизить актив- поверхностное затвердевание происходит в меньшей ность разрушающих витаминов, ферменты и сохра- степени. нить вкус, но при этом потери витаминов могут до- стигнуть 30%. Интенсификация процесса сушки способствует улучшению качества сушеного продукта и снижению Как показывают исследования, если применять потерь витаминов и других ценных питательных ве- предварительную обработку с ИК-нагревами, то в ществ. Однако толщина кусочка может быть умень- бланшировании нет особой необходимости. Объем- шена до определенной величины (2 мм), так как резка ный и глубокий ИК-нагревы на первом этапе сушки на более тонкие кусочки приводит к образованию создают достаточно высокую температуру (70-80°С) большого количества крошек. внутри частиц продукта. Такая температура уже дает бланширующий эффект. По мере снижения скорости Предпочтение отдают сушеным фруктам, наре- сушки снижается, и температура внутри продукта и занным в виде кубиков, пластинок и стружки, по- он досушивается до характерной корочки подсыха- скольку такой продукт имеет большую насыпную ния. массу, равномерно смешивается в смесях, хорошо дозируется в мягкую упаковку на автоматах и имеет Авторами использована усовершенствованная более привлекательный вид. установка для предварительной обработки сель- хозпродуктов. Содержание мелочи в нарезанном сырье не должно превышать 5-8%. Повышенное содержание На рис. 1 представлена установка для предвари- крошки ухудшает условия сушки и приводит к из- тельной обработки фруктов в электромагнитном лишним потерям, так как при этом уменьшается вы- поле ИК-диапазона. Рабочая камера 1 выполнена в ход стандартного сушеного продукта и увеличива- виде изолированного корпуса 2 с перемещающимся с ется расход сырья. помощью винта 3 в вертикальной плоскости экраном. На экране в рабочей камере закреплён набор трубча- Неравномерная резка по ширине и толщине, тых источников ИК-излучения 4 марки КГ-220-1000. наличие слипшихся или не полностью разрезанных Установка подключается к сети посредством щита 5, частиц также недопустимы, так как нарушается пра- камера экранирована алюминиевой фольгой 6 для от- вильный режим сушки, продукт неравномерно обез- ражения ИК-лучей, установка снабжена электрон- воживается, в результате чего требуются затраты до- ным потенциометром 7 для измерения температуры полнительного труда на сортировку и досушку круп- на образце продукта 13, в камеру введена хромель- ных кусочков, выходящих из сушилки с повышенной копелевая (ХК) термопара 8 потенциометра, образец влажностью. Поверхность среза должна быть ров- размещается на поддоне (чашке) 11 весов 9, убыль ной, гладкой, при этом клетки сырья разрушаются влаги в материале регистрируется прибором 10. Рас- меньше и потери витамина С тоже небольшие. положение по высоте поддона устанавливается регу- лятором 12. В настоящее время среди ученых и специалистов в области сушки сельхозпродуктов существуют про- При экспериментировании варьирования уста- тиворечивые мнения о целесообразности примене- новка ламп от поверхности продукта осуществлялась ния перед сушкой сырья предварительной тепловой в диапазоне 100-200 мм за счёт перемещения экрана обработки. с лампами с помощью винта 3. С одной стороны, считается, что бланширование Температура тела нагрева излучателей регулиру- паром или горячей водой является необходимым ется изменением напряжения. Температура, создава- условием сохранения цвета, вкуса, запаха, витамин- емая ИК-излучением в продукте, измеряется с помо- ной активности, ускорения восстановления, главным щью хромель-копелевых термоэлектрических преоб- образом, разрушения окислительных ферментов – разователей 7 и регистрируется самопишущим элек- оксидаз и предотвращения гидролиза или окисления тронным потенциометром типа КСП-4 (класс точно- липидов в целях предупреждения потери уровня ка- сти 0,25). чества продуктов в процессе обезвоживания и осо- бенно последующего хранения. С другой стороны, Исследована кинетика сушки фруктов путём ре- многочисленные исследования показывают, что в ализации радиационно-вакуумного способа сушки в бланшировании нет особой необходимости. установке собственной разработки и изготовления, параметры которых приведены в табл.1. Отрицательным фактором тепловой обработки фруктов перед сушкой является частичное выщела- чивание из них растворимых веществ (сахаров, мине- ральных веществ, кислот и др.) и потери водораство- римых витаминов. 72
№ 10 (67) октябрь, 2019 г. Рисунок 1. Схема экспериментальной установки Обоснование выбора предварительной тепловой для предварительной обработки обработки базируется на исследованиях, проводи- мых в консервной и сушильной промышленности. 1-поддон; 2-обрабатываемый материал; Исследования и разработка терморадиационных спо- 1-рабочая камера; 3-отражатель из алюминиевой собов обезвоживания фруктов актуальны как с точки зрения повышения качества сушеных продуктов, так фольги; 4-корпус; 5- ИК-лампа; 6-щит для и энергосбережения при их производстве. подключения установки; 7-термодатчик ; 8-весы; Разработанная авторами технология конвектив- 9-регулятор расположения поддона; ной сушки с предварительной обработкой включает в себя процессы конвективного энергоподвода на стадии постоянной скорости сушки и ИК-энергопод- вода на стадии падающей скорости сушки. Исследования проводили на специально для этого созданной установке, содержащей ИК-лампы мощностью в 500 Вт (2 шт), термометр, барометр и вакуум-компрессор. Исследовательский комплекс позволяет включать отдельно предварительную об- работку с ИК большой мощности в течение 90 с, ва- рьировать параметры и режимы сушки. Определение количества сахаров осуществля- лось методом высокоэффективной жидкостной хро- матографии по ГОСТу 53883-2010, определение ас- корбиновой кислоты по ГОСТу 53693-2009, влаж- ность материала и активность воды по ГОСТ 28561- 90 в лаборатории Госстандарта Республики Узбеки- стан. Таблица 1. Параметры предварительных обработок Параметры процесса Виды предварительных обработок Мощность источников ИК СВЧ ИК+СВЧ Плотность теплового потока Продолжительность обработки 1 кВт 1 кВт 1 -1 кВт 25-30 кВт/м2 25-30 кВт/м2 25-30 кВт/м2 90 с 90 с 45 с / 45 с (45х60х45) (45х60х45) (45х60х45) Плотность поверхностного теплового потока из- мерена при помощи измерительного блока ИПП-2 по ГОСТ 25380-92. Список литературы: 1. Норкулова К.Т., Сафаров Ж.Э., Султанова Ш.А., Ахмедов Ш.И., Жумаев Б.М. Перераспределения биоактив- ных веществ в процессах сушки. // Международный сельскохозяйственный журнал. –Москва, 2016. №5. С.51- 52. 2. Сафаров Ж.Э., Парпиев З.Т., Жумаев Б.М. Сушка моркови с сохранением биологически активных веществ. // Научная конференция молодых ученых «Актуальные проблемы химии природных соединений». -Ташкент, 2015. С.194. 73
ДЛЯ ЗАМЕТОК
ДЛЯ ЗАМЕТОК
Научный журнал UNIVERSUM: ТЕХНИЧЕСКИЕ НАУКИ № 10(67) Октябрь 2019 Часть 1 Свидетельство о регистрации СМИ: ПИ № ФС 77 – 66236 от 01.07.2016 Свидетельство о регистрации СМИ: ЭЛ № ФС 77 – 54434 от 17.06.2013 Подписано в печать 25.10.19. Формат бумаги 60х84/16. Бумага офсет №1. Гарнитура Times. Печать цифровая. Усл. печ. л. 4,75. Тираж 550 экз. Издательство «МЦНО» 123098, г. Москва, улица Маршала Василевского, дом 5, корпус 1, к. 74 E-mail: [email protected] www.7universum.com Отпечатано в полном соответствии с качеством предоставленного оригинал-макета в типографии «Allprint» 630004, г. Новосибирск, Вокзальная магистраль, 3 16+
