[На предыдущий раздел]

Субмодель (С1), описывающая явление ТП. Данная субмодель предназначена для описания температурного поля и ТП, ее включение необходимо при моделировании процесса, протекающего при неизотермических условиях. Подробное описание данной субмодели представлено в работах [12–15]. Она построена на основе транспортного интеграла, описывающего продольный перенос тепла:

Для упрощения математических преобразований использованы безразмерные переменные e и h, позволяющие использовать математические формулировки для профиля скорости, u(h), и температуры, t(e), в безразмерном виде.

Перенос в зоне А. Интегрирование и последующие преобразования выражения, описывающего транспорт в пределах зоны А (см. первое слагаемое в правой части выражения 19), осуществлялось в безразмерных координатах, e:

Промежуточные выражения X_A и Y_A представлены в форме единичных матриц, являющихся произведением соответствующих векторов.

Перенос в зоне В. Второе слагаемое в правой части выражения (19) отражает транспорт в пределах зоны В, лежащей между верхней границей температурного слоя и осевой линией канала, где наблюдается изменение скорости, в то время как температура остается постоянной, то есть равной температуре в ядре t₁:

Выражения (21) и (23) представляют собой количественную связь между продольным транспортом тепла и параметрами температурного и гидродинамического поля, а именно, осевой скоростью, U_MAX, и объемной температурой, t₁. Решенные относительно степени ТП (b₁), данные выражения позволяют прогнозировать изменение степени ТП по длине канала. Серия расчетных профилей степени ТП для плоскорамного мембранного модуля при различных значениях реологического индекса m показана на рис. 3. Данная субмодель дает связь между значениями температуры на поверхностях мембраны, величиной трансмембранного потока, параметрами мембраны и теплофизическими свойствами среды.

Субмодель (С2), описывающая явление КП. Данная субмодель, см. рис. 2, позволяет моделировать явление КП и прогнозировать профиль степени КП по длине канала. Подробное описание данной субмодели представлено в работах [16–18]. Она построена на основе транспортного интеграла, описывающего продольный перенос массы:

Для упрощения математических преобразований контрольное сечение потока подразделяется на две зоны: зона С, лежащая между Z = 0 и Z = d_C, и зона D, лежащая между Z = d_C и Z = Н, при этом транспортный интеграл (24) можно представить в виде суммы двух слагаемых:

Перенос в зоне С. Для дальнейшего моделирования переноса в пределах зоны С (см. первое слагаемое в правой части выражения 25) были использованы безразмерные координаты, q:

Перенос в зоне D. Второе слагаемое в правой части выражения (25) отражает транспорт в пределах зоны D, , лежащей между верхней границей диффузионного слоя и осевой линией канала, где наблюдается изменение скорости, в то время как концентрация остается постоянной, равной концентрации в ядре C₁:

Выражения (27) и (29) представляют собой количественную связь между продольным транспортом массы и параметрами концентрационного и гидродинамического поля. Выражения (27) и (29), решенные относительно степени КП (a), позволяют количественно прогнозировать изменение степени КП по длине канала. Серия расчетных профилей степени КП для плоскорамного мембранного модуля показана на рис. 4.

Данная субмодель содержит также блок уравнений, описывающих концентрационное поле в процессах микрофильтрации для более сложных случаев, когда коэффициент диффузии является функцией скорости деформации [11]. В основе субмодели лежит транспортное уравнение, описывающее перенос, где коэффициент диффузии, , и нормальная скорость, V(l), представляют функциональные зависимости от поперечной координаты, в безразмерных координатах, l, данное уравнение имеет вид:

где D(l) – коэффициент диффузии, зависящий от скорости деформации и, как следствие, изменяющийся от точки к точке, V(l) – поперечная скорость. Для описания поперечной составляющей скорости, V(l), может быть использована функция, предложенная Берманом в работе [6]. Она имеет вид:

Для описания зависимости коэффициента диффузии, , от скорости деформации, может быть рекомендовано выражение, предложенное Зидни и Колтоном [9, 10, 11], которое имеет вид: , где K_D = 0.03r².

Подстановка соответствующих выражений для скорости, V(l), и коэффициента диффузии, D(l), в выражение (30) дает выражение концентрационного профиля в процессах микрофильтрации, когда коэффициент диффузии является функцией скорости деформации, . Функция представлена в l-координатах, где l = 1 – z/H (см. рис. 1):

где K_D = 0.03r².

Представленная субмодель (32) позволяет анализировать влияние нормальной и тангенциальной гидродинамической составляющих, Re_Z и Re_X, геометрической полувысоты канала, Н, и радиуса суспендированных частиц, r, на конфигурацию концентрационного поля. Серия расчетных профилей степени КП для плоскорамного мембранного модуля показана на рис. 5.

Субмодель (С3), описывающая формирование и перемещение слоя геля вдоль фильтрующей поверхности. При моделировании перемещения слоя учитывается его неньютоновское поведение. Данная субмодель построена на основе следующего балансового уравнения:

где w_G(z) – профиль скорости перемещения слоя геля, C_G – концентрация геля в слое, принята постоянной и равной концентрации гелеобразования. Транспортные составляющие в выражении (33) проиллюстрированы на рис. 6. В результате преобразований данного выражения получен ряд соотношений в аналитическом виде для количественной оценки скорости формирования слоя и нестационарности процесса. Подробное описание данной субмодели представлено в работе [19].

Группа субмоделей (С4–С5) для описания транспорта через мембрану. Данная группа субмоделей основана на следующих физических предпосылках:

А. Мембрана представляет собой пористое тело, численно характеризуемое такими параметрами, как пористость, удельная внутренняя поверхность, извилистость капилляра, число пор на единице поверхности, а также кривой распределения по размерам f(R). Для аппроксимации плотности вероятности неоднородного распределения пор была использована функция Максвелла f(R) = A_MR² exp(B_MR²), где f(R) – плотность вероятности, пропорциональная числу пор, имеющих радиус в интервале между R и R + dR; А_M и В_M – параметры функции Максвелла, полученные на основе анализа гистограммы.

В. Жидкость представляет собой неньютоновское тело, при моделировании ее поведения учитываются зависимость вязкости, m_N, от скорости деформации (dV/dr), которая аппроксимирована следующей степенной функцией: m(r) = K₁(dV/dr)^m–1. Для описания течения в индивидуальной поре радиуса R использовано следующее выражение для профиля скорости: .

С. Трансмембранный перенос диспергированной фазы. В модели учитывается полидисперсность разделяемой суспензии. Для описания распределения частиц по размерам использована следующая функция [20]:

где f_SUSP(r) – плотность вероятности распределения диспергированной фазы в диапазоне от r до r + dr.

Данная субмодель базируется на “ситовом” механизме разделения, предполагающем размер проникающей частицы меньше размера капилляра. При проникновении через мембрану отсутствует физико-химическое взаимодействие между диспергированной фазой и материалом мембраны. Адсорбция и механическое задержание диспергированной фазы в объеме мембраны пренебрежимо мало. Представленная в работе субмодель позволяет прогнозировать фракционный состав фильтрата и концентрата на основании кривых распределения исходной суспензии и мембраны. Расчетные кривые распределения суспендированной фазы в фильтрате и концентрате приведены на рис. 7. Описание данной субмодели представлено в работе [21].

где dF_MEM = dxdy и dF_FLOW = dydz – контрольные сечения для продольного и поперечного потоков. В транспортные интегралы заложена математическая формулировка профиля скорости, учитывающего неньютоновское поведение фазы (см. выражение 6). В результате интегрирования и последующих преобразований получены следующие выражения для оценки соответствующих потоков механической энергии. Поток механической энергии через выходное контрольное сечение (см. левую часть выражения 35) представлен следующим образом:

где F_MEM = BDL.

Субмодель, описанная уравнениями (36–38) и решенная относительно давления в выходном контрольном сечении (Ре), позволяет количественно прогнозировать изменение давления и движущей силы по длине. Серия расчетных значений тангенциального изменения давления при различных реологических показателях степени (m) представлена на рис. 8.

[На следующий раздел] [На Содержание] [Условные обозначения] [Подстрочные индексы] [Безразмерные комплексы]

Copyright ©