Энергоемкость и энергетические потери индуктивного накопителя энергии на основе композитных ВТСП лент

Полный текст

ВВЕДЕНИЕ

Использование возобновляемых источников электроэнергии в последнее время вызывает интерес во всем мире [1, 2]. Однако, из-за прерывистого характера энергопроизводства большинства существующих возобновляемых источников энергии, таких как ветряные генераторы и солнечные батареи, накопители энергии являются важным компонентом любой устойчивой и надежной электросистемы. Наиболее широко исследованные системы накопления энергии включают в себя кинетические накопители энергии, индуктивные накопители энергии, гидроаккумулирующие системы, и суперконденсаторы [3, 4]. Они используются для обеспечения стабильности электросетей, благодаря их высокой способности к передаче и накоплению энергии [5–7]. Одними из перспективных устройств для накопления электроэнергии являются индуктивные накопители энергии на основе высокотемпературных сверхпроводников (ВТСП) [8].

Устройство СПИН представляет собой комбинацию трех технологий, а именно сверхпроводимости, криогеники и силовой электроники. Сверхпроводящая катушка – сердце СПИН. Электрически это чистый индуктор (без внутреннего сопротивления), и постоянный ток может течь через него без каких-либо омических потерь [9]. В результате сверхпроводящая катушка может сохранять ток или энергию в течение многих лет с плотностью энергии до 100 МДж/м³.

Накопленная энергия передается в электросеть в режиме разрядки. Индуктивные накопители способны быстро срабатывать и могут переходить из режима зарядки в режим разрядки за времена порядка нескольких микросекунд. Дополнительными преимуществами таких систем являются отсутствие движущихся компонентов и высокая циклическая эффективность [10]. Наиболее распространенными конструкциями катушек СПИН являются соленоид и тороид. При изготовлении катушек индуктивных ВТСП накопителей энергии используются такие технологии как CORC® кабель и рейстрек-катушки [11, 12].

Чтобы провод оставался в сверхпроводящем состоянии, ВТСП элементы индуктивного накопителя энергии всегда должны поддерживаться при температурах, плотностях тока и магнитном поле ниже критических значений. Для достижения и поддержания рабочей температуры используется криокулерная техника, а также жидкие хладогенты, такие как жидкий азот, гелий, неон и водород [13]. Критический ток в сверхпроводнике и, как следствие, энергоемкость СПИН зависят от температуры. Целью представленной работы является выбор оптимальной конфигурации индуктивного накопителя энергии, а также исследование зависимости максимально возможной накопленной энергии и потерь СПИН от температуры.

ТЕОРИЯ Н-ФОРМАЛИЗМА

Построение и расчет модели проводились в программной среде COMSOL Multiphysics®. В модели решается управляющее уравнение Фарадея методом конечных элементов:

$\nabla \times E+\frac{\partial B}{\partial t}=\nabla \times E+\frac{d\left({\mu }_{0}H\right)}{dt}=0$                                                       (1)

где  – напряженность электрического поля,  – индукция магнитного поля,  – напряженность магнитного поля,  – магнитная постоянная. Вольт-амперная характеристика сверхпроводника задается в виде степенного закона для нелинейной зависимости сопротивления от тока:

$\rho =\frac{E_c}{J_c}{\left(\frac{\left|J\right|}{J_c}\right)}^n$,                                                                                                      (2)

где E_c – критерий, принятый равным 1 мкВ/см, J_c – критическая плотность тока, n – показатель степени, характеризующий крутизну ВАХ сверхпроводника. Для корректного расчета в выражение (2) вводится зависимость критической плотности тока от величины приложенного магнитного поля J_c(B). Компоненты тока J=[J_x,J_y,J_z] вычисляются с помощью уравнений Максвелла:

$J=\nabla \times H=\left\{\begin{array}{c}\frac{d{H}_z}{dy}-\frac{d{H}_y}{dz}\\ -\frac{d{H}_z}{dx}+\frac{d{H}_x}{dz}\\ \frac{d{H}_y}{dx}-\frac{d{H}_x}{dy}\end{array}\right.$                                                                                    (3)

МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ТЕПЛОВЫХ ПРОЦЕССОВ В ВТСП

Тепловые процессы в ВТСП описываются управляющим уравнением теплопередачи, в котором температура Т является зависимой переменной:

$\rho \left(T\right)C_p\left(T\right)\frac{\partial T}{\partial t}+\nabla \cdot (-k(T)\nabla T)=Q+\rho \left(T\right)C_p\left(T\right)w\nabla T$,                         (4)

где ρ – плотность, k – коэффициент теплопроводности, Q – все источники тепла, C_p – теплоемкость при постоянном давлении, w – поле скорости температур. Термоэлектрический локальный источник тепла Q определяется выражением:

$Q=\left|\right|E\left|\right|\cdot \left|\right|J\left|\right|$,                                                                               (5)

где E и J определяются из уравнений (1)–(3). В выражении (2) в присутствии магнитного поля зависимость плотности критического тока  от температуры T и индукции магнитного поля B в ВТСП описывается экспериментальной интерполированной критической поверхностью, изображенной на Рис. 1 [14].

 

Рис. 1. Экспериментальная критическая поверхность

Fig. 1. Experimental critical surface

 

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

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

${\alpha }_{conv}=C_{conv}\mathrm{\Delta }{T}^{\frac{1}{3}}$,                                                                           (6)

где $C_{conv}$ – коэффициент, определяемый размером образца, $\mathrm{\Delta }T$ – разность температур на границе между ВТСП композитом и жидким азотом [15]. Для пузырькового кипения коэффициент теплоотдачи задается следующим образом:

${\alpha }_{{}_{boil}}=C_h{q}^{0.624}{\left(\rho C_{p}k\right)}^{0.117}$,                                                               (7)

где C_h – коэффициент, зависящий от материала поверхности ВТСП композита (в нашем случае медь) и параметров жидкого хладогента, С_р – теплоемкость, k – теплопроводность жидкого азота, ρ – плотность, q – тепловой поток. Выбор того или иного коэффициента теплоотдачи определяется, используя гистерезисную зависимость теплового потока q от разности температур $\mathrm{\Delta }T$ [15].

ОПИСАНИЕ ЭЛЕМЕНТОВ ИНДУКТИВНОГО НАКОПИТЕЛЯ

 

Рис. 2. Архитектура ВТСП-ленты SuperOx

Fig. 2. Architecture of HTSC tape SuperOx

 

Для создания замкнутых элементов в форме рейстрек и CORC® кабеля использовались композитные ВТСП-ленты REBa₂Cu₃O_7-x 2G фирмы SuperOx. На Рис. 2 приведена архитектура ВТСП-ленты. Такие ленты обладают высокими токонесущими характеристиками: критический ток ) в собственном поле при температуре кипения жидкого азота (77 К составляет 550 A. Толщина сверхпроводящего слоя составляет 1–3 мкм.

 

Рис. 3. Принципиальная схема CORC® кабеля

Fig. 3. Cable schematic diagram CORC®

 

На Рис. 3 представлена принципиальная схема CORC® кабеля. Конструкция каждого провода соответствует одной и той же базовой архитектуре: сверхпроводящие ленты наматываются спиральными слоями вокруг цилиндрического сердечника из мельхиорового сплава. На Рис. 4 представлена принципиальная схема рейстрек-катушки, в которой ВТСП ленты наматываются на формер послойно.

 

Рис. 4. Принципиальная схема рейстрек-катушки

Fig. 4. Schematic diagram of a racetrack coil

 

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Сравнение характеристик тороидальной и соленоидальной конфигураций

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

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

 

Рис. 5. Геометрия анализируемых конфигураций.

Слева – СПИН с тороидальной геометрией; Справа – СПИН с соленоидальной геометрией

Fig. 5. Geometry of analyzed configurations. Left: SPIN with toroidal geometry; Right – SPIN with solenoidal geometry

 

В обоих случаях сверхпроводящие катушки обладают следующими параметрами: диаметр провода – 5 мм, радиус кольца – 100 мм, рабочий ток – 5 кА. Длина ВТСП кабеля, необходимого для намотки: 7.34 м – для тороидальной конфигурации, 7.43 м – для соленоидальной конфигурации.

 

Рис. 6. Норма индукции магнитного поля и магнитные линии:

a – в соленоидальном СПИН; b – в тороидальном СПИН при рабочем токе 5кА

Fig. 6. Magnetic field induction rate and magnetic lines: a – in solenoidal SPIN; b – in toroidal SPIN at an operating current of 5 kA

 

На Рис. 6 и 7 представлены распределения величины индукции магнитного поля и линии магнитного поля в центральном сечении для обоих конструкций. Количество уровней для магнитных линий в обоих случаях одинаково. Охлаждение здесь производилось жидким азотом при температуре кипения 77 К.

Анализ Рис. 6 и 7 показал, что для соленоида рассеянное магнитное поле значительно больше, чем для тора (линии уровней в 1 случае отстоят от границы сверхпроводника максимум на 200 мм, а во 2 случае на 50 мм). Максимальное поле больше в торе почти на 60%. При этом максимум индукции магнитного поля в обеих конфигурациях приходится на сами ВТСП обмотки. Указанное обстоятельство может привести к ухудшению критических характеристик сверхпроводящей ленты в соответствии с Рис. 1 и к дополнительным потерям, при этом чем выше напряженность магнитного поля на обмотках ВТСП кабеля, тем выше будут гистерезисные потери. При разработке СПИН необходимо использовать высокополевую ВТСП ленту, способную сохранять критические характеристики в сильных магнитных полях.

На основе разработанных численных моделей был проведен расчет энергоемкости, механических напряжения и индуктивности тороидальной и соленоидальной конфигураций катушек СПИН. В Табл. 1 представлены основные расчетные характеристики индуктивных ВТСП накопителей энергии. Для секционной тороидальной конфигурации энергоемкость практически в 2 раза меньше, чем в соленоидальной конфигурации. Такая значительная разница в этом показателе объясняется более высокими гистерезисными потерями в случае использования тороида.

 

Таблица. Расчетные характеристики

Table. Design characteristics

	Тороид	Соленоид
Энергоемкость, Дж	157	301
Растягивающая нагрузка, кПа	221	198
Индуктивность, Гн	0,3	0,8

 

Полученные результаты позволяют заключить, что конфигурация соленоидной катушки СПИН выбирается для малого и среднего масштаба, поскольку она подвержена более низким механическим напряжениям и имеет большую энергоемкость. Для крупномасштабных СПИН критически важным является показатель рассеянного магнитного поля, поэтому для реализации такого устройства, несмотря на более высокие механические нагрузки, тороидальная конструкция катушки предпочтительнее в сравнении с соленоидом. Кроме того, разница в полной длине ВТСП ленты, использованной в этих конфигурациях оказывается значительной только на крупномасштабных системах, где тороидальная конфигурация оказывается более экономически выгодной [16]. В нашем случае длина сверхпроводящей ленты отличается незначительно. Таким образом, для реализации маломасштабного прототипа спин выбирается соленоидальная конфигурация ВТСП катушки, и дальнейшие расчеты энергоемкости и потерь в системе выполняются для нее.

Сравнение потерь в CORC® кабеле и рейстрек кабеле

Для ВТСП CORC® кабеля и рейстрек катушек в соленоидальной конфигурации ВТСП накопителя был проведен расчет общих потерь при температуре кипения жидкого азота. Результаты расчетов представлены на Рис. 7.

 

Рис. 7. Полные потери в рейстрек и CORC® кабеле в зависимости от плотности тока при 77 К

Fig. 7. Total loss in racetrack and CORC® cable as a function of current density at 77 K

 

Видно, что при малых плотностях тока (J/J_c < 0,5) потери для обеих конфигураций катушек слабо отличаются (менее чем на 10%). При более высоких токах крутизна зависимости потерь от плотности тока увеличивается и CORC® кабель становится значительно более выгодным. Этот факт объясняется тем, что рейстрек кабель из-за своей слоистости греется изнутри, что приводит к большим потерям. Более того, в случае рейстрек кабеля рост потерь усиливается начиная с J/J_c = 0,6, а у CORC® кабелей с 0,7.

Исходя из вышеизложенного можно сделать следующий вывод: для реализации прототипа индуктивного ВТСП накопителя энергии оптимальным является выбор CORC® кабеля в конфигурации соленоида. Оптимальное отношение J/J_c при охлаждении жидким азотом составляет 0,7.

 

Рис. 8. Полные потери в CORC® кабеле в зависимости от плотности тока при различных температурах

Fig. 8. Total loss in a CORC® cable as a function of current density at different temperatures

 

На Рис. 8 представлены результаты расчетов потерь для случая CORC® кабеля в соленоидальной геометрии при различных температурах. Как и ожидалось, потери снижается уменьшением температуры. При этом следует отметить, что для 77 К при J/J_c > 0,5 параметры теплоотвода жидким азотом ухудшаются и криокулер в данном случае лучше подходит для охлаждения. Однако, в то же самое время, криокулер тратит запасенную в СПИН энергию на охлаждение. Учет этой связи станет предметом последующих исследований.

ЗАКЛЮЧЕНИЕ

На основе метода конечных элементов разработана численная модель индуктивного ВТСП накопителя энергии в тороидальной и соленоидальной геометрии. Рассмотрены случаи конструирования ВТСП обмоток на основе рейстрек катушек и CORC кабеля. Проведена сравнительная оценка характеристик тороидальной и соленоидальной конфигураций СПИН с использованием CORC® и рейстрек кабелей, в результате которой, была выбрана оптимальная с точки зрения минимизации потерь соленоидальная конструкция маломасштабного индуктивного накопителя с проводом из CORC® кабеля. В случае охлаждения жидким азотом для выбранной конфигурации было получено оптимальное, с точки зрения энергоэффективности, отношение плотности тока к плотности критического тока, равное 0.7. Проведены расчеты потерь в геометрии соленоида, охлаждаемого при помощи криокулерной установки, при различных температурах.

БЛАГОДАРНОСТИ

Работа выполнена в рамках Госзадания (проект FSWU-2022-0013) при поддержке Министерства науки и высшего образования РФ.

Авторы заявляют, что:

1. У них нет конфликта интересов;
2. Настоящая статья не содержит каких-либо исследований с участием людей в качестве объектов исследований.

The authors state that:

1. They have no conflict of interest;
2. This article does not contain any studies involving human subjects.

Об авторах

Дмитрий Александрович Александров

Национальный исследовательский ядерный университет «МИФИ»

Автор, ответственный за переписку.
Email: cfrfcfrfdima123@gmail.com
ORCID iD: 0009-0001-7383-0094
SPIN-код: 5365-6190

инженер-исследователь

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

Ирина Валерьевна Мартиросян

Национальный исследовательский ядерный университет «МИФИ»

Email: mephizic@gmail.com
ORCID iD: 0000-0003-2301-1768
SPIN-код: 3368-8809

кандидат физико-математических наук, инженер-исследователь

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

Сергей Владимирович Покровский

Национальный исследовательский ядерный университет «МИФИ»

Email: sergeypokrovskii@gmail.com
ORCID iD: 0000-0002-3137-4289
SPIN-код: 6643-7817

кандидат физико-математических наук, заведующий научно-исследовательской лаборатории

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

Виктория Витальевна Залеткина

Национальный исследовательский ядерный университет «МИФИ»

Email: viktoriazaletkina@gmail.com
ORCID iD: 0009-0009-9854-5028

инженер-исследователь

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

Игорь Анатольевич Руднев

Национальный исследовательский ядерный университет «МИФИ»

Email: iarudnev@mephi.ru
ORCID iD: 0000-0002-5438-2548
SPIN-код: 2070-5265

доктор физико-математических наук, профессор, ведущий научный сотрудник

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

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

Дополнительные файлы