CC BY
53
Научный журнал КубГАУ, №103(09), 2014 года
1
УДК 517.958: 531.32, 544.6
ПЕРЕНОС ИОНОВ СОЛИ В ЭЛЕКТРОХИМИЧЕСКОЙ ЯЧЕЙКЕ С ВРАЩАЮЩИМСЯ МЕМБРАННЫМ ДИСКОМ С УЧЕТОМ ЭЛЕКТРОКОНВЕКЦИИ. ЧАСТЬ 1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ1
Коваленко Анна Владимировна к.э.н., доцент
Уртенов Махамет Али Хусеевич д.ф.-м.н., профессор
Казаковцева Екатерина Васильевна аспирантка
Кубанский государственный университет, Краснодар, Россия
Данная статья является продолжением работ [1,2], которые были посвящены исследованию гидродинамики и переноса ионов соли в экспериментальной электрохимической ячейке с вращающейся дисковой катионообменной мембраной при допредельных токовых режимах, когда выполняется условие локальной электронейтральности. В данной работе приведена математическая модель переноса ионов соли в ячейке с вращающейся дисковой катионообменной мембраной при запредельных токовых режимах, с учетом электроконвекции. При этих условиях гидродинамика зависит от процесса переноса ионов соли и описывается системой уравнений Навье-Стокса в цилиндрической системе координат с учетом электрической (Кулоновской) силы
Ключевые слова: ОБЕССОЛИВАНИЕ, ВРАЩАЮЩАЯСЯ ДИСКОВАЯ МЕМБРАНА, РАВНОДОСТУПНАЯ ПОВЕРХНОСТЬ, ЭЛЕКТРОДИАЛИЗ, УРАВНЕНИЯ НАВЬЕ-СТОКСА, УРАВНЕНИЯ НЕРНСТА-ПЛАНКА-ПУАССОНА, ЭЛЕКТРОКОНВЕКЦИЯ, ЦИЛИНДРИЧЕСКАЯ СИСТЕМА КООРДИНАТ
UDC 517.958: 531.32, 544.6
THE TRANSFER OF SALT IONS IN AN ELECTROCHEMICAL CELL WITH ROTATING MEMBRANE DISK WITH ELECTRO CONVECTION. PART 1. MATHEMATICAL MODEL
Kovalenko Anna Vladimirovna Cand.Econ.Sci., associate professor
Urtenov Makhamet Ali Khuseevich Dr.Sci.Phys.-Math., professor
Kazakovtseva Ekaterina Vasilyevna postgraduate student
Kuban State University, Krasnodar, Russia
This article is a continuation of the works [1,2], which were devoted to the study of hydrodynamics and transport of salt ions in the experimental electrochemical cell with a rotating disk with a cation exchange membrane of exact current modes, when the condition of local electroneutrality. This article presents a mathematical model of transport of salt ions in a cell with a rotating disk with a cation exchange membrane exorbitant current regimes, taking into account electroconvection. Under these conditions, fluid dynamics depends on the ion transport process salt and described by the system of Navier-Stokes equations in cylindrical coordinate system with the electric forces
Keywords: DESALTING, ROTATING DISK MEMBRANE, FAIRNESS SURFACE, ELECTRO DIALYSIS, NAVIER-STOKES EQUATION, NERNST-PLANK-POISSON EQUATION, ELECTRO CONVECTION, CYLINDRICAL COORDINATE SYSTEM
Введение
В работе [1] была исследована гидродинамика экспериментальной электрохимической ячейки с вращающейся дисковой мембраной (ВДМ) с учетом ее реальных размеров при допредельных токовых режимах, когда
1 Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 13-08-00464 а.
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
2
выполняется условие локальной электронейтральности. Работа [2] была посвящена исследованию переноса ионов соли в экспериментальной электрохимической ячейки с вращающейся дисковой катионообменной мембраной, также при допредельных токовых режимах. В работе [3] был исследован перенос ионов соли в запредельных токовых режимах с использованием уравнения Пуассона вместо условия
электронейтральности, но без учета электроконвекции. При этих условиях гидродинамика не зависит от процесса переноса ионов соли и можно воспользоваться результатами работы [1]. Данная работа является продолжением работ [1,2,3]. В ней приведена математическая модель переноса ионов соли в закрытой ячейке при запредельных токовых режимах, с учетом электроконвекции, что приводит к существенному изменению гидродинамики.
1. Постановка задачи
Рассматривается задача о переносе ионов соли при вращении мембранного диска внутри вертикально стоящей цилиндрической ячейки вокруг центральной оси [1] при запредельных токовых режимах с учетом электроконвекции.
1.1 Область
При моделировании и численном решении используется осевая симметрия модели, поэтому описывается половина сечения цилиндрической области, где и определяются уравнения и граничные условия (рис.1).
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
3
Рисунок 1. Исследуемое сечение области и ее границы: 1 - глубина раствора, где выполняется условие электронейтральности, 2 - ось симметрии, 3 - катионообменная мембрана, 4- открытая граница
При интерпретации результатов нужно иметь в виду, что фигура изображенная на рис. 1 вращается вокруг оси симметрии 2. Граница 1 моделирует бесконечно удаленную от катионообменной мембраны часть пространства, где выполняется условие электронейтральности,
концентрация раствора постоянная (С0). Граница 1 считается также анодом и открытой границей (входом) для раствора. Граница 4 считается открытой границей (выходом) для раствора. Скорость течения раствора на входе и выходе определяется по ходу решения. В этой работе мы рассматриваем ячейку, вначале полностью заполненную раствором с концентрацией С0.
1.2 Уравнения
1) Моделирование переноса ионов соли
При указанных выше условиях моделирование переноса ионов соли может быть осуществлено с использованием системы уравнений Нернста-Планка для катионов и анионов, уравнения Пуассона для потенциала электрического поля. Векторная запись этой системы для бинарного
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
4
электролита при отсутствия химических реакций, в декартовой системе координат, имеет вид:
N, =Fz,D,C,E - D,V C, + CV, i = 1,2 RT (1)
dC r 2 1 = divNi, i = 1,2 (2)
єАФ = — F (z,C, + z 2C 2) (3)
1 = F (z1j1 + Z 212 ) (4)
где V - градиент, А - оператор Лапласа, V - скорость течения раствора, N1, N2, C1, C2 — потоки и концентрации катионов и анионов в растворе,
соответственно, z1, z2 - зарядовые числа катионов и анионов, 1 -плотность тока, D1, D2 - коэффициенты диффузии катионов и анионов,
соответственно, Ф - потенциал электрического поля, E = -VF -напряженность электрического поля, є - диэлектрическая проницаемость электролита, F - постоянная Фарадея, R - газовая постоянная, T -абсолютная температура, t - время, n - коэффициенты кинематической вязкости.
Уравнения Нернста-Планка (1) описывают поток растворенных компонентов, обусловленный миграцией в электрическом поле, диффузией и конвекцией; (2) - уравнение материального баланса в точке (малом элементе объема); (3) - уравнение Пуассона для потенциала электрического поля; (4) - плотность тока в растворе электролита,
обусловленая движением заряженных компонентов.
Замечание 1. Если подставить (1) в (2), то уравнения (2) запишутся в
виде
дСг = F dt ~ RT
ziDidiv(CiVF) + ЦА Сг - div(ClV), i = 1,2
Для решения задачи система уравнений Нернста-Планка и Пуассона записывается в цилиндрической системе координат.
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
5
2) Моделирование течения раствора
Для моделирования течения жидкости используются уравнения Навье-Стокса с объемной электрической силой f . В уравнении Навье-Стокса и обозначает скорость, р плотность, ц динамическая вязкость, и P давление:
Эи r
р-----hDU + р(и ■ V)u + VP = f
dt
(5)
V ■ и = 0, (6)
Влиянием гравитационной конвекции можно пренебречь, поскольку катионообменная мембрана расположена горизонтально и под ней образуется запирающий слой обессоленного раствора. Таким образом, при запредельных токовых режимах можно считать, что объемные силы в системе уравнений (5) являются только электрическими:
f = рЕ = -eDFE = eDFVF = eEdivE .
Модель включает все три скоростных компонента, однако для 3D осесимметричного потока циркулирующий поток является уже 2D, т.е. аргументов у неизвестных функций будут только две после перехода к цилиндрической системе координат (х, y, z) ® (r, z). Уравнения (5) в цилиндрических координатах примут вид:
Эи ( ^ -2 э.Л
р*+р
Эи v2 Эи
и------+ w—
Эг r Эz
Эр + — = л
Эг '
1 Э ( ЭиЛ
г
г
Эг
V
Эг
и
Э2 и
Эz7
+ bi
(7)
Эи
Рt+р
Эw
рэ7+р
ии
и
V
Эг
+ w-r Эz
\ Г1 э (эи v Э 2vl
=h У г Эг V Эг у r2 Эz2 _
+ bn,
^ Эw Эw^
и— + w— Эг Эz
Эр
+ — = л Эz '
і э
г Эг
^ Эw^ V Эг у
+
Э2
w
Эz2
+ b
3 5
где и - радиальная скорость, v - вращательная скорость, и w - осевая скорость (м\с).
2
r
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
6
Для использования системы уравнений необходимо выразить электрическую силу в цилиндрических координатах, т.е. найти
b = (b^ b2 , b3).
3) Расчет электрической силы в цилиндрических координатах
1) Вычислим силу в цилиндрической системе координат. Так как:
f(t,x, y, z) = fi( t, x, y, z)i + /2 (t, x, y, z)j + f3(t, x, y, z)k =
= b(r,j,z) = bi(r, j, z) re + b2 ( r, j, z) ij + Ьз( r, j, z)k,
bi(t,r,j,z) = fi( t,x,y ,z) cos j + /2 (t,x,y,z) sin j, где b2(t, r, j, z) = - fx(t,x, y, z) sin j + /2(t,x, y, z) cos j,
Ьз(t,r,j, z) = /3 (t,x,y,z),
b\( t,r, j,z) = fi(t ,x,y ,z) cos j + /2 (t,x,y,z) sin j, причем b2 (t,r,j,z) = -/l(t,x,y,z)sinj + /2 (t,x,y,z)cosj,
Ьз( t,r, j,z) = /з( t ,x,y ,z),
Эф Эф Эф
то /i(t,x,y,z) = eDF —, /2(t,x,y,z) = eDF —, /,(t,x,y,z) = eDF —,
Эx Эу Эz
x = r cosj, y = r sinj, z = z,
/ (x, y, z) = /i(x,y,z)i + /2 (x,y,z)j + /3 ( x, y, z )k =
= F( r, j, z) = Fi( r, j, z)?e + F2 (r, j, z)ij + F3 (r, j, z)k.
/ = eAFVF,
F(t, x, y, z) = F(t, r,j, z),
э2ф э2ф э2Ф 1 Э ґ эф 1 э2ф э2ф АФ = = (r ) ^ +
2
3x2 Эу2 Эz2 r Эг Эг r2 Э j2 Эz
„ _ Эф r Эф r Эф г Эф r 1 Эф r эф г
VF = — 1 +— 1 +—к = — г +---iq +—к .
Эx Эу Эz Эг r Э j Эz
2) Для задач с осевой симметрией, в том числе, для ячейки с ВМД функция Ф не зависит от j, поэтому:
1 э Эф э 2 ф
АФ =(r—) + (8)
r Эг Эг Эz2 W
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
7
„ _ эф r эф r эф ,
УФ = і +-т— J +к
Эх Эу
Следовательно,
Эх
ЭФ r _ . ЭФ г r + 0 • iq + —— к
Эг
Эх
(9)
r Л э / ЭФЧ э2Ф ЭФ r ,1 Э , ЭФЧ Э2ФЧЭФ r
J V Эг Эг Эх2) Эг г Эг Эг Эх2) Эх
т.к.
1 Э ЭФ Э2 Ф Э2 Ф Э2 Ф 1 ЭФ
----(г —) \------=--------\-----\------, то
г Эг Эг Эх2 Эг2 Эх2 г Эг
1 Э , Эф Э2Ф, ЭФ А2Ф Э2Ф, ЭФ 1 (ЭФЛ
^=£(aA+aAV+£-
г Эг 4 Эг 7 Эх27 Эг
/1 =£(~-(^г-) +
f2 = 0
г А Э , ЭФ Э2 Ф ЭФ
/з =£(- — () + ^тт)
г
Э2 ф э 2 ф эф
£(^ ^гтЬг-
+ £-
V Эг у 1 ЭФ ЭФ
г Этх Эт' Эх2 Эх Эг2 Эх27 Эх г Эг Эх Замечание 2. Из (8) следует, что уравнение Пуассона (3) в цилиндрической системе координат имеет вид:
2
1 Э , ЭФ Э2Ф F
_7Г( г^Г~) ^Г"Г =-----
г Эг Эг Эх £
1.3 Г раничные условия
Опишем граничные условия на каждой из границ.
Г раница 1 считается входом. На границе 1 концентрации катионов и анионов считаются постоянным: С10 = C0,i = 1,2. Для скорости ставится
(х1С1 + х 2 C 2 )
условие отсутствия нормального напряжения (Ум + (Ум )T )п = 0, давление при этом считается равным нулю. Граница 1 также считается эквипотенциальной поверхностью, причем Ф = 0 .
Граница 2 соответствует оси вращения, поэтому на этой границе используется условие симметрии.
Граница 3 соответствует вращающейся идеально селективной катионообменной мембране, поэтому она считается выходом для катионов, концентрация которых постоянна и равна емкости мембраны: С10 = Ckm.
Для анионов используется условие непроницаемости (отсутствия потока):
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
8
- n ■ N2 = 0. Поверхность катионообменной мембраны считается эквипотенциальной: Ф = dF. Для радиальной скорости используется условие: v = cor.
На открытой границе 4 для ионов ставятся условие выноса
выходом и для скорости ставится такое же граничное условие, как и для границы 1.
1.4 Начальные значения и свойства раствора.
Будем рассматривать водный раствор хлористого натрия. При проведении экспериментов возможны две различные методики:
1) Перед экспериментом ячейка полностью заполняется идеально перемешенным раствором хлористого натрия и в нее через границу №2 подается идеально перемешанный раствор. В качестве начального условия тогда берется постоянная концентрация С0 во всем объеме ячейки,
например С0 = 0.01 моль / м .
2) Ячейка заполнена чистой водой и в нее через границу №2 подается идеально перемешанный раствор. В этом случае в качестве начального условия берется нулевая концентрация.
В данной статье мы исследуем случай 1, когда ячейка полностью заполняется идеально перемешенным раствором хлористого натрия.
1.5. Метод численного решения
Для решения используется метод конечных элементов, реализованный в среде COMSOL Multiphysics 4.4 с неравномерной сеткой, количество элементов 6325.
2. Анализ численных результатов при постоянном начальном условии
http://ej.kubagro.ru/2014/09/pdf/80.pdf
конвективным потоком Ni =-u ■ Ct, i = 1,2. Для потенциала используется
условие
)T = 0. Граница 4 считается
Научный журнал КубГАУ, №103(09), 2014 года
9
Рассмотрим некоторые результаты моделирования переноса ионов соли в электрохимической ячейке с вращающимся мембранным диском при фиксированном скачке потенциала в 0.3В и изменение угловой скорости.
2.1 Анализ численных результатов при w = p/2
Из приведенных ниже на рис.2 линий тока раствора видно, что электроконвективный вихрь образуется около оси симметрии. В
пространстве этот вихрь имеет вид тора.
Скорость движения раствора в вихре очень мала и конвективный перенос сопоставим с диффузионным переносом, и как показывает рис.3, на распределение концентраций катионов и анионов практически не влияет.
а) б)
Рисунок 2. Линии тока раствора при угловой скорости w = p/2 рад/сек (15 оборотов в минуту) в момент времени t = 1000 с совместно с : а) - радиальной координатой скорости, б) - концентрацией катионов
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
10
Рисунок 3. Графики концентрация при угловой скорости w = p / 2 рад/сек (15 оборотов в минуту) в момент времени t = 1000 с: а) -катионов, б) - анионов
Как видно из рис. 3 возле мембраны наблюдается резко выраженный диффузионный слой, толщина которого практически не меняется вдоль радиуса мембраны. Интерес представляет время стабилизации процесса переноса.
Как следует из рис. 4 процесс достаточно быстро выходит на стационарный режим. Сопоставление рис.4 и рис.2 показывается, что время стабилизации примерно равна 200с.
2.2 Анализ численных результатов при w = 2p
При увеличении угловой скорости вращения, при том же скачке потенциала, течение раствора может качественно измениться.
Из рис.5 видно при увеличении угловой скорости электроконвективный вихрь исчез. Это связано с тем, что при увеличении угловой скорости радиальная скорость также увеличивается, а образование электроконвективных вихрей зависит от соотношения падения потенциала и линейной скорости (см. пороговую кривую в [4]).
Как следует из рис.6 при отсутствии электроконвекции ширина диффузионного слоя остается постоянной вдоль радиуса мембраны с большой точностью, что согласуется с теорией Левича [5].
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
11
О 10 О 10 О 10 О 10
О 10 О 10 О 10 О 10 ▼ 2.45x10
Рисунок 4. Изменение концентрации катионов и линий тока жидкости по времени при угловой скорости w = p/2 рад/сек (15 оборотов в минуту)
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
12
Surface: Velocity magnitude (m/s) A 6,28x10' Streamline: Velocity field х10’Л
_i_________________________________i_T 2,63x10'
0 10
Time=10Q0 s Surface: Concentration (mol/m3) Streamline: Velocity field
а) б)
Рисунок 5. Линии тока раствора при угловой скорости w = 2p
рад/сек (60 оборотов в минуту) в момент времени t = 1000 с совместно с : а) - радиальной координатой скорости, б) - концентрации катионов
Рисунок 6. Графики концентрация при угловой скорости скорости w = 2p рад/сек (60 оборотов в минуту) в момент времени t = 1000 с : а) -катионов, б) - анионов
Из сопоставления рис.5 и рис.7 следует, что с увеличение угловой скорости время установления (выхода на стационарный режим) уменьшается. Например, при угловой скорости w = 2p рад/сек, время
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
13
стабилизации примерно равно 10-12 с ( ср. с 200 с при угловой скорости w = p/2 рад). Это связано с отсутствием вихреобразования.
Рисунок 7. Изменение концентрации катионов и линий тока жидкости во времени при угловой скорости w = 2p рад/сек (60 оборотов в минуту)
Выводы
В работе рассчитана электрическая сила в цилиндрической системе координат и построена математическая модель переноса ионов соли в ячейке с вращающейся дисковой катионообменной мембраной при запредельных токовых режимах, с учетом электроконвекции в виде
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
14
краевой задачи для связанной системы уравнений Навье-Стокса и Нернста-Планка-Пуассона в цилиндрической системе координат.
Проведены численные исследования краевой задачи и показано, что образование электроконвективных вихрей начинается возле оси симметрии и зависит от соотношения угловой скорости вращения и скачка потенциала.
В дальнейшем планируется исследовать влияние электроконвективных вихрей на зависимость толщины диффузионного слоя от угловой скорости вращения и скачка потенциала и на равнодоступность поверхности мембраны при различных граничных условиях на ее поверхности с использованием вольтамперной характеристики.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 13-08-00464 а.
Литература
1. Коваленко А. В. Математическое моделирование и численное исследование гидродинамики в экспериментальной электрохимической ячейке с вращающимся мембранным диском / А.В. Коваленко, В.И. Заболоцкий, М.Х. Уртенов, Е.В. Казаковцева, М.В. Шарафан // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №10 (094). - Режим доступа: http://ei.kubagro.ru/2013/10/pdf/24.pdf
2. Коваленко А. В. Исследование переноса ионов соли в экспериментальной электрохимической ячейки с вращающейся дисковой мембраной / А. В. Коваленко,
В.И. Заболоцкий, М.Х. Уртенов, Е.В. Казаковцева, М.В. Шарафан // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №10 (094). - Режим доступа: http://ei .kubagro.ru/2013/10/pdf/25 .pdf
3. Kovalenko A. Mathematical modeling of transfer of salt ions in the electrochemical cell with rotating membrane disk in view of electroconvection / Urtenov M.,
http://ej.kubagro.ru/2014/09/pdf/80.pdf
Научный журнал КубГАУ, №103(09), 2014 года
15
Kovalenko A., Kazakovtseva E. // Ion transport in organic and inorganic membranes: proceeding international conference. 2014. p. 235-236
4. Коваленко А.В. Критериальные числа образования нестабильных
электроконвективных вихрей в канале обессоливания электродиализного аппарата / Коваленко А.В., Уртенов М.Х., Узденова А.М., Никоненко В.В. // Сорбционные и хроматографические процессы: научный журнал. Т. 14. № 2. 2014. с. 260-269
5. Левич В.Г. Физико-химическая гидродинамика. М.: Физматгиз, 1959, 700с.
References
1. Kovalenko A.V. Matematicheskoe modelirovanie i chislennoe issledovanie
gidrodinamiki v jeksperimental'noj jelektrohimicheskoj jachejke s vrashhajushhimsja membrannym diskom / A.V. Kovalenko, V.I. Zabolockij, M.H. Urtenov, E.V. Kazakovceva, M.V. Sharafan // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs]. - Krasnodar: KubGAU, 2013. - №10(094). - Rezhim dostupa:
http://ej.kubagro.ru/2013/10/pdf/24.pdf
2. Kovalenko A.V. Issledovanie perenosa ionov soli v jeksperimental'noj jelektrohimicheskoj jachejki s vrashhajushhejsja diskovoj membranoj / A.V. Kovalenko, V.I. Zabolockij, M.H. Urtenov, E.V. Kazakovceva, M.V. Sharafan // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs]. - Krasnodar: KubGAU, 2013. -№10(094). - Rezhim dostupa: http://ej.kubagro.ru/2013/10/pdf/24.pdf
3. Kovalenko A. Mathematical modeling of transfer of salt ions in the electrochemical cell with rotating membrane disk in view of electroconvection / Urtenov M., Kovalenko A., Kazakovtseva E. // Ion transport in organic and inorganic membranes: proceeding international conference. 2014. p. 235-236
4. Kovalenko A.V. Kriterial'nye chisla obrazovanija nestabil'nyh jelektrokonvektivnyh vihrej v kanale obessolivanija jelektrodializnogo apparata / Kovalenko A.V., Urtenov M.H., Uzdenova A.M., Nikonenko V.V. // Sorbcionnye i hromatograficheskie processy: nauchnyj zhurnal. T. 14. № 2. 2014. s. 260-269
5. Levich V.G. Fiziko-himicheskaja gidrodinamika. M.: Fizmatgiz, 1959, 700s.
http://ej.kubagro.ru/2014/09/pdf/80.pdf
