УДК 663.479.1 (075)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОБЕССОЛИВАНИЯ ВОДЫ МЕТОДОМ ЭЛЕКТРОДИАЛИЗА
Т.А. Краснова, А.Г. Семёнов
Кемеровский технологический институт пищевой промышленности бульвар Строителей, 47, Кемерово, Россия, 650056
Описано применение электродиализа для очистки природной воды от избыточного количества неорганических солей. Предложены принципы построения математической модели электро-диализного аппарата, получены уравнения модели, учитывающей динамику движения воды и изменение характеристик мембран по длине аппарата и с течением времени.
Ключевые слова: электродиализ, неорганические соли, динамика движения воды.
Проблема опреснения соленых и обессоливания пресных вод в настоящее время является актуальной. Одним из перспективных способов очистки природной воды от содержащихся в ней неорганических солей является электродиализ. В аппарате для электродиализа очищаемая вода проходит через систему параллельных каналов, образованных чередующимися катионо- и анионообменными мембранами. Находясь в воде, мембраны приобретают определенный заряд и получают способность обмениваться с ней ионами, имеющими заряд противоположного знака («противоионами»). Под действием приложенного электрического поля ионы растворенных солей мигрируют. При этом каждая мембрана пропускает через себя в основном противоионы. В результате в одних (рассольных) камерах происходит концентрирование растворенных веществ, а в других (ди-люатных) раствор освобождается от них.
Сочетание процессов конвективного, диффузионного и электромиграцион-ного переноса ионов в потоках, протекающих через камеры, приводит к вознино-вению концентрационной поляризации — перераспределению содержания растворенных солей и появлению у поверхностей мембран концентрационных пограничных слоев, называемых также поляризационными или диффузионными. Концентрация примесей в диффузионном слое резко отличается от их концентрации в основном потоке через камеру. Оптимизация конструкции электродиализатора и определение оптимальных режимных параметров процесса электродиализа требуют построения адекватных математических моделей, учитывающих как перенос ионов в мембранах, так и конвективный перенос в камерах.
Наиболее известная гидродинамическая модель электродиализа [1] основана на следующих основных допущениях: движение воды в камерах является ламинарным; мембраны обладают идеальной селективностью; рассматривается разбавленный раствор бинарного электролита, диссоциация воды не учитывается; выполняется условие локальной электронейтральности; рассматривается стационарный процесс. Полученные уравнения решены приближенным методом [1]. Получены численные решения для той же модели [2]. Аналогичная модель, построенная для нестационарного процесса, использована для моделирования уста-
новок порционной обработки воды [3]. Модель обобщена за счет введения в качестве постоянных параметров чисел переноса ионов в мембранах и отказа от условия бинарности электролита [4]. Полученные уравнения также решены приближенным методом.
В данной статье предлагается гидродинамическая модель процесса элект-родиализного обессоливания природной воды и растворов минеральных солей. Приняты следующие допущения: движение воды в камерах является ламинарным и профиль скорости описывается уравнением Хагена-Пуазейля; выполняется условие локальной электронейтральности; диссоциация воды не учитывается. Перенос ионов соли в сплошной среде описывается уравнениями
§1 = ёУ - Д.гД.У\|Т- , (1)
Л = (2)
где: § — плотность потока ионов .-сорта (т.е., катионов или анионов электролита), — плотность парциального тока, с. — концентрация ионов, — заряд иона (с учетом знака), Д — коэффициент диффузии иона в среде, Г — константа Фарадея, у = Г ф / ЯТ —
безразмерный электрический потенциал (ф — потенциал, Я — газовая постоянная, Т — температура).
После подстановки (1, 2) в уравнения сохранения числа ионов
+V -ус = ду2ё
д* 1 1 1
и сохранения заряда У - ] = 0 (] — суммарный ток), после ряда преобразований получаем:
дс- + V -Уё = ДУ 2ё, (3)
ді
У- (сУу) + 0У2с = 0, (4)
где с = 2сес = 2аеа — приведенная концентрация электролита (индексы с, а обозначают характеристики соответственно катионов и анионов), 2с, 2а — абсолютные величины зарядов ионов.
В уравнениях введены параметры
О О (2 + 2 ) О — О
О = с а^ с-------а±, 0 =-----£----. (5)
2 О + 2 О 2 О + 2 О
с с а а с с а а
Поскольку длина аппарата значительно превышает толщину камер и мембран, используется допущение тонкого слоя: в камерах пренебрегаем конвективным переносом вдоль координатной оси 0у, направленной по нормали к поверхности мембран и продольным (вдоль оси 0х) переносом электрического заряда. В мембранах учитывается перенос и изменение всех характеристик только вдоль оси 0у. В этом случае уравнение (3) принимает вид
дс 3G
- + -
ді 4mhd
(Л у2 1 дс „д 2с
1—
. d2 У
Граничные условия к этому уравнению учитывают, что мембраны пропускают как противоионы, так и некоторое количество коионов (имеющих заряд того же знака, что и мембрана). Этот факт учитывается путем введения чисел переноса ионов в мембране і{ = / ] . Граничные условия имеют вид
| = ( — Іар ) <7>
на поверхностях анионообменной мембраны и
| = ( — ^)' - V <8)
на поверхностях катионообменной мембраны. В этих условиях введены параметры (миграционные числа переноса ионов в растворе)
і = 2аОа і = 2сОс
аР 2 О + 2 Б ’ сР 2 О + 2 Б'
с с а а с с а а
Если в растворе числа переноса можно считать постоянными, то в мембранах они зависят от концентраций прилегающих растворов. По мере нарастания разности примембранных концентраций рассола и дилюата как с течением времени, так и вдоль поверхности мембран числа переноса противоионов в мембранах уменьшаются, что снижает эффективность процесса. Для оценки чисел переноса необходимо проинтегрировать уравнения Нернста—Планка, которые получаются из уравнений (1, 2), рассмотренных в проекции на 0у при отбрасывании конвективных членов. При этом следует учесть, что из закона сохранения числа ионов и уравнения Пуассона при выполнении условия локальной электронейтральности в мембране следует, что парциальные токи ионов и напряженность электрического поля в мембране вдоль оси 0у оказываются неизменными. В результате получаем:
1
і„ =-
1 +
V 2а У
V с а у
ЄХр(в) — ( "с/с ’с ) ^ {1 — ЄХр(—Р)
'У
ехр( В) — 1
(9)
,(<5 "а/ё 'а )- еХР(-Р)^
*с = 1 - *а •
Здесь одним штрихом отмечены концентрации соответствующих ионов у поверхности мембраны с той стороны, с которой в нее входит вектор напряженности электрического поля, двумя штрихами — поверхностные концентрации с противоположной стороны. Параметры:
в = -2 л —, Р = -гЛ —, (10)
с т 1 ’ а т » ’
лу лу
где Лт — толщина мембраны.
Поверхностные концентрации ионов в мембране ёт1 связаны с концентрациями тех же ионов в растворе у поверхности мембраны ёр1 соотношениями Доннана
^ = ехр (-2 ). (11)
с рі
Кроме того, в мембране должно выполняться условие электронейтральности
2сёс + 0 = 2аёа , (12)
где о — плотность фиксированного заряда мембраны (с учетом знака).
Совместное решение уравнений (11; 12) позволяет определить доннанов-ский скачок потенциала Аул и поверхностные концентрации ионов в мембране.
Система замыкается интегральным соотношением, выражающим условие постоянства напряжения, приложенного к ячейке, по длине аппарата. Оно получается путем интегрирования (4) в растворе и мембранах вдоль направления 0у с учетом (11):
RT
F
j =----
------0
V zc У
—+ 0
(-dd) Vza У cd (dd)
ln-
■Фо
—, (13)
R + RT (f dp dl + fdd dy}
m F2 (zaDa + zcDc )VJ-dp~cp J-dd~cdj
где p, d —характеристики рассольных и дилюатных камер соответственно; d — полу-толщина камеры; Rm — общее омическое сопротивление пары мембран; ф0 — напряжение в расчете на ячейку.
Соотношения (6—13) образуют замкнутую систему уравнений, позволяющую моделировать процессы электродиализа и оценить влияние геометрических и режимных параметров.
ЛИТЕРАТУРА
[1] Sonin A.A., Probstein R.F. Hydrodynamic theory of desalination by electrodialysis // Desalination. — 1968. — V. 168. — P. 293—329.
[2] Заболоцкий В.И. и др. Конвективно-диффузионная модель процесса электродиализно-го обессоливания: Распределение концентраций и плотности тока // Электрохимия. — 1985. — Т. 21. — № 3. — С. 296—302.
[3] Краснова Т.А., Семенов А.Г. Математическое моделирование электродиализных установок с замкнутыми контурами // Теоретические основы химической технологии. — 1994. — Т. 28 — № 2. — С. 158—163.
[4] Кузьминых В.А., Григорчук О.В., Шапошник В.А. Гидродинамическая модель электродиализа с ионообменными мембранами разной селективности // Электрохимия. — 1994. — Т. 30. — № 9. — С. 1101—1108.
MATHEMATICAL MODEL OF WATER DESALTING BY ELECTRODIALYSIS
T.A. Krasnova, A.G. Semenov
Kemerovo Technological Institute of Food Industry Boulevard Builders, 47, Kemerovo, Russia, 650056
Application of an electrodialysis for the natural water treating from superfluous quantity of inorganic salts is described. Principles of the modeling of electrodialysis apparatus are offered. The equations are obtained for the model, that take into account fluid dynamics and changes of membrane characteristics with the length of apparatus and in time.
Key words: water treatment, electrodialysis, mathematical model.