Теоретические основы химической технологии, 2023, T. 57, № 5, стр. 612-616

Математическое описание электро-массобменных процессов на основе уравнения переноса субстанции

Л. В. Равичев a*, С. И. Ильина a, В. Я. Логинов a, В. И. Быков a, А. А. Титов a

a Российский химико-технологический университет им. Д.И. Менделеева
Москва, Россия

* E-mail: ravichevlv@list.ru

Поступила в редакцию 02.07.2023
После доработки 11.07.2023
Принята к публикации 23.07.2023

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

Аннотация

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

Ключевые слова: электро-массобменные процессы, электродиализ, математическое описание, уравнение переноса субстанции

ВВЕДЕНИЕ

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

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

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

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

• механические;

• гидромеханические;

• теплообменные;

• массообменные.

Необходимо отметить:

1. Большинство рассматриваемых процессов являются сопряженными так, например, на протекание теплообменных и массообменных процессов существенно влияет гидродинамика потоков.

2. В приведенном выше перечне не выделены процессы, протекающие под действием разности электрических потенциалов, однако, такие процессы (или сопряженные с ними) также нашли свое применение в химической промышленности.

МЕТОДЫ РАСЧЕТА ХИМИКО-ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ

В настоящее время существует множество различных методов расчетов протекания химико-технологических процессов.

Развитие компьютерной техники позволяет быстро решать сложные задачи моделирования. Однако в основе построения математических моделей используются те же принципы, что и несколько десятилетий назад. К таким методам относится и использование теории подобия. На основании базовых уравнений выделяются критерии подобия, определяющие протекание процесса и далее составляются уравнения, устанавливающие соотношения между рассматриваемыми критериями. В данном методе выделяются следующие этапы:

1. Нахождение базового уравнения для выделения безразмерных критериев.

2. Установление зависимости между критериями.

3. Определение граничных условий адекватности критериальных уравнений.

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

Анализ вышеупомянутых уравнений переноса показал, что все они могут быть получены из уравнения переноса субстанции

(1)
$\frac{{\partial \varphi }}{{\partial \tau }} = - {\text{div}}\bar {q} + \gamma .$

Плотность притока может быть обусловлена, например, влиянием поля или образованием продуктов химической реакции [1, 2].

Дальнейшие выводы сводятся к определению вида переноса субстанции (материи), где плотность потока $\bar {q}~$ складывается из переноса на молекулярном и конвективном уровне, а также при выводе уравнения Навье–Стокса учитывается поле сил тяжести (учет этот происходит за счет составляющей внешних сил – плотность притока γ). Необходимо отметить, что при выводе уравнения Фурье–Кирхгофа, поток тепла за счет теплового излучения не учитывается, а при выводе уравнений конвективной диффузии величина $\gamma $ принимается равной 0. В то же время основные постулаты гласят о корпускулярно-волновом характере материи, из чего следует, что в целом перенос субстанции должен иметь обязательно три составляющие, а именно: перенос на молекулярном уровне, перенос на конвективном уровне и влияние поля.

Итак, если схематично изобразить виды переноса субстанции, то получится следующая картина (рис. 1):

Рис. 1.

Составляющие переноса субстанции.

Теперь рассмотрим известные уравнения переноса субстанции.

Перенос импульса – уравнение Навье–Стокса отражает все составляющие. Поле обусловлено силой тяжести – гравитационное поле.

Перенос тепла – уравнение Фурье–Кирхгоффа объединяет молекулярный и конвективный перенос, но математически не затрагивает описание поля. Тепловое излучение обусловлено электромагнитными волнами.

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

Однако субстанция (или материя) помимо массы и энергии имеет и другие виды – одним из которых является электрический заряд.

Так как заряд является видом материи, то, с учетом того, что потенциалом переноса является плотность заряда, для его переноса уравнение (1) может быть преобразовано в

(2)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} = - {\text{div}}{{\vec {J}}_{{\text{q}}}} + {{\gamma }_{{\text{q}}}}~.$

Перенос заряда – на молекулярном уровне перенос заряда описывается законом Ома и имеет вид, аналогичный законам Ньютона, Фурье, Фика. Влияние электромагнитного поля описывается уравнениями Максвелла. Математическое описание конвективного переноса заряда в отечественной и зарубежной литературе представлено фрагментарно. Однако очевидно, что явление конвективного переноса заряда необходимо учитывать при описании процессов электромассо- и электротеплообмена.

Исходя из утверждения о единстве подходов к описанию переноса всех видов субстанции, можно получить уравнение конвективной электропроводности по аналогии с уравнениями конвективной теплопроводности и конвективной диффузии.

Итак, выразим вышеизложенное в виде соответствующих уравнений:

(3)
${{\vec {J}}_{{\text{q}}}} = {{\vec {J}}_{{{{{\text{q}}}_{{\text{м}}}}}}} + {{\vec {J}}_{{{{{\text{q}}}_{{\text{к}}}}}}}.$

На молекулярном уровне плотность потока заряда описывается уравнением электропроводности:

(4)
${{\vec {J}}_{{{{{\text{q}}}_{{\text{м}}}}}}} = - \sigma \frac{{\partial U}}{{\partial n}}.$

На конвективном уровне плотность заряда описывается следующим уравнением:

(5)
${{\vec {J}}_{{{{{\text{q}}}_{{\text{к}}}}}}} = \vec {\omega } \cdot {{\rho }_{v}}$

${{\vec {J}}_{{{{{\text{q}}}_{{\text{к}}}}}}}$ имеет размерность $\left[ {\frac{{\text{м}}}{{\text{с}}}\frac{{{\text{Кл}}}}{{{{{\text{м}}}^{{\text{3}}}}}}} \right]$ = $\left[ {\frac{{{\text{Кл}}}}{{{\text{с}}\,\,{{{\text{м}}}^{{\text{2}}}}}}} \right]$.

Принимая ${{\gamma }_{{\text{q}}}}$ = 0 уравнение (2) можно преобразовать к виду

(6)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} = - {\text{div}}\left( { - \sigma \frac{{\partial U}}{{\partial n}} + \bar {\omega }{{\rho }_{v}}} \right).$

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

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Нами было проведено преобразование уравнения (6) по аналогии с уравнением Фурье–Кирхгофа.

С учетом того, что отношение заряда к напряжению – есть емкость системы:

(7)
$C = \frac{q}{U}.$

Напряжение можно выразить через плотность заряда:

(8)
$U = \frac{{{{\rho }_{v}}V}}{C} = \frac{{{{\rho }_{v}}m}}{{C\rho }} = \frac{{{{\rho }_{v}}}}{{{{C}_{{\text{m}}}}\rho }}.$

Далее, подставляя полученное уравнение (8) в уравнение электропроводности (4), получаем уравнение плотности потока заряда на молекулярном уровне

(9)
${{\vec {j}}_{{{{{\text{q}}}_{{\text{м}}}}}}} = - \sigma {\text{grad}}\left( {{{{{\rho }_{v}}} \mathord{\left/ {\vphantom {{{{\rho }_{v}}} {{{C}_{{\text{m}}}}\rho }}} \right. \kern-0em} {{{C}_{{\text{m}}}}\rho }}} \right) = ~ - \frac{\sigma }{{{{C}_{{\text{m}}}}\rho }}{\text{grad}}{{\rho }_{v}}.$

Тогда запишем

(10)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} + {\text{div}}\bar {\omega }{{\rho }_{v}} = b~{\text{div\;grad}}{{\rho }_{v}},$
а с учетом математических преобразований

(11)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} + {{\rho }_{v}}{\text{div}}\bar {\omega } + \bar {\omega }{\text{grad}}{{\rho }_{v}} = b{{\nabla }^{2}}{{\rho }_{v}}.$

Полученное уравнение учитывает перенос заряда в различных средах (как в неразрывных, так и при наличии кавитации).

Если поток неразрывен (${\text{div}}\bar {\omega } = 0$), уравнение (11) имеет вид

(12)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} + {{\omega }_{x}}\frac{{\partial {{\rho }_{v}}}}{{\partial x}} + {{\omega }_{y}}\frac{{\partial {{\rho }_{v}}}}{{\partial y}} + {{\omega }_{z}}\frac{{\partial {{\rho }_{v}}}}{{\partial z}} = \frac{\sigma }{{{{C}_{{\text{m}}}}\rho }}{{\nabla }^{2}}{{\rho }_{v}},$
или
(12а)
$\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} = b{{\nabla }^{2}}{{\rho }_{v}},$
где, обозначив $\frac{\sigma }{{{{C}_{{\text{m}}}}\rho }}$ = b, можно определить b как “коэффициент зарядопроводности”, названный по аналогии с коэффициентом температуропроводности, $\left[ {\frac{{{{{\text{См}}} \mathord{\left/ {\vphantom {{{\text{См}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}}}{{\frac{{{\text{Кл}}}}{{{\text{В\;кг}}}}\frac{{{\text{кг}}}}{{{\text{\;}}{{{\text{м}}}^{3}}}}}}} \right] = \left[ {\frac{{\frac{{{\text{Кл}}}}{с}}}{{\frac{{{\text{Кл}}}}{{{{{\text{м}}}^{2}}}}}}} \right] = \left[ {\frac{{{{{\text{м}}}^{2}}}}{{\text{с}}}} \right]$.

Полученное уравнение (12) является уравнением конвективной электропроводности, которое выражает в общем виде распределение плотности заряда в движущемся потоке. Нетрудно заметить, что вид уравнения (12) аналогичен уравнению конвективной теплопроводности Фурье–Кирхгофа и уравнению конвективной диффузии, что соответствует приведенной схеме (рис. 1) и еще раз подтверждает единство явлений переноса субстанций различного вида. Это подтверждается также тем, что коэффициент зарядопроводности b имеет ту же размерность, что и коэффициент температуропроводности, коэффициенты диффузии и кинематической вязкости.

Как и другие дифференциальные уравнения (Навье–Стокса, Фурье–Кирхгофа, конвективной диффузии) уравнение (12) может быть решено численными методами или использовано как основное для получения математических моделей. На практике для описания процессов переноса заряда чаще используют напряжение. Так как плотность тока ${{\rho }_{{\text{v}}}}$ зависит от напряжения U, то можно сделать вывод, что слагаемое $\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }}$, включающее частную производную напряжения по времени $\frac{{\partial U}}{{\partial \tau }}$, так как $\frac{{\partial {{\rho }_{v}}}}{{\partial \tau }} = {{C}_{{\text{m}}}}\rho \frac{{\partial U}}{{\partial \tau }}$, характеризует параметры переменного тока.

Из уравнения (12) можно получить числа Пекле и Фурье. Электрическое число Пекле можно выразить как

(13)
$\frac{{b{{\rho }_{v}}}}{{{{l}^{2}}}}\sim \frac{{\omega {{\rho }_{v}}}}{l} \Rightarrow \frac{{\omega {{\rho }_{v}}{{l}^{2}}}}{{lb{{\rho }_{v}}}} = \frac{{\omega l}}{b} = {\text{P}}{{{\text{e}}}_{{\text{q}}}}.$

Электрическое число Фурье соответственно будет иметь вид

(14)
$\frac{{b{{\rho }_{v}}}}{{{{l}^{2}}}}\sim \frac{{{{\rho }_{v}}}}{\tau } \Rightarrow \frac{{b{{\rho }_{v}}\tau }}{{{{l}^{2}}{{\rho }_{v}}}} = \frac{{b\tau }}{{{{l}^{2}}}} = {\text{F}}{{{\text{o}}}_{{\text{q}}}}.$

На практике широко используемым критерием является число Прандтля, которое рассчитывается как отношение числа Пекле к числу Рейнольдса

(15)
${\text{P}}{{{\text{r}}}_{{\text{q}}}} = {{{\text{P}}{{{\text{e}}}_{{\text{q}}}}} \mathord{\left/ {\vphantom {{{\text{P}}{{{\text{e}}}_{{\text{q}}}}} {{\text{Re}}}}} \right. \kern-0em} {{\text{Re}}}} = \frac{{\omega l}}{b}\frac{\mu }{{\omega l\rho }} = \frac{\mu }{{b\rho }} = \frac{{\mu {{C}_{{\text{m}}}}\rho }}{{\sigma \rho }} = \frac{{\mu {{C}_{{\text{m}}}}}}{\sigma }.$

Анализируя размерности (табл. 1), можно заметить сопоставимость размерностей переменных, входящих в критерий Прандтля.

Таблица 1.  

Сопоставление размерностей классического и электрического критерия Прандтля

Классический критерий Прандтля $\Pr = \frac{{{{C}_{{\text{p}}}}\mu }}{\lambda }$ $ \leftrightarrow $ Электрический критерий Прандтля ${{\Pr }_{{\text{q}}}} = \frac{{{{C}_{{\text{m}}}}\mu }}{\sigma }$
обозначение размерность   обозначение размерность
${{C}_{{\text{p}}}}$ $\frac{{{\text{Дж}}}}{{{\text{кг\;К}}}} \to \frac{{{\text{Вт\;с}}}}{{{\text{кг\;К}}}}$ $ \leftrightarrow $ ${{C}_{{\text{m}}}}$ $\frac{{\text{Ф}}}{{{\text{кг}}}} \to \frac{{{\text{Кл}}}}{{{\text{В\;\;кг}}}}$
$\lambda $ $\frac{{{\text{Вт}}}}{{{\text{м\;К}}}}$ $ \leftrightarrow $ $\sigma $ $\frac{1}{{{\text{Ом\;м}}}} \to \frac{{{\text{Кл}}}}{{{\text{В\;с\;м}}}}$
$\mu $ ${\text{Па\;с}} \to \frac{{{\text{кг}}}}{{{\text{м\;с}}}}$ $ \leftrightarrow $ $\mu $ ${\text{Па\;с}} \to \frac{{{\text{кг}}}}{{{\text{м\;с}}}}$

Таким образом, выражения для расчета электрических чисел, полученные из уравнения конвективной электропроводности, совпадают с выражениями, полученными методом аналогий и применяемыми при составлении математических описаний [35]. Кроме того, в этом случае более ясен физический смысл получаемых критериев. Электрическое число Прантдля (Prq) характеризует проводящие свойства материала, а электрическое число Фурье характеризует влияние переменного тока на процесс переноса заряда.

Полученные критерии электрического подобия могут быть использованы как для описания процессов переноса заряда, так и для сопряженных процессов, таких как электромассообменных или электротеплообменных. Например, для электромассообменных процессов можно записать

(16)
$\begin{gathered} {\text{N}}{{{\text{u}}}_{{\text{D}}}} = {\text{Sh}} = \\ = {\text{f}}\left( {{\text{Fo}},{\text{F}}{{{\text{o}}}_{{\text{q}}}},{\text{F}}{{{\text{o}}}_{{\text{D}}}},{\text{Re}},{\text{Pr}},{\text{P}}{{{\text{r}}}_{{\text{q}}}},{\text{P}}{{{\text{r}}}_{{\text{D}}}},{i \mathord{\left/ {\vphantom {i {{{i}_{{{\text{lim}}}}}}}} \right. \kern-0em} {{{i}_{{{\text{lim}}}}}}},{l \mathord{\left/ {\vphantom {l d}} \right. \kern-0em} d}, \ldots } \right). \\ \end{gathered} $

Также, в зависимости от режима течения жидкости, в уравнении (16) может быть включено число Фруда (Fr).

ЗАКЛЮЧЕНИЕ

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

ОБОЗНАЧЕНИЯ

$C$ емкость, Ф
${{C}_{{\text{m}}}}$ удельная емкость, Ф/кг
${{C}_{{\text{p}}}}$ удельная теплоемкость, Дж/кг К
$\overrightarrow {{{J}_{{\text{q}}}}} $ поток заряда, Кл/м3 с
$U$ напряжение, В
$b$ коэффициент зарядопроводность, м/с
$d$ диаметр, м
$i$ плотность тока, А/м2
$l$ линейный размер, м
m масса, кг
$q$ заряд, Кл
$\bar {q}~$ плотность потока
$\gamma ~$ плотность притока
${{\gamma }_{{\text{q}}}}$ плотность притока заряда, Кл/м3 с
$\mu $ вязкость, Па с
$\rho $ плотность, кг/м3
${{\rho }_{v}}$ плотность заряда, Кл/м3
σ удельная электропроводность, См/м
$\tau $ время, с
$\varphi $ потенциал переноса
$\vec {\omega }$ скорость, м/с
${\text{Fo}}$ число Фурье
${\text{Nu}}$ число Нуссельта
${\text{Pr}}$ число Прандтля
${\text{Re}}$ число Рейнольдса
${\text{Sh}}$ число Шервуда

ИНДЕКСЫ

D диффузионное число
K конвективный уровень
M молекулярный уровень
q электрическое число
lim предельное значение

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

  1. Жакин А.И. Электрогидродинамика // Успехи физических наук, 2012. Т. 182. № 5. С. 495–520. https://doi.org/10.3367/UFNr.0182.201205b.0495

  2. Саранин В.А. Электрогидродинамика: Равновесия, зарядка и конвекция жидких масс в электрических полях: дис. Док. Тех. Наук: 01.02.05 // Саранин Владимир Александрович. – Глазов, 1999.

  3. Ильин В.А., Пономарева Л.А. Электроконвекция слабопроводящей жидкости в высокочастотном электрическом поле // Вестник пермского университета. 2013. № 3(25). С. 28–36.

  4. Ильин В.А., Куршина Е.В. Исследование модели электроконвекции идеального диэлектрика в конденсаторе с твердыми границами // Вестник пермского университета. 2009. № 1(27). С. 3–6.

  5. Гросу Ф.П., Болога М.К., Болога А.М. Особенности теплообмена в условиях электрической конвекции // Электронная обработка материалов. 2010. № 4. С. 41–55.

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