Статьи

 




dx.doi.org/ 10.18577/2307-6046-2017-0-6-9-9
УДК 66.017
С. В. Кондрашов, К. А. Шашкеев, О. В. Попков
СПОСОБ ВЫЧИСЛЕНИЯ ЭФФЕКТИВНЫХ ПАРАМЕТРОВ ПЕРИОДИЧЕСКОЙ СРЕДЫ

Предложен способ вычисления эффективных электродинамических параметров (тензоров восприимчивостей) частично разупорядоченной периодической среды с прямоугольной ячейкой, размеры которой существенно меньше длины волны. Предложенный итеративный алгоритм имеет вычислительную сложность итерации O(NxNyNzlnNxlnNylnNz), а его потребность в объеме памяти ЭВМ составляет Θ(NxNyNz), где Nx, Ny и Nz– размер сетки по соответствующим направлениям. Алгоритм может быть расширен на случай кристаллов с произвольной геометрической формой ячейки путем введения дополнительных линейных преобразований.

Ключевые слова: эффективные электромагнитные параметры, диэлектрическая проницаемость, магнитная проницаемость, периодическая среда, разупорядоченная среда, параллельные вычисления, effective electromagnetic parameters, permittivity, permeability, periodic medium, disordered medium, parallel computations.

Введение

Уровень современной авиационной техники в большой степени определяется уровнем радиоэлектронного оборудования, а именно – радиоэлектронной части летательного аппарата: систем, отвечающих за интерфейс «человек–машина», пространственное положение и пилотирование, средств телекоммуникации, энергообеспечения, безопасности, управления, контроля и обеспечения жизненного цикла всех агрегатов и систем летательного аппарата и др. В последние десятилетия электроника, в том числе и авиационная, именуемая авионикой, является одной из наиболее динамично развивающихся отраслей человеческой деятельности, в особенности – применительно к созданию авиационных комплексов с беспилотными летательными аппаратами [1, 2].

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

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

Одним из наиболее перспективных радиопоглощающих наполнителей для решения поставленной задачи являются углеродные нанотрубки (УНТ) [4–9]. К преимуществам УНТ относятся низкий порог проводимости и возможность достижения равномерного ориентационного распределения нанотрубок внутри тонких слоев композита благодаря их малым размерам.

В ряде работ, например [10–13], определены частотные зависимости действительной и мнимой частей диэлектрической проницаемости различных нанокомпозитов с углеродными нанотрубками. Характерной особенностью всех исследованных композитов является резкое изменение величин действительной и мнимой частей диэлектрической проницаемости (ε′ и εʺ) при увеличении проводимости после достижения порога перколяции. Так, авторы работы [13] показали, что при изменении концентрации одностенных УНТ с 0,1 до 1% величина ε′ изменяется от 4 до 8. При этом величина тангенса угла диэлектрических потерь (tgδ) изменяется с 0,1 до 0,24, а проводимость нанокомпозита увеличивается с 10-12 до 10-9 См/см. Увеличение концентрации УНТ до 3% (проводимость нанокомпозита составляет 10-3 См/см) приводит к возрастанию величины ε′ до 16, а tgδ – до 0,7.

Использование плоских пластин из композитов с углеродными нанотрубками позволяет получить узкополосный радиопоглощающий материал. Положение максимума поглощения и его ширину можно регулировать путем изменения концентрации УНТ. Например, согласно результатам работы [14], максимум поглощения и ширина пика поглощения по уровню -10 дБ изменяется с 6,4–8,2 (1,8 ГГц) до 7,5–10,1 (2,6 ГГц) и 12,0–15,1 (3,1 ГГц) при изменении концентрации одностенных УНТ в полиуретане. На ширину и положение пика максимума поглощения оказывает влияние тип матрицы [15] нанокомпозита, качество диспергирования УНТ [16] и другие факторы, которые влияют на проводимость нанокомпозита.

Возможным способом увеличения ширины полосы поглощения нанокомпозитов с углеродными нанотрубками является интегрирование наночастиц ферромагнитного металла или оксида внутрь углеродной нанотрубки. Образующийся карбид со структурой «ядро–оболочка» наряду с диэлектрическими потерями обладает высоким уровнем магнитных потерь. Такая комбинация позволяет получать нанокомпозиты с шириной полосы поглощения 16 ГГц по уровню -10 дБ [17] или даже 11 ГГц – по уровню -20 дБ [18].

Одна из основных проблем, которую необходимо решить при разработке радиопоглощающего композита, заключается в построении модели, описывающей диэлектрические и магнитные свойства среды с разупорядоченным распределением нанонаполнителя.

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

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

Работа выполнена в рамках реализации комплексного научного направления 15. «Наноструктурированные, аморфные материалы и покрытия» («Стратегические направления развития материалов и технологий их переработки на период до 2030 года») [2].

 

Теоретическая часть

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

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

– в этом случае понятие эффективных параметров среды остается актуальным независимо от строения ячейки;

– распространение электромагнитного излучения от близко расположенных (оказывающих существенное влияние) ячеек к рассматриваемой ячейке можно считать мгновенным (в противном случае в рамках предложенного метода возникают сложности при вычислении);

– это позволяет считать внешнее электромагнитное поле однородным.

Предлагаемая схема расчета основывается на итеративном вычислении полей поляризации и намагниченности периодической среды в заданном внешнем однородном монохроматическом электромагнитном поле. На каждой итерации сначала вычисляются поправки к полям поляризации и намагниченности, а затем – вызванные ими поправки к электрическому и магнитному полям.

Изначально задается однородное внешнее поле  с круговой частотой ω (где    – комплекснозначные векторы напряженности электрического поля и индукции магнитного поля), а также комплекснозначные параметры среды, устанавливающие линейную зависимость между параметрами   (где    – комплекснозначные векторы поляризации и намагниченности). Дальнейший расчет строится на следующих рекуррентных соотношениях:

где  – комплекснозначный тензор второго ранга, описывающий восприимчивость величины α к величине β в линейном приближении (α и β – индексы, указывающие на тип связанных величин – см. выражения (7) и (8)); φ – скалярный электрический потенциал;  – векторный потенциал; , k0=ω/c, c – скорость света в вакууме; i=√-1; j – номер итерации; векторы со штрихом – найденная на текущей итерации поправка к полю;    – радиус-векторы; r – абсолютное значение радиус-вектора;  – элемент объема, по которому производится интегрирование.

 

Решением данной системы уравнений является поле 

Зададим прямое и обратное преобразования Фурье следующим образом:

                        

где  – Фурье-образ некоторой функции F.

 

Дискретное преобразование Фурье (ДПФ) зададим в следующем виде:

                            

 

Избавимся от «сверток» и дифференцирования в выражениях (1)–(12), переписав их с привлечением операций в обратном пространстве:

                  

Рассмотрим кристаллическую ячейку, имеющую форму параллелепипеда со сторонами ах, ау, аz. Введем функцию f′, которая учитывает вклад от одинаковых элементов всех таких ячеек:

                    

где  – вектор смещения рассматриваемой ячейки,

тогда

                     

 

Можно показать, что функция  отлична от нуля только при kαaα=2πN, где α – индекс компоненты вектора; aα – параметр решетки в направлении a, . Переходя к дискретному представлению ряда Фурье, можно записать:

                         

где  V– объем ячейки.

 

Согласно выражениям (17) и (18), достаточно рассматривать параметры  и  с теми же  , для которых отличаются от нуля коэффициенты Фурье-разложения полей поляризации и намагниченности. Далее будем считать, что поля целиком описываются набором синусоид с пространственными частотами, не превышающими частот Найквиста по соответствующим направлениям ячейки.

Коэффициенты Фурье-разложения   для k=0 не влияют на векторы   и на любом этапе вычислений могут быть приняты равными нулю, благодаря чему можно также принять  . В случае k0<<k выражения (21) и (22) могут быть упрощены:

                      

 

При этом условии также упрощаются выражения (17)–(20):

 

             

Перепишем уравнение (29) для напряженности магнитного поля:

                     

Таким образом, в случае нулевых параметров    оказывается возможным произвести разделение задачи на независимые подзадачи для электрического и магнитного полей.

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

Как следует из вычислительной сложности быстрого преобразования Фурье (БПФ), предложенный итеративный алгоритм имеет вычислительную сложность итерации О(NxNyNzlnNxlnNylnNz), а его потребность в объеме памяти ЭВМ составляет Θ(NxNyNz), где Nx, Ny и Nz – размер сетки по соответствующим направлениям. Таким образом, предложенный способ так же хорошо масштабируется, как и БПФ.

 

Выводы

Предложенный способ расчета эффективных электродинамических параметров (тензоров восприимчивостей) периодической среды (кристалла с ячейкой в виде параллелепипеда) обладает хорошей масштабируемостью по отношению к размеру сетки.

В случае равенства нулю тензоров   задача может быть разделена на независимые подзадачи для электрического и магнитного полей, что позволяет уменьшить потребность в объеме памяти ЭВМ и ускорить расчет.

 

Данная работа выполнена при финансовой поддержке РФФИ в рамках научного проекта №15-38-20717 мол_а_вед.

 


ЛИТЕРАТУРА REFERENCE LIST
1. Каблов Е.Н. Конструкционные и функциональные материалы – основа экономического и научно-технического развития России // Вопросы материаловедения. 2006. №1. С. 64–67.
2. Каблов Е.Н. Инновационные разработки ФГУП «ВИАМ» ГНЦ РФ по реализации «Стратегических направлений развития материалов и технологий их переработки на период до 2030 года» // Авиационные материалы и технологии. 2015. №1 (34) С. 3–33. DOI:10.18577/2071-9140-2015-0-1-3-33.
3. Кондрашов С.В., Гуревич Я.М., Попков О.В., Шашкеев К.А. и др. Широкополосные радиопоглощающие материалы на основе пористых композитов с углеродными нанотрубками // Все материалы: Энциклопедический справочник (в печати).
4. Qinand F., Brosseau C. A review and analysis of microwave absorption in polymer composites filled with carbonaceous particles // J. Appl. Phys. 2012. Vol. 111. Paper No. 061301. DOI: 10.1063/1.3688435.
5. Кондрашов С.В., Шашкеев К.А., Попков О.В., Соловьянчик Л.В. Перспективные технологии получения функциональных материалов конструкционного назначения на основе нанокомпозитов с УНТ (обзор) // Труды ВИАМ: электрон. науч.-технич. журн. 2016. №3. Ст. 07. URL: http://www.viam-works.ru (дата обращения: 05.07.2016). DOI: 10.18577/2307-6046-2016-0-3-7-7.
6. Кондрашов С.В., Шашкеев К.А., Попков О.В., Соловьянчик Л.В. Физико-механические свойства нанокомпозитов с УНТ (обзор) // Труды ВИАМ: электрон. науч.-технич. журн. 2016. №5. Ст. 08. URL: http://www.viam-works.ru (дата обращения: 05.07.2016). DOI: 10.18577/2307-6046-2016-0-5-8-8.
7. Каблов Е.Н., Кондрашов С.В., Юрков Г.Ю. Перспективы использования углеродсодержащих наночастиц в связующих для полимерных композиционных материалов // Российские нанотехнологии. 2013. Т. 8. №3–4. С. 28–42.
8. Кондрашов С.В., Гуняева А.Г., Шашкеев К.А. и др. Электропроводящие гибридные полимерные композиционные материалы на основе нековалентно функционализированных углеродных нанотрубок // Труды ВИАМ: электрон. науч.-технич. журн. 2016. №2. Ст. 10. URL: http://www.viam-works.ru (дата обращения: 05.07.2016). DOI: 10.18577/2307-6046-2016-0-2-10-10.
9. Акатенков Р.В., Аношкин И.В., Беляев А.А., Битт В.В., Богатов В.А., Дьячкова Т.П., Куцевич К.Е., Кондрашов С.В. и др. Влияние структурной организации углеродных нанотрубок на радиоэкранирующие и электропроводящие свойства нанокомпозитов // Авиационные материалы и технологии. 2011. №1. С. 35–42.
10. Challa R.K., Kajfez D., Demir V., Gladden J.R., Elsherbeni A.Z. Characterization of Multiwalled Carbon Nanotube (MWCNT) Composites in a of Square Cross Section Waveguide // IEEE Microwave and Wireless Components Letters. 2008. Vol. 18. No. 3. P. 161–163. http://ieeexplore.ieee.org/document/4459260/ (дата обращения: 23.01.2017). DOI: 10.1109/LMWC.2008.916776.
11. De Vivo B., Guagno L., Lamberti P. et al. Electromagnetic properties of Carbon NanoTube/epoxy nanocomposites // International Symposium on Electromagnetic Compatibility – EMC Europe. Athens, 2009. P. 1–4. http://ieeexplore.ieee.org/document/5189674/ (дата обращения: 23.01.2017). DOI: 10.1109/EMCEUROPE.2009.5189674.
12. Huang Q., Holland T.B., Mukherjee A.K. et al. Carbon nanotube RF absorbing materials // Proc. SRF2009. Berlin, 2009. P. 648–651.
13. Huang Yi., Li N., Ma Y. et al. The influence of single-walled carbon nanotube structure on the electromagnetic interference shielding efficiency of its epoxy composites // Carbon. 2007. Vol. 45. P. 1614–1621. DOI: 10.1016/j.carbon.2007.04.016.
14. Liu Z.F., Bai G., Huang Y. et al. Microwave Absorption of Single-Walled Carbon Nanotubes/Soluble Cross-Linked Polyurethane Composites // J. Phys. Chem. C. 2007. Vol. 111. No. 37. P. 13696–13700. DOI: 10.1021/jp0731396.
15. Fan Z., Luo G., Zhang Z. et al. Electromagnetic and microwave absorbing properties of multi-walled carbon nanotubes/polymer composites // Materials Science and Engineering B. 2006. Vol. 132. P. 85–89.
16. Zhang Y., Zhang A., Ding L. et al. The effect of polymer spatial configuration on the microwave absorbing properties of non-covalent modified MWNTs // Composites. Part A. 2016. Vol. 81. P. 264–270.
17. Che R.C., Liang C.Y., Shi H.L. et al. Electron energy-loss spectroscopy characterization and microwave absorption of iron-filled carbon-nitrogen nanotubes// Nanotechnology. Vol. 18. No. 35. Paper No. 355705.
18. Qi X., Hu Q., Xu J. et al. Enhanced microwave absorption properties and mechanism of core/shell structured magnetic nanoparticles/carbon-based nanohybrids // Materials Science and Engineering B. 2016. Vol. 211. P. 53–60.
1. Kablov E.N. Konstrukcionnye i funkcionalnye materialy – osnova ekonomicheskogo i nauchno-tehnicheskogo razvitiya Rossii [Constructional and functional materials are the base of economic and scientific and technical development of Russia] // Voprosy materialovedeniya. 2006. №1. S. 64–67.
2. Kablov E.N. Innovacionnye razrabotki FGUP «VIAM» GNC RF po realizacii «Strategicheskih napravlenij razvitiya materialov i tehnologij ih pererabotki na period do 2030 goda» [Innovative developments of FSUE «VIAM» SSC of RF on realization of «Strategic directions of the development of materials and technologies of their processing for the period until 2030»] // Aviacionnye materialy i tehnologii. 2015. №1 (34). S. 3–33. DOI: 10.18577/2071-9140-2015-0-1-3-33.
3. Kondrashov S.V., Gurevich Ya.M., Popkov O.V., Shashkeev K.A. i dr. Shirokopolosnye radiopogloshhayushhie materialy na osnove poristyh kompozitov s uglerodnymi nanotrubkami [Broadband radio absorbing materials on the basis of porous composites with carbon nanotubes] // Vse materialy: Enciklopedicheskij spravochnik (v pechati).
4. Qinand F., Brosseau C. A review and analysis of microwave absorption in polymer composites filled with carbonaceous particles // J. Appl. Phys. 2012. Vol. 111. Paper No. 061301. DOI: 10.1063/1.3688435.
5. Kondrashov S.V., Shashkeev K.A., Popkov O.V., Solovyanchik L.V. Perspektivnye tehnologii polucheniya funkcionalnyh materialov konstrukcionnogo naznacheniya na osnove nanokompozitov s UNT (obzor) [Prospective producing methods for functional structural materials based on CNT-filled nanocomposites (review)] // Trudy VIAM: elektron. nauch.-tehnich. zhurn. 2016. №3. St. 07. Available at: http://www.viam-works.ru (accessed: July 05, 2016). DOI: 10.18577/2307-6046-2016-0-3-7-7.
6. Kondrashov S.V., Shashkeev K.A., Popkov O.V., Solovyanchik L.V. Fiziko-mehanicheskie svojstva nanokompozitov s UNT (obzor) [Mechanical properties of CNT nanocomposites (review)] // Trudy VIAM: elektron. nauch.-tehnich. zhurn. 2016. №5. St. 08. Available at: http://www.viam-works.ru (accessed: July 05, 2016). DOI: 10.18577/2307-6046-2016-0-5-8-8.
7. Kablov E.N., Kondrashov S.V., Yurkov G.Yu. Perspektivy ispolzovaniya uglerodsoderzhashhih nanochastic v svyazuyushhih dlya polimernyh kompozicionnyh materialov [Perspectives of use of carbon-containing nanoparticles in binding for polymeric composite materials] // Rossijskie nanotehnologii. 2013. T. 8. №3–4. S. 28–42.
8. Kondrashov S.V., Gunyaeva A.G., Shashkeev K.A., Barinov D.Ya., Soldatov M.A., Shevchenko V.G., Muzafarov A.M. Elektroprovodyashhie gibridnye polimernye kompozicionnye materialy na osnove nekovalentno funkcionalizirovannyh uglerodnyh nanotrubok [Electrically-conductive hybrid polymer composite materials on the basis of noncovalent functional carbon nanotubes] // Trudy VIAM: elektron. nauch.-tehnich. zhurn. 2016. №2. St. 10. Available at: http://www.viam-works.ru (accessed: July 05, 2016). DOI: 10.18577/2307-6046-2016-0-2-10-10.
9. Akatenkov R.V., Anoshkin I.V., Belyaev A.A., Bitt V.V., Bogatov V.A., Dyachkova T.P., Kutsevich K.E., Kondrashov S.V., Romanov A.M., Shirokov V.V., Horobrov N.V. Vliyanie strukturnoj organizacii uglerodnyh nanotrubok na radiojekraniruyushhie i elektroprovodyashhie svojstva nanokompozitov [Influence of the structural organization of carbon nanotubes on radio shielding and electrocarrying-out properties of nanocomposites] // Aviacionnye materialy i tehnologii. 2011. №1. S. 35–42.
10. Challa R.K., Kajfez D., Demir V., Gladden J.R., Elsherbeni A.Z. Characterization of Multiwalled Carbon Nanotube (MWCNT) Composites in a of Square Cross Section Waveguide // IEEE Microwave and Wireless Components Letters. 2008. Vol. 18. No. 3. P. 161–163. http://ieeexplore.ieee.org/document/4459260/ (accessed: December 23, 2017). DOI: 10.1109/LMWC.2008.916776.
11. De Vivo B., Guagno L., Lamberti P. et al. Electromagnetic properties of Carbon NanoTube/epoxy nanocomposites // International Symposium on Electromagnetic Compatibility – EMC Europe. Athens, 2009. P. 1–4. http://ieeexplore.ieee.org/document/5189674/ (accessed: December 23, 2017). DOI: 10.1109/EMCEUROPE.2009.5189674.
12. Huang Q., Holland T.B., Mukherjee A.K. et al. Carbon nanotube RF absorbing materials // Proc. SRF2009. Berlin, 2009. P. 648–651.
13. Huang Yi., Li N., Ma Y. et al. The influence of single-walled carbon nanotube structure on the electromagnetic interference shielding efficiency of its epoxy composites // Carbon. 2007. Vol. 45. P. 1614–1621. DOI: 10.1016/j.carbon.2007.04.016.
14. Liu Z.F., Bai G., Huang Y. et al. Microwave Absorption of Single-Walled Carbon Nanotubes/Soluble Cross-Linked Polyurethane Composites // J. Phys. Chem. C. 2007. Vol. 111. No. 37. P. 13696–13700. DOI: 10.1021/jp0731396.
15. Fan Z., Luo G., Zhang Z. et al. Electromagnetic and microwave absorbing properties of multi-walled carbon nanotubes/polymer composites // Materials Science and Engineering B. 2006. Vol. 132. P. 85–89.
16. Zhang Y., Zhang A., Ding L. et al. The effect of polymer spatial configuration on the microwave absorbing properties of non-covalent modified MWNTs // Composites. Part A. 2016. Vol. 81. P. 264–270.
17. Che R.C., Liang C.Y., Shi H.L. et al. Electron energy-loss spectroscopy characterization and microwave absorption of iron-filled carbon-nitrogen nanotubes// Nanotechnology. Vol. 18. No. 35. Paper No. 355705.
18. Qi X., Hu Q., Xu J. et al. Enhanced microwave absorption properties and mechanism of core/shell structured magnetic nanoparticles/carbon-based nanohybrids // Materials Science and Engineering B. 2016. Vol. 211. P. 53–60.
Вы можете оставить комментарий к статье. Для этого необходимо зарегистрироваться на сайте.