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

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

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

где

Г²(i)=z(i)y(i);

где *, **, *** - соответственно три значения частоты, при которых обеспечивается одинаковое значение А_а-А_n=Р

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

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

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

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

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

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

Р_r - число Прандтля;

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

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

- коэффициент теплопроводности;

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

С_р - удельная теплоемкость при постоянном давлении;

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

а₂, b₂ -заданные const;

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

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

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

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

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

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

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

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

Наиболее близким к предлагаемому изобретению является способ определения активной пористости материалов (Патент №2180743, кл. G 01 N 15/02, 2002 г.), включающий ламинарную фильтрацию газа через заключенный в обечайку из газонепроницаемого материала образец из замкнутого пространства перед образцом в замкнутое пространство после него путем создания в замкнутом пространстве перед образцом гармонических колебаний давления газа изменяющейся частоты и определенной амплитуды и измерение амплитуды колебаний давления газа в замкнутом пространстве после образца. Измеряют изменение амплитуды давления газа в замкнутом пространстве до и после образца, вычисляют разность амплитуд и определяют три значения частоты из условия равенства разности амплитуд, определяя активную пористость материалов, эквивалентный радиус пор и коэффициент, учитывающий извилистость пор и потери давления на входе в образец и выходе из него в результате совместного решения на ЭВМ трех уравнений:

,

где

у(i )=-

L= Lм

где *, **, *** соответственно три значения частоты, при которых обеспечивается одинаковое значение А_а-А_n= Р (или

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

;

;

где

;

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

, , , L= L_м

где *, **, *** соответственно три значения частоты, при которых обеспечивается одинаковое значение А -А_n= Р (или

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

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

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

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

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

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

Р_r - число Прандтля;

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

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

- коэффициент теплопроводности;

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

С_р - удельная теплоемкость при постоянном давлении;

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

а₂, b₂ - заданные const;

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

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

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

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

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

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

- уравнение Навье-Стокса:

- уравнение неразрывности:

- уравнение сохранения энергии:

- уравнение Ван-дер-Ваальса:

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

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

где р - давление;

Т - температура;

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

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

- плотность;

- коэффициент теплопроводности;

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

х, r, - цилиндрические координаты;

t - время;

С_р - удельная теплоемкость при постоянном давлении;

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

₀, Т₀ - значения , Т при установившемся режиме;

а₂, b₂ - заданные const.

При а₂=b₂=0 получаем уравнение состояния Менделеева-Клапейрона идеального термодинамического газа:

Р=R₀T

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

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

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

- число Прандтля,

К системе (1, 3, 6, 7) добавим граничные условия, которые получаются при следующих допущениях.

1) Давление во входном сечении равняется давлению Р_а, а давление на выходе канала равняется давлению Р_n.

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

Итак, граничные условия

где F - площадь поперечного сечения образца.

Применим к системе (1), (3), (6), (7) преобразование Лапласа по переменной t, тогда получим систему уравнений

где

Первые два уравнения - обыкновенные линейные неоднородные дифференциальные уравнения.

При их решении общее решение уравнения есть сумма общего решения соответствующего однородного и частного решения данного неоднородного уравнения. Соответствующее однородное уравнение как для (11), так и для (12) есть уравнение Бесселя. Учитывая его, а также определяя частное решение, входящее в правую часть, мы получили общие решения первых двух уравнений

где J₀ - функция Бесселя, f(x,s), F(x,s) - произвольные функции, причем

Прежде чем решать третье и четвертое уравнение системы, т.е. (13), (14), мы исключили из них производные используя найденные решения (15) и (16).

После ряда преобразований получим:

Подставляя (18) и (19) в (13) и (14) преобразуем третье и четвертое уравнения, которые примут вид

где

Исключая из (20) и (21) поочередно и получили уравнения

где

Решая уравнение (24), (25) методом Эйлера получим:

где

,

а А(s) и В(s) пока неизвестные функции.

Итак, решением системы (11), (12), (13) и (14) являются функции (15), (16), (27), (28), куда входят неизвестные функции F(x,s), A(s), B(s).

Для определения А(s) и В(s) используем условия (9) и (10). После применения к ним преобразования Лапласа они примут вид:

при х=0

при х=0

Подставляя (27), (28), (29) и (30) и (31) и (32) получим:

А(s)+В(s)= (s),

А(s)е^Г(s)L+В(s)e^-Г(s)L= (s),

Имеем три уравнения с двумя неизвестными А(s) и В(s). (t) можно считать известным, т.к. это давление в измерительной камере, которое можно измерить.

Из первого уравнения выразим В(s) через А(s). Из второго и третьего уравнений найдем А(s), и, наконец, подставив А(s) и В(s) во второе уравнение выразим давление в измерительной камере (s) через внешнее давление (s).

где

chГ(s)L= shГ(s)L=

Далее предполагая, что изменение давления в рабочей камере происходит по гармоническому закону, осуществляем замену переменной s на i .

В результате получаем

где

Если в ячейке плотно установить пористое тело с порами в виде ориентированных вдоль оси х цилиндрических капилляров, то с учетом того, что FL=П_а, (где F - суммарная площадь поперечных сечений капилляров, L - длина капилляров, П_а - активная пористость исследуемого материала) уравнение (34) преобразуется к виду

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

При фильтрации газа через пористую среду, будут наблюдаться потоки газа, которые можно привести к рассмотренной нами физической модели (вытянутым вдоль оси Х цилиндрическим капиллярам различного диаметра) через коэффициент поровой извилистости а_изв. При этом L в уравнении (38) будет равна

где L_м - толщина пористого материала в направлении фильтрации.

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

где _вх и _вых - коэффициент входных и выходных потерь;

Р_вх - давление газа до входа в пористое тело;

Р_вых - давление газа за пористым телом;

и - соответственно давление газа на входном и выходном сечении пор;

и - скоростной напор во входном и выходном сечении пор.

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

Он заключается в искусственном увеличении длины входных и выходных участков пористого тела в направлении фильтрации до тех пор, пока не выровняются относительные (отнесенные к единице длины) потери давления на входных, выходных и внутренних участках пористого материала. С учетом этого длина капилляров L в уравнении (38) будет равна

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

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

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

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

Для определения активной пористости материалов в камере 3 сильфоном 4 создают колебания газа. В результате фильтрации газа на образце возникает давление Р=А_а-А_n. Меняя частоту колебаний давления газа и фиксируя разницу амплитуд, определяют три значения частоты, *, **, ***, при которых перепад давления на образце постоянен Р=А_а-А_n=const. Значения частот *, **, *** и отношения амплитуд А_n*/А_a*, А _n**/А_a**, А_n***/А_a*** подаются на ЭВМ 7.

Данные о номинальном объеме измерительной камеры V_n, температуре Т, площади F образца, разделяющего камеры, газовой постоянной R₀, коэффициенте кинематической вязкости , плотности газа , толщине образца L, удельной темплоемкости газа С_р , введены в ЭВМ заранее. Поэтому ЭВМ 7 осуществляет определение активной пористости материалов, эквивалентного радиуса пор и коэффициента, учитывающего извилистость пор и потери давления на входе в образец и выходе из него из совместного решения трех уравнений.

;

;

Пример: осуществляли измерение активной пористости материалов, эквивалентного радиуса пор и коэффициента, учитывающего извилистость пор сложной керамической системы (бетон). Образец выполнен в виде цилиндра диаметром 30 мм и толщиной 15 мм. Расход компонента связки, К_с 300 кг на 1 м³ раствора. Отношение затворителя к компоненту связки З/К_с=0,63. ЭВМ 7 зафиксировал значения частот *=2,49 Гц, **=4,11 Гц, ***=8,56 Гц и отношение амплитуд А_n*/А_a *=0,109, А_n**/А_a**=0,081, А_n ***/А_a***=0,058 и выдал значение П_а=11,5311%. R=13,9875 мкм, =5,1124.