УДК 538.955.35
МОДЕЛИРОВАНИЕ СИСТЕМ МАГНИТНОЙ ЛЕВИТАЦИИ НА ОСНОВЕ СВЕРХПРОВОДЯЩИХ КОЛЕЦ
И. В. Анищенко, С. В. Покровский, И. А. Руднев
Разработана комплексная модель магнита на основе ВТСП колец. Для расчетов использованы магнитные и транспортные характеристики промышленных ВТСП лент, учтены особенности их слоистой структуры. Проведен расчет параметров магнитной системы во внешних магнитных полях различных конфигураций. Выполнен расчет магнитных полей, токов, намагниченности, силы левитации для различных геометрических параметров намагничиваемых колец.
Ключевые слова: сверхпроводящие кольца, ВТСП ленты, сила левитации.
Введение. Высокотемпературные сверхпроводники (ВТСП) нашли применение в качестве элементов левитационных подвесов, магнитных подшипников, магнитов захваченного потока (МЗП). Размеры первых МЗП на основе объемных сверхпроводящих (СП) керамик [1, 2] ограничены ввиду хрупкости образцов. Размеры предложенных относительно недавно МЗП на основе стопок ВТСП лент [3, 4] также ограничены максимальной шириной промышленных лент. Один из вариантов реализации магнитной системы, позволяющий снять эти ограничения - замкнутые кольца. Кольца, изготовленные из ВТСП лент, более гибкие в размерах и позволяют изготавливать МЗП с большими габаритными размерами для промышленных применений. Кроме того, при помощи МЗП на основе ВТСП колец можно получить более высокие захваченные магнитные поля, что делает этот тип МЗП перспективным для применения в левитационных системах. Характеристики ВТСП колец зависят не только от критических параметров используемых материалов, но и от их геометрии (количества, диаметра и ширины колец), условий намагничивания и охлаждения. В данной работе приводятся результаты моделирования системы магнитной левитации на основе ВТСП колец из ленточного высокотемпературного сверхпроводника YBCO.
НИЯУ МИФИ, 115409 Россия, Москва, Каширское ш., 31; e-mail: mephizic@gmail.com.
Расчетная модель. Для того чтобы изготовить замкнутый СП контур, способный поддерживать созданное магнитное поле без внешних источников питания, предлагается выполнить разрез в центре каждой ВТСП ленты, как показано на рис. 1, и затем уложить ленты вместе, сформировав таким образом стопку ВСТП колец. Число колец варьировалось от 1 до 9. Выбирая различную длину разреза лент I, можно регулировать диаметр магнита й и, следовательно, его размер.
(а) (Ь)
Рис. 1: Процесс построения МЗП на основе ВТСП колец и геометрия левитационной системы. (а) Ь - длина ВСТП ленты, I - длина разреза, к-размер краевой области, т -ширина ленты, й - внутренний диаметр магнита; (Ь) общий вид расчётной области. На правой половине показана расчётная сетка элементов.
ВТСП ленты, из которых изготовлены СП кольца, в рамках модели обладают следующими характеристиками: ширина - 12 мм, толщина слоя подложки - 0.1 мм, толщина СП слоя - 1 мкм, толщина стабилизирующих слоев - 3 мкм, критический ток в собственном поле - 317 А при температуре 77.4 К. ВТСП-магнит по аналогии с работой [5] был изготовлен из 9 СП колец, сделанных из ленты длиной Ь = 150 мм, с длиной разреза I = 120 мм и размерами краевой области к =15 мм с каждой стороны от разреза. Внутренний диаметр магнита й при этом составил 50 мм. Кольцевой МЗП помещается в градиентное поле постоянного магнита (ПМ) с характеристиками, соответствующими №ЕеВ магниту марки N42 (максимальная индукция магнитного поля на поверхности ПМ составляет 0.4 Т). Намагничивание ВСТП колец в градиентном поле магнита рассмотрено в двух режимах охлаждения: в присутствии магнитного поля
(ЕС) и в нулевом магнитном поле ^ЕС). Геометрия расчетной системы представлена на рис. 1. Иллюстрация граничных условий для режимов ЕС и ZЕC приводится на рис. 2.
Рис. 2: Различие условий наложения внешнего магнитного поля для режимов охлаждения ЕС и ЕЕС. Я - комнатная температура.
Ранее проверенная нами формулировка в терминах компонент магнитного поля [6] была реализована в среде численного моделирования Сошяо1 МиШрЬувюБ для описания поведения сверхпроводника во внешнем магнитном поле. В рамках модели учитываются экспериментально измеренные угловые зависимости критического тока от магнитного поля, температурная зависимость критического тока лент [6], а также реальные тепловые характеристики всех слоев ленты (температурные зависимости теплоемкости и теплопроводности). Для соответствия реальным геометрическим размерам всех слоев ленты при построении модели применялись инструменты протяжки сетки элементов и многомасштабное структурирование (см. рис. 1).
Результаты и обсуждение. На рис. 3 представлено расчетное распределение г-компоненты индукции магнитного поля МЗП на высоте г = 10 мм в случае, когда все ВСТП кольца полностью намагничены и имеют одинаковое распределение токов. Так же как и в работе [5], в нашем случае наблюдается асимметрия в распределении поля, обусловленная геометрическим фактором.
Обозначим максимальное захваченное поле магнита при условии полной намагниченности всех 9 колец Втах = 26 мТ. Рассмотрим в качестве примера случай охлаждения сверхпроводника в отсутствие внешнего магнитного поля ^ЕС). Расчеты показали, что процесс намагничивания МЗП можно условно разделить на три участка в зависимости от величины внешнего приложенного магнитного поля. 1) В случае, когда внешнее приложенное магнитное поле меньше чем Втах, наблюдается очень низкая сте-
н °-02 , г 0.015 03 0.01 0.005 0
-0.005 -0.01 -0.015
0.04 0.035 0.03 0.025
-50 -40 -30 -20 -10 0 10 20 30 40 50
у-соогсШше, тт
Рис. 3: Профиль индукции магнитного поля МЗП на основе СП колец на высоте г =10 мм.
пень захвата магнитного поля, поскольку токи, возникающие на поверхности внешних ВТСП колец, практически полностью экранируют внешнее магнитное поле. Назовем этот участок зоной экранировки. 2) Когда приложенное поле больше чем Втах, начинается линейная зона, в которой захваченный магнитный поток линейно возрастает при увеличении внешнего магнитного поля. 3) Когда приложенное внешнее поле в четыре раза превышает Втах, все ВСТП кольца становятся полностью намагниченными, и дальнейшего роста захваченного потока при увеличении внешнего магнитного поля не происходит. Это так называемая зона насыщения. В случае ЕС для перехода МЗП в зону насыщения внешнее магнитное поле должно превысить Втах в два раза. Кроме того, было уставлено, что в краевой области ленты существует некоторая область "растекания" токов (рис. 4), размер которой зависит от характеристик конкретной ленты. Если размер краевой области к будет меньше области растекания токов, это приведет к локальному увеличению плотности тока и нагреву. Размер области растекания токов оценивался как расстояние от начала краевой области вдоль оси х (на высоте г, равной четверти ширины ВТСП ленты), на котором максимальный сверхпроводящий ток затухает в е раз, и для одного СП кольца размер области растекания токов составляет ~80% ширины ленты.
Также было проведено исследование влияния ширины ВТСП ленты и диаметра кольца на величину захваченного потока и силу левитации. Как правило, чем больше ширина ленты и чем меньше диаметр кольца, тем больше степень захвата потока (степень захвата потока определяется как отношение максимального захваченного поля к
Рис. 4: Распределение х-компоненты токов 3 в краевой области СП слоя ВТСП кольца.
приложенному внешнему магнитному полю). Для одиночного ВТСП кольца и ВТСП диска того же объема сверхпроводника и толщины 1 мкм (диаметр диска при этом составил ~48 мм, тогда как диаметр ВТСП кольца был равен 50 мм) были построены зависимости силы левитации от величины левитационного зазора между сверхпроводником и ПМ (рис. 5, слева). Сила левитации для ВТСП кольца существенно больше,
0.6-1-1-1-1-г
Displacement, mm Number of rings
Рис. 5: Зависимости силы левитации от величины левитационного зазора для ВТСП диска и ВТСП кольца (слева) и нормированные на максимальное значение зависимости отталкивающей силы и намагниченности от числа колец в МЗП.
следовательно, для достижения той же отталкивающей силы при использовании ВТСП дисков потребуется больший объем ВТСП материала. Зависимости намагниченности и силы левитации от количества колец, составляющих МЗП, носят линейный характер (рис. 5, справа). Аналогичные результаты наблюдались при исследовании магнитных и левитационных характеристик стопок ВТСП лент [7]. Приведенные зависимости являются линейными при варьировании числа СП колец от i до 9. Поведение данных зависимостей при дальнейшем увеличении числа колец требует дальнейшего изучения.
Таким образом, на основе метода конечных элементов была разработана комплексная модель магнитной левитационной системы на основе ВТСП колец. Выявлено существование трех характерных зон в процессе намагничивания лент. Установлено влияние геометрических параметров колец на величину захваченного магнитного потока и силу левитации. Для ВСТП диска и кольца проведено сравнение силы левитации, показаны преимущества МЗП на основе ВТСП колец.
Исследование выполнено при поддержке Российским научным фондом (проект № 17-19-01527).
ЛИТЕРАТУРА
[1] Z. Huang, M. Zhang, W. Wang, and T. A. Coombs, IEEE Trans. Appl. Supercond. 24, 4602605 (20i4).
[2] H. Matsuzaki, Y. Kimura, I. Ohtani, et al., IEEE Trans. Appl. Supercond. 15, 2222 (2005).
[3] A. Patel, K. Filar, V. Nizhankovskii, et al., Appl. Phys. Lett. 102, i0260i ^З).
[4] T. Tamegai, T. Hirai, Y. Sun, and S. Pyon, Physica C 530, 20 (20i6).
[5] Jie Sheng et al., Supercond. Sci. Technol. 30, 094002 (20i7).
[6] I. V. Anischenko, S. V. Pokrovskii, and I. A. Rudnev, IOP Conf. Series: Journal of Physics: Conf. Series 945, 0i20i5 (20i8).
[7] M. A. Osipov, D. A. Abin, S. V. Pokrovskiy, et al., Progress in Superconductivity and Cryogenics 17(i), 2i (20i5).
Поступила в редакцию i8 сентября 20i8 г.
После доработки i8 октября 20i8 г.
Принята к публикации 7 ноября 20i8 г.
Публикуется по результатам VII Международной молодежной научной школы-конференции "Современные проблемы физики и технологий" (16-21 апреля 2Q1S года, НИЯУ МИФИ, Москва).