2D моделирование переноса 1:1 электролита в электромембранных системах при выполнении условия электронейтральности

Научный журнал КубГАУ, №110(06), 2015 года

1

УДК 519.85

01.00.00 Физико-математические науки

2D МОДЕЛИРОВАНИЕ ПЕРЕНОСА 1:1 ЭЛЕКТРОЛИТА В ЭЛЕКТРОМЕМБРАННЫХ СИСТЕМАХ ПРИ ВЫПОЛНЕНИИ УСЛОВИЯ ЭЛЕКТРОНЕЙТРАЛЬНОСТИ1

Коваленко Анна Владимировна к.э.н., доцент

РИНЦ SPIN-код автора: 3693-4813

Scopus Author ID: 55328224000

savanna-05@mail.ru

Кубанский государственный университет,

Россия,350040, Краснодар,Ставропольская, 149

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

Ключевые слова: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, 2D МОДЕЛИРОВАНИЕ, УРАВНЕНИЯ НЕРНСТА-ПЛАНКА-ПУАССОНА, КОНВЕКТИВНАЯ ДИФФУЗИЯ

UDC 519.85

Physics and Mathematical sciences

2D MODELING OF TRANSPORT 1: 1 ELECTROLYTES IN ELECTRO-MEMBRANE SYSTEMS WHEN HAVING THE CONDITION OF ELECTRONEUTRALITY

Kovalenko Anna Vladimirovna

Cand.Econ.Sci., associate professor

RSCI SPIN-code: 3693-4813

Scopus Author ID: 55328224000

savanna-05@mail.ru

Kuban State University, Krasnodar, Russia

The article presents a new approach to 2D modeling of transport of salt ions in EMC (electro systems: electrodialysis devices, electro-cells, etc.) under the condition of electrical neutrality with limiting and overlimiting current density. For definiteness as seen half of EMS channel EDA desalting (electrodialysis apparatus), the right border, which serves as a CEM (cation exchange membrane). The new approach in the use of partial differential equations of the first order, instead of equations of convective diffusion. A common method of transport modeling binary electrolyte in the EMS under the condition of electrical neutrality, is to use the equation of convective diffusion (partial differential equations of the second order). The article presents a new approach to modeling 2D transfer binary electrolyte in EMS under the same conditions, using partial differential equation of the first order for the decision, which does not require a boundary condition for concentration on the membrane surface. This allows you to simulate the transport of salt ions, as in prelimit and exorbitant current density and to determine the boundaries of the field of electrical neutrality

Keywords: MATHEMATICAL MODELING, 2D MODELING, NERNST-PLANCK-POISSON EQUATION, CONVECTIVE DIFFUSION

1 Работа выполнена при финансовой поддержке грантов РФФИ № 13-08-93105-НЦНИЛ_а, 13-08-96519 р_ юг_а и администрации Краснодарского края.

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

2

Введение

Для моделирования переноса бинарного электролита в электромембранных системах, как правило, используется система уравнений Нернста-Планка [1]. При выполнении условия электронейтральности из этих уравнений для бинарного электролита несложно получить уравнение конвективной диффузии [1] для определения концентрации. Однако использование этого уравнения при запредельных плотностях тока вызывает ряд сложностей. Ниже показано, что в этом случае удобнее не переходить к уравнению конвективной диффузии, а использовать для определения концентрации дифференциальные уравнения с частными производными первого порядка. Это позволяет при запредельных плотностях тока оценить область пространственного заряда.

1. Постановка задачи

1.1. Уравнения

Система уравнений Нернста-Планка и условия электронейтральности имеет безразмерный вид [1]:

ji = ziDiCiE - DtVC, + PeCV, i = 1,2

Pe = -divj, i = 1,2

at

z1C1 + z2C 2 = 0

1 = Z1j1 + Z2j2

Для простоты изложения рассмотрим стационарный перенос 1:1 бинарного электролита с D1 = D2 = 1, тогда уравнения (1)-(4) примут вид:

j = QE -V Cj + PeC V

j2 = -C2E — V C2 + PeC 2V

div j1 = 0

0 II . 1 ^

C1 - C 2 = 0

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

3

1 = j1 - j2

1.2. Геометрия

Рассмотрим половину канала обессоливания. Пусть х = 0 соответствует глубине раствора, где концентрация сохраняется постоянной, х = 1 соответствует условной границе раствор/катионообменная мембрана.

1.3. Г раничные условия

Ниже приведены общепринятые граничные условия для системы уравнений (5)-(10).

1) В глубине раствора при х = 0: С1 х=0 = С2| х=0 = 1.

2) На поверхности КОМ при х = 1: Q| х=1 = C1m, предполагается также, что КОМ идеально селективна, т.е. поток анионов j2 через катионообменную мембрану равен нулю j2| х=1 = 0.

3) На входе в канал при у = 0: cj у=0 = C2

у=0

4) На выходе из канала при у = L

ЭС

Эу

у=L

= 1.

0,i = 1,2.

Замечание 1. Кроме граничных условий 1)-4) на концентрацию, обычно задают следующие условия на потенциал р: р|х=0 = 0, — х=1 = р0 (т.е. прямая х = 0 и поверхность КОМ ( х = 1 ) являются эквипотенциальными), на входе в канал задается распределение р у=0 = р0х, на выходе «мягкое»

условие — у=L = 0 . Можно показать, что при допредельных плотностях

Эу у L

тока величины р0 и C1m связаны между собой (см.п.4.1). Указанные выше условия на р задают потенциостатический режим, однако в дальнейшем, когда будет удобно будем переходить к гальваностатическому режиму, поскольку между этими режимами существует однозначная связь. В гальваностатическом режиме задается средняя плотность тока icp по формуле

1 L

icp =—J Ь(х,уМу. L 0

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

4

2. Декомпозиция

1) Сложим первые два уравнения, тогда, обозначая S0 = C1 + C2 = 2C1, получим: v (C1 + C2) = (C1 - C2 )E + Pe(Cj + C2 )V - (j + j2) . С учетом условия электронейтральности, имеем: v S = PeSV - J, где J = j1 + j2.

2) Вычитаем первые уравнения, тогда:

v (C1 - c2) = (C1 + c2)E + Pe(Cj - c2)V -1 или se-1 = 0, т.е.:

E = 1J

S

Таким образом, для решения системы уравнений (1)-(10) достаточно решить следующую систему уравнений:

v S = PeSV - J .

3) Из уравнений divji = 0, i = 1,2, следует, что j., J, I - соленоидальные вектора.

3. Решение в одномерном случае

Если V = 0, то распределение концентрации не зависит от длины канала, поэтому достаточно рассмотреть распределение концентрации в некотором фиксированном сечении канала, т.е решать одномерное уравнение

(12). В одномерном случае имеем —S = -J, где J постоянная величина.

dx

Откуда S = -Jx + a, и, соответственно E =

1

- J, причем J = I = ic

Jx + a cp

для идеально селективной катионообменной мембраны. Находим постоянное а, используя условие при x =0, т.е. S|x=0 = 2. Таким образом, а = 2 = I [2], и функция S: S = -Ix +1 положительна лишь на отрезке

I

1

[0, ——) . Только на этом отрезке определена и функция E = — I.

IS

4. Анализ 2D задачи

Уравнение v S = PeSV - J, а также формула E = I /S, зависят от соле-ноидальных векторов I и J. Уравнения, не зависящие от соленоидальных

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

5

векторов I и J, можно получить, если взять div от обеих частей (11), (12), тогда получим уравнение конвективной диффузии для функции S:

AS = Pediv

и уравнение переноса для E : div(sE) = 0. Это уравнение обычно записывают для потенциала с учетом Е = -V j, тогда:

SAj +

Для решения уравнений (13), (14) необходимо ставить дополнительные граничные условия из-за того, что получились дифференциальные уравнение более высокого порядка.

4.1 Допредельный режим

Общепринятая практика, наряду с переходом к уравнениям (13),(14), заключается в том, что ставятся граничные условия при х = 1, например, условие:

SU = Sm

причем S1m(y) > 0, поскольку S является суммарной концентрацией. Из этого сразу следует, что такая задача моделирует лишь допредельный режим. Если же положить S|x=1 < 0, то задача для пары функций (S,E) не

имеет решение.

Рассмотрим для наглядности одномерный аналог уравнения (13):

d2 dx 2

S = 0, S|x=0

= 2,

S( 1) = Sm, причем S1m > 0.

Откуда S(x) = 2 - (2 - S1m)x или S(x) = lnp - (lnp - S1m)x. Это решение согласуется с решением s = -lx + inp, если положить inp - S1m = I. Отсюда получаем, что I < Inp, при S1m > 0. Если же S1m < 0, то существует точка х = Inp/I, где S(x) = 0 и функция S(x) = Inp - (Inp - S1m)x положительна лишь на отрезке [0,х). Только на этом отрезке определена и функция e = I/S. Таким образом, при запредельных плотностях тока отрезок [0,х) является

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

6

областью электронейтральности, а отрезок (х,1 ] областью пространственного заряда.

4.2. Запредельный режим в 2D. Сведение к задаче с неизвестной границей.

Дифференциальные уравнения (13), (14) при запредельных плотностях тока можно рассматривать как 2D краевые задачи в области с неизвестной правой границей и необходимо при этом доопределить краевые условия. Пусть нули функции S описывается формулой x = q(y), У е [0,L], т.е. S(q(y),y) = 0 для любого y e [0,L].

Обозначим S решение краевой задачи: AS = Pediv(sV) в области

W = {(х,У):У e (0,L),° < x < q(y)} : S x=0 = 2, S(q(y),y) = 0, AS(q(y),y) = 0,

Решения переопределенных задач с неизвестной границей является достаточно сложной задачей [3].

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

5. Решение вспомогательного уравнения v s = f в w при известном векторе f.

Рассмотрим решение уравнения v S = f в 2D односвязной области W при известном гладком векторе f .

5.1. Необходимое условие разрешимости

Определим двумерный аналог rot, - оператор r(a) = A^ -A;1. По-

S

y=0

y=L

= 0.

dx dy

скольку r

Э Э

S-А_Э

S = 0 для SeC2(W ), то необходимым услови-

Эx Эу Эу Эx

ем разрешимости является условие r( f ) =

r ) = f-f ’ Эx Эу

Л ЭЛ Э^

= 0 или = ^ll

Эx Эу

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

7

Если еще вектор f является и соленоидальным, то существует

Эh

dh

функция h, что — = -f1, — = f2, т.е. r(f) = ^ —2i = —2 + —2 = Dh, необ-

= Э/2 f = Э 2h + Э 2h =

Эу Эх

ходимым условием будет Dh = 0.

Эх Эу Эх2 Эу2

5.2. Достаточность и формула для решения

Покажем, что при известной вектор-функции f, удовлетворяющей условию разрешимости, решение v S = f определяется однозначно.

Запишем уравнение V S = f в координатной форме:

S = f1(x,y) , ^ S = f2 (х,у) .

Эх Эу

х

Решая первое уравнение, получаем: S = J f1(s,y)ds + у(у).

0

Подставляя это решение во второе уравнение имеем для у(у) урав-

х Э

нение: J—f^s^ds + У (у) = f2 (х,у). Используя условие разрешимости

0 Эу

f = Э1, получаем: Jdf2(s,у) + У(у) = f2(х,у).

Эх Эу 0

Поэтому f2 (х,у) - f2 (0,у) + у (у) = f2 (х, у). Таким образом, для у( у) имеем уравнение: У (у) = f2(0,у). Откуда следует, что:

у - х у

у(у) = J f2 (0,s)ds + const. Таким образом S = J f^s^)ds + J f2 (0,s)ds + const.

0 0 0

ух

Аналогично можно получить формулу: S = Jf2(х,s)ds + Jf/s,0)ds + const.

00

Формулу для S можно записать и симметрично:

1 х у 1 х у

S = - (J f/s^ds + J f2 (х,s)ds) + - ( J fx(s,0)ds + J fj (0,s)ds) + S o, где S0 = S( 0,0) -

2 0 0 2 0 0

некоторая постоянная. Таким образом, для однозначного решения системы уравнений v S = f достаточно условия S(0,0 ) = S0 .

6. Непосредственное решение системы уравнений

Перейдем к решения системы уравнений v S = PeSV - J, где J неко-

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

8

торый соленоидальный вектор. Поскольку J соленоидальный вектор, т.е.

divJ = 0, то существует функция g, что — = - J1, — = J2.

ду дх

6.1. Условие разрешимости

Обозначим f = PeSV - J . Необходимое и достаточное условие разрешимости уравнения V S = f для известной гладкой функции f в односвязной области WcR2, как показано выше, имеет вид r(f) = 0, где

r(f) = f-Ml . Таким образом, условие разрешимости уравнения (12) име-

дх ду

ет вид: r(PeSV - J) = 0. С учетом свойств r(f) имеем:

PeSr(V) + Pe(VS,V)-r(J) = 0. Так как r(J) = J-J = Dg, то получаем

дх ду

уравнение для функции g : Dg = Pe(VS,V)l + Pe  S  r(V).

Подставим в правую часть вместо VS правую часть (12), тогда Dg = Pe(PeSV - J,V) + Pe  S  r(V) или Dg= (Pe)2S(V,V)1 -Pe(J,V)L + Pe  S r(V) . Так как (V,V)l = 0 и (J,V)l = (Vg,V) , то Dg = -Pe(Vg,V) + Pe  S  r(V). К этому уравнению добавляем уравнение (12) и получаем систему из двух уравнений с двумя неизвестными g и S:

V S = PeSV - J ,

Dg=-Pe

6.2. Граничные условия

Наряду с граничными условиями на функцию S:

S|x=0 = 2, S|x=q(y) = 0, S|y=0 = 2, дУ

у=L

=0

необходимы граничные условия для функции g. Как следует из п. 5 условия (18) избыточны для нахождения решения уравнения (16), они должны рассматриваться как граничные условия для функции g. С учетом (16) и

(18) для вектора J получаем: J1 х=0 = 2 PeVi\ х=0

—I А

дх Iх=0 5 Jох=q

д ^

дх Iх=q(y) ,

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

9

J2 x=0 = 2PeV2x = 0

as, Эу

x = 0 ’

J

2x = q

as

Эу

j I = 2Pev\ a S1

x = q(yn j 1 у=0 2pevl\у = 0 dx

у = 0 ’

a s.

a s

J1\у=L = 2PeVl\у =L ax |y=L , J2\у = 0 = 2PeV2\у = 0 Эу

у=0 , J2 у=L = 2PeV2\у=L

a s a s

В силу условий (18) имеем — I = 0 = 0,

ax

эу

x=0

= 0. Откуда, с учетом

— = -J1, — = j2, получаем избыточное количество граничных условий

Эу dx

для функции у. ^1 x=0 = -2 PeVJ x=0 + a S' ag =a S

Эу

ay I

axU=0 2PeV2x=o, axix=q

ax |x=0 ’ Эу |x=q(}>) ax

ay a Si

x=q(y) ’

эу

x = q( у) ’

ay

ay i

эу

a Si

a Si

Эу |у=0 = 2PeV1 у=0, Эу Iу=L = 2PeV1 у=L + ax ^=L’

^ | = 2Pev I_______i

ax Iу=0 2Iу=0 Эу'у=05

у=L = 2 PeV2| у=L

Эх|у=L Д у=L

Из этих условий выбирается необходимое число граничных условий исходя из целей конкретного исследования. Как видно, из уравнения (17) и граничных условий (19) функция у определяется с точностью до постоянной величины. Определим конкретное значение функции у и свяжем ее с заданной плотностью тока. Для этого воспользуемся соотношением I = SE. Так как, по физическому смыслу I = SE > 0, и I = j\ - j2 , то j1 > j2 . Предположим, что j1 >> j2 . Это предположение выполняется в одномерном случае в силу предположения об идеальной селективности катионообменной мембраны, т.к. j2 = 0 и j = icp . В 2D оно будет заведомо выполняться,

например, при небольших скоростях прокачки раствора. При условии j1 >> j2 средний ток icp может быть определен следующей формулой

1 L 1 L

icp = — Jhix^d « — J J^x^d L 0 L 0

I j g рУ)ру = -1 (y(x,L) -y(x,0)) или L 0 Эу L

y(x,L) -y(x,0) = -Licp. Поскольку функция у определяется с точностью до постоянной величины, то можно положить, например, y(x, L) = 0, тогда

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

10

g(x,0) = Licp . Дополним эти условия следующими двумя условиями из (19):

|g|x=0 = 2PeV2\x_0 и g=q(y) = -|S|x=q(y) . ИтаK, для нахождения функциЙ g

и s имеем краевую задачу:

V S = PeSV - J ,

Dg=-Pe(Vg,V) + Pe • S • r(V) ,

S

x=0

= 2.

S

y=0 = 2, S x = q(y) = 0,

g(x, L) = 0, g(x,0) = Licp,

dg

dx

x=0

= 2 PeVJ x=0

dg I

axlx=q

d S,

dy <x=q

7. Решение краевой задачи (20)-(23).

7.1. Решение краевой задачи методом последовательных приближений.

Для решения краевой задачи (20) - (23) можно использовать следующий метод последовательных приближений:

1) Возьмем некоторый соленоидальный вектор J( 0);

2) Найдем решение S(0), уравнения V S = PeSV - J(0), удовлетворяющее граничным условиям при x = 0 и y = 0;

3) Найдем нули функции S(0), т.е. такую функцию x = q(0}(y), y e [0,L],

что S(0}(q(0}(y), y) = 0, y e [0, L].

4) Найдем решение g( 1 уравнения Dg = -Pe(Vg,V) + Pe  S(0) • r(V), в области W0 = {(x,y):0 < y < L,0 < x < q(0}(y)j, удовлетворяющее соответствующим

граничным условиям (23);

5) Найдем j( 1) и J 21) по формулам j( 1) = -

dg

(1) __ dg

и J21)

dg1

dy dx

6) Найдем решение S( 1) уравнения V S = PeSV - J( 1), удовлетворяющее соответствующим граничным условиям;

7) Проверим условие сходимости S( 1) -S(0) <d, где 8 заданная точность.

Если условие сходимости выполняется, S( 1) принимаем за решения, иначе

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

11

полагаем J( 0) = J( 1 и идем к п.2

7.2. Решение краевой задачи асимптотическим методом при небольших числах Пекле.

Если Pe << 1, т.е. скорость прокачки раствора небольшая, то функции g и S можно разложить по степеням Pe: g = g0) + Peg1) +...,

S = S( 0) + PeS( 1) +...

Подставляя эти разложения в систему уравнений (20), (21) и краевые условия (22), (23) и, приравнивая коэффициенты при одинаковых степенях в каждом уравнений и граничном условии, получаем краевые задачи для коэффициентов разложения. Для простоты изложения ограничимся приведем краевые задачи лишь для первых двух приближений.

Для g0), S( 0) получаем краевую задачу:

Л/0) = 0,

v s( 0) =- J( 0),

s( 0 г=0=2

S

у=0

g0)(x,L) = 0 , g0)(х,0) = Li ,

x=q

dg0;

Эх х 0 " Эх

Для g1, S( 1) получаем краевую задачу:

Dg1) =-(Vg0),V) + S(0) • r(V),

v s( 1) = - J( 1) + s( 0 V,

= 0.

Эg0)

Э S

x=q

Эу

x=q

S'1 )x=0 = 0, S<1 '|y=0 = 0, S'11

g1 '(x,L) = 0, g1'(x,0) = 0,

x=q

= 0,

Эу( 1

x = 0 _  2 x = 0 -1 ix=0 '

Эx Эx Эx

1) Решение краевой задачи для начального приближения g0), S

Рассмотрим два разных метода решения краевой задачи.

2V2 x=0 -

ЭУ0 -1, эг( 1'

x=0 ,

x=q

Э S ЭУ

x=q

0

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

12

Первый метод. Сначала решаем краевую задачу для уравнения V S = - J аналитически и находим S = G(J), где Г интегральный оператор (см. выше). Далее находим решение уравнения Dg = 0 в области W, с заранее неизвестной границей, определяемой дополнительным условием r(J) = 0 и граничным условием V r(J) = -J на этой границе.

Второй метод. Если для решения использовать метод последовательных приближений п.7.1, то получаем, что сначала задается некоторый соленоидальный вектор J(0). Затем находим S(0) =r(J(0)) и определяем градиент Vr(J(0)) и нули x = q(0}(y) функции S(0). Далее решаем уравнение Dg = 0 в области W0 = {(x,y):0 < y < L,0 < x < q(0)(y)j, удовлетворяющее соответствующим граничным условиям, и далее по алгоритму п.7.1.

2) Решение краевой задачи для начального приближения /1, S( 1).

Краевая задача для /1), S( 1) отличается от краевой задачи для g( 0), S(0) лишь тем, что уравнение для /1) неоднородное, а для S( 1 граничные условия однородные. Поэтому краевая задача для /1), S( 1) решается аналогично краевой задачи для g(0) , S(0).

8. Случай течения Пуазейля

8.1 Краевая задача для течения Пуазейля

Рассмотрим частный случай, когда течение раствора в канале является Пуазейлевским, т.е. V = (0,V2(x))T, где y = V2(x) = 3(1 - x2) - полупарабо-

ла Пуазейля. В этом случае имеем:

д_

dx

S = -J x(x,y),

Л

— S = Pe  SV2 (x) - J2 (x,y),

Лу

Dg = -PeV2(x)J 1 + PeSV2(x) .

8.2 Сведение краевой задачи к краевой задаче для инте-

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

13

гро-дифференциального уравнения

Из первого уравнения с учетом граничного условия S(0,y) = 2 для любого y е [0,L], получаем:

S = 2 - J J 1

0

С другой стороны, применяя общую формулу п.5, с учетом

f1 = - J1 ( х, у) и f2 = Pe  SV2 (x) - J2 (X,y), получаем

x y

S = - J J 1(s,y )ds + J (Pe  S(0,y )V2 (0) - J2 (0, y ))dy + a или

0 0

X y

S = - J J 1( s,y)ds + 2 PeV2 ( 0)y - J J 2 (0,y)dy + a,

00

где a постоянная, не зависящая от x,y. Это формула согласуется с формулой (37), поскольку в силу граничных условий:

y

a = 2 и 2PeV2 (0 )y - J J 2 (0,y)dy = 0.

0

Действительно, второе соотношение, как легко видеть после дифференцирования по y, эквивалентно граничному условию: 2PeV2\X=0 = X=0.

Подставляя (37) в уравнение (36) получаем для функции g краевую задачу для интегро-дифференциального уравнения:

X

Dg = -PeV2

0

с соответствующими граничными условиями из (19), (20), (21).

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

Краевую задачу для интегро-дифференциального уравнения (38) можно решать, например, методом последовательных приближений

X

Dg

0

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

14

в области Wк = {(x,y):0 < y < L,0 < x < q(k)(y)}, где J(0) некоторый соленои-дальный вектор, S(0) рассчитывается по формуле (37), x = q(0gy) - нули функции S( 0). Решая краевую задачу для (39) при к = 0 находим функции g1, J( 1), а затем функцию S( 1 и т.д.

8.4 Пример вычисления последовательных приближений

1) Положим J(0; = (J(0gy),J20}(x))T, где j(0)(y),J20}(x) некоторые

1 L

функции удовлетворяющие условиям: — J J 1( 0 gy)=i

L 0

cp :

и т. д.

2) Вычислим S

x.

3) Найдем область W0 = {(x,y):0 < y < L,0 < x < q(0J(y)} . Для этого положим 2 - J{0)(y)x = 0, следовательно x = 2/J(0)(y). Таким образом:

x = q( 0 }(y) = 2/J(i0)(y) . Область W 0 = {(x,y): 0 < y < L,0 < x < q(0)(y)} является начальным приближение области электронейтральности при заданной запредельной плотности тока icp > Inp = 2.

4) Решим уравнение: А/1) = -PeV2(x)j{0) + PeV2L(x)(2 - j(0)(y)x в области W0 = {(x,y):0 < y < L,0 < x < q(0}(y)} с заданными краевыми условиями:

g(x,L) = 0g(x,0) = Licp ,g=0 = 2PeV2\x=0 ,

dg

dx

d S

x=q

9y

x=q

например, численным методом конечных разностей находим g1), J( 1), а затем функцию S( 1).

Замечание 2. Если V2(x) = const = V0, уравнение (36) не зависит от S и

dg т~.

имеет вид Ag = - PeJ 1 или Ag = Pe—. Решая это уравнение с соответствую-

щими условиями, находим J и по формуле (37) находим функцию S .

Замечание 3. Метод 2D моделирования развитый здесь, несложно обобщается на случай произвольного бинарного электролита, и, более того, на случай произвольного число ионов.

Заключение. Общепринятый метод моделирования переноса би-

http://ej.kubagro.ru/2015/06/pdf/25.pdf

Научный журнал КубГАУ, №110(06), 2015 года

15

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

Работа выполнена при финансовой поддержке РФФИ и администрации Краснодарского края грантов: №13-08-93105-НЦНИЛ_а и 13-0896519 р_юг_а.

Литература

1. Ньюмен Дж. Электрохимические системы. М.: Мир. 1977. - 463с.

2. Коваленко А.В., Уртенов М.Х. Краевые задачи для системы электродиффузионных уравнений. Часть 1. Одномерные задачи. Germany, Saarbrucken: LAP LAMBERT Academic Publishing GmbH & Co. KG. 2011. - 281 c.

3. Вабищевич П.Н. Численные методы решения задач со свободной границей. М.: Изд-во МГУ, 1987. -164с.

References

1. N'jumen Dzh. Jelektrohimicheskie sistemy. M.: Mir. 1977. - 463s.

2. Kovalenko A.V., Urtenov M.H. Kraevye zadachi dlja sistemy jelektro-diffuzionnyh uravnenij. Chast' 1. Odnomernye zadachi. Germany, Saarbrucken: LAP LAMBERT Academic Publishing GmbH & Co. KG. 2011. - 281 c.

3. Vabishhevich P.N. Chislennye metody reshenija zadach so svobodnoj grani-cej. M.: Izd-vo MGU, 1987. -164s.

http://ej.kubagro.ru/2015/06/pdf/25.pdf