Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

THE MODEL OF THE PROCESS OF EVAPORATION AND GROWTH OF MOVING DROPLETS

Fedorov V.N. 1 Kulishenko R.Yu. 1 Prostitenko O.V. 1 Bobtenkov A.V. 1
1 St. Petersburg State Institute of Technology
The model proposed in this article describes the process of evaporation and growth of droplets moving in a gas. Creation of a strict model of droplet growth can be based on time-dependent kinetic equations. The growth of spherical shape single-particle from the supercooled vapor or gas-saturated liquid vapor obeys the general laws of hydrodynamics and heat and mass transfer in continuous medium, which can accurately predict its growth rate. Also note that in droplet movement in medium its rate of evaporation may change, since the vapor molecules, which are close to the surface of the droplet, are carried away by the environment. For sufficiently large droplets, compared to the size of atoms, molecules and clusters, they will have a developed surface, so the solution of these problems is complicated by the presence of interfaces. Raises difficult problems of formation of the boundary conditions and solutions of boundary value problems for kinetic equation. In the model supposed by authors of this article considered the effect of all the above factors.
evaporation
hydrodynamics
gas atmosphere
aerosols
fluid mechanics
dedusting
Представим конденсационное пылеулавливание как процесс эволюции во времени большой системы капель. В принципе, вероятно, возможно рассмотреть поставленную задачу, решив уравнение Лиувилля при надлежащем выборе начальных и граничных условий. Детальный анализ такого решения должен выявить особенности, наблюдаемые при конденсации в дисперсных системах. Сказанное основано на следующей фундаментальной идее: если задана некоторая система уравнений и начальное ее состояние при , то ее эволюцию в последующие моменты времени можно объяснить посредством задания точных законов микрокинетики изучаемого явления. Однако подход, основанный на решении уравнения Лиувилля, требует глубоких знаний методов статистической физики и законов микрокинетики конденсации в дисперсных системах. В связи с этим воспользуемся более традиционным подходом [1] и уже на самом первом этапе рассуждений откажемся от точных законов микрокинетики, умышленно заменим их вероятностным законом. Такая операция обладает, по меньшей мере, одним достоинством - ясностью.

Итак, рассмотрим аэрозольную частицу (каплю) объемом V, помещенную в переохлажденный пар. Как уже отмечалось выше, состояние каждой капли в аппарате определяется координатами ее центра массы, размером и формой. Чтобы упростить анализ, примем, что основной характеристикой капли является ее объем, что дает возможность рассматривать ее как шар, объем которого равен объему капли.

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

.

Под влиянием каждой отдельной флуктуации, вызванной изменением условий взаимодействия отдельных частиц со средой, происходит малое отклонение скорости роста от ее среднего значения. Если мы не хотим входить в детали динамики взаимодействия системы многих частиц со средой (газом), то единственное утверждение, которое можно высказать относительно флуктуаций, заключается в том, что они весьма многочисленны и чрезвычайно нерегулярны по своей величине. Это утверждение дает нам необходимую основу для применения закона больших чисел и теории вероятностей в описании процесса. Мы не можем считать величину  заданной функцией времени, однако можем сделать разумные предположения о влиянии флуктуации при усреднении по большому числу одинаковых частиц (то есть по их ансамблю). Аналогично мы не можем предсказать скорость роста или объем капли в каждый момент времени τ, но можем предсказать средний результат большого числа экспериментов, выполненных в одинаковых условиях. Следовательно, подход к решению данного уравнения должен отличаться от традиционной детерминированной задачи роста частицы. Данное является типичным представителем класса так называемых стохастических (или случайных) уравнений и относится к теории стохастических процессов.

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

,

где: ,  - масса и теплоемкость капли;  - коэффициент теплоотдачи;  -разность между температурой  пара на поверхности капли и средней температурой капли . Величина , определяется по известной формуле:

,

где: второй член учитывает частичное восстановление температуры в паровом слое вблизи поверхности капли;

 - число Прандтля;

 - теплоемкость пара.

Во многих случаях поправкой, учитывающей частичное восстановление температуры торможения на поверхности капли, можно пренебречь.

Коэффициент теплоотдачи для режима континуума, когда длина свободного пробега молекул значительно меньше диаметра частицы, можно определять, используя эмпирическое соотношение:

. (1)

В частности известны формулы: Ранца - Маршала (цит. по [2]):

,

а также Дрейка:

.

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

,   (2)

где, как и ранее, ;  - исходное число Нуссельта (1). Для режима течения со скольжением, охватывающего и переход к свободномолекулярному течению, можно использовать формулу [3]:

 . (3)

Расчетами [4] установлено, что частицы размером  м  движутся с сохранением температурного равновесия. В работе [5] подробно проанализированы теоретические и экспериментальные данные по этому вопросу и сделан вывод, что согласно теории, в области Стокса при , возрастание интенсивности испарения на передней поверхности частицы компенсируется его уменьшением на задней поверхности частицы. Таким образом, полная скорость испарения не меняется. К тому же, при очень высоких числах  вылетающие из частицы молекулы уносятся от нее так быстро, что процесс напоминает испарение в вакууме. Этот процесс несущественен для частиц диаметром меньше 40 мкм, так как их движение относительно среды происходит при малых числах Re, или они быстро замедляются до малых чисел Re. То есть, при оценках скорости испарения и конденсации, движением капли можно пренебречь.

Выводы

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

            (4)

где  - вероятность нахождения частицы в -м состоянии в момент времени ;

 - условная вероятность перехода частицы из состояния  в состояние  за время ,

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

Рецензенты:

  • Веригин Александр Николаевич, доктор технических наук, профессор, заведующий кафедрой машин и аппаратов химической промышленности ФГОУ ВПО «Санкт-Петербургский государственный технологический институт (технический университет)», г. Санкт-Петербург.
  • Холоднов Владислав Алексеевич, доктор технических наук, профессор, заведующий кафедрой математического моделирования и оптимизации химико-технологических процессов ФГОУ ВПО «Санкт-Петербургский государственный технологический институт (технический университет)», г. Санкт-Петербург.