Введение. Фильтрование жидкостей с тонкодисперсными примесями встречается во многих отраслях народного хозяйства [1, 2, 3]. В качестве примера можно назвать проблему восстановления технических масел, подготовку технической воды промышленности, очистку масел растительного происхождения и сгущения фруктовых соков в пищевой промышленности. Одним из методов фильтрования тонкодисперсных суспензий является использование вспомогательных фильтрующих материалов. Вспомогательный фильтрующий материал может быть добавлен в исходную суспензию с целью повышения проницаемости слоя образуемого осадка. Тогда рассматривается процесс фильтрования трехфазной гетерогенной среды. Вспомогательный фильтрующий материал может быть использован в качестве фильтрующей перегородки. В этом случае он намывается отдельно в виде нижнего слоя в начале цикла фильтрования. Тонкодисперсная суспензия фильтруется через нижний слой, образуя слой тонкодисперсных частиц. При этом возникает задача расчета фильтрации среды через двойной слой осадка, когда толщины верхнего слоя является переменной величиной. В настоящее время изучение процесса фильтрования неньютоновских сред с тонкодисперсными частицами на многослойных проницаемых подложках остается актуальной задачей [4, 5, 6].Цель исследований - математическое моделирование и численный расчет процесса фильтрования тонкодисперсных неньютоновских суспензий, подчиняющихся степенной модели реологического состояния среды, в намывных фильтрах периодического действия.Условия, материалы и методы. Многие аномально вязкие среды хорошо описываются моделью степенной жидкости. При этом обобщенное уравнение фильтрации имеет вид [7, 8, 9] (1) В данном уравнении параметр n указывает на степень нелинейности среды. При значении n=1 получаем классическое уравнение фильтрации Дарси. Через коэффициенты Ki – обозначены проницаемости слоев пористого осадка. Величина Ψ называется эквивалентной вязкостью. Она зависит не только от реологических параметров жидкости, но и от пористости и проницаемости соответствующих слоев осадка. Для вычисления эквивалентной вязкости известны разные формулы, которые справедливы при определенных диапазонах изменения скорости и конкретных моделей пористой среды [9]. Существуют разные конструкции намывных фильтров периодического действия – пластинчатые, трубчатые [10, 11, 12]. Возможны разные режимы их работы: фильтрование при постоянной разности давлений, фильтрование при постоянной скорости, а также фильтрование при переменных разности давлений и скорости. В этих аппаратах предварительно намывается нижний фильтрующий слой с достаточно высокой проницаемостью, через которого фильтруется тонкодисперсная суспензия. По мере увеличения толщины второго слоя возрастает сопротивление, оказываемое потоку жидкой фазы. В режиме фильтрования при постоянной разности давлений уменьшается производительность, в режиме фильтрования при постоянной скорости увеличивается разность давлений. После достижения их предельных значений цикл работы аппарата прекращается. Далее выполняется промывка фильтрующего элемента, и далее цикл повторяется заново. Результаты и обсуждение. Пусть на поверхности трубчатого фильтра образован слой вспомогательного материала толщины δ = const. В процессе фильтрования тонкодисперсной суспензии происходит накопление второго слоя осадка h(t). Уравнения фильтрации (1) для соответствующих слоев в цилиндрической системе координат имеют вид (2)  (3) Запишем уравнения сохранения массы (4)  (5) Для решения уравнений фильтрации и сохранения массы необходимо задавать граничные условия. В рабочей зоне фильтра создается давление P, которое определяется с учетом режима работы аппарата. В математической модели это давление задается на поверхности верхнего слоя осадка (6) На границе раздела двух слоев задаются условия сшивания для давления и скорости. Поскольку пористости разных слоев осадка различаются, на границе раздела произойдет скачкообразное изменение давления и скорости фильтрации. Эти условия могут быть записаны в виде [13] (7) где εδ, εh- пористости нижнего и верхнего слоев осадка.На поверхности фильтрующего трубчатого элемента задается выходное давление (8) Решение системы уравнений (2), (3), (4) и (5) при заданных граничных условиях (6), (7), (8) не представляют труда. После интегрирования уравнения (4) имеем  (9) После подстановки зависимости (9) можно проинтегрировать уравнения фильтрации (2)    Постоянного интегрирования  С2 определим из условия (8). Тогда последнее соотношение примет вид   Отсюда найдем  (10) Аналогично проинтегрируем уравнения (3) и (5). Если интеграл уравнения (5) представим в виде   то решения уравнения (3) ,с учетом граничного условия (6), запишется так  (11) Условие сшивания скоростей  εhVh=εδVδ на границе r=R+δ дает зависимость εhC3=εδC1. При этом постоянная интегрирования C3 выражается через константу C1 в виде:   Тогда зависимость (11) можно представить следующим образом (12) Полученные соотношения для давлений (10) и (12) должны выполнять граничное условие (7), с помощью которого определяется постоянное интегрирования С1. После определения С1, в результате несложных преобразований, для расчета скоростей фильтрации окончательно получим  (13)   (14) Уравнения (13), (14) позволяют вычислять скорость фильтрации сплошной фазы суспензии в аппаратах трубчатого типа. Полученные уравнения могут быть преобразованы для расчета фильтровальных аппаратов пластинчатого типа. Для этого необходимо выполнить предельный переход R→∞. Предварительно проведем необходимые преобразования.Несложно показать, что Выполним предельный переход R→∞, используя замечательный предел  после замены переменных 1/R=x  или  1/(R+δ)=x. С учетом того, что при условии R→∞ имеет место предел r/R→1, получим      Тогда, уравнение (13) для фильтровальных аппаратов пластинчатого типа примет вид (15) а уравнение(14) остается в силе.Если фильтруемая среда подчиняется линейному реологическому закону состояния, то уравнения (13) также необходимо преобразовывать, используя предельный переход  n→1. С помощью замечательного предела    можно показать, что имеет место равенства        Тогда уравнение (13) преобразуется к виду (16) Построим уравнение для расчета толщины слоя осадка. При заданной концентрации дисперсных частиц άT, их приток на элементарную площадку ∆S=(R+δ+h)∆l составляет   Продифференцируем данное соотношение с учетом условия (1-άT)V=εhVh на границе осадок-суспензия (17) С другой стороны за время ∆t объем слоя осадка на единичной площадке ∆S изменится на величину ∆q=(1—εh)∆S∆hРазделим обе части данного равенства на ∆t и выполним предельный переход ∆t→0: (18) Приравнивая правых частей соотношений (17) и (18), для расчета изменения толщины верхнего слоя осадка получим дифференциальное уравнение (19) Построенное уравнение решается с начальным условием h=0 при t=0.Численные расчеты по построенным математическим моделям были проведены в безразмерных переменных   где V* - характерная скорость, которая вычисляется через число Рейнольдса. В частности, после перехода к безразмерным переменным, уравнение (13) примет следующий вид (знаки тильда ~ опущены),   где  При этом вид уравнения (19) не изменится.Некоторые результаты расчетов, для режима фильтрования с постоянной скоростью, приведены на рисунках 1-3. В режиме фильтрования с постоянной скоростью, по мере увеличения толщины осадка, происходит возрастание давления в аппарате.На процесс фильтрации сильно влияют реологические параметры суспензии Ψi и n. Эквивалентная вязкость среды обратно пропорциональна числу Рейнольдса Re. На рис. 1 показано влияние параметра Re на интенсивность возрастания давления. Как видим, более интенсивный рост давления наблюдается при малых значениях числа Re. Следовательно, при увеличении эквивалентной вязкости среды, предельное значение давления в аппарате достигает раньше, что приведет к сокращению времени фильтровального цикла. Влияние степени нелинейности среды n показано на рис. 2. При увеличении этого параметра давление в аппарате растет быстрее, следовательно, сокращается время достижения критического давления.При фильтровании суспензий с повышенной концентрации нарастание слоя осадка происходит более интенсивно. В результате этого быстрее растет давление в аппарате и сокращается длительность фильтровального цикла (рис.3).     Выводы. Построена математическая модель фильтрования тонкодисперсных суспензий, подчиняющихся степенной модели реологического состояния среды, в намывных фильтрах периодического типа. Для режима фильтрования с постоянной скоростью приведены результаты численных расчетов. Установлены закономерности влияния эквивалентной вязкости, степени нелинейности среды, а также исходной концентрации дисперсной фазы, на длительность фильтровального цикла.         .