В материаловедении приближения эффективной среды ( EMA ) или теория эффективной среды ( EMT ) относятся к аналитическому или теоретическому моделированию, которое описывает макроскопические свойства композитных материалов . EMA или EMT разрабатываются путем усреднения множественных значений компонентов, которые непосредственно составляют композитный материал. На уровне компонентов значения материалов различаются и неоднородны . Точный расчет многих значений компонентов практически невозможен. Однако были разработаны теории, которые могут производить приемлемые приближения, которые, в свою очередь, описывают полезные параметры, включая эффективную диэлектрическую проницаемость и проницаемость материалов в целом. В этом смысле приближения эффективной среды являются описаниями среды (композитного материала) на основе свойств и относительных долей его компонентов и выводятся из расчетов, [1] [2] и теории эффективной среды . [3] Существует две широко используемые формулы. [4]
Эффективная диэлектрическая и магнитная проницаемость являются усредненными диэлектрическими и магнитными характеристиками микронеоднородной среды. Они обе были выведены в квазистатическом приближении, когда электрическое поле внутри частицы смеси можно считать однородным. Поэтому эти формулы не могут описать эффект размера частиц. Было предпринято много попыток улучшить эти формулы.
Существует много различных приближений эффективной среды, [5] каждое из которых более или менее точно в определенных условиях. Тем не менее, все они предполагают, что макроскопическая система однородна, и, что типично для всех теорий среднего поля, они не могут предсказать свойства многофазной среды вблизи порога перколяции из -за отсутствия в теории дальнодействующих корреляций или критических флуктуаций.
Рассматриваемыми свойствами обычно являются проводимость или диэлектрическая проницаемость [6] среды. Эти параметры взаимозаменяемы в формулах в целом ряде моделей из-за широкой применимости уравнения Лапласа. Задачи, выпадающие из этого класса, в основном относятся к области упругости и гидродинамики, что обусловлено тензорным характером эффективных констант среды более высокого порядка.
EMA могут быть дискретными моделями, например, применяемыми к сетям резисторов, или теориями континуума, применяемыми к упругости или вязкости. Однако большинство современных теорий испытывают трудности в описании просачивающихся систем. Действительно, среди многочисленных приближений эффективной среды только симметричная теория Бруггемана способна предсказать порог. Эта характерная черта последней теории ставит ее в ту же категорию, что и другие теории среднего поля критических явлений . [ необходима цитата ]
Для смеси двух материалов с диэлектрическими проницаемостями и с соответствующими объемными долями и Д.А.Г. Брюггеман предложил формулу следующего вида: [7]
Здесь положительный знак перед квадратным корнем должен быть изменен на отрицательный знак в некоторых случаях, чтобы получить правильную мнимую часть эффективной комплексной диэлектрической проницаемости, которая связана с затуханием электромагнитных волн. Формула симметрична относительно перестановки ролей 'd' и 'm'. Эта формула основана на равенстве
где — скачок потока электрического смещения по всей поверхности интегрирования, — компонента микроскопического электрического поля, нормальная к поверхности интегрирования, — локальная относительная комплексная диэлектрическая проницаемость, которая принимает значение внутри выбранной металлической частицы, значение внутри выбранной диэлектрической частицы и значение снаружи выбранной частицы, — нормальная компонента макроскопического электрического поля. Формула (4) вытекает из равенства Максвелла . Таким образом, в подходе Бруггемана рассматривается только одна выбранная частица. Взаимодействие со всеми остальными частицами учитывается только в приближении среднего поля, описываемом выражением . Формула (3) дает разумную резонансную кривую для плазмонных возбуждений в металлических наночастицах , если их размер составляет 10 нм или меньше. Однако она не способна описать размерную зависимость для резонансной частоты плазмонных возбуждений, которые наблюдаются в экспериментах [8]
Не теряя общности, рассмотрим изучение эффективной проводимости (которая может быть как постоянного, так и переменного тока) для системы, состоящей из сферических многокомпонентных включений с различной произвольной проводимостью. Тогда формула Бруггемана примет вид:
В системе евклидовой пространственной размерности, имеющей произвольное число компонентов, [9] сумма производится по всем компонентам. и являются соответственно долей и проводимостью каждого компонента, а является эффективной проводимостью среды. (Сумма по 's равна единице.)
Это обобщение уравнения (1) на двухфазную систему с эллипсоидальными включениями проводимости в матрицу проводимости . [10] Доля включений равна и система является размерной. Для случайно ориентированных включений,
где 's обозначают соответствующий дублет/триплет факторов деполяризации, который регулируется соотношениями между осями эллипса/эллипсоида. Например: в случае круга ( , ) и в случае сферы ( , , ). (Сумма по 's равна единице.)
Наиболее общий случай, к которому был применен подход Бруггемана, включает бианизотропные эллипсоидальные включения. [11]
На рисунке показана двухкомпонентная среда. [9] Рассмотрим заштрихованный объем проводимости , примем его за сферу объемом и предположим, что он вложен в однородную среду с эффективной проводимостью . Если электрическое поле вдали от включения равно , то элементарные соображения приводят к дипольному моменту , связанному с объемом
Эта поляризация производит отклонение от . Если среднее отклонение должно исчезнуть, общая поляризация, суммированная по двум типам включений, должна исчезнуть. Таким образом
где и — соответственно объемная доля материала 1 и 2. Это можно легко распространить на систему размерности , которая имеет произвольное число компонентов. Все случаи можно объединить, чтобы получить уравнение (1).
Уравнение (1) также можно получить, потребовав, чтобы отклонение тока исчезло. [12] [13] Оно было выведено здесь из предположения, что включения являются сферическими, и его можно модифицировать для форм с другими факторами деполяризации, что приводит к уравнению (2).
Также доступен более общий вывод, применимый к бианизотропным материалам. [11]
Основное приближение заключается в том, что все домены находятся в эквивалентном среднем поле. К сожалению, это не так вблизи порога перколяции, где система управляется крупнейшим кластером проводников, который является фракталом, и дальними корреляциями, которые полностью отсутствуют в простой формуле Бруггемана. Пороговые значения в целом не предсказаны правильно. Это 33% в EMA, в трех измерениях, далеко от 16%, ожидаемых из теории перколяции и наблюдаемых в экспериментах. Однако в двух измерениях EMA дает порог 50% и, как было доказано, моделирует перколяцию относительно хорошо. [14] [15] [16]
В приближении Максвелла Гарнетта [17] эффективная среда состоит из матричной среды с и включений с . Максвелл Гарнетт был сыном физика Уильяма Гарнетта и был назван в честь друга Гарнетта, Джеймса Клерка Максвелла . Он предложил свою формулу для объяснения цветных картин, которые наблюдаются в стеклах, легированных металлическими наночастицами. Его формула имеет вид
где - эффективная относительная комплексная диэлектрическая проницаемость смеси, - относительная комплексная диэлектрическая проницаемость фоновой среды, содержащей малые сферические включения относительной диэлектрической проницаемости с объемной долей . Эта формула основана на равенстве
где — абсолютная диэлектрическая проницаемость свободного пространства , а — электрический дипольный момент одиночного включения, индуцированный внешним электрическим полем E. Однако это равенство справедливо только для однородной среды и . Кроме того, формула (1) игнорирует взаимодействие между одиночными включениями. В силу этих обстоятельств формула (1) дает слишком узкую и слишком высокую резонансную кривую для плазмонных возбуждений в металлических наночастицах смеси. [18]
Уравнение Максвелла Гарнетта гласит: [19]
где — эффективная диэлектрическая проницаемость среды, включений и матрицы; — объемная доля включений.
Уравнение Максвелла Гарнетта решается следующим образом: [20] [21]
пока знаменатель не исчезнет. Простой калькулятор MATLAB, использующий эту формулу, выглядит следующим образом.
% Этот простой калькулятор MATLAB вычисляет эффективную диэлектрическую проницаемость % смеси материала включений в базовой среде % в соответствии с теорией Максвелла Гарнетта % ВХОДЫ: % eps_base: диэлектрическая проницаемость основного материала; % eps_incl: диэлектрическая проницаемость материала включений; % vol_incl: объемная доля материала включений; % ВЫХОД: % eps_mean: эффективная диэлектрическая проницаемость смеси.функция eps_mean = МаксвеллГарнеттФормула ( eps_base, eps_incl, vol_incl ) small_number_cutoff = 1e-6 ; if vol_incl < 0 || vol_incl > 1 disp ( 'ВНИМАНИЕ: объемная доля материала включения выходит за пределы диапазона!' ); end factor_up = 2 * ( 1 - vol_incl ) * eps_base + ( 1 + 2 * vol_incl ) * eps_incl ; factor_down = ( 2 + vol_incl ) * eps_base + ( 1 - vol_incl ) * eps_incl ; if abs ( factor_down ) < small_number_cutoff disp ( 'ВНИМАНИЕ: эффективная среда является единственной!' ); eps_mean = 0 ; else eps_mean = eps_base * factor_up / factor_down ; end end
Для вывода уравнения Максвелла Гарнетта мы начнем с массива поляризуемых частиц. Используя концепцию локального поля Лоренца, мы получаем соотношение Клаузиуса-Моссотти :
Где - число частиц в единице объема. Используя элементарную электростатику, мы получаем для сферического включения с диэлектрической проницаемостью и радиусом поляризуемость :
Если мы объединим с уравнением Клаузиуса-Моссотти , мы получим:
Где - эффективная диэлектрическая проницаемость среды включений; - объемная доля включений.
Поскольку модель Максвелла Гарнетта представляет собой композицию матричной среды с включениями, мы усиливаем уравнение:
В общих чертах, ожидается, что ЭМА Максвелла Гарнетта будет действителен при низких объемных долях , поскольку предполагается, что домены пространственно разделены и электростатическое взаимодействие между выбранными включениями и всеми другими соседними включениями пренебрегается. [22] Формула Максвелла Гарнетта, в отличие от формулы Бруггемана, перестает быть верной, когда включения становятся резонансными. В случае плазмонного резонанса формула Максвелла Гарнетта верна только при объемной доле включений . [23] Была изучена применимость приближения эффективной среды для диэлектрических многослойных структур [24] и металл-диэлектрических многослойных структур [25] , что показало, что существуют определенные случаи, когда приближение эффективной среды не выполняется, и необходимо проявлять осторожность при применении теории.
Уравнение Максвелла Гарнетта описывает оптические свойства нанокомпозитов, которые состоят из набора идеально сферических наночастиц. Все эти наночастицы должны иметь одинаковый размер. Однако из-за эффекта ограничения оптические свойства могут зависеть от распределения размеров наночастиц. Как показали Батти и др. [26] , уравнение Максвелла Гарнетта можно обобщить, чтобы учесть это распределение.
и — радиус наночастиц и распределение по размерам, соответственно. и — средний радиус и объемная доля наночастиц, соответственно. — первый электрический коэффициент Ми. Это уравнение показывает, что классическое уравнение Максвелла Гарнетта дает ложную оценку объемной доли наночастиц, когда распределением по размерам нельзя пренебречь.
Уравнение Максвелла Гарнетта описывает только оптические свойства набора идеально сферических наночастиц. Однако оптические свойства нанокомпозитов чувствительны к распределению формы наночастиц. Чтобы преодолеть это ограничение, Y. Battie и др. [27] разработали теорию эффективной среды с распределенной формой (SDEMT). Эта теория эффективной среды позволяет рассчитать эффективную диэлектрическую функцию нанокомпозита, который состоит из набора эллипсоидальных наночастиц, распределенных по форме.
с
Факторы деполяризации ( ) зависят только от формы наночастиц. — распределение факторов деполяризации. f — объемная доля наночастиц.
Теория SDEMT использовалась для извлечения распределения формы наночастиц из спектров поглощения [28] или эллипсометрических спектров. [29] [30]
Была предложена новая формула, описывающая эффект размера. [18] Эта формула имеет вид
где a — радиус наночастицы, а — волновое число. Здесь предполагается, что временная зависимость электромагнитного поля задается множителем В данной работе использовался подход Бруггемана, но электромагнитное поле для электрической дипольной моды колебаний внутри выбранной частицы вычислялось без применения квазистатического приближения . Таким образом, функция обусловлена неоднородностью поля внутри выбранной частицы. В квазистатической области ( , т.е. для Ag эта функция становится постоянной и формула (5) становится идентичной формуле Бруггемана.
Формула эффективной проницаемости смесей имеет вид [18]
Здесь — эффективная относительная комплексная проницаемость смеси, — относительная комплексная проницаемость фоновой среды, содержащей мелкие сферические включения с относительной проницаемостью с объемной долей . Эта формула выведена в дипольном приближении. Магнитная октупольная мода и все другие магнитные моды колебаний нечетных порядков здесь не учитывались. При и эта формула имеет простой вид [18]
Для сети, состоящей из высокой плотности случайных резисторов, точное решение для каждого отдельного элемента может быть непрактичным или невозможным. В таком случае случайную сеть резисторов можно рассматривать как двумерный граф , а эффективное сопротивление можно моделировать в терминах мер графа и геометрических свойств сетей. [31] Предполагая, что длина ребра намного меньше расстояния между электродами, а ребра должны быть равномерно распределены, можно считать, что потенциал падает равномерно от одного электрода к другому. Поверхностное сопротивление такой случайной сети ( ) можно записать в терминах плотности ребер (проводов) ( ), удельного сопротивления ( ), ширины ( ) и толщины ( ) ребер (проводов) как:
{{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link){{cite journal}}
: CS1 maint: multiple names: authors list (link)