СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ Российский патент 2006 года по МПК G01N15/08 

Описание патента на изобретение RU2275617C2

Изобретение относится к измерительной технике и может быть использовано при оценке качества пористых материалов, например керамики, металлокерамики.

Известны способы определения коэффициента диффузии (А.С. № 652473 кл G 01 N 15/00, 1979 г.), коэффициента фильтрации (А.С. № 1056001 А, кл G 01 N 15/08, 1983 г.), способ измерения газопроницаемости материалов (А.С. №750346 кл G 01 N 15/08, 1980 г.).

Недостатком известных способов является то, что они позволяют измерить только один из параметров пористости материалов.

Наиболее близким по технической сущности к предлагаемому решению является способ определения параметров пористости материалов (Патент РФ № 2186365 кл G 01 N 15/08, 2002 г.).

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

где Vn - номинальный объем измерительной камеры;

М - молярная масса газа;

F - площадь образца, разделяющего камеры;

R0 - газовая постоянная;

Т - температура газа в измерительной камере;

Pc - давление в измерительной камере при равенстве предыдущего и последующего приращений давления;

Pa - атмосферное давление;

υ - коэффициент кинематической вязкости, м2/с;

tc - время в измерительной камере при равенстве предыдущего и последующего приращений давления;

Рn - давление во время tn, где tn>tc;

l - толщина исследуемого образца.

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

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

Это достигается тем, что в способе, заключающемся в вакуумировании герметично соединенной с материалом измерительной камеры, регистрации температуры и измерении величины давления в измерительной камере за фиксированное время, регистрируют изменение объема газа в измерительной камере с момента поступления газа из окружающей среды через площадь контролируемого материала, ограниченную измерительной камерой, определяют давление и время в измерительной камере, соответствующие переходу от кнудсеновского к ламинарному режиму течения газа - Рс, tc, а коэффициенты фильтрации (Кф), проницаемости (К0), диффузии (D) растворимости (константа Генри - Кг) и пористости (П) определяют из соотношений:

- коэффициент растворимости:

- коэффициент фильтрации:

- коэффициент пористости:

- коэффициент газопроницаемости

- коэффициент диффузии:

где или

где Vu - объем газа, протекающего в единицу времени tu при давлении Рu в измерительную камеру, м3;

F - площадь контролируемого материала, ограниченная измерительной камерой, м2;

R0 - газовая постоянная, Дж/(кг·К);

μ - коэффициент динамической вязкости, Па·с;

Рc, tc - давление (Па) и время (с) в измерительной камере, соответствующие переходу от кнудсеновского к ламинарному режиму течения газа;

Рu, tu - текущее значение давления (Па) и времени (с), tu>tc;

Pa - атмосферное давление, Па;

ν - коэффициент кинематической вязкости, м2/с;

l - толщина исследуемого материала, м;

d - диаметр материала, ограниченный измерительной камерой, м.

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

Основными уравнениями, характеризующими процесс проникновения воздуха через материал являются:

- уравнение Дарси:

- уравнение Фика для стационарного потока:

- уравнение Фика для нестационарного потока:

- уравнение Лейбензона для изотермической фильтрации:

- добавочное уравнение:

где ρ - плотность газа, кг/м3; q - поток газа в единицу времени, кг/(с·м2); Кф - коэффициент фильтрации, с; l - линейный размер в направлении фильтрации, м; ΔР - перепад давления в материале (потери напора), Па; D - коэффициент диффузии, м2/с; ▿Р - модуль градиента давления в исследуемом материале, Па/м, ▿ρ - модуль градиента концентрации в исследуемом материале, кг/м4; К0 - коэффициент проницаемости, м2; П - коэффициент пористости, %; μ - коэффициент динамической вязкости, Па·с; Ра - атмосферное давление, Па; F - площадь материала, м2; R0 - газовая постоянная, Дж/(кг·К); Vu - объем, м3; Рu - давление, Па; Т - температура, К; t - время, с.

Для двумерного случая (движение газа через пористый материал длиной l и диаметром d) уравнения (2)-(4) примут вид:

- уравнение Фика для стационарного потока:

- уравнение Фика для нестационарного потока:

- уравнение Лейбензона:

Решение уравнения (7) ρ(x,y,t) должно удовлетворять следующим граничным и начальным условиям:

0<x<l, 0<y<d, t>0,

Для решения системы уравнений (6)-(8) рассмотрим уравнение (7) в следующем виде:

где материал начинаем рассматривать с противоположного торца:

Решением уравнения (10), удовлетворяющим граничным и начальным условиям (11), является функция:

где

В качестве начального условия возьмем функцию:

Тогда

Полученное значение аij подставляем в решение (12), после чего имеем

Учитывая граничные условия (11), т.e. то, что , получаем решение уравнения (7), удовлетворяющее условиям (9).

При d=∞, из выражения (17), как частный случай, вытекает известное решение уравнения , а именно:

Введем обозначения:

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

Следовательно, вместо выражения (17) получаем

Граничные условия (9) имеют вид:

Используя уравнения (20) и (6) при х1=l, y1=d, определим поток q:

Для решения системы уравнений (6)-(8) используем добавочное уравнение (5), что позволяет получить следующее выражение:

или

При интегрировании от 0 до t, уравнение (23) примет вид:

Если d=∞, имеем одномерный случай:

Если l=∞, получаем результат при двумерном случае:

Определим величину двойной суммы, раскладывая в числовой ряд:

тогда при l=∞ и t>tc

где

Полагая, что в уравнении (8) Р2=Ф(х,у,t), вместо уравнения (8) получаем:

где

Решение уравнения (29) удовлетворяет следующим граничным и начальным условиям:

Фy(x,0,t)=0(x,d,t)=0,

Ф(х,у,0)=ƒ(х,у), 0<х<1, 0<y<d.

Вместо уравнения (29) рассмотрим следующее уравнение:

где материал начинаем рассматривать с противоположного торца, т.е. из вакуумированной полости:

Решением (31), удовлетворяющим условиям (32) является функция:

где

По аналогии с (14) в качестве начального условия возьмем функцию:

Тогда

Следовательно:

Учитывая граничные условия (32), что , получаем:

Введем обозначения:

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

Учитывая, что Ф(х,у,t)=Р2(х,y,t), вместо выражения (37) с учетом (38), имеем:

Граничные условия (30) выглядят так:

Для определения давления на входе вакуумированной полости согласно уравнению Лейбензона предварительно найдем массу газа (воздуха), прошедшую в вакуумированную полость за время t проведения эксперимента.

Тогда давление на входе вакуумированной полости равно:

Если d=∞, получаем одномерный случай:

Пологая, что , уравнение (43) примет вид:

При выходе на линейный участок в уравнении (44) заменяем Рu на Рuc, т.е.

Если l=∞, то

Так как

то при t>tc

где

Рассмотрим связи между параметрами. Из сравнения левых частей уравнений (28) и (48) получаем:

где

Сравнивая уравнение Дарси с уравнением Фика для стационарного потока, получаем:

ΔP≈Pa-Pu≈Pa.

Отсюда,

Из уравнения (28) находим коэффициент диффузии:

где Pu, tu - текущее значение давления и времени.

Тогда коэффициент проницаемости:

Найденное значение К0 подставляем в (49), получаем коэффициент пористости:

повышая тем самым функциональные возможности способа.

Известно, что

тогда коэффициент фильтрации:

Из равенства (52) имеем коэффициент растворимости газа:

Таким образом, для двумерной модели истечения газа по порам материала, получаем следующие зависимости:

- коэффициент растворимости:

- коэффициент фильтрации:

- коэффициент пористости:

- коэффициент газопроницаемости

- коэффициент диффузии:

где или .

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

На чертеже представлена схема устройства для реализации данного способа.

Устройство содержит измерительный стакан 1, соединенный с контролируемым материалом 2, который при помощи скобы (условно не показано) прижимается к материалу 2 усилием Q, образуя измерительную камеру 3 устройства. Кольцо 4 из вакуумной резины служит для герметизации системы.

Измерительная камера 3 посредством клапана 5 соединена с вакуумным насосом (условно не показан), который позволяет создавать разряжение в измерительной камере 3.

Для измерения изменения давления в камере 3 перед клапаном 5 установлены датчик давления 6 и датчик времени 7, выходы которых подсоединены к ЭВМ 8. Датчик давления 6 соединен с усилителем сигнала 9, выход которого подключен к запоминающему устройству 10. Усилитель сигнала 9 и запоминающее устройство 10 соединены со сравнивающим устройством 11, выход которого соединен с преобразователем сигнала 12, второй выход которого подключен к ЭВМ 8.

Способ осуществляется следующим образом.

После открытия клапана 5 и создания вакуума при помощи вакуумного насоса в измерительной камере 3 клапан 5 закрывается. Так как с внешней плоскости материала атмосферное давление, а в измерительной камере 3 вакуум, будет наблюдаться фильтрационный поток газа (Фф) через материал. Давление в камере 3 будет повышаться. Изменение давления с течением времени будет отслеживаться датчиком давления 6 и датчиком времени 7, данные с которых передаются на вход ЭВМ 8. Концентрация газа в камере 3 при повышении давления будет изменяться. Усилитель сигнала 9 при изменении давления газа будет усиливать сигнал об изменении приращения давления в камере 3, а запоминающее устройство 10 - фиксировать каждое последующее приращение давления соответствующее определенному промежутку времени. Сравнивающее устройство 11 осуществляет сравнение последующего приращения давления с предыдущим и подает информацию о новом значении давления в измерительной камере 3 на преобразователь сигнала 12. Преобразователь сигнала преобразовывает приращение давления газа в камере 3 в приращение объема газа и подает данные на ЭВМ 8, которая строит графическую зависимость изменения давления в измерительной камере 3 с течением времени, определяет численные значения давления Рc, Рu, времени tc и tu (фиксированные промежутки времени), где tu>tc, Рuс и объема газа Vu.

Данные о площади F материала, ограниченной измерительной камерой, атмосферном давлении Рa, газовой постоянной R0, коэффициенте кинематической вязкости ν, коэффициенте динамической вязкости μ, толщине l и диаметре d материала введены в ЭВМ 8 заранее. Поэтому ЭВМ 8 осуществляет определение коэффициента фильтрации из формулы:

коэффициента проницаемости из формулы:

коэффициента диффузии из формулы:

коэффициента растворимости (константа Генри) из формулы:

а коэффициента пористости из формулы:

Похожие патенты RU2275617C2

название год авторы номер документа
СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ 2006
  • Житников Юрий Захарович
  • Иванов Алексей Николаевич
  • Матросова Юлия Николаевна
  • Матросов Александр Евгеньевич
RU2305828C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОНИЦАЕМОСТИ И ПОРИСТОСТИ МАТЕРИАЛОВ 2003
  • Житников Ю.З.
  • Иванов А.Н.
  • Матросова Ю.Н.
  • Матросов А.Е.
RU2247358C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ 1999
  • Можегов Н.А.
  • Житников Ю.З.
  • Матросова Ю.Н.
RU2186365C2
СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ 2003
  • Житников Ю.З.
  • Иванов А.Н.
  • Матросова Ю.Н.
  • Матросов А.Е.
RU2235308C1
СПОСОБ ОПРЕДЕЛЕНИЯ МИНИМАЛЬНЫХ ДИАМЕТРОВ ПОР ОБРАЗЦА 1999
  • Можегов Н.А.
  • Житников Ю.З.
  • Матросова Ю.Н.
RU2186364C2
Способ определения режима движения потокагАзА 1979
  • Бондаренко Анатолий Георгиевич
  • Костина Наталья Владимировна
  • Британ Игорь Михайлович
  • Голубев Илья Федорович
SU834521A1
СПОСОБ ОПРЕДЕЛЕНИЯ АКТИВНОЙ ПОРИСТОСТИ МАТЕРИАЛОВ 2003
  • Житников Ю.З.
  • Матросова Ю.Н.
  • Матросов А.Е.
RU2244915C1
СПОСОБ ИЗМЕРЕНИЯ ПОРИСТОСТИ ЧАСТИЦ СЫПУЧИХ МАТЕРИАЛОВ 2012
  • Медведева Анна Валерьевна
  • Мордасов Денис Михайлович
  • Мордасов Михаил Михайлович
RU2527656C2
СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА МАССООТДАЧИ ПОРИСТЫХ ПРОНИЦАЕМЫХ МАТЕРИАЛОВ 2014
  • Светлов Юрий Валентинович
  • Никифоров Юрий Борисович
RU2566726C1
СПОСОБ ЗАПОЛНЕНИЯ ТЕПЛОВЫХ ТРУБ 2010
  • Серяков Аркадий Владимирович
RU2431101C1

Реферат патента 2006 года СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ

Использование: изобретение относится к измерительной технике и может быть использовано при оценке качества различных пористых материалов. Сущность способа заключается в том, что измерительный стакан 1 соединяют с контролируемым материалом 2, образуя измерительную камеру 3. Измерительную камеру 3 вакуумируют. Регистрируют датчиком давления 6 и датчиком времени 7 изменение давления в камере 3 с течением времени. Усилителем сигнала 9 усиливают сигнал об изменении приращения давления в камере 3 и фиксируют каждое последующее приращение давления запоминающим устройством 10. Осуществляют сравнение последующего приращения давления с предыдущим и преобразователем сигнала 12 преобразуют приращение давления газа в камере 3 в приращение объема газа, а коэффициенты фильтрации (Кф), проницаемости (К0), диффузии (D), растворимости (константа Генри - Кг) и пористости (П) определяют из соотношений: - коэффициент растворимости: коэффициент фильтрации: - коэффициент пористости: - коэффициент газопроницаемости: - коэффициент диффузии: где или ; где Vu - объем газа, протекающего в единицу времени tu при давлении Рu в измерительную камеру, м3; F - площадь контролируемого материала, ограниченная измерительной камерой, м2; R0 - газовая постоянная, Дж/(кг·К); μ - коэффициент динамической вязкости, Па·с; Рс, tc - давление (Па) и время (с) в измерительной камере, соответствующие переходу от кнудсеновского к ламинарному режиму течения газа; Рu, tu - текущее значение давления (Па) и времени (с), tu>tc; Pa - атмосферное давление, Па; ν - коэффициент кинематической вязкости, м2/с; l - толщина исследуемого материала, м; d - диаметр материала, ограниченный измерительной камерой, м. Технический результат изобретения заключается в повышении точности. 1 ил.

Формула изобретения RU 2 275 617 C2

Способ определения параметров пористости материалов, при котором герметично соединенную с материалом измерительную камеру вакуумируют, регистрируют температуру и измеряют величину давления в измерительной камере за фиксированное время, отличающийся тем, что регистрируют изменение объема газа в измерительной камере с момента поступления газа из окружающей среды через площадь контролируемого материала, ограниченную измерительной камерой, определяют давление и время в измерительной камере, соответствующие переходу от кнудсеновского к ламинарному режиму течения газа Рс, tc, а в качестве параметров пористости материала определяют коэффициент фильтрации (Кф), коэффициент проницаемости (К0), коэффициент диффузии (D), коэффициент растворимости (константа Генри - Кг) и коэффициент пористости (П) из соотношений

коэффициент растворимости

коэффициент фильтрации

коэффициент пористости

коэффициент газопроницаемости

коэффициент диффузии

где или

где Vu - объем газа, протекающего в единицу времени tu при давлении Рu в измерительную камеру, м3;

F - площадь контролируемого материала, ограниченная измерительной камерой, м2;

μ - коэффициент динамической вязкости, Па·с;

Рc, tc - давление (Па) и время (с) в измерительной камере, соответствующие переходу от кнудсеновского к ламинарному режиму течения газа;

Рu, tu - текущее значение давления (Па) и времени (с), tu>tc;

Pa - атмосферное давление, Па;

ν - коэффициент кинематической вязкости, м2/с;

l - толщина исследуемого материала, м;

d - диаметр материала, ограниченный измерительной камерой, м.

Документы, цитированные в отчете о поиске Патент 2006 года RU2275617C2

СПОСОБ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПОРИСТОСТИ МАТЕРИАЛОВ 1999
  • Можегов Н.А.
  • Житников Ю.З.
  • Матросова Ю.Н.
RU2186365C2
СПОСОБ ОПРЕДЕЛЕНИЯ МИНИМАЛЬНЫХ ДИАМЕТРОВ ПОР ОБРАЗЦА 1999
  • Можегов Н.А.
  • Житников Ю.З.
  • Матросова Ю.Н.
RU2186364C2
СПОСОБ ОПРЕДЕЛЕНИЯ АКТИВНОЙ ПОРИСТОСТИ МАТЕРИАЛОВ 1999
  • Можегов Н.А.
  • Житников Ю.З.
  • Матросова Ю.Н.
RU2180743C2
Печь для непрерывного получения сернистого натрия 1921
  • Настюков А.М.
  • Настюков К.И.
SU1A1

RU 2 275 617 C2

Авторы

Житников Юрий Захарович

Иванов Алексей Николаевич

Матросова Юлия Николаевна

Матросов Александр Евгеньевич

Даты

2006-04-27Публикация

2004-06-28Подача