CC BY
14
Труды БГТУ, 2016, № 6, с. 147-151
147
УДК 531.19
Р. Н. Ласовский, Д. В. Гапанюк
Белорусский государственный технологический университет
МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ ЭЛЕКТРОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ТОКОПРОВОДЯЩЕЙ КЕРАМИКИ
Модифицированная система уравнений Нернста - Планка - Пуассона применяется для описания распределения электрического потенциала, напряженности электрического поля и заряда по пространству твердого электролита при стационарном переносе заряда. Свободная энергия рассматривается как функционал поля плотности и используется его разложение в ряд до квадратичных по градиенту слагаемых, вследствие чего выражение для химического потенциала дополняется градиентным слагаемым, и в уравнениях переноса появляется градиентная составляющая, вносящая определяющий вклад в области существования двойного слоя. Сформулирована система трех дифференциальных уравнений четвертого порядка для задачи о контакте трех твердотельных фаз и получено ее аналитическое решение, описывающее немонотонное распределение заряда. Для определения постоянных интегрирования используются условия сшивки, заключающиеся в требовании непрерывности функции, ее первой, второй и третьей производных на границах областей. Несмотря на наличие трех областей, распределения потенциала и заряда являются гладкими функциями ввиду принятых условий сшивки. Следует отметить, что полученное решение записывается сразу во всей области твердого электролита без его искусственного деления на отдельные зоны.
Ключевые слова: электропроводящая керамика, уравнения Нернста - Планка - Пуассона, химический потенциал, электропроводность, коэффициент диффузии, плотность заряда.
R. N. Lasovsky, D. V. Gapanjuk
Belarusian State Technological University
SIMULATION OF ELECTROPHYSICAL CHARACTERISTICS OF ELECTROCONDUCTING CERAMICS
The modified system of Nernst - Planck - Poisson equations is used to describe the distribution of electric potential, electric field and space charge of the solid electrolyte under a stationary charge transfer. The free energy is considered as the density field functional and its series decomposition up to quadratic in the gradient terms is used, resulting in the expression for the chemical potential that is complemented by a gradient term. The gradient component appears in the transport equations, making a decisive contribution to the area of a double layer. A system of three fourth-order differential equations for the problem of the contact of three solid phases is formulated, and its analytic solution describing nonmonotonic distribution of charge is obtained. To determine the constants of integration the cross-linking conditions requiring continuity of the electric field, its first, second and third derivatives on the boundary are used. Despite the presence of the three areas, the potential distribution and the charge density are smooth functions due to the conditions adopted by the cross-linking. It should be noted that the solution obtained is recorded in the entire region of the solid electrolyte without artificial division into separate zones.
Key words: conductive ceramics, Nernst - Planck - Poisson equations, chemical potential, electrical conductivity, diffusion coefficient, charge density.
Введение. В настоящее время наблюдается значительный интерес к области энергетики, связанной с твердотельными химическими источниками тока [1, 2]. Это обусловлено возросшей потребностью в портативных источниках электрической энергии различного назначения. Отмечается динамичное развитие рынка батарей для нужд мобильных устройств [3, 4].
В современных литиевых батареях в основном применяются жидкие полимерные материалы. Переход к твердым электролитам позволит повысить экологичность и безопасность источников энергии. Важной задачей
является теоретическое описание равновесного распределения и движения ионов в таких системах.
Одним из подходов считается решение системы уравнений Нернста - Планка - Пуассона (НПП), позволяющее устанавливать механизмы электропереноса в различных случаях, опираясь на методы электрохимического импеданса, и исследовать кинетические характеристики электрохимических систем при различных термодинамических условиях в широком диапазоне характеристических частот или характерных масштабов времени.
148 Моделирование распределения электрофизических характеристик токопроводяшей керамики
Ранее рассматривались варианты решения линеаризованных уравнений. В соответствии с ними суммарный поток частиц в среде определяется градиентом электрохимического потенциала и состоит из двух составляющих, порождаемых движущими силами, пропорциональными градиентам электрического потенциала и концентрации числа частиц. Этот подход оказался успешным при анализе миграции частиц в пределах одной фазы и при малой неоднородности распределения плотности в системе.
В данной работе рассмотрен вариант модификации уравнений НПП [5, 6] для описания распределения заряда по пространству твердого электролита при стационарном переносе заряда. При этом химический потенциал рассматривается не как функция, а как функционал поля плотности. Вследствие этого в уравнениях переноса появляется градиентная составляющая, вносящая определяющий вклад в области существования двойного слоя.
Система определяющих уравнений. Рассмотрим процесс стационарного переноса заряда в электрохимической среде при условии суммарной электронейтральности среды. При этом полагаем, что исследуется фрагмент среды, моделирующий токопроводящую керамику, состоящую из трех областей: приэлек-тродная, объемная фаза и межзеренная.
Неравновесное состояние среды описывается электрохимическим потенциалом [6]:
х) = M*ch + <?Ф + а
д 25р
дх2 '
(1)
где Цсь - химический потенциал; д - заряд частицы; ф - электрический потенциал; а - параметр взаимодействия; 5р - отклонение плотности заряда от равновесного значения.
Рассматривается случай одномерной неоднородной вариации плотности 5р, возникающей в среде под действием внешней разности потенциалов 2Дф, приложенной к электродам, расположенным на расстоянии 2Ь друг от друга вдоль оси х. Ширину межфазной области примем равной 21, а начало координат выберем по центру этой области. В условиях стационарно -сти, учитывая симметрию, можно утверждать, что распределение потенциала и плотности заряда вдоль оси х антисимметрично.
Для получения плотности распределения заряда выразим электрохимический потенциал через отклонение плотности:
Мх) = ^eq +Y5P( х) и используем уравнение Пуассона:
д 2ф
дх'
■ = -- 5р( х),
8
(2)
(3)
где - значение химического потенциала, соответствующее выбранному равновесному состоянию среды; у = дцеч / ф - термодинамический фактор; е - диэлектрическая проницаемость среды.
Соотношения (1)-(3) позволяют записать вариацию электрохимического потенциала в форме
ае д4ф ув д2 ф
5ц( х) =--—^ + дф-—(4)
д дх д дх
Уравнение для плотности тока имеет вид
J = -Ф—
дх
(5)
где D - коэффициент диффузии.
При стационарных условиях (J = const) соотношение (5) после интегрирования по х с учетом выражения (4) приобретает вид
Jx а8 (4) у8 (2) .
— =-ф()+^ф(2) - дф + A, (6)
qD q q
где A - постоянная, учитывающая граничные условия на левой границе участка (ось х предполагается направленной вправо, при этом рассматривается область х > 0, поскольку распределение потенциала электрического поля считается антисимметричным), а верхний индекс при Ф указывает порядок производной по пространственной переменной.
Из уравнения (6) следует, что ф удовлетворяет неоднородному линейному дифференциальному уравнению со специальной правой частью, которое для трех участков запишется в виде системы трех уравнений в стандартной форме:
Ф(4) + а1ф(2) + Ь1ф = d1 х, х < l, (7)
Ф(4) + а2ф(2) + Ь2ф = d1l + d2( х -1), х > l, (8) Ф(4) + а2ф(2) + Ь2ф = d1l + d2(L -1), х > L, (9)
где
a1 <0; b1 = -0_> 0, 8i
Y 2
(10)
a2 = < 0; b2 =—> 0, (11)
а
J
J
Djaj8j ^D 2 a 282
(12)
Индексы 1 и 2 обозначают номер участка в порядке следования уравнений (7)-(9), при-
2
Р. Н. Ласовский, Д. В. Гапанюк
149
чем параметры среды предполагаются одинаковыми для участков 2 и 3. Правые части этих уравнений сформированы определением константы А последовательным интегрированием левой части исходного уравнения (6) при условии, что 5ц(0) = 0:
А1 = 0 для 0 < х < 1,
Л
А2 =- для I < х < Ь,
яА
(13)
(14)
. Л Л(Ь -1) А3 =-+ —-- для х > Ь. (15)
яА яЦ2
Кроме того, в системе уравнений (7)-(9) учтено, что ток протекает между электродами, и, следовательно, на третьем участке ток равен нулю.
Ввиду характера правой части и граничных условий уравнение (7) должно удовлетворять следующему условию симметрии:
ф( х) = -ф(-х).
(16)
Корни характеристических уравнений (значения величин а и Ь выбираются согласно уравнениям (10) с параметрами, соответствующими рассматриваемому участку):
2
-Ь >0 при -а <"7"у. (17)
2 V 4
4я2
Следовательно, с учетом свойства симметрии (16) решение уравнения (7) должно содержать лишь нечетные функции координаты х:
ф1( х) = С^Ь^ х) + С2$Щ2 х) +—х. (18)
Ь1
В объемной области решение уравнения (8) выглядит следующим образом:
ф2( х) = С3е-к3( х-1) + С4вк3( х-1) + С5в-к4( х-1) +
. С ек4(х-1) + ^^ + ^2 (х - 1)
Ь1 Ь2
(19)
Также с учетом ограниченности потенциала на бесконечности решение уравнения (9) записывается в виде
ф3 (х) = С7 е~к3( х-Ь) + С8е-к4( х-Ь) +
+ + d2(Ь -1) Ь2 Ь2
(20)
Для определения постоянных интегрирования С (/ = 1, ..., 8) используются условия сшивки, заключающиеся в требовании непре-
рывности функции, ее первой, второй и третьей производных на границах областей, что дает 8 условий для 8 коэффициентов.
Решение в явном виде достаточно громоздкое и здесь не приводится, а конкретные вычисления выполнены в программном пакете МаШса«!
В результате распределения потенциала ф, напряженности поля Е, плотности заряда р выражаются зависимостями
ф( х) =
ф1 (х) при 0 < х < 1, ф2 (х) при 1 < х < Ь, ф3 (х) при х > Ь,
Е (х) = ^, р( х) =
dх dх
(21)
(22)
Функции (21) и (22) описывают распределение электрических характеристик при перетекании постоянного тока из точки х = -Ь в точку х = Ь.
Результаты вычислений и их анализ. На рис. 1 и 2 представлены зависимости электрического потенциала, напряженности электрического поля и плотности заряда от координаты, полученные по формулам (21) и (22) при к1 = 0,1, к2 = 0,3, к3 = 0,5, к4 = 0,6, ^ = 0,5, Ь\ = Ь2 = 2, d2 = 0,3dl, с различной шириной межфазной области.
На рисунках приведена лишь область х > 0, поскольку распределение потенциала электрического поля и заряда антисимметрично, а напряженности поля симметрично по пространственной координате.
Графики указанных зависимостей показывают, что, несмотря на наличие трех областей, распределения потенциала и заряда являются гладкими функциями ввиду принятых условий сшивки.
Следует также отметить немонотонный характер изменения плотности заряда в каждой из фаз при приближении к границе, их разделяющей (появление «крыльев»).
Как и следовало ожидать, ширина межфазной прослойки не влияет на ширину «крыльев» распределения заряда (которая определяется волновыми числами к) или быстроту изменения характеристик электрического поля по длине системы, а меняет лишь их местоположение, если только «крылья» не перекрываются. При этом напряженность поля имеет четко выраженные горизонтальные участки в областях нулевого заряда.
Локальные плотности заряда отличаются от нуля на границах участков вследствие скачков электрофизических характеристик среды либо присутствия электродов.
1 50 Моделирование распределения электрофизических характеристик токопроводяшей керамики
р 0,005 0,000 -0,005 -0,010 -0,015 -0,020
0 40 80 120 160 х
р 0,005 0,000 -0,005 -0,010 -0,015 -0,020 -0,025
0 40 80 120 160
х
б
б
в
Рис. 1. Зависимости электрического потенциала ф (а), напряженности электрического поля Е (б) и плотности заряда р (в) от координаты при I = 5 и Ь = 80
В целом система остается электронейтральной, поскольку в области межзеренной прослойки формируется двойной электрический слой, причем при достаточной ширине слоя положительная и отрицательная части заряда пространственно разделены (рис. 2), тогда как при малой ширине (по сравнению с обратными величинами волновых чисел к1 и к2) они образуют компактный двойной слой.
в
Рис. 2. Зависимости электрического потенциала ф (а), напряженности электрического поля Е (б) и плотности заряда р (в) от координаты при I = 50 и Ь = 80
Заключение. Выполнен анализ распределения плотности заряда в неоднородных многофазных системах, которыми являются керамические ионные проводники. В таких системах изменение распределения плотности заряда в каждой из фаз при приближении к границе, их разделяющей, должно быть немонотонным. Однако изучение решений линеаризованной системы уравнений (НПП) при
P. H. Ласовский, Д. В. Гапанюк
151
стационарных условиях показало, что они носят монотонный характер и, следовательно, эта система неприменима для описания вышеназванных систем.
Разработан вариант модификации уравнений НПП, устраняющий отмеченный недостаток. В результате в уравнениях переноса появляется градиентная составляющая, вносящая определяющий вклад в области межфазных слоев.
Сформулирована система трех дифференциальных уравнений четвертого порядка
для задачи о контакте трех твердотельных фаз и получено ее аналитическое решение, описывающее немонотонное распределение заряда.
Следует отметить, что полученное решение записывается сразу во всей области твердого электролита без его искусственного деления на отдельные зоны (приэлектродная область, область объема зерна, прослойка между зернами и т. п.) и четко отражает наличие неоднородно-стей в системе.
Литература
1. Armand M., Tarascon J.-M. Building better batteries // Nature. 2008. Vol. 451. P. 652-657.
2. Rechargeable Batteries: Grasping for the Limits of Chemistry / E. J. Berg [et al.] // Journ. Electrochem. Soc. 2015. Vol. 162, no. 14. P. A2468-A2475.
3. The Industry Review 2013. URL: http://fuelcelltoday.com/analysis/industry-review/2013/ the-industry-review-2013 (date of access: 27.09.2015).
4. Tarascon J.-M. Key challenges in future Li-battery research // Phil. Trans. Royal Soc. A. 2010. Vol. 368. P. 3227-3241.
5. Бокун Г. С. Асимптотическое решение уравнений Нернста - Планка - Пуассона для ионного кристалла // Труды БГТУ. 2014. № 6: Физ.-мат. науки и информатика. С. 32-34.
6. Бокун Г. С., Вихренко В. С. Химический потенциал в форме функционала плотности для описания импеданса электрохимической системы // Труды БГТУ. 2015. № 6: Физ.-мат. науки и информатика. С. 43-47.
References
1. Armand M., Tarascon J.-M. Building better batteries. Nature, 2008, vol. 451, pp. 652-657.
2. Berg E. J., Villevieille C., Streich D., Trabesinger S., Novak P. Rechargeable Batteries: Grasping for the Limits of Chemistry. Journ. Electrochem. Soc., 2015, vol. 162, no. 14, pp. A2468-A2475.
3. The Industry Review 2013. Available at: http://fuelcelltoday.com/analysis/industry-review/2013/the-industry-review-2013 (accessed 27.09.2015).
4. Tarascon J.-M. Key challenges in future Li-battery research. Phil. Trans. Royal Soc. A., 2010, vol. 368, pp. 3227-3241.
5. Bokun G. S. Asymptotic solution of the Nernst - Planck - Poisson equation for the ionic crystal. Trudy BGTU [Proceedings of BSTU], 2014, no. 6: Physical-mathematical sciences and informatics, рр. 32-34 (In Russian).
6. Bokun G. S., Vikhrenko V. S. Chemical potential in the form of density functional for description of inhomogeneous electrochemical system. Trudy BGTU [Proceedings of BSTU], 2015, no. 6: Physical-mathematical sciences and informatics, рр. 43-47 (In Russian).
Информация об авторах
Ласовский Руслан Николаевич - кандидат физико-математических наук, доцент кафедры теоретической механики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: [email protected]
Гапанюк Дмитрий Владимирович - кандидат физико-математических наук, доцент кафедры теоретической механики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13 а, Республика Беларусь). E-mail: [email protected]
Information about the authors
Lasovsky Ruslan Nikolaevich - PhD (Physics and Mathematics), Assistant Professor, the Department of Theoretical Mechanics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Gapanjuk Dmitry Vladimirovich - PhD (Physics and Mathematics), Assistant Professor, the Department of Theoretical Mechanics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Поступила 10.03.2016
