Электродинамический расчет комплексной эффективной диэлектрической проницаемости нанокомпозитов на основе массивов углеродных нанотрубок в диапазоне сверхвысоких частот Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 535.32

О. А. Голованов, Г. С. Макеева, Д. Н. Ширшиков, Г. Г. Горлов

ЭЛЕКТРОДИНАМИЧЕСКИЙ РАСЧЕТ КОМПЛЕКСНОЙ ЭФФЕКТИВНОЙ ДИЭЛЕКТРИЧЕСКОЙ ПРОНИЦАЕМОСТИ НАНОКОМПОЗИТОВ НА ОСНОВЕ МАССИВОВ УГЛЕРОДНЫХ НАНОТРУБОК В ДИАПАЗОНЕ СВЕРХВЫСОКИХ ЧАСТОТ

Аннотация.

Актуальность и цели. Наноструктурные материалы на основе углеродных нанотрубок (УНТ) благодаря электронным свойствам УНТ играют важную роль в создании материалов для поглощения и экранирования электромагнитного излучения СВЧ. Целью данной работы является теоретическое исследование взаимодействия электромагнитных волн с массивами ориентированных УНТ и расчет эффективной диэлектрической проницаемости анизотропных нанокомпозитов на УНТ в СВЧ-диапазоне на основе математического моделирования электродинамического уровня строгости.

Материалы и методы. Методом автономных блоков с каналами Флоке проведено математическое моделирование распространения электромагнитных волн в периодических 3D-решетках ориентированных УНТ в диапазоне СВЧ. Вычислительный алгоритм решения краевой 3D-задачи дифракции для определения матрицы проводимости автономных блоков с каналами Флоке в виде прямоугольного параллелепипеда, содержащего УНТ, построен проекционным методом Галеркина. Разработана методика электродинамического расчета эффективной диэлектрической проницаемости нанокомпозита на основе УНТ, базирующаяся на теории эффективной среды.

Результаты. Получены результаты электродинамического расчета частотных зависимостей действительной и мнимой частей комплексной диэлектрической проницаемости нанокомпозита на основе УНТ в диапазоне частот 2-5 ГГц для различной ориентации вектора электрического поля Е волны к оси УНТ при различном соотношении объема УНТ к общему объему композита.

Выводы. Показано, что степень ослабления СВЧ-излучения за счет поглощения массивом УНТ, существенно возрастает при увеличении проводимости УНТ, фактора заполнения композита УНТ, изменяется в зависимости от взаимной ориентации вектора электрического поля Е волны и осей УНТ и при параллельной ориентации Е в сравнении с ортогональной увеличивается.

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

O. A. Golovanov, G. S. Makeeva, D. N. Shirshikov, G. G. Gorlov

ELECTRODYNAMIC CALCULATION OF THE COMPLEX EFFECTIVE PERMEABILITY OF NANOCOMPOSITES BASED ON THE ARRAYS OF CARBON NANOTUBES IN THE MICROWAVE FREQUENCY RANGE

Abstract.

Background. The nanostructured materials, based on carbon nanotubes (CNTs), have attracted a special interest due applications thereof in the fields of microwave absorption and shielding. The goal of the present work is the theoretical research of

interaction of electromagnetic waves with arrays of CNTs and the calculation of the effective permeability for the anisotropic CNT nanocomposites in the microwave frequency range using mathematical modeling at the electrodynamic level of rigor.

Materials and methods. Using the numerical method of autonomous blocks with Floquet channels (FABs) the authors carried out mathematical modeling of propagation of electromagnetic waves in the 3D- periodic arrays of CNTs at microwaves.

The computational algorithm for determining the conductivity matrix of FAB in the form of a rectangular parallelepipeds, containing CNT, was developed to solve 3D-diffraction boundary problem using the Galerkin’s projection method. The methodology of calculation of the effective permeability of nanocomposites based on the arrays of CNTs was developed based on the effective medium theory.

Results. The authors obtained the results of electrodynamic calculation of the frequency dependencies of the complex effective permeability for the CNT nanocomposites, at different orientation of the rf electric field of the wave with respect to CNT axis depending on the CNT filling factor of materials at microwave frequencies of 2-5 GHz.

Conclusions. It is shown that the dielectric losses and the absorption in arrays of CNTs increase due to the conductivity of CNTs, the CNT filling factor of materials, change depending on the orientation of the rf electric field of wave with respect to CNT axis and increase with parallel orientation in comparison with the orthogonal one.

Key words: carbon nanotubes, nanocomposite, effective permeability, electromagnetic waves, propagation, autonomous blocks

Введение

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

Особый интерес для создания высокоэффективных систем защиты от электромагнитного излучения СВЧ представляют нанокомпозиты на основе углеродных нанотрубок (УНТ) [1-3], благодаря их низкому удельному весу и соответствующим электронным свойствам.

Углеродные нанотрубки характеризуются высокими значениями удельной поверхностной проводимости и аспектного отношения размеров (длины и диаметра УНТ l/d), они сильно анизотропны, их характеристики (включая и электродинамические) отличаются в продольном и поперечном направлении УНТ [1].

В работе [3] представлены результаты экспериментальных исследований взаимодействия электромагнитного излучения с массивом вертикально ориентированных многослойных УНТ (рис. 1). Исследованы характеристики экранирования массивом вертикально ориентированных УНТ при различной ориентации нанотрубок относительно вектора электрического поля Е электромагнитной волны и УНТ, изотропно распределенными в эпоксидной матрице. Экспериментальные исследования коэффициента ослабления выполнялись при различной ориентации УНТ в волноводе, т.е. при распространении электромагнитной волны вдоль оси УНТ и в направлении, перпендикулярном оси УНТ. Обнаружено, что величина ослабления электромагнитного излучения массивом УНТ достаточно высока и изменяется в пределах от 25 до 38 дБ

(при толщине образца 1 мм в волноводе) в зависимости от взаимной ориентации вектора электрического поля волны и осей УНТ [3].

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

Нанокомпозит на основе массива периодической 3Б-решетки ориентированных УНТ рассматриваем как периодическую 3Б-наноструктуру (рис. 2,а). Элементарная ячейка периодической 3Б-наноструктуры моделируется автономным блоком, содержащим углеродную нанотрубку, и виртуальными каналами Флоке на гранях (ФАБ) (рис. 2,6).

Сложная структура наноматериала требует определения дескриптора ФАБ, содержащего УНТ. Дескриптор (в линейном приближении в виде матрицы рассеяния R или проводимости Y) ФАБ определяем из решения краевой 3Б-задачи дифракции проекционным методом Галеркина.

Построим вычислительный алгоритм решения краевой 3Б-задачи дифракции для определения матрицы проводимости Y ФАБ в виде прямоугольного параллелепипеда, содержащего УНТ (рис. 3).

Запишем систему уравнений для областей ФАБ, используя кусочнонеоднородную функцию заполнения полости ФАБ (прямоугольного параллелепипеда):

Рис. 1. Массив углеродных нанотрубок [1]

1. Математическая модель

rot H = i roeE, rot E = — ro|i0 H,

(1)

где

£ - диэлектрическая проницаемость УНТ; E, H - векторы напряженности электрического и магнитного поля; £д, М-о - электрическая и магнитная постоянные.

-Л-V 1

і I

/'

* і Г I I

и

Я

' I I I

I.

• Т "Л"

„К: ✓1 і ■ і

4-і. і

✓ I

I /

-У-

✓ I I I

I ✓ I '

✓ I I I

_ -I .

... l

УГ ■

I

-a'-

а)

rc

Рис. 2. Расчетная схема математической модели композита на основе углеродных нанотрубок: а - периодическая 3Б-структура композита; б - моделирование ячейки периодической 3Б-наноструктуры автономным блоком с каналами Флоке: 1 - углеродные нанотрубки; a, Ь, c - геометрические размеры ФАБ

x

у 4

o

Уб , „

\Г*

1 ' z

71

/

o„

I /

I /

I /

V z6 z4 _!/

v.:s "гг-----Г--7І—

У5

o

унт

/

S

S4

x

/ z5! /---

‘//V0 S S‘z

___________S________________

/1 / I / I I

x

У1

1*71

o

z1 o

Уз V °З

x

5\

Виртуальные каналы Флоке

Рис. 3. Автономный блок в виде прямоугольного параллелепипеда, содержащего УНТ, и виртуальными каналами Флоке на гранях:

V - основная область ФАБ, Vунт - область УНТ,

oaxaу0za (а = 1,2,..., 6) - локальные системы координат для сечений S0

1

Построим проекционную модель [4] для системы дифференциальных уравнений (4). Сформулируем вспомогательную краевую задачу на собственные значения (частоты) [4] для прямоугольного резонатора с геометрическими размерами основной области ФАБ (область V) на рис. 3):

1°; Йк = і щ Є) Ёк; |

^ ^ > в области Vо , (2)

1°Ёк =-іЩ Цо Йк >

Ёк(Sl) = Ёк(S4), Йк(Sl) = Йкф); '

Ёк (S2) = Ёк (S5), Йк ^2) = Йк (S5); > на гранях,

Ёк(Sз) = Ёк(S6), Йк^) = Йк(S6) ,

где Щк - собственные частоты резонатора; Ёк, Йк - собственные электрические и магнитные поля (собственные функции) резонатора.

Система собственных функций {Ёк, Йк} состоит из соленоидальной

{Ёк,Йр} и потенциальной подсистем {Ё^,Й^-} [4]. Индекс к определен на

множестве индексов к' и к" .

Собственные функции ортогональны и нормированы:

Цо \ Й\- ЙndV = £) \ Ёк • Ё^ = дк„. (3)

V) V)

Запишем выражения для собственных электромагнитных волн виртуальных каналов Флоке ФАБ [5]:

Ч(а) =(а) ±\а)) exp( Tk(а) z

а

Ч(а) =(±((а) ±%)) exp(±iTk(а) -а); (4)

H,

(а)

k = 1,2,..., ^; а = 1,2, ...,б;

где к - номера мод собственных волн каналов Флоке; (а), (а) - поперечные электрические и магнитные компоненты собственных волн; ,

Нк(а) - продольные электрические и магнитные компоненты собственных волн; Г к (а) - постоянные распространения собственных волн.

Поперечные электрические и магнитные компоненты собственных волн каналов Флоке образуют полную систему ортогональных функций

{вк (а), Нк (а)} [2]. Любое поперечное электромагнитное поле на входных сечениях £0 автономного блока представляем в виде разложения по этим системам {вк(0),Нк(0)} в ортогональные ряды Фурье:

Ёа 2 ak (a)ek (а)

k=1

Йа = 2 bk (a)hk(a), k=1

а = 1, 2,... ,6,

где с^, ек(а) - амплитуды падающих и отраженных волн. Из рядов Фурье (5) и их нормировки [2]

[0, к ■£ п,

|1, к = п,

J (k(а) x ^г(а) ) • dSa

(5)

(6)

следуют следующие интегральные выражения:

ak(а) = J (Ёа xhk(a)) • dSa ; Sa

bk(а) = J (ek(а) xЙа) •dSа .

(7)

(S)

Выражения (7), (S) являются интегральными краевыми условиями на гранях автономного блока и известны в электродинамике [4] как условия неасимптотического излучения.

Используя краевую задачу на собственные значения (2), тождество векторного анализа brot a - a rot b = rot (a x b), формулу Остроградского - Гаусса и условие неасимптотического излучения (S), запишем для системы уравнений (1) проекционную интегральную модель:

^J (Ё x Щ) • dSp=- i ®k eo J P=1 Sr

Vo

Ё• ЁkdV-iro^o JЙ • ЙkdV,

vo

О

2 J (Й x Kk) • dSa = i® J еЁ • KkdV + i ®k ^o J Й • Йk dV,

(9)

P=1 Sr

Vo

Vo

bq(a) J (eq(а) xЙа ) • dSa,

к = 1,2,...N0, Ч = 1,2,..., а = 1,2,...,6, где N - количество базисных функций в области автономного блока; Ыа - количество базисных функций на

гранях автономного блока.

Решение краевой задачи ищем в виде линейной комбинации по системам функций {Еп }, {Нп } (собственные функции прямоугольного резонато-

ра), {%)}, {'

/(Р)

(собственные функции каналов Флоке).

а

а

а

В области V ФАБ (рис. 3):

No No

ё = 2 anЁn, й = 2 ЬпЙп . n=1 n=1

(1o)

На гранях Sp(Р = 1,2,...,6) ФАБ:

Nr

Nr

ёр = 2a/ (в)е/ (в), йр=2b/ (в)h/ (в).

/=1 /=г

(11)

Подставляя (10) и (11) в (9), получаем следующую систему алгебраических уравнений:

А • а + В • Ь = -Ь • а,

Б • а + и • Ь = 0,

W • Ь = Ь,

где А, В, Б, и, W, Ь - матрицы с элементами:

Акп = * ®к ^ кп , Вкп = * ®^кп :

(12)

Dkn =

i ® J е En • Ёр dV i ® J е Ёп • Ё|* dV Л

Vo Vo

i ® J е En • Ёп* dV i ® J е ё1П" • Ef dV

7И рп*

V Vo

Vo

Ukn =

'®k'h'n' - 2 J (Hcx ЁС*) • dSp i ®k*kV - 2 J (Йп x ЁС*)

• dSR

P=1 SR

P=1 SR

i ®k'SkV - 2 J (ЙС x Ёп ) • dSp i ®k hn - 2 J (Йп x Ёп*)

• dSR

P=1 SR

P=1 Sr

Wq

q(a)n

J (eq(a) xЙ%) • dSa J (eq(a) xЙп) • dS0

V Sa

L

k/ (P)

J (e/(P) x ЙС*) • dSp

SP

J (e/(P) x Hk*) • dSp

v SP

, k = 1,2,...No, q = 1,2,...Na, a = 1,2,...,6.

Компонентами векторов а, Ь, а, Ь являются коэффициенты рядов Фурье (10) и (11), которые равны |йп } , [Ь } , |й(р)}, |Ь(р)} .

Исключая векторы а, Ь из системы линейных алгебраических уравнений (12), получаем

\-1

b = W-(а - D-1 - U - б) - L -

откуда следует матрица проводимости ФАБ:

Y = W

•(А - D-1

-1

Б I - L .

(1З)

(14)

2. Методика и результаты электродинамического расчета эффективной диэлектрической проницаемости нанокомпозита на основе углеродных нанотрубок

Геометрия задачи - направление распространения волнового процесса и электродинамическая модель нанокомпозита на основе массива - периодической 3Б-решетки ориентированных УНТ, который рассматриваем как ква-зипериодическую 3Б-наноструктуру с геометрическими размерами ячейки a, b, с (рис. 3).

Свободный электромагнитный процесс, в 3Б-периодической наноструктуре (рис. 4) подчинен теореме Флоке [4] в форме

E(x + a, y, z) = E(x, y, z) ехр(-/фх), H(x + a, y, z) = H(x,y, z) exp(-iqx),

E(x,y + b, z) = E(x, y, z)exp(-i'9y), H(x,y + b, z) = H(x, y, z) exp(-i'9y),

E(x, y, z + c) = E(x,y, z)exp(-i'9z), H(x, y, z + c) = H(x, y, z) exp(-i'9z), (15)

где фх =r„acospx, = rnbcospy, = r„ccospz ; (3X,Py,Pz - углы ори-

ентации направления распространения волнового процесса.

а)

z

б)

Рис. 4. Электродинамическая модель нанокомпозита на основе периодической 3Б-решетки ориентированных УНТ: а - направление распространения волнового процесса с волновым вектором к; б - периодическая 3Б-наноструктура

Ячейку периодической ЗБ-наноструктуры моделируем ФАБ с дескриптором в виде матрицы проводимости У. Для матрицы проводимости ФАБ (14) справедливо выражение

b = Y • a.

(1б)

где векторы а и Ь составлены соответственно из коэффициентов ап(р) и

Ь^и) (а,в = 1,2,...6; k,п = 1,2,...Naв) рядов Фурье (5). Тогда условия теоре-

мы Флоке (15) принимают следующий вид:

ап(4) = ап(1) ехр(-/фx), aп(5) = ап(2) ехр(-^у ),

aп(6) = aп(3) ехР(-% X Ь(4) = Ч(1) ехР(-ФхX (17)

Ч(5) = Ч^М-ЩуX Ч(6) = bk(3^)QxР(-iФzX 1с,п = 1,2,...^,р.

Подставляя (17) в (16), получаем характеристическое уравнение для определения постоянных распространения Гп волн в периодической 3Б-наноструктуре на основе периодической 3Б-решетки ориентированных УНТ:

А(Гп ) =

Глл - H-1 • Гвл + Глв • H - H-1 • Гвв • H

= о.

(18)

где А(Гп ) - определитель матрицы; Глл, Гвл, Глв, Гвв - клетки матрицы

проводимости Г =

Г Г,

лл Глв Гвл Гвв

(л, в - индексы входных сечений ФАБ

S1, S2, Sз и S4, S3, S6 соответственно); H =

Г hx о о

о hy о

о о h

- диагональная мат-

z

рица с элементами: hx(l j) = - 8l j Гпа cos Рx , hy(l j) = - 8l j rn b cos Py , hz(l j) = - 8l j rn c cos Pz •

Предметом исследования в периодической 3Б-наноструктуре является нулевая пространственная гармоника, которой соответствует постоянная распространения волны Г о •

Методика определения комплексной эффективной диэлектрической проницаемости нанокомпозитов на основе массивов УНТ базируется на теории эффективной среды: постоянные распространения волн в неограниченной сплошной среде совпадают с аналогичными постоянными распространения волн в периодической 3Б-наноструктуре. В этом случае комплексная эффективная диэлектрическая проницаемость определяется следующим образом:

еХ

Го2

(19)

На рис. 5, 6 показаны результаты расчета частотных зависимостей комплексной диэлектрической проницаемости нанокомпозита на основе периодической 3Б-решетки ориентированных УНТ для ориентации вектора напряженности электрического поля Е волны параллельно (рис. 5) и ортогонально (рис. 6) оси УНТ при различном соотношении V объема УНТ к общему объему композита (выбор периода решетки а,Ь - выполнение этого соотношения).

Рис. 5. Частотные зависимости действительной и мнимой частей комплексной диэлектрической проницаемости нанокомпозита на основе периодической 3Б-решетки ориентированных УНТ при различном соотношении V объема УНТ к общему объему композита для ориентации вектора электрического поля Е волны параллельно оси УНТ: А = 3 нм, 2г = 25 нм, I = 1000 нм, с = 1025 нм;

Ет , Нт - ориентация электрического и магнитного полей квази-ТЕМ волны;

кривая 1 - Еунт = 62, Оунт = 2,5 Ом-1м-1; кривая 2 - Еунт = 62,

о унт = 5,0 Ом 1м 1; кривая 3 - Еунт = 62, о унт = 7,5 Ом 1м 1;

кривая 4 - Еунт = 62, Оунт = 10,0 Ом м

2.5

3.5

4.5

f, ГГц

Рис. 5. Окончание

2

3

4

Как следует из результатов математического моделирования (рис. 5, 6), величина ослабления излучения в СВЧ-диапазоне достаточно высока, что связано с хорошей электропроводностью УНТ.

Рис. 6. Частотные зависимости действительной и мнимой частей комплексной диэлектрической проницаемости нанокомпозита на основе периодической 3Б-решетки ориентированных УНТ при различном соотношении V объема УНТ к общему объему композита для ориентации вектора электрического поля Е волны ортогонально оси УНТ: А = 3 нм, 2г = 25 нм,

I = 1000 нм, с = 1025 нм; кривая 1 - е унт = 62, о унт = 2,5 Ом-1м-1; кривая 2 - еунт = 62, оУНТ = 5,0 Ом-1м-1; кривая 3 - Еунт = 62, о унт = 7,5 Ом 1м 1; кривая 4 - Еунт = 62, о унт = 10,0 Ом 1м 1

2.5

3.5

4.5

f, ГГц

Рис. 6. Окончание

2

3

4

Как видно из графиков (рис. 5, 6), чем выше проводимость УНТ и соотношение V объема УНТ к общему объему композита, тем больше мнимая часть комплексной диэлектрической проницаемости нанокомпозита , и, следовательно, выше степень ослабления электромагнитного излучения за счет поглощения. Для массива УНТ ослабление зависит от взаимной ориентации вектора напряженности электрического поля Em волны и осей углеродных нанотрубок: при параллельной (рис. 5) ориентации в сравнении с ортогональной (рис. 6) ослабление СВЧ излучения увеличивается.

Список литературы

1. Дьячков, П. Н. Углеродные нанотрубки: строение, свойства, применения / П. Н. Дьячков. - М. : БИНОМ. Лаборатория знаний, 2006. - 294 с.

2. Decrossas, E. Carbon Nanotubes for Electromagnetic Compatibility Applications / E. Decrossas, M. A. EL Sabbagh, V. F. Hanna and S. M. El-Ghazaly // Electromagnetic Compatibility (EMC), IEEE International Symposium (6-10 Aug. 2012). - Pittsburgh, PA, 2012. - P. 428-433.

3. Во вченко, Л. Л. Резонансный характер взаимодействия многослойных углеродных нанотрубок с излучением миллиметрового диапазона волн / Л. Л. Вовченко, Л. Ю. Мацуй, В. В. Олейник, В. Л. Лаунец, В. В. Загородний, Ф. Ле Норманд // Наносистеми, наноматеріали, нанотехнологп. - 2011. - Т. 9, № 4. - С. 759-769.

4. Никольский, В. В. Проекционные методы в электродинамике / В. В. Никольский // Сборник научно-методических статей по прикладной электродинамике. -М. : Высшая школа, 1977. - С. 4-23.

5. Голованов, О. А. Автономные блоки с виртуальными каналами Флоке и их применение для решения прикладных задач электродинамики / О. А. Голованов // Радиотехника и электроника. - 2006 - Т. 51, № 12. - С. 1423-1430.

References

1. D'yachkov P. N. Uglerodnye nanotrubki: stroenie, svoystva, primeneniya [Carbonic nanotubes: structure, properties, application]. Moscow: BINOM. Laboratoriya znaniy, 2006, 294 p.

2. Decrossas E., EL Sabbagh M. A., Hanna V. F. and El-Ghazaly S. M. Electromagnetic Compatibility (EMC), IEEE International Symposium (6-10 Aug. 2012). Pittsburgh, PA, 2012, pp. 428-433.

3. Vovchenko L. L., Matsuy L. Yu., Oleynik V. V., Launets V. L., Zagorodniy V. V., Le Normand F. Nanosistemi, nanomatenali, nanotekhnologu [Nanosystems, nanomaterials, nanotechnologies]. 2011, vol. 9, no. 4, pp. 759-769.

4. Nikol'skiy V. V. Sbornik nauchno-metodicheskikh statey po prikladnoy elektrodinamike [Collection of scientifical and methodological articles on applied electrodynamics]. Moscow: Vysshaya shkola, 1977, pp. 4-23.

5. Golovanov O. A. Radiotekhnika i elektronika [Radio engineering and electronics]. 2006, vol. 51, no. 12, pp. 1423-1430.

Голованов Олег Александрович доктор физико-математических наук, профессор, кафедра общепрофессиональных дисциплин, Филиал Военного учебно-научного центра сухопутных войск «Общевойсковая академия Вооруженных сил РФ» (Россия, г. Пенза-5)

E-mail: golovanovol@mail.ru

Макеева Галина Степановна

доктор физико-математических наук, профессор, кафедра радиотехники и радиоэлектронных систем, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

E-mail: radiotech@pnzgu.ru

Ширшиков Дмитрий Николаевич

адъюнкт, филиал военной академии материально-технического обеспечения (Россия, г. Пенза-5)

E-mail: shirshikov1981@mail.ru

Горлов Геннадий Геннадьевич

адъюнкт, филиал военной академии материально-технического обеспечения (Россия, г. Пенза-5)

E-mail: mitsubisi-gor82@mail.ru

Golovanov Oleg Aleksandrovich Doctor of physical and mathematical sciences, professor, sub-department of interprofessional disciplines, branch of the Military Educational and Research Center of the land forces “Combined Arms Academy of the Armed Forces of RF” (Penza-5, Russia)

Makeeva Galina Stepanovna Doctor of physical and mathematical sciences, professor, sub-department of radio engineering and radio electronic systems, Penza State University (40 Krasnaya street, Penza, Russia)

Shirshikov Dmitriy Nikolaevich Postgraduate student, branch of the Military academy of maintenance supply (Penza-5, Russia)

Gorlov Gennadiy Gennad'evich Postgraduate student, branch of the Military academy of maintenance supply (Penza-5, Russia)

УДК 538.95 Голованов, О. А.

Электродинамический расчет комплексной эффективной диэлектрической проницаемости нанокомпозитов на основе массивов углеродных нанотрубок в диапазоне сверхвысоких частот / О. А. Голованов, Г. С. Макеева, Д. Н. Ширшиков, Г. Г. Горлов // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. -№ 1 (29). - С. 141-155.