CC BY
40
И. И. Гильмутдинов, И. М. Гильмутдинов, И. В. Кузнецова,
Р. И. Гиззатов, А. Н. Сабирзянов
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗАРОДЫШЕОБРАЗОВАНИЯ И РОСТА ЧАСТИЦ
ИБУПРОФЕНА В ПРОЦЕССЕ БЫСТРОГО РАСШИРЕНИЯ СВЕРХКРИТИЧЕСКОГО
ФЛЮИДНОГО РАСТВОРА В МИКРОННОМ КАНАЛЕ И СВОБОДНОЙ СТРУЕ
Ключевые слова: ибупрофен, зародышеобразование, рост частиц, расширение сверхкритического флюидного раствора.
Проведено описание процессов зародышеобразования и роста частиц в процессе расширения двухмерного, стационарного, осесимметричного, вязкого и сжимаемого потока «сверхкритический СО2- ибупрофен» в канале постоянного сечения и в свободной струе. Обнаружено, что поток является сверхзвуковым с наличием ударных волн и турбулентных потоков. Результаты моделирования процесса зародышеобразования и роста частиц ибупрофена совпадают с полученными экспериментальными данными.
Keywords: ibuprofen, nucleation, particle growth, the expansion of supercritical fluid solution.
The description of the processes of nucleation and growth ofparticles in the process of expanding the two-dimensional, stationary, axially symmetric, and the viscous compressible flow "supercritical CO2 ibuprofen" in a channel of constant cross section and in a free jet. Found that the supersonic flow is the presence of shock waves and turbulence. The simulation results of the process of nucleation and growth ofparticles of ibuprofen coincide with the experimental data.
Введение
При быстром расширении сверхкритического раствора достигаются большие перенасыщения растворенного вещества, что способствует образованию субмикронных и наночастиц с узким распределением по размерам. Особенности процессов заро-дышеобразования, конденсации и коагуляции во флюидной среде обусловлены уникальным набором теплофизических свойств сверхкритического флюида (например, плотность флюида как у жидкости, а сжимаемость больше чем у идеального газа, кинематическая вязкость сопоставима с соответствующим параметром идеального газа). Кроме того, в околокритической области растворителя физические параметры очень чувствительны к малым изменениям температуры и давления. Отдельно друг от друга эти явления достаточно изучены, но их совокупное влияние на процесс RESS делает эту проблему особенно интересной.
Для проведения энерготехнологической оптимизации, автоматизации RESS процесса и управления размером и морфологией получаемых частиц требуется полная математическая модель гидродинамики, зародышеобразования, роста частиц в канале и свободной струе. Адекватное математическое описание требует введения эмпирических параметров, получаемых из экспериментальных исследований. Для определения величины зародышеобразования и роста частиц необходимо определение профилей температуры, плотности, давления и скорости в направлении движения потока расширяющегося сверхкритического раствора. Отсутствие полного математического описания обуславливает необходимость применения численных подходов в моделировании процесса. В настоящей работе проведено описание зародышеобразования и роста частиц ибупрофена в процессе быстрого расширения диоксида углерода в канале постоянного сечения и свободной струе с использованием пакета программ Ansys Fluent.
Математическая модель и результаты
Образование и рост частиц в пределах устройства расширения и в свободной струе в результате перенасыщения сверхкритического раствора происходит за счет двух явлений: образования критических зародышей, способных к дальнейшему росту и конденсации одиночных молекул на поверхности критических зародышей и на поверхности растущих частиц. Предполагается, что зародышеобразование и конденсация частиц сферической формы протекает равномерно в каждой точке расширяющегося потока «сверхкритический С02 - ибупрофен».
Для определения скорости образования критических зародышей используем уравнение [1]:
°к Г3
кТ
I = 2N,
РУ2
exp
Ют
" 3
1
LnS-КУ?
(1)
где т2 - молекулярная масса субстанции; у2 - молекулярный объем субстанции в твердой фазе (у23=т2/р3Ь; Ь - число Авогадро; р3 - плотность субстанции в конденсированном состоянии); N -концентрация растворенного во флюидной фазе метилпарабена, ибупрофена; у2 - фактическая мольная доля растворенного метилпарабена, ибупрофена во флюидной фазе; у2еч - равновесная мольная доля растворенного метилпарабена, ибупрофена во флюидной фазе; 8 - величина перенасыщения
(5 = у 2/ У т ); К - коэффициент кристаллизации; ст - поверхностное натяжение на границе флюид -фармацевтическая субстанция; к - константа
Больцмана.
Конденсация одиночной молекулы на §- мерной сферической частице определяется уравнением [1]:
' лк____________>л
1+1.7Жи+1.33Жи ) (2)
F(g) = (4&г2 gvS ) - N? (g)f
где g - количество молекул в частице, К2еч^) - равновесная концентрация на поверхности сферической g размерной частицы, Л - корреляционный параметр, характеризующий кинетику роста частиц.
Коэффициент диффузии Б растворенного вещества в сверхкритическом растворителе находится по уравнению:
Б = 7,4 х10
ТМ
1/2
(3)
где Б измеряется в м2 сек-1, М (г моль-1) - молекулярный вес растворителя; ц (кг м-1 сек-1) - вязкость чистого растворителя; V (см-3 моль-1) - молярный объем твердого растворенного вещества.
Число Кнудсена находится по уравнению:
(4)
где X - средний путь свободного пробега молекул,
который определяется по уравнению:
1
^ 2
Х = — р
ц( 7ЇЇПХ 2кТ
(5)
Vі'1-' /
где р - плотность растворителя; т1 - молекулярная масса растворителя.
Размер критических ядер g рассчитывается по уравнению [1]:
_ .
-|2
д =-
32п
3
4^ Г3 " 1 "
кТ !пБ -
(6)
Радиус критических ядер находится по уравнению [1]:
г* = 2
о^2
~кТ
1
!пв - КУ?
(7)
При расчете динамики роста частиц путь поток сверхкритического раствора разбивается на элементарные объемы, в пределах которых остаются постоянными температура, давление, плотность и скорость. В каждом элементарном объеме в связи с перенасыщением одновременно происходит образование новых частиц и рост существующих частиц за счет конденсации. Переход к следующему элементарному объему приводит к образованию новой фракции частиц одинаковых размеров. Если, 1 - порядковый номер фракции частиц, ] - порядковый номер элементарных объемов, то уравнение динамики образования новых частиц запишется следующим образом [2]:
N
і(і+1)
N' Т. рі > г«? ‘
(8)
д,(,+1) = д * (Т - р) , (9)
где, N - количество частиц, Т - температура в объёме _|, PJ - давление в объёме Т*ес - время прохождения элементарного объёма ], и **ес - >й объем.
Рост частиц в пределах элементарного объёма определяется количеством молекул, осаждаемых за счет конденсации:
д, = е д у-1),) (10)
Из уравнений (6) и (7) можно получить зависимость диаметра частицы от количества молекул:
3ду 2
4п
(11)
Из последнего элементарного объёма, когда) принимает максимальное значение, рассчитывается итоговое распределение по размерам частиц:
о 1 = о, (д,) и N = к]г (12)
Средний размер частиц рассчитывается следующим образом:
I І____
п
(13)
Для моделирования зародышеобразования и роста частиц во всем интервале расширяющегося потока необходимо знать поля температуры, давления, плотности и скорости. Для этого поток сверх-критического раствора в расширительном устройстве постоянного сечения (рис. 1) рассматривается в настоящей работе двухмерным, стационарным, осесимметричным, вязким и сжимаемым. Расчёт ведется в области капиллярного канала (1-2) и свободной струи (2-3).
Рис. 1 - Упрощенная схема расширительного устройства: 1-2 - капиллярный канал, 2-3 - свободная струя
Гидродинамика процесса расширения из резервуара бесконечного объёма через микроканал с постоянным сечением с учетом трения описывается в виде системы дифференциальных уравнений сохранения: массы, импульса, энергии и уравнения состояния. В гидродинамических расчётах использовалось уравнение состояния для чистого С02, так как раствор «сверхкритический СО2-ибупрофен» является разбавленным.
Для двухмерной (2Б) осесимметричной геометрии потока уравнение сохранения массы в цилиндрических системах координат имеет вид [2]:
-0 р (рЭг ) + р- = 0 (14)
ОХ ОГ г
где Р - плотность, х - осевая координата, г - радиальная координата, Зх - осевая скорость и Эг - радиальная скорость.
Уравнение сохранения импульса описывается в общем виде [2]:
О - - - - -
— (р9 + У- (рЗЗ) = -Ур + р + ¥ (15)
Ог
15
где р - плотность, р - давление, Т - тензор напряжений и ~g и ~ Б являются силой тяжести тела и внешней силой, зависящий от условий задачи.
Для 2Б осесимметричной геометрии потока, осевые и радиальные уравнения сохранения импульса определяются в виде [2]:
- [грЗхЭх ) + 1 {трЗгЗх ) =
г дх г дг
др + 1 д Ox r Ox
гц
1 _o
r dr
1 d
r
1 d
гц
( дЗ
+
L дЗх 2 2—x dx
дЗг I
M
,Y
+
(1б)
dr dx
1 d^p33 )+10 (рза )=-dx +1 Ox гЦц
г dx г O dx г dx V
г 0г
где
, с3> 2 -1
гЦЧ -зН
-2ц-2 +3 - (v3)-+p-^:
3 г г
3 д3Х 03Г 3Г V-3 =—- + —- + — дХ дГ Г
(18)
Зх - осевая скорость, Зг - радиальная скорость, ^ -
молекулярная вязкость.
Закон сохранения энергии в цилиндрических системах координат описывается в виде:
д
Ot
(pE) + V-( 3(pE + p) | = -V• IhjJ
(
\
(19)
V ] / ;
В нашем случае уравнение энергии решается следующей форме:
(рЕ + •( • Т + (Т-V}) (20)
где
Е = h - * + 3
p
2
где к, у - молекулярные коэффициенты теплопроводности, вязкости, Е - полная энергия, Т - температура. Слагаемые в правой части уравнения (20) выражают передачу энергии теплопроводностью и диссипацией энергии за счёт трения внутри потока. Течение в канале Ь/0<10 в настоящей работе рассматривается как ламинарное, поэтому в уравнениях (15)-(20) используются молекулярные коэффициенты переноса. Ввиду пренебрежительно малого значения теплового потока на границе стенка-среда по сравнению с энтальпией расширяющегося сверх-критического раствора процесс истечения сверхзвукового потока принимается адиабатным.
Для описания термодинамических свойств чистого диоксида углерода в гидродинамической модели было выбрано уравнение состояния идеального газа:
Вязкость диоксида углерода рассчитывается по формулам [3]:
П = Пd * exp(In n*),
где Пid -вязкость идеального газа,
4 ( 1 a
In n*=i ar
V j-0
wJ
(22)
(23)
где \siij} - коэффициенты, представленные в работе
[156]. Здесь ю=р/рс - приведённая плотность.
Теплопроводность находится по следующей зависимости [3]:
X = Хш * ехр(1п X)
где Х^- теплопроводность идеального газа
(24)
W3 0311 + In Л*=Х ( 2 a 1 (і v I
^ ox дг J (17) /=1 Vj=o T J
w1
(25)
дель к-є второго порядка.
Граничная температура на входе в канал (Тх=0) приравнивается к температуре предрасширения, начальное давление на входе в канал (Рх=0) равняется давлению в системе до начала процесса расширения. На стенке канала принимается их=0,
8Т„те„,_, _ о Термодинамические условия на выхо-
дх
де из камеры расширения равны нормальным физическим условиям. Геометрические условия определяются соотношениями L/D. Уравнения сохранения импульса, массы и энергии решены численно методом конечных объемов. В начальном приближении задается температура по всей расчетной области равной 320 К, давление 250 бар, а плотность 50 кг/м3.
Построение расчетной области выполнялось в графическом приложении программы Fluent - Gambit. В нашем случае для двухмерной (2D) осесимметричной геометрии потока была построена сетка для канала с отношением L/D=1000/300 мкм, а также для области свободной струи для камеры расширения L/D=100/40 мм. Входные граничные условия капиллярного канала (область 1 в рис.2) определяются в соответствии с экспериментальными данными и равны: T=343 К, P=250 бар. Результаты моделирования поля скорости приведен на рис. 2.
Pv=RT
(21)
Рис. 2 - Профиль осевой скорости вдоль сопла и свободной струи
в
б5
Из полученного поля можно выяснить наличие ударной волны поперек потока на расстоянии 2 мм от выхода из капиллярного канала, как отмечается в теоретических и экспериментальных работах посвященных ИЕ88 процессу. Ударная волна влечет за собой скачкообразное изменение термодинамических параметров и скорости поток. Как видно из рис. 4.6, после ударной волны осевая скорость в определенной области приобретает отрицательные значения, можно предположить о сильно турбулентности в этой области. Скорость потока до ударной волны достигает значений в 700 м/с, что является нормой для подобного сопла и данных термодинамических параметров, а расчетное значение числа Маха достигает до 6, что свидетельствует о сверхзвуковом характере потока.
Профили температуры, давления, плотности потока расширяющегося сверхкритического флюида в канале и свободной струе позволяют оценивать значение равновесной концентрации растворенного вещества в сверхкритическом флюиде по уравнению состояния Пенга-Робинсона [4]. Так как уравнения (1)-(13) являются функциями от температуры, давления, мольной доли растворённого компонента и перенасыщения, то можно рассчитать зародыше-образование и рост частиц вдоль канала и свободной струи в соответствии с профилями термодинамических параметров состояния.
В результате моделирования процесса заро-дышеобразования и роста частиц получаем зависимости показателя перенасыщения, показателя нук-леации, радиуса критического зародыша, среднего размера частиц вдоль пути расширения по осевой координате (рис. 3).
чкхяяясмяяюм*
Рис. 3 - Профиль среднего размера частиц по осевой координаты. Линия - расчет; ■ - эксперимент
Наилучшее описание среднего размера получено при значении корреляционного параметра в уравнении (4.28) Л=300. По полученным профилям можно отметить скачкообразное изменение ключевых параметров зародышеобразования и роста частиц в окрестности ударной волны. Так же из рис. 3 следует, что частицы в пределах канала имеют наноразмеры, и достигают микронных размеров после
выхода из капиллярного канала. Увеличение среднего размера частиц в эксперименте относительно расчетных значений можно объяснить учетом в расчетах образовавшихся в свободной струе нанораз-мерных зародышей, которые в экспериментальном исследовании, в связи с особенностями методики, не учитывались. На рис. 4 представлено распределение частиц по размерам.
40 -
и —* * ■■ ■ ■ г Щ ■ ■ ■
2,5 3,0 3,5 4,0 4,5 5,0 5,5 6,0
Средний размер, мкм
Рис. 4 - Сравнение экспериментального и расчётного распределения частиц по размерам ибу-профена
Как видно из результатов, в расчетах зафиксированы, преимущественно, частицы более мелкого размера, а в эксперименте более крупного, хотя пики обеих распределений находятся в одной области. Это можно объяснить, так же агломерацией частиц в эксперименте после осаждения и механическим уносом более мелких частиц. В целом математическая модель адекватно описывает процессы за-родышеобразования и роста частиц в канале и в свободной струе.
Благодарность
Работа выполнена при поддержки Федеральной целевой программы «Научные и научнопедагогические кадры инновационной России» на 2009-2013 годы Государственный контракт № 14.B37.21.0944 от 5 сентября 2012 года.
Литература
1. Debenedetti P.G. Homogeneous nucleation in supercritical fluids. ALChE J. - 1990. - № 36. - 1289.
2. Гидродинамика и зародышеобразование в канале и свободной струе в процессе быстрого расширения сверх-критического раствора / И.В. Кузнецова, И.И. Гильмутдинов, И.М. Гильмутдинов, А.А.Мухамадеев, А.Н. Са-бирзянов // Вестник Казанского технологического университета. - Казань. - 2012. - №1. - С.111-118.
3. Алтунин В.В. Теплофизические свойства двуокиси углерода/ В.В. Алтунин. - М., Издательство стандар-тов,1975. - С.546.
4. Описание растворимости ибупрофена с использованием уравнения состояния Пенга-Робинсона / И.В.Кузнецова, И.И. Гильмутдинов, И.М. Гильмутдинов, А. А. Мухама-диев, А.Н. Сабирзянов// Вестник Казанского технологического университета. - Казань. - 2011.- №19.- С.7-11.
© И. И. Гильмутдинов - асп. каф. теоретических основ теплотехники КНИТУ, i1nur1988@inbox.ru; И. М. Гильмутдинов -к.т.н., асс. той же кафедры, gi1mutdinov@kstu.ru; И. В. Кузнецова - асс. той же кафедры, Irina301086@ramb1er.ru; Р. И. Гизза-тов - асп. той же кафедры, rami1.bua@mai1.ru; А. Н. Сабирзянов - д.т.н., проф. той же кафедры, sabirz@kstu.ru.
