Научная статья на тему 'Оптимизация ВТСП кабелей постоянного и переменного тока с учетом эффекта продольного магнитного поля'

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

CC BY
20
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВТСП-кабели / эффект продольного магнитного поля / численный анализ / HTS cables / longitudinal magnetic field effect / numerical analysis

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Зубко В. В., Занегин С. Ю., Фетисов С. С., Высоцкий В. С., Носов А. А.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Зубко В. В., Занегин С. Ю., Фетисов С. С., Высоцкий В. С., Носов А. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Optimization of HTS DC and AC cables with take account longitudinal magnetic field effect

It is known that in magnetic fields parallel to the transport current, an increase of the critical current is observed in most low-temperature and high-temperature superconductors — the effect of a longitudinal magnetic field (LMFE). LMFE has been theoretically predicted and experimentally confirmed in many studies. This article presents various methods for optimizing power coaxial cables based on high-temperature superconductors (HTS), considering LMFE. LMFE provides a higher current-carrying capacity of HTS DC cables and higher stability for HTS power AC cables, including lower AC losses. The results of optimization of HTS cables for both DC and AC applications with consideration of LMFE are presented and discussed. Two methods for calculating the magnetic field in coaxial HTS cables are presented: analytical and using the finite element method. The results of AC loss calculations and measurements for HTS cables are also presented.

Текст научной работы на тему «Оптимизация ВТСП кабелей постоянного и переменного тока с учетом эффекта продольного магнитного поля»

УДК 538.945 Статья

Оптимизация ВТСП кабелей постоянного и переменного тока с учетом эффекта продольного магнитного поля

В.В. Зубко1'2'*, С.Ю. Занегин1,2, С.С. Фетисов1, В.С. Высоцкий1, А.А. Носов1, Э.С. Отабе3, Т. Акасака 4

1 ОАО «Всероссийский научно-исследовательский проектно-конструкторский и технологический

институт кабельной промышленности» («ВНИИКП»), 111024, Москва, Россия

2 Московский авиационный институт (национальный исследовательский университет), 125993,

Москва, Россия

3 Kyushu Institute of Technology, 680 - 4, Iizuka, Japan

4 Superconductivity and Cryogenics, Maglev Systems Technology Division, Railway Technical Research

Institute, 2-8-38 Hikaricho

* e-mail: [email protected]

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

Ключевые слова: ВТСП-кабели; эффект продольного магнитного поля; численный анализ.

DOI: 10.62539/2949-5644-2024-0-1-53-62

1. Введение

При разработке силовых ВТСП-кабелей основной целью является достижение его максимально высокой токонесущей способности при минимальном использовании основных сверхпроводников. Разработано множество методов оптимизации для увеличения критических токов в ВТСП-кабелях. Одной из возможных идей является использование так называемого эффекта продольного магнитного поля [1].

С конца 1960-х годов известно, что в низкотемпературных сверхпроводниках, наблюдается заметное увеличение плотности критического тока в продольном магнитном поле (поле параллельно транспортному току) в несколько раз больше, чем её величина в поперечном магнитном поле (поле перпендикулярно транспортному току) [2, 3]. Такое увеличение критического тока было объяснено бессиловой конфигурацией для магнитного поля и тока в сверхпроводнике в продольном магнитном поле [4, 5]. Это условие предполагает, что магнитное поле и ток практически параллельны друг другу.

ЭПМП для ВТСП-кабелей был теоретически предсказан позднее в работе [6]. Конструкция силового кабеля постоянного тока на основе высокотемпературных сверхпроводников второго поколения (ВТСП-2), использующая ЭПМП была предложена в работе [7].

Недавно были разработаны два коаксиальных кабеля постоянного тока с токонесущими элементами из ВТСП-2 лент с учетом ЭПМП [1, 8]. Эти кабели прошли успешные испытания и продемонстрировали увеличение критического тока и, следовательно, увеличение токонесущей способности за счет ЭПМП. Чтобы добиться увеличения критического тока в этих кабелях, их

конструкция (шаг скрутки в зависимости от диаметра повива) должна обеспечивать параллельное магнитное поле для каждой ВТСП ленты в каждом повиве. В этом случае критический ток лент будет увеличиваться за счет ЭПМП с увеличением суммарного критического тока кабеля.

В данной работе мы продолжаем исследования, начатые нами ранее в работе [1], распространяя наш анализ на возможные кабели переменного тока с использованием ЭПМП. Для нашего анализа мы также используем данные из работы [8].

Очевидно, что способы оптимизации ВТСП-кабелей постоянного и переменного тока с использованием ЭПМП должны быть разными, и это обсуждается в данной статье. Поскольку потери на переменном токе очень важны, мы также оценили потери на переменном токе в оптимизированном кабеле с использованием ЭПМП из работы [1].

2. Зависимости критического тока ВТСП-2 лент

Экспериментальные зависимости критических токов современных ВТСП-2 лент от магнитного поля для различных углов между транспортным током и внешним магнитным полем 1с(В,р) приведены в работах [1,8]. Зависимости критического тока для ВТСП-2 лент из [1,8] для параллельного 1с(В,0°) и перпендикулярного 1с(В,90°) магнитных полей по отношению к току при температуре 77.4 К показаны на рис. 1.

Рис. 1. (а) Экспериментальная зависимость критического тока от магнитного поля для ВТСП-2 лент из [1] и [8] для углов р =0° (продольное магнитное поле (By)) и @ = 90° (поперечное магнитное поле).

Зависимости Ic(B,fi) представленные в работах [1,8] могут быть аппроксимированы следующим уравнением:

1С(В,р) = 0.5(7С (В,0) + 1С(В,90)) + 0.5(Ic(В,0) - Ic(В,90)) cos (2р), (1)

Ici, приведенную на рис. 1, можно рассматривать как зависимость критического тока от магнитного поля при отсутствии ЭПМП (Ic(B,90°)), а Ic2 — при наличии максимального ЭПМП (Ic(B,0°)). Из рис. 1 видно Ic (собственное поле) = 150 А для лент, представленных в [1], и Ic

(собственное поле) = 162 А, для лент, представленных в [8]. На рис. 1 также показаны направления магнитного поля и тока.

3. Оптимизация коаксиальных ВТСП-кабелей постоянного тока с учетом эффекта продольного магнитного поля

Токонесущая способность ВТСП кабеля постоянного тока равна сумме критического тока каждого повива проводника (жилы) кабеля:

=^л=\1сл> (2)

где т количество повивов в проводнике.

Для кабеля постоянного тока задачей оптимизации является разработка конструкции кабеля с максимальной токонесущей способностью, т.е. целевая функция — это максимум 1С ш. Для оптимизации кабеля постоянного тока (с точки зрения получения максимума 1С ш.с учетом ЭПМП) был использован итерационный алгоритм:

Р.

1 с,г'

(и (3)

сначала для каждого () повива ВТСП-кабеля задавалось начальное значение тока 1С!° = 1С(Вг- = 0), далее вычислялось магнитное поле (В{) и угол (Зг) в каждом повиве кабеля, затем вычислялось !С г = Затем на каждом новом шаге ток в каждом повиве изменяется до тех пор, пока

ток в каждом повиве будет удовлетворять условию:

тк+1_ тк

c,i c,i

C,l

где £ = 1 0 4 - остановочный критерий.

Для вычисления магнитного поля В¿, аксиальной В2 г- и азимутальной Вд ^ составляющих магнитного поля в повиве коаксиального ВТСП кабеля можно использовать аналитические уравнение [7]:

m

в Vohsin (вк) | в

I

И ohcos(e к)

]=ti2 п Rkcos(ek) ext'

m

2 nR, '

k=i+1 1

M ВЬ + Bef2* (5)

где Ik — ток в повиве, Rk — внутренний радиус повива, вк — угол скрутки (угол наложения) ВТСП лент в повиве, Вext — внешнее магнитное поле, создаваемое повивами экрана кабеля.

Кроме того, мы можем использовать 3D-модель на основе метода конечных элементов (МКЭ) для определения В,-, Bz^Bq ¡. Данная модель с использованием ANSYS (EMAG software) [9] подробно представлена в [10,11]. В этой модели каждый, повив проводника и экрана моделируются тонким цилиндром, имеющим толщину ВТСП слоя ВТСП-2 ленты. Угол скрутки

ВТСП лент в повиве моделируется анизотропией электропроводности каждого цилиндра. Электропроводность анизотропного тонкого проводника представляет собой тензор 2-го порядка, где все недиагональные элементы равны нулю, данный тензор имеет вид:

о =

О 0 0

0 о 0

0 0 '

(6)

где ох' = о ^ ^ охг — электропроводность в соответствующих координатах, г/ - координата,

параллельная направлению намотки ВТСП ленты в повиве.

После определения В¿, Ви В ^ ¿, магнитное поле параллельное току в лентах (В \\ Д и угол для каждого повива определяется уравнениями:

Ву = В ^ ■ сов (в {) + В 0 б т (в{),

(7)

(8)

Далее был проведен новый анализ ВТСП кабелей постоянного тока разработанных с учетом ЭПМП: компактного кабеля с четырьмя ВТСП-повивами в жиле и двумя повивами в экране [1] и кабеля класса 10 кА с четырьмя повивами в жиле и четырьмя повивами в экране [8].

На рис. 2 показаны расчетные значения критического тока кабелей с учетом и без учета ЭПМП в зависимости от угла скрутки ВТСП лент в крайнем повиве ( втах) проводника,

о — о

остальные углы определялись уравнением в1 = 5°, в^ = —■ Угол скрутки ВТСП-лент в

повивах экрана одинаковы и равны 30о. Расчеты проводились для кабелей, когда ВТСП ленты в жиле скручены в одну сторону, а ленты экрана — в противоположную; например, направление намотки кабеля [1] а = [1,1 ,1,1, - 1, - 1]. Внутренние радиусы (Я) повивов кабелей взяты из [1,8]. Из рис. 2 видно, что увеличение критического тока кабеля (около 1 кА) можно получить, учитывая ЭПМП. Расчет ограничен вд^ах <;35° , поскольку общая длина ВТСП ленты увеличивается по мере увеличения, и использование большого количества сверхпроводника нецелесообразно. Согласно [7], максимальный критический ток кабеля получается при значениях угла ниже 50о.

< г*

а

с р

>к к ы

и

О! г

г —

г о.

Кабель [8] е- учет ЭПМП нет

Кабель [1] — учет ЭПМП 1— нет

10 15 20 25 Угол скрутки, втах, град

30

Рис. 2. Токонесущая способность ВТСП к<беая (1С 1о1) в зависимости от

Для примера, на рис. 3 показано распределение магнитного поля для вшах = 3 4° при критическом токе в кабеле [8]. Видно, что максимальное магнитное поле получается в сверхпроводящем повиве и составляет около 0.18 Тл. Для кабеля из [1] рассчитанные значения магнитных полей примерно такие же.

Таким образом, нам удалось разработать метод оптимизации для кабеля постоянного тока с ЭПМП.

Рис. 3. Распределение плотности магнитного потока при критическом токе в кабеле с 4 повивами в проводнике и 4 повивами в экране (0max = 34°).

4. Оптимизация коаксиальных ВТСП-кабелей переменного тока с учетом эффекта продольного магнитного поля

Для ВТСП кабелей переменного тока задача оптимизации состоит в обеспечении равномерного распределения токов между повивами проводника и между повивами экрана кабеля и может быть описана как задача минимизации целевой функции:

т-1 т N-\ N

F(X) = min/(X) = £ £ \1,(Х)-1к(Х)\+ £ £ (9)

А=1 i=i+l t=m+l i=k+l

где Ij(X) и 1к{Х) - значение тока в -м и -м повиве, т- количество повивов в проводнике, N -общее количество повивов (включая проводник и экран) в кабеле. Вектор управляющих переменных:

X [R1,ei,a1,^,Ri,ei,ai,^,RN,eN,aN], (10)

где К,- - внутренний радиус -го повива, яг- and в — направление и угол скрутки ВТСП-лент в -м повиве.

Для расчета токов в каждом - м повиве мы можем использовать уравнения Кирхгофа

Кроме того, для полученного токов в каждом повиве мы можем использовать 3D-МКЭ модель, как описано в разделе 3, см. также [10, 11], где управляющими переменными являются направление электропроводности каждого цилиндра (повива) (см. уравнение (6)).

Ограничением на управляющие переменные является условие сохранения сверхпроводящих свойств ВТСП лент в повиве.

Для надежности ВТСП кабелей переменного тока необходим запас по току, для этого отношение I,- / 1сг- выбирается около 0.7-0.8. На практике при заданных максимальных токах в повивах проводника запас по току рассчитывается без учета ЭПМП и, следовательно, это значение оказывается завышенным. Для более точного расчета запаса по току следует учитывать ЭПМП.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Если заданы максимальные токи в повивах проводника, то можно увеличить запас по току (или увеличить 1с г) повивов проводника с помощью ЭПМП. Для этого на значения переменных 0г- для каждого повива экрана необходимо добавить ограничения, которое должно быть > 30°, что обеспечит максимальное продольное магнитное поле, создаваемое повивами экрана в повивах проводника.

5. Расчет и измерение потерь на переменном токе в коаксиальном ВТСП кабеле

Проверка влияния ЭПМП на потери на переменном токе проведена на ВТСП кабеле из работы [1], который оптимизирован с учетом ЭПМП.

Измерения потерь на переменном токе проводились электрическим методом, который основан на расчете активной мощности (потери переменного тока) для каждого периода тока для проводника по формуле:

(11)

где и, г - мгновенные значения напряжения и тока, £т — время измерения нескольких периодов сигнала для повышения точности.

Измеренные потери в жиле ВТСП кабеля и в каждом повиве в отдельности показаны на рис. 4 [1]. На рис. 4 повивы отсчитываются от внутреннего повива, в каждом повиве имеется по 14 лент [1]. Измеренные потери в проводнике превышают измеренные потери в каждом повиве в отдельности, поскольку в проводнике наблюдается незначительная неравномерность распределения тока между повивами проводника, поскольку проводник не был оптимизирован для использования на переменном токе [1], также повивы проводника имеют разные критические токи [1].

Рис. 4. Измеренные потери на переменном тока не ленту в зависимости от относительного тока в проводнике и в каждом повиве отдельно.

На рис. 5 показано сравнение измеренных и рассчитанных МКЭ моделью потерь в проводнике. В расчете учитывалось распределение измеренного тока между повивами и измеренный критический ток каждого повива [1].

Для расчета потерь в ВТСП-2 кабелях на переменном токе была разработана 2D МКЭ модель (ANSYS EMAG software, A/A-V формулировка) [9, 11, 12]. В модели рассматривается полное сечение кабеля, учитывая его полигональную структуру и зазоры между лентами, рассчитываются только радиальная и азимутальная составляющие магнитного поля в кабеле. Поэтому мы попытались смоделировать влияние ЭПМП на потери следующим образом, расчет проводился для критических токов ВТСП-2 ленты, представленных на рис. 1. Jcl (без ЭПМП) и Ic2 (максимальное ЭПМП).

Расчетные потери на ленту ниже измеренных потерь, если относительный ток (I/Ic) в проводнике меньше 0.4. По мере увеличения I/Ic из меренные по тери лу чше со впадают с расчетными если в расчетах учитывать ЭПМП. Как видно из рис. 5, вследствие ЭПМП потери в кабеле переменного тока незначительно снижаются. Причиной малого влияния на потери ЭПМП в данном кабеле является неравномерность распределения тока между повивами проводника и небольшое изменение распределения поля в кабеле с ЭПМП по сравнению с кабелем без ЭПМП. Для проверки этого эффекта предполагается изготовить более подходящий (оптимизированный с учетом ЭПМП) кабель переменного тока.

Следует отметить, что существуют две модели, объясняющие структуру внутри сверхпроводника в продольном магнитном поле: первая модель представлена в работе [13], а вторая модель представлена в работе [14]. И вокруг этого вопроса ведется очень много дискуссий. Расчеты МКЭ моделью не учитывают структуру сверхпроводника. В расчетах мы использовали экспериментальные результаты, приведенные на рис. 1 с помощью которых оценили величину потерь на переменном токе.

10

В

а

10'

Я 10 1

а ю-а>

ю"

- • 11роводнвк(пзм,

—Г I nun

• ¿г для Irl для Ii 2

fV i i

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 I

I/Ic

Рис. 5. Измеренные потери на переменном токе на ленту в зависимости от относительного тока в жиле кабеля по сравнению с расчетными потерями в проводнике.

5. Выводы

Разработаны различные методы оптимизации с учетом ЭПМП коаксиальных ВТСП кабелей постоянного и переменного тока. Учет ЭПМП позволяет: для кабелей постоянного тока увеличить критический ток (токонесущую способность), для кабелей переменного тока при заданных максимальных токах в повивах проводника повысить его надежность за счет увеличения критического тока повивов проводника. Анализ экспериментальных данных с помощью МКЭ модели показал, что, вследствие ЭПМП потери в кабеле на переменном токе незначительно снижаются.

Сделан вывод, что 10% увеличение общего критического тока ВТСП кабеля может быть достигнуто, если при оптимизации учитывать ЭПМП в современных ВТСП-2 лентах.

Литература

[1] S. S. Fetisov, et al, IEEE Trans. Appl. Supercond. 32,4 (2022). DOI: 10.1109/TASC.2021.3137763

[2] Yu. F. Bychkov, et al, JETP Lett. 9, 12, 652 (1969).

[3] B. ten Haken, et al, IEEE Trans. Magn. 28, 755 (1992). DOI: 10.1109/20.119989

[4] C.J. Bergeron, Appl. Phys. Lett. 3, 171 (1963). DOI: 10.1063/1.1753920

[5] T. Matsushita Jpn. J. Appl. Phys. 51, 010111 (2012). DOI 10.1143/JJAP.51.010111

[6] T. Matsushita, et al, Supercond. Sci. Technol., 25, 125009 (2012). DOI 2048/25/12/125009

[7] V.S. Vyatkin, et al, IEEE Trans. Appl. Supercond., 25, 4, 6606207 10.1109/TASC.2015.2427357

[8] Y. Kinoshita, et al, J. Phys., Conf. Ser. 1975, 012037 (2021). DOI: 6596/1975/1/012037

[9] Ansys Multiphysics, www.ansys.com

[10] S.S. Fetisov, et al, IEEE Trans. Appl. Supercond. 28, 4, 5400905 10.1109/TASC.2018.2821708

[11] V.V. Zubko, et al, J. Phys., Conf. Ser. 2043, 012004 (2021). DOI 6596/2043/1/012004

[12] V.V. Zubko, et al, IEEE Trans. Appl. Supercond. 26, 3, 8202005 10.1109/TASC.2016.2535151

[13] J.R. Cave, et al, J. Phys. Colloques, 39 (1978). DOI: 10.1051/jphyscol:19786276

[14] T. Matsushita, J. Phys. Soc. Jpn, 54, 3, 1066 (1985). DOI: 10.1143/JPSJ.54.1066

10.1088/0953-

(2015). DOI: 10.1088/1742-

(2018). DOI: 10.1088/1742-

(2016). DOI:

Optimization of HTS DC and AC cables with take account longitudinal magnetic field effect

V.V. Zubko1'2'*, S.Yu. Zanegin1,2, S.S. Fetisov1, V.S. Vysotsky1, A.A. Nosov1,2, E.S. Otabe3, T. Akasaka3

1 All-Russian scientific research and development cable institute («VNIIKP»), 111024, Moscow, Russia

2 Moscow Aviation Institute (MAI), 125993, Moscow, Russia

3 Kyushu Institute of Technology, 680 - 4, Iizuka, Japan * e-mail: [email protected]

It is known that in magnetic fields parallel to the transport current, an increase of the critical current is observed in most low-temperature and high-temperature superconductors — the effect of a longitudinal magnetic field (LMFE). LMFE has been theoretically predicted and experimentally confirmed in many studies. This article presents various methods for optimizing power coaxial cables based on high-temperature superconductors (HTS), considering LMFE. LMFE provides a higher current-carrying capacity of HTS DC cables and higher stability for HTS power AC cables, including lower AC losses. The results of optimization of HTS cables for both DC and AC applications with consideration of LMFE are presented and discussed. Two methods for calculating the magnetic field in coaxial HTS cables are presented: analytical and using the finite element method. The results of AC loss calculations and measurements for HTS cables are also presented.

Keywords: HTS cables, longitudinal magnetic field effect, numerical analysis.

Зубко Василий Васильевич — доктор техн. наук, главный научный сотрудник ОАО «ВНИИКП», Москва, Россия, ведущий научный сотрудник МАИ (НИУ), Москва, Россия. Dr. Vasily Zubko — Chief Researcher «VNIIKP», MAI.

Занегин Сергей Юрьевич — кандидат техн. наук, научный сотрудник ОАО «ВНИИКП», Москва, Россия, научный сотрудник МАИ (НИУ), Москва, Россия.

Dr. Sergey Zanegin — Researcher, «VNIIKP», MAI.

Фетисов Сергей Сергеевич — доктор техн. наук, заведующий отделением,

ОАО «ВНИИКП», Москва, Россия.

Dr. Sergey Fetisov — Head of Department, «VNIIKP».

Высоцкий Виталий Сергеевич — доктор техн. наук, заместитель заведующего, научный руководитель отделения сверхпроводящих и криорезистивных кабелей, ОАО. «ВНИИКП», Москва, Россия. Dr. Vitaly Vysotsky — Ph.D; Dr. of Sci. (habil.) Deputy head of Department, «VNIIKP».

Носов Александр Анатольевич — кандидат техн. наук, старший научный

сотрудник ОАО «ВНИИКП», Москва, Россия.

Dr. Aleksandr Nosov — Senior Researcher «VNIIKP»

Эдмунд Соджи Отабе, профессор, Кюсю Технологический Институт, Иизука. Япония.

Dr. Edmund Soji Otabe — Professor, Kyushu Institute of Technology, Iizuka, Japan.

Томоюки Акасака — научный сотрудник, Кюсю Технологический Институт, Иизука. Япония.

Tomoyuki Akasaka — Researcher, Kyushu Institute of Technology, Iizuka, Japan.

i Надоели баннеры? Вы всегда можете отключить рекламу.