Москва, г. Москва и Московская область, Россия
Тверь, Тверская область, Россия
Тверь, Тверская область, Россия
Москва, г. Москва и Московская область, Россия
Москва, г. Москва и Московская область, Россия
Москва, г. Москва и Московская область, Россия
ВАК 01.04.2007 Физика конденсированного состояния
ВАК 01.04.2015 Физика и технология наноструктур, атомная и молекулярная физика
ВАК 05.02.2000 Машиностроение и машиноведение
УДК 62 Инженерное дело. Техника в целом. Транспорт
ГРНТИ 29.17 Физика газов и жидкостей. Термодинамика и статистическая физика
ГРНТИ 55.03 Машиноведение и детали машин
Рассмотрен перенос электромагнитного излучения в системах сферических и цилиндрических дисперсных частиц с различными оптическими характеристиками. Методом конечных элементов произведен расчет температурного поля для системы сферических частиц в двумерном приближении. Рассмотрена возможность управления процессом теплопереноса в таких системах. Предложен численный метод определения плотности тепловых источников для двух сферических частиц. Метод применим как для частиц с постоянными оптическими характеристиками, так и для частиц, диэлектрическая проницаемость которых следует закону Керра или зависит от четных степеней модуля электрического вектора.
дисперсные системы, электромагнитное излучение, теплоперенос, метод конечных элементов, управление, наночастицы
Введение
В современной действительности все большую актуальность приобретает разработка эффективных способов воздействия на дисперсные системы (в том числе и наносистемы) с различными физико-химическими свойствами [1- 4]. Это связано как с глобальной экологической проблемой загрязнения атмосферы, когда аэрозоли с их огромной суммарной поверхностью оказывают существенное влияние на все физико-химические процессы в атмосфере, так и с развитием современного производства, на различных технологических стадиях которого используются дисперсные вещества [5]. В качестве управляющего фактора рассматриваются электромагнитные поля и исследуются процессы переноса электромагнитного излучения и тепла в дисперсных системах, подвергающихся интенсивному электромагнитному воздействию. При управлении технологическим процессом необходимо учитывать нелинейный характер зависимости различных характеристик дисперсных частиц от температуры, векторов электромагнитного поля и др. Так, исследования взаимодействия лазерного излучения с дисперсными частицами необходимы для решения таких практических задач как:
- получение веществ с заданными свойствами;
- конструирование двигателей (жидкостно-реактивных и внутреннего сгорания);
- конструирование аппаратов химических производств;
- моделирование процессов в атмосфере (таких, как создание зон просветления лазерным излучением).
При рассмотрении наносистем необходимо учесть, что по форме наночастицы условно можно разделить на следующие группы:
- трехмерные (три размера в наноинтервале), например коллоидные растворы, эмульсии, кристаллы, капли, пузырьки;
- двухмерные (один размер находится в наноинтеравале, другие могут быть сколь угодно велики) тонкие плёнки;
- одномерные частицы (поперечные размеры - в наноинтервале, а длина может быть велика), например, к ним можно отнести игольчатые аэрозоли, нанотрубки и др.
На основе изучения данных о реальных дисперсных системах и наносистемах использованы такие модели дисперсных частиц как сфера и цилиндр, являющиеся предельными случаями их форм. В работах [6, 10-12] рассмотрены задачи взаимодействия электромагнитных волн с поглощающими излучение аэрозольными частицами и наночастицами со следующими характеристиками:
- две однородные взаимодействующие сферические частицы;
- уединенная сферическая частица, имеющая неоднородное включение;
- коллективы сферических однородных частиц одинаковых радиусов;
- коллективы сферических однородных частиц различных радиусов;
- коллективы сферических неоднородных частиц одинаковых радиусов;
- коллективы сферических неоднородных частиц различных радиусов;
- пара однородных по составу цилиндрических частиц;
- пара неоднородных по составу цилиндрических частиц с коаксиальными и не коаксиальными включениями.
В работе [6] рассматривались дисперсные системы «капли воды в воздухе», при этом их форма полагалась сферической (отметим, что при определенном воздействии некоторых факторов - взаимодействия капель, электромагнитного излучения - капли становятся сфероидальными). Задача взаимодействия электромагнитного излучения с двумя каплями сферической формы, даже при условии их взаимодействия, может быть решена аналитически [6]. Для дисперсных частиц цилиндрической формы при определенных условиях тоже могут быть получены аналитические решения [6,7]. Отметим, что реальные аэрозольные системы содержат полидисперсные частицы в твердой фазе, имеющие неоднородную структуру. Аналитическое решение систем уравнений переноса (Максвелла (если дисперсная система находится в электромагнитном поле), теплопереноса, Навье-Стокса) оказывается чрезвычайно громоздким, и расчеты необходимо проводить численно (задача, как правило, является нелинейной). В разделе 1 настоящей работы описывается способ расчета тепловых источников внутри поглощающих частиц с учетом влияния соседней частицы, полученный на основе теории Ми. В общем случае частицы полагаются оптически нелинейными, а именно, подчиняющимися закону Керра. Раздел 2 посвящен решению уравнения Лапласа в длинноволновом приближении, нахождению напряженности электромагнитного поля внутри вытянутых цилиндрических частиц и определению инициированных электромагнитным полем источников тепла. В разделе 3 проведен расчет температуры в неоднородных дисперсных частицах в полидисперсных системах методом конечных элементов. Тепловые источники определяются на основе решения задачи о взаимодействии электромагнитного поля с частицами.
Перенос электромагнитного излучения в системе сферических частиц
Рассматриваются парные взаимодействия сферических дисперсных частиц, находящихся под воздействием монохроматического электромагнитного излучения. Поглощенная частицами энергия электромагнитного поля, выделяющаяся в единице объема j-ой частицы в единицу времени, представляет собой тепловой источник , инициированный излучением и определяется значением квадрата модуля амплитуды вектора напряженности электрического вектора в данной точке [1]. Ранее разработан метод расчета компонент векторов напряженности электромагнитного поля на основе теории Ми [8] и теории представления групп [9]. Для решения задачи связываем с центрами частиц O1 и O2 радиусами R1
и R2 и направлениями векторов (напряженности внешнего электромагнитного поля) и
(волнового вектора) декартовой системой координат (рис.1) [9].
|
Рис. 1. Система двух частиц в поле электромагнитного излучения |
Электрический вектор падающей волны зависит от времени по гармоническому закону и для его амплитуды можно записать [10]:
(1)
где k(j) - волновое число в веществе j,
- циклическая частота,
- диэлектрическая проницаемость среды,
- удельная проводимость, с - скорость света. Число j принимает значения: j=0 (внешняя среда), j=1 (частица № 1), j=2 (частица №2). Для магнитного вектора можно записать аналогичное уравнение. В общем случае диэлектрическая проницаемость зависит от электрического поля, то есть
,
,
где - параметр нелинейности.
Учесть влияние k-й частицы можно, если представить внешнее по отношению к j-й частице поле в виде [8]
, (2)
где представляет собой вектор падающей волны,
- волны, рассеянной на соседней частице. Векторы
и
можно представить через электрический U и магнитный V потенциалы Дебая [6]. Потенциалы Дебая поля падающего излучения выражены в системе координат х’у’z’, а потенциалы Дебая поля, рассеянного на k-й частице в системе координат xkykz (k=1,2). Далее полученные выражения представляются в виде разложений по собственным сферическим функциям j-й сферы
. Углы поворота Эйлера p/2 - b, θ, 3p/2 полностью определяют вращение g относительно точки O1, переводящее систему координат x’y’z’ в систему x1y1z. Данный подход может быть применен и для частиц с нелинейными свойствами, однако уравнения для определения полей имеют другой вид. Полагая, что параметр линеаризации можно рассматривать как малый параметр, получим в нелинейном случае неоднородные уравнения Гельмгольца для определения электрического вектора. Действительно, рассматривая нелинейную зависимость Керра для амплитуды вектора электрической напряженности в каждой частице, уравнения для
в каждой из частиц будут иметь вид:
, (3)
где - комплексная диэлектрическая проницаемость,
. Система (3) линеаризуется, если величина
является малым параметром и соответственно,
. При этом основная часть электрического уравнения определяется на основе уравнения (1), а для
получается неоднородное уравнение Гельмгольца с правой частью:
.
Частное решение такого уравнения может быть найдено с использованием разложения правой части по собственным функциям решения уравнения (1). Вклад полученного частного решения отражается на величине коэффициентов рассеяния и поглощения.
Отметим, что полученные решения могут быть использованы и при рассмотрении наночастиц относительно больших размеров. Однако при рассмотрении наночастиц с диаметром порядка нескольких наноразмеров необходимо учитывать вклад продольных волн в общую величину электрического поля.
Для определения результатов воздействия электромагнитного поля на частицу или каплю, прежде всего её нагрев в результате поглощения воздействующей электромагнитной волны, необходимо решать уравнения теплопроводности с соответствующим тепловым источником электромагнитной природы.
Для расчета источника тепла частица разбивается на малые объемы, внутри которого практически постоянно. Используя такую схему, мы рассчитываем средние по объему частицы значения плотности источников тепла, отнесенные к интенсивности падающей волны q* для различных ориентаций системы частиц относительно векторов поля и различных параметров нелинейности при выбранном законе нелинейности, например закон Керра. Такой метод ввиду выше приведенной правой части неоднородного уравнения Гельмгольца оказывается весьма эффективным и для численного определения теплового источника в каждой из частиц. В таблице в качестве примера приведены расчеты для частиц водного аэрозоля радиусами R1 = R2 =1 мкм. Расстояние между центрами частиц R = 20 мкм. Длина волны инициирующего теплоперенос лазерного излучения l = 10.6 мкм. Температура окружающей среды T = 273 К. Коэффициент теплопроводности среды cе = 24,1 мВт/(мК), коэффициенты теплопроводности вещества частиц равны c1 = c2 = 566 мВт/мК×. Значения комплексных показателей преломления вещества частиц (воды) 1.173 + i ×0.0823. Показатель преломления среды ne = 1. В данном случае коэффициент нелинейности очень мал (10-5) и почти не влияет на величину тепловых источников.
Таблица
Зависимость q* от углов Эйлера
|
θ, град |
b, град |
qj*× 10-4, м-1 |
|
90 |
0 |
8,7 |
|
90 |
90 |
5,7 |
|
45 |
45 |
6,3 |
Из расчетов следует [11], что средняя плотность источников тепла зависит от расположения линии, соединяющей центры двух частиц относительно вектора напряженности электрического поля и волнового вектора. Даже на больших расстояниях между частицами это расположение является значимым.
Ранее в работах [10-12] было показано, что не только электромагнитное, но и тепловое взаимодействие дисперсных частиц может быть заметным и при больших расстояниях между их центрами. Проведенные указанным методом расчеты показывают, что вблизи границы дисперсной частицы распределение носит более неоднородный характер, чем в центральных областях, что позволяет говорить о существовании теплового поверхностного слоя частицы, который характеризуется сильно неоднородным температурным полем в его объеме.
Распространение электромагнитного излучения в системе двух вытянутых одномерных частиц (цилиндров)
Другой распространенной геометрической моделью дисперсных частиц является цилиндр или вытянутый сфероид. Разработанный нами метод расчета квадрата амплитуды для данной геометрии [11] применим при условии, когда радиус цилиндров много меньше длины волны падающего излучения. Такое ограничение позволяет использовать для решения бицилиндрическую систему координат, что значительно упрощает решение и расчеты [13, 14].
Для указанных систем частиц проведены расчеты отношения величин 
(то есть отношения квадрата модуля амплитуды светового вектора в выбранной точке внутри частицы к квадрату модуля амплитуды вектора напряженности электрического поля падающей волны), характеризующего плотность источников тепла в зависимости от радиуса и угловой координаты в сечениях, перпендикулярных оси цилиндра.
Рассмотрены системы неоднородных по составу цилиндров с различными значениями комплексной диэлектрической проницаемости вещества частицы. При расчетах использовались следующие значения: l =2 мм; e1 =1 (для воздуха); (для льда).
Проведены вычислительные эксперименты для двухслойных цилиндров с концентрическими (ранее в работе [11]) и неконцентрическими слоями (внутренняя область частицы представляет собой структуру со случайно распределенным радиусом в зависимости от угла по сечению цилиндра) различного химического состава.
|
Рис. 2. Система цилиндров с неконцентрическими слоями |
, координатная поверхность цилиндра в бицилиндрической системе координат – τ1. Данное решение формально сводится к приведенному в работе [11], но здесь радиус внутренней границы задается заново для каждой рассматриваемой точки.
Найдены распределения величины 
внутри цилиндрической частицы радиуса R0=1 мкм при радиусе окружности выбранного сечения 0.95 мкм со случайным распределением внутренней границы внешнего слоя. В этом случае точка окружности может принадлежать как внешнему, так и внутреннему слою. Были выбраны следующие значения комплексной диэлектрической проницаемости: для внутреннего слоя – диэлектрическая проницаемость для льда; для внешнего – модельные значения, соответствующие сильному поглощению.
|
а) б)
Рис. 3. Распределение величины b по сечению цилиндра
|
Рис. 3а,б соответствуют различным конфигурациям внутренней границы внешнего слоя, выбранной случайным образом. Полученные распределения имеют форму, зависящую от толщины внешнего слоя в выбранной точке, и вписываются в диаграмму, рассчитанную для радиуса внутренней границы внешнего слоя, соответствующей минимальной толщине внешнего слоя.
Теплоперенос в дисперсных системах с различной структурой и конфигурацией
В разделе 1 был предложен метод расчета плотности теплового источника, обусловленного поглощением при воздействии на частицу электромагнитной волны. Метод применим как для частиц с постоянными оптическими характеристиками, так и для частиц, диэлектрическая проницаемость которых следует закону Керра или зависит от четных степеней модуля электрического вектора. Для определения температуры нагреваемых частиц под действием тепловых источников электромагнитной природы использовался метод конечных элементов. Тепловые источники определялись с использованием результатов задачи взаимодействия электромагнитных волн с частицами. В модельных схемах рассмотрены сферические частицы реальных для природных аэрозолей размеров, однако разработанная программа позволяет проводить расчеты и для частиц произвольной формы.
Использовалась схема, представляющая собой прямоугольную площадку (90 на 70 мкм) с размещенными на ней сферическими частицами радиуса 5 и 10 мкм. Количество частиц варьировалось от одной до девяти, расположение частиц на площадке также варьировалось.
Основные результаты исследования процесса теплопереноса в дисперсных системах с однородными по составу частицами под действием произвольно заданных тепловых источников опубликованы в работе [12]. Проводились расчеты температурного поля и для систем, содержащих частицы различного химического состава.
Из анализа результатов вычислительных экспериментов следует, что присутствие частиц различных по размерам и теплофизическим свойствам приводит к увеличению неоднородности температурного поля в системе, причем при линейной зависимости источника от температуры проявляются нелинейные эффекты, а при нелинейной зависимости эти эффекты усиливаются. Если система полидисперсная, то существенное влияние на такое распределение оказывает присутствие более крупных частиц, а при граничных условиях третьего рода на внешней границе системы влияние коллективных эффектов может изменять температуру в 1,5 и более раз. Это влияние тем больше, чем система более неоднородна (по размеру частиц и составу).
На основе общего подхода к решению задачи оптимизации процесса теплопереноса в системе из N частиц, обусловленного тепловыми источниками в каждой частице, получены расчетные формулы для определения оптимальных параметров конфигурации системы из двух сферических частиц и оптимального внешнего воздействия. Характерные результаты численных модельных расчетов по полученным формулам приведены в [11,12]. При решении задачи оптимизации использовались решения задачи теплопереноса под действием электромагнитного излучения для двух сферических частиц [10,11]. В расчетах использовались такие управляющие параметры, как мощность внешнего электромагнитного источника, определяемого из задачи электродинамики, и конфигурация расположения частиц.
Найден минимум функционала для соответствующей системы с распределенными параметрами. Так как конфигурации расположения частиц выбирались произвольно, что существенно усложняло расчеты, то в качестве характерных размеров системы рассматривались продольный и поперечный диаметры скопления частиц, расстояние между центрами двух наиболее крупных частиц и подобные естественные геометрические характеристики. Варьировались размеры частиц, вещество частиц и некоторые другие величины.
Мы полагали, что процесс теплопереноса в такой системе является квазистационарным, и рассматривали последовательную оптимизацию, считая, что в этом случае необходимо минимизировать следующий функционал [6]
(3)
|
Рис. 4. Оптимальная зависимость радиуса частицы от отношения тепловых источников при η = π/3, m = 10, ν = 5, k = 22; 44 |
где - заданная температура. Другие варианты минимизации функционала представлены в работах [11,12]. Полагая, что поверхности всех частиц имеют одинаковую температуру (на основе второго закона термодинамики), последнее выражение преобразуем к виду:
(4)
Для системы «две сферические частицы - континуальная среда» в [11,12] было получено аналитическое решение для распределения температуры в системе из двух частиц в бисферической системе координат. Дифференцируя это решение, мы получили уравнения для экстремальных параметров задачи: (где
средние по объему частицы плотности тепловых источников, 
- радиусы частиц), что соответствует решению задачи оптимизации (4). В качестве управляющего параметра можно рассматривать среднее по объёму значение теплового источника (в частном случае, мощность внешнего воздействия, например, лазерного воздействия) и расстояние между центрами сфер или осями цилиндров (в ортогональной плоскости к осям). Решение для случая двух сферических частиц, где в качестве управляющего параметра выбрана мощность источника и расстояние между двумя центрами сфер
, было получено в работах [11,12].
Решая системы уравнений
, (5)
можно определить оптимальные параметры. При численных модельных расчетах полагалось, что (на поверхности частицы). Для уменьшения числа параметров вводились величины [6]
. (6)
Параметры варьировались в пределах: .
В процессе численного моделирования было установлено, что оптимальные значения отношений характерных размеров частиц зависят от вещества частиц, их оптических свойств и геометрии системы. Характерные диаграммы таких зависимостей приведены на рис. 4.
Из полученных результатов расчета отношения источников для системы «вода-воздух» следует, что зависимость источника от координаты является немонотонной.
В последние годы проводилось интенсивное исследование кластеров, имеющих прямое отношение к процессам ассоциации твердых частиц (образованию кластеров в облаках, коагуляции частиц в дымах, образованию структур при релаксации металлического пара и т.д.) численными методами. Так, было получено, что при воздействии на крупную аэрозольную частицу размером порядка 10 - 100 мкм лазерным излучением большой мощности, частица распадается на несколько осколков (образуется кластер). Если внутренняя энергия таких осколков близка к внутренней энергии систем, описываемых уравнениями сплошной среды, то предложенная теория может быть использована для описания теплопереноса в системах, включающих кластеры.
Заключение
1. Из проведенных расчетов для системы двух сферических частиц следует, что среднее значение плотности источников тепла существенно зависит от расположения системы двух частиц относительно волнового вектора и вектора напряженности электрического поля. Показано, что даже на больших расстояниях между частицами влияние второй частицы на распределение плотности источников тепла по объему первой значительно. Не только электромагнитное, но и тепловое взаимодействие дисперсных частиц может быть заметным при больших расстояниях между их центрами. Расчеты показывают, что вблизи границы дисперсной частицы распределение температуры носит более неоднородный характер, чем в центральных областях, что позволяет нам говорить о существовании теплового поверхностного слоя частицы, который характеризуется сильно неоднородным температурным полем в его объеме.
2. Предложен численный метод определения плотности тепловых источников для двух сферических частиц. Метод применим как для частиц с постоянными оптическими характеристиками, так и для частиц, диэлектрическая проницаемость которых следует закону Керра или зависит от четных степеней модуля электрического вектора.
3. Проведенные вычислительные эксперименты для двухслойных цилиндров с концентрическими и неконцентрическими слоями различного химического состава со случайным образом выбранной внутренней границей между слоями показывают, что распределение плотности внутренних источников тепла (пропорциональных квадрату амплитуды электрического вектора поглощенного поля) в зависимости от угла в полярной системе координат, имеют на выбранном сечении форму, определяемую размером внешнего слоя и вписывается в диаграмму, рассчитанную для радиуса внутренней границы внешнего слоя, соответствующей минимальной толщине внешнего слоя.
4. Из анализа результатов численных экспериментов, выполненных на основе метода конечных элементов и проведенных для изучения переноса тепла в дисперсных системах, содержащих частицы, отличающиеся по составу, конфигурации и размерам следует, что в таких неоднородных системах возрастает неоднородность температурного поля, причем даже при линейной зависимости источника от температуры проявляются нелинейные эффекты, а при нелинейной зависимости эти эффекты усиливаются. Если система полидисперсная, то существенное влияние на такое распределение оказывает присутствие более крупных частиц, а если на поверхности частиц реализуются граничные условия третьего рода, влияние коллективных эффектов может изменять температуру в 1,5 и более раз. Это влияние тем больше, чем система более неоднородна (по размеру частиц и составу).
1. Пришивалко, А.П. Оптические и тепловые поля внутри светорассеивающих частиц / А.П.Пришивалко. - Мн.: Наука и техника, 1983. - 190 с.
2. Лушников, А.А. Физико-химические процессы образования атмосферных аэрозолей / А.А. Лушников, А.А. Загайнов, И.Е. Аграновский, Ю.С. Любовцева // Журнал физической химии. - 2008. - Т.82. - № 10. - С. 1950-1958.
3. Лушников, А.А. Начальная стадия образования аэрозоля из пересыщенных паров / А.А. Лушников, В.А. Загайнов, Ю.С. Любовцева // Журнал физической химии. - 2018. - Т.92. - № 3. - С. 501- 507.
4. Лушников, А.А. Образование наноаэрозолей в тропосфере под действием космического излучения / А.А. Лушников, В.А. Загайнов, Ю.С. Любовцева, А.Д. Гвишиани // Известия РАН. Физика атмосферы и океана. - 2014. - Т.50. - № 2. - С.175-184.
5. Barakov, A. V. Study of heat and mass transfer in the dryer with a pseudo liquefied layer of a disperse material / A. V. Barakov, A. P. Byrdin , A. A. Nadeev // Scientific Herald of the Voronezh State University of Architecture And Civil Engineering. Construction And Architecture. - 2016. - Vol.3. - № 31. - P.18-26.
6. Кривенко, И.В. Моделирование воздействия электромагнитного излучения на природные аэрозоли/ И.В. Кривенко, М.А. Смирнова // Математическая физика и компьютерное моделирование. - 2019. - Т. 22. - № 4. - С. 64-79.
7. Смирнова, М.А. Теплоперенос в дисперсных системах различного состава и конфигураций: автореферат диссертации на соискание ученой степени канд.физ.-мат.наук: 02.00.04 / М. А. Смирнова. - Тверь, 2000. - 24 с.: ил. библиогр.: С. 22-24.
8. Борн, М. Основы оптики / М. Борн, Э. Вольф. - М.: Наука, 1970. - 850 с.
9. Виленкин, Н.Я. Специальные функции и теория представления групп / Н.Я. Виленкин - М.: Наука, 1956. - 588 с.
10. Гамаюнов, Н.И. Особенности распространения электромагнитного излучения и инициированного им теплопереноса в системе «аэрозольные частицы - окружающая среда» / Н.И. Гамаюнов, И.В. Кривенко, Л.А. Уварова, Ю.З. Бондарев // Журнал физической химии. - 1997. - Т.71. - № 2. - С.2270-2274.
11. Uvarova, L.A. Modeling of Heat Transfer in the System of Small Spherical and Cylindrical Particles under the Action Electromagnetic Radiation / L.A. Uvarova, I.V. Krivenko, M.A. Smirnova, A.F. Ivannikov //Nonlinearity: Problems, Solutions and Applications / ed. L.A. Uvarova. - NY: Nova Science Publishers, Inc, 2017. - С. 519-541.
12. Уварова, Л.А. Малые поглощающие частицы в электромагнитном поле и теплоперенос в дисперсных системах / Л.А. Уварова [и др]. - М.: Янус-К, 2014. - 192 с.
13. Корн, Г. Справочник по математике для научных работников и инженеров / Г. Корн. - М.: Наука, 1984. - С. 831.
14. Морс, Ф.М. Методы теоретической физики / Ф.М. Морс, Г. Фешбах. - М.: Изд-во иностр. литературы, 1958. - С. 415.
