Протекание через фрактальную пористую среду Текст научной статьи по специальности «Физика»

УДК 621.891

ПРОТЕКАНИЕ ЧЕРЕЗ ФРАКТАЛЬНУЮ ПОРИСТУЮ СРЕДУ

© 2011 В.П. Тихомиров, О.А. Горленко, М.А. Измеров Брянский государственный технический университет

Поступила в редакцию 10.11.2011

Контакт шероховатых поверхностей представлен в виде пористой среды, через которую происходит утечка уплотняемой среды. Рассмотрена методика определения параметров структуры пористой среды. Ключевые слова: фрактал, пористая среда, протекание.

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

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

Задача, как наилучшим образом описать просачивание уплотняемой жидкости или газа в пористой среде, представляет интерес в разных областях знания, в том числе и при проектировании металл-металлических герметичных устройств. Интерес также представляет существующее расхождение результатов, полученных экспериментально с теоретическими количественными оценками. Это связано, как показано ранее, с большим разбросом и недостаточной обоснованностью назначения и выбора значений факторов, входящих в уравнение Дарси [1]. Вот почему при проектировании конкретного герметизирующего устройства конструктор и технолог часто не имеют достоверных данных, чтобы на начальном этапе создания этого устройства обеспечить его работоспособное состояние. Так, коэффициент извилистости канала протекания, состоящего из цепочки связанных между собой пор, по литературным данным изменяется от 2 до 10. Не совсем обоснован выбор коэффициента формы канала протекания и т.д.

Тихомиров Виктор Петрович, доктор технических наук, профессор кафедры детали машин.

Горленко Олег Александрович, доктор технических наук, профессор кафедры управление качеством производствен-

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

Для фрактальной пористой среды (Пуазейль) расход может быть записан в виде

д <х — Ар сс ФАр.

и

Здесь г - радиус поры; l - длина поры; п - динамическая вязкость; Ар - перепад давления (напор).

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

Обозначим проницаемость и порог проницаемости соответственно через Ф и Фс (наименьшее значение проницаемости в бесконечном кластере пор). Вероятность протекания через пору - W(Ф).

Тогда порог протекания (концентрация пор, приводящая к формированию цепочки связанных между собой пор, через которую происходит утечка) определяется формулой р

С фс

Плотность распределения пор [2] по радиусам г

представим в виде W ^ — I)

-r

>2-Л max

ных и технических систем.

Измеров Михаил Александрович кандидат технических наук, доцент кафедры детали машин. E-mail: maik51@yandex.ru.

Здесь D - фрактальная размерность (1<D<2); rmax -радиус наибольшей поры.

Представление пористой среды в виде набора пор сферической формы является частным случаем набора вытянутых пор (в виде цилиндрических трубок). Пусть трубки имеют радиус r и длину Cr. Тогда объем трубки равен V=nr2Cr=Cnr3. Для пор сферической формы объем V=(4/3)nr3.

В этой модели [3] плотность распределения пор определяется зависимостью

- Г1'0 * С -л -гъ~в

^ п 'max

Вероятность нахождения трубки радиуса r равна

\Г ос г

или

-1-я

С -л - к

3-Я

ф.

Пористость определим из соотношения г

тах ____

т = =

-1-15

= I

(- ' 71 ' Vу

С -л -гъс1г.

Интегрируя это выражение, найдем

Л

т = 1-

г

V тах у

Связь между пористостью и фрактальной размерностью определяется выражением

£> = 3 —

1§<

-т

(1)

^т/'шах^

Обычно при определении пористости учитывают только связанные между собой поры и не учитывают изолированные поры, не участвующие в перемещении уплотняемой среды внутри пористой среды.

Используя фрактальные представления, используемые в работе [4], пористость может быть опреде-

1 и'Ю-Е"

лена из соотношения: т = 1 — о где Ь, г - коэффициенты; Б - фрактальная размерность профиля зазора (1<Б<2); Е - эвклидова размерность (Е=2).

Фрактальная размерность в этом случае будет

равна

М§6

(2)

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

Расположение пор (белые клетки на рис. 1) формируется с помощью распределения Бернулли. Остальные черные клетки указывают на отсутствие пор и невозможность протекания через них.

Проведенные многими исследователями численные эксперименты показали, что протекание наступает при концентрации пор (отношение числа пор к общему числу клеток, равному, например 70х70),

п

пор

оцениваемой величинои р = ——.

По

Здесь ппор - количество пор; по - общее число клеток квадратной решетки.

Анализ рис. 1 показывает, что при малой концентрации пор, не достигшей порогового значения,

протекание невозможно, а при превышении порогового значения наблюдается протекание. Пороговое значение концентрации:

• для 2-0 перколяционной модели /)6=0.5927:

• для 3-Б модели концентрация пор рс=0,3116.

б

Рис. 1. Заполнение квадратной решетки порами (белые квадраты):а) концентрация пор р=0,21б) концентрация пор р=0,73.

На рис. 2 представлен канал протекания в трёхмерном пористом слое, смоделированном с вероятностью заполнения 65% непроницаемыми кубами, при этом независимый сквозной канал один длиной 115 кубов при длине пористого слоя в 50 кубов. Однако эти теоретические модели не рассматривают анизотропию или негомогенность пористой структуры и пористость разного масштаба измерения.

Рис. 2. Канал протекания.

г

г

г

а

3-0 модель требует сравнения с реальными процессами протекания (фильтрации) в пористой среде (через межконтактный зазор инженерных металл-металлических поверхностей. Фрактальные представления о пористой среде позволяют дать более точные количественные оценки параметров, входящих в исходное уравнение А. Дарси.

Согласно абстрактным 2-0 и 3-0 моделям предполагается, что пористость численно равна критическому значению концентрации пор, при которой осуществляется переход от состояния герметичности к протеканию т =рс.

Так, для 2-0 модели при Ь=3; /=1 фрактальная размерность пористого слоя при т=0,5927 определяется по формуле

а

^ С-0,5927

м§3

■2 = 1,182.

Для 3-0 модели при Ь=3; /=1 фрактальная размерность пористого слоя при »7=0,3116 равна

в = х^-т\Е =

/1

^<-0,3116^

м§3

■3 = 2,662.

Сравним полученную оценку фрактальной размерности с расчетом по формуле (1). При отношении гтах/гт,„=3 и т=0,3116 найдем 0=2.66.

При малых значениях соотношения (гтп / гтах)3'° , где гт,„=тс - критическая величина радиуса поры, ниже которой просачивание невозможно, справед-

ливо выражение

= 1 - т.

г

\ шах /

Тогда отношение гс / гт

тах имеет вид

' о'

■т

о

Отношение гс / гтах зависит от пористости и фрактальной размерности (критическое значение радиуса изменяется в пределах [3] от гс=гтах/5 до гс=гтах/2). Отметим, что минимальный радиус поры, через которую еще возможно просачивание, является для рассматриваемого конкретного случая постоянной величиной. Минимальный (критический) радиус поры, в свою очередь, зависит от многих факторов: вязкости, числа Рейнольдса, наличия граничного слоя прилипших к стенкам поры молекул уплотняемой среды и др.

В ряде случаев удается установить связь между пористостью и фрактальной размерностью, которую можно записать в первом приближении в виде как т = 3-Д 2<£><3.

Зависимость отношения гт1„/гтах от фрактальной размерности Б представлена на рис.3.

Полученные данные о фрактальной размерности можно считать предварительными. Окончательную оценку фрактальной размерности и пористости (межконтактного зазора) дадим с помощью компьютерного моделирования.

2 2.2 2.4 2.6 2.3 В

Рис. 3. Зависимость отношения гтп / гтах от фрактальной размерности Б.

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

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

б

Рис. 4. Поверхности 3-Б: а - модель, б - реальная поверхность.

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

з-и

г

шт

Г

с

Г

а

друга. Для упрощения модели материалы двух сопряжённых поверхностей считаются одинаковыми. Сближение поверхностей автоматически продолжается до момента, когда протекание жидкости в направлении сверху вниз становится невозможным -наступает полная герметичность.

Рис. 5. Моделирование контакта поверхностей.

Картина распределения пор, пятен касания и каналов представлена на рис. 6. Эта картина меняется при увеличении сближения. В этом случае площадь, занятая металлическим контактом, растет.

Таким образом, существующие модели протекания (фильтрационная - Дарси, капиллярная - Гаге-на-Пуазейля, щелевая) содержат для согласования с экспериментальными данными ряд коэффициентов, которые в ряде случаев либо трудно определить, либо принимаемые в расчетах значения недостаточно обоснованы.

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

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

Фрактальная модель позволяет рассматривать пористость на разных масштабных уровнях. Так, для стыка мультифрактальных поверхностей пористость

равна сумме, определяющей раздельно особенности пористого слоя на разных уровнях

т=Шц + тг,

Здесь ггц =

3-/Х

I

Рис. 6. Картина зазора между шероховатыми поверхностями приразличном сближении (а - 7 мкм, б - 20 мкм, с - 27 мкм; светлые пятна -поры, темные - площадки металлического контакта, серые - капилляры).

3- Д

I

Дополнительным условием подобного разделения является следующее утверждение: радиус САмой наименьшей «структурной» поры всегда больше радиуса самой большой «текстурной» поры.

Для определения пористости необходимо определить фрактальные размерности этих особенностей пористой среды.

Итак, фрактальная модель в определенной степени объединяет фильтрационную и капиллярную

г

2-й

г

г

$ тах

г

$ тт

а

б

в

г

2- Д

I =

г

г

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

СПИСОК ЛИТЕРАТУРЫ

1. Суслов, А.Г. Технологическое обеспечение параметров состояния поверхностного слоя деталей.-М.: Машиностроение, 1987.-208 с.

2. Hunt, A. Percolation Theory for Flow in Porous Media/A. Hunt, R. Ewing.-Springer, Berlin Heidelberg, 2009.-306 p.

3. Hunt, A.G. Effects of Bubbles on the Hydraulic Conductivity of Porous Materials. Theoretical Results/ A. G. Hunt, M. Manga// Transport in Porous Media, 2003.-V. 52.-P. 51-65.

4. Sukop, M.C. Percolation Thresholds in 2-Dimensional Prefractal Models of Porous Media/ M.C. Sukop, G-J. Van Dijk,E. Perfect, W.K.P. Loon// Transport in Porous Media: Kluwer Academic Publishers. Printed in the Netherlands, 2002.-48. -P. 187-208.

PERCOLATION THROUGH FRACTAL POROUS MEDIA

© 2011 V.P. Tikhomirov, O.A. Gorlenko, M.A. Izmerov Bryansk State Technical University

Contact of rough surfaces is presented in the form of the porous environment through which there is a leak of the condensed environment. The technique of definition of parameters of structure of the porous environment is considered. Key words: a fractal, the porous media, percolation.

Tikhomirov Victor Petrovich, Dr.Sci.Tech., the professor of chair of a detail of cars.

Gorlenko Oleg Aleksandrovich, Dr.Sci. Tech., the professor of chair quality management ofproizvodstven th and technical systems. Izmerov Michael Aleksandrovich Cand.Tech.Sci., the senior lecturer of chair of a detail of cars. An E-mail: maik51@yandex.ru.