Физикохимия поверхности и защита материалов, 2020, T. 56, № 4, стр. 386-394

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

Ю. Г. Чирков 1*, В. И. Ростокин 2**, В. Н. Андреев 1***, В. А. Богдановская 1****, О. В. Корчагин 1*****

1 Институт физической химии и электрохимии им. А.Н. Фрумкина РАН
119071 Москва, Ленинский пр., 31, корп. 4, Россия

2 Национальный исследовательский ядерный университет (МИФИ)
115409 Москва, Каширское шоссе, 31, Россия

* E-mail: olga.nedelina@gmail.com
** E-mail: viktor.rostockin@yandex.ru
*** E-mail: vandr@phyche.ac.ru
**** E-mail: bogd@elchem.ac.ru
***** E-mail: oleg-kor83@mail.ru

Поступила в редакцию 13.03.2019
После доработки 27.09.2019
Принята к публикации 03.10.2019

Полный текст (PDF)

Аннотация

Активный слой (АС) катода литий-кислородного аккумулятора (ЛКА) должен иметь сложную структуру из пор двух типов (макропор и мезопор). Для успешной работы АС катода при разряде ЛКА стараются создать два кластера пор: кластер из макропор обеспечивает транспорт кислорода в зону, где происходит образование конечного продукта – пероксида лития, и кластер из мезопор, который гарантирует доставку ионов лития. В модельных расчетах проводится оптимизация этой структуры АС катода, состоящей из двух кластеров. Однако даже после реального воплощения этих теоретических рекомендаций, как свидетельствуют данные экспериментов, габаритные характеристики ЛКА при разряде (в частности, плотность тока i, мА/см2 и емкость катода С, Кл/см2) остаются низкими. В данной работе предложен новый тип структуры АС катода: регулярная бипористая модель. В такой модели каналы подачи в АС катода кислорода и ионов лития отделены друг от друга. Это позволяет одновременно и независимо улучшать функционирование каждого из двух каналов. Проведенные расчеты показали явное преимущество АС с новой структурой. В частности, плотность тока i и емкость катода С выросли до величин i ~ десятков мА/см2 и С ~ примерно до тысячи Кл/см2.

Ключевые слова: литий-кислородный аккумулятор, процесс разряда, активный слой катода, регулярная бипористая модель, компьютерное моделирование, константа процесса расхода кислорода k

ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ

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

Характерной особенностью процесса разряда катода ЛКА с неводным электролитом является “закупорка” пор положительного электрода нерастворимым в электролите и неэлектропроводным веществом – пероксидом лития Li2O2 [3, 4], конечным продуктом сложной реакции, происходящей в процессе восстановления кислорода [5, 6]:

(1)
$\begin{gathered} {{{\text{О}}}_{2}} + {\text{е}} \to {\text{О}}_{2}^{ - },\,\,\,\,{\text{О}}_{2}^{ - } + {\text{L}}{{{\text{i}}}^{ + }} \to {\text{Li}}{{{\text{O}}}_{{2~}}}, \\ 2{\text{Li}}{{{\text{O}}}_{2}} \to {\text{L}}{{{\text{i}}}_{2}}{{{\text{O}}}_{2}} + {{{\text{O}}}_{2}}. \\ \end{gathered} $

При заряжении ЛКА на аноде идет процесс регенерации кислорода:

(2)
${\text{L}}{{{\text{i}}}_{2}}{{{\text{O}}}_{2}} \to {\text{ }}2{\text{L}}{{{\text{i}}}^{ + }} + {{{\text{O}}}_{2}} + {\text{ }}2{\text{e}}{\text{.}}$

Существует много проблем, которые необходимо преодолеть для коммерческого применения ЛКА. Эти проблемы включают низкую эффективность процессов заряда и разряда, малую удельную емкость при разряде и реализуемую плотность энергии, невысокую циклируемость, большую длительность времени разряда, малую величину плотности тока. Все это в основном результат низкой эффективности работы АС катода ЛКА. Поэтому необходимо искать новые катодные материалы и специальные структуры АС, которые могли бы способствовать уменьшению катодного перенапряжения, особенно во время процесса разряда. К настоящему времени исследователями, разрабатывающими новые материалы для положительного электрода ЛКА, было рассмотрено огромное число композиций. Подробный анализ этих работ можно найти в обзорах [7, 8] и в работе [9].

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

В структуре пористого АС катода ЛКА имеет смысл выделить пять основных элементов, перечислим их. Во-первых, в катоде требуется присутствие макропор с радиусами в сотни нм. Они практически не “зарастают” пероксидом лития и потому не препятствуют прохождению через них кислорода. В [10] было показано, как быстро с ростом радиусов пор катода возрастает время окончания разряда, обусловленное “зарастанием” устья этих пор. Ввиду этого макропоры, если их сравнивать с микропорами и мезопорами, могут стать “газовым” каналом для подачи кислорода вглубь АС катода.

Во-вторых, в АС катода также необходимо присутствие мезопор. Они обладают высокой удельной поверхностью, необходимой для протекания процесса (1). При этом в отличие от микропор мезопоры имеют значительный собственный объем, необходимый для запасания в них значительного количества пероксида лития.

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

В-четвертых, в структуре АС катода должен быть сформирован с точки зрения теории перколяции [11] перколяционный кластер, способный обеспечивать путем диффузии в макропорах транспорт молекул кислорода в мезопоры. Кластер из макропор дальше будем называть “кислородным кластером”.

В-пятых, в структуре АС катода требуется создание и перколяционного “ионного кластера” из мезопор. Крайне важно подчеркнуть, что “зарастание” устьев мезопор ограничивает не только поступление кислорода из макропор, но также и поступление из них ионов лития. Поэтому можно с большой степенью точности полагать, что в основном ионы лития должны поступать в мезопоры не из макропор, а с “противоположной” стороны, из кластера заполненных электролитом мезопор, которые питаются в свою очередь ионами лития из межэлектродного пространства ЛКА.

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

Поскольку объем АС катода конечен и на 100% заполнен структурными элементами вышеперечисленных двух типов пор (макропорами и мезопорами), то при увеличении концентрации макропор характеристики кислородного кластера могут улучшиться, но при этом автоматически уменьшится доля объема АС катода, приходящаяся на материал с мезопорами. Это приводит к тому, что характеристики ионного кластера становятся хуже.

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

РЕГУЛЯРНАЯ БИПОРИСТАЯ МОДЕЛЬ

Рассмотрим предложенную нами модель структуры АС катода ЛКА. Она регулярна и бипориста. Изучаемая модель представлена на рис. 1. Она максимально проста. Имеется регулярный набор равновеликих крупных щелей толщиной δ (эти “макропоры” не “зарастают” пероксидом лития, поэтому это постоянно действующий канал для подачи в катод кислорода на большую глубину). Щели разделены между собой монопористыми слоями толщиной 2l, в мезопорах которых и происходит образование пероксида лития. С одной из сторон, назовем ее фронтальной, АС катода граничит с газовой фазой, с поступающим из нее в катод кислородом, с противоположной, тыльной стороны АС катода граничит с электролитной камерой.

Рис. 1.

Регулярная бипористая модель структуры активного слоя катода литий-кислородного аккумулятора. Имеет место регулярное чередование пористых слоев толщиной 2l, разделенных щелями размером δ. Толщина активного слоя L, его высота H = 1 см.

Все поры в катоде – и крупные, и мелкие – полностью заполнены электролитом – растворителем с ионами лития. В начальный момент процесса разряда ЛКА концентрация кислорода в электролите, находящемся в порах, нулевая. Затем начинается диффузионное проникновение кислорода в щели с последующим поступлением его в мезопоры пористых слоев. Будем полагать, что толщина пористого катода L, его высота H = 1 см (рис. 1).

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

Еще можно по желанию развивать поверхность встречи кислорода в щелях с мезопорами в монопористых слоях. АС катода теперь представляет собой набор однотипных элементарных ячеек – щели и два примыкающих к ним монопористых полуслоя.

Размер каждой ячейки δ + 2l, а поверхность в ней, на которой молекулы кислорода из щели могут попасть в поры полуслоя 2HL, таким образом, полная удельная поверхность, приходящаяся на 1 см2 видимой поверхности АС катода, равна, очевидно, 2HL/(δ + 2l). Если, допустим, толщина АС катода L = 1 см, H = 1 см, ширина ячеек δ = 1 мкм, толщина полуслоя l = 5 мкм, то удельная приведенная поверхность встречи “газовой” среды (пространство щелей) с мезопорами будет равным 2/(11 × 10–4) = 1818 см.

ОПИСАНИЕ ПРОЦЕССОВ В АКТИВНОМ СЛОЕ КАТОДА

В представленной на рис. 1 модели АС катода совершаются два процесса – диффузия растворенного в электролите кислорода по щелям и процесс образования пероксида лития в мезопорах монопористых слоев. Вывод уравнений, описывающих эти два процесса по отдельности, был выполнен в [12] и основан на системе уравнений.

Полагаем, что в начальный момент (при t* = 0) все щели в катоде и все поры в пористых слоях свободны от растворенного в электролите кислорода. Обозначим через а = а*/с0 – приведенную концентрацию кислорода в электролите в щелях, где с0 – растворимость кислорода в электролите. Тогда в области приведенных координат 0 ≤ х = = х*/L ≤ 1 в начальный момент в щелях

(3)
$а = {\text{ }}0\,\,\,\,{\text{при}}\,\,\,\,0 < х \leqslant 1.$

При этом на фронтальной поверхности АС катода в любой момент времени имеет место условие

(4)
$а = 1\,\,\,\,{\text{при}}\,\,\,\,х = 0.$

А на тыльной стороне АС катода (при х* = L) в щелях положим

(5)
${{dа} \mathord{\left/ {\vphantom {{dа} {dх}}} \right. \kern-0em} {dх}} = 0\,\,\,\,{\text{при}}\,\,\,\,х = 1.$

После проведения процесса обезразмеривания диффузионное уравнение примет вид

(6)
${{dа} \mathord{\left/ {\vphantom {{dа} {dt}}} \right. \kern-0em} {dt}} = {{{{d}^{2}}а} \mathord{\left/ {\vphantom {{{{d}^{2}}а} {{{d}^{2}}x}}} \right. \kern-0em} {{{d}^{2}}x}}--\tau I\left( {x,t} \right),$
где t = t*/t** – приведенное время, а t** = L2/D – характерное время для процесса диффузии в щелях, τ = 2L2/FδDc0, D – коэффициент диффузии кислорода в электролите, I(x, t), А/см2 – габаритный ток, характеризующий расход кислорода в полуслоях (он будет определен позднее). Таким образом, процесс диффузии кислорода в щелях существенным образом зависит от процессов, идущих в двух пористых полуслоях.

Теперь рассмотрим процессы, идущие в пористых слоях. Будем рассматривать монопористую модель пористых слоев в ее простейшем варианте. А именно полагаем, что в пористом слое толщиной 2l с начальной пористостью g0 имеется набор извилистых непересекающихся капилляров с постоянным начальным радиусом r0. Кислород (или воздух) поступает в два пористых полуслоя из щелей. Далее не очевидное, но важное замечание: считаем, что концентрация ионов лития в заполняющем поры катода электролите остается практически постоянной [1315].

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

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

Теперь введем для пор набор приведенных величин: приведенной координаты z = y/l, при этом координата z находится в пределах 0 ≤ z = y/l ≤ 1. Так как в поры кислород может поступать только из щели, то если в данном сечении щели, к которому примыкает пора, концентрация кислорода a = 0 (кислород отсутствует), то во всех сечениях рассматриваемой поры выполняются условия:

(7)
$r = {{r{\text{*}}} \mathord{\left/ {\vphantom {{r{\text{*}}} {{{r}_{{\text{0}}}}}}} \right. \kern-0em} {{{r}_{{\text{0}}}}}} = 1\,\,\,\,{\text{и}}\,\,\,\,с = {{с{\text{*}}} \mathord{\left/ {\vphantom {{с{\text{*}}} {{{с}_{{\text{0}}}}}}} \right. \kern-0em} {{{с}_{{\text{0}}}}}} = 0\,\,\,\,{\text{при}}\,\,\,\,{\text{ }}0 < z \leqslant 1,$
где r – приведенный радиус мезопоры (радиус поры r* с накоплением в поре пероксида лития уменьшается), с – приведенная концентрация кислорода. При этом на фронтальной (примыкающей к щели) поверхности поры для концентрации имеем условие

(8)
$c = а(х,t)\,\,\,\,{\text{при}}\,\,\,\,z = 0.$

А на тыльной стороне активного пористого полуслоя катода примем условие

(9)
${{dc} \mathord{\left/ {\vphantom {{dc} {dz}}} \right. \kern-0em} {dz}} = 0\,\,\,\,{\text{при}}\,\,\,\,z = 1.$

Итак, считаем, что

(10)
$t = {{t{\text{*}}} \mathord{\left/ {\vphantom {{t{\text{*}}} {t{\text{***}}}}} \right. \kern-0em} {t{\text{***}}}},\,\,\,\,r = {{r{\text{*}}} \mathord{\left/ {\vphantom {{r{\text{*}}} {{{r}_{{\text{0}}}}}}} \right. \kern-0em} {{{r}_{{\text{0}}}}}},\,\,\,\,z = {y \mathord{\left/ {\vphantom {y l}} \right. \kern-0em} l},\,\,\,\,c = {{с{\text{*}}} \mathord{\left/ {\vphantom {{с{\text{*}}} {{{с}_{{\text{0}}}}}}} \right. \kern-0em} {{{с}_{{\text{0}}}}}},$
где t*** = l2/D[g0]1/2 – характерное время для процесса образования пероксида в порах.

В системе двух взаимно согласованных процессов – диффузия с оттоком кислорода из щели и притоком кислорода в мезопоры, ведущего к образованию пероксида лития – мы имеем два характерных времени: большое для щели t** = L2/D и малое для пор t*** = l2/D[g0]1/2. Однако процессы в щели и в порах идут одновременно, поэтому естественно взять за основу более значительное по величине характерное время L2/D, полагая процесс образования пероксида лития в пористых полуслоях гораздо более быстрым процессом, нежели диффузия кислорода в щелях.

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

(11)
$d({{r}^{{4/3}}}{{c)} \mathord{\left/ {\vphantom {{c)} {dt}}} \right. \kern-0em} {dt}} = \alpha d({{{{r}^{2}}dc} \mathord{\left/ {\vphantom {{{{r}^{2}}dc} {dz}}} \right. \kern-0em} {dz}})--2\beta {{r}^{{1/3}}}c,$
(12)
${{dr} \mathord{\left/ {\vphantom {{dr} {dt}}} \right. \kern-0em} {dt}} = --\gamma c$
в них параметры α, β, γ имеют вид
(13)
$\begin{gathered} \alpha = {{({L \mathord{\left/ {\vphantom {L l}} \right. \kern-0em} l})}^{2}}g_{0}^{{1/2}},\,\,\,\,\beta = {{k{{L}^{2}}} \mathord{\left/ {\vphantom {{k{{L}^{2}}} {D{{r}_{{\text{0}}}}}}} \right. \kern-0em} {D{{r}_{{\text{0}}}}}}, \\ \gamma = k({{3{{L}^{2}}{{с}_{{\text{0}}}}} \mathord{\left/ {\vphantom {{3{{L}^{2}}{{с}_{{\text{0}}}}} {2D{{r}_{{\text{0}}}}}}} \right. \kern-0em} {2D{{r}_{{\text{0}}}}}})({{{{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}} \mathord{\left/ {\vphantom {{{{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}} {{{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}}}} \right. \kern-0em} {{{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}}}), \\ \end{gathered} $
где ${{{{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}} \mathord{\left/ {\vphantom {{{{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}} {{{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}}}} \right. \kern-0em} {{{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}}}$ – молярный объем, приходящийся на пероксид лития, ${{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}$, г/моль – его мольная масса и ${{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}$, г/см3 – плотность Li2O2.

Конечный итог расчетов – оценка изменений во времени t в зависимости от концентрации кислорода в щели а(х, t) величины габаритного тока I(х, t), его следует рассчитывать по формуле

(14)
$\begin{gathered} I(х,t) = k[{{4F{{g}_{0}}{{с}_{{\text{0}}}}l} \mathord{\left/ {\vphantom {{4F{{g}_{0}}{{с}_{{\text{0}}}}l} {{{r}_{{\text{0}}}}}}} \right. \kern-0em} {{{r}_{{\text{0}}}}}}]\int {({\text{от}}\,\,z = 0\,\,{\text{до}}\,\,z = 1)} \\ [{{r}^{{1/3}}}(x,z,t)c(х,z,t)]dz. \\ \end{gathered} $

А зависимость от времени количества накопленного в каждой элементарной ячейке заряда Q при окончании процесса разряда ЛКА можно рассчитать по формуле

(15)
$\begin{gathered} Q = {{2LH{{l}^{2}}} \mathord{\left/ {\vphantom {{2LH{{l}^{2}}} D}} \right. \kern-0em} D}{{[{{g}_{0}}]}^{{1/2}}}\iint {I\left( {х,\tau } \right)}{\text{ }}({\text{от}}\,\,t = 0\,\,{\text{до}}\,\,T) \\ d\tau \left( {{\text{от}}\,\,х = 0\,\,{\text{до}}\,\,х = 1} \right)dх, \\ \end{gathered} $
где Т – момент окончания процесса разряда, приведенное время разряда t* = Т/(l2/D[g0]1/2).

Также следует рассчитать, как со временем меняется распределение концентрации а (х) по толщине щели и как меняется со временем по толщине щели распределение габаритного тока I(х, t).

МЕТОДИКА РАСЧЕТОВ

При расчетах, проведенных в [12], в качестве растворителя был выбран тетраглим (ТГМ). Растворимость кислорода в нем с0 = 4.43 × 10–6 моль см–3, коэффициент диффузии D = 2.17 × 10–6 см2 с–1. Другие параметры: начальный радиус мезопор в пористых слоях r0 = 2 нм, толщина АС катода L = 0.07 см, начальная пористость слоев 0.73, ${{M}_{{{\text{L}}{{{\text{i}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{2}}}}}}}$, г/моль = = 45.8768 г моль–1, плотность пероксида лития ${{\rho }_{{{\text{Li}}}}}_{{_{{\text{2}}}{{{\text{O}}}_{{\text{2}}}}}}$, г/см3 = 2.3, ширина щелей δ = 1 мкм, половина толщины пористого слоя l = 5 мкм.

Описание того, как во времени ведут себя два сопряженных диффузионных процесса (система уравнений (6) и (11), (12)) – задача сложная. Характерное время процесса диффузии в щелях катода t** = L2/D на порядки больше характерного времени процесса диффузии в пористых слоях l2/D[g0]1/2. Поэтому при решении указанной системы уравнений для получения оценочных результатов можно прибегнуть к приближенному методу расчетов.

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

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

1. Выбирали достаточно малый интервал Δt1 безразмерного времени, ток полагали I = 0 (согласно начальным условиям), численно решали уравнение (6) и определяли распределение концентрации a(x, t1) вдоль щели, удовлетворяющее граничным условиям (4), (5) и начальному условию (3).

2. Пространство щели разделено на N областей, каждое шириной Δx = 1/N, в пределах каждой из которых a(xi, t1), взятое в середине данной области, задает граничное условие (17) на фронтальной поверхности полуслоя.

3. Для каждой цилиндрической поры, имеющей контакт со щелью в точке xi (i = 1…N), решаем уравнения (11)(12), предварительно преобразовав их к виду:

(16)
$\begin{gathered} \frac{\partial }{{\partial t}}{{c}_{i}}\left( {z,t} \right) = {{r}_{i}}{{\left( {z,t} \right)}^{{{2 \mathord{\left/ {\vphantom {2 3}} \right. \kern-0em} 3}}}}\frac{{{{\partial }^{2}}{{c}_{i}}(z,t)}}{{\partial {{z}^{2}}}} + \frac{2}{{{{r}_{i}}{{{\left( {z,t} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0em} 3}}}}}} \times \\ \times \,\,\frac{{\partial {{r}_{i}}\left( {z,t} \right)}}{{\partial z}}\frac{{\partial {{c}_{i}}\left( {z,t} \right)}}{{\partial z}} + \frac{{{{c}_{i}}\left( {z,t} \right)\left( {\beta - \frac{4}{3}\gamma {{c}_{i}}\left( {z,t} \right)} \right)}}{{{{r}_{i}}\left( {z,t} \right)}}, \\ \frac{\partial }{{\partial t}}{{r}_{i}}\left( {z,t} \right) = - \gamma {{c}_{i}}\left( {z,t} \right),\,\,\,\,\frac{\partial }{{\partial z}}{{r}_{i}}\left( {z,t} \right) = {{v}_{i}}(z,t) \\ \end{gathered} $
с граничными условиями:
(17)
$\begin{gathered} {{с}_{i}}(0,t) = a{\text{*}}(x_{i}^{*},t),\,\,\,\,\frac{\partial }{{\partial z}}{{с}_{i}}(1,t) = 0, \\ {{v}_{i}}(1,t) = 0\,\,\,\,(0 \leqslant t \leqslant {{t}_{1}}) \\ \end{gathered} $
и начальными условиями:

(18)
${{c}_{i}}(z,0) = 0;\,\,\,\,{{r}_{i}}(z,0) = 1;\,\,\,\,(0 \leqslant z \leqslant 1).$

4. Найдя распределения величин сi(z, t1) и ri(z, t1) в поре, вычисляем значение габаритного тока I(xi, t1) в каждой поре и путем сплайн-интерполяции определяем функцию I(x, t1) согласно формуле (14). На этом 1-й шаг вычислений завершается.

5. Далее, на каждом очередном интервале Δtn = = tntn – 1 решается уравнение (6) уже с учетом габаритного тока, а начальными условиями для амплитуд служат значения a(xi, 0) = a(xi, tn – 1), полученные на предыдущем шаге.

6. Затем находим распределение концентрации сi (z, tn) и радиусов пор ri (z, tn) в момент времени tn, решая систему уравнений (16) для каждой i-й поры с начальными условиями:

(19)
${{c}_{i}}(z,0) = {{c}_{i}}(z,{{t}_{{n - 1}}});\,\,\,\,{{r}_{i}}(z,0) = {{r}_{i}}(z;{{t}_{{n - 1}}}).$

Вычисляем габаритный ток I(x, tn) и повторяем численное решение для следующего интервала времени.

Когда радиус ri(0, t) устья поры становится меньше 0.1, считаем ее неактивной, не дающей вклад в габаритный ток, и потому исключаем ее из рассмотрения. Процесс вычисления габаритного тока и остальных параметров заканчивается, когда “закрывается” последняя (N-я) система пор. Это означает, что процесс разряда ЛКА завершен.

Численное решение дифференциальных уравнений в частных производных (6), (16) проводили с помощью математического пакета Maple [16, 17].

РЕЗУЛЬТАТЫ РАСЧЕТОВ

В работе [18] для АС катода ЛКА с растворителем пропиленкарбонатом была проведена оценка порядка величин параметра k. Результат – k = 1.025 × 10–8 см/с. Однако следует оговориться, что, скорее всего, приемы, примененные в [18] для оценки величины параметра k, непригодны для АС катода с регулярной бипористой моделью. Исходя из этого, для катода с регулярной бипористой структурой в работе [12] достаточно произвольно было выбрано три значения для параметра k: 10–7, 10–6 и 10–5 см/с. Результаты проведенных в [12] расчетов представлены в табл. 1. В первом столбце табл. 1 выбранный набор значений параметра k, во втором – величины соответствующих удельных емкостей, в третьем – продолжительность времени разряда ЛКА, примерно от 2 до 6 ч, очевидно, что при уменьшении k это время может только расти. И в четвертом столбце даны примерные плотности тока, необходимые для проведения разряда.

Таблица 1.  

Зависимость габаритных характеристик при разряде ЛКА с регулярной бипористой моделью активного слоя катода от константы процесса расхода кислорода k

k, см/с C, Кл/см2 t, с i, мА/см2
10–7 240.1 21136 16.1
10–6 397.7 9900 82.7
10–5 537.7 6500 152.7

Плотности тока в межэлектродном пространстве ЛКА i оценивали следующим образом. На рис. 2 показано, как во времени для катодов с тремя значениями параметра k в каждой отдельной элементарной ячейке АС катода с регулярной бипористой структурой растет габаритный ток I (количество образующегося в ячейках пероксида лития, формула (14)). Видно, что в начальные моменты разряда габаритный ток стремительно растет, так как по мере быстрого проникновения кислорода в щель в создание пероксида лития вовлекается все большее и большее количество пор в пористых слоях. Затем габаритный ток достигает своего максимума. А из-за последующего возрастающего “зарастания” пор в пористых слоях начинается неизбежный спад габаритного тока до нуля.

Рис. 2.

Зависимость от времени разряда катода ЛКА величины габаритного тока I в отдельной элементарной ячейке (щель плюс два примыкающих к ней пористых полуслоя). k, см/с: 1 – 10–5, 2 – 10–6, 3 – 10–7.

Выберем в качестве ориентира для оценки величины плотности тока i максимумы на кривых рис. 2. Затем принимаем, что i = Imax, толщина АС катода L = 0.07 см, а также то, что на квадратном сантиметре АС размещается 1/(δ + 2l) = 1/(1 + 10) × × 10–4 ячеек. С учетом этих допущений получим, например, для кривой 1 на рис. 2 (у нее максимум плотности тока равен ~2.4 мА/см2) такой результат: i = 2.4 × 10–3 × 0.07 × 1/11 × 10–4 = 152.7 мА/см2. Данные последнего столбца в таблице 1 показывают, что АС катода со структурой регулярной бипористой модели способен значительно (на порядки) увеличить плотности тока в ЛКА.

Все главные параметры регулярной бипористой модели – толщина АС L = 0.07 см, ширина щели δ = 1 мкм, толщины пористых полуслоев l = = 5 мкм, начальные радиусы пор в пористых слоях r0 = 2 нм, а также значения параметра k при расчетах, результаты которых были собраны в табл. 1, выбраны нами произвольно. Представляется, что наибольшее для функционирования АС катода ЛКА значение имеют два параметра структуры АС – радиус пор в монопористых слоях r0 и толщина пористых полуслоев l. Будем варьировать эти два параметра, оставив неизменными значения для величин L, δ и положив вначале для определенности, что k = 10–6 см/с. Это необходимо для оценки влияния радиусов пор и толщины слоя на величину удельной емкости С, Кл/см2 при разряде ЛКА.

Таким образом, так же, как это было при расчетах, выполненных в [12], вначале принимаем l = 5 мкм, а r0 будет произвольно варьироваться в окрестностях значения r0 = 2 нм. Затем проведем расчеты для емкости катода С, Кл/см2. Результат этих расчетов представлен на рис. 3. Видно, что варьирование радиуса пор при фиксированном значении параметра l = 5 мкм дает возможность для роста емкости. Для радиуса r0 = 2 нм на рис. 3 емкость примерно такая, как и в табл. 1 – 395.36 Кл/см2, а вот для r0 = 1.3 нм емкость возрастает до величины С = 463.94. Видно, что при отклонении радиуса пор от значения r0 = 1.3 вправо или влево емкость катода начинает быстро снижаться.

Рис. 3.

Регулярная бипористая модель. Зависимость емкости активного слоя катода ЛКА от размеров начального радиуса пор r0 в пористых слоях. Остальные параметры активного слоя катода: L = 0.07 см, δ = 1 мкм, k = = 10–6 см/с, l = 5 мкм.

Будем по-прежнему варьировать величину радиуса и дополнительно изменим толщину АС, пусть l будет на этот раз не 5 мкм, как в [12], а 15 мкм (рис. 4). В этом случае на С, r0 – кривой максимум наблюдается при r0 = 2.5 нм. При достижении этого максимума имеется более высокое значение емкости С = 738.23 Кл/см2, чем С = = 463.94 Кл/см2, как это показано на рис. 3.

Рис. 4.

Регулярная бипористая модель. Зависимость емкости активного слоя катода ЛКА от размеров начального радиуса пор r0 в пористых слоях. Остальные параметры активного слоя катода: L = 0.07 см, δ = 1 мкм, k = = 10–6 см/с, l = 15 мкм.

Рассмотрим, к чему приведет увеличение толщины пористого полуслоя. Увеличим ее до l = 25 мкм. И вновь будем варьировать величину начального радиуса пор. Результат этих новых расчетов представлен на рис. 5. Кривая на рис. 5 имеет максимум, на этот раз уже в точке r0 = 6 нм. Однако максимальная величина емкости на рис. 5 по сравнению с данными рис. 4 стала ниже и составляет С = 685.58 Кл/см2 (С = 738.23 Кл/см2 на рис. 4). Чтобы исключить негативное снижение емкости, изменим тактику подбора параметров r0 и l. Зафиксируем величину начального радиуса пор r0 и будем наращивать толщины пористых слоев l.

Рис. 5.

Регулярная бипористая модель. Зависимость емкости активного слоя катода ЛКБ от размеров начального радиуса пор r0 в пористых слоях. Остальные параметры активного слоя катода: L = 0.07 см, δ = 1 мкм, k = = 10–6 см/с, l = 25 мкм.

Будем считать, что r0 = 2.5 нм (точка максимума кривой на рис. 4), а толщину пористых слоев будем постепенно увеличивать. Результат этих расчетов представлен на рис. 6. В этом случае кривая имеет максимум в точке l = 30 мкм, который свидетельствует, что дальнейшее увеличение толщины АС катода ЛКА уже не имеет смысла. Кроме того, видно, что на рис. 6 при l = 15 мкм мы имеем, как это и должно было быть, то же значение емкости, которое было максимумом на рис. 4. Главная особенность рис. 6 состоит в том, что для емкости АС катода ЛКА на этот раз получено наибольшее значение – С = 1117.0 Кл/см2.

Рис. 6.

Регулярная бипористая модель. Зависимость емкости активного слоя катода ЛКА от размеров полутолщины пористых слоев l. Остальные параметры активного слоя катода: L = 0.07 см, δ = 1 мкм, k = = 10‒6 см/с, r0 = 2.5 нм.

До сих пор при расчетах предполагали, что параметр k = 10–6 см/с. Увеличим эту величину на порядок, до значения k = 10–5 см/с, и проведем расчет на одном примере при l = 25 мкм, а величину начального радиуса пор r0 будем произвольно варьировать. Результаты расчетов представлены на рис. 7. Вновь наблюдается кривая с максимумом, который свидетельствует, что нет необходимости для дальнейшего увеличения радиуса пор. Удивительным на рис. 7 является то, что максимум емкости теперь приходится на очень большие радиусы мезопор r0 = 45 нм! И величина удельной емкости опять принимает очень большое значение. При r0 = 45 нм и l = 25 мкм емкость составляет 956.42 Кл/см2. Величины получаемых при разряде ЛКА емкостей, видимо, можно еще значительно увеличить. Действительно, не существует никаких ограничений на толщину АС катода L. До сих пор мы полагали L = 700 мкм. Видимо, если взять величину L = 1 см, скорее всего, величина емкости при этом возрастет.

Рис. 7.

Регулярная бипористая модель. Зависимость емкости активного слоя катода ЛКА от размеров начального радиуса пор r0 в пористых слоях. Остальные параметры активного слоя катода: L = 0.07 см, δ = 1 мкм, k = = 10–5 см/с, l = 25 мкм.

Проведенные расчеты указывают на необходимость оптимизации набора из пяти параметров – L, δ, l, r0 и k. Это самостоятельная большая и трудная задача. Не только потому, что она потребует большого расхода машинного времени, но может возникнуть ряд ограничений на величину параметров, подлежащих оптимизации. К примеру, для обеспечения больших потоков ионов лития будут необходимы значительные величины толщин пористых слоев (регулярная бипористая модель), а получаемая при оптимизации величина l будет рекомендовать не столь значительные толщины пористых слоев для получения больших величин емкости катода при разряде.

ЗАКЛЮЧЕНИЕ

Для создания АС катодов ЛКА испробовано огромное количество композиций различных материалов. Однако пока экспериментально полученные габаритные характеристики ЛКА все же продолжают оставаться низкими. Ведутся и попытки компьютерного моделирования процессов в ЛКА. Результаты модельных расчетов показывают, что в целях оптимизации процесса образования пероксида лития необходимо выполнение ряда условий. В [9], к примеру, прямо указывается, что для структуры АС катода надо стремиться к увеличению общей пористости. Общая удельная поверхность пор должна составлять величину 3.5–4.5 × 106 см–1. В распределении пор по размерам должны присутствовать два пика групп пор с размерами 7–8 нм для мелких пор и 33–110 нм для пор крупных.

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

Поэтому и предлагается для АС катода ЛКА новая модель структуры – регулярная бипористая модель. Ее основное достоинство состоит в том, что в ней каналы подачи в АС катода кислорода и ионов лития разъединены, никак не зависят друг от друга. И это позволяет совершенствовать характеристики сразу двух каналов, независимо. Что все это дает? В данной статье после проведения серии компьютерных расчетов показано, что удается увеличить плотность тока от величин порядка десятых долей мА/см2 (при прежних методах создания структур АС катодов) до десятков мА/см2. И поднять емкость катодов ЛКА при разряде в АС катода с регулярной бипористой моделью от прежних десятков Кл/см2 до величин С ~ ~ 1000 Кл/см2.

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

  1. Abraham K.M. // J. Electrochem. Soc. 1996. V. 143. P. 1.

  2. Bruce P.G. // 2012. V. 11. P. 19.

  3. Tran C. // J. Power Sources. 2010. V. 195. P. 2057.

  4. Yang X.-h., He P., Xia Y.-y. // 2009. V. 11. P. 1127.

  5. Laoire C.O., Mukerjee S., Abraham K.M. // J. Physical Chemistry C. 2009. V. 113. P. 20127.

  6. Laoire C.O., Mukerjee S., Abraham K.M. // J. Physical Chemistry C. 2010. V. 114. P. 9178.

  7. Ma Z., Yuan X., Li L., Ma Z-F., Wilkinson D.P., Zhang L., Zhang J. // Energy Environ. Sci. 2015. V. 8. P. 2144.

  8. Тарасевич М.Р., Андреев В.Н., Корчагин О.В., Трипачев О.В. // Физикохимия поверхности и защита материалов, 2017. Т. 53. С. 3. [Tarasevich M.R., Andreev V.N., Korchagin O.V., Tripachev O.V. // Prot. Met. Phys. Chem. Surf. 2017. V. 53. P. 1].

  9. Pan W., Yang X., Bao J., Wang M. // J. Electrochem. Soc. 2017. V. 164. P. E3499.

  10. Чирков Ю.Г., Андреев В.Н., Ростокин В.И., Богдановская В.А. // Альтернативная энергетика и экология. 2018. № 4–6. С. 95. [Chirkov Yu.G., Andreev V.N., Rostokin V.I., Bogdanovskaya V.A. // Alternative Energy and Ecology (ISJAEE) (In Russian). 2018. № 4–6. P. 95].

  11. Тарасевич Ю.Ю. Перколяция: теория, приложения, алгоритмы, М.: Эдиториал УРСС, 2011. 112 с. [Tarasevich Yu.Yu. Percolation: Theory, Applications, Algorithms (In Russian). Moscow: Editorial URSS, 2011].

  12. Чирков Ю.Г., Андреев В.Н., Ростокин В.И., Богдановская В.А. // Электрохимия. 2019, в печати. [Chirkov Yu.G., Andreev V.N., Rostokin V.I., Bogdanovskaya V.A. // Russ. J. Electrochem., 2019, In press.]

  13. Sandhu S., Fellner J., Brutchen G. // J. Power Sources. 2007. V. 164. P. 365.

  14. Dabrowski T., Struck A., Fenske D., Maaß P., Colombi Ciacchi L. // J. Electrochem. Soc. 2015. V. 162. P. A2796–A2804.

  15. Read J., Mutolo K., Ervin M., Behl W., Wolfenstine J., Driedger A., Foster D. // J. Electrochem. Soc. 2003. V. 150. P. A1351.

  16. Эдвардс Ч.Г., Пенни Д.Э. Дифференциальные уравнения и краевые задачи: моделирование и вычисление с помощью Mathematica, Maple и MATLAB, 3-е издание. Киев: Диалектика-Вильямс, 2007 [Edwards C.H., Penny D.E. Differential Equations and Boundary Value Problems: Computing and Modeling (3rd Ed.), Мoscow, 2008].

  17. Голоскоков Д.П. Уравнения математической физики. Решение задач в системе Maple, С-Пб: Питер, 2004. 539 с. [Goloskokov D.P. Equations of mathematical physics. Solving problems in the Maple system (In Russian). St. Petersburg: Peter, 2004. 539 p.].

  18. Богдановская В.А., Андреев В.Н., Чирков Ю.Г., Ростокин В.И., Емец В.В., Корчагин О.В., Трипачев О.В. // Физикохимия поверхности и защита материалов. 2018. Т. 54. С. 549. [Bogdanovskaya V.A., Chirkov Yu.G., Rostokin V.I., Yemetz V.V., Korchagin O.V., Andreev V.N., Tripachev O.V. // Prot. Met. Phys. Chem. Surf. 2018. V. 54. P. 1015].

Дополнительные материалы отсутствуют.