УДК 621.594
Анализ процесса кристаллообразования диоксида углерода в проточной части центростремительного турбодетандера
Канд. техн. наук Данилов М.М. [email protected]
Университет ИТМО 191002, Россия, Санкт-Петербург, ул. Ломоносова, 9
Разработанная автором математическая модель вымораживания диоксида углерода в низкотемпературных турбодетандерах позволяет проанализировать основные стадии процесса кристаллообразования. Сопоставляется скорость образования кристаллов-зародышей и растущая удельная масса кристаллов в единице объема расширяющегося газового потока. Прослеживается изменение размеров образующихся и растущих кристаллов, а также среднего размера получаемых кристаллов. Указывается на взаимосвязь роста размеров и удельной массы кристаллов. Приводится сравнение изменений удельной массы кристаллов и размеров основной массы кристаллов при различных начальных параметрах расширяемой газовой смеси: давления и концентрации паров диоксида углерода. Важным результатом проведенного анализа стал вывод о том, что основной рост кристаллов происходит за зоной Вильсона.
Ключевые слова: кристаллы, удельная масса, размер кристалла, диоксид углерода, кристаллообразование, переохлаждение потока.
The analysis of process of grain formation of carbon dioxide in flowing a chastitsentrostremitelny turbodetander
Ph.D. Danilov M.M. [email protected]
ITMO University 191002, Russia, St. Petersburg, Lomonosov St., 9
The mathematical model developed by the author of freezing carbon dioxide in the low-temperature expansion turbines allows to analyze the main stages of the process of crystal formation. Matches rate of crystal nucleation and crystal growing specific weight per unit volume of the expanding gas stream. Track changes the image size and the growing crystals, as well as the average size of the resulting crystals. It indicates the relationship size growth and specific mass of crystals. The comparison of changes in the specific gravity of the crystals and the size of the bulk crystals with different initial parameters expandable gas mixture: pressure and vapor concentration of carbon dioxide. An important result of the analysis was the conclusion that the main crystal growth occurs in the Wilson area.
Keywords: crystals, specific gravity, the size of the crystal, carbon dioxide, crystal formation, hypothermia flow.
Получение твердого диоксида углерода (сухого льда) в объеме расширяющегося в турбо-детандере газовой холодильной машины потока газовой смеси экспериментально исследовано и предложена математическая модель, которая позволяет отслеживать изменение вдоль проточной части турбодетандера скорости образования кристаллов-зародышей диоксида углерода, массовой концентрации образующихся и растущих кристаллов, а также их размеров [1].
Разбивка проточной части центростремительного турбодетандера на участки, а также изменение степени переохлаждения потока вдоль проточной части, место расположения и протяженность зоны кристаллообразования рассмотрены в работе [2]. Изменение скорости образования кристаллов-
зародышей и массовой концентрации образующихся в каждом из последовательных сечений (участков) и растущих кристаллов диоксида углерода можно проследить на примере расширения газовой смеси с начальным давлениемРсм,0 = 2,0 • 105 Па и начальной концентрацией паров диоксида углерода гсо2,о=0,05 (рис. 1), при котором протяженность процесса кристаллообразования будет наибольшей. Скорость образования кристаллов-зародышей в большей степени зависит от размеров образующихся кристаллов.
В работе [3] определен возможный диапазон этих размеров (параметров кристалла а) при их дискретном увеличении от а = 2ам до а = 5ам (ам - параметр молекулярного кристалла). Для процесса кристаллообразования при вымораживании диоксида углерода в низкотемпературных турбодетандерах характерно образование кристаллов-зародышей с параметром кристалла а = 2ам и а = 3ам [4], причем с увеличением степени переохлаждения потока вдоль проточной части турбодетандера [2] вначале образуются более крупные кристаллы (а = 3ам), количество, а следовательно и удельная масса которых невелики, а затем более мелкие (а = 2ам), количество и масса которых возрастает на несколько порядков.
Расположение стартовых точек роста массовой концентрации, т.е. удельной массы кристаллов-зародышей, образовавшихся на каждом из последующих участков, начиная с участка начала кристаллообразования (/ = х), адекватно характеру изменения скорости образования центров кристаллизации, которая, в пределах одного и того же параметра а, сначала уменьшается вследствие понижения температуры потока при его расширении, а затем возрастает, следуя повышению температуры вследствие подвода к потоку теплоты фазового перехода (рис.1).
и 3
Рис. 1. Изменение скорости образования кристаллов-зародышей/, 1/(с^м )
и удельной массы кристаллов g вдоль проточной части центростремительного турбодетандера
19
(/-номер участка) при начальных параметрах Рсм,0 = 2,0 • 105 Па; гсо _о=0,05
Изменение характера кривых (/—номер участка, на котором образовались кристаллы-зародыши) для] = 6 на г = 13, для] = 7 на г = 14 и для] = 8 на г =Ьу связано с переходом приращения удельной массы кристаллов с микроскопического на макроскопический уровень, когда размер кристаллов становится больше длины свободного пробега молекул газовой смеси [5]. Приращение удельной массы кристаллов заканчивается на участке Ъу, когда переохлаждение потока АТ полностью снимается за счет подвода к потоку теплоты фазового перехода. Так как участки проточной части имеют определенную длину [2], то поток может достигнуть состояния насыщения (ДГ=0) еще до окончания этого участка (г=у), поэтому вводится коэффициент Ь, который уменьшает длину участка (Ь = 0,1.. .0,9). В рассматриваемом случае г = 15, а Ь = 0,2, т.е. процесс кристаллообразования заканчивается на двадцати процентах длины пятнадцатого участка. Основное влияние на величину суммарной получаемой массы кристаллов (удельной массы gi на участке Ьу) оказывает масса кристаллов, образовавшихся на участке г = 8, т.е. на том участке, где впервые начали появляться кристаллы-зародыши с параметром а = 2ам, чему способствует достижение потоком второго критического переохлаждения.
Аналогичный вывод можно сделать, если на том же примере проследить за изменением размера растущих кристаллов (рис. 2).
Рис. 2. Изменение параметра кристалла а вдоль проточной части центростремительного турбодетандера (/-номер участка) при начальных параметрах Pсм,0 = 2,0 • 105 Па; г со ^=0,05
Средний размер получаемых кристаллов а будет наиболее близок к размеру кристаллов, образовавшихся на участке г = 8, т.к. на этом участке образуется наибольшее количество кристаллов (около 60 %). Взаимосвязь роста размеров кристаллов и роста их удельной массы отражена в математической модели анализируемого процесса [1]. Наибольшего размера достигают кристаллы, образо-
вавшиеся в самом начале процесса кристаллизации (при достижении потоком первого критического переохлаждения), но количество таких кристаллов ничтожно мало (на 13 порядков меньше, чем кристаллов, образовавшихся на участке г = 8).
Теперь рассмотрим процесс кристаллообразования при расширении газовых смесей с различными начальными концентрациями паров диоксида углерода и различными начальными давлениями газовой смеси [2].
На характер изменения удельной массы всех образующихся кристаллов по мере их роста (рис. 3) основное влияние оказывает стартовая масса кристаллов, т.е. масса зародышей, образующихся при втором критическом переохлаждении потока, которая пропорциональна скорости их образования, и которая будет тем больше, чем выше концентрация паров диоксида углерода в газовой смеси [6].
Рис. 3. Изменение суммарной удельной массы ^вдоль проточной части центростремительного турбодетандера (/-номер участка) при различных начальных параметрах газовой смеси Рсм,0; гС02,о (заштрихованный участок - зона Вильсона)
Именно эта тенденция оказывает решающее влияние на величину максимального переохлаждения потока АТтох и его месторасположение [2], т.к. резкое повышение температуры потока вследствие подвода к нему теплоты фазового перехода, вызывающее снижение степени переохлаждения, зависит от величины массы растущих кристаллов. Наступление этого момента (вход потока в зону Вильсона) происходит тогда, когда удельная масса кристаллов достигает величины 1.. .8
3 „
г/м (рис.3). Величина максимальной степени переохлаждения АТтах = 11,4... 19,6 К [2], причем более высокие значения этой величины соответствуют более низкой начальной концентрации паров диоксида углерода в газовой смеси. Можно сделать вывод и о том, что, если до зоны Вильсона рост массы кристаллов зависит от степени переохлаждения потока, то затем наблюдается обратная зависимость, т.е. дальнейший рост массы вызывает снижение переохлаждения (зоной Вильсона считается зона, где разность между равновесной и текущей концентрациями диоксида углерода достигает максимума).
Увеличение удельной массы кристаллов зависит и от характера роста их размера (рис. 4), т.к. приращение массы на любом из рассматриваемых участков будет зависеть от размера кристаллов на предшествующем участке [1], т.е. рост размера будет первичен по отношению к росту массы.
Рис. 4. Изменение параметра кристалла аосновной массы кристаллов вдоль проточной части центростремительного турбодетандера (/-номер участка, /-участок зарождения) при различных
начальных параметрах газовой смеси Рсм,0; тС02,о
При одинаковом параметре кристаллов-зародышей, например а = 2ам, более высокая начальная концентрация паров диоксида углерода способствует образованию большей массы кристаллов, поэтому, даже при одинаковом приросте массы на последующем участке, общая масса кристаллов на нем, а следовательно, и выросший размер кристаллов, будут также тем больше, чем выше концентрация диоксида углерода. Это приводит к тому, что и следующее приращение массы будет следовать той же тенденции. Вышеперечисленные зависимости находят свое отражение в характере наклона кривых на рис. 3 и 4.
Проведенные на рис. 4 линии Вильсона соединяют точки, соответствующие размеру кристаллов в момент достижения потоком максимального переохлаждения для газовых смесей с различной начальной концентрацией паров диоксида углерода с учетом места зарождения этих кристаллов. Размеры кристаллов в зоне Вильсона - от 0,002 • 10-6 м для высоких концентраций паров диоксида углерода до 0,032 • 10-6 м для низких концентраций.
Если сравнить эти величины с конечными размерами кристаллов, то можно сделать вывод о том, что основной рост кристаллов происходит уже за зоной Вильсона, т.е. в период резкого снижения переохлаждения потока.
Список литературы
1. Данилов М.М. Моделирование процесса вымораживания диоксида углерода в объеме расширяющегося газового потока // Известия СПбГУНиПТ. 2007. № 1. С. 6 - 8.
2. Данилов М.М. Вымораживание диоксида углерода из газовых смесей в проточной части низкотемпературных турбодетандеров // Известия СПбГУНиПТ. 2009. № 3/4. С. 10 - 13.
3. Данилов М.М., Смирнов А.С. Основные особенности образования твердой фазы диоксида углерода // Вестник Международной академии холода. 2014. № 2. С. 37 - 40.
4. Данилов М.М. Скорость образования центров кристаллизации при вымораживании диоксида углерода // Известия СПбГУНиПТ. 2006. № 1. С. 21 - 23.
5. Данилов М.М. Рост кристаллов диоксида углерода в объеме расширяющегося газового потока // Известия СПбГУНиПТ. 2006. № 1. С. 18 - 20.
6. Данилов М.М., Смирнов А.С. Анализ влияния параметров газовой смеси на скорость кристаллообразования при вымораживании диоксида углерода // Современные проблемы холодильной техники и технологий: Сб. тез.докл. ЕХМеждунар. науч.-техн. конф. - Одесса, 2013. С. 64 - 65.
Статья поступила в редакцию 25.10.2015 г.