ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Публикационная активность

(сведения по итогам 2017 г.)
2-летний импакт-фактор РИНЦ: 0,500
2-летний импакт-фактор РИНЦ без самоцитирования: 0,405
Двухлетний импакт-фактор РИНЦ с учетом цитирования из всех
источников: 0,817
5-летний импакт-фактор РИНЦ: 0,319
5-летний импакт-фактор РИНЦ без самоцитирования: 0,264
Суммарное число цитирований журнала в РИНЦ: 6012
Пятилетний индекс Херфиндаля по цитирующим журналам: 404
Индекс Херфиндаля по организациям авторов: 338
Десятилетний индекс Хирша: 17
Место в общем рейтинге SCIENCE INDEX за 2017 год: 527
Место в рейтинге SCIENCE INDEX за 2017 год по тематике "Автоматика. Вычислительная техника": 16

Больше данных по публикационной активности нашего журнале за 2008-2017 гг. на сайте РИНЦ

Вход


Забыли пароль? / Регистрация

Добавить в закладки

Следующий номер на сайте

4
Ожидается:
16 Декабря 2018

Моделирование вакуумной сублимационной сушки

Modeling of vacuum freeze-drying
Статья опубликована в выпуске журнала № 4 за 2011 год. [ на стр. 170 – 173 ][ 13.12.2011 ]
Аннотация:Описывается программный продукт, позволяющий рассчитывать процесс вакуумной сублимационной сушки и подбирать параметры его ведения. Используя разработанную программу, пользователь может рассчитать время про-ведения процесса, изменяя температурные режимы. Программа создана на основе математической модели, описывающей процессы сублимационной сушки, протекающие в локальном объеме, и смещение границы раздела фаз.
Abstract:The paper describes a software product which enables calculation of the vacuum freeze-drying process and to select the parameters of the process. Using the developed program, the user can calculate the time of the process of changing temperature regimes. The program is based on a mathematical model describing the processes occurring in the local volume, the displacement of the interface.
Авторы: Диденко А.А. (didenko_al@bk.ru) - Российский химико-технологический университет им. Д.И. Менделеева, г. Москва, , , Ершова А.Н. (didenko_al@bk.ru) - Российский химико-технологический университет им. Д.И. Менделеева, г. Москва, , , Гордиенко М.Г. (chemcom@muctr.ru) - Международный учебно-научный центр Российского химико-технологического университета им. Д.И. Менделеева, Москва, Россия, кандидат технических наук, Меньшутина Н.В. (chemcom@muctr.ru) - Российский химико-технологический университет им. Д.И. Менделеева, г. Москва, Россия
Ключевые слова: кинетика сушки, температурный профиль, температурный режим, вакуумная сублимационная сушка
Keywords: kinetics of drying, temperature profile, temperature, vacuum freeze-drying
Количество просмотров: 6490
Версия для печати
Выпуск в формате PDF (5.83Мб)
Скачать обложку в формате PDF (1.28Мб)

Размер шрифта:       Шрифт:

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

Перед разработчиками стояла задача создать ПО, позволяющее смоделировать процесс ВСС и описать его в единичном объеме. Для ее решения была разработана модель, рассматривающая влияние температурных режимов на ход сублимации. Актуальность задачи обусловлена тем, что режим теплоподвода является важным фактором ведения процесса: повышение температуры позволяет интенсифицировать процесс и сократить время сушки, однако превышение критического значения может привести к расплавлению материала.

Описание модели

На рисунке 1 представлена схема процесса сублимационной сушки, на которой показан так называемый флакон с двухфазовой субстанцией – льдом и высушенным продуктом. Между слоями находится граница раздела двух фаз, смещающаяся в ходе сушки вниз. Теплота к материалу подводится двумя способами: тепловым излучением от верхней полки и кондуктивным теплоподводом от нижней.

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

1)     все потоки тепла и массы в системе имеют одно измерение и направлены вдоль нормали к поверхностям;

2)     сушка вымораживанием протекает под воздействием теплового излучения от полки сверху и кондуктивного теплоподвода от полки снизу;

3)     во время сушки замороженный продукт условно можно разделить на две части: верхнюю – слой высушенного материала (R) и нижнюю – слой льда (С); граница между слоями в процессе сушки движется сверху вниз;

4)     влагосодержание замороженного продукта принимается равным начальному, влагосодержание высушенного продукта – эквивалентному;

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

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

Модель состоит из нестационарных балансовых энергетических уравнений и соответствующих граничных и начальных условий, записанных для областей замороженного слоя С и высушенного R [2].

Дифференциальные уравнения, начальные и граничные условия могут быть записаны в следующем виде.

1.     Для слоя С (замороженный слой – лед).

Подпись:  

Рис. 1. Схема процесса
Уравнение теплопроводности для замороженного материала:

                                                   (1)

в диапазоне 0£x0; X(t) – подвижная граница; граничные условия на дне флакона:

TC=T0, где x=0, t>0.                                         (2)

Уравнение изменения температуры пограничного слоя на границе с сухим слоем за счет испарения:

,                 (3)

где x=XC(t), t>0;

начальные условия: TC=TS,                          (4)

0£x

Здесь λ – теплопроводность; T0 – температура полки; TS – равновесная температура сублимации; TH – температура излучающей поверхности; TL – температура поверхности материала, нагреваемой излучением; L – толщина слоя материала; x – координата; используемые в модели индексы: ef – замороженный слой, ed – высушенный слой, eq – равновесный слой, i – исходный (замороженный) слой; TC – температура в слое замороженного материала; α – коэффициент температуропроводности; t – время; W – среднее влагосодержание по исходному материалу; ρbu – насыпная плотность высушенного материала; ∆HS – равновесная теплота сублимации; XС – положение границы между замороженным и высушенным слоями материала при наличии нагрева теплопроводностью.

2.     Для слоя R (высушенный слой).

Уравнение теплопроводности для высушенного материала:

                                                    (5)

в диапазоне XC(t)£x£L, t>0;

граничные условия: TR=TL,                                (6)

x=L, t>0.

Уравнение изменения температуры пограничного слоя сухого материала на границе с замороженным слоем:

;            (7)

x=XC, t>0;

начальные условия:

TR=TS, 0£x£L, t=0.                                                (8)

Здесь TR – температура в слое высушенного материала; XR – положение границы между замороженным и высушенным слоями материала при наличии нагрева излучением.

На границе раздела фаз должны выполняться равенства температур ТR=ТC, а также потоков тепла и массы.

Для нахождения температурного профиля воспользуемся разностной схемой Кранка–Николсо­на. Она является абсолютно устойчивой, поэтому шаг по времени может быть выбран достаточно большой для ускорения расчета, но это может привести к потере участка высушенного материала за период Δt. Чтобы избежать неточности, необходимо учитывать скорость удаления влаги и в соответствии с этим вести перерасчет шага по времени, уменьшая или увеличивая его в процессе расчета [3].

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

          (9)

Уравнение для изменения влагосодержания в материале:

     (10)

Алгоритм расчета

Для решения уравнений (1)–(10) был разработан алгоритм, графически представленный на рисунке 2. Он позволяет рассчитать температурные профили и изменение влагосодержания в материале. Порядок действий алгоритма следующий.

1.     Задаются начальные и граничные условия.

2.     Оценивается величина шага в замороженном (С) и высушенном (R) слоях.

3.     Одновременно рассчитываются температурные профили в замороженном (С) и высушенном (R) слоях.

4.     Сравниваются конечные значения температуры на границе раздела фаз, и при необходимости уточняется шаг по времени и координате.

5.     Полученное значение фактической температуры сравнивается с критической температурой.

6.     Рассчитывается изменение влагосодержания образца и нового положения границы раздела фаз.

7.     Фактическое значение влажности сравнивается с требуемым значением остаточного влагосодержания.

8.     Делается вывод о завершении процесса сушки.

ПО и результаты расчета

Для описания процесса ВСС разработано ПО на С++ с применением блочного подхода. В состав ПО входят блок ввода начальных условий и ограничений, блок математического моделирования и блок отображения результатов.

Разработанная программа позволяет рассчитать все возможные температурные режимы проведения процесса, варьировать массу загрузки. На основании полученных данных могут быть выбраны параметры проведения процесса.

Результат работы ПО – графические зависимости температурного профиля по высоте флакона, а также кинетика вакуумной сублимационной сушки (рис. 3). Программа рассчитывает три стадии ВСС: сублимационную сушку в изотермичных условиях, неизотермичную сушку и тепловую изотермичную сушку.

Графики сублимационной сушки показывают характер изменения температуры в процессе ВСС. На первой (рис. 3а) и третьей (рис. 3в) стадиях сушки температура не меняется, следовательно, во всем объеме флакона она остается неизменной.

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

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

Литература

1.     Семенов Г.В. Вакуумная сублимационная сушка. Основные понятия и определения // Сублимационная сушка в фармацевтической и пищевой промышленности: матер. Междунар. науч.-технич. конф. М.: МГУПБ, 2005. С. 92–98.

2.     Nastaj J.F. Some aspects of freeze drying of dairy biomaterials // Drying technology. 1996. № 14 (9), pp. 1967–2002.

3.     Кольцова Э.М., Гордеев Л.С. Методы синергетики в химии и химической технологии. М.: Химия, 1999. 256 с.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=2941
Версия для печати
Выпуск в формате PDF (5.83Мб)
Скачать обложку в формате PDF (1.28Мб)
Статья опубликована в выпуске журнала № 4 за 2011 год. [ на стр. 170 – 173 ]

Назад, к списку статей

Хотите оценить статью или опубликовать комментарий к ней - зарегистрируйтесь