Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

3702

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
11.31 Mб
Скачать

где aj – расстояние между гранулами, а интеграл взят по всей площади грану-

лы.

Из второго уравнения Лондонов легко получить, что поле в грануле

H(x,y) H(x) H0

 

ch(x / )

 

 

 

 

.

(2)

ch(ax / 2 )

Тогда

ax aj

 

 

 

 

 

H0

 

 

 

He .

(3)

aj

 

sh(ax

/ )

 

 

ch(ax

/2 )

 

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

Из (3) следует, однако, что Н0 должно на порядки превосходить Не, их равенство возможнотолько при выполнении соотношения

ax / sh(ax / )/ch(ax /2 ).

(4)

Это не согласуется с результатами эксперимента, где не обнаружено резкого уменьшения Bc1 в каком-либо материале, что практически невозможно при таком жестком условии. Действительно, условие (4) утверждает, что глубина проникновения много больше толщины гранул. Но реально размеры гранул на один-два порядка превосходят глубину проникновения (десятки – сотни нанометров для иттриевых ВТСП).

Модель прямоугольных гранул, очевидно, является грубым приближением. Более строгий подход может быть построен, если рассмотреть структуру цилиндрических гранул во внешнем магнитном поле, меньшем первого критического. Согласно [5], сверхпроводящий цилиндр в продольном магнитном поле удовлетворяет уравнениям

B 2 rotrotB 0,

divB 0

(5)

с граничными условиями равенства полей на границе сверхпроводник – вакуум. Решение системы в цилиндрических координатах при дополнительном условии малости глубины проникновения есть

B Be

R/ rexp((r R)/ ).

(6)

Таким образом, индукция в гранулах практически отсутствует. Рассмотрим ячейку такой структуры, включающую в себя четыре цилиндрические гранулы. Общая площадь сечений сверхпроводящих цилиндров 4πr2 = 0.25πa2, где a – сторона охватывающего эту ячейку квадрата. Тогда из условия сохранения магнитного потока легко получить, что Н0 = Не(1–0.25π)–1. Внутреннее и внешнее поля в гранулированных сверхпроводниках могут отличаться только в несколько раз. Такой подход к структуре гранул учитывает наличие пор в кера-

50

мических ВТСП и позволяет объяснить близость значений поля зарождения вихрей Абрикосова Bc1 различных материалов.

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

Автор выражает благодарность Н.А. Андреевой, И.И. Соломахиной и М.А. Авдееву за помощь в изготовлении образцов и проведении некоторых измерений.

Литература

1.И.М. Шушлебин // Известия АН, сер. физическая. 57, 11,178 (1993).

2.В.Е. Милошенко, И.М. Шушлебин, О.В. Калядин. ФТТ, 48, 3, 403 (2006).

3.И. М. Шушлебин, О.В. Калядин // Вестник ВГТУ, 3, 78 (2007).

4.К.И. Кугель, А.Л. Рахманов // СФХТ, 4, 11, 2072 (1991).

5.В.В. Ломакин, В.В. Томчик // Физико-математическое моделирование систем: Материалы II Междунар. семинара, Воронеж, ВГТУ, 2005. C. 93-98.

THE FIRST CRITICAL FIELD OF A GRANULAR SUPERCONDUCTOR

I.M. Shushlebin

Voronezh State Technical University

Demagnetization effect in the high-Tc superconductor Y-Ba-Cu-O on the low critical field Bc1 was studied

Keywords: superconductor, ceramics, Abrikosov vortex, demagnetization, magnetic flux

51

УДК 538.915

РЕНТГЕНОВСКИЕ АБСОРБЦИОННЫЕ Sn L3 СПЕКТРЫ ОЛОВА И ОКСИДОВ ОЛОВА ПО ДАННЫМ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ

М.Д. Манякин, О.И. Дубровский, Е.Р. Лихачев, О.А. Чувенкова, С.Ю. Турищев

Воронежский государственный университет manyakin@phys.vsu.ru

Методом линеаризованных присоединенных плоских волн с помощью программного пакета Wien2k проведено исследование электронной структуры объемных кристаллов металлического олова и полупроводниковых моно- и диоксида олова. Вычислены Sn L3 спектры рентгеновского поглощения. Проведен сравнительный анализ результатов расчета с данными синхротронного эксперимента

Ключевые слова: олово, оксиды, электронная структура, моделирование, Wien2k, ЛППВ, XANES

Процесс окисления олова и образования оксидных фаз SnOx привлекает к себе значительное внимание и часто становится объектом научного исследования [1], поскольку получающиеся в результате материалы обладают большим потенциалом промышленного применения [2, 3]. При этом важно контролировать состав и структуру образовавшихся материалов на различных этапах окисления. С этой целью, как было показано нами ранее [4], может успешно применяться метод линейной комбинации для интерпретации спектров XANES [5] в ультрамягкой области, соответствующей Sn M4,5 краю поглощения. В работе [6] авторами была предпринята попытка проведения аналогичного анализа, но с использованием Sn L3 спектров поглощения, которая, однако, не позволила им сделать однозначного вывода о структуре изучаемого материала. Причиной этого является недостаточно полный набор эталонных спектров, применяемый при проведении процедуры построения линейной комбинации спектров.

Хотя Sn L3 спектры поглощения оксидов олова ранее становились предметом экспериментальных исследований, расчеты подобных спектров до настоящего момента, по-видимому, не проводились. В связи с этим целью данной работы было проведение первопринципного расчета XANES Sn L3 спектров для эталонных материалов – металлического олова, его моно- и диоксида. На рис. 1 представлены элементарные ячейки исследованных в работе материалов. При моделировании металлического олова β-Sn, тетрагонального монооксида олова SnO и тетрагональной фазы диоксида олова SnO2 (T) использовались кристаллические параметры, приведенные в базе данных [7]. Для орторомбического диоксида олова SnO2 (O) были взяты значения кристаллических параметров, рассчитанные в [8].

52

Рис. 1. Элементарные ячейки: a) Sn, b) SnO, c) SnO2 (T), d) SnO2 (O)

Расчет электронной структуры выполнялся методом ЛППВ, реализованным в программном пакете Wien2k. Метод ЛППВ позволяет проводить вычисления рентгеноспектральных характеристик твердых тел и надежно интерпретировать экспериментальные данные [9, 10]. Непосредственно для вычисления спектров XANES использовалось приближение остовной дырки [4, 11]. Результаты моделирования, сопоставленные с данными, полученными в ходе синхротронного эксперимента, приведены на рис. 2. Экспериментальные рентгеновские спектры L3 края поглощения олова были получены на канале DCM (Double-Crystal Monochromator) синхротрона SRC (Synchrotron Radiation Center)

Университета Висконсин – Мэдисон, г. Стоутон, США. Исследовались объёмные эталонные материалы фирмы Alfa Aesar – фольга металлического олова, порошки SnO и SnO2.

Видно, что для металлического олова и тетрагонального диоксида олова наблюдается очень хорошее совпадение результатов теории и эксперимента. Расчетный спектр Sn L3 для монооксида олова также достаточно хорошо согласуется с экспериментальным по форме и положению основных особенностей тонкой структуры. Однако, положение главного максимума в расчете (3947.4) и эксперименте (3950.5) различается на 3.1 eV, что, очевидно, связано с доокислением поверхности частичек, формирующих порошок SnO, исследованный экспериментально. Особый интерес представляет рассмотрение рассчитанного Sn L3-спектра для орторомбической модификации диоксида олова. Спектр орторомбического диоксида олова в области энергий ниже 3949 eV оказывается близок к спектру тетрагональной фазы SnO2. В спектре SnO2 (O) присутствуют характерные и для SnO2 (T) пики A при 3926.5 и B при 3935.4 eV. Однако далее

53

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

Рис. 2. Сравнение экспериментальных (exp.) и рассчитанных (calc.) XANES Sn L3-спектров проанализированных образцов

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

54

при анализе сложных реальных многофазных объектов системы оловокислород.

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-32-00860 мол_а, при финансовой поддержке РФФИ и правительства Воронежской области в рамках научного проекта № 16-42-360612 р_а.

Литература

1.Рябцев С.В., Чувенкова О.А., Канныкин С.В., Попов А.Е., Рябцева Н.С., Воищев С.С., Турищев С.Ю., Домашевская Э.П. // Физика твердого тела. – 2016. – №2. – С. 180–184.

2.Шапошник А.В., Сизаск Е.А., Корчагина С.Н., Звягин А.А., Рябцев С.В., Назаренко И.Н. // Сорбционные и хроматографические процессы. – 2014. – Т. 14. – №. 4. – С. 674–677.

3.Журбина И.А., Цетлин. О.И., Тимошенко В.Ю. // ФТП. – 2011. – Т. 45. – № 2. – С.

241–244.

4.С.И. Курганский, М.Д. Манякин, О.И. Дубровский, О.А. Чувенкова, С.Ю. Турищев, Э.П. Домашевская. // Физика твердого тела. – 2014. – T. 56, №9. – С. 1690–1695.

5.Manyakin M.D., Kurganskii S.I., Dubrovskii O.I., Chuvenkova O.A., Domashevskaya E.P., Ryabtsev S.V., Ovsyannikov R., Turishchev S.Yu. // Computational Materials Science – 2016.

V. 121. – P. 119–123.

6.H. Hulme, F. Baxter, R. P. Babu, M. A. Denecke, M. Gass, A. Steuwer, K. Noren, S. Carlson, M. Preuss // Corrosion Science. – 2016. – V. 105. – P. 202 -208.

7.http://database.iem.ac.ru/mincryst/

8.F. E. H. Hassan, S. Moussawi, W. Noun, C. Salameh, A.V. Postnikov. // Computational Materials Science. – 2013. – V. 72. – P. 86-92.

9.Переславцева Н.С., Курганский С.И. // Физика твердого тела. – 1999. – Т. 41, № 11.

С. 2075–2080.

10.Kurganskii S.I., Pereslavtseva N.S. // Physica Status Solidi (b). – 2000. – V. 218, № 2. – P. 431–439.

11.Манякин М.Д., Курганский С.И., Дубровский О.И., Чувенкова О.А., Домашевская Э.П., Турищев С.Ю. // Физика твердого тела. – 2016. – T. 58, № 12. – С. 2294–2298.

X-RAY ABSORPTION Sn L3 SPECTRA OF TIN AND TIN OXIDES

BY COMPUTATIONAL MODELLING

M.D. Manyakin, O.I. Dubrovskii, E.R. Likhachev,

O.A. Chuvenkova, S.Y. Turishchev

Voronezh State University

Electronic structure of Sn bulk crystal and tin oxides has been theoretically investigated by the WIEN2k package in the network of the linearized augmented plane-wave method. Sn L3 X-ray absorption spectra have been calculated. A comparative analysis of the calculation results with data from a synchrotron experiment have been performed

Keywords: tin, oxides, electronic structure, calculation, Wien2k, LAPW, XANES spectra

55

УДК 621.316:176.681.3

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

Ю.И. Дёмин

Северо-Кавказский государственный технический университет, г. Ставрополь yurij.dm12@mail.ru

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

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

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

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

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

56

туальные электронные устройства (ИЭУ), используемые в составе автоматизированных систем технологического управления сетей и АСУ ТП подстанций, должны обеспечивать быстрые синхронные измерения параметров режима электрической сети [1 - 3].

Перечисленным требованиям в максимальной степени соответствуют ИЭУ на основе технологии векторных измерений. Подобные устройства в последние годы наряду с традиционной областью применения WAMS все чаще стали применяться в АСУ ТП подстанций, в автоматизированных системах технологического управления нового поколения WACS, а также в качестве основных элементов управления и регулирования активно-адаптивных сетей [1, 3]. Более доступной альтернативой по отношению к ИЭУ на основе технологии векторных измерений являются ИЭУ с функциями синхронных измерений параметров режима электрической сети, например, ЭНИП - 2 и ЭНИП - 3 [2].

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

В качестве исходной, для рассмотрения, принята обобщенная идеализированная структура радиальной электрической сети (рис. 1). Выделено 4 характерных структур радиальных схем. Предположим, что во всех приемных пунктах (подстанциях) и источнике питания (ИП) установлены синхронизированные векторные измерители параметров режима нагрузок: тока (Ij), напряжения (Uj), активной (Pj) и реактивной (Qj) мощностям. Здесь j 0,1, 2,...,n – теку-

щие номера нагрузок и ИП при j = 0.

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

1.Суммируются измеренные мощности нагрузок (отдельно по активным

иреактивным мощностям)

S

Sk

 

 

3 Ik Uk.

(1)

 

k \ 0

 

k \ 0

 

 

2. В первом приближении находятся суммарные: потери мощности и сопротивления отдельно по активным и реактивным составляющим (рис. 2, б)

S

 

 

 

 

 

S0 S

Sxxmpk ;

 

 

 

 

k \{0}

 

(2)

 

S U02

 

 

*

 

 

 

Z

 

 

.

 

 

 

2

2

 

 

 

P0

Q0

 

 

 

57

где Sxxmpk – потери в стали трансформатора k-й нагрузки. Напряжение виртуальных эквивалентных шин

 

 

 

 

*

 

 

*

 

 

 

 

 

 

P X Q R

 

 

 

U0*

0

 

0

 

.

 

(3)

 

 

U0

 

 

 

 

 

 

 

 

 

 

 

 

3. Находим в первом приближении потери напряжения и сопротивления

от каждой нагрузки до виртуальных шин

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

U U0* U

 

Z S

 

 

,

 

 

 

 

 

 

 

 

 

 

U0*

(4)

 

 

U

U0*

 

 

 

 

*

 

 

 

 

 

 

Z

 

 

 

 

,

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где \{0}.

4. Уточняем сопротивление от k-й нагрузки на φ-й итерации. Для этого

проверим соотношение:

 

 

 

 

 

 

0 подключение -йнагрузкик виртуальной

 

 

 

 

шине через Z* ,

 

 

 

 

 

 

 

 

Z*

(Z*

Zтрмеди) 0 подключениек ИП,

 

(5)

 

 

 

 

0(с точностью ) подключение -йнагрузки

 

 

 

 

 

тр

 

 

 

 

 

к виртуальнойшине через Z

,

 

 

 

 

меди

где Zmpмеди

– потери в меди трансформатора -й нагрузки.

 

 

5.В соответствии с (5) производится: деление схемы на 4 группы, возврат

кначалу алгоритма (п. 1) и пересчет схемы без учета нагрузок 2 группы (подключенных непосредственно к ИП), см. рис. 1 (структура I).

6.Заменяем U0 в(2) на U0* из (3) на каждой итерации.

7.При совпадении с заданной точностью полученного напряжения с измеренным на µ-й нагрузке можно считать, что она не участвует в дальнейших разветвлениях и является «кустом»: фрагменты схем I и IV на рис. 1.

U0* U ,

(6)

где – заданная точность.

8.Условие невыполнения (6) потребует дальнейшей работы алгоритма: п.п. 1–6 при отсутствии нагрузок, для которых выполнены условия (6).

9.Реализация (2) предполагает расчеты параметров режима для активных

иреактивных составляющих:

*

Q U2

 

*

P U2

 

Х

 

0

;

R

 

0

.

(7)

P2

Q2

P2

 

 

 

 

Q2

 

 

0

0

 

 

0

0

 

 

58

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

Вкачестве иллюстрации работоспособности алгоритма проведены аналитические и машинные расчеты актуализации структуры и параметров для четырех характерных схем радиальных сетей (рис. 1, 3-6). Проверочные компьютерные расчеты выполнялись с помощью программы RERS PС [4]. Результаты расчетов показали высокую достоверность получаемых структур и параметров режима. В качестве примера рассмотрим расчеты одной итерации теста 1 (рис. 3). Параметры режима теста приведены в табл. 1.

 

 

 

 

Расчетные данные для теста 1 (рис. 3)

 

Таблица 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U,

Тип транс-

QН,

Рн ,

ΣΔPХХТР,

 

ΣΔQХХ,

ΣΔPМ,

ΣΔQМ,

 

 

 

кВ

форматора

квар

кВт

 

кВт

 

квар

кВт

 

квар

 

1

10,5

158,685

224,513

2,063

13,159

2,165

 

5,332

 

2

10,49

ТМ-250

65

 

 

100

 

 

 

3

10,483

ТМН-100

30

 

 

50

 

 

 

 

 

 

 

 

 

 

5

10,479

ТМ-1810

45

 

 

70

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

1

 

2

 

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

I I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

ИП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

III

2

1

1

2

IV

1

2

2

Рис. 1. Возможные структурные особенности радиальных разомкнутых сетей

59

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]