способ определения активной пористости материалов

Формула изобретения

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





;

где



Г(i) = [z(i)y(i)]^1/2;

^*, ^**, ^***-

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

А_n - амплитуда колебаний давления в замкнутом пространстве после образца;

P - перепад давлений на образце;

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

соответственно отношение амплитуд в замкнутом пространстве после образца и до образца при ^*, ^**, ^***;

L = L_м, где - коэффициент, учитывающий извилистость пор и потери давления на входе и выходе из образца;

R - эквивалентный радиус пор;

П - значение активной пористости;

₀ - плотность газа в установившемся состоянии;

i - комплексное число;

К - коэффициент адиабаты;

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

J₀, J₂ - функции Бесселя;

ch(i), sh(i)- эллиптические синус и косинус;

V - газовый объем пространства после образца;

а₀ - скорость звука.

Описание изобретения к патенту

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

Известен способ определения пористости материалов, заключающийся в помещении исследуемого тела в измерительную емкость, соединении измерительной емкости с калиброванной емкостью избыточного давления, перепуске газа из калиброванной емкости в измерительную и измерении давлений в емкостях до и после перепуска (А.С. 1368720 кл. G 01 N 15/08, 1988 г.).

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

Наиболее близким к предлагаемому изобретению является способ определения пористости материалов, включающий ламинарную фильтрацию газа через заключенный в обечайку из газонепроницаемого материала образец, создании гармонических колебаний газа определенной частоты и амплитуды в замкнутом пространстве перед образцом и измерения колебаний давления газа в замкнутом пространстве после образца (А.С. 1679287 кл. G 01 N 15/02, 1991 г.).

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

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

Поставленная цель достигается тем, что в известном способе, включающем ламинарную фильтрацию газа через заключенный в обечайку из газонепроницаемого материала образец из замкнутого пространства перед образцом в замкнутое пространство после него, путем создания в замкнутом пространстве перед образцом гармонических колебаний давления газа определенной частоты и амплитуды, и измерения амплитуды колебаний газа в замкнутом пространстве после образца, гармонические колебания давления газа в замкнутом пространстве перед образцом осуществляют с изменяющейся частотой, измеряют изменение амплитуды давления газа в замкнутом пространстве перед образцом и после образца, вычисляют разность амплитуд и определяют три значения частоты из условия равенства разности амплитуд, определяя активную пористость материалов в результате совместного решения на ЭВМ трех уравнений:







где







^*, ^**, ^*** - соответственно три значения частоты, при которых обеспечивается одинаковое значение



где А_а - амплитуда колебаний давления в замкнутом пространстве перед образцом;

A_n - после образца;

P - перепад давлений на образце;

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

соответственно отношение амплитуд в замкнутом пространстве после образца и до образца при ^*, ^**, ^***.

₀ - плотность газа в установившемся состоянии;

i - комплексное число;

R - эквивалентный радиус пор;

Р₀ - давление в установившемся состоянии;

К - коэффициент адиабаты;

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

J₀, J₂ - функции Бесселя;

ch(i), sh(i) - эллиптические синусы и косинусы;

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

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

V - газовый объем пространства после образца;

П - значение активной пористости материала;

а₀ - скорость звука.

В системе цилиндрических координат х, r, , при допущении осевой симметрии емкостей 4 и 3 и обечайки 1 (при условии, что образец 2 имеет поры в виде капилляров, длина которых намного превышает диаметр) математическая модель перетекания газа по этим порам может быть представлена следующими уравнениями:

- уравнение Навье-Стокса или уравнение, вытекающее из закона сохранения количества движения:



- уравнение неразрывности или уравнение, вытекающее из закона сохранения массы



- уравнение адиабаты:

PV^K=const (3)

Линеаризируем уравнение (3). Будем иметь



где - плотность газа;

Р - давление газа;

u - осевая составляющая скорости течения газа;

радиальная составляющая скорости течения газа;

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

r - текущий радиус капилляра;

К - показатель адиабаты.

Решив совместно уравнения 2 и 4, исключим плотность и, умножив полученное выражение на r, выразим

Будем иметь:

(5)

Проинтегрировав последнее от 0 до R по r при условии, что при r=0 и r=R, а также получим:



Аналогично проинтегрируем уравнение 1. Используя условие u/r/_r=0 = 0 и предварительно умножив на r, будем иметь:



где коэффициент кинематической вязкости.

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



где среднее значение осевой составляющей скорости в поперечном сечении капилляра образца 2.

В результате уравнения 6 и 7 преобразуются к виду:





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

Из уравнения (1) будем иметь:



Последнее уравнение представляет собой обыкновенное линейное неоднородное дифференциальное уравнение, решением которого будет уравнение вида:



где J₀ - функция Бесселя, f(х,s) - неизвестная функция.

Применим преобразование Лапласа к уравнениям 9 и 10. Получим:





Исключим из (13) производную

Будем иметь:



Определим неизвестную f(x,S)



Поставив полученное выражение в (15), будем иметь



А из уравнения 14



Примем





Г²(S)=z(S)y(S).

С учетом этого из уравнений (17) и (18) получим следующее соотношение между давлением в емкости 4 и в емкости 3:



Так как давление в емкости 4 является гармоническим незатухающим, то согласно теории управления возможна замена S на i.

С учетом этого будем иметь:



Произведение FL=П - активная пористость образца. Отсюда



Если пористое тело 2 имеет множество ориентированных вдоль оси Х цилиндрических капилляров различного диаметра, то R - эквивалентный диаметр капилляров.

Реальные открытые поры материалов в большинстве своем очень далеки от идеальных, вытянутых вдоль оси Х капилляров различного диаметра. Тем не менее, при фильтрации газа через реальную пористую среду, как и через идеальную, будут наблюдаться трубки тока, которые можно привести с рассмотренной нами физической модели (вытянутым вдоль оси Х цилиндрическим капиллярам различного диаметра) через коэффициент извилистости _изв. Необходимо также учитывать входные и выходные потери газа на входе и выходе из пористого тела. Если эти потери учитывать традиционным методом через уравнения Вейсбаха, то это приведет к нелинейности уравнений. Чтобы избежать этого, применим также широко известный способ учета входных и выходных сопротивлений, заключающийся в искусственном увеличении длины входных и выходных участков пористого тела в направлении фильтрации. С учетом этого длина L в уравнении 24 будет равна

L = L_мизвL_вх+L_вых = L_м, (25)

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

_изв- коэффициент, учитывающий извилистость пор;

L_вх - искусственное удлинение участка пористого тела для того, чтобы учесть потери на входе;

L_вых - искусственное удлинение участка пористого тела для того, чтобы учесть потери на выходе;

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

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

Вследствие этого при определении П, и R необходимо совместно решить три уравнения 24 при различных и А_а/A_n.

Так как потери давления на входе и выходе образца согласно уравнениям Вейсбаха пропорциональны квадрату скорости фильтрации газа в порах образца, а скорость фильтрации для одних и тех же образцов пропорциональна перепаду давления на образце и обратно пропорциональна толщине образца L_м, в трех уравнениях 24 отношение амплитуд А_а/A_n должны соответствовать тем значениям и L_м, при которых отношение перепада давления на образце к толщине образца постоянно (P/L_м = const, или при постоянной P = const).

На фиг.2 представлены характерные осциллограммы измерения А_а/A_n и P/L_м в зависимости от частоты при L_м=const.

Характер изменения P/L_м позволяет, как это видно из графика, найти три значения частоты ^*, ^**, ^***, при которых P/L_м = const. Извилистость кривой P/L_м = f() объясняется наличием функций Бесселя и эллиптических синусов и косинусов в уравнении 24. Поэтому при любых материалах и любых L_м всегда можно найти три значения частоты , при которых P/L_м = const и P = const.. Так как L_м в процессе измерения остается неизменной, то для упрощения находим три значения частоты при P = A_a-A_n =const.

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

Устройство состоит из цилиндрической, выполненной из газонепроницаемого материала обечайки 1, в которую плотно установлен цилиндрический пористый образец 2 таким образом, что длина образца 2 равна длине обечайки 1. Один открытый конец обечайки 1 соединен, например, с помощью резьбового соединения и уплотнителя (условно не показаны) с калиброванной емкостью 3, а другой с емкостью 4, являющейся генератором гармонических колебаний. Для возбуждения колебаний в емкость 4 вмонтирован, например, сильфон 5. В стенках емкостей 4 и 3 установлены датчики 6 и 7, выходы которых соединены с входом ЭВМ 9 через электронное согласующее устройство 8.