УДК 539.219.3:303.725.23
ВЫБОР ОПТИМАЛЬНЫХ УСЛОВИЙ ИЗМЕРЕНИЯ ХАРАКТЕРИСТИК ВЛАГОПЕРЕНОСА В ТОНКОЛИСТОВЫХ КАПИЛЛЯРНО-ПОРИСТЫХ МАТЕРИАЛАХ МЕТОДОМ «МГНОВЕННОГО» ИСТОЧНИКА ВЛАГИ
С.Н. Мочалин, И.Н. Исаева, С.В. Пономарев
Кафедра «Управление качеством и сертификация», ГОУВПО «ТГТУ»;
kafedra@asp.tstu.ru
Ключевые слова и фразы: влагоперенос; влагосодержание исследуемого материала; коэффициент влагопроводности; коэффициент диффузии; относительные среднеквадратичные погрешности; плотность материала; тонколистовой капиллярно-пористый материал.
Аннотация: Предложены расчетные зависимости и алгоритмы для вычислений значений коэффициента диффузии влаги, плотности сухого материала и коэффициента влагопроводности тонколистового капиллярно-пористого материала. Разработанные математические модели погрешностей измерений использованы при определении оптимальных режимных параметров метода измерения и рациональных конструкционных размеров измерительного устройства для экспериментального измерения характеристик переноса влаги.
Введение
Целью данной работы является выбор оптимальных режимных параметров и конструкционных размеров устройства для проведения измерений коэффициента диффузии влаги am, плотности р0 сухого материала и коэффициента влагопроводности Xm тонколистового капиллярно-пористого материала методом «мгновенного» источника влаги.
1. Измерительное устройство и основные измерительные операции, осуществляемые при экспериментальном определении характеристик влагопереноса в тонколистовых капиллярно-пористых материалах
При измерении характеристик переноса влаги разработанным методом используют измерительное устройство (рис. 1), основой которого являются две пластины 1 и 2, между которыми размещают тонколистовой капиллярно-пористый материал 3.
В центре верхней пластины 1 просверлено отверстие 4, предназначенное для подведения необходимого количества влаги к исследуемому тонколистовому капиллярно-пористому материалу 3. На нижней поверхности верхней пластины 1 через специально просверленные небольшие отверстия установлены датчики вла-госодержания в виде пары электродов, изготовленных из медной и оцинкованной проволочек. При этом возможны 3 варианта (I, II и III) выполнения этих электродов. В I варианте электроды 5 и 6 размещают вдоль одной окружности радиуса г так, чтобы расстояние между ближними концами этих электродов 5 и 6 было порядка 1...2 мм. Во II варианте основной электрод 7 размещают на окружности
радиуса г, а электрод 8 - на окружности радиуса г’ < г, причем расстояние между электродами 7 и 8 также находятся в пределах 1...2 мм. В III варианте оба электрода 9 и 10 выполняют в виде иголок, размещенных на окружности радиуса г, причем при сборке измерительного устройства эти электроды-иголки прокалывают исследуемый материал 3, а заостренные концы электродов-иголок 9 и 10 входят в полость 11, профрезерованную со стороны верхней поверхности нижней пластины 2.
В процессе эксперимента обеспечивают хороший контакт верхней 1 и нижней 2 пластин с исследуемым материалом 3 за счет использования специального прижимающего механизма, который в простейшем случае может быть выполнен в виде груза 12, равномерно распределенного на верхней поверхности пластины 1 .
Измерительные операции согласно предлагаемому методу осуществляют в основном в следующем порядке:
а) путем выдерживания исследуемого образца в эксикаторе (в который добавлено необходимое количество влаги) задают необходимое начальное влагосо-держание По каждого исследуемого образца;
б) после извлечения из эксикатора каждый образец взвешивают и регистрируют его массу М;
в) образец исследуемого материала 3, имеющий начальное влагосодержание П0, размещают между верхней 1 и нижней 2 пластинами измерительного устройства;
г) за счет использования прижимного устройства создают постоянную для всех экспериментов силу, обеспечивающую хороший контакт медных и цинковых электродов 5 и 6 (либо 7 и 8, либо 9 и 10) датчика влагосодержания с исследуемым тонколистовым капиллярно-пористым материалом 3;
д) до начала активной стадии эксперимента измеряют сигнал датчика влагосодержания; этот сигнал Ео представляет собой электродвижущую силу, разви-
ваемую парой электродов 5 и 6 (либо 7 и 8, либо 9 и 10), и зависит от начального влагосодержания По образца;
е) в момент начала активной стадии эксперимента через канал 13 в грузе 12 и отверстие 4 в центре верхней пластины 1 вводят заданное количество влаги (в виде капли воды из дозирующего шприца);
ж) на протяжении всей активной стадии эксперимента с заданным шагом во
времени Дтш в моменты времени тг = хг-_1+Дхш (/ = 1, 2, ..., п; = 0) измеряют и
регистрируют значения Е, (/ = 1, 2, ., п) сигнала датчика влагосодержания;
з) эксперимент прекращают в момент времени, когда сигнал Е, датчика вла-госодержания снижается до величины
Е < в (ЕШах - Ео) + Ео,
где Етах - максимальное значение сигнала датчика влагосодержания, достигнутое в ходе эксперимента в момент времени ттах; Ег- - текущее значение сигнала датчика влагосодержания в момент времени тг-; в - заданное значение из диапазона
0 < в < 1, принятое в качестве критерия окончания эксперимента;
и) искомое значение коэффициента диффузии влаги в тонколистовом капиллярно-пористом материале вычисляют на основе полученной в эксперименте первичной информации путем ее обработки по приведенным ниже зависимостям и алгоритмам.
2. Расчетные соотношения метода «мгновенного» источника влаги
Исходя из математической модели рассматриваемого метода, изложенной в [1], легко получается решение, имеющее вид
„2
/ ч W
U (r, t)-U0 =--------exp
4пХтт
4атт
(1)
где U(r, т) - влагосодержание в цилиндрическом сечении исследуемого материала на расстоянии r (от места подвода влаги к образцу) в момент времени т; Uo - начальное влагосодержание образца (в момент времени т = 0); W - мощность импульсного источника влаги, подействовавшего в момент времени т в начале координат r = 0, вычисляемая как отношение количества влаги (подведенного из шприца-дозатора) к толщине листа исследуемого образца капиллярно-пористого материала.
В работах [1, 2] показано, что определение коэффициента диффузии am можно осуществлять с использованием формулы
r 2
am = , (2)
4т
^ ‘"max
где ттах - момент времени, в который достигается максимальное значение Umax влагосодержания.
Использование формулы (2) приводит к большим погрешностям, поскольку сложно достаточно точно определить значение момента времени ттах. Попробуем определить такие моменты времени т' и т", используя которые можно минимизировать погрешность определения коэффициента диффузии.
Для этого измерим значения влагосодержания в точке с координатой r во времени т и зарегистрируем эту кривую (рис. 2).
Рис. 2. Изменение во времени т влагосодержания [Щ(г, т) - Щ0] в точке с координатой г и иллюстрация непостоянства абсолютной погрешности Ат определения моментов времени в ходе эксперимента при постоянной абсолютной погрешности АЩ измерения влагосодержания
г 2
Введем безразмерную переменную г =----------. Если принять, что величина
4атт
начального влагосодержания По ф 0, то решение (1) примет вид
U (r, т) -Uo =-
-ехр
4amт
W
r 2nPo
rexp[- z] = U (r, z) -Uo,
(3)
где P0 -■
- плотность абсолютно сухого исследуемого материала.
Введем переменный параметр
U (r, т)-U о
Umax - U0
(4)
где П - начальное влагосодержание исследуемого материала, условно принимаемое за начало шкалы влагосодержания в каждом эксперименте; Птах - максимальное значение влагосодержания П(г, ттах), достигаемое в ходе эксперимента в сечении радиуса г образца в момент времени ттах.
Запишем решение (3) для моментов времени т' = т'(а) и ттах:
U (r, т')-U о =
W
nr 2Po
■z(x'(a))exp[-z(x'(a))].
(3a)
U,
тах П0
W
z(l:max)exp[ z(l:max)]
W
2 -max
nr Po
max 2
nr Poe
(3b)
2
r
m
a
m
где принято во внимание, что согласно [3] максимальное значение влагосодержа-ния Umax достигается при z(xmax) = 1.
Поделив (3а) на (3b), получим
а = U(^Т (а))—— = ez(T(a))exp[—z(T(a))]. (4a)
Umax — U 0
Преобразовав выражение (4a) с учетом того, что eexp(z) = exp(1 — z), получим уравнение
z exp(1—z) = a . (5)
Обозначим z ' = z(x'(a)) и z '' = z(x'' (a)) соответственно больший и меньший корни уравнения (5). После преобразований легко получаем формулы для вычисления искомого коэффициента диффузии am по экспериментально измеренным моментам времени т' и х":
2 2 r r
из z' = z(x' (a)) =--- следует, что am =-------------------------------; (6)
4amx' 4z V
rr из z ' ' = z(x' (a)) =--- следует, что am =------------------------------------, (6a)
4amx' ' 4z ' 'x' '
где х' и х" - меньший и больший моменты времени, соответствующие большему і' и меньшему і” корням уравнения (5), при которых достигается заданное значение параметра а, определенное формулой (4).
Для вычисления плотности ро сухого исследуемого материала на основе зависимости (3 а) легко получается формула
W г(т' (а)) ехр[—г(т '(а))]
ро =-----Гг-----;----------—,
п г [и(г, т (а)) — ио ]
которая с учетом (4а) принимает вид
Wa (7)
Ро =—гг---------;----------? (7)
п е г [и(г, т (а)) — ио ]
а при a = 1
W
Ро =------Г. (7а)
п e r [Umax U0]
Отметим, что для вычисления коэффициента влагопроводности Хт = атр0 можно воспользоваться формулами:
r2 W W
Xm =-----------Гг---------ч =-------;----------7 , (8)
4z'x'n e r (Umax-U0) 4nez’x' (Umax-U0)
r2 W W
x m =-------------г-------------ч=----------;------------ч. (8a)
4 z "x'' п e r (Umax - U0 ) 4п e z "x " (Umax - U0 )
3. Вывод соотношений для вычисления погрешностей
По методике, описанной в [3, 4], легко получить выражения для вычисления погрешностей для случая определения (вычисления) коэффициента диффузии по формулам (6) и (6а).
На основании (6) получено следующее выражение
5ат =,14(5г)2 + (5г')2 + (5т1)2 , (9)
~ Аа
где 5а = — - среднеквадратическая оценка относительной погрешности измере-
а
Аг
ния коэффициента диффузии; 5г =--------относительная погрешность определения
г
расстояния г между началом координат (местом ввода влаги из шприца-дозатора)
5 ' Аг'
и измерителем влагосодержания; 5г =— - относительная погрешность опреде-
2 '
' 5 • Ат'
ления величины г ; 5т'=---- - относительная погрешность определения момента
т
времени т'; Аат, Аг, А2, Ат ' - абсолютные погрешности измерения (определения) величин ат, г, 2 и т .
Определим составляющую погрешности 5г'
, йг 1 йг 1 йг 52 и— =--------йаи----Аа .
2 2 йа 2 йа
Аналогично определим абсолютную погрешность Да по методике, изложенной в работах [3, 4]. Учитывая, что а рассчитывается по формуле (4), получаем
Аа = ад/52(и-ио) + 52(итах -V<,) , (10)
5 (и-V о) = Аи—и°1, 5(итах - и о) = А(итах—и°2 - относительные погрешно-
V-ио итах-ио
сти измерений разностей влагосодержаний (V - ио) и (итах - ио).
Предполагая, что А(и -и о) = А(итах -ио) = Аи, после преобразования (1о) получаем
Аа = а5(итах -и о)^+1 , (11)
откуда следует
1 dz’ (a) 1
5z = a^Uma^ 2 +
z'(a) da V a
где 8 Umax = 8(Umax — U0 ) .
Повторив вышеизложенное, легко получаем, что
1 dz” (a) 1
8z =-----------a8UmaxJ—+1 .
z '' (a) da V a
,2
Определим входящую в (9) составляющую 5т', поскольку ясно, что эта составляющая не представляет собой просто число, а изменяется во времени при постоянной абсолютной погрешности АП измерения влагосодержания.
Погрешности определения моментов времени ттах, т' и т'' связаны не только непосредственно с измерением времени. Из представленного на рис. 2 графика видно, что в реальных условиях эксперимента пик на кривой, соответствующий максимуму влагосодержания, заметно зашумлен и поэтому несколько размыт. Если определение значения максимального значения влагосодержания Птах можно провести достаточно точно, то соответствующий момент времени ттах определяется с заметно большей абсолютной погрешностью Аттах, чем погрешности Ат' и Ат'' определения моментов времени т' и т"
Определение моментов времени т' и т'' осуществляется через определение влагосодержания, поэтому погрешность определения этих моментов времени можно выразить через погрешность определения влагосодержания. Учитывая, что
dU(r,x) _ AU dx Ax
----, получаем Ax";
AU
dU Ua) dx
(12)
Ax"
AU
dU Ua) dx
(12a)
Из формулы (12) следует, что абсолютная погрешность Дт' определяется моментом времени т', зависит от абсолютной погрешности Ди измерения влагосо-
СЮ (а)
держания U(a) и от производной
dU Ua)
dx
в момент времени т'. При этом значения
производной
dx
и погрешности A^(a) зависят от выбора конкретного значе-
ния безразмерного параметра a =
U (r, x")-U о
Umax U о
представляющего собой отноше-
' max u 0
ния разностей влагосодержаний (U(r,т') - Uо ) и (Umax - Uо). Если предположить, что абсолютные погрешности измерения влагосодержаний A(Umax - Uo) = = A(U - Uo) = AU = const остаются постоянными, то с учетом (12) получим соотношения:
Ат' AU .(Umax-Uo) SUmax .
= 8x " = -
dU
dx
_ " г dU
1 dx ) x "
AU
dU
dx
(Umax -Uо) SU„
ГdU ^
1 dx ) x " "
где принято во внимание, что AU = (U
Ur, x
dU Ur, x) ' "
Значения производных------------—- для моментов времени т' и т'':
dx
x
x
i
I
X
x
Ф (r, x)-U о]
dx
д[ (r, x)-U о]
W
r 2ПРо
(l- z )exp(- z )
W
r 2про
(l- z' (a))exp(- z ' (a));
dx
W r 2пРо
(l- z )exp(- z )
W
r 2nP0
(l - z '' (a))exp(- z '' (a)),
а для момента времени ттах, когда г = 1:
дП (г, т)
dx
= 0:
Umax U 0
W
nr2РоЄ
Проведем следующие преобразования:
d [U (r,x) -U о ]
dx
W
Umax U 0
Wn r ро
-[l- z' (a)]exp[- z '(a)]:
d (U (r, x)-U о) dx
e n r 2p0W [l-z' (a)]exp[- z ' (a)] e [l-z' (a)]exp[- z ' (a)]
x "
(Umax - U0 )S Umax =
SUm
т ' Г д[ (r, t)-U o]
V дт Действуя аналогично, получаем
(Umax - U0 )S Umax
x " e [l - z ' (a )]exp [- z ' (a )]
SUm
, " Г т (a)1 x"" e [l-z'" (a)]exp[-z '' (a)]
1 дг ) x=x " "
В итоге получаются следующие формулы для определения среднеквадратичных погрешностей измерения коэффициента диффузии:
(Sam ) =•
4S2r +
SUm
x" e[l—z ' (a)]exp(-z ' (a))
1 dz’ (a)
-------------—-a
z' (a) da
1+_2 |SUmax
T=X
I
X=T
x=x
тах
l
2
2
+
(Sam) =
4S2r +
SUm
x "" e [l — z '' (a)]exp(-z '' (a))
1 dz'' (a)
-------------——a
" (a) da
1+_2 |s2umax a2
По аналогии с изложенным выше, на основе формулы (7) была получена зависимость для вычисления среднеквадратичной оценки относительной погрешности 5ро измерения плотности ро сухого материала:
2
2
+
z
8(р0) = ^52Ж + 452г + 52а + (5[и(г,т '(а)) - и0 ])2 .
Принимая во внимание, что на основании (11) 5а=^а=5итах I——+1, а
а Уа2
8 [и(г,т')-и0] = 5^тах
получаем, что
а
8(Ро) = ^52Ж + 452г + (5 итах )2 ^0? +1^)
где 5Ж- относительная погрешность измерения мощности источника влаги Ж, которая мгновенно выделяется в единице толщины исследуемого материала в начальный момент времени т = 0; 5итах - относительная погрешность измерения максимальной разности влагосодержаний (итах- ио); 5г - относительная погрешность измерения расстояния г между центром отверстия 4 (местом введения влаги) и измерителем влагосодержания (см. рис. 1).
4. Оценка погрешностей измерений
По полученным формулам были рассчитаны среднеквадратичные погрешности для различных значений а, Ж, г и выбраны оптимальные параметры, обеспечивающие наименьшую погрешность.
При выполнении расчетов были использованы следующие значения: ат = 5-10-9м2/с; ро = 300 кг/м3; Ж = 0,01 кг/м; АЖ = 0,0001 кг/м; г = 4...6 мм; Аг = 0,2 мм; Аитах = 0,01 кг/кг. Пример зависимости среднеквадратичных погрешностей 5а'т, 5а'т, 5(р0) от параметра а при г = 4 мм и Ж = 0,01 кг/м приведен на рис. 3.
Рис. 3. Зависимость относительных погрешностей измерения коэффициента диффузии Ъа'т, Ъа'т и относительных погрешностей измерения плотности сухого материала 8(р0) от значения а для Ж = 0,01 кг/м, г = 4 мм
Как видно из графиков, приведенных на рис. 3, оптимальные значения а (для измерения коэффициента диффузии ат) находятся в диапазоне 0,4.0,6. После проведения более точных расчетов при изменении а с шагом 0,01 были определены минимальные значения погрешностей и соответствующие им значения а:
(5ат ) = 12,4 % при а' = 0,48 и (5ат ) = 18,0 % при а" = 0,64.
Из рис. 3 видно, что значения относительных погрешностей 5а'т измеренного коэффициента диффузии ат при его вычислении по большему моменту времени т" оказывается заметно выше по сравнению со значениями 5а'т в случае вычисления коэффициента диффузии по величине меньшего момента времени т'.
Это объясняется тем, что абсолютная погрешность Ат '' *
AU
dU (а) dx
определения
момента времени т" оказывается значительно больше абсолютной погрешности
Ат '!
AU
dU (а) dx
так как производная
dU (а) dx
при т" > ттах заметно меньше
производной
dU (а) dx
при т' > ттах, что вытекает из формул (12), (12а) и из
характера изменения влагосодержания во времени (см. рис. 2).
Понятно, что от использования большего по величине момента времени т'' следует отказаться, а при вычислении коэффициента диффузии при обработке экспериментальных данных надо использовать меньшее значение момента времени т', соответствующее оптимальному значению параметра а ~ 0,48.
В процессе выполнения данного исследования стало очевидно, что наряду с выбором оптимального значения параметра а = аопт = 0,48, следует определить
т
т
О 0.2 0.4 а 0;6 0.8 1
Рис. 4. Линии равного уровня погрешностей измерения коэффициента диффузии ат и плотности сухого материала р0
(при использовании метода «мгновенного» источника влаги)
величину расстояния г, при котором погрешность измерения коэффициента диффузии будет минимальной. В рамках этой части данного исследования были рассчитаны значения погрешностей измерения коэффициента диффузии при значениях расстояния г из диапазона 2 мм < г < 6 мм.
В результате были построены линии равных уровней погрешностей измерения коэффициента диффузии ат и плотности сухого материала р0, представленные на рис. 4.
Заключение
Из материалов, изложенных в данной статье, следует, что минимальное значение относительных среднеквадратичных погрешностей Ъат измерения коэффициента диффузии достигается при а = 0,48, г = 4 мм, приемлемые значения погрешностей имеют место при 0,3 < а < 0,65; 3,75 мм < г < 4,25 мм.
Для минимизации среднеквадратичных погрешностей 5р0 измерения плотности р0 сухого материала ее нужно измерять в тот момент времени, когда влаго-содержание достигает максимума, то есть при а = 1. Тогда при Ж = = 0,01 кг/м; г = 4 мм и Аг = 0,2 мм; Аитах = 0,01 кг/кг; АЖ = 0,0001 кг/м, получаем 5р0тш = 12,3 % при аопт = 1 и гопт = 4 мм.
Список литературы
1. Мочалин, С.Н. Математическая модель экспрессного метода измерения коэффициента диффузии влаги в капиллярно-пористых материалах / С.Н. Мочалин // Теплофизика в энергосбережении и управлении качеством : материалы шестой междунар. теплофиз. школы : в 2 ч / Тамб. гос. техн. ун-т. - Тамбов, 2007. -Ч. II. - С. 160-161.
2. Мочалин, С. Н. К вопросу о выборе оптимальных режимных параметров метода измерения коэффициента диффузии влаги в тонколистовых капиллярнопористых материалах / С.Н. Мочалин, С.В. Пономарев // Междунар. науч.-техн. семинар «Актуальные проблемы сушки и термовлажностной обработки материалов» : материалы семинара / Федер. агентство по образованию, ГОУ ВПО «Воронеж. гос. лесотехн. акад.». - Воронеж, 2010. - С. 63-67.
3. Теоретические и практические основы теплофизических измерений / С.В. Пономарев [и др.] ; под ред. С.В. Пономарева. - М. : Физматлит, 2008. - 408 с.
4. Пономарев, С.В. Теоретические и практические аспекты теплофизических измерений : монография. В 2 кн. Кн. 2 / С.В. Пономарев, С.В. Мищенко, А.Г. Дивин. - Тамбов : Изд-во Тамб. гос. техн. ун-та, 2006. - 216 с.
The Choice of Optimal Measuring Conditions of Moisture Transfer in Thin Sheet Capillary Porous Materials by the Technique of “Prompt” Moisture Source
S.N. Mochalin, I.N. Isaeva, S.V. Ponomarev
Department “Quality Management and Certification”, TSTU
Key words and phrases: moisture transfer; moisture content of the examined material; moisture transfer coefficient; diffusion coefficient; relative mean-square errors; material density; thin sheet capillary porous material.
Abstract: The paper offers the calculations, dependencies and algorithms for computation of moisture diffusion coefficient in a thin sheet capillary porous material. The mathematical models of measuring errors are used to determine of the optimal mode parameters for the measuring technique and rational design sizes of the device for moisture transfer characteristics measurement.
Auswahl der optimalen Bedingungen der Messung der Charakteristiken der Feuchtetibertragung in den feinen kapillarporosen Stoffen durch die Methode der “augenblicklichen” Feuchtigkeitsquelle
Zusammenfassung: Es sind die berechneten Abhangigkeiten und die Algorithmen fur die Berechnung der Werte des Koeffizientes der Feuchtigkeitsdiffusion, der Dichte des Trockenstoffes und des Koeffizientes der Feuchteleitfahigkeit des feinen kapillarporosen Stoffes vorgeschlagen. Die erarbeiteten matematischen Modelle der Meftfehler wurden bei der Bestimmung der optimalen Regimeparameter der Meftmethode und der rationalen Konstruktionsgroften der Mefteinrichtung fur die experimentellen Messung der Charakteristiken der Feuchteubertragung verwendet.
Choix des conditions optimales de la mesure des caracteristiques du transfert de la chaleur dans les toles fines capillaires et poreuses par la methode de la source «instantanee» de l’humidite
Resume: Sont proposes les calculs, les dependances et les algorithmes pour la mesure de la grandeur du coefficient de la diffusion de l’humidite dans une tole fine capillaire et poreuse. Les modeles mathematiques des erreurs des mesures elabores sont utilises lors de la definition des regimes optimaux des parametres de regime de la methode de mesure et des dimensions rationnelles de construction du dispositif pour la mesure des caracteristiques du transfert de la chaleur.
Авторы: Мочалин Сергей Николаевич - аспирант кафедры «Управление качеством и сертификация»; Исаева Ирина Николаевна - аспирант кафедры «Управление качеством и сертификация»; Пономарев Сергей Васильевич - доктор технических наук, профессор кафедры «Управление качеством и сертификация», ГОУ ВПО «ТГТУ».
Рецензент: Коновалов Виктор Иванович - доктор технических наук, профессор кафедры «Химическая инженерия», ГОУ ВПО «ТГТУ».