УДК 539.21 Вестник СПбГУ. Математика. Механика. Астрономия. 2023. Т. 10 (68). Вып. 4 МБС 74Е05
Моделирование неидеальных контактов
при определении эффективных коэффициентов
диффузии
К. П. Фролова, Е. Н. Вильчевская, В. А. Полянский
Институт проблем машиноведения Российской академии наук, Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61
Для цитирования: Фролова К. П., Вильчевская Е. Н., Полянский В. А. Моделирование неидеальных контактов при определении эффективных коэффициентов диффузии // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2023. Т. 10(68). Вып. 4. С. 650-664. https://doi.org/10.21638/spbu01.2023.405
В работе развивается комплексный подход к учету неидеальных контактов при определении эффективных свойств различной природы: диффузионных, тепло- и электропроводящих. Под неидеальными контактами понимаются такие, при которых поля на микроуровне не являются непрерывными. Возможность создания единого подхода к учету таких контактов обусловливается подобием уравнений диффузии, тепло- и электропроводности. При этом с физической точки зрения появление неидеальных контактов может иметь разные причины, связанные как с особенностями микроструктуры материала, так и со спецификой самого процесса. Для конкретности в работе определяются эффективные коэффициенты диффузии, поскольку данный случай позволяет рассматривать в качестве причин появления неидеальных контактов как образование дефектов структуры, так и наличие характерного для массопереноса эффекта сегрегации, под которым понимается скопление примеси на границе раздела фаз. В работе обобщается и сравнивается два подхода к учету неидеальных контактов при построении моделей гетерогенного материала. В рамках первого подхода предполагается, что такие контакты могут быть учтены посредством задания скачка поля концентрации примеси или нормальной компоненты потока. При втором подходе в рассмотрение вводится неоднородность с тонким покрытием, обладающим экстремальными свойствами. Всесторонний анализ проводится на примере материала со сферическими неоднородностями. Приводятся аналитические выражения для определения тензора вклада эквивалентной неоднородности в искомое свойство, что позволяет существенно упростить процедуру обобщения разных методов гомогенизации на случай микрогетерогенного материала с неидеальными контактами.
Ключевые слова: эффективные свойства, неоднородность, неидеальные контакты, сегрегация, тензор диффузии.
1. Введение. Микроструктура материала оказывает непосредственное влияние на его макроскопические характеристики. При определении эффективных свойств в случае, когда в материале имеются ограниченные области, свойства которых отличны от свойств окружающего материала, считается, что материал состоит из матрицы и локализованных неоднородностей [1, 2]. При этом, как правило, делается предположение о непрерывности полей на границе раздела фаз. С физической точки зрения это означает наличие «идеальных контактов». В то же время при опи-
© Санкт-Петербургский государственный университет, 2023
сании ряда явлений необходимо учитывать наличие неидеальных контактов. Такие контакты могут возникать как минимум по двум разным причинам: в связи с особенностями микроструктуры и в связи со спецификой самого процесса. К первому случаю можно отнести растрескивание материала или образование отслоений, наличие проскальзывающего слоя, либо суперпроводящего/изолирующего покрытия у частиц. Ко второму — эффект сегрегации при диффузии (оседание диффундирующего вещества в дефектах). Вопросы учета неидеальных контактов при построении моделей рассматривались в литературе по отдельности в контексте описания различных процессов.
На необходимость учета неидеальных контактов при определении эффективных диффузионных свойств указали Ю. Занг и Л. Лиу (Y. Zhang and L. Liu) [3], обратив внимание, что концентрация не является непрерывной и может испытывать скачок на границе раздела фаз (данное явление известно как эффект сегрегации [4, 5]). Выразив поток массы через непрерывный градиент химического потенциала, они получили оценки Фойгта — Рейсса и Хашина — Штрикмана для эффективной подвижности диффузанта, после чего выразили через нее непосредственно эффективные коэффициенты диффузии. Иной способ учета сегрегации был предложен в работах [6, 7], где авторы ввели постоянный параметр сегрегации, задающий скачок концентрации, в модифицированный метод эффективной среды и определили макроскопический коэффициент диффузии примеси в материале со сферическими неоднородностями. Учет постоянного параметра сегрегации в методе без учета взаимодействия неоднородностей, методах Мори — Танаки и Максвелла был реализован в работе [8]. При этом наличие скачка концентрации не было учтено на этапе осреднения полей, что привело к необходимости дальнейшего уточнения использованных методов. Такой поэтапный учет был реализован авторами настоящей статьи в работе [9], дополнительно был рассмотрен метод гомогенизации Канауна — Левина [1].
Таким образом, в контексте задачи диффузии под неидеальным контактом понимался исключительно эффект сегрегации. Иная интерпретация появилась при рассмотрении задач теплопроводности. Предполагалось, что появление неидеальных контактов вызвано реальным нарушением идеальных (шероховатостью поверхности, отслоениями и т.д.). В работах [10-12] неидеальные контакты вводились в модель посредством задания постоянной величины скачка поля. Иной подход был использован в работах [13, 14], где решалась задача об определении эффективных свойств материала с неоднородностями, состоящими из конфокальных эллипсоидов с различными свойствами, после чего делалось предположение об устремлении к нулю толщины слоя и к бесконечности — его свойств (проводимости, либо сопротивляемости).
Подобие уравнений диффузии, теплопроводности и электропроводности позволяет сделать предположение о принципиальной возможности разработки единого комплексного подхода к моделированию неидеальных контактов. В связи с этим представляет интерес сравнение подходов, разработанных в литературе в контексте описания разных физических процессов. Целью настоящей работы является обобщение и сравнение двух подходов (посредством задания скачка поля и рассмотрения неоднородности с покрытием) к учету неидеальных контактов при определении эффективных свойств аналитическим способом. Для конкретности рассматривается задача об определении эффективной диффузионной проницаемости, под которой понимается эффективный тензор диффузии примеси в материале, поскольку данный случай представляется наиболее общим в связи с принципиальной возможно-
стью образования неидеальных контактов по разным причинам: в связи с наличием реальных дефектов и эффекта сегрегации.
2. Задача гомогенизации. Для определения эффективного тензора диффузии необходимо решить задачу гомогенизации для гетерогенного материала. Процедура гомогенизации реализуется для репрезентативного объема V. Решается стационарная задача диффузии. Закон сохранения в отсутствие внутренних источников имеет следующий вид:
V I (х) = 0, (1)
где Л — диффузионный поток; х — радиус-вектор точки внутри V.
Предполагается, что поток в каждой точке репрезентативного объема определяется в соответствии с линейным законом Фика (влияние напряжений не учитывается):
Л (х) = -Б (х) V (х), (2)
где е — концентрация; Б — тензор диффузионной проницаемости материала (тензор диффузии примеси в данном материале).
Эффективный тензор диффузии Б6® связывает между собой средние по репрезентативному объему диффузионный поток и градиент концентрации, которые соответствуют значениям в точке сплошной среды на макроуровне:
Л = -Б6® -^е), (3)
где (...} = у /у (...) ¿V — осреднение по объему.
Одной из основных проблем при определении эффективных свойств материалов является идентификация микроструктурных параметров, через которые будут выражены искомые макроскопические характеристики. Эти параметры должны адекватным образом отражать вклад отдельных неоднородностей в оцениваемое свойство. В работе используется подход, разработанный И. Севостьяновым и М. Кача-новым [2], согласно которому роль микроструктурного параметра играет сумма так называемых тензоров вклада неоднородностей в искомое свойство, которая позволяет учесть влияние свойств и формы отдельных неоднородностей, а также их распределения по ориентациям. Учет взаимодействия неоднородностей при таком подходе реализуется отдельно в рамках самосогласованных схем. При этом наличие неидеальных контактов достаточно учесть только на этапе нахождения тензоров вклада неоднородностей. В этом случае можно ввести в рассмотрение эквивалентную неоднородность с идеальными контактами, дающую такой же вклад в искомое свойство, что и исходная. Поскольку целью настоящей работы является моделирование неидеальных контактов, сконцентрируемся на нахождении и анализе обобщенных на случай неидеальных контактов выражений для тензоров вклада неоднородностей, тогда как для определения эффективных свойств материала со множественными неоднородностями ограничимся построением границ области допустимых значений таких характеристик.
3. Тензор вклада неоднородности в диффузионную проницаемость.
Определим вклад отдельной неоднородности в искомое эффективное свойство. Рассмотрим репрезентативный объем V, состоящий из матрицы с диффузионной проницаемостью Бо и неоднородности объемом VI ^ V с проницаемостью Б1. На границе £ репрезентативного объема удобно задавать однородное условие Хилла
с (х)|Е = О0 • х, поскольку в таком случае при отсутствии внутренних источников среднее по объему значение градиента концентрации определяется вектором О0 [15]:
{Ус)у = О0. (4)
Средний поток, зависящий от микроструктуры, будет определяться следующим выражением [2]:
= + (5)
где Ис — тензор вклада неоднородности в диффузионную проницаемость (далее — тензор вклада неоднородности).
Тензор вклада эллипсоидальной неоднородности может быть найден на основании решения задачи Эшелби для изолированной неоднородности в бесконечной матрице. В этом случае он может быть выражен через тензор концентрации градиента концентрации Л (далее — тензор концентрации), связывающий постоянное поле внутри неоднородности с приложенным:
(УсК = Л • О0. (6)
Таким образом, для определения тензора вклада неоднородности необходимо, во-первых, решить задачу об осреденении полей и, во-вторых, найти тензор концентрации. Стоит отметить, что наличие неидеальных контактов, вообще говоря, должно быть учтено на обоих этапах. Тем не менее в литературе зачастую предполагается, что наличие таких контактов никак не влияет на этап осреднения, дополнительные проверки такого предположения не проводятся. Рассмотрим разные способы введения в модель неидеальных контактов и учтем их на обоих этапах.
3.1. Моделирование неидеального контакта посредством задания скачка поля. Неидеальный контакт может быть учтен в модели непосредственно с помощью задания скачка поля. В литературе рассмотрен случай задания скачка концентрации [6-9], моделирующий наличие сегрегации, тогда как схожий случай задания скачка нормальной компоненты потока, который может образовываться в результате растрескивания, не рассматривался. Исследуем оба случая.
Обратимся сперва к случаю неидеального контакта, при котором имеет место нарушение непрерывности потока, в частности образование дополнительных путей ускоренной/замедленной диффузии. Для моделирования такого контакта посредством задания скачка нормальной компоненты потока на границе раздела фаз Г изотропной матрицы (+) и изотропной неоднородности (—) поставим следующие условия:
дс (х)
•Лг+ — —> -Со" п
+ ' дпГ
дпГ
х^Г-
1Г
с (х)|х^Г+ = с (х)|х^Г- > (7)
где в^ играет роль параметра сегрегации для потока (термин введен по аналогии с известным термином «параметр сегрегации» в контексте описания оседания примеси).
Средний поток может быть определен через значения нормальной компоненты на границе Я репрезентативного объема V как
(Л) = -Л(х)хйЕ, (8)
где пт — внешняя нормаль. Здесь учтено, что в каждой точке внутри репрезентативного объема выполняется закон сохранения (1), а также что Ух = I.
Выражение (8) удобно переписать, дополняя поверхностный интеграл по внешней границе £ репрезентативного объема интегралом по поверхности неоднородности Г (добавляя и вычитая интегралы по поверхности Г от компонент потока .1п+ и 1п_) и используя теорему Остроградского— Гаусса:
(Л)
1
V
1
V
/ пт • Л (х) xd£ ^ ,1п+ xdГ ± 1п_ х^Г ит Jг Jг
(У пЕ • Лxd£ — ^ Зп+ xdГ^ + ^ Зп_ xdГ + ^ {1п+ - п) xdГ / ЛdVo + / ЛdVl + / (1п+ - П) xdГ
= (9)
1
V
где пг — внешняя нормаль к поверхности Г; V) = 1 — VI — объем матрицы.
Тогда с учетом скачка, задаваемого верхним условием (7), выражение (9) сведется к следующему:
VI
VI
V
(10)
где (.. .)уо = ± /уо (• • •) с1У0, (• • .)у! = Iу1 (• • •) ~ осреднение по матрице и по неоднородности соответственно.
Концентрация в рассматриваемом случае не испытывает скачка на границе раздела фаз, а потому выполняется следующее равенство:
= ( 1 - у) + у (Vc)Уl = С0.
(11)
Принимая во внимание, что и материал матрицы, и материал неоднородности подчиняются линейному закону Фика, получим из (10) с учетом (11)
(Л) = -Б0 • С0 - £ (8/Б! - Б0) • (Ус)У1 . С учетом равенства (6) выражение (12) сведется к
(Л) = —
Бо + у(5/Б1-Бо)-А/
Сс
(12)
(13)
где Лf — тензор концентрации для неоднородности с неидеальным контактом, моделируемым скачком нормальной компоненты потока.
Сравнивая выражения (5) и (13), получим следующее представление для тензора вклада неоднородности:
Н
о
(sfБ! — Бс) • А,.
(14)
Задача Эшелби для электропроводности в случае идеальных контактов на границе раздела фаз была решена Х. Фрике (Н. Епске) [16]. Следуя [16] и ставя условия
(7) на границе раздела матрица/неоднородность, получим следующее выражение
для тензора концентрации в случае эллипсоидальной неоднородности:
3 2
Л/= -ТТЪ-\eiei' (15)
i= 1 2 + Ljaia2a3 i- 1J где a и ej (г = 1, 2, 3) — полуоси и орты главных осей эллипсоида,
у Г , +
./0 V (a2 + s)(a2 + s)(a3 + s)
В случае сферической неоднородности выражение (15) сводится к
Л/ = --1, (16)
f 2D0 + sfD1 ' V ^
а тензор вклада неоднородности принимает вид
= 3 {sfD\ — Dp) = 3(Д/а-1)
2,00 + s/ßi 2 + S/a ' V ;
где а = Di/Do.
Перейдем к обсуждению случая неидеального контакта, при котором имеет место нарушение непрерывности концентрации, в частности оседание частиц диффундирующего вещества на границе раздела фаз. При моделировании неидеального контакта посредством задания скачка концентрации ставятся следующие условия [6, 7]:
дс{х)
дпГ
= D дс(х)
дп
х^Г-
С (Ж)1х^Г+ = всС , (18)
где вс — параметр сегрегации. Случай вс > 1 соответствует скоплению примеси снаружи неоднородности, а случай вс < 1 — оседанию внутри. Осреднение полей было осуществлено нами в работе [9]:
= (1 - У) + *су {Ус)У1 = С0. (20)
Тензор вклада неоднородности в рассматриваемом случае принимает вид [9]
= (Б1 - всБс) • Ас, (21)
где тензор концентрации Лс определяется следующим образом:
3 1 2 Лс = V---—-^е(22)
<=1 8°2 + Ьа1а2а3 (^-1)
Вестник СПбГУ. Математика. Механика. Астрономия. 2023. Т. 10(68). Вып. 4
Для сферической неоднородности справедливы следующие выражения для тензоров концентрации и вклада в диффузионную проницаемость соответственно:
Лс = 2 П^П1* (23)
2всОо + О
2 / _-у
в = ЦВ1-асВ0) = у
0 о „ п I п 0 о I а V/
3.2. Моделирование неидеального контакта посредством рассмотрения неоднородности с покрытием. Другим способом учета неидеального контакта на границе раздела фаз гетерогенного материала является рассмотрение неоднородности с представляющим собой суперпроводник или изолятор покрытием, толщина которого стремится к нулю [14, 17]. Рассмотрим оба случая.
Обратимся к задаче об осреднении полей для материла, состоящего из матрицы и неоднородности, состоящей в свою очередь из конфокальных эллипсоидов. Поскольку на границах раздела трех фаз имеют место только идеальные контакты, средний градиент концентрации может быть найден следующим образом:
(Vе) = (1 - у) (Vс)у0 + У ^с)уа + ^у^ {Ус)уь_уа = С0, (25)
где нижние индексы а и Ь относятся к внутреннему и внешнему эллипсоиду. Средний поток с учетом закона Фика и выражения (25) находится как
(Л) = _в0 • С0 - ^ (Б! - Б0) • (Ус)к - (Вя - Б0) • &с)уь_уа, (26)
где Б1, Бя — тензоры диффузионной проницаемости внутреннего эллипсоида и внешнего слоя неоднородности соответственно.
Отметим, что выражения (25) и (26) справедливы для обоих случаев моделирования неидеального контакта.
Поля во внутреннем эллипсоиде и во внешнем слое могут быть выражены через приложенное через тензоры концентрации Ла и Ая соответственно:
{Уе)Уа = Аа • С0, {Ус)уъ_уа = Ая • С0. (27)
В результате сравнения выражения (5) при VI = V и выражения (26) с учетом (27) получим следующее представление для тензора вклада неоднородности:
нл- 1 Н Й
Va (Б1 - Бо) • Аа + (V - Va) (бв - Бо) • Ля] . (28)
Заметим, что такое же выражение было получено в работе [14], где тензор вклада одиночной неоднородности сразу определялся через тензор концентрации по стандартной формуле (без явного учета неидеальных контактов).
Для описания поверхностного явления, связанного с наличием реального неидеального контакта, необходимо рассмотреть предельный переход, когда толщина слоя
6 стремится к нулю и эллипсоидальная неоднородность становится однородной. Тензоры концентрации и в окончательном виде тензоры вклада эллипсоидальных неод-нородностей были найдены В. Левиным и М. Марковым [14]. Так, согласно [14], тензор вклада неоднородности с проводящим покрытием имеет вид
где А = 6П8,
С = -^77, Н=-°-, с=*21°2+*214 + 4+4г Л = 1а1а2азМ< = 1,2,3)
4пС а\а2аз 2а\а2аз 2
F =
а1а2аз 3(s+a2 ) (s+a3 ) + (s+al ) (s+a2 ) + («+ai ) («+Дз ) ^
/ v^3 — 2 А л й1а2аз f°° s+fl9 ) ( s+а^ ) + ( s+a^ ) ( s+a, ) + ( s+a^ )( s+a
\2^k=1 ak ) A1--2— Je
^(S+af)5(S+al)3(S+al)3
2
здесь ^ — площадь поверхности эллипсоида. Выражения для £2 и £3 могут быть определены по последней формуле посредством замены индекса 1 на 2 и 3 соответственно.
В случае сферической неоднородности выражение для тензора вклада в диффузионную проницаемость сводится к
3^-1 + 2^
£>1+2Д, + 2£ « + 2 + 2Воа
= -и-- = £)g — -Т^-. (30)
Тензор вклада неоднородности с изолирующим покрытием имеет вид [14]
„D л ^_Dl-Do- DoDiRHGAj_
Н = Dn у ---ejej, (31)
¿i A^i + (1 - Ai) D0 + (l- Ai) DoDiRHG - f)
здесь R = -ß-.
В случае сферической неоднородности выражение (31) сводится к
п 3 - D0 - £>0£>i-) 3 (a — 1 — £>i—) H D=D0^ ° V=A) У n nтл1?¡,. (32)
aa
Отметим, что представленные в разделах 3.1 и 3.2 выражения для тензоров вклада неоднородностей могут в общем случае включать зависимость от температуры T как от параметра посредством рассмотрения зависимостей Do (T), Di (T), Ds (T).
3.3. Анализ выражений для тензоров вклада неоднородности в диффузионную проницаемость. Введем в рассмотрение так называемую эквивалентную неоднородность с идеальными контактами на границе раздела фаз, которая дает такой же вклад в макроскопические свойства материала, что и неоднородность с неидеальными контактами, и определим, какими диффузионными свойствами D* она должна обладать. Введение эквивалентной неоднородности позволяет существенно упростить анализ влияния неидеального контакта на эффективные свойства.
Обратимся сперва к случаю неидеального контакта, при котором имеет место нарушение непрерывности потока. Такой случай можно моделировать либо непосредственно заданием скачка нормальной компоненты потока, либо, согласно [17], наличием суперпроводящего слоя. В первом случае из равенства тензоров вклада исходной и эквивалентной неоднородности следует, что
Увеличение параметра в $ приводит к увеличению диффузионной проницаемости эквивалентной неоднородности. В случае = 1, соответствующем отсутствию скачка, тензоры диффузионной проницаемости эквивалентной и исходной неоднородности совпадают. В предельном случае, когда ^ ж>, эквивалентная неоднородность характеризуется бесконечной проницаемостью независимо от значения Ох. Величина а может принимать значение, равное единице, и тогда О*/Оо = а = 1, что в конечном счете означает, что свойства гомогенизированного материала будут совпадать со свойствами матрицы. То есть происходит «нивелирование» плохой/хорошей проницаемости неоднородности наличием скачка нормальной компоненты потока.
Можно показать, что при моделировании неидеального контакта посредством рассмотрения суперпроводящего покрытия у неоднородности из равенства тензоров вклада такой неоднородности и эквивалентной следует, что последняя характеризуется в общем случае ортотропным тензором диффузионной проницаемости. Анализ существенно упрощается в случае сферических неоднородностей. Тогда
Увеличение параметра Л приводит к увеличению диффузионной проницаемости эквивалентной неоднородности.
Сравнивая выражения (33) и (34), легко заметить, что два подхода к учету рассматриваемого типа неидеального контакта дают один и тот же результат при выполнении следующего равенства:
В случае Л = 0 (при отсутствии суперпроводящего покрытия), = 1 (отсутствует скачок нормальной компоненты потока). Увеличение Л приводит к увеличению скачка. Поскольку параметр Л > 0, параметр также должен быть больше или равен единице, т. е. величина потока со стороны матрицы может быть только больше потока со стороны неоднородности, что в случае диффузии соответствует случаю образования каналов ускоренной диффузии при растрескивании/отслоениях.
Возвращаясь к общему случаю эллипсоидальной неоднородности, заметим, что попытки вычисления параметра через величину Л приводят к разным результатам при сравнении разных компонент ортотропного тензора вклада эквивалентной неоднородности. Два подхода дают качественно разные результаты, интерпретация которых требует проведения отдельных исследований, выходящих за рамки работы.
Перейдем к обсуждению случая неидеального контакта, при котором имеет место нарушение непрерывности поля концентрации. Такой случай можно моделировать либо непосредственно заданием скачка концентрации, либо согласно [17]
Б* = Ох1.
(33)
(34)
в} = 1 + (2Л) / (Оха).
(35)
наличием изолирующего слоя. Тогда из равенства тензоров вклада исходной и эквивалентной неоднородностей следует, что в случае учета наличия неидеального контакта посредством задания скачка концентрации
D = (D1/sc) I. (36)
Увеличение параметра sc приводит к уменьшению диффузионной проницаемости эквивалентной неоднородности. В случае sc = 1, соответствующему отсутствию скачка, тензоры диффузионной проницаемости эквивалентной и реальной неоднородности совпадают. В случае sc > 1 диффузионная проницаемость эквивалентной неоднородности меньше соответствующей проводимости реальной неоднородности с неидеальным контактом, что соответствует физическим представлениям, поскольку наличие скачка приводит к тому, что диффузант проникает в неоднородность в меньшей степени. Соответственно, чтобы добиться такого же эффекта при рассмотрении неоднородности с идеальным контактом, необходимо уменьшить ее мас-сопроницаемость. В случае sc < 1 диффузионная проницаемость эквивалентной неоднородности больше соответствует проницаемости реальной неоднородности с неидеальным контактом, что также соответствует физике процесса, поскольку наличие скачка в этом случае ведет к увеличению концентрации внутри неоднородности, чего можно добиться посредством увеличения ее проницаемости. При a/sc = 1, D*/D0 = 1.
При моделировании неидеального контакта посредством рассмотрения изолирующего покрытия, эквивалентная неоднородность характеризуется в общем случае ортотропным тензором диффузии, компоненты которого зависят от формы неоднородности. Тензор становится изотропным в случае сферической неоднородности, рассмотрением которой далее и ограничимся. Тогда
(37)
Увеличение параметра R приводит к уменьшению диффузионной проницаемости эквивалентной неоднородности.
Сравнивая выражения (36) и (37), нетрудно убедиться, что два подхода к учету рассматриваемого типа неидеального контакта дают один и тот же результат, если
sc = 1 + (DiR/a). (38)
При R = 0 (в отсутствие изолирующего покрытия), sc = 1 (скачок концентрации отсутствует). Увеличение R приводит к увеличению скачка концентрации. Положительность R влечет за собой требование к величине параметра сегрегации: sc > 1. Такое условие накладывает ограничения на применимость моделей, а именно: они могут использоваться только при описании оседания примеси с внешней стороны границы раздела фаз, но не внутри неоднородностей. Таким образом, такое характерное для диффузионного процесса явление, как оседание диффузанта в порах необходимо моделировать либо иными способами, либо исключительно в рамках подхода, в котором скачок концентрации задается напрямую (см. раздел 3.1).
4. Обобщенные оценки Хашина — Штрикмана. Обобщим оценки Хаши-на — Штрикмана [2], представляющие собой верхнюю и нижнюю оценку значений
эффективных свойств изотропного на макроуровне материала. Заменив неоднородность с неидеальным контактом эквивалентной неоднородностью с идеальными контактами, преобразуем оценки следующим образом:
1 - Ф О6® ф
+ Ф ,1 ^ — 1-* , 1 (39)
Зsfa 1 — Sfa 3 Sfa—l
в случае материала с неидеальными контактами, при которых имеет место образование ускоренных путей диффузии, и
а 1 - ф О6® ф
<-- д се I _ се и д I се _
в случае материала с неидеальными контактами, при которых наблюдается скопление примеси на границе раздела фаз. Здесь ф — объемная доля неоднородностей.
В неравенствах (39), (40) для однозначности рассмотрен случай, когда О* < О0. В случае О* > О о выражения останутся прежними, но верхняя и нижняя оценки поменяются местами. Подчеркнем, что сравниваться между собой должны именно значения коэффициентов диффузии для матрицы и для эквивалентной, а не исходной, неоднородности, поскольку одновременно могут выполняться следующие условия: О* > Оо и Ох < Оо. В частности, именно такая ситуация может иметь место при моделировании неидеального контакта либо посредством задания скачка нормальной компоненты потока, либо посредством рассмотрения суперпроводящего покрытия (см. пояснения в тексте в разделе 3).
Аналогично посредством замены неоднородности с неидеальным контактом эквивалентной могут быть обобщены и различные методы гомогенизации.
5. Эффективная диффузионная проницаемость поликристалла. Проведем анализ влияния неидеальных контактов, моделируемых в рамках разных подходов, на эффективную диффузионную проницаемость поликристалла. Рассмотрим двухфазный материал, матрица которого представляет собой границы зерен, а изолированные неоднородности — зерна, при этом а = Ох/Оо < 1 (примем для определенности а = 0.2). В поликристаллах могут возникать неидеальные контакты обоих рассмотренных в статье типов. Так, за счет растрескивания могут образовываться дополнительные ускоренные пути диффузии, моделируемые либо посредством задания скачка нормальной компоненты потока, либо за счет рассмотрения супер-проводящего покрытия. Характерное для массопереноса явление оседания примеси в несплошностях, в частности, на границах зерен, может моделироваться либо посредством задания скачка концентрации, либо за счет рассмотрения изолирующего покрытия.
Обобщенные по формуле (39) оценки представлены на рис. 1. При отсутствии неидеальных контактов (в$ = 1) увеличение объемной доли неоднородностей ведет к уменьшению эффективной диффузионной проницаемости. В случае наличия неидеального контакта поведение функции зависит от величины а. Когда а < 1, эквивалентная неоднородность характеризуется коэффициентом диффузии О* < Оо, и ее присутствие приводит к уменьшению эффективных свойств (см. график при = 4). В случае а > 1, О* > Оо, в результате чего эффективные свойства гомогенизированного материала увеличиваются (см. графики при = 10, = 20).
яе70о
4
0
■ — Sf = 20
Sf = 10
Sf = 4
Sf = 1
■■■■•■■.■::::игг1111........
0.0
0.2
1.0
0.4 0.6 0.8 Ф
Рис. 1. Оценки Хашина — Штрикмана при наличии растрескивания.
0.4 0.6 Ф
Рис. 2. Оценки Хашина — Штрикмана при наличии сегрегации.
3
2
1
й^/й
о
Чем больше в^а отличается от единицы (т. е. чем больше О* отличается от Оо), тем шире границы области допустимых значений эффективных свойств.
Обобщенные по формуле (40) оценки представлены на рис. 2. Поскольку отношение а/вс всегда меньше единицы, О* всегда меньше О о, и присутствие такой эквивалентной неоднородности приводит только к уменьшению эффективной диффузионной проницаемости. Увеличение вс приводит к большему контрасту свойств эквивалентной неоднородности и матрицы и расширению границ области допустимых значений эффективных свойств.
6. Заключение. В работе обобщены способы учета неидеальных контактов при определении эффективных свойств материала. На примере материала со сферическими неоднородностями проведено сравнение двух подходов, предложенных в литературе в контексте описания разных физических процессов: подхода, в рамках которого напрямую задается скачок поля, и подхода, в котором в рассмотрение вводится неоднородность с покрытием, обладающим экстремальными свойствами. Проведена интерпретация полученных результатов на примере задачи диффузии, в которой неидеальные контакты могут возникать по разным причинам: как в связи с наличием растрескивания, приводящего к появлению обходных ускоренных путей
диффузии, так и вследствие характерного для массопереноса эффекта сегрегации, заключающегося в оседании примеси на границе раздела фаз. Показана ограниченность применимости подходов к описанию эффекта оседания диффузанта внутри неоднородностей.
Литература
1. Kanaun S.K., Levin V. M. Effective field method in mechanics of matrix composite materials. Advances in mathematical modelling of composite materials 1 (58) (1994). https://doi.org/10.1142/9789814354219_0001
2. Kachanov M., Sevostianov I. Micromechanics of materials, with applications. Cham, Springer (2018).
3. Zhang Y., Liu L. On diffusion in heterogeneous media American Journal of Science 312 (9), 1028-1047 (2012). https://doi.org/10.2475/09.2012.03
4. Бокштейн Б. С., Магидсон И. А., Светлов И. Л. О диффузии в объеме и по границам зерен. Физика металлов и металловедение 6 (6), 1040-1052 (1958).
5. Kaur I., Mishin Y., Gust W. Fundamentals of grain and interphase boundary diffusion. John Wiley (1995).
6. Kalnin J.R., Kotomin E.A., Maier J. Calculations of the effective diffusion coefficient for inhomogeneous media. Journal of physics and chemistry of solids 63 (3), 449-456 (2002). https://doi.org/10.1016/S0022-3697(01)00159-7
7. Belova I. V., Murch G.E. Calculation of the effective conductivity and diffusivity in composite solid electrolytes. Journal of Physics and Chemistry of Solids 66 (5), 722-728 (2005). https://doi.org/10.1016/j .jpcs.2004.09.009
8. Knyazeva A. G., Grabovetskaya G. P., Mishin I. P., Sevostianov I. On the micromechanical modelling of the effective diffusion coefficient of a polycrystalline material. Philosophical Magazine 95 (19), 2046-2066 (2015). https://doi.org/10.1080/14786435.2015.1046965
9. Frolova K. P., Vilchevskaya E.N. Effective diffusivity of transversely isotropic material with embedded pores. Materials Physics & Mechanics 47 (6), 937-950 (2021). https://doi.org/10.18149 /MPM.4762021_12
10. Miloh T., Benveniste Y. On the effective conductivity of composites with ellipsoidal inhomogeneities and highly conducting interfaces. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 455 (1987), 2687-2706 (1999). https://doi.org/10.1098/rspa.1999.0422
11. Duan H.L., Karihaloo B. L. Effective thermal conductivities of heterogeneous media containing multiple imperfectly bonded inclusions. Physical Review B. 75 (6), 064206 (2007). https://doi.org/10.1103/PhysRevB.75.064206
12. Kushch V. I., Sevostianov I., Belyaev A. S. Effective conductivity of spheroidal particle composite with imperfect interfaces: Complete solutions for periodic and random micro structures. Mechanics of Materials 89, 1-11 (2015). https://doi.org/10.1016Zj.mechmat.2015.05.010
13. Endres A. L., Knight R. J. A model for incorporating surface phenomena into the dielectric response of a heterogeneous medium. Journal of colloid and interface science 157(2), 418-425 (1993). https://doi.org/10.1006/jcis.1993.1204
14. Levin V., Markov M. Effective thermal conductivity of micro-inhomogeneous media containing imperfectly bonded ellipsoidal inclusions. International Journal of Engineering Science 109, 202-215 (2016). https://doi.org/10.1016/j.ijengsci.2016.09.012
15. Hill R. Elastic properties of reinforced solids: some theoretical principles. Journal of the Mechanics and Physics of Solids 11 (5), 357-372 (1963). https://doi.org/10.1016/0022-5096(63)90036-X
16. Fricke H. Mathematical treatment of the electric conductivity and capacity of disperse systems I. The electric conductivity of a suspension of homogeneous spheroids. Physical Review 24 (5), 575 (1924). https://doi.org/10.1103/physRev.24.575
17. Markov K. Z. Elementary micromechanics of heterogeneous media. Birkhauser, Boston (2000).
Статья поступила в редакцию 23 марта 2023 г.;
доработана 17 мая 2023 г.; рекомендована к печати 18 мая 2023 г.
Контактная информация:
Фролова Ксения Петровна — канд. физ.-мат. наук, ст. науч. сотр.; fkp@ipme.ru Вильчевская Елена Никитична — д-р физ.-мат. наук, доц., вед. науч. сотр.; ven@ipme.ru Полянский Владимир Анатольевич — д-р техн. наук; vapol@mail.ru
Modeling of imperfect contacts in determining the effective diffusion permeability
K. P. Frolova, E.N. Vilchevskaya, V. A. Polyanskiy
Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences, 61, Bolshoy pr. V. O., St. Petersburg, 199178, Russian Federation
For citation: Frolova K.P., Vilchevskaya E. N., Polyanskiy V. A. Modeling of imperfect contacts in determining the effective diffusion permeability. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2023, vol. 10(68), issue 4, pp. 650-664. https://doi.org/10.21638/spbu01.2023.405 (In Russian)
The work develops a universal approach to acounting for imperfect contacts in determining the effective properties of various nature, namely, effective diffusivity, thermal and electrical conductivity. Imperfect contacts take place when fields at the microlevel are not continuous. The possibility of creating a unified approach is due to the similarity of the governing equations. At the same time, the appearance of imperfect contacts can be caused by microstructural features and by the specifics of the process itself. For concreteness, the effective diffusion permeability is determined, since various reasons for the appearance of imperfect contacts can be considered. The reasons can be associated both with the formation of structural defects and with the presence of the specific segregation effect. The paper generalizes and compares two approaches to accounting for imperfect contacts. In the first case, a field jump is set. In the second case, an inhomogeneity with a thin coating possessing extreme properties is introduced. A comprehensive analysis is carried out on the example of a material with spherical inhomogeneities. Analytical expressions for contribution tensor of the equivalent inhomogeneity are obtained, which results in simplification of generalization of various homogenization methods.
Keywords: effective properties, inhomogeneity, imperfect contacts, segregation, diffusivity. References
1. Kanaun S.K., Levin V. M. Effective field method in mechanics of matrix composite materials. Advances in mathematical modelling of composite materials 1 (58) (1994). https://doi.org/10.1142/9789814354219_0001
2. Kachanov M., Sevostianov I. Micromechanics of materials, with applications. Cham, Springer (2018).
3. Zhang Y., Liu L. On diffusion in heterogeneous media American Journal of ¡Science 312 (9), 1028-1047 (2012). https://doi.org/10.2475/09.2012.03
4. Bokstein B. S., Magidson I. A., Svetlov I. L. On diffusion in volume and along grain boundaries. Physics of metals and metallography 6 (6), 1040—1052 (1958). (In Russian)
5. Kaur I., Mishin Y., Gust W. Fundamentals of grain and interphase boundary diffusion. John Wiley (1995).
6. Kalnin J.R., Kotomin E.A., Maier J. Calculations of the effective diffusion coefficient for inhomogeneous media. Journal of physics and chemistry of solids 63 (3), 449—456 (2002). https://doi.org/10.1016/S0022-3697(01)00159-7
7. Belova I. V., Murch G.E. Calculation of the effective conductivity and diffusivity in composite solid electrolytes. Journal of Physics and Chemistry of Solids 66 (5), 722—728 (2005). https://doi.org/10.1016/j .jpcs.2004.09.009
8. Knyazeva A. G., Grabovetskaya G. P., Mishin I. P., Sevostianov I. On the micromechanical modelling of the effective diffusion coefficient of a polycrystalline material. Philosophical Magazine 95 (19), 2046-2066 (2015). https://doi.org/10.1080/14786435.2015.1046965
9. Frolova K. P., Vilchevskaya E.N. Effective diffusivity of transversely isotropic material with embedded pores. Materials Physics & Mechanics 47 (6), 937-950 (2021). https://doi.org/10.18149 /MPM.4762021_12
10. Miloh T., Benveniste Y. On the effective conductivity of composites with ellipsoidal inhomogeneities and highly conducting interfaces. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 455 (1987), 2687-2706 (1999). https://doi.org/10.1098/rspa.1999.0422
11. Duan H.L., Karihaloo B. L. Effective thermal conductivities of heterogeneous media containing multiple imperfectly bonded inclusions. Physical Review B. 75 (6), 064206 (2007). https://doi.org/10.1103/PhysRevB.75.064206
12. Kushch V. I., Sevostianov I., Belyaev A. S. Effective conductivity of spheroidal particle composite with imperfect interfaces: Complete solutions for periodic and random micro structures. Mechanics of Materials 89, 1-11 (2015). https://doi.org/10.1016/j.mechmat.2015.05.010
13. Endres A.L., Knight R.J. A model for incorporating surface phenomena into the dielectric response of a heterogeneous medium. Journal of colloid and interface science 157(2), 418-425 (1993). https://doi.org/10.1006/jcis.1993.1204
14. Levin V., Markov M. Effective thermal conductivity of micro-inhomogeneous media containing imperfectly bonded ellipsoidal inclusions. International Journal of Engineering Science 109, 202-215 (2016). https://doi.org/10.1016/j.ijengsci.2016.09.012
15. Hill R. Elastic properties of reinforced solids: some theoretical principles. Journal of the Mechanics and Physics of Solids 11 (5), 357-372 (1963). https://doi.org/10.1016/0022-5096(63)90036-X
16. Fricke H. A mathematical treatment of the electric conductivity and capacity of disperse systems I. The electric conductivity of a suspension of homogeneous spheroids. Physical Review 24 (5) 575 (1924). https://doi.org/10.1103/PhysRev.24.575
17. Markov K. Z. Elementary micromechanics of heterogeneous media. Birkhauser, Boston (2000).
Received: March 23, 2023 Revised: May 17, 2023 Accepted: May 18, 2023
Authors' information:
Ksenia P. Frolova — fkp@ipme.ru Elena N. Vilchevskaya — ven@ipme.ru Vladimir A. Polyanskiy — vapol@mail.ru